key: cord-0722715-j3rw85i5 authors: Alesawy, Mohamed S.; Abdallah, Abdallah E.; Taghour, Mohammed S.; Elkaeed, Eslam B.; H. Eissa, Ibrahim; Metwaly, Ahmed M. title: In Silico Studies of Some Isoflavonoids as Potential Candidates against COVID-19 Targeting Human ACE2 (hACE2) and Viral Main Protease (M(pro)) date: 2021-05-10 journal: Molecules DOI: 10.3390/molecules26092806 sha: 6c8a7dabc79a4230107930614fd8f592e6bb6994 doc_id: 722715 cord_uid: j3rw85i5 The Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) caused the “COVID-19” disease that has been declared by WHO as a global emergency. The pandemic, which emerged in China and widespread all over the world, has no specific treatment till now. The reported antiviral activities of isoflavonoids encouraged us to find out its in silico anti-SARS-CoV-2 activity. In this work, molecular docking studies were carried out to investigate the interaction of fifty-nine isoflavonoids against hACE2 and viral M(pro). Several other in silico studies including physicochemical properties, ADMET and toxicity have been preceded. The results revealed that the examined isoflavonoids bound perfectly the hACE-2 with free binding energies ranging from −24.02 to −39.33 kcal mol(−1), compared to the co-crystallized ligand (−21.39 kcal mol(–1)). Furthermore, such compounds bound the M(pro) with unique binding modes showing free binding energies ranging from −32.19 to −50.79 kcal mol(–1), comparing to the co-crystallized ligand (binding energy = −62.84 kcal mol(–1)). Compounds 33 and 56 showed the most acceptable affinities against hACE2. Compounds 30 and 53 showed the best docking results against M(pro). In silico ADMET studies suggest that most compounds possess drug-likeness properties. In December 2019, an outbreak of severe pneumonia caused by the novel severe SARS-CoV-2 originated in Wuhan, China. The infection spread all over the world causing coronavirus disease [1, 2] . By October 2020, COVID-19 caused more than 33 million infections and more than 1 million deaths according to the WHO [3] . Unfortunately, till now there is no specific antiviral drug available for the treatment of COVID-19-infected people. However, some drugs such as remdesivir showed modest activity through decreasing the mortality rate and treatment time [4] . M pro is an essential non-structural chymotrypsin-like cysteine proteases enzyme for the replication of coronavirus. It works on two large polyproteins (PP1a and PP1ab) releasing 16 essential non-structural proteins (NSPs 1-16) [5, 6] . Angiotensin-converting enzyme (ACE-2) is a crucial enzyme in the renin-angiotensin system. It is a significant target for antihypertensive drugs [7] . It is primarily expressed in renal tubular epithelium and vascular endothelium cells [8] . It was also reported to be Drug-likeness properties were calculated using Lipinski's rule of five, which suggested that the absorption of an orally administered compound is more likely to be better if the molecule satisfies at least three of the following rules: (i) H bond donors (OH, NH, and SH) ≤5; (ii) H bond acceptors (N, O, and S atoms) ≤10; (iii) molecular weight <500; (iv) logP <5. Compounds violating more than one of these rules could not have good oral bi- Drug-likeness properties were calculated using Lipinski's rule of five, which suggested that the absorption of an orally administered compound is more likely to be better if the molecule satisfies at least three of the following rules: (i) H bond donors (OH, NH, and SH) ≤5; (ii) H bond acceptors (N, O, and S atoms) ≤10; (iii) molecular weight <500; (iv) logP <5. Compounds violating more than one of these rules could not have good oral bioavailability [38] . The pharmacokinetic properties (ADMET) of isoflavonoids and adherence with Lipinski's rule of five were calculated using Discovery studio 4.0 software(Accelrys software Inc., San Diego, CA, USA) [39] . The title molecules were investigated with the aid of docking studies using Discovery At first, water molecules were removed from the complex. Using the valence monitor method, the incorrect valence atoms were corrected. The energy minimization was then accomplished through the application of force fields CHARMM and MMFF94 [42] [43] [44] [45] . The binding sites were defined and prepared for docking processes. Structures of the tested isoflavonoids and the co-crystallized ligands were sketched using ChemBioDraw Ultra 14.0 (PerkinElmer, Waltham, MA, USA) [46] and saved as MDL-SD files. Next, the MDL-SD files were opened, 3D structures were protonated, and energy was minimized by implementing force fields CHARMM and MMFF94, then adjusted for docking. CDOCKER protocol was used for docking studies using CHARMM-based molecular dynamics (MD) to dock the co-crystallized ligands into a receptor binding site [47, 48] . In the docking studies, a total of 10 conformers were considered for each molecule. Finally, according to the minimum free energy of binding interaction, the most ideal pose was chosen. The toxicity parameters for the examined compounds were calculated using Discovery studio 4.0 software (Accelrys software Inc., San Diego, CA, USA). Simeprevir was used as a reference drug. At first, the CHARMM force field was applied then the compounds were prepared and minimized according to the preparation of small molecule protocol. Then different parameters were calculated from the toxicity prediction (extensible) protocol. In the present study, an in silico computational study of compounds (1-59) was performed to determine their physicochemical properties according to the directions of Lipinski's rule of five [38] (Table 1) . Moreover, all the tested isoflavonoids were predicted to have good absorption behavior better than that of remdesivir. Also, the solubility levels of most compounds were expected to be in the good range (Table 2 and Figure 2 ). a BBB level, blood brain barrier level, 0 = very high, 1 = high, 2 = medium, 3 = low, 4 = very low. b Absorption level, 0 = good, 1 = moderate, 2 = poor, 3 = very poor. c PBB, plasma protein binding, 0 means less than 90%, 1 means more than 90%, 2 means more than 95%. d solubility level, 0 = extremely low, 1 = very low, 2 = low, 3 = good, 4 = optimal. Validation of the docking procedures was achieved via re-docking of the co-crystallized ligands against the active pocket of hACE2 and M pro . The calculated RMSD values Coronavirus spike receptor-binding domain complexed with its receptor hACE-2 (PDB: 6LZG) used as a target for the docking studies of selected isoflavonoids. The results demonstrated that all isoflavonoids bound strongly to hACE-2 with binding energies bitter than that of the co-crystallized ligand (NAG). This indicated that the affinity of the tested isoflavonoids toward hACE-2 is higher than that of the co-crystallized ligand ( Table 3) . Moreover, almost all the tested isoflavonoids exhibited binding modes similar to that of NAG. Coronavirus spike receptor-binding domain complexed with its receptor hACE-2 (PDB: 6LZG) used as a target for the docking studies of selected isoflavonoids. The results demonstrated that all isoflavonoids bound strongly to hACE-2 with binding energies bitter than that of the co-crystallized ligand (NAG). This indicated that the affinity of the tested isoflavonoids toward hACE-2 is higher than that of the co-crystallized ligand (Table 3) . Moreover, almost all the tested isoflavonoids exhibited binding modes similar to that of NAG. Table 3 . Free binding energies of the selected isoflavonoids and the co-crystallized ligand (NAG) against hACE-2 and amino acid residues involved in H. bonds and hydrophobic interaction. Binding Energy (kcal mol −1 ) The binding pattern of co-crystallized ligand (NAG) demonstrated single hydrogen bonding interaction with Ser371 residue (Figure 4) . NAG showed binding energies of −21.39 kcal mol −1 . It was found that most of the tested isoflavonoids exhibited binding modes similar to the reference molecule. Compounds 1 (−30.90 kcal mol -1 ) and 8 (−27.41 kcal mol −1) demonstrated an additional hydrogen bond with Asn343 residue (Figures 5 and 6 ). This extra hydrogen bond may account for the relatively high binding affinity of both compounds. Furthermore, compounds 33 ( Figure 7 ) and 56 ( Figure 8 ) were found to have good binding energy values of −36.35 and −34.90 kcal mol −1 , respectively. Compound 33 formed a binding mode like that of the reference ligand as it formed one hydrogen bond with Ser371 and seven hydrophobic interactions with Phe374, Phe342, Ser371, Asn343, Cys336, Glu340, and Ser373. Interestingly, compound 56 formed two hydrogen bonds with Ser371 and Cys336 in addition to seven hydrophobic interactions with Phe374, Phe338, Ser371, Val367, Cys336, Leu368, and Ser373. Such results indicate the significance of the tested isoflavonoids as potential inhibitors for hACE-2. Consequently, such compounds may inhibit the entrance of coronavirus into human cells. The docking results of isoflavonoids into the active site of coronavirus M pro (PDB: 6LU7) were listed in Table 4 . The results showed that all tested isoflavonoids can bind to M pro with one or more hydrogen bonds. At the same time, the tested compounds bound to the receptor with free binding energies ranging from −32.19 to −50.79 kcal mol −1 , compared to the co-crystallized (binding energy = −62.84 kcal mol −1 ). These results revealed that the affinities of the presented isoflavonoids against M pro are lower than that of N3. Despite that, the binding energies are still considerable, and their binding modes are great which making these isoflavonoids seem to be biologically active ligands to some extent. Compound 30 formed a binding mode like that of the reference ligand as it formed three hydrogen bonds with Glu166, Tyr54, and Asp187. Furthermore, it formed eight hydrophobic interactions with His41, Gln189, His163, Met165, Tyr54, Asp187, Leu167, and Glu166. For compound 53, it formed two hydrogen bonds with Glu166, Phe140. Besides, it formed six hydrophobic interactions with Glu166, Gln189, Leu141, Met165, His172, and Phe140. The docking results of isoflavonoids into the active site of coronavirus M pro (PDB: 6LU7) were listed in Table 4 . The results showed that all tested isoflavonoids can bind to M pro with one or more hydrogen bonds. At the same time, the tested compounds bound to the receptor with free binding energies ranging from −32.19 to −50.79 kcal mol −1 , compared to the co-crystallized (binding energy = −62.84 kcal mol −1 ). Glu166. For compound 53, it formed two hydrogen bonds with Glu166, Phe140. Besides, it formed six hydrophobic interactions with Glu166, Gln189, Leu141, Met165, His172, and Phe140. Figure 9 . the co-crystallized ligand (N3) docked into M pro , forming four H. bonds with Gln189, Tyr 54, Asp 142, Asp187. Based on the binding affinities of the tested compounds against hACE-2, we can obtain valuable SAR. Generally, the tested compounds showed decreased affinity against hACE-2 in descending order of isoflavone derivatives (compounds 34, 33, 35 and 37) > isoflavane derivatives (compounds 50, 53, 57 and 59) > isoflavone derivatives (compounds 19, 20 and 23) > isoflavone derivatives (compounds 1-3) > isoflava-3-ene derivatives (compounds 6-8) > isoflavane derivatives (compounds 4 & 5) > pterocarpanes derivatives (compounds 9-11) (Figure 15 ). For isoflavone derivatives (compounds 34, 33, 35 and 37) , it was found that compound 34 incorporating 3-hydroxy-3-methylbutyl moiety at 6-position was more active that compound 33 incorporating 3-methylbut-2-en-1-yl moiety at the same position. The Based on the binding affinities of the tested compounds against hACE-2, we can obtain valuable SAR. Generally, the tested compounds showed decreased affinity against hACE-2 in descending order of isoflavone derivatives (compounds 34, 33, 35 and 37) > isoflavane derivatives (compounds 50, 53, 57 and 59) > isoflavone derivatives (compounds 19, 20 and 23) > isoflavone derivatives (compounds 1-3) > isoflava-3-ene derivatives (compounds 6-8) > isoflavane derivatives (compounds 4 & 5) > pterocarpanes derivatives (compounds 9-11) ( Figure 15 ). member with free OH group at position-1 of phenyl ring (compound 8) which was more potent than compound 7 incorporating a dioxolan ring. Observing binding affinities of isoflavane derivatives. It was found that compound 4 incorporating an additional methoxy group at 6 position of phenyl ring showed better binding affinity against hACE-2 than the unsubstituted derivative (compound 5). Such a result may be attributed to the electron donating effect of the methoxy group. Concerning the activity of different pterocarpan derivatives, it was noted that compound 11, which contained an additional tetrahydrofuran ring attached to the chromene ring, showed better binding affinity inside the hACE-2 than compounds 9 and 10, which contained free OH groups at the chromene ring. For isoflavone derivatives (compounds 34, 33, 35 and 37) , it was found that compound 34 incorporating 3-hydroxy-3-methylbutyl moiety at 6-position was more active that compound 33 incorporating 3-methylbut-2-en-1-yl moiety at the same position. The latter was more active than compound 35 incorporating 2-hydroxy-3-methylbut-3-en-1-yl moiety at the same position. Compound 37 incorporating 3-methylbut-2-en-1-yl moiety at 8-position was less active than the corresponding members. With regard to isoflavane derivatives (compounds 50, 53, 57 and 59), it was found that 3-methylbut-2-en-1-yl moiety is critical for binding affinity. compound 50 incorporating this moiety at 3-and 6-positions of 4-chromanone nucleus was more active than compound 53 incorporating this moiety at 5-position of phenyl ring. The latter was more active than compound 57 incorporating 3-methylbut-2-en-1-yl moiety at both 6-position position of 4-chromanone nucleus and 3-position of phenyl ring. Compound 59 incorporating 3-methylbut-2-en-1-yl moiety at both 8-position position of 4-chromanone nucleus and 3-position of phenyl ring was less active than the corresponding members. Regarding isoflavone derivatives (compounds 19, 20, and 23), it was found that the presence of 1,3-dioxole moiety can affect the affinity depending on it position. compound 19 incorporating 1,3-dioxole moiety at 3,4-positions of phenyl ring was more active than compound 20 incorporating this moiety at 7,8-position of 4H-chromen-4-one nucleus. The latter was more active than compound 23 incorporating 1,3-dioxole moiety at 4,5-position position of phenyl ring. Then, we investigated the effect of substitutions at isoflavone derivatives on the binding affinity. It was found that the substitutions at 5-position with hydroxyl (compound 1) and methoxy (compound 3) group, increase the binding of isoflavones against hACE-2, with an increased affinity of hydroxyl derivative. Regarding the effect of substitutions at isoflava-3-ene, it was found that the derivative with additional pyran ring (compound 6) was more active than the corresponding member with free OH group at position-1 of phenyl ring (compound 8) which was more potent than compound 7 incorporating a dioxolan ring . Observing binding affinities of isoflavane derivatives. It was found that compound 4 incorporating an additional methoxy group at 6 position of phenyl ring showed better binding affinity against hACE-2 than the unsubstituted derivative (compound 5). Such a result may be attributed to the electron donating effect of the methoxy group. Concerning the activity of different pterocarpan derivatives, it was noted that compound 11, which contained an additional tetrahydrofuran ring attached to the chromene ring, showed better binding affinity inside the hACE-2 than compounds 9 and 10, which contained free OH groups at the chromene ring. Toxicity prediction was carried out based on the validated and constructed models in Discovery studio 4.0 software [49, 50] as follows. (i) FDA rodent carcinogenicity test which computes the probability of a compound to be a carcinogen. (ii) Carcinogenic potency TD 50 which predicts the tumorigenic dose rate 50 (TD 50 ) of a drug in a rodent chronic exposure toxicity test of carcinogenic potency [51] . (iii) Rat maximum tolerated dose (MTD) [52, 53] . (iv) Rat oral LD 50 which predicts the rat oral acute median lethal dose (LD 50 ) of a chemical [54] . (v) Rat chronic LOAEL which predicts the rat chronic lowest observed adverse effect level (LOAEL) value [55, 56] . (vi) Ocular irritancy predicts whether a particular compound is likely to be an ocular irritant and how severe the irritation is in the Draize test [57] . (vii) Skin irritancy predicts whether a particular compound is likely to be a skin irritant and how severe it is in a rabbit skin irritancy test [57] . As shown in Table 5 , most compounds showed in silico low toxicity against the tested models. FDA rodent carcinogenicity model indicated that most of the tested compounds are non-carcinogens. Only compounds 6, 9, and 10 were predicted to be carcinogens so that, these compounds do not have the likeness to be used as drugs. COVID-19 agents. Generally, the tested compounds showed decreased affinity against hACE-2 in descending order of isoflavone derivatives (compounds 33-35 and 37) > isoflavane derivatives (compounds 50, 53, 57 and 59) > isoflavone derivatives (compounds 19, 20 and 23) > isoflavone derivatives (compounds 1-3) > isoflava-3-ene derivatives (compounds 6-8) > isoflavane derivatives (compounds 4 and 5) > pterocarpan derivatives (compounds 9-11) Finally, compounds 33 and 56 showed the most acceptable affinity against hACE2; compounds 30 and 53 showed the best docking results against M pro . In addition, these compounds showed good physicochemical and cytotoxicity profiles. Moreover, in silico investigation of physicochemical properties, ADMET and toxicity studies revealed good properties and general low toxicity. Consequently, this study strongly suggests in vitro and in vivo studies for the most active isoflavonoids against COVID-19. Coronaviruses -drug discovery and therapeutic options Genetic Recombination, and Pathogenesis of Coronaviruses WHO. WHO Coronavirus Disease (COVID-19) Dashboard. Available online Remdesivir-An Important First Step Crystal structure of SARS-CoV-2 main protease provides a basis for design of improved α-ketoamide inhibitors Drug targets for corona virus: A systematic review Angiotensin-I-converting enzyme and its relatives Preparation of fatty acid methyl esters for gas-liquid chromatography A novel coronavirus associated with severe acute respiratory syndrome Quantitative mRNA expression profiling of ACE 2, a novel homologue of angiotensin converting enzyme Structural basis for the recognition of SARS-CoV-2 by full-length human ACE2 Structural insights into coronavirus entry Molecular interaction and inhibition of SARS-CoV-2 binding to the ACE2 receptor Angiotensin-converting enzyme 2 (ACE2) as a SARS-CoV-2 receptor: Molecular mechanisms and potential therapeutic target ACE2 receptor expression and severe acute respiratory syndrome coronavirus infection depend on differentiation of human airway epithelia Black ginseng and its saponins: Preparation, phytochemistry and pharmacological effects Chemical Constituents of Stems and Leaves of Tagetespatula L. and Its Fingerprint Antimicrobial peptides from marine invertebrates: Challenges and perspectives in marine antimicrobial peptide discovery Comprehensive Virtual Screening of the Antiviral Potentialities of Marine Polycyclic Guanidine Alkaloids against SARS-CoV-2 (Covid-19) New α-Pyrone derivatives from the endophytic fungus Embellisia sp New antimalarial benzopyran derivatives from the endophytic fungus Alternaria phragmospora Comparative biological evaluation of four endophytic fungi isolated from nigella sativa seeds Induction of apoptosis in human cancer cells through extrinsic and intrinsic pathways by Balanites aegyptiaca furostanol saponins and saponin-coated silvernanoparticles New combination approaches to combat methicillin-resistant Staphylococcus aureus (MRSA) Two new antileishmanial diketopiperazine alkaloids from the endophytic fungus Trichosporum sp Antileukemic α-pyrone derivatives from the endophytic fungus Alternaria phragmospora Nigrosphaerin A a new isochromene derivative from the endophytic fungus Nigrospora sphaerica In Vitro and In Silico Cytotoxic and Antibacterial Activities of a Diterpene from Cousinia alata Schrenk Biological evaluation and molecular docking study of metabolites from Salvadora Persica L. Growing in Egypt Biological and chemical evaluation of some African plants belonging to Kalanchoe species: Antitrypanosomal, cytotoxic, antitopoisomerase I activities and chemical profiling using ultra-performance liquid chromatography/quadrupoletime-of-flight mass spectrometer Antifungal prenylated isoflavonoids from Maclura aurantiaca The Flavonoids: Advances in Research since 1980 Antiviral isoflavonoid sulfate and steroidal glycosides from the fruits of Solanum torvum Soybean Saponin and Isoflavonoids. In Food Phytochemicals for Cancer Prevention I Induction of a 5-lipoxygenase product by daidzein is involved in the regulation of influenza virus replication Antiviral activity of substituted homoisoflavonoids on enteroviruses Synthesis and anti-picornavirus activity of homo-isoflavonoids Experimental and computational approaches to estimate solubility and permeability in drug discovery and development settings Design, synthesis, molecular modeling and anti-hyperglycemic evaluation of novel quinoxaline derivatives as potential PPARγ and SUR agonists Design, synthesis, molecular modeling and biological evaluation of novel 2, 3-dihydrophthalazine-1, 4-dione derivatives as potential anticonvulsant agents Design, synthesis, molecular modeling and anti-hyperglycemic evaluation of quinazolin-4 (3H)-one derivatives as potential PPARγ and SUR agonists Discovery and antiproliferative evaluation of new quinoxalines as potential DNA intercalators and topoisomerase II inhibitors ChemBioDraw Ultra 14.0; PerkinElmer Synthesis, docking, QSAR, ADMET and antimicrobial evaluation of new quinoline-3-carbonitrile derivatives as potential DNA-gyrase inhibitors Design, synthesis, molecular modeling and anti-hyperglycemic evaluation of phthalimide-sulfonylurea hybrids as PPARγ and SUR agonists Classification of kinase inhibitors using a Bayesian model Structure-Activity Relationships for Carcinogenic Potential Comparison of the dependence of the TD50 on maximum tolerated dose for mutagens and nonmutagens Correlation Between Carcinogenic Potency and the Maximum Tolerated Dose: Implications for Risk Assessment Comparison of in silico tools for evaluating rat oral acute toxicity In silico models for repeated-dose toxicity (RDT): Prediction of the no observed adverse effect level (NOAEL) and lowest observed adverse effect level (LOAEL) for drugs Assessment of the oral rat chronic lowest observed adverse effect level model in TOPKAT, a QSAR software package for toxicity prediction The Draize eye test The authors declare no conflict of interest.Sample Availability: Samples of the compounds are not available from the authors. Potency TD 50 (Rat) a For the carcinogenic potency TD 50 rat model, the examined compounds showed TD 50 values ranging from 0.44 to 322.42 mg Kg −1 body weight/day which are higher than simeprevir (0.280 mg Kg −1 body weight/day).Regarding the rat MTD model, the compounds showed MTD with a range of 0.061 to 0.764 g kg −1 body weight higher than simeprevir (0.003 g kg −1 body weight).Concerning the rat oral LD 50 model, the tested compounds showed oral LD 50 values ranging from 0.10 to 4.66 mg Kg −1 body weight/day), while simeprevir exhibited an oral LD 50 value of 0.21 mg Kg -1 body weight/day. About the rat chronic LOAEL model, the compounds showed LOAEL values ranging from 0.004 to 0.865 g kg −1 body weight. These values are higher than simeprevir (0.002 g kg −1 body weight). Moreover, most of the compounds were predicted to be irritant against the ocular irritancy model. On the other hand, the tested compounds were predicted to be mild or non-irritant against the skin irritancy model. There is an urgent global need to find a cure for COVID-19. The present work is an attempt to find some natural compounds with potential activity against COVID-19. Accordingly, docking studies were carried out for fifty-nine isoflavonoid derivatives against two essential targets (hACE-2 and M pro ). The obtained results showed that the tested isoflavonoids can strongly bind the hACE-2 and M pro with great binding modes. Based on in silico studies, SARs were established. SAR studies afforded an insight into the pharmacophoric groups which may serve as a guide for the design of new potential anti-