- Research article
- Open Access
Molecular docking and simulation of Zika virus NS3 helicase
© The Author(s) 2019
- Received: 7 November 2017
- Accepted: 2 May 2019
- Published: 17 May 2019
The Zika virus (ZIKV) has gained attention for the last few years due to the congenital microcephaly and Guillain–Barre Syndrome that resulted in humans. The non-structural protein-3 (NS3) helicase of ZIKV play an important role in viral RNA replication. In this article, we performed hundred nanosecond molecular dynamics simulation and molecular docking of the NS3 helicase of ZIKV with 1,4-benzothiazine derivatives. The root mean square deviation (RMSD) analyses showed the stability of the NS3 helicase. The simulation showed that the flexible and rigid domains of the protein play a crucial role during the RNA replication process. All such domains with ligand binding pockets can be targeted for drug design. The molecular docking showed that the strong hydrogen bonding and arene-cation interactions are responsible for the binding between NS3 and 1,4-benzothiazine derivatives, which provides a new dimension for potent drug design for ZIKV.
- Nonstructural protein-3 Helicase
- Molecular dynamics simulation
- Molecular docking
Zika virus (ZIKV) is a mosquito-borne flavivirus like yellow fever virus, dengue virus (DENV), West Nile virus (WNV) and chikungunya . It contains a single-stranded positive RNA with a total of 10,794-nucleotides in its genome . The ZIKV transmission occurs through numerous Aedes spp. mosquitoes, including Aedes africanus, Aedes luteocephalus, Aedes hensilli, and Aedes aegypti [3–6]. Initially the ZIKV was observed in rhesus monkey amid observation of yellow fever in the zika woodland in Uganda in the year 1947 and it was initially reported in people in 1952 [7, 8]. Later on in 1954, it was observed in humans with febrile diseases in West Africa [7, 8]. In 2014 several cases of ZIKV were recorded in Indonesia, Micronesia, Thailand, Philippines, French Polynesia and Easter Island in South Pacific [9, 10]. From May, 2015 onward there was a boom in the spread of ZIKV in Brazil and from there it spread to other countries [11, 12].
In early 2016, the World Health Organization declared the ZIKV spread as a international public health emergency due to the presence of congenital neurological disorders like microcephaly, Guillain–Barre syndrome and cranial nerve dysfunction . The exact mechanism of action of ZIKV that how it causes microcephaly is not known but research efforts are in progress to establish the link between ZIKV and microcephaly [10, 14, 15]. Generally, the mosquito borne flavivirus pathogenesis begins with its replication in dendritic cells close to the site of inoculation from where it move to the lymph centers and blood circulation system . After the initial viral attack, the virions can be observed in the blood, while the viral RNA can been seen after eleven days of ailment . The helicases are an omnipresent exceedingly diverse group of proteins that carry out an astonishing variety of functions in cells . These are ATP-ases which depend on nucleic acid and have the potential to unwind DNA or RNA duplex substrates . Because of this unwinding properties, they are significantly involved in every process of cells linked with DNA like replication, repair, transcription, translation, synthesis of ribosomes, RNA maturation and splicing and the procedures related with the nuclear export . Recently a number of structural as well as non-structural proteins from ZIKV that have enzymatic functions are resolved at high resolution through various techniques and they are useful target for drug design . The NS3 helicase is one of the most thoroughly studied antiviral drug target [20, 21]. In NS3 of flaviviruses, the enzymatic activity is coupled with the C-terminal region, namely an RNA helicase (NS3-Hel) concerned with genome replication and RNA synthesis [22, 23]. The NS3-Hel is part of the superfamily helicases, and its inhibition in DENV makes the virus unable to replicate [24, 25]. The two motor domains formed in the helicase part of the NS3 proteins are identical with all other helicases of superfamily 1 and 2 but the third domain is very different from other helicases. The binding site for ATP is located between the two domains while the RNA binds between the third and motor domain . The inchworm like movement of NS3 helicase in a 3′ to 5′ directions along one strand of RNA was observed in detailed mechanistic studies which involved the displacing of the complementary strand in an ATP utilized reaction. The same type of structural and enzymatic studies of DENV helicase suggest that there is no difference between the mechanism of hepacivirus and flavivirus [27, 28]. The non-nucleoside based viral inhibitors are an active area of research and benzothiazine based compounds have been studied against several flaviviruses [21, 29, 30]. They contains a nitrogen and a sulfur atom in their basic nucleus that makes them a good interacting agent inside the active site of an enzyme . Different research groups have shown that benzothiazine are suitable inhibitors of helicases [31, 32]. Currently there are two MD simulation studies available that explore the dynamics of ZIKV NS3 helicase but to get further insight docking novel inhibitors and further simulation studies are required . In this study we targeted the ZIKV NS3 helicase with 1, 4-benzothiazine analogues (Additional file 1: Figure S1) using molecular dynamics simulation and molecular docking methods. The reason for choosing 1,4-benzothiazine based derivatives was that a number of research groups have used 1,4-benzothiazine analogues both computationally and experimentally against different viral helicase enzymes [31, 34]. Therefore, we design these derivatives and used it for molecular docking and simulation studies. We hope that this study will also be vital in the effort to find suitable inhibitors against the spread of this devastating viral disease.
The crystal structure of the NS3 helicase of ZIKV was obtained from Protein Data Bank with the PDB ID 5JRZ . PyMOL software version 1.7 was used for visualization and checking the protein . The ligands were searched from the literature and newly reported compounds were selected for molecular docking . For molecular docking, we used the Molecular Operation Environment (MOE) 2014 docking software program . The prediction of the active site in the protein was carried out with MOE 2014 software, site finder tool which identified the various sites in which one was selected for docking procedure. The removal of water molecules was done from the crystal structure. The addition of missing hydrogen atoms was done, correction of charges and assigning of the correct hybridization state of each residue was done through the preparation program of MOE. The corrected protonation was done using 3D protonate module present in MOE with Generalized Born/Volume Integral (GB/VI) electrostatic function. To explore the potential binding site, the complete structure of the enzyme was used as a receptor. For each ligand, generation of multiple conformations were done by applying a selected torsion angles to all the rotatable bonds in each ligand. For each drug candidate about thirty conformations were generated. For each ligand and the receptor, the accepted conformations were scored by using the London Dock scoring function which calculates the free energy for the binding ligand in a given conformation .
The binding affinities were calculated using generalized-Born volume integral/weighted surface area (GBVI/WSA) method present in MOE 2014. Generalized Born interaction energy is the non-bonded interaction energy between the receptor molecule and the ligand. It is composed of Coulomb electrostatic interaction, Van der Waals, and implicit solvent interaction energies. The binding affinity was calculated for each hit after energy minimization, and reported in unit of kcal/mol . The PCA was performed as reported previously .
Molecular dynamics simulations
The crystal structure of NS3 helicase of ZIKV, PDB ID 5JRZ having 1.62 Å resolution was used for atomic coordinates for both Apo and in docked complex with derivative 7 saved as PDB file, which was constructed through SYBYLs-X2.1.1 . For molecular dynamics simulation the AMBER12 software package [43, 44] was utilized adjusting all the parameters. The addition of hydrogen atoms was done through LEaP module of AMBER-12. For the maintenance of neutral condition of the system the counter-ions were utilized. A truncated octahedral box of TIP3P molecules of water was used for the solvation of both the systems. The cutoff distance is kept 8 Å for the determination of pairwise interactions (van der Waals and direct coulombic interactions). In order to compute the long range electrostatic interactions the particle mesh Ewald (PME) method  of AMBER12 was used while the ff12SB force field was selected for the finding of the intramolecular interactions. For preparations runs a Langevin thermostat was utilized with friction constant of 1 ps−1 while for the production runs a Berendsen thermostat was used . All MD simulations were accelerated using the CUDA version of PMEMD in GPU cores of NVIDIAs Tesla K20. The 500-step steepest descent minimization and 2000-step conjugated gradient minimization were performed for initial minimization with macromolecules frozen in order to avoid further structural clash in the solvated system. The whole system was again passed through the minimization process at the end, keeping 1000-step steepest descent minimization and 19,000-step conjugate gradient minimization. The adjustment of parameters after energy minimization at 400 ps heating up and 200 ps equilibration in the NVT ensemble at 310 K were done before the start of MD simulations in the NPT ensemble at 310 K. In order to carry out a comparison of Apo and complex various simulated trajectories at different nanoseconds were collected at 310 K during the 100 ns simulation timescale.
Molecular docking of 1,4-benzothiazine derivatives
Stability of the simulated systems
Structural fluctuation of the NS3 protein
Drug designing involved two opposite views about the drug targets in literature; one is that inhibitors inhibit the most rigid residues and the second view is that inhibitors target highly flexible residues . The rigid residues and flexible regions in enzymes have specified role in various biological reactions . The molecular recognition and catalytic properties are linked with the mobility of specific residues while production of various structures in beta-folds is due to the rigid residues, that are essential for interaction with the substrate and its catalysis.
Principal component analysis
Superposition of Apo and complex before and after simulation
In case of the complex form, before the simulation time, the inhibitors was represented by green color while the portions of three domains were represented by dark brown color and after the simulation the ligand was shown by blue color and the protein mobile spheres were represented by cyan color (Fig. 6b, c). When we performed the superposition of the NS3 helicase with the inhibitor, there is a difference in the position of the mobile residues which showed that the inhibitor produced conformational changes when it binds inside the active site of the NS3 helicase which was previously confirmed from RMSF analysis. Secondly, again the amount of displacement and orientation in the spheres is different from one another in the same protein. While an internal motion towards the center of axis from left and right side is also observed in the complex form after the simulation time.
Helicase enzyme of various viruses is an important drug target and in this study, we reported the utilization of 1-4 benzothiazine derivatives as inhibitor of ZIKV helicase for the docking study. The docking results indicated that these compounds were in interactions mainly with residues Thr201, Lys200, Arg462 and Arg459 of the NS3 which was reported hotspots by making hydrogen bonds and arene-cation interactions. The 1,4-benzothiazine derivative 7 showed maximum interactions with best docking scores and binding energies while its validation was further carried out with molecular dynamics simulation for 100 ns and the results were found excellent including, RMSD, RMSF values, principal component analysis and superposition expressed the stability of the drug inside the binding pocket. From the above study, we can conclude that this class of compounds can be used as possible novel inhibitors of ZIKV helicase (NS3) and it’s a way forward toward virtual screening and pharmacophore mapping tools to explore the treasures of biological important class of compounds to inhibit the genome replication of ZIKV. The molecular dynamic simulation confirmed that the docked conformation is reliable. Binding energy calculations through the MOE docking program showed that van der Waal interaction and hydrogen bonding provided the most substantial force for the binding of the inhibitor.
The authors extend their appreciation to the Deanship of Scientific Research at King Khalid University for its funding this prolific research Group No. (R.G.P.2/23/40/2019).
King Khalid University grant awarded to Dr. Yahia Mabkhot.
SB and NA contributed equally. SB, NA and AR performed the molecular docking and simulation. NA, KH, SB, AU, AM, AA and YM wrote the manuscript. All authors read and approved the final manuscript.
All the authors declare that they have no competing interests.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Moulin E, Selby K, Cherpillod P et al (2016) Simultaneous outbreaks of dengue, chikungunya and Zika virus infections: diagnosis challenge in a returning traveller with nonspecific febrile illness. New Microbes New Infect. 11:6–7View ArticleGoogle Scholar
- Cunha MS, Esposito DLA, Rocco IM et al (2016) First complete genome sequence of Zika virus (Flaviviridae, Flavivirus) from an Autochthonous transmission in Brazil. Genome Announc 4:2015–2016. https://doi.org/10.1128/genomeA.00032-16 View ArticleGoogle Scholar
- Gardner LM, Chen N, Sarkar S (2016) Global risk of Zika virus depends critically on vector status of Aedes albopictus. Lancet Infect Dis. 16:522–523View ArticleGoogle Scholar
- Wong PSJ, Li M, Zhi I, Chong CS et al (2013) Aedes (Stegomyia) albopictus (Skuse): a potential vector of Zika virus in Singapore. PLoS Negl Trop Dis 1:1. https://doi.org/10.1371/journal.pntd.0002348 View ArticleGoogle Scholar
- Leandro H, Dutra C, Rocha MN et al (2016) Wolbachia blocks currently circulating Zika virus isolates in Brazilian Aedes aegypti mosquitoes. Cell Host Microbe 19:1–4. https://doi.org/10.1016/j.chom.2016.04.021 View ArticleGoogle Scholar
- Bajpai S, Nadkar MY (2016) Zika virus infection, the recent menace of the Aedes mosquito. J Assoc Physicians India 64:42–45PubMedGoogle Scholar
- Dick GWA (1952) Zika virus (I). Isolations and serological specificity. Trans R Soc Trop Med Hyg 46:509–520. https://doi.org/10.1016/0035-9203(52)90042-4 View ArticlePubMedGoogle Scholar
- Wikan N, Smith DR (2016) Zika virus: history of a newly emerging arbovirus. Lancet Infect Dis 16:e119–e126View ArticleGoogle Scholar
- Musso D (2015) Zika virus transmission from French Polynesia to Brazil. Emerg Infect Dis 21:1887–1889View ArticleGoogle Scholar
- Cauchemez S, Besnard M, Bompard P et al (2016) Association between Zika virus and microcephaly in French Polynesia, 2013–2015: a retrospective study. Lancet 387:2125–2132. https://doi.org/10.1016/s0140-6736(16)00651-6 View ArticlePubMedPubMed CentralGoogle Scholar
- Bonaldo MC, Ribeiro IP, Lima NS et al (2016) Isolation of infective Zika virus from urine and saliva of patients in Brazil. PLoS Negl Trop Dis 1:1. https://doi.org/10.1371/journal.pntd.0004816 View ArticleGoogle Scholar
- Syed S, Mabkhot YN, Naeem A, et al (2018) Zika virus, microcephaly and its possible global spread. In: Current topics in ZikaGoogle Scholar
- Gulland A (2016) Zika virus is a global public health emergency, declares WHO. BMJ 352:i657. https://doi.org/10.1136/bmj.i657 View ArticlePubMedGoogle Scholar
- De Góes Cavalcanti LP, Tauil PL, Alencar CH et al (2016) Zika virus infection, associated microcephaly, and low yellow fever vaccination coverage in Brazil: is there any causal link. J Infect Dev Ctries 10:563–566. https://doi.org/10.3855/jidc.8575 View ArticlePubMedGoogle Scholar
- Chan JFW, Choi GKY, Yip CCY et al (2016) Zika fever and congenital Zika syndrome: an unexpected emerging arboviral disease. J Infect 72:507–524View ArticleGoogle Scholar
- Diamond MS, Shrestha B, Mehlhop E et al (2003) Innate and adaptive immune responses determine protection against disseminated infection by West Nile encephalitis virus. Viral Immunol 16:259–278. https://doi.org/10.1089/088282403322396082 View ArticlePubMedGoogle Scholar
- Filipe AR, Martins CMV, Rocha H (1973) Laboratory infection with Zika virus after vaccination against yellow fever. Archiv für die gesamte Virusforschung 43:315–319. https://doi.org/10.1007/bf01556147 View ArticlePubMedGoogle Scholar
- Gorbalenya AE, Koonin EV (1993) Helicases: amino acid sequence comparisons and structure-function relationships. Curr Opin Struct Biol 3:419–429. https://doi.org/10.1016/S0959-440X(05)80116-2 View ArticleGoogle Scholar
- Badshah S, Naeem A, Mabkhot Y (2017) The new high resolution crystal structure of NS2B-NS3 protease of Zika virus. Viruses 9:7. https://doi.org/10.3390/v9010007 View ArticlePubMed CentralGoogle Scholar
- De Clercq E (2014) Current race in the development of DAAs (direct-acting antivirals) against HCV. Biochem Pharmacol 89:441–452View ArticleGoogle Scholar
- Fang J, Li H, Kong D et al (2016) Structure-based discovery of two antiviral inhibitors targeting the NS3 helicase of Japanese encephalitis virus. Sci Rep 6:34550. https://doi.org/10.1038/srep34550 View ArticlePubMedPubMed CentralGoogle Scholar
- Frick DN, Lam AMI (2006) Understanding helicases as a means of virus control. Curr Pharm Des 12:1315–1338. https://doi.org/10.2174/138161206776361147 View ArticlePubMedPubMed CentralGoogle Scholar
- Luo D, Vasudevan SG, Lescar J (2015) The flavivirus NS2B-NS3 protease-helicase as a target for antiviral drug development. Antiviral Res 118:148–158View ArticleGoogle Scholar
- Matusan AE, Pryor MJ, Davidson AD, Wright PJ (2001) Mutagenesis of the Dengue virus type 2 NS3 protein within and outside helicase motifs: effects on enzyme activity and virus replication. J Virol 75:9633–9643. https://doi.org/10.1128/JVI.75.20.9633 View ArticlePubMedPubMed CentralGoogle Scholar
- Singleton MR, Dillingham MS, Wigley DB (2007) Structure and mechanism of helicases and nucleic acid translocases. Annu Rev Biochem 76:23–50. https://doi.org/10.1146/annurev.biochem.76.052305.115300 View ArticlePubMedGoogle Scholar
- Luo D, Xu T, Watson RP et al (2008) Insights into RNA unwinding and ATP hydrolysis by the flavivirus NS3 protein. EMBO J 27:3209–3219. https://doi.org/10.1038/emboj.2008.232 View ArticlePubMedPubMed CentralGoogle Scholar
- Frick DN (2007) The hepatitis C virus NS3 protein: a model RNA helicase and potential drug target. Curr Issues Mol Biol 9:1–20PubMedPubMed CentralGoogle Scholar
- Gu M, Rice CM (2010) Three conformational snapshots of the hepatitis C virus NS3 helicase reveal a ratchet translocation mechanism. Proc Natl Acad Sci USA 107:521–528. https://doi.org/10.1073/pnas.0913380107 View ArticlePubMedGoogle Scholar
- de Vicente J, Hendricks RT, Smith DB et al (2009) Non-nucleoside inhibitors of HCV polymerase NS5B. Part 2: synthesis and structure-activity relationships of benzothiazine-substituted quinolinediones. Bioorganic Med Chem Lett 19:3642–3646. https://doi.org/10.1016/j.bmcl.2009.05.004 View ArticleGoogle Scholar
- Badshah Syed L, Naeem A (2016) Bioactive thiazine and benzothiazine derivatives: green synthesis methods and their medicinal importance. Molecules 21:1054View ArticleGoogle Scholar
- Lin K (2010) Development of novel antiviral therapies for hepatitis c virus. Virol Sin 25:246–266View ArticleGoogle Scholar
- Rehman S, Ashfaq UA, Javed T (2011) Antiviral drugs against hepatitis C virus. Genet Vaccines Ther 9:11. https://doi.org/10.1186/1479-0556-9-11 View ArticlePubMedPubMed CentralGoogle Scholar
- Mottin M, Braga RC, da Silva RA et al (2017) Molecular dynamics simulations of Zika virus NS3 helicase: insights into RNA binding site activity. Biochem Biophys Res Commun. https://doi.org/10.1016/j.bbrc.2017.03.070 View ArticlePubMedGoogle Scholar
- Hendricks RT, Fell JB, Blake JF et al (2009) Non-nucleoside inhibitors of HCV NS5B polymerase. Part 1: synthetic and computational exploration of the binding modes of benzothiadiazine and 1,4-benzothiazine HCV NS5b polymerase inhibitors. Bioorg Med Chem Lett 19:3637–3641. https://doi.org/10.1016/j.bmcl.2009.04.119 View ArticlePubMedGoogle Scholar
- Jain R, Coloma J, Garcia-Sastre A, Aggarwal AK (2016) Structure of the NS3 helicase from Zika virus. Nat Struct Mol Biol 23:752–754. https://doi.org/10.1038/nsmb.3258 View ArticlePubMedPubMed CentralGoogle Scholar
- DeLano WL (2002) The PyMOL molecular graphics system. Schrödinger LLC version 1. http://www.pymol.org
- Chu J-J, Hu B-L, Liao Z-Y, Zhang X-G (2016) Copper-catalyzed three-component tandem cyclization for one-pot synthesis of 1,4-benzothiazines. J Org Chem 81:8647–8652. https://doi.org/10.1021/acs.joc.6b01668 View ArticlePubMedGoogle Scholar
- Molecular Operating Environment (MOE), Version 2009.10, Chemical Computing Group Inc., Montreal, Quebec, CanadaGoogle Scholar
- Kumar YN, Kumar PS, Sowjenya G et al (2012) Comparison and correlation of binding mode of ATP in the kinase domains of Hexokinase family. Bioinformation 8:543–547. https://doi.org/10.6026/97320630008543 View ArticlePubMedPubMed CentralGoogle Scholar
- Labute P (2008) The generalized born/volume integral implicit solvent model: estimation of the free energy of hydration using London dispersion instead of atomic surface area. J Comput Chem 29:1693–1698. https://doi.org/10.1002/jcc.20933 View ArticlePubMedGoogle Scholar
- ur Rahman M, Liu H, Wadood A, Chen H-F (2016) Allosteric mechanism of cyclopropylindolobenzazepine inhibitors for HCV NS5B RdRp via dynamic correlation network analysis. Mol BioSyst 12:3280–3293. https://doi.org/10.1039/C6MB00521G View ArticlePubMedGoogle Scholar
- Ash S, Cline MA, Homer RW et al (1997) SYBYL line notation (SLN): a versatile language for chemical structure representation. J Chem Inf Model 37:71–79. https://doi.org/10.1021/ci960109j View ArticleGoogle Scholar
- Salomon-Ferrer R, Case DA, Walker RC (2013) An overview of the Amber biomolecular simulation package. Wiley Interdiscip Rev Comput Mol Sci 3:198–210. https://doi.org/10.1002/wcms.1121 View ArticleGoogle Scholar
- Pearlman DA, Case DA, Caldwell JW et al (1995) AMBER, a package of computer programs for applying molecular mechanics, normal mode analysis, molecular dynamics and free energy calculations to simulate the structural and energetic properties of molecules. Comput Phys Commun 91:1–41. https://doi.org/10.1016/0010-4655(95)00041-D View ArticleGoogle Scholar
- Darden T, York D, Pedersen L (1993) Particle mesh Ewald: an N·log(N) method for Ewald sums in large systems. J Chem Phys 98:10089. https://doi.org/10.1063/1.464397 View ArticleGoogle Scholar
- Ruhle V (2007) Berendsen and Nose-Hoover thermostats. Thesis 1–4Google Scholar
- Ramharack P, Ogundate S, Soliman M (2017) Delving into Zika virus structural dynamics—a closer look at NS3 helicase loop flexibility and its role in drug discovery. RSC Adv 7:22133View ArticleGoogle Scholar
- Badshah SL, Naeem A, Mabkhot Y (2016) Molecular dynamics simulation of cholera toxin A-1 polypeptide. Open Chem 14:188–196. https://doi.org/10.1515/chem-2016-0021 View ArticleGoogle Scholar
- Tian H, Ji X, Yang X et al (2016) The crystal structure of Zika virus helicase: basis for antiviral drug design. Protein Cell 7:450–454. https://doi.org/10.1007/s13238-016-0275-4 View ArticlePubMedPubMed CentralGoogle Scholar