Molecular dynamics simulations of Zika virus NS3 helicase: Insights into RNA binding site activity

Comparative evaluation of methods for the prediction of protein-ligand binding sites

The accurate identification of protein-ligand binding sites is of critical importance in understanding and modulating protein function. Accordingly, ligand binding site prediction has remained a research focus for over three decades with over 50 methods developed and a change of paradigm from geometry-based to machine learning. In this work, we collate 13 ligand binding site predictors, spanning 30 years, focusing on the latest machine learning-based methods such as VN-EGNN, IF-SitePred, GrASP, PUResNet, and DeepPocket and compare them to the established P2Rank, PRANK and fpocket and earlier methods like PocketFinder, Ligsite and Surfnet. We benchmark the methods against the human subset of our new curated reference dataset, LIGYSIS. LIGYSIS is a comprehensive protein-ligand complex dataset comprising 30,000 proteins with bound ligands which aggregates biologically relevant unique protein-ligand interfaces across biological units of multiple structures from the same protein. LIGYSIS is an improvement for testing methods over earlier datasets like sc-PDB, PDBbind, binding MOAD, COACH420 and HOLO4K which either include 1:1 protein-ligand complexes or consider asymmetric units. Re-scoring of fpocket predictions by PRANK and DeepPocket display the highest recall (60%) whilst IF-SitePred presents the lowest recall (39%). We demonstrate the detrimental effect that redundant prediction of binding sites has on performance as well as the beneficial impact of stronger pocket scoring schemes, with improvements up to 14% in recall (IF-SitePred) and 30% in precision (Surfnet). Finally, we propose top-N+2 recall as the universal benchmark metric for ligand binding site prediction and urge authors to share not only the source code of their methods, but also of their benchmark.Scientific contributionsThis study conducts the largest benchmark of ligand binding site prediction methods to date, comparing 13 original methods and 15 variants using 10 informative metrics. The LIGYSIS dataset is introduced, which aggregates biologically relevant protein-ligand interfaces across multiple structures of the same protein. The study highlights the detrimental effect of redundant binding site prediction and demonstrates significant improvement in recall and precision through stronger scoring schemes. Finally, top-N+2 recall is proposed as a universal benchmark metric for ligand binding site prediction, with a recommendation for open-source sharing of both methods and benchmarks.

Classification of likely functional class for ligand binding sites identified from fragment screening

Fragment screening is used to identify binding sites and leads in drug discovery, but it is often unclear which binding sites are functionally important. Here, data from 37 experiments, and 1309 protein structures binding to 1601 ligands were analysed. A method to group ligands by binding sites is introduced and sites clustered according to profiles of relative solvent accessibility. This identified 293 unique ligand binding sites, grouped into four clusters (C1-4). C1 includes larger, buried, conserved, and population missense-depleted sites, enriched in known functional sites. C4 comprises smaller, accessible, divergent, missense-enriched sites, depleted in functional sites. A site in C1 is 28 times more likely to be functional than one in C4. Seventeen sites, which to the best of our knowledge are novel, in 13 proteins are identified as likely to be functionally important with examples from human tenascin and 5-aminolevulinate synthase highlighted. A multi-layer perceptron, and K-nearest neighbours model are presented to predict cluster labels for ligand binding sites with an accuracy of 96% and 100%, respectively, so allowing functional classification of sites for proteins not in this set. Our findings will be of interest to those studying protein-ligand interactions and developing new drugs or function modulators.

Learning from COVID-19: How drug hunters can prepare for the next pandemic

Over 3 years, the SARS-CoV-2 pandemic killed nearly 7 million people and infected more than 767 million globally. During this time, our very small company was able to contribute to antiviral drug discovery efforts through global collaborations with other researchers, which enabled the identification and repurposing of multiple molecules with activity against SARS-CoV-2 including pyronaridine tetraphosphate, tilorone, quinacrine, vandetanib, lumefantrine, cetylpyridinium chloride, raloxifene, carvedilol, olmutinib, dacomitinib, crizotinib, and bosutinib. We highlight some of the key findings from this experience of using different computational and experimental strategies, and detail some of the challenges and strategies for how we might better prepare for the next pandemic so that potential antiviral treatments are available for future outbreaks.

