key: cord-0878187-lu7pofuq authors: Imran, Muhammad; Iqbal, Sana; Hussain, Ajaz; Uddin, Jalal; Shahzad, Mohsin; Khaliq, Tanwir; Razzaq Ahmed, Abdul; Mushtaq, Laiba; Kashif, Muhammad; Mahmood, Khalid title: In silico screening, SAR and Kinetic studies of Naturally occurring flavonoids against SARS CoV-2 main Protease date: 2021-10-11 journal: Arabian journal of chemistry DOI: 10.1016/j.arabjc.2021.103473 sha: 41c679864dd191cf4788e26bb0ecc50442205218 doc_id: 878187 cord_uid: lu7pofuq The Severe Acute Respiratory Syndrome Coronavirus (SARS-CoV-2) pandemic has become a global challenge based on its replication within the host cells that relies on non-structural proteins, protease (Mpro). Flavonoids, an important class of naturally occurring compounds with medicinal importance, are frequently available within fruits and vegetables. Herein, we report the in silico studies on naturally occurring flavonoids consisting of molecular docking studies and evaluation of theoretical kinetics. In this study, we prepared a library of nine different classes of naturally occurring flavonoids and screened them on Autodock and Autodockvina. The pharmacokinetic properties of most promising compounds have been predicted through ADMET SAR, inhibition constants, ligand efficiency and ligand fit quality have been worked out theoretically. The results revealed that naturally occurring flavonoids could fit well in the receptor's catalytic pocket, interact with essential amino acid residues and could be useful for future drug candidates through in vitro and in vivo studies. Moreover, MD simulation studies were conducted for two most promising flavonoids and the protein-ligand complexes were found quite stable. The selected natural flavonoids are free from any toxic effects and can be consumed as a preventive measure against SARS CoV-2. The severe acute respiratory syndrome corona virus-2 (SARS-CoV-2) belongs to the family Coronaviridae and is known to cause respiratory disorders and digestive tract infections (Hui et al., 2020) in both humans and animals. In humans, coronavirus infection caused acute lung injury that changed to ARDS (acute respiratory distress syndrome) and resulted in the COVID-19, a condition associated with person-to-person transmission. Emerging and reemerging new coronavirus strains have become a global challenge for public health concerns and become a global pandemic. Different countries adopt various strategies to treat COVID-19, including social distancing, self-quarantine, and anti-HIV and anti-malarial drugs (Colson et al., 2020 , Musarrat et al., 2020 . The high rate of recombination and novelty of the 2019 n-CoV are major reasons for its abrupt damage and widespread. Hence to combat the COVID-19, the search for new antiviral agents is of great interest rather than vaccine therapy. The SARS-CoV-2 is an enveloped virus containing a non-segmented single-stranded positive-sense RNA, and the envelope possesses the glycoprotein projections. The genome of n- CoV revealed six open reading frames (OFRs), and the first OFRs encode the two large proteins, polyprotein 1a, and polyprotein 1b (pp1b). These polyproteins being responsible for viral replication and transcription, are processed by Chymotrypsin like Protease (3CL pro ) or Main Protease (M pro ) and Papain-like Protease into sixteen non-structural proteins (16 nsps). The remaining OFRs encode the four structural proteins, including spike (S), membrane (m), envelope (E), and nucleocapsid (N) proteins (Ziebuhr et al., 2000 , Tufan et al., 2020 . The spike protein binds to the host cell, namely ACE2 (angiotensin-converting enzyme 2), and leads the entry of SARS-CoV-2 to the host cell (Walls et al., 2020) . Precisely speaking, papain-like protease (PL pro ) and 3-chymotrypsin-like Protease (3CL pro ) or M pro is indispensable for CoV replication and also play a role in the inhibition of immune responses (Báez-Santos et al., 2015 , Amin et al., 2020 . Thus, any interruption in the main Protease is of prime significance for blocking viral replication. The proteases have very little similarity with human proteases and make the drugs less toxic to humans (Anand et al., 2002) . Hence targeting the main Protease appears crucial for COVID-19, including SARS (Drosten et al., 2003 , Zhong et al., 2003 and MERS-CoV (Zaki et al., 2012) treatment. The absence of specific therapy for COVID-19 leads to medicinal herbs in several regions of the world. Flavonoids are natural phytochemical compound with several biological activities such as anti-inflammatory, analgesic, antibacterial, antifungal, antiviral, anticarcinogenic, carbonic anhydrase and cholinesterase inhibitory properties (Havsteen, 1983 , Pietta, 2000 , Imran et al., 2020 . Although flavonoids have been reported to have antiviral activity against several types of viruses, the molecular mechanism was rarely known. The compounds like Kaempferol, Narigenin, Epigenin, and Quercetin show a pharmacophore similar to nelfinavir, a protease inhibitor, effective against HIV (human immunodeficiency virus) (Dabeek and Marra, 2019 , Salehi et al., 2019 , Yamamoto et al., 2004 . ACE2 is a receptor of the coronavirus; thus, compounds that can inhibit the activity of the ACE2 enzyme are also considered helpful in COVID-19 treatment . Some citrus flavonoids have been reported to bind the ACE2 enzyme [18] . In most severe COVID-19 patients, a cytokine storm was observed (Huang et al., 2020b) . The host immunity's antiviral response involves producing several pro-inflammatory cytokines and activating the T-cells, which are essential for controlling the viral replication, inflammation, and cleaning the cells (Li et al., 2020, Ivashkiv and Donlin, 2014) . However, the excessive immune responses lead to cytokine stormcytokine storm formation, known as macrophage activation syndrome (MAS), to further tissue damage (Tufan et al., 2020) . The cytokine storm in COVID-19 patients blocks airways, changes vascular permeability and is a cause of exacerbation and even death (Huang et al., 2020a , Xu et al., 2020 . Recent studies on flavonoids showed that flavonoids could be the active ingredients to suppress cytokine release, eliminate inflammation and alleviate the immune response . In this study, a flavonoids library consisting of nine different subcategories of flavonoids (consisting of 51 natural flavonoids) was prepared to evaluate in silico interactions and binding affinities against SARS-CoV-2 main Protease and structure-activity relationship (SAR). The physicochemical ADMET (absorption, distribution, metabolism, excretion, and toxicity) properties were also worked out. The kinetic aspects of these compounds have also been worked out. Moreover, the food sources consumed as a preventive measure against SARSCoV-2 have been figured out from our results. The docking studies were carried out through Vina (Trott and Olson, 2010) interlinked with MGL tools (Morris et al., 2009) . The binding pattern was explored using visualization tools, including Biovia discovery studio (Studio, 2008) and Pymol (DeLano, 2002) . The compounds possessing efficient binding affinity were selected for further analysis of absorption, distribution, metabolism, excretion profiling through AdmetSAR (Cheng et al., 2012) , http://lmmd.ecust.edu.cn/admetsar2/, a ligand property evaluating tool. The ligands three-dimensional structures were obtained from the PubChem database (National centre for biochemistry Information, U.S. National Library of Medicine) https://pubchem.ncbi.nlm.nih.gov/ in .sdf format, and structures were drawn on Marvin sketch (https://chemaxon.com/products/marvin) after the literature survey. The molecular energy optimization of ligands was performed using UCSF Chimera (Pettersen et al., 2004) . AutoDock tools (ADT1.5.6) [28] were used to set the number of torsions for the flavonoids and COVID-19 3CL pro /Mpro, and vina files were prepared and saved in pdbqt format. The crystal structure of COVID-19 3CL pro /Mpro complexed with inhibitor N3 (PDB ID: 6lu7) was downloaded from the protein data bank and used as a target in the screening of selected flavonoids. In preparing the receptor co-crystallized ligands, chain B and water molecules were excluded, and Gasteiger charges were added for a protein using UCSF chimera (Pettersen et al., 2004) . The grid size was set to 60 x.y.z points having grid spacing of 0.375 A , and × 60 × 60°g rid recenter (x, y, z): -9.732, 12.695, 69.958, respectively was designated. The NCBI conserved domain (https://www.ncbi.nlm.nih.gov/Structure/cdd/wrpsb.cgi database was used to showed the conserved residues of main protease structure. Aminoacids that are present in the vicinity of 6.5Å from the N3 inhibitor are considered as main binding residues. The residues such as Thr24, His41, Cys44, Phe140, Thr26, Gly123, His164, Leu27, Glu166, His 172, Leu141, Thr190, Asn142, Gln189 and Gln192 are present in the catalytic pocket of main protease (Figure2). Molecular docking studies between the retrieved compounds and main protease were carried out by using AutoDock Vina interlinked with AutoDock MGL tools. The binding pocket of the receptor was well selected using the grid box generation process of AutoDock MGL tools. All the optimized flavonoids are then subjected to docking process, to evaluate their binding pose and binding energy in the active site of main protease. The flavonoids with best confirmations and binding affinity with main protease are evaluated by intermolecular interactions. The compounds with best docking score (i.e. Amentoflavone and Morusin) were further used for molecular dynamic simulation studies. Molecular dynamic was performed for Amentoflavone-M pro complex and Morusin-M pro complex using Desmond module of Schrodinger suite. The complexes were simulated for 100ns of timescale for understanding the interaction pattern and dynamic behavior (Guo et al., 2010) . The prepared complexes were solvated using orthorhombic solvent box conditioned for 10Å buffer region (Jorgensen et al., 1983) . The system in each case was neutralized by adding desired number of counter ions, and the total system consisted of 39326 atoms and 7697 water molecules. An ensemble (NPT) of Nose-Hoover thermostat used to maintain the constant temperature of 300K and Martyna-Tobias-Klein barostat was applied to maintain the pressure (Martyna et al., 1992) . 1000 steps of steepest descent energy minimization followed by conjugate gradient algorithms were utilized. Smooth particle mesh Ewald method was used to evaluate the short range interactions with cut-off radius of 9 Å and long range coulombic interactions were estimated. Reference system propagator algorithms (RESPA) integration was used in the dynamic study of non-bonded interactions. Trajectories were documented for every 100pc, and the obtained results was analyzed with Maestro graphical interphase. NCBI CD results showed that catalytic dyad of M pro composed of two residues like His41 and The flavonoids in nine different scaffolds having the lowest binding energies and highest bonding interactions in the main Protease's catalytic pocket are presented in Table 1 . Moreover, Lipinski's rule of five (RO5) and Veber's rule have been employed to determine the drug-likeness and oral activity of compounds under study, respectively. (Lipinski et al., 1997 , Lipinski, Giménez et al., 2010 . It has been found that twelve flavonoids show no violation of RO5 and represent the druglike molecular nature. Their log P values are within the range of -0.79 to 4.54. Molecular weight, hydrogen bond donor, and hydrogen bond acceptor are within the accepted range except for Amentoflavone and Naringin, which showed two and three violations. Amentoflavone showed a violation of molecular weight and the number of hydrogen bond donors, while Naringin showed a violation of molecular weight, the number of hydrogen bonds, and the acceptor. All the selected flavonoids except Amentoflavone fulfilled the criteria of Veber's rule (Veber et al., 2002) with the number of rotatable bonds and TPSA (total polar surface area) within the range for oral bioavailability. The number of rotatable bonds for drug-like candidates should be less than seven to improve oral bioavailability (Veber et al., 2002) . A polar surface area equal to or less than 140 Å is essential for good oral bioavailability. Oral bioavailability is the fractional extent of administered drug that finally reach the site of action (Benet et al., 1996 , Kim et al., 2014 The flavonoid's pharmacokinetic properties of all the compounds were calculated using SwissADME (http://www.swissadme.ch/) as shown in Table 2 . All flavonoids showed good gastrointestinal absorption, with the probability of the absorption of 96 to 98%. The low absorption rate of the molecules is correlated with their high molecular weight. However, despite its high molecular weight, the Amentoflavone showed an increased absorption of 97%. The solubility range lies between -1 to -5 (Tsaioun and Kates, 2011) , and the aqueous solubility of the flavonoids ranges from -3.78 to -2.77 presented moderate to high solubility. The high solubility of flavonoids is due to hydroxy groups that participate in hydrogen bonding with water molecules. The blood-brain barrier is a semi-permeable border present within endothelial cells. It prevents the passage of pathogens and large or hydrophilic molecules into the extracellular fluid of the CNS (central nervous system) passage of small polar and hydrophobic molecules (Obermeier et al., 2016) . BBB (blood-brain barrier) is the layer of endothelial cells that prevents the entry or exit of solutes from blood stream to the extracellular space. It may be easily evaluated by considering the compound's concentration in the brain's concentration in blood (Obermeier et al., 2016) . In the central nervous system, active drugs and high penetrations of the blood-brain barrier are necessary. In contrast, in non-central nervous system active drugs, low penetration reduces the CNS side effects. The selected flavonoids are not good blood-brain barrier permeant. ADMET properties of a useful drug-like compound accepted the following parameters: such as i) absorb the CaCO 3 permeability ii) non-toxic iii) non-carcinogenic iv) it must absorb in the human intestine (HI) it must pass the Blood-brain barrier. In the cytochrome p450 inhibitory activity, almost all the selected flavonoids are an inhibitor of CYP to show the drug-drug interactions. CytochromeP450 (CYP) is a family of an enzyme involved in the metabolism of polyunsaturated fatty acids, metabolism of endogenous and exogenous chemicals, and metabolism of toxic compounds, including drugs bilirubin (Berka et al., 2011 , Waxman et al., 1991 . The apigenin, wogonin, and morusin show the high human colonic adenocarcinoma cell line Caco-2 permeability than other selected compounds. In the case of acute oral Toxicity Apegenin, wogonin, and morusin fall into class III (classification given by the Center for Drug Evaluation and Research) labeled as slightly toxic (Osterberg and See, 2003, Onawole et al., 2018) . All the selected flavonoids are predicted to be not readily biodegradable, but all are noncarcinogenic. All the chosen flavonoids showed no AMES toxicity. The structure-activity relationship of flavonoids is derived from the interactions with the main protease active site. The flavonoids of different categories have similar molecular weight values, several rotatable bonds, total polar surface area, and hence activity can be related to hydrogen bond acceptor and donor hydroxyl groups are essential in demonstrating the SARs of flavonoids, especially at positions 3 of ring C, 3' of ring B, and 7 of ring A (Figure 4) . Herein, we applied a library of flavonoids composed of nine different subcategories. Apegenin, wogonin and scutellarein, cirsilol, tangeritin, baiclain belongs to subclass flavones having the binding energies of -8.3, -8.2, -8.3, -7.4 and -6.4, -8.1 Kcal/mol, respectively. The lower the binding affinities of scutellarein, apegenin, baiclain, and wogonin with higher interactions were due to hydroxyl group at carbon seven. Scutellarein, baiclain, apegenin, and wogonin have four, three, three, and two hydroxyl groups, respectively. These groups participate in forming a conventional hydrogen bond with active residues (Thr24, His41, Cys44, Phe140, Thr26, Gly123, His164, Leu27, Glu166, His172, Leu141, Thr190, Asn142, Gln189, and Gln192). Hydrogen bonds have an essential role in stability and macromolecular recognition (Klaholz and Moras, 2002) . Hydrogen bonds also play a crucial role in several biological processes; hence, hydrogen bonds have a primary role in determining pharmacological receptors. The presence of bulky methyl groups in tangeritin and cirsilol reduces the interactions and increases the binding energy. The binding modes of some representative compounds are presented in Figure 5 . The binding pose of baiclain in the active site of the main protease and 3D interactions are shown in Figure 6 . According to docking results of prenylated flavonoids, the binding energies of kuwanon C, morusin, Artonin E were -8.1, -8.3 and -9.0 Kcal/mol. Their lower binding energies and higher efficient in silico are demonstrated due to the presence of five types of interactions. These interactions are the conventional hydrogen bond, pi-donor hydrogen bond, pi-sigma, pi-alkyl and alkyl type. Although Kuwanon C formed the five conventional hydrogen bonds with active residues, the binding energy of Artonin E was much lower than Kuwanon C. It may be due to the presence of other types of interaction that stabilized the ligand-receptor complex. Moreover, from docking results it was cleared that hydrophobic substituent may enhanced the bioactivity of the compounds. The binding pattern of prenylated flavonoids is represented in Figure 7 . Naringenin, Eriodictoyl, kaempferide, Hesperetin belong to the subcategory flavanones. The binding energy of no one compound was less than -8.0Kcal/mol, although they have a hydroxyl group at position seven. The binding affinity of Eriodictoyl was -7.4Kcal/mol, and it formed the conventional bonds with critical residues. Naringenin formed the six conventional hydrogen bonds with Thr190, Glu166, His163, and Leu141, having the binding affinity of -7.7 Kcal/mol. Naringin is a disaccharide derivative of naringenin that presented the binding affinity of -8.2Kcal/mol. It means that sugar residue by substituting hydroxyl group resulted in a decrease in binding energy with higher interactions. The binding modes of Hesperetin, Naringenin are depicted in Figure 8 . In bioflavonoids, namely Amentoflavone, robustaflavone presented the lowest binding affinity in the whole library of flavonoids. The binding affinity of Amentoflavone and pi-alkyl interactions with Met165, Met49, Cys145. Fisetin also has strong hydrophobic contact with Asp187, Gly143, Asn142, Ser144, His172, Glu166, His164, Gln189 Arg18 8. The predicted 2D interaction was depicted in Figure 10 . The docking results of isoflavones showed that Daidzein's two hydroxyl groups formed the four hydrogen bonds with Arg188, Leu141, Cys145 and Gly143 with the binding affinity of -7.5Kcal/mol Figure 10 . In Chalcones and Flavan-3-ols, the lower binding energy (-7.1Kcal/mol) was observed in Chalconaringenin has four free hydroxyl groups that formed hydrogen bonds with His163, Leu141, Gly143, and Arg188 ( Figure 10 ). While in Flavan-3-ols the lower binding energy (-7.6Kcal/mol) of Catechin may be due to presence of five free hydroxyl groups and in these four conventional hydrogen bonds was observed with Phe140, Glu166, and Gln189 key residues. The type and amount of bonding exhibited by the inhibitor in the protein's active site showed the compounds' affinity toward the drug. In flavonoids, the number and position of hydroxyl groups are essential in determining the inhibitory potential. In addition to conventional hydrogen bonds, other types of interactions like van der waals, pi-sigma, pi-sulfur, pi-alkyl, and carbon-hydrogen bonds are also important in accelerating the inhibitory potential of flavonoids. The AutoDcok vina score presented in Table1 has been further used to calculate the inhibition constant (Ki) using equation 1. The compounds with low inhibition constant values, usually in the micromolar range had a higher likelihood of pharmacological promiscuity (Onawole et al., 2018) . All the flavonoids with efficient binding affinity showed the inhibition constant within the micromolar range and qualified as lead compounds. The inhibition constant was calculated using equation 1. Other derivable ligand efficiency metrics were estimated using equation 2-5. The other metric ligand efficiency (LE) is used to calculate potency's efficiency per heavy atom in a ligand. The threshold value of ligand efficiency and fit quality (FQ) for drug-like candidates should be between 0.3 and 0.8, respectively (Murray et al., 2014 , Keserü and Makara, 2009 , Leeson and Springthorpe, 2007 , Reynolds et al., 2007 , Schultes et al., 2010 . The ligand efficiency in lipophilicity is known as ligand efficiency lipophilic price, denoted as LELP, is expected to be within the range of -10 to 10. All the selected flavonoids presented in The 100ns MD simulations of M pro complexed with Amentoflavone and Morusin were performed using Desmond. The fluctuation and stability of the receptor-ligand complexes during simulations in each case was analyzed and resulting trajectory was made with backbone root mean square deviation. It was observed that overall protein ligand was stable during entire simulation. Some fluctuations in protein were observed for first 10ns then it remained stable for next 90ns. Some ligand fluctuations were observed that might be due to hydrogen bond formation (figure 11 A). Root mean square fluctuation (RMSF) was studied to recognize the major fluctuations in the domain of protein. Peaks in the plot indicated that the major fluctuation in the loop region of protein was observed after 300 residues. This might be due to that alpha helices and beta strands were more rigid than unstructured (loop region) of the protein (figure 11B). A total of 27 ligand contacts were observed with residues of protein, from Thr24 to Leu50, Asn142 to Cys145, His163 to Thr169 and Val186 to Gln192. These residues formed four types of interactions including Hydrogen Bonds, Hydrophobic, Ionic and Water Bridges. It was observed that ligand achieved stability by forming polar and hydrophobic interaction with the active residues ( figure 11D ). The schematic details of ligand atoms two dimensional interactions with protein residues is presented in figure Figure 13 presents the MD trajectory analysis data for the morusin-Mpro complex. Ligand-protein complex remained almost stable throughout the simulation and the RMSD value ranged from 0.8 to 3.0Å for both morusin and main protease complex (figure13A). The flexibility of protein was also checked by evaluating the Root mean square fluctuations (RMSF) of individual residues of main protease. The entire main protease had an RMSF value between 0.5 to 1.5Å except for the terminal amino acid residues this might be due to that N and C terminal move rapidly than any other part of the protein (figure13B). The interaction between morusin and active residues of main protease is depicted in figure 13C and it presented the actual percentage of interaction. It was observed that Glu166 (72% and 46%), Gln189 (64%) and Gln192 (79% and 46%) interact through hydrogen bonding. The figure 13 D showed the fraction analysis of morusin with main protease during 100ns simulation period. The major active residues were found to interact through hydrogen bonding and water bridges with ligand. Finally, a timeline of the morusin-main protease interactions is plotted in the form of hydrogen-bonds, hydrophobic contact, ionic interaction and water bridges during 100ns of simulation period. The top panel presents the total contacts of morusin-main protease complex and the bottom panel indicates residue level interaction of the morusin. Overall, the morusin interacted well with the main protease binding pocket, especially, Glu166 Gln192 showed the efficient contact with ligand throughout the simulation period (figure14). The experimental results recently reported for some of the flavonoids against M pro have been presented in Table 5 which supported our findings through molecular docking and MD simulation studies. The experimental data also showed that bioactivity of the compounds against SARS MPro depended on the modifications in their molecular structures. Amentoflavone, a biflavone in proved to be more potent than apegenin (parent compound). The compounds with substituted groups exerted stronger inhibition than parent compound. Moreover prenylated flavonoids (Morusin, Kuwanon C) having the hydrophobic substituent were presented in literature with good inhibitory potential against influenza virus. Generally the prenylation of flavonoids generally enhanced the bioactivity and bioavailability of the compounds (Grienke et al., 2016) . The food sources for some common dietary flavonoids belonging to nine different subcategories are presented in Table 4 . These dietary sources are rich in flavonoids contents and may be used without any special precautions or trials for the prevention from the risk of COVID-19. The present study aims to identify the inhibitory activity of various flavonoids against SARS- Table 4 may serve as an essential preventive measure against SARS CoV-2 and helpful in adjuvant therapy of SARS CoV-2 patients. Induction of biologically active flavonoids in cell cultures of Morus nigra and testing their hypoglycemic efficacy Protease targeted COVID-19 drug discovery and its challenges: Insight into viral main protease (Mpro) and papain-like protease (PLpro) inhibitors Structure of coronavirus main proteinase reveals combination of a chymotrypsin fold with an extra α-helical domain The SARS-coronavirus papain-like protease: structure, function and inhibition by designed antiviral compounds Pharmacokinetics: the dynamics of drug absorption, distribution, metabolism, and elimination. Goodman and Gilman's the pharmacological basis of therapeutics Membrane position of ibuprofen agrees with suggested access path entrance to cytochrome P450 2C9 active site. The journal of physical chemistry A admetSAR: a comprehensive source and free tool for assessment of chemical ADMET properties Chloroquine and hydroxychloroquine as available weapons to fight COVID-19 Dietary quercetin and kaempferol: Bioavailability and potential cardiovascular-related bioactivity in humans Pymol: An open-source molecular graphics tool Duke's Phytochemical and Ethnobotanical Databases. Artemisia annua Identification of a novel coronavirus in patients with severe acute respiratory syndrome Which sources of flavonoids: complex diets or dietary supplements? Evaluation of blockbuster drugs under the rule-of-five Discovery of prenylated flavonoids with dual activity against influenza virus and Streptococcus pneumoniae Probing the α-helical structural stability of stapled p53 peptides: molecular dynamics simulations and analysis Flavonoids, a class of natural products of high pharmacological potency Content of potentially anticarcinogenic flavonoids of 28 vegetables and 9 fruits commonly consumed in the Netherlands Clinical features of patients infected with 2019 novel coronavirus in Wuhan, China. The lancet The continuing 2019-nCoV epidemic threat of novel coronaviruses to global health-The latest 2019 novel coronavirus outbreak in Wuhan Carbonic anhydrase and cholinesterase inhibitory activities of isolated flavonoids from Oxalis corniculata L. and their first-principles investigations Regulation of type I interferon responses Inhibition of SARS-CoV 3CL protease by flavonoids Comparison of simple potential functions for simulating liquid water The influence of lead discovery strategies on the properties of drug candidates Development of chemical inhibitors of the SARS coronavirus: viral helicase as a potential target Anatomy, micromorphology, and physiochemical analysis of Rhus succedanea var. himalaica root. Microscopy Research and Technique Critical evaluation of human oral bioavailability for pharmaceutical drugs by using various cheminformatics approaches C-H··· O Hydrogen Bonds in the Nuclear Receptor RARγ-a Potential Tool for Drug Selectivity Chemistry and biological activities of flavonoids: an overview. The scientific world journal The influence of drug-like concepts on decision-making in medicinal chemistry Coronavirus infections and immune responses Angiotensin-converting enzyme 2 is a functional receptor for the SARS coronavirus Drug Discovery Today: Technol Lead-and drug-like compounds: the rule-of-five revolution Experimental and computational approaches to estimate solubility and permeability in drug discovery and development settings Nosé-Hoover chains: The canonical ensemble via continuous dynamics Flavonoid (myricetin, quercetin, kaempferol, luteolin, and apigenin) content of edible tropical plants AutoDock4 and AutoDockTools4: Automated docking with selective receptor flexibility Validity of ligand efficiency metrics The anti-HIV drug nelfinavir mesylate (Viracept) is a potent inhibitor of cell fusion caused by the SARSCoV-2 spike (S) glycoprotein warranting further evaluation as an antiviral against COVID-19 infections The blood-brain barrier. Handbook of clinical neurology Structure based virtual screening of the Ebola virus trimeric glycoprotein using consensus scoring Toxicity of excipients-a food and drug administration perspective UCSF Chimera-a visualization system for exploratory research and analysis Flavonoids as antioxidants The role of molecular size in ligand efficiency Quantitative survey of narirutin, naringin, hesperidin, and neohesperidin in citrus The therapeutic potential of naringenin: a review of clinical trials Ligand efficiency as a guide in fragment hit selection and optimization Occurrence of flavonols in tomatoes and tomato-based products Phytoestrogen content of foods consumed in Canada, including isoflavones, lignans, and coumestan AutoDock Vina: improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading COVID-19, immune system response, hyperinflammation and repurposing antirheumatic drugs Molecular properties that influence the oral bioavailability of drug candidates Structure, function, and antigenicity of the SARS-CoV-2 spike glycoprotein Steroid hormone hydroxylase specificities of eleven cDNA-expressed human cytochrome P450s Preliminary exploration of the mechanism of Qingfei Paidu decoction agains novel coronavirus pneumonia based on network pharmacology and molecular docking technology Pathological findings of COVID-19 associated with acute respiratory distress syndrome. The Lancet respiratory medicine HIV protease inhibitor nelfinavir inhibits replication of SARS-associated coronavirus Chemical composition and pharmacological mechanism of Qingfei Paidu Decoction and Ma Xing Shi Gan Decoction against Coronavirus Disease 2019 (COVID-19): in silico and experimental study Isolation of a novel coronavirus from a man with pneumonia in Saudi Arabia Angiotensin-converting enzyme 2 (ACE2) as a SARS-CoV-2 receptor: molecular mechanisms and potential therapeutic target Urinary disposition of the soybean isoflavones daidzein, genistein and glycitein differs among humans with moderate fecal isoflavone degradation activity Epidemiology and cause of severe acute respiratory syndrome (SARS) in Guangdong, People's Republic of China A Review of Qing-Fei-Pai-Du-Tang and its Major Herbs Virus-encoded proteinases and proteolytic processing in the Nidovirales 2020) Amentoflavone -9 Submission of work requires that the piece to be reviewed has not been previously published. Upon acceptance, the Author assigns to the Arabian Journal of Chemistry (ARABJC) the right to publish and distribute the manuscript in part or in its entirety. The Author's name will always be included with the publication of the manuscript. The Author has the following nonexclusive rights: (1) to use the manuscript in the Author's teaching activities;(2) to publish the manuscript, or permit its publication, as part of any book the Author may write;(3) to include the manuscript in the Author's own personal or departmental (but not institutional) database or on-line site; and (4) to license reprints of the manuscript to third persons for educational photocopying. The Author also agrees to properly credit the Arabian Journal of Chemistry (ARABJC) as the original place of publication. The Author hereby grants the Arabian Journal of Chemistry (ARABJC) full and exclusive rights to the manuscript, all revisions, and the full copyright. The Arabian Journal of Chemistry (ARABJC) rights include but are not limited to the following:(1) to reproduce, publish, sell, and distribute copies of the manuscript, selections of the manuscript, and translations and other derivative works based upon the manuscript, in print, audio-visual, electronic, or by any and all media now or hereafter known or devised;(2) to license reprints of the manuscript to third persons for educational photocopying;(3) to license others to create abstracts of the manuscript and to index the manuscript; (4) to license secondary publishers to reproduce the manuscript in print, microform, or any computer-readable form, including electronic on-line databases; and (5) to license the manuscript for document delivery. These exclusive rights run the full term of the copyright, and all renewals and extensions thereof. I hereby accept the terms of the above Author Agreement. Author :-Date :-Editor in Chief:-AbdulRahman Al-Warthan Date:-10-06-2021 The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.