Pharmacognosy Magazine

: 2018  |  Volume : 14  |  Issue : 57  |  Page : 481--487

Molecular docking and density function theory studies of compounds from Euphorbia hirta and Bacopa monnieri to zika virus structural and nonstructural proteins

Sangeetha Kothandan1, Indu Purushothaman2, Stella Puthur George1, Sasikala Reddy Purushothaman1, Meena Karunakaran Sulochana3,  
1 Bioinformatics Infrastructure Facility Centre of DBT, Queen Mary's College (Autonomous), Chennai, Tamil Nadu, India
2 Department of Microbiology and Biotechnology, Presidency College (Autonomous), Chennai, Tamil Nadu, India
3 Bioinformatics Infrastructure Facility Centre of DBT, Queen Mary's College (Autonomous); Department of Chemistry, Queen Mary's College (Autonomous), Chennai, Tamil Nadu, India

Correspondence Address:
Meena Karunakaran Sulochana
Queen Mary's College (Autonomous), Chepauk, Chennai, Tamil Nadu


Background: Zika virus is an arbovirus belongs to the genus flavivirus and pose a serious global threat. The recent 2015 outbreak in Brazil was associated with a significant increase in microcephaly cases and other neurological complications in newborn babies and WHO declared Zika to be an international public health emergency. Currently, there is no specific treatment or Vaccine available for the Zika virus, and thus due to the unavailability of the antiviral drugs, the need for the identification of novel drugs is paramount. Materials and Methods: The compounds from two medicinal plants (Bacopa monnieri and Euphorbia hirta) were selected for the in silico molecular docking studies against the structural and nonstructural proteins of Zika virus. Quantum–chemical parameters density functional theory and absorption, distribution, metabolism, and excretion-toxicity (ADMET) was performed to identify the drug-likeliness properties. Results: Among the tested compounds, galloylquinic acid, Bacopaside III, and Bacopaside A were identified as leads against multiple targets of Zika virus. The identified compounds also exhibited desirable quantum chemical and ADMET properties. Conclusion: Hence, the compounds hampering the active site of the three different proteins playing a prime role in replication and fusion with desirable pharmacokinetic properties could be suggested for further in vitro and in vivo analysis of Zika virus. Abbreviations used: WHO: World Health Organization; DFT: Density function theory; ADMET: Absorption, distribution, metabolism, and excretion-toxicity; HBV: Hepatitis B virus; HCV: Hepatitis C virus; HIV: Human immunodeficiency virus; HSV: Herpes simplex virus; PDB: Protein data bank; OPLS: Optimized potential for liquid simulations; RMSD: Root-mean-square deviation; IFD XP: Induced fit docking extra precision; XP: Extra precision; PSA: Polar surface area; BBB: Blood brain barrier; A LogP 98: Atom-based Log P98; HOMO: Highest occupied molecular orbital; LUMO: Lowest unoccupied molecular orbital; DS V 4.0: Discovery Studio version 4.0; NS3: Nonstructural protein 3; NS1: Nonstructural protein 1

How to cite this article:
Kothandan S, Purushothaman I, George SP, Purushothaman SR, Sulochana MK. Molecular docking and density function theory studies of compounds from Euphorbia hirta and Bacopa monnieri to zika virus structural and nonstructural proteins.Phcog Mag 2018;14:481-487

How to cite this URL:
Kothandan S, Purushothaman I, George SP, Purushothaman SR, Sulochana MK. Molecular docking and density function theory studies of compounds from Euphorbia hirta and Bacopa monnieri to zika virus structural and nonstructural proteins. Phcog Mag [serial online] 2018 [cited 2022 Nov 30 ];14:481-487
Available from:

Full Text



About 32 compounds from two plants namely, Bacopa monneri and Euphorbia hirta were subjected for molecular docking against Non Structural proteins (NS1 and NS3) and the structural protein (Envelope protein Domain III) in Discovery studio V 4.0. Pharmacokinetic properties were predicted by ADMET. The compounds Galloylquinic acid, Bacopaside III and Bacopaside A were identified as leads against multiple targets and hence could be suggested for further in vitro and in vivo analysis.