Discovery of New Zika Protease and Polymerase Inhibitors through the Open Science Collaboration Project OpenZika

The Zika virus (ZIKV) is a neurotropic arbovirus considered a global threat to public health. Although there have been several efforts in drug discovery projects for ZIKV in recent years, there are still no antiviral drugs approved to date. Here, we describe the results of a global collaborative crowdsourced open science project, the OpenZika project, from IBM's World Community Grid (WCG), which integrates different computational and experimental strategies for advancing a drug candidate for ZIKV. Initially, molecular docking protocols were developed to identify potential inhibitors of ZIKV NS5 RNA-dependent RNA polymerase (NS5 RdRp), NS3 protease (NS2B-NS3pro), and NS3 helicase (NS3hel). Then, a machine learning (ML) model was built to distinguish active vs inactive compounds for the cytoprotective effect against ZIKV infection. We performed three independent target-based virtual screening campaigns (NS5 RdRp, NS2B-NS3pro, and NS3hel), followed by predictions by the ML model and other filters, and prioritized a total of 61 compounds for further testing in enzymatic and phenotypic assays. This yielded five non-nucleoside compounds which showed inhibitory activity against ZIKV NS5 RdRp in enzymatic assays (IC50 range from 0.61 to 17 μM). Two compounds thermally destabilized NS3hel and showed binding affinity in the micromolar range (Kd range from 9 to 35 μM). Moreover, the compounds LabMol-301 inhibited both NS5 RdRp and NS2B-NS3pro (IC50 of 0.8 and 7.4 μM, respectively) and LabMol-212 thermally destabilized the ZIKV NS3hel (Kd of 35 μM). Both also protected cells from death induced by ZIKV infection in in vitro cell-based assays. However, while eight compounds (including LabMol-301 and LabMol-212) showed a cytoprotective effect and prevented ZIKV-induced cell death, agreeing with our ML model for prediction of this cytoprotective effect, no compound showed a direct antiviral effect against ZIKV. Thus, the new scaffolds discovered here are promising hits for future structural optimization and for advancing the discovery of further drug candidates for ZIKV. Furthermore, this work has demonstrated the importance of the integration of computational and experimental approaches, as well as the potential of large-scale collaborative networks to advance drug discovery projects for neglected diseases and emerging viruses, despite the lack of available direct antiviral activity and cytoprotective effect data, that reflects on the assertiveness of the computational predictions. The importance of these efforts rests with the need to be prepared for future viral epidemic and pandemic outbreaks.

Antiviral activity of ouabain against a Brazilian Zika virus strain

Zika virus (ZIKV) is an emerging arbovirus associated with neurological disorders. Currently, no specific vaccines or antivirals are available to treat the ZIKV infection. Ouabain, a cardiotonic steroid known as Na⁺/K⁺-ATPase inhibitor, has been previously described as an immunomodulatory substance by our group. Here, we evaluated for the first time the antiviral activity of this promising substance against a Brazilian ZIKV strain. Vero cells were treated with different concentrations of ouabain before and after the infection with ZIKV. The antiviral effect was evaluated by the TCID₅₀ method and RT-qPCR. Ouabain presented a dose-dependent inhibitory effect against ZIKV, mainly when added post infection. The reduction of infectious virus was accompanied by a decrease in ZIKV RNA levels, suggesting that the mechanism of ZIKV inhibition by ouabain occurred at the replication step. In addition, our in silico data demonstrated a conformational stability and favorable binding free energy of ouabain in the biding sites of the NS5-RdRp and NS3-helicase proteins, which could be related to its mechanism of action. Taken together, these data demonstrate the antiviral activity of ouabain against a Brazilian ZIKV strain and evidence the potential of cardiotonic steroids as promising antiviral agents.

Quantum Biochemistry and MM-PBSA Description of the ZIKV NS2B-NS3 Protease: Insights into the Binding Interactions beyond the Catalytic Triad Pocket

