Research Article | Open Access
Zekun Liu, Junpeng Zhao, Weichen Li, Xinkun Wang, Jingxuan Xu, Jin Xie, Ke Tao, Li Shen, Ran Zhang, "Molecular Docking of Potential Inhibitors for Influenza H7N9", Computational and Mathematical Methods in Medicine, vol. 2015, Article ID 480764, 8 pages, 2015. https://doi.org/10.1155/2015/480764
Molecular Docking of Potential Inhibitors for Influenza H7N9
As a new strain of virus emerged in 2013, avian influenza A (H7N9) virus is a threat to the public health, due to its high lethality and pathogenicity. Furthermore, H7N9 has already generated various mutations such as neuraminidase R294K mutation which could make the anti-influenza oseltamivir less effective or ineffective. In this regard, it is urgent to develop new effective anti-H7N9 drug. In this study, we used the general H7N9 neuraminidase and oseltamivir-resistant influenza virus neuraminidase as the acceptors and employed the small molecules including quercetin, chlorogenic acid, baicalein, and oleanolic acid as the donors to perform the molecular docking for exploring the binding abilities between these small molecules and neuraminidase. The results showed that quercetin, chlorogenic acid, oleanolic acid, and baicalein present oseltamivir-comparable high binding potentials with neuraminidase. Further analyses showed that R294K mutation in neuraminidase could remarkably decrease the binding energies for oseltamivir, while other small molecules showed stable binding abilities with mutated neuraminidase. Taken together, the molecular docking studies identified four potential inhibitors for neuraminidase of H7N9, which might be effective for the drug-resistant mutants.
As the WHO reported on June 27, 2014, there were 450 laboratory-confirmed cases of avian influenza A (H7N9) virus human infection including 165 deaths (http://www.who.int). Most patients presented with severe pneumonia, acute respiratory distress syndrome (ARDS), and multiple organ failure [1, 2]. Sequence analysis indicated that H7N9 gene descended from avian species, and what is more, the other six genes originated from H9N2 [1–3]. Since the H7N9 virus could bind to both the avian and human receptors, the potential risk of spreading between avian species and human was highly increased . Previously, it was observed that the Asp701Asn mutation in PB2 and the Ser31Asn mutation in hemagglutinin increased the adaptability to mammals, and the Ser31Asn mutation of M2 contributed at least partly to the resistance of antivirus drugs including amantadine and rimantadine [3, 4]. These mutations made the treatment of H7N9 virus infected patients more difficult and increased the threats to public health. The neuraminidase inhibitor oseltamivir was reported to be effective and safe for therapy of influenza A  and employed in a number of H7N9 cases. Although a number of cases were cured [6, 7], oseltamivir was still invalid for a part of cases and the efficiency is still debatable [8, 9]. Recently, Wu et al. solved the structures for neuraminidases of A/Anhui/1/2013 (N9, R294) and A/Shanghai/1/2013 (N9, K294), which showed that the mutation of R294K could induce conformation changes and generate oseltamivir resistance through interfering with the binding of oseltamivir carboxylate . Furthermore, the authors proposed that the drug resistance caused by mutation of R294K in H7N9 was more serious than that caused by mutation of H274Y in H7N1 . Taken together, the development of new inhibitors for mutated neuraminidase is urgently needed.
Recently, molecular modeling and computational chemistry based computer-aided drug design provided great help for modern drug development . Software programs such as AutoDock were widely used to search potential inhibitor for protein targets . Previously, traditional Chinese medicine has been demonstrated to have antivirus ability according to the clinical data, and a number of small molecules such as quercetin, chlorogenic acid, oleanolic acid, and baicalein were regarded as active molecules . For example, quercetin could serve as the inhibitor of Bcl-2 and Bcl-xL through binding with the BH3 domain, which promoted cancer cell apoptosis [14, 15]. Chlorogenic acid could protect the ApoE knockout mice against atherosclerosis through accelerating the cholesterol efflux from macrophages [16, 17]. Baicalein was identified as the inhibitor of neuraminidase for pandemic 2009 H1N1 and also inhibited interleukin-1 beta and TNF-alpha mediated inflammation [18, 19]. Through inducing autophagy, oleanolic acid could inhibit the proliferation and invasiveness in K-ras transformed cells and induce Nrf-2 related antioxidant to reduce ethanol-induced liver injury [20, 21]. However, whether these small molecules could inhibit the neuraminidase of H7N9 is still unclear.
In this study, we employed the strategy of molecule docking to explore potential H7N9 neuraminidase inhibitors from small molecules including quercetin, chlorogenic acid, oleanolic acid, and baicalein, which were bioactive components in traditional Chinese medicine. It was observed that these small molecules showed high binding energies with neuraminidases from A/Anhui/1/2013. The binding energies for these small molecules were comparable with oseltamivir, which indicated that the molecules were potential inhibitors for neuraminidase. Furthermore, the results showed a weak binding energy between oseltamivir and neuraminidase from A/Shanghai/1/2013, which was consistent with the oseltamivir resistance of H7N9 A/Shanghai/1/2013. However, the four molecules quercetin, chlorogenic acid, oleanolic acid, and baicalein showed robust binding energies with neuraminidase from A/Shanghai/1/2013. Taken together, these molecules were potential neuraminidase inhibitors which might be helpful for anti-influenza drug development.
2. Materials and Methods
2.1. Molecular Structure Preparing
To perform the molecular docking between neuraminidase and potential inhibitors, we employed the H7N9 neuraminidase structures of A/Anhui/1/2013 (PDB code: 4MWQ) and A/Shanghai/1/2013 (PDB code: 4MWW), which were solved by Wu et al. at the resolutions of 2.00 Å and 1.90 Å, respectively . The PDB files for these proteins were downloaded from the RCSB Protein Data Bank (available at http://www.rcsb.org) . The water and nonprotein molecules in the PDB files were removed, while only neuraminidase monomer structures were reserved. All the visualizations of molecular structures were performed with PyMOL, which is state-of-the-art software, to present the supercomplexes of proteins and small molecules (Version 1.7.2, available at http://pymol.org) . Four small molecules including quercetin, chlorogenic acid, oleanolic acid, and baicalein from traditional Chinese medicine were selected for molecular docking based on the following principles: (1) the molecules should be identified from traditional Chinese medicine; (2) the molecules should have been studied previously and be commercially available. The structures in mol2 file format of the five small molecules, including oseltamivir carboxylate (ZINC03929509), quercetin (ZINC03869685), chlorogenic acid (ZINC02138728), oleanolic acid (ZINC03785416), and baicalein (ZINC03943903), were retrieved from the ZINC database, which is a free public resource for ligand molecules (available at http://zinc.docking.org) [24, 25].
2.2. Molecular Docking
In this study, we employed AutoDock software (Version 4.2, available at http://autodock.scripps.edu) to perform molecular docking. The AutoDock is an automated docking program designed for prediction of the binding among small molecules such as substrates or drug candidates and the receptor with known 3D structure . AutoDockTools software (Version 1.5.6, also available at http://autodock.scripps.edu), the graphical user interface (GUI) program for setting up AutoDock software, was used to prepare the procedures of docking . The volume of the grid box was set as 60 × 50 × 60 with the default 0.375 Å spacing, and the number of docking runs was set as 100. The neuraminidase structures were used as the receptor for docking with the small molecules. Since the active sites of neuraminidase in influenza A were highly conserved, the docking pocket was set as the eight key residues including Arg-118, Asp-151, Arg-152, Arg-224, Glu-276, Arg-292, Arg-371, and Tyr-406 [26–28]. The chemoinformatic analyses of potential toxicities and promiscuities for the five small molecules were performed with FAF-Drugs3 software (available at http://fafdrugs3.mti.univ-paris-diderot.fr) .
3. Result and Discussion
3.1. Characterizing the Structures of the Small Molecules for Molecular Docking
In this study, five small molecules including oseltamivir carboxylate (Figure 1(a)), quercetin (Figure 1(b)), baicalein (Figure 1(c)), chlorogenic acid (Figure 1(d)), and oleanolic acid (Figure 1(e)) were employed for docking with neuraminidase to explore the potential inhibitor. Previously, it was observed that the oseltamivir phosphate is the prodrug without antivirus activity while its product of metabolism oseltamivir carboxylate is active . Thus, in this study, the oseltamivir carboxylate was employed for molecular docking. The chemoinformatic analyses were performed with the FAF-Drugs3 software  to further understand the properties of the small molecules. Since the hydrogen bonds (H-bonds) were critical for the interactions and binding affinity between the protein and small molecules, we analyzed the H-bond donors and acceptors of these small molecules. It was observed that oseltamivir carboxylate and quercetin, chlorogenic acid and baicalein have five and six H-bond donors, respectively, while oleanolic acid has two (Table 1). Baicalein has 11 H-bond acceptors, while chlorogenic acid, quercetin, oseltamivir carboxylate, and oleanolic acid have 9, 7, 6, and 3 H-bond acceptors, respectively (Table 1). These results showed that quercetin, chlorogenic acid, oleanolic acid, and baicalein had the potential to form H-bonds with proteins. Furthermore, the potential toxicities and promiscuities were analyzed by the FAF-Drugs3 software . It was observed that baicalein, chlorogenic acid, quercetin, and oseltamivir carboxylate have low toxicities and promiscuities according to the threshold of < 3 and tPSA > 75 (Table 1). Furthermore, oleanolic acid has the of 7.49 and tPSA of 60.36 (Table 1). However, oleanolic acid was proved as a drug for human liver disorders treatment . Thus, it seemed that four small molecules had limited toxicities and promiscuities and were proper for further drug development.
3.2. Molecular Docking between NA and Small Molecules
The molecular docking was performed among the neuraminidases from H7N9 A/Anhui/1/2013 and A/Shanghai/1/2013. Oseltamivir carboxylate could serve as a control to evaluate the neuraminidase binding ability of other molecules. The conformation variations, binding energies, and inhibition constants for the molecular docking among neuraminidase and small molecules were presented in Figure 2 and Tables 2 and 3. The calculated root mean square deviations (RMSDs) indicated that the conformation variations among neuraminidases binding with different small molecules were limited, while A/Shanghai/1/2013 neuraminidase had moderate large variations (Tables 2-3). It was observed that the A/Anhui/1/2013 neuraminidase binding energies for oseltamivir carboxylate and quercetin were almost the same, which were −9.38 kcal/mol and −9.41 kcal/mol, respectively, while chlorogenic acid, oleanolic acid, and baicalein have higher binding energies as −12.23 kcal/mol, −11.20 kcal/mol, and −11.45 kcal/mol, separately. The inhibition constants calculated for oseltamivir carboxylate, quercetin, chlorogenic acid, oleanolic acid, and baicalein were 132.26, 126.87, 1.09, 6.19, and 4.02 nM, respectively (Table 2). For the A/Shanghai/1/2013 neuraminidase, the binding energies for oseltamivir carboxylate, quercetin, chlorogenic acid, oleanolic acid, and baicalein were −6.91 kcal/mol, −9.24 kcal/mol, −9.82 kcal/mol, −9.51 kcal/mol, and −10.32 kcal/mol. It was obvious that oseltamivir carboxylate showed weaker binding ability for A/Shanghai/1/2013 neuraminidase than A/Anhui/1/2013, which was consistent with the observation that A/Shanghai/1/2013 has resistance for oseltamivir . Other molecules including chlorogenic acid, oleanolic acid, and baicalein also presented the decrease of binding energies; quercetin showed robust binding abilities for neuraminidases from H7N9 A/Anhui/1/2013 and A/Shanghai/1/2013. The inhibition constants calculated for oseltamivir carboxylate, quercetin, chlorogenic acid, oleanolic acid, and baicalein were , 168.62, 63.6, 106.11, and 27.11 nM, respectively (Table 2), which also indicated that H7N9 neuraminidase had significant oseltamivir resistance. Furthermore, the detailed results for the top 3 scored poses for each small molecule were presented in Supplementary Table S1 (see Supplementary Table S1 in Supplementary Material available online at http://dx.doi.org/10.1155/2015/480764).
3.3. Structural Insights of the Docking Complex
To further dissect the interaction between neuraminidase and the small molecules, the structures in complex derived from the best dock results with neuraminidase for oseltamivir carboxylate (Figure 3(a)), quercetin (Figure 3(b)), chlorogenic acid (Figure 4(a)), oleanolic acid (Figure 4(b)), and baicalein (Figure 5) derived from molecular docking were shown. From the molecule docking among A/Anhui/1/2013 neuraminidase and small molecules, it was observed that there were 18 H-bonds in the neuraminidase-chlorogenic acid complex (Figure 4(a)), while quercetin (Figure 3(b)) has 14 H-bonds with neuraminidase. Oseltamivir carboxylate (Figure 3(a)) and baicalein (Figure 5) both have 12 H-bonds with neuraminidase, while oleanolic acid (Figure 4(b)) has 5. The H-bond residues for neuraminidase from A/Anhui/1/2013 and A/Shanghai/1/2013 were provided in Tables 4 and 5, respectively. These results indicated that chlorogenic acid, quercetin, and baicalein could intensively bind neuraminidase through H-bonds. For A/Shanghai/1/2013 neuraminidase, chlorogenic acid could form 17 H-bonds (Figure 4(a)), which indicated that chlorogenic acid has high potential of neuraminidase inhibition. It was interesting that the H-bonds for baicalein (Figure 5) and oleanolic acid (Figure 4(b)) were increased to 15 and 8 for A/Shanghai/1/2013 neuraminidase compared to those from A/Anhui/1/2013, while there were 9 and 11 H-bonds for oseltamivir carboxylate (Figure 3(a)) and quercetin (Figure 3(b)), respectively. These results indicated that the oseltamivir resistance caused by R294K mutation for neuraminidase in A/Shanghai/1/2013 might be generated from the H-bond loss, while other molecules might overcome the resistance. Furthermore, other known inhibitors of neuraminidase including zanamivir, peramivir, and laninamivir had similar structures and presented similar inhibition abilities for neuraminidase from A/Anhui/1/2013 . However, they also showed significantly decreased inhibition for A/Shanghai/1/2013 . Since the small molecules including quercetin, chlorogenic acid, oleanolic acid, and baicalein had longer structures and there was considerable room for improvements, we anticipated that they had great potential to overcome the drug resistance of H7N9 mutant.
Zekun Liu and Junpeng Zhao should be regarded as joint first authors.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
This study was supported by grants from the Innovative Training Center for College Student of Life Sciences (201301012), the National Natural Science Foundation of China (81100054), and the Construct Program of the Key Discipline of Basic Medicine in Hunan Province.
The Supplementary Material of Table S1 provide the detailed results of the top 3 scored poses for the complexes of neuraminidases from A/Anhui/1/2013 and A/Shanghai/1/2013 docking with each small molecule.
- R. Gao, B. Cao, Y. Hu et al., “Human infection with a novel avian-origin influenza A (H7N9) virus,” The New England Journal of Medicine, vol. 368, no. 20, pp. 1888–1897, 2013.
- Y. Wang, Z. Dai, H. Cheng et al., “Towards a better understanding of the novel avian-origin H7N9 influenza A virus in China,” Scientific Reports, vol. 3, article 2318, 2013.
- Y. Chen, W. Liang, S. Yang et al., “Human infections with the emerging avian influenza A H7N9 virus from wet market poultry: Clinical analysis and characterisation of viral genome,” The Lancet, vol. 381, no. 9881, pp. 1916–1925, 2013.
- Q. Liu, L. Lu, Z. Sun, G.-W. Chen, Y. Wen, and S. Jiang, “Genomic signature and protein sequence analysis of a novel influenza A(H7N9) virus that causes an outbreak in humans in China,” Microbes and Infection, vol. 15, no. 6-7, pp. 432–439, 2013.
- J. J. Treanor, F. G. Hayden, P. S. Vrooman et al., “Efficacy and safety of the oral neuraminidase inhibitor oseltamivir in treating acute influenza: a randomized controlled trial,” The Journal of the American Medical Association, vol. 283, no. 8, pp. 1016–1024, 2000.
- S.-Y. Chang, P.-H. Lin, J.-C. Tsai, C.-C. Hung, and S.-C. Chang, “The first case of H7N9 influenza in Taiwan,” The Lancet, vol. 381, no. 9878, p. 1621, 2013.
- Q. Li, L. Zhou, M. Zhou et al., “Epidemiology of human infections with avian influenza A(H7N9) virus in China,” The New England Journal of Medicine, vol. 370, pp. 520–532, 2014.
- X. Liu, T. Li, Y. Zheng, K.-W. Wong, S. Lu, and H. Lu, “Poor responses to oseltamivir treatment in a patient with influenza A (H7N9) virus infection,” Emerging Microbes & Infections, vol. 2, article e27, 2013.
- S. Lu, Y. Zheng, T. Li et al., “Clinical findings for early human cases of influenza A(H7N9) virus infection, Shanghai, Shina,” Emerging Infectious Diseases, vol. 19, no. 7, pp. 1142–1146, 2013.
- Y. Wu, Y. Bi, C. J. Vavricka et al., “Characterization of two distinct neuraminidases from avian-origin human-infecting H7N9 influenza viruses,” Cell Research, vol. 23, no. 12, pp. 1347–1355, 2013.
- G. R. Marshall, “Computer-aided drug design,” Annual Review of Pharmacology and Toxicology, vol. 27, pp. 193–213, 1987.
- G. M. Morris, H. Ruth, W. Lindstrom et al., “Software news and updates AutoDock4 and AutoDockTools4: automated docking with selective receptor flexibility,” Journal of Computational Chemistry, vol. 30, no. 16, pp. 2785–2791, 2009.
- X. Wang, W. Jia, and A. Zhao, “Anti-influenza agents from plants and traditional Chinese medicine,” Phytotherapy Research, vol. 20, no. 5, pp. 335–341, 2006.
- T. Thangasamy, S. Sittadjody, S. Lanza-Jacoby, P. R. Wachsberger, K. H. Limesand, and R. Burd, “Quercetin selectively inhibits bioreduction and enhances apoptosis in melanoma cells that overexpress tyrosinase,” Nutrition and Cancer, vol. 59, no. 2, pp. 258–268, 2007.
- A. Primikyri, M. V. Chatziathanasiadou, E. Karali et al., “Direct binding of Bcl-2 family proteins by quercetin triggers its pro-apoptotic activity,” ACS Chemical Biology, vol. 9, no. 12, pp. 2737–2741, 2014.
- C. Wu, H. Luan, X. Zhang et al., “Chlorogenic acid protects against atherosclerosis in ApoE−/− mice and promotes cholesterol efflux from RAW264.7 macrophages,” PLoS ONE, vol. 9, no. 9, Article ID e95452, 2014.
- E. Heitman and D. K. Ingram, “Cognitive and neuroprotective effects of chlorogenic acid,” Nutritional Neuroscience, 2014.
- M.-J. Hour, S.-H. Huang, C.-Y. Chang et al., “Baicalein, ethyl acetate, and chloroform extracts of Scutellaria baicalensis inhibit the neuraminidase activity of pandemic 2009 H1N1 and seasonal influenza A viruses,” Evidence-Based Complementary and Alternative Medicine, vol. 2013, Article ID 750803, 11 pages, 2013.
- C.-J. Hsieh, K. Hall, T. Ha, C. Li, G. Krishnaswamy, and D. S. Chi, “Baicalein inhibits IL-1beta- and TNF-alpha-induced inflammatory cytokine production from human mast cells via regulation of the NF-kappaB pathway,” Clinical and Molecular Allergy, vol. 5, article 5, 2007.
- J. Liu, L. Zheng, L. Ma et al., “Oleanolic acid inhibits proliferation and invasiveness of Kras-transformed cells via autophagy,” The Journal of Nutritional Biochemistry, vol. 25, no. 11, pp. 1154–1160, 2014.
- J. Liu, X. Wang, R. Liu et al., “Oleanolic acid co-administration alleviates ethanol-induced hepatic injury via Nrf-2 and ethanol-metabolizing modulating in rats,” Chemico-Biological Interactions, vol. 221, pp. 88–98, 2014.
- P. W. Rose, A. Prli , C. Bi et al., “The RCSB Protein Data Bank: views of structural biology for basic and applied research and education,” Nucleic Acids Research, vol. 43, no. D1, pp. D345–D356, 2015.
- W. L. DeLano, The PyMOL Molecular Graphics System, DeLano Scientific, Palo Alto, Calif, USA, 2002.
- J. J. Irwin and B. K. Shoichet, “ZINC—a free database of commercially available compounds for virtual screening,” Journal of Chemical Information and Modeling, vol. 45, no. 1, pp. 177–182, 2005.
- J. J. Irwin, T. Sterling, M. M. Mysinger, E. S. Bolstad, and R. G. Coleman, “ZINC: a free tool to discover chemistry for biology,” Journal of Chemical Information and Modeling, vol. 52, no. 7, pp. 1757–1768, 2012.
- X. Xu, X. Zhu, R. A. Dwek, J. Stevens, and I. A. Wilson, “Structural characterization of the 1918 influenza virus H1N1 neuraminidase,” Journal of Virology, vol. 82, no. 21, pp. 10493–10501, 2008.
- P. M. Mitrasinovic, “Advances in the structure-based design of the influenza A neuraminidase inhibitors,” Current Drug Targets, vol. 11, no. 3, pp. 315–326, 2010.
- J. Gong, W. Xu, and J. Zhang, “Structure and functions of influenza virus neuraminidase,” Current Medicinal Chemistry, vol. 14, no. 1, pp. 113–122, 2007.
- D. Lagorce, J. Maupetit, J. Baell et al., “The FAF-Drugs2 server: a multistep engine to prepare electronic chemical compound collections,” Bioinformatics, vol. 27, no. 14, Article ID btr333, pp. 2018–2020, 2011.
- B. E. Davies, “Pharmacokinetics of oseltamivir: an oral antiviral for the treatment and prophylaxis of influenza in diverse populations,” The Journal of Antimicrobial Chemotherapy, vol. 65, supplement 2, pp. ii5–ii10, 2010.
- J. Liu, “Pharmacology of oleanolic acid and ursolic acid,” Journal of Ethnopharmacology, vol. 49, no. 2, pp. 57–68, 1995.
Copyright © 2015 Zekun Liu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.