Zika virus, an arbovirus belongs to the genus flavivirus poses a serious global threat. Zika virus was first isolated from a Macaca monkey in 1947 in Zika Forest near Entebbe, Uganda.[1] Zika virus generally causes mild disease with most common symptoms of fever, rash, joint pain, and conjunctivitis, however, microcephaly and Guillain-Barre syndrome are reported in the fetuses of the infected mother through prenatal transmission.[2],[3],[4],[5]

During the Zika outbreak in French Polynesia, an unusual increase in the number of neurological and autoimmune complications was identified. The French Polynesian outbreak spread to other Pacific islands, and autochthonous cases have been reported in New Caledonia (1400 confirmed cases), Cook Islands (932 suspected cases, 50 confirmed) Fiji, Samoa, and Solomon Islands.[6],[7],[8],[9],[10] The lack of antiviral drugs and the licensed vaccine to treat the disease necessitates the development of drugs to ZIKA virus.

Medicinal plants in the treatment of specific ailments have been in existence for several centuries. The novel scaffolds of the plants and their wide complex chemical constituents provide the source for the synthesis of new drugs. Antiviral studies of certain Indian plants specific to particular viruses, namely, Phyllanthus niruri (HBV), Glycyrrhiza glabra (HCV), Phyllanthus niruri, Aristolochia indica, Cassia occidentalis, Phyllanthus niruri, Withania somnifera, Tinospora cordifolia, Camellia sinenis, Calophyllum spp, Glycyrrhiza glabra, Phytolacca americana, trichosanthes, Kirilowii, Calophyllum, lanigerum (HIV), Sarracenia purpurea, Glycyrrhiza glabra, Rhus chinensis and Rhus javanica, Punica granatum (HSV), Phyllanthus emblica, and Sophora spp. (Coxsackie B virus) has been carried out for few viruses, however, this cannot be concluded as exhaustive.[11],[12],[13],[14],[15],[16],[17],[18],[19],[20],[21],[22],[23] Thus, it is constructive to identify the antiviral compounds against emerging Zika virus where the approved drug is unavailable. Hence, the compounds from two medicinal plants already in medicinal practice in India were selected for the in silico molecular docking studies against the structural and nonstructural proteins of Zika virus.

Bacopa monnieri is a perennial plant and possesses triterpenoid saponins called bacosides. It had proven antioxidant, hepatoprotective, and neuroprotective activity. Furthermore, earlier studies had demonstrated mechanism of action in acetylcholinesterase inhibition, choline acetyltransferase activation, β-amyloid reduction, increased cerebral blood flow, and monoamine potentiation.[24]

Euphorbia hirta belongs to the genus Euphorbiaceae and possess wide pharmacological activities such as antibacterial, antimalarial, anti-inflammatory, antiasthmatic, antidiarrheal, antioxidant, antifungal, and antiamoebic activity.[25] Thus, in this study, the compounds reported in the plants of B. monnieri and E. hirta were subjected for in silico molecular docking against the structural and nonstructural targets of Zika virus and furthermore absorption, distribution, metabolism, and excretion-toxicity (ADMET), density function theory (DFT) were undertaken to study the suitability of the compounds both in the biological and chemical standpoint.

 Materials and Methods


Structures of Bacopaside I, Bacopaside II, Bacopaside III, Bacopaside IV, Bacopaside V, bacoside, pseudojujubogenin reported in B. monnieri and the structures of β-Sitosterol, campesterol, stigmasterol, geranin, 1,2,3,4,6-penta-O-galloyl-β-D glucose, euphorbin-B, euphorbin–D, heptacosane, nonacosane, shikimic acid, choline, camphol, alpha-amyrin, beta-amyrin, campesterol, cycloartenol, euphorbol hexacosaote, friedelin, galloylquinicacid, geranin, leukocyanidin, quercitol, taxaxerol, taxaxerone, tetramethyl-2 hexadecen-ol from E. hirta were retrieved from PubChem, ChemSpider in Mol format. For each ligands, conformational search and optimization were carried out using CharmMM Force field.

Molecular docking

Receptor grid generation