The Zika virus protease NS2B-NS3 has a binding site formed with the participation of a H51-D75-S135 triad presenting two forms, active and inactive. Studies suggest that the inactive conformation is a good target for the design of inhibitors. In this paper, we evaluated the co-crystallized structures of the protease with the inhibitors benzoic acid (5YOD) and benzimidazole-1-ylmethanol (5H4I). We applied a protocol consisting of two steps: first, classical molecular mechanics energy minimization followed by classical molecular dynamics were performed, obtaining stabilized molecular geometries; second, the optimized/relaxed geometries were used in quantum biochemistry and molecular mechanics/Poisson-Boltzmann surface area (MM-PBSA) calculations to estimate the ligand interactions with each amino acid residue of the binding pocket. We show that the quantum-level results identified essential residues for the stabilization of the 5YOD and 5H4I complexes after classical energy minimization, matching previously published experimental data. The same success, however, was not observed for the MM-PBSA simulations. The application of quantum biochemistry methods seems to be more promising for the design of novel inhibitors acting on NS2B-NS3.

In silico testing of flavonoids as potential inhibitors of protease and helicase domains of dengue and Zika viruses

Background

Dengue and Zika are two major vector-borne diseases. Dengue causes up to 25,000 deaths and nearly a 100 million cases worldwide per year, while the incidence of Zika has increased in recent years. Although Zika has been associated to fetal microcephaly and Guillain-Barré syndrome both it and dengue have common clinical symptoms such as severe headache, retroocular pain, muscle and join pain, nausea, vomiting, and rash. Currently, vaccines have been designed and antivirals have been identified for these diseases but there still need for more options for treatment. Our group previously obtained some fractions from medicinal plants that blocked dengue virus (DENV) infection in vitro. In the present work, we explored the possible targets by molecular docking a group of molecules contained in the plant fractions against DENV and Zika virus (ZIKV) NS3-helicase (NS3-hel) and NS3-protease (NS3-pro) structures. Finally, the best ligands were evaluated by molecular dynamic simulations.

Methods

To establish if these molecules could act as wide spectrum inhibitors, we used structures from four DENV serotypes and from ZIKV. ADFR 1.2 rc1 software was used for docking analysis; subsequently molecular dynamics analysis was carried out using AMBER20.

Results

Docking suggested that 3,5-dicaffeoylquinic acid (DCA01), quercetin 3-rutinoside (QNR05) and quercetin 3,7-diglucoside (QND10) can tightly bind to both NS3-hel and NS3-pro. However, after a molecular dynamics analysis, tight binding was not maintained for NS3-hel. In contrast, NS3-pro from two dengue serotypes, DENV3 and DENV4, retained both QNR05 and QND10 which converged near the catalytic site. After the molecular dynamics analysis, both ligands presented a stable trajectory over time, in contrast to DCA01. These findings allowed us to work on the design of a molecule called MOD10, using the QND10 skeleton to improve the interaction in the active site of the NS3-pro domain, which was verified through molecular dynamics simulation, turning out to be better than QNR05 and QND10, both in interaction and in the trajectory.

Discussion

Our results suggests that NS3-hel RNA empty binding site is not a good target for drug design as the binding site located through docking is too big. However, our results indicate that QNR05 and QND10 could block NS3-pro activity in DENV and ZIKV. In the interaction with these molecules, the sub-pocket-2 remained unoccupied in NS3-pro, leaving opportunity for improvement and drug design using the quercetin scaffold. The analysis of the NS3-pro in complex with MOD10 show a molecule that exerts contact with sub-pockets S1, S1', S2 and S3, increasing its affinity and apparent stability on NS3-pro.

Molecular dynamics simulations of allosteric motions and competitive inhibition of the Zika virus helicase

The 2015 Zika outbreak sparked major global concern and emphasized the reality and dangers still posed by mosquito borne pathogens. While efforts have been made to develop a vaccine and other therapeutics, there is still a great demand for antiviral drugs targeting Zika and other flaviviruses. The non-structural protein 3 (NS3) helicase is a vital component of the viral replication complex, tasked with unwinding the viral dsRNA molecule into single strands. Given this critical function, the Zika virus helicase is a potential therapeutic target and the focus of many ongoing research efforts. Using a combination of drug docking and molecular dynamics simulations, we have identified a list of competitive helicase inhibitors targeting the ATP hydrolysis site and have discovered a potential allosteric site capable of distorting both of the protein's active sites.

Discovery of human coronaviruses pan-papain-like protease inhibitors using computational approaches

