Computational repositioning of ethno medicine elucidated gB-gH-gL complex as novel anti herpes drug target
© Basha et al.; licensee BioMed Central Ltd. 2013
Received: 15 February 2012
Accepted: 11 March 2013
Published: 15 April 2013
Herpes viruses are important human pathogens that can cause mild to severe lifelong infections with high morbidity. They remain latent in the host cells and can cause recurrent infections that might prove fatal. These viruses are known to infect the host cells by causing the fusion of viral and host cell membrane proteins. Fusion is achieved with the help of conserved fusion machinery components, glycoproteins gB, heterodimer gH-gL complex along with other non-conserved components. Whereas, another important glycoprotein gD without which viral entry to the cell is not possible, acts as a co-activator for the gB-gH-gL complex formation. Thus, this complex formation interface is the most promising drug target for the development of novel anti-herpes drug candidates. In the present study, we propose a model for binding of gH-gL to gB glycoprotein leading from pre to post conformational changes during gB-gH-gL complex formation and reported the key residues involved in this binding activity along with possible binding site locations. To validate the drug targetability of our proposed binding site, we have repositioned some of the most promising in vitro, in vivo validated anti-herpes molecules onto the proposed binding site of gH-gL complex in a computational approach.
Hex 6.3 standalone software was used for protein-protein docking studies. Arguslab 4.0.1 and Accelrys® Discovery Studio 3.1 Visualizer softwares were used for semi-flexible docking studies and visualizing the interactions respectively. Protein receptors and ethno compounds were retrieved from Protein Data Bank (PDB) and Pubchem databases respectively. Lipinski’s Filter, Osiris Property Explorer and Lazar online servers were used to check the pharmaceutical fidelity of the drug candidates.
Through protein-protein docking studies, it was identified that the amino acid residues VAL342, GLU347, SER349, TYR355, SER388, ASN395, HIS398 and ALA387 of gH-gL complex play an active role in its binding activity with gB. Semi flexible docking analysis of the most promising in vitro, in vivo validated anti-herpes molecules targeting the above mentioned key residues of gH-gL complex showed that all the analyzed ethno medicinal compounds have successfully docked into the proposed binding site of gH-gL glycoprotein with binding energy range between -10.4 to -6.4 K.cal./mol.
Successful repositioning of the analyzed compounds onto the proposed binding site confirms the drug targetability of gH-gL complex. Based on the free binding energy and pharmacological properties, we propose (3-chloro phenyl) methyl-3,4,5 trihydroxybenzoate as worth a small ethno medicinal lead molecule for further development as potent anti-herpes drug candidate targeting gB-gH-gL complex formation interface.
KeywordsHerpes simplex virus gB-gH-gL complex Ethnomedicine Docking (3-chloro phenyl) methyl-3, 4, 5 trihydroxybenzoate
Role of ethnomedicine in finding novel antiherpes drugs
The nucleoside analog acyclovir, a widely used anti-herpes drug, is unable to cure the chronic infections and recurrence. Use of this drug often leads to the development of viral resistance coupled with side effects and viral latency . Therefore discovery of novel antiherpes drugs with high specificity to inhibit the viral infection deserves great effort. On the other hand, a wide array of antimicrobial compounds with antiviral activity is described in traditional medical literature. Traditional medicine enjoys new found interest among microbiologists because of increasing antibiotic resistance and limited life span of conventional antimicrobials. It is a common practice in various fields of Biomedicine to choose plants described in traditional medicine for finding new successful therapeutic agents . Most of the ethno medicinal compounds were ignored by medical practitioners despite of their proven in vitro and in vivo antiviral properties due to lack of evidence on their mode of action, effective dosage, chemical composition and toxicity .
Recent studies have reported promising anti Herpes viral activities of several plant extracts, some of the active constituents of this plant extracts are reported to have potential to interfere with RNA and protein synthesis, DNA replication, cellular fusion, attachment and penetration, virus entry and target cell binding, whereas some are thought to interfere in more than one step of herpes virus life cycle, resulting in complementary mechanisms of action to the existing antiviral drugs, as reviewed elsewhere in the literature . In this scenario, this present work is an attempt to elucidate the possible mode of action of various promising in vitro, in vivo validated anti herpes ethno medicinal compounds to inhibit the viral fusion mechanism by attenuating the gB-gH-gL complex formation interface in a computational approach.
Preparation of protein receptor
The crystal structures of glycoprotein B (gB) (PDB ID: 2GUM)  and gH-gL complex (PDB ID: 3M1C)  of Herpes simplex virus, which were resolved using X-Ray diffraction method with a resolution factor of 2.10 Å and 3.00 Å respectively were obtained from Protein Data Bank (PDB) . Methodology adapted for the preparation of protein receptor gH-gL complex for the semi flexible docking studies was described elsewhere in detail . The modified structure so obtained was saved in .pdb format and used for all docking studies.
Selection and preparation of ligands
Structures of the ethno medicinal compounds with proven in vitro and in vivo activity against herpes as reviewed elsewhere in literature  were retrieved from Pubchem database . All the compounds retrieved were screened with ADME constraints according to Lipinski’s rule. Finally, 29 compounds were selected for further studies, which were in accordance with Lipinski’s rule. Preparation of the selected ligands was done using Arguslab 4.0.1 software , addition of missing hydrogen atoms and fidelity of all bonds was checked using “add hydrogens” and “Clean Hybridization” options respectively. Geometry optimization was done using UFF Molecular Mechanics (MM) method [16–20]. Finally, all the compounds were saved in .mol format for further docking studies.
Protein-protein docking was performed using Hex 6.3 software  between Domain II of glycoprotein gB (receptor) and H2 Domain of gH-gL heterodimer complex (ligand). Out of 2000 possible solutions, the best solution binding mode is selected to identify the key residues involved in their binding activity. Following parameters are used in Docking controls of Hex 6.3 to get the best results: Correlation type-Shape only; FFT mode-3D; Post Processing-MM minimization; Grid dimension-0.4; Solutions-2000; Receptor range-180; Step size-7.5; Ligand range-180; Step size-7.5; Twist range-360; Step size-5.5; Distance range-40; Scan step-0.8; Substeps-0; Steric scan-18; Final search-25, whereas the parameters used in clustering controls to get the best results are as follows: Max Clusters-500; Sort solutions by-Cluster; Display clusters-Best; Cluster window-200; RMS threshold-3.0; Bumps threshold-0.
Docking between receptor gH-gL protein and ligand was performed using “Dock a Ligand” option of Arguslab 4.0.1 software. A spacing of 0.4 Å between the grid points was used. “ArgusDock” was selected as docking engine. Search was performed by enabling “High precision” option in Docking precision menu, “Dock” was chosen as the calculation type, “Flexible” for the ligand and “AScore” was used as the scoring function. A maximum of 150 poses were allowed to be analyzed, binding site box was set to 25 × 25 × 25 angstroms to encompass the entire active site, and each docking run was repeated three times to get the best results.
Analysis and visualization of docking simulation results
The docking results of Arguslab 4.0.1 were saved in .pdb format. Out of 150-model results, one with the least free energy bonding data was selected. Analysis, visualization of the interactions and rendering of all the docking images was done using Accelrys® Discovery Studio 3.1 Visualizer software .
Pharmacological properties of the compounds
Lipinski’s Filter (http://www.scfbio-iitd.res.in/utility/LipinskiFilters.jsp), Osiris Property Explorer (http://www.organic-chemistry.org/prog/peo/) and Lazar (http://lazar.in-silico.de/predict) online servers were used to check the pharmaceutical fidelity of the drug candidates.
Molecular descriptors, such as Log P, the number of hydrogen bond donors, the number of hydrogen bond acceptors, and the molecular mass of the compounds were analyzed by utilizing Lipinski Filter. Osiris Property Explorer and Lazar were used in analyzing various attributes of the drugs, such as toxicity and drug score.
Results and discussion
Active site identification and validation
Interaction of the ethno compounds with active site
The binding free energy docking simulation result of ethno compounds toward gH-gL complex of HSV and the number of hydrogen bonds formed along with the amino acid residues involved in the interactions
Name of the compound
Binding energy in K.cal/mol
No. of H-Bonds formed with receptor
Amino acid residues involved in H-bonding
ARG 318, TYR 319, SER 349 and ALA 392
TYR 319 and ASN 395
PRO 346 and ASN 395
ASN 395, PRO 346 and PRO 186
ASN 395 and TYR 319
ARG 318, HIS 398 and THR 251
THR 251, ARG 318, TYR 319, ASN 395 and HIS 398
ARG 318, GLY 270 and THR 350
ARG 318, TYR 319 and ASN 395
GLY 268 and VAL 394
ARG 318, ASN 395 and ALA 392
(3-chloro phenyl) methyl-3,4,5 trihydroxybenzoate
ARG 318, TYR 319, ASP 187 and ASN 395
THR 251, GLY 270 and HIS 398
ASN 395 and SER 349
GLY 268 and HIS 398
ARG 318, ASN 395, PRO 346, SER 349 and ASP 187
ARG 318, GLY 270 and ALA 392
GLY 268 and HIS 398
ARG 318 and ASP 187
GLY 270, ARG 318 and PRO 346
ARG 318, TYR 319 and ASN 395
ARG 318, TYR 319 and HIS 398
Binding free energy (ΔG binding)
The argusdock values of ΔG binding given in Table 1 shows that every ligand has negative ΔG, which shows that all the ethno medicinal compounds analyzed for binding with gH-gL complex are highly stable. These negative binding free energies also evidenced that the reaction is spontaneous.
Prediction of pharmacological properties
Lipinski Filter and Osiris Property Explorer were utilized to screen the drug candidate based on Lipinski’s Rule of Five and Oral Bioavailability. The pharmacological attributes prediction results are displayed in (Additional file 1: Table S1) and the structures of all the ethno medicinal compounds used for the study are shown in Additional file 1: Figure S1.
Based on the experimental values, it was inferred that all the ethno medicinal compounds successfully satisfied all the parameters of Lipinski’s Rule of Five with the exception of the hydrogen donor for Caffeoylquinic acid. The parameters of the Lipinski’s rule are as follows: the molecular weight must be < 500 Da, Log P < 5, the number of hydrogen donors must be < 5, the number of acceptor hydrogens must be < 10, and the refractivity molar range must be between 40–130.
As per the veber’s rule, oral bioavailability of drugs could be measured by the molecular weight, number of rotatable bonds (n rotb), number of hydrogen bonds, and the expanse of the drug’s polar surface (TPSA). The oral bioavailability was marked by small molecular weight (less than 500 Da); also, the number of rotatable bond must be less than 10, the number of hydrogen bond donors and acceptors must be less than 12, and TPSA values less than 140. Additional file 1: Table S1 shows that all the compounds have a good oral bioavaibility, with the exception of Caffeoylquinic acid and Rosmarinic acid.
The hydrophobicity of drugs could be inferred from Log P value. Log P Values are directly proportional to the oral hydrophobicity of the drug. More hydrophobic the drug is, higher is the ability of the drug to circulate longer in our body. It wouldn’t be easy to excrete such a drug. In the present investigation, the Log P values of the drug molecules were observed to be in the range of -0.4 to 6. The highest Log P value of 6 was observed for Bisbenzylisoquinoline, followed by 3.8 for Lignan and Yatein, which shows that Bis-benzylisoquinoline was the most hydrophobic molecule amongst all the analyzed. However, Anthraquinone, Galangin, Isoborneol, Lignan, Linalool, (3-chloro phenyl) methyl-3,4,5 trihydroxybenzoate, Oxyresveratrol, Physcion, Resveratrol, Rhein, Rosmarinic acid, Wogonin and Yatein also showed good hydrophobicity with a range of log P between 2–5.
The drug score value is one such useful practical value which is a combination of variables such as drug likeness, Log P, solubility, molecular weight and toxicity risk, which could be used for the evaluation of a potential drug candidate. A drug molecule with a better drug score has a greater chance to be a drug candidate. As per the values shown in Additional file 1: Table S1, Catechin with a drug score 0.87 showed the highest possibility of being a potential drug and Zeatin with a drug score of 0.1 shown to be the least possible potential drug.
Prediction of toxicity
Osiris Property Explorer and Lazar were used to predict toxicity of the compounds in the present study. The predictions were based on the functional group similarity for the query molecule with the in vitro and in vivo validated compounds present in the database of these two online programs.
Osiris Property Explorer results were visualized using color codes. The result of toxicity analysis of all the analyzed ethno compounds is shown in Additional file 1: Table S2. Green color shows low toxic tendency, yellow shows the midcore and red color shows high tendency of toxicity.
Additional file 1: Table S2 shows that Anthraquinone, Apigenin, Caffeic acid, Galangin, Kaempferol, Linalool, Morin, Phloroglucinol, Quercetin, Resveratrol, Zeatin have high toxicity. Whereas Allicin, Bis-benzylisoquinoline, Caffeoylquinic acid, Catechin, Cochinolide, Epiafzelechin, Isoborneol, Lignan, (3-chloro phenyl) methyl-3,4,5 trihydroxybenzoate, Methylgallate, Oxyresveratrol, Rhein, Rosmarinic acid, Wogonin and Yatein had no toxic effect, whereas other compounds had midcore toxicity.
Lazar is a software package which is used to detect mutagenic and/or carcinogenic properties based on the similarities in functional group with mutagenic and/or carcinogenic ones present in its database. Additional file 1: Table S3 shows the results of compounds toxicity analysis by lazar. The training dataset of molecules present in Lazar program were verified by conducting an assay with Salmonella typhimurium.
Additional file 1: Table S3 shows that Aloe-emdoin, Anthraquinone, Baicalein, Bisbenzylisoquinoline, Galangin, Morin, Physcion, Quercetin, Rhein, Wogonin and Zeatin have mutagenic properties. Whereas Caffeic acid, Caffeoylquinic acid, Catechin, Cochinolide, Epiafzelechin, Lignan, Linalool, (3-chloro phenyl) methyl-3,4,5 trihydroxybenzoate, Methylgallate, Oxyresveratrol, Phloroglucinol, Resveratrol, Rosmarinic acid and Yatein have no mutagenic properties. Moreover, the carcinogenicity of compounds was verified by in silico animal testing studies with Mouse. Additional file 1: Table S3 shows that all the compounds analyzed have no carcinogenicity, with the exception of Zeatin.
Screening for the best compound based on docking and drug likeliness results
It was found that (3-chloro phenyl) methyl-3,4,5 trihydroxybenzoate was the best anti herpes viral ethno medicinal compound based on its ΔG binding compared to other compounds followed by Bis-benzylisoquinoline and Rosmarinic acid with a binding energy of around -10 K.cal/mol. However, Bis-benzylisoquinoline and Rosmarinic acid have bad pharmacological and toxicity attributes, based on pharmacological properties. Rosmarinic acid has a TPSA value more in excess of threshold Veber’s rule parameter, whereas Bisbenzylisoquinoline showed considerable mutagenicity. (3-chloro phenyl) methyl-3,4,5 trihydroxybenzoate showed good pharmacological attributes, since it satisfied the Lipinski’s Rule, Veber’s Rule, Log P values, toxicity and drug score threshold.
However, based on Pharmacological properties, Caffeoylquinic acid, Catechin, Cochinolide, Epiafzelechin, Lignan, (3-chloro phenyl) methyl-3,4,5 trihydroxybenzoate, Methylgallate, Oxyresveratrol, Rosmarinic acid and Yatein showed good pharmacological attributes. These compounds were found to comply with Lipinski’s rule, Veber’s rule and oral bioavailability parameters with the exception of Caffeoylquinic acid and Rosmarinic acid, because of their TPSA value exceeded the threshold of 140.
The binding energy analysis showed that (3-chloro phenyl) methyl-3,4,5 trihydroxybenzoate had a least binding energy value of -10.4 K.cal/mol compared to other compounds followed by Bisbenzylisoquinoline and Rosmarinic acid with binding energy values of -10.1 K.cal/mol each, among all the compounds analyzed. Visualization of the interactions showed that O atom on C = O and -OH groups of all tested compounds played a major role in hydrogen bonds formation with amino acid residues of receptor molecule. Conclusively, our results strongly suggest that (3-chloro phenyl) methyl-3, 4, 5 trihydroxybenzoate is worth a lead ethno compound for developing further as a potent inhibitor against Herpes simplex virus targeting its fusion protein complex. On the other hand, successful repositioning of the analyzed compounds onto the proposed binding site of gH-gL complex shed light on its drug targetability for further investigations. However, as this present work is only a step forward towards understanding the mechanistic insights of a potent inhibitor for herpes viral entry into the host cell targeting gB-gH-gL complex interface, further in vitro and in vivo validations are required.
SHB specially thank and dedicate this work to his parents Syed Abdul Khader Basha and Syed Khursheed Begum for providing the required facilities to carry out this whole work at home and for their constant encouragement and support, without which this research work wouldn’t have been possible. The authors wish to thank Dr. N. Ramesh, Principal, REVA ISM for his much needed help. Dr. Atuluri Venu, Principal, Montessori Siva sivani Institute of Science and Technology-College of Pharmacy for extending his help and motivation. Dr.Vidyavathi N, Prof. Dept. of Biotechnology, REVA ISM for reading the earlier version of this manuscript. We gratefully acknowledge Doina Atanasiu and Ashoke Sharon for generously donating their valuable time in providing peer review for this manuscript.
- Prevention and control of Herpes virus diseases. Clinical and laboratory diagnosis and chemotherapy. Volume 63. 1985, Bulletin of the WHO, 182-185.Google Scholar
- Whitley RJ: Herpes simplex viruses. Fields Virology. Edited by: Knipe DM, Howley PM. 2001, Philadelphia: Lippincott, Williams and Wilkins, 2461-2509.Google Scholar
- Wild K, Bohner T, Folkers G, Schulz GE: The structures of thymidine kinase from Herpes simplex virus type 1 in complex with substrates and a substrate analogue. Protein Sci. 1997, 6: 2097-2106.View ArticlePubMedPubMed CentralGoogle Scholar
- Heldwein EE, Lou H, Bender FC, Cohen GH, Eisenberg RJ, Harrison SC: Crystal structure of glycoprotein B from Herpes simplex virus-1. Science. 2006, 313: 217-220. 10.1126/science.1126548.View ArticlePubMedGoogle Scholar
- Chowdary TK, Cairns TM, Atanasiu D, Cohen GH, Eisenberg RJ, Heldwein EE: Crystal structure of the conserved herpesvirus fusion regulator complex gH-gL. Nat Struct Mol Biol. 2010, 17 (7): 882-888. 10.1038/nsmb.1837.View ArticlePubMedPubMed CentralGoogle Scholar
- Krummenacher C, Supekar VM, Whitbeck JC, Lazear E, Connolly SA, Eisenberg RJ, Cohen GH, Wiley DC, Carfí A: Structure of unliganded HSV gD reveals a mechanism for receptor-mediated activation of virus entry. EMBO J. 2005, 24: 4144-4153. 10.1038/sj.emboj.7600875.View ArticlePubMedPubMed CentralGoogle Scholar
- Eisenberg RJ, Atanasiu D, Cairns TM, Gallagher JR, Krummenacher C, Cohen GH: Herpes Virus Fusion and Entry: A Story with Many Characters. Review, Viruses. 2012, 4: 800-832. 10.3390/v4050800.View ArticleGoogle Scholar
- Chatis PA, Crumpacker CS: Resistance of herpesviruses to antiviral drugs. Antimicrob Agents Chemother. 1992, 36 (8): 1589-1595. 10.1128/AAC.36.8.1589.View ArticlePubMedPubMed CentralGoogle Scholar
- Cox PA, Balick M: The ethnobotanical approach to drug discovery. Sci Am. 1994, 270: 82-87.View ArticlePubMedGoogle Scholar
- Grifo F, Rosenthal J: Biodiversity and human health. 1997, Washington DC, USA: Island PressGoogle Scholar
- Debprasad C, Sonali D, Sekhar C, Bhattacharya SK: Ethnomedicines for the development of anti-herpesvirus agents. Ethnomedicine: A Source of Complementary Therapeutics. 2010, 117-147.Google Scholar
- Bernstein FC, Koetzle TF, Williams GJ, Meyer EE, Brice MD, Rodgers JR, Kennard O, Shimanouchi T, Tasumi M: The Protein Data Bank: A Computer-based Archival File For Macromolecular Structures. J Mol Biol. 1977, 112: 535-10.1016/S0022-2836(77)80200-3.View ArticlePubMedGoogle Scholar
- Hussain Basha S, Prasad RN: In-Silico screening of Pleconaril and its novel substituted derivatives with Neuraminidase of H1N1 Influenza strain. BMC Research Notes. 2012, 5: 105-10.1186/1756-0500-5-105.View ArticlePubMedPubMed CentralGoogle Scholar
- Bolton E, Wang Y, Thiessen PA, Bryant SH: PubChem: Integrated Platform of Small Molecules and Biological Activities. Chapter 12. Annual Reports in Computational Chemistry, Volume 4. 2008, Washington, DC: American Chemical SocietyGoogle Scholar
- Thompson MA: ArgusLab 4.0.1. Seattle, WA: Planaria Software LLC, [http://www.arguslab.com]
- Rappe AK, Casewit CJ, Colwell KS, Goddard WA, Skiff WM: UFF, a full periodic table force field for molecular mechanics and molecular dynamics simulations. J Am Chem Soc. 1992, 114: 10024-10035. 10.1021/ja00051a040.View ArticleGoogle Scholar
- Rappe AK, Colwell KS, Casewit CJ: Application of a Universal force field to metal complexes. Inorg Chem. 1993, 32: 3438-3450. 10.1021/ic00068a012.View ArticleGoogle Scholar
- Rappe AK, Goddard WA: Charge Equilibration for molecular dyanamics simulations. J Phys Chem. 1991, 95: 3358-3363. 10.1021/j100161a070.View ArticleGoogle Scholar
- Casewit CJ, Colwell KS, Rappe’ AK: Application of a universal force field to organic molecules. J Am Chem Soc. 1992, 114: 10035-10046. 10.1021/ja00051a041.View ArticleGoogle Scholar
- Casewit CJ, Colwell KS, Rappe’ AK: Application of a universal force field to main group compounds. J Am Chem Soc. 1992, 114: 10046-10053. 10.1021/ja00051a042.View ArticleGoogle Scholar
- Ritchie DW: Recent progress and future directions in protein-protein docking. Curr Prot Pep Sci. 2008, 9 (1): 1-15. 10.2174/138920308783565741.View ArticleGoogle Scholar
- Accelrys Software Inc: Discovery studio Visualizer 3.1. 2011,http://accelrys.com/products/discovery-studio/visualization-download.php,Google Scholar
- Doina A, Whitbeck JC, de Leon MP, Lou H, Hannah BP, Cohen GH, Eisenberg RJ: Bimolecular Complementation Defines Functional Regions of Herpes Simplex Virus gB that Are Involved with gH/gL as a Necessary Step Leading to Cell Fusion. J Virol. 2010, 84 (8): 3825-3834. 10.1128/JVI.02687-09.View ArticleGoogle Scholar
- Hussain Basha S, Naresh Kumar K: Ligand and Structure based virtual screening studies to identify potent inhibitors against herpes virus targeting gB-gH-gL complex interface as a novel drug target. Open access sci rep. 2012, 1 (12): 566-Google Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1472-6882/13/85/prepub
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.