The proteins for NSP3 (PDB ID: 5JRZ), NSP1 (PDB ID: 5K6K), and envelope protein Domain III (PDB ID: 5KVD) were retrieved from PDB and prepared. The prepared proteins were calculated for the receptor grids so that the ligands could bind within the predicted active site. The parameters of van der Waals scaling factor 1.00 and charge cutoff 0.25 were kept as default subjected to optimized potential for liquid simulations 2001 force field. A cubic box of specific dimensions was set around the centroid of the active site residues. The bounding box dimensions of 14 Å ×14 Å ×14 Å was placed for docking experiments.

Induced fit docking extra precision

Induced fit docking (IFD) extra precision (XP) was performed using the module-Induced Fit Docking of Schrodinger-Maestrov 9.1, LLC, New York. The entire receptor molecule constrained minimized with a root-mean-square deviation cutoff of 0.18 Å was selected for generation of centroid of the residues, and the box size was generated automatically. The initial Glide docking for each ligand was carried out. Side chains were trimmed automatically based on B-factor, with receptor and ligand van der Waals scaling of 0.70 and 0.50, respectively, and the number of poses generated were set to be 20. Prime side chain prediction and minimization were carried out in which residues were refined within 5.0 Å of ligand poses and side chains were optimized. This leads to a ligand structure and confirmation that is induced fit to each pose of the receptor structure. Finally, Glide XP redocking was carried out into structures within 30.0 kcal/mol of the best structure, and within the top 20 structures overall. The ligand was rigorously docked into the induced-fit receptor structure and obtained the XP score.[26],[27]

Absorption, distribution, metabolism, and excretion-toxicity

ADME and toxicity studies was performed in Discovery studio(DS) v4.0, Dassault systems, BIOVIA by considering the parameters such as Atom-based Log P98 (A LogP 98), polar surface area (PSA), blood-brain barrier (BBB), cytochrome P450, and hepatotoxicity.[28],[29],[30],[31],[32]

Density functional theory studies

DFT was carried out for the screened compounds to study the orbital energies (highest occupied molecular orbital [HOMO] and lowest unoccupied molecular orbital [LUMO]), energy gap and the dipole moment using DMOL3/B3 LYP in DS V 4.0.[33],[34],[35],[36]


Molecular docking

Among the compounds that were subjected for molecular docking using Schrodinger, 1,2,3,4,6-penta-O-galloyl-b-D-glucose showed highest Glide energy to the targets nonstructural protein 3 (NS3), nonstructural protein 1 (NS1), and E1 domain and XP bond followed by galloylquinic acid to NS3 and NS1. Likewise, compounds bacoside A and leukocyanidin showed higher binding energy and XP bond to envelope protein. Besides, the other compounds also exhibited equally higher binding energies, and the amino acid interaction profile of the compounds are represented in [Table 1] and [Table 2].{Table 1}{Table 2}

Interaction analysis of compound -1,2,3,4,6-penta-O-galloyl-b-D-glucose virus proteins

The compound 1,2,3,4,6-penta-O-galloyl-b-D-glucose interacted with a promising docking score of −11.32 and the Glide binding energy of −72.277 against the target NS3 of Zika virus protein. The backbone amino acids Arg 211, Glu 213 interacted by hydrogen bond interactions with the galloyl moiety of the compound. Side-chain hydrogen bond interactions occurred with the OH group of other galloyl moieties at Glu 187 and Glu 123. The Glide energy against the target envelope protein was −94.77 and Glide score of −13.24 and backbone hydrogen bond interaction was observed between the galloyl moiety and Lys 277. Furthermore, other six side chain interactions occurred between the OH group of galloyl moieties and the amino acids Glu 224, Lys 206, Asp 208, Ser 176, and Ser 228. Likewise, for the target NSP1, two hydroxyl groups and one oxygen group of galloyl moiety interacted with aminoacid residue Leu 430. Other side chain interactions were observed between hydroxyl groups and aminoacids Asp410, Glu 413, Glu 392, Asp 540, ASP 291, and the oxygen group with Arg 388. π–cation interaction was observed between the galloyl moiety and Arg 388 [Figure 1].{Figure 1}

Interaction Analysis of compound - galloylquinic acid with Zika virus proteins