The papain-like protease (PLpro) is vital for the replication of coronaviruses (CoVs), as well as for escaping innate-immune responses of the host. Hence, it has emerged as an attractive antiviral drug-target. In this study, computational approaches were employed, mainly the structure-based virtual screening coupled with all-atom molecular dynamics (MD) simulations to computationally identify specific inhibitors of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) PLpro, which can be further developed as potential pan-PLpro based broad-spectrum antiviral drugs. The sequence, structure, and functional conserveness of most deadly human CoVs PLpro were explored, and it was revealed that functionally important catalytic triad residues are well conserved among SARS-CoV, SARS-CoV-2, and middle east respiratory syndrome coronavirus (MERS-CoV). The subsequent screening of a focused protease inhibitors database composed of ∼7,000 compounds resulted in the identification of three candidate compounds, ADM_13083841, LMG_15521745, and SYN_15517940. These three compounds established conserved interactions which were further explored through MD simulations, free energy calculations, and residual energy contribution estimated by MM-PB(GB)SA method. All these compounds showed stable conformation and interacted well with the active residues of SARS-CoV-2 PLpro, and showed consistent interaction profile with SARS-CoV PLpro and MERS-CoV PLpro as well. Conclusively, the reported SARS-CoV-2 PLpro specific compounds could serve as seeds for developing potent pan-PLpro based broad-spectrum antiviral drugs against deadly human coronaviruses. Moreover, the presented information related to binding site residual energy contribution could lead to further optimization of these compounds.

Identification of novel human USP2 inhibitor and its putative role in treatment of COVID-19 by inhibiting SARS-CoV-2 papain-like (PLpro) protease

Human ubiquitin carboxyl-terminal hydrolase-2 (USP2) inhibitors, such as thiopurine analogs, have been reported to inhibit SARS-CoV papain-like proteases (PLpro). The PLpro have significant functional implications in the innate immune response during SARS-CoV-2 infection and considered an important antiviral target. Both proteases share strikingly similar USP fold with right-handed thumb-palm-fingers structural scaffold and conserved catalytic triad Cys-His-Asp/Asn. In this urgency situation of COVID-19 outbreak, there is a lack of in-vitro facilities readily available to test SARS-CoV-2 inhibitors in whole-cell assays. Therefore, we adopted an alternate route to identify potential USP2 inhibitor through integrated in-silico efforts. After an extensive virtual screening protocol, the best compounds were selected and tested. The compound Z93 showed significant IC50 value against Jurkat (9.67 μM) and MOTL-4 cells (11.8 μM). The binding mode of Z93 was extensively analyzed through molecular docking, followed by MD simulations, and molecular interactions were compared with SARS-CoV-2. The relative binding poses of Z93 fitted well in the binding site of both proteases and showed consensus π-π stacking and H-bond interactions with histidine and aspartate/asparagine residues of the catalytic triad. These results led us to speculate that compound Z93 might be the first potential chemical lead against SARS-CoV-2 PLpro, which warrants in-vitro evaluations.

Discovery of novel Hepatitis C virus inhibitor targeting multiple allosteric sites of NS5B polymerase

HCV is a viral infection posing a severe global threat when left untreated progress to end-stage liver disease, including cirrhosis and HCC. The NS5B polymerase of HCV is the most potent target that harbors four allosteric binding sites that could interfere with the HCV infection. We present the discovery of a novel synthetic compound that harbors the potential of NS5B polymerase inhibition. All eight compounds belonging to the benzothiazine family of heterocycles displayed no cellular cytotoxicity in HepG2 cells at nontoxic dose concentration (200 μM). Subsequently, among eight compounds of the series, merely compound 5b exhibited significant inhibition of the expression of the HCV NS5B gene as compared to DMSO control in semi-quantitative PCR. Based on our western blot result, 5b at the range of 50, 100 and 200 μM induced 20, 40, and 70% inhibition of NS5B protein respectively. To estimate the binding potential, 5b was docked at respective allosteric sites followed by molecular dynamics (MD) simulations for a period of 20 ns. In addition, binding free energy calculation by MM-GB/PBSA method revealed a conserved interaction profile of residues lining the allosteric sites in agreement with the reported NS5B co-crystallized inhibitors. The presented results provide important information about a novel compound 5b which may facilitate the the discovery of novel inhibitors that tends to target multiple sites on NS5B polymerase.

Déjà vu: Stimulating open drug discovery for SARS-CoV-2

In the past decade we have seen two major Ebola virus outbreaks in Africa, the Zika virus in Brazil and the Americas and the current pandemic of coronavirus disease (COVID-19), caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). There is a strong sense of déjà vu because there are still no effective treatments. In the COVID-19 pandemic, despite being a new virus, there are already drugs suggested as active in in vitro assays that are being repurposed in clinical trials. Promising SARS-CoV-2 viral targets and computational approaches are described and discussed. Here, we propose, based on open antiviral drug discovery approaches for previous outbreaks, that there could still be gaps in our approach to drug discovery.

Structure-based screening and validation of bioactive compounds as Zika virus methyltransferase (MTase) inhibitors through first-principle density functional theory, classical molecular simulation and QM/MM affinity estimation

Recent Zika virus (ZIKV) outbreak and association with human diseases such as neurological disorders have raised global health concerns. However, in the absence of an approved anti-ZIKV drug has generated urgency for the drug development against ZIKV infection. Here, structure-based virtual screening of 8589 bioactive compounds, screened at the substrate-binding site of ZIKV nonstructural 5 (NS5)-based structure N-terminal methyltransferase (MTase) domain followed by ADMET (absorption, distribution, metabolism, excretion and toxicity) profiling concluded the four potential lead inhibitors, i.e. (4-acetylamino-benzenesulfonylamino)-acetic acid (F3342-0450), 3-(5-methylfuran-2-yl)-N-(4-sulfamoylphenyl)propanamide (F1736-0142), 8-(2-hydroxy-ethylamino)-1,3-dimethyl-7-(3-methyl-benzyl)-3,7-dihydro-purine-2,6-dione (F0886-0080) and N-[4-(aminosulfonyl)phenyl]-2,3-dihydro-1,4-benzodioxine-2-carboxamide (F0451-2187). Collectively, extra precision docking and Density Functional Theory(DFT) calculations studies identified the F3342-0450 molecule, having strong interactions on the active site of MTase, further supported by molecular dynamics simulation, binding affinity and hybrid QM/MM calculations, suggest a new drug molecule for the antiviral drug development against ZIKV infection. Communicated by Ramaswamy H. Sarma.

RNA-Dependent Structures of the RNA-Binding Loop in the Flavivirus NS3 Helicase

The flavivirus NS3 protein is a helicase that has pivotal functions during the viral genome replication process, where it unwinds double-stranded RNA and translocates along the nucleic acid polymer in a nucleoside triphosphate hydrolysis-dependent mechanism. Crystallographic and computational studies of the flavivirus NS3 helicase have identified the RNA-binding loop as an interesting structural element that may function as a component of the RNA-enhanced NTPase activity observed for this family of helicases. Microsecond-long unbiased molecular dynamics and extensive replica exchange umbrella sampling simulations of the Zika NS3 helicase have been performed to investigate the RNA dependence of this loop's structural conformations. Specifically, the effect of the bound single-stranded RNA (ssRNA) oligomer on the putative "open" and "closed" conformations of this loop is studied. In the Apo substrate state, the two loop conformations are nearly isoergonic (ΔAO→C = -0.22 kcal mol-1), explaining the structural ambiguity observed in Apo NS3h crystal structures. The bound ssRNA is seen to stabilize the "open" conformation (ΔAO→C = 1.97 kcal mol-1) through direct protein-RNA interactions at the top of the loop. Interestingly, a small ssRNA oligomer bound over 13 Å away from the loop is seen to affect the free energy surface to favor the "open" structure, while minimizing barriers between the two states. Both the mechanism of the "open" to "closed" transition and important residues of the RNA-binding loop structures are characterized. From these results, point mutations that are hypothesized to stabilize the "closed" RNA-binding loop and negatively impact RNA-binding and the RNA-enhanced NTPase activity are posited.

Perspectives towards antiviral drug discovery against Ebola virus