Galloylquinic acid interacted with high energy score of −8.292 and Glide energy of −55.274 against NSP3. It extended 5 hydrogen bond interaction with NSP3 of Zika virus. Hydroxyl group of quinic acid moiety interacted both by backbone and side chain interaction with Pro 120, Glu 123, and Tyr 186, respectively. Galloylquinic acid interacted with envelope protein Domain III of ZIKA virus with Glide score of −7.972 and Glide energy of −42.958 and showed hydrogen bond interactions at three different sites. Hydroxyl groups of the quinic acid formed side chain hydrogen bond interactions with amino acids Ser 228 and Glu 178. Galloylquinic acid interacted with NS1 of Zika virus with low Glide energy of −39.026 and Glide score of −7.368. Three hydroxyl groups of quinic acid interacted with Thr 290 and GLN 413 and an oxygen group of gallic acid extended its interaction with HIE 288 [Figure 2].{Figure 2}

Interaction analysis of compound - leukocyanidin with Zika virus proteins

Leukocyanidin interacted with NS3 protein of Zika virus with Glide score of −6.593 and Glide energy of 42.12. In addition, it interacted with 5 hydrogen bond interactions. 2-hydroxyl groups of phenyl moiety interacted with Arg 211 through hydrogen bond side chain interaction whereas the two hydroxyl groups of chromene interacted by both backbone and side chain interaction with Glu 213 and Val 229, respectively. In addition, oxygen group of chromene interacted with Arg228.

Similarly, it exhibited higher binding energy and Glide score of −94.775 and −13.24, respectively, against envelope protein Domain III. The hydroxyl group of phenyl moiety interacted by hydrogen bond backbone with Leu 430. Furthermore, a pi-pi stacking was observed at HIS 486. Similarly, 3 hydroxyl groups of chromene interacted with the aminaocid Thr 290 by side chain interaction and with Asp 410 by backbone hydrogen bond interaction and another pi-pi stacking was observed with Phe at 299 [Figure 3].{Figure 3}

Interaction analysis of bacopasides with Zika virus proteins

Bacopaside III interacted with NSP3 of Zika virus with higher interaction energy of −36.625 and Glide score of −5.909. The hydroxyl and oxygen group of pseudojujubogenin moiety of the compound interacted with Asn 212 both by backbone and side chain interactions. Furthermore, the OH group of β-D-glucopyranosyl–(1-3)-α-L-arabinopyranosyl interacted with Val 229 [Figure 4].{Figure 4}

The hydroxyl group of jujubogenin of Bacopaside A interacted with Glu 213 of NSP3 and ASP 540 of NSP1 by side chain interactions. Likewise, the hydroxyl groups of sugar moieties interacted with Ser 208, Asn 201, Val 127 of NS3 and with ASP 410, Leu 430, Ser 452, and Thr 449 of NSP1 by both side chain and hydrogen bond interaction. Bacopaside IV interacted with both NS3 and NS1whereas Bacopaside V interacted with NS1 and envelope protein Domain III of ZIKA virus and the detailed interaction has been tabulated in [Table 1] and [Table 2].

Thus, in Bacopasides, both the jujobogenin and pseudojujubogenin moiety along with the sugar moieties extended interactions with the multitargets of ZIKA virus.

Interaction analysis of compound quercitol with Zika virus proteins

Quercitol exhibited Glide energy of −26.77 and Glide score of −5.316 against the target NS1 and Glide energy of −27.604 and Glide score of −5.904 to NS3, however, has not shown any interaction with the envelope protein Domain 3 of ZIKA virus. Specific hydrogen bond interactions were observed both at the backbone and side chain between the amino acid residues Arg 211, Tyr 186, Arg 228 of NS1, and OH group of quercitol. Furthermore, both the Arg 211 and 228 interactions were positively charged. Similar interactions were observed between Trp210, Thr 256, and HIE 253 of NSP3 and the hydroxyl group of quercitol.

Density function theory

DFT studies were carried out for the compounds that showed higher binding energies. The Energy gap (ΔE) represents the function of reactivity and thus lower separation energy signifies the higher reactivity of the compounds. The top-ranked or lead compounds 1,2,3,4,6-penta-O-galloyl-b-D-glucose, Bacopaside III, IV, A, galloylquinic acid showed very low separation energies in comparison to other molecules (0.1665, 0.099, 0.192,0.150, 0.156) whereas quercitol and leukocyanidin showed a modest increased energy (0.2981, 0.2022). Hence, the results obtained from molecular orbital energies were highly in association to binding energies of the compounds obtained during molecular docking. The parameters such as HOMO, LUMO, energy gap, dipole movement of the compounds are presented in [Table 3].{Table 3}

Absorption, distribution, metabolism, and excretion-toxicity

ADME properties such as PSA, Alog P98, absorption, aqueous solubility, BBB level, hepatotoxicity, and CYP2D6 were studied for the compounds. The intestinal absorption models are represented by 95% and 99% confidence ellipses in the ADMET_PSA_2D and ADMET_Alogp98 plane in which the upper limit are 131.62 and 148.12, respectively. In addition, the absorption of compounds is represented as good (0), moderate (1), poor (2), and very poor (3). Plasma surface area influences the drug transportation and permeability, and the lipophilicity is represented as the logarithm of the partition coefficient between n-octanol and water. Compounds with PSA <140 A °2 and A logp98 <5 showed optimum cell permeability. Thus, the compounds leukocyanidin, quercitol, shikimic acid showed the value of PSA (133.82, 104.07, and 100.56, respectively). Similarly, all the compounds showed good log P value. Leukocyanidin and shikimic acid showed increased absorption activity while the other compounds showed poor absorption. Admet aqueous solubility predicts the solubility of each compound in water at 25°C. The tested compounds were found to be hydrophilic except quercitol and shikimic acid exhibiting lipophilic properties hence suggesting the good bioavailability of the compounds.

Except 1,2,3,4,6-penta-O-galloyl-b-D-glucose, other compounds were non hepatotoxic and showed the value <1. The BBB penetration levels of the compounds were 4 without any violations and also non-inhibitors of CYP2D6 suggesting the possible good metabolization by CYP2D6. The tested compounds were both noncarcinogenic and nonmutagenic and exhibited moderate to no ocular irritancy except for the compounds Bacopaside IV, V, leukocyanidin, and shikimic acid that exhibited severe ocular irritancy. The ADME and toxicity properties of the compounds are represented in [Table 4].{Table 4}


Of the 32 compounds subjected for screening against the ZIKA virus structural and nonstructural protein targets, 10 compounds showed promising interaction with at least one of the targets or to all the three targets. The compounds 1,2,3,4,6-penta-O-galloyl-b-D-glucose, Bacopaside III, IV, V, bacoside A, geranin, galloylquinic acid, leukocyanidin, quercitol, and shikimic acid interacted with the ZIKA targets.

Furthermore, the quantum chemical parameters of the compounds were estimated based on the three components that the polar molecules dissociate better than nonpolar molecules which is expressed by dipole movement. Second, the Dipole movement is also an index of lipophilicity and ability of drug molecule to cross various biological membranes. Third, the smaller energy gap (ΔE) between HOMO and LUMO permits the transfer and exchange of electron which leads to an increase in the reactivity of the compounds.

Based on that the compounds 1,2,3,4,6-penta-O-galloyl-b-D-glucose, galloylquinic acid, shikimic acid, and leukocyanidin showed the lowest dipole movement (1.16491, 1.323, 1.385, 2.453) when compared to the other subjected compounds. This suggested that the inhibitors were found to be more hydrophobic (lipophilic) and hence could possibly show increased biological activities on further in vitro evaluation.

Although the compound 1,2,3,4,6-penta-O-Galloyl-β-D-glucose showed higher binding energy, lowest dipole moment, the ADMET analysis revealed it as hepatotoxic. Thus, based on the molecular docking, quantum chemical parameters and the ADMET properties, galloylquinic acid from E. hirta, Bacopaside III, Bacopaside A from B. monnieri were identified as potential lead molecules that could be tested for further in vitro analysis.


The study has identified the compound galloylquinic acid, Bacopaside III, and Bacopaside A as potential inhibitors of Zika virus protein targets. The screened compounds were effective against multiple targets (NS1, NS3, and envelope protein domain III of ZIKA virus) involved in genome replication, RNA synthesis and the fusion with the host cell. Hence, the compounds hampering the active site of the three different protein playing prime roles in replication with acceptable pharmacokinetic properties could be suggested for further in vitro and in vivo analysis.