Ebola virus disease (EVD), caused by Ebola viruses, resulted in more than 11 500 deaths according to a recent 2018 WHO report. With mortality rates up to 90%, it is nowadays one of the most deadly infectious diseases. However, no Food and Drug Administration‐approved Ebola drugs or vaccines are available yet with the mainstay of therapy being supportive care. The high fatality rate and absence of effective treatment or vaccination make Ebola virus a category‐A biothreat pathogen. Fortunately, a series of investigational countermeasures have been developed to control and prevent this global threat. This review summarizes the recent therapeutic advances and ongoing research progress from research and development to clinical trials in the development of small‐molecule antiviral drugs, small‐interference RNA molecules, phosphorodiamidate morpholino oligomers, full‐length monoclonal antibodies, and vaccines. Moreover, difficulties are highlighted in the search for effective countermeasures against EVD with additional focus on the interplay between available in silico prediction methods and their evidenced potential in antiviral drug discovery.

PARP-1 mediated cell death is directly activated by ZIKV infection

Zika virus (ZIKV) has emerged as a severe health threat due to its association with microcephaly. It has been reported that the strong cytopathic effects, including cell-cycle arrest and cell death are responsible for the nervous system disease. However, the mechanisms by which ZIKV infection induced cell death were largely unknown. Here, we reported that cell death is readily detected after ZIKV infection as indicated by PI staining and the reduction of cell viability. Importantly, cell death can be induced by overexpression of ZIKV NS3 protein alone but not the other non-structure proteins. Mass spectrometry analysis revealed that NS3 bond to and activated PARP-1. In agreement with these observations, we found that PARP-1 was massively activated during ZIKV infection and the intracellular ATP and NAD⁺ concentrations rapidly declined. Finally, PARP-1 knockdown simultaneously restrained ZIKV infection-induced cell death and ablated host restriction of virus infection. Our finding indicates that PARP-1 activation is an important cellular event during ZIKV infection, which contributes to the cell death.

Structure Based Identification of Potential Inhibitors of NS3 Protein of Zika Virus

Background:

The re-emerging Zika virus has posed a serious threat to human health due to its association with the neurological disorders. The NS3 protein of Zika virus plays a pivotal role in the genome replication and thus may prove to be a critical target for the drug designing studies.

Objective:

The present study was conceptualized to analyze the crystal structure of NS3 protein of Zika virus followed by the identification of it’s potential inhibitors.

Methods:

Crystal structure of the NS3 protein was evaluated in detail. Docking of the NS3 protein was done with 130 different ligands including dengue virus inhibitors and their similar compounds along with some approved drugs. The drug likeliness properties were checked for non drug compounds.

Results:

Structural analysis of the NS3 protein revealed three important sites namely ATP- and RNAbinding sites as well as a central cavity. The selected ten ligands (ZINC05487635, ZINC0092398, ZINC13345444, 4-methoxyphenyl 4-chloro-3-nitrobenzoate, Luteolin, Ivermectin, Suramin, Dasatinib, Panduratin A, and ARDP0009) showed a higher binding affinity for the NS3 protein and good drug likeliness properties.

Conclusion:

These inhibitors could possibly act as potential lead molecules for future drug designing studies. Our present computational data is envisaged to be useful for gathering experimental evidences towards the development of potential therapeutic molecules against this arthropod mediated pathogen.

Zika virus circumvents host innate immunity by targeting the adaptor proteins MAVS and MITA

Recently, Zika virus (ZIKV) has generated extraordinary concern because of its severe neurotoxicity. Disturbingly, there is no vaccine or specific drug to prevent or treat the diseases caused by ZIKV infection. Thus, it is extremely urgent to characterize the pathogenesis of ZIKV. It has been documented that ZIKV can evade antiviral responses of host cells. Here, we demonstrate that ZIKV strain SZ-WIV01 down-regulates the production of type I IFN and IFN-stimulated genes along with the expression of mitochondrial antiviral signaling protein (MAVS) and mediator of IFN regulatory factor 3 activation (MITA). In the mechanism, ZIKV nonstructural (NS) 3 and NS2B3 negatively regulate IFN-related retinoic acid-inducible gene I-like receptor signaling pathway by targeting MAVS and MITA, respectively. Overexpression of ZIKV NS3 and NS2B3 dramatically inhibits expression of IFN-β. ZIKV NS3 interacts with MAVS, and NS2B3 interacts with MITA, which catalyzes K48-linked polyubiquitination of MAVS and MITA for degradation. Further investigations suggest that ZIKV NS2B3 impairs polyinosinic:polycytidylic acid-triggered K63-linked polyubiquitination of MITA, thereby subverting the activation of downstream sensors. Our study reveals an undiscovered mechanism for ZIKV to escape the innate immune response, providing new insights into clinical study of vaccines or effective drugs.-Li, W., Li, N., Dai, S., Hou, G., Guo, K., Chen, X., Yi, C., Liu, W., Deng, F., Wu, Y., Cao, X. Zika virus circumvents host innate immunity by targeting the adaptor proteins MAVS and MITA.

Molecular docking and simulation of Zika virus NS3 helicase

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.

In silico structural elucidation of RNA-dependent RNA polymerase towards the identification of potential Crimean-Congo Hemorrhagic Fever Virus inhibitors

The Crimean-Congo Hemorrhagic Fever virus (CCHFV) is a segmented negative single-stranded RNA virus (-ssRNA) which causes severe hemorrhagic fever in humans with a mortality rate of ~50%. To date, no vaccine has been approved. Treatment is limited to supportive care with few investigational drugs in practice. Previous studies have identified viral RNA dependent RNA Polymerase (RdRp) as a potential drug target due to its significant role in viral replication and transcription. Since no crystal structure is available yet, we report the structural elucidation of CCHFV-RdRp by in-depth homology modeling. Even with low sequence identity, the generated model suggests a similar overall structure as previously reported RdRps. More specifically, the model suggests the presence of structural/functional conserved RdRp motifs for polymerase function, the configuration of uniform spatial arrangement of core RdRp sub-domains, and predicted positively charged entry/exit tunnels, as seen in sNSV polymerases. Extensive pharmacophore modeling based on per-residue energy contribution with investigational drugs allowed the concise mapping of pharmacophoric features and identified potential hits. The combination of pharmacophoric features with interaction energy analysis revealed functionally important residues in the conserved motifs together with in silico predicted common inhibitory binding modes with highly potent reference compounds.

Discovery and Computational Analyses of Novel Small Molecule Zika Virus Inhibitors

Zika virus (ZIKV), one of the flaviviruses, has attracted worldwide attention since its large epidemics around Brazil. Association of ZIKV infection with microcephaly and neurological problems such as Guillain–Barré syndrome has prompted intensive pathological investigations. However, there is still a long way to go on the discovery of effective anti-ZIKV therapeutics. In this study, an in silico screening of the National Cancer Institute (NCI) diversity set based on ZIKV NS3 helicase was performed using a molecular docking approach. Selected compounds with drug-like properties were subjected to cell-based antiviral assays resulting in the identification of two novel lead compounds (named Compounds 1 and 2). They inhibited ZIKV infection with IC₅₀ values at the micro-molar level (8.5 μM and 15.2 μM, respectively). Binding mode analysis, absolute binding free energy calculation, and structure–activity relationship studies of these two compounds revealed their possible interactions with ZIKV NS3 helicase, suggesting a mechanistic basis for further optimization. These two novel small molecules may represent new leads for the development of inhibitory drugs against ZIKV.

A comparative study of the efficiency of HCV NS3/4A protease drugs against different HCV genotypes using in silico approaches

AIMS

To investigate the efficacy of Direct Acting Antivirals (DAAs) in the treatment of different Hepatitis C Virus (HCV) genotypes.

MAIN METHODS

Homology modeling is used to predict the 3D structures of different genotypes while molecular docking is employed to predict genotype - drug interactions (Binding Mode) and binding free energy (Docking Score).

KEY FINDINGS

Simeprevir (TMC435) and to a lesser degree MK6325 are the best drugs among the studied drugs. The predicted affinity of drugs against genotype 1a is always better than other genotypes. P2-P4 macrocyclic drugs show better performance against genotypes 2, 3 and 5. Macrocyclic drugs are better than linear drugs.

SIGNIFICANCE

HCV is one of the major health problems worldwide. Until the discovery of DAAs, HCV treatment faced many failures. DAAs target key functional machines of the virus life cycle and shut it down. NS3/4A protease is an important target and several drugs have been designed to inhibit its functions. There are several NS3/4A protease drugs approved by Food and Drug Administration (FDA). Unfortunately, the virus exhibits resistance against these drugs. This study is significant in elucidating that no one drug is able to treat different genotypes with the same efficiency. Therefore, treatment should be prescribed based on the HCV genotype.

In silico approaches to Zika virus drug discovery

INTRODUCTION

After the WHO declared Zika virus (ZIKV) as a public health emergency of international concern, intense research for the development of vaccines and drugs has been undertaken, leading to the development of several candidates. Areas covered: This review discusses the developments achieved so far by computational methods in the discovery of candidate compounds targeting ZIKV proteins, i.e. the envelope and capsid structural proteins, the NS3 helicase/protease, and the NS5 methyltransferase/RNA-dependent RNA polymerase. Expert opinion: Research for effective drugs against ZIKV is still in a very early discovery phase. Notwithstanding the intense efforts for the development of new drugs and the identification of several promising candidates by using different approaches, including computational methods, so far only a few candidates have been experimentally tested. An important caveat of anti-flavivirus drug development is represented by the difficult of reproducing the in vivo microenvironment of the replication complex, which may lead to discrepancies between in vitro results and experimental evaluation in vivo. Moreover, anti-ZIKV drugs have the additional requirement of an excellent safety profile in pregnancy and ability to diffuse to different tissues, including the central nervous system, the testis, and the placenta.

The A-Z of Zika drug discovery

Despite the recent outbreak of Zika virus (ZIKV), there are still no approved treatments, and early-stage compounds are probably many years away from approval. A comprehensive A-Z review of the recent advances in ZIKV drug discovery efforts is presented, highlighting drug repositioning and computationally guided compounds, including discovered viral and host cell inhibitors. Promising ZIKV molecular targets are also described and discussed, as well as targets belonging to the host cell, as new opportunities for ZIKV drug discovery. All this knowledge is not only crucial to advancing the fight against the Zika virus and other flaviviruses but also helps us prepare for the next emerging virus outbreak to which we will have to respond.

Allostery in the dengue virus NS3 helicase: Insights into the NTPase cycle from molecular simulations

The C-terminus domain of non-structural 3 (NS3) protein of the Flaviviridae viruses (e.g. HCV, dengue, West Nile, Zika) is a nucleotide triphosphatase (NTPase) -dependent superfamily 2 (SF2) helicase that unwinds double-stranded RNA while translocating along the nucleic polymer. Due to these functions, NS3 is an important target for antiviral development yet the biophysics of this enzyme are poorly understood. Microsecond-long molecular dynamic simulations of the dengue NS3 helicase domain are reported from which allosteric effects of RNA and NTPase substrates are observed. The presence of a bound single-stranded RNA catalytically enhances the phosphate hydrolysis reaction by affecting the dynamics and positioning of waters within the hydrolysis active site. Coupled with results from the simulations, electronic structure calculations of the reaction are used to quantify this enhancement to be a 150-fold increase, in qualitative agreement with the experimental enhancement factor of 10–100. Additionally, protein-RNA interactions exhibit NTPase substrate-induced allostery, where the presence of a nucleotide (e.g. ATP or ADP) structurally perturbs residues in direct contact with the phosphodiester backbone of the RNA. Residue-residue network analyses highlight pathways of short ranged interactions that connect the two active sites. These analyses identify motif V as a highly connected region of protein structure through which energy released from either active site is hypothesized to move, thereby inducing the observed allosteric effects. These results lay the foundation for the design of novel allosteric inhibitors of NS3.

Non-structural protein 3 (NS3) is a Flaviviridae (e.g. Hepatitis C, dengue, and Zika viruses) helicase that unwinds double stranded RNA while translocating along the nucleic polymer during viral genome replication. As a member of superfamily 2 (SF2) helicases, NS3 utilizes the free energy of nucleotide triphosphate (NTP) binding, hydrolysis, and product unbinding to perform its functions. While much is known about SF2 helicases, the pathways and mechanisms through which free energy is transduced between the NTP hydrolysis active site and RNA binding cleft remains elusive. Here we present a multiscale computational study to characterize the allosteric effects induced by the RNA and NTPase substrates (ATP, ADP, and P_i) as well as the pathways of short-range, residue-residue interactions that connect the two active sites. Results from this body of molecular dynamics simulations and electronic structure calculations are highlighted in context to the NTPase enzymatic cycle, allowing for development of testable hypotheses for validation of these simulations. Our insights, therefore, provide novel details about the biophysics of NS3 and guide the next generation of experimental studies.