This work was supported by Bioinformatics Infrastructure Facility Centre of Department of Biotechnology, Ministry of Science and Technology, Government of India vide Grant No. BT/BI/25/068/2012/2015.

Financial support and sponsorship

This work was supported by Bioinformatics Infrastructure Facility Centre of Department of Biotechnology, Ministry of Science and Technology, Govt. of India vide Grant No.BT/BI/25/068/2012/2015.

Conflicts of interest

There are no conflicts of interest.


1Kuno G, Chang GJ, Tsuchiya KR, Karabatsos N, Cropp CB. Phylogeny of the genus Flavivirus. J Virol 1998;72:73-83.
2Buathong R, Hermann L, Thaisomboonsuk B, Rutvisuttinunt W, Klungthong C, Chinnawirotpisan P, et al. Detection of Zika virus infection in Thailand, 2012-2014. Am J Trop Med Hyg 2015;93:380-3.
3Simpson DI. Zika virus infection in man. Trans R Soc Trop Med Hyg 1964;58:335-8.
4Oehler E, Watrin L, Larre P, Leparc-Goffart I, Lastere S, Valour F, et al. Zika virus infection complicated by Guillain-Barre syndrome – Case report, French Polynesia, December 2013. Euro Surveill 2014;19. pii: 20720.
5Carteaux G, Maquart M, Bedet A, Contou D, Brugières P, Fourati S, et al. Zika virus associated with meningoencephalitis. N Engl J Med 2016;374:1595-6.
6Dirlikov E, Major CG, Mayshack M, Medina N, Matos D, Ryff KR, et al. Guillain-Barré syndrome during ongoing Zika virus transmission-Puerto Rico, January 1-July 31, 2016. MMWR Morb Mortal Wkly Rep 2016;65:910-4.
7Dos Santos T, Rodriguez A, Almiron M, Sanhueza A, Ramon P, de Oliveira WK, et al. Zika virus and the Guillain-Barré syndrome – Case series from seven countries. N Engl J Med 2016;375:1598-601.
8Rozé B, Najioullah F, Fergé JL, Apetse K, Brouste Y, Cesaire R, et al. Zika virus detection in urine from patients with Guillain-Barré syndrome on Martinique, January 2016. Euro Surveill 2016;21:30154.
9Besnard M, Eyrolle-Guignot D, Guillemette-Artur P, Lastère S, Bost-Bezeaud F, Marcelis L, et al. Congenital cerebral malformations and dysfunction in fetuses and newborns following the 2013 to 2014 Zika virus epidemic in French Polynesia. Euro Surveill 2016;21:1560-7917.
10Cauchemez S, Besnard M, Bompard P, Dub T, Guillemette-Artur P, Eyrolle-Guignot D, et al. Association between Zika virus and microcephaly in French Polynesia, 2013-15: A retrospective study. Lancet 2016;387:2125-32.
11Venkateswaran PS, Millman I, Blumberg BS. Effects of an extract from Phyllanthus niruri on hepatitis B and woodchuck hepatitis viruses: In vitro and in vivo studies. Proc Natl Acad Sci U S A 1987;84:274-8.
12Thyagarajan SP, Subramanian S, Thirunalasundari T, Venkateswaran PS, Blumberg BS. Effect of Phyllanthus amarus on chronic carriers of hepatitis B virus. Lancet 1988;2:764-6.
13Wang M, Cheng H, Li Y, Meng L, Zhao G, Mai K, et al. Herbs of the genus Phyllanthus in the treatment of chronic hepatitis B: Observations with three preparations from different geographic sites. J Lab Clin Med 1995;126:350-2.
14van Rossum TG, Vulto AG, Hop WC, Schalm SW. Glycyrrhizin-induced reduction of ALT in European patients with chronic hepatitis C. Am J Gastroenterol 2001;96:2432-7.
15Da Nagao Y, Sata M, Suzuki H, Tanikawa K, Itoh K, Kameyama T, et al. Effectiveness of glycyrrhizin for oral lichen planus in patients with chronic HCV infection. J Gastroenterol 1996;31:691-5.
16Naik AD, Juvekar AR. Effects of alkaloidal extract of Phyllanthus niruri on HIV replication. Indian J Med Sci 2003;57:387-93.
17Pengsuparp T, Serit M, Hughes SH, Soejarto DD, Pezzuto JM. Specific inhibition of human immunodeficiency virus type 1 reverse transcriptase mediated by soulattrolide, a coumarin isolated from the latex of Calophyllum teysmannii. J Nat Prod 1996;59:839-42.
18Spino C, Dodier M, Sotheeswaran S. Anti-HIV coumarins from Calophyllum seed oil. Bioorg Med Chem Lett 1998;8:3475-8.
19Taylor PB, Culp JS, Debouck C, Johnson RK, Patil AD, Woolf DJ, et al. Kinetic and mutational analysis of human immunodeficiency virus type 1 reverse transcriptase inhibition by inophyllums, a novel class of non-nucleoside inhibitors. J Biol Chem 1994;269:6325-31.
20Kurokawa M, Hozumi T, Basnet P, Nakano M, Kadota S, Namba T, et al. Purification and characterization of eugeniin as an anti-herpesvirus compound from Geum japonicum and Syzygium aromaticum. J Pharmacol Exp Ther 1998;284:728-35.
21Jadhav P, Kapoor N, Thomas B, Lal H, Kshirsagar N. Antiviral potential of selected Indian medicinal (ayurvedic) plants against herpes simplex virus 1 and 2. N Am J Med Sci 2012;4:641-7.
22Liu Q, Wang YF, Chen RJ, Zhang MY, Wang YF, Yang CR, et al. Anti-coxsackie virus B3 norsesquiterpenoids from the roots of Phyllanthus emblica. J Nat Prod 2009;72:969-72.
23Mishra KP, Sharma N, Diwaker D, Ganju L, Singh SB. Plant derived antivirals: A potential source of drug development. Virol Antivir Res 2013;2:1-9.
24Aguiar S, Borowski T. Neuropharmacological review of the nootropic herb Bacopa monnieri. Rejuvenation Res 2013;16:313-26.
25Kim YH, Yang X, Yamashita S, Kumazoe M, Huang Y, Nakahara K, et al. 1,2,3,4,6-penta-O-galloyl-β-D-glucopyranose increases a population of T regulatory cells and inhibits IgE production in ovalbumin-sensitized mice. Int Immunopharmacol 2015;26:30-6.
26Halgren TA. Identifying and characterizing binding sites and assessing druggability. J Chem Inf Model 2009;49:377-89.
27Halgren TA, Murphy RB, Friesner RA, Beard HS, Frye LL, Pollard WT, et al. Glide: A new approach for rapid, accurate docking and scoring 2. Enrichment factors in database screening. J Med Chem 2004;47:1750-9.
28Egan WJ, Merz KM Jr., Baldwin JJ. Prediction of drug absorption using multivariate statistics. J Med Chem 2000;43:3867-77.
29Cheng A, Merz KM Jr. Prediction of aqueous solubility of a diverse set of compounds using quantitative structure-property relationships. J Med Chem 2003;46:3572-80.
30Susnow RG, Dixon SL. Use of robust classification techniques for the prediction of human cytochrome P450 2D6 inhibition. J Chem Inf Comput Sci 2003;43:1308-15.
31Ertl P, Rohde B, Selzer P. Fast calculation of molecular polar surface area as a sum of fragment-based contributions and its application to the prediction of drug transport properties. J Med Chem 2000;43:3714-7.
32van de Waterbeemd H, Smith DA, Jones BC. Lipophilicity in PK design: Methyl, ethyl, futile. J Comput Aided Mol Des 2001;15:273-86.
33Xia X, Maliski EG, Gallant P, Rogers D. Classification of kinase inhibitors using a Bayesian model. J Med Chem 2004;47:4463-70.
34Atlam FM, Awad MK, El-Bastawissy EA. Computational simulation of the effect of quantum chemical parameters on the molecular docking of HMG-CoA reductase drugs. J Mol Struct 2014;1075:311-26.
35Becke AD. Density-functional thermochemistry. III. The role of exact exchange. J Chem Phys 1993;98:5648
36Lee C, Yang W, Parr RG. Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density. Phys Rev B 1988;37:785.