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.

NS2B-NS3 protease inhibitors as promising compounds in the development of antivirals against Zika virus: A systematic review.

Zika virus (ZIKV) infections are associated with severe neurological complications and are a global public health concern. There are no approved vaccines or antiviral drugs to inhibit ZIKV replication. NS2B-NS3 protease (NS2B-NS3 pro), which is essential for viral replication, is a promising molecular target for anti-ZIKV drugs. We conducted a systematic review to identify compounds with promising effects against ZIKV; we discussed their pharmacodynamic and pharmacophoric characteristics. The online search, performed using the PubMed/MEDLINE and SCOPUS databases, yielded 56 articles; seven relevant studies that reported nine promising compounds with inhibitory activity against ZIKV NS2B-NS3 pro were selected. Of these, five (niclosamide, nitazoxanide, bromocriptine, temoporfin, and novobiocin) are currently available on the market and have been tested for off-label use against ZIKV. The 50% inhibitory concentration values of these compounds for the inhibition of NS2B-NS3 pro ranged at 0.38-21.6 µM; most compounds exhibited noncompetitive inhibition (66%). All compounds that could inhibit the NS2B-NS3 pro complex showed potent in vitro anti-ZIKV activity with a 50% effective concentration ranging 0.024-50 µM. The 50% cytotoxic concentration of the compounds assayed using A549, Vero, and WRL-69 cell lines ranged at 0.6-1388.02 µM and the selectivity index was 3.07-1698. This review summarizes the most promising antiviral agents against ZIKV that have inhibitory activity against viral proteases.

Chikungunya nsP4 homology modeling reveals a common motif with Zika and Dengue RNA polymerases as a potential therapeutic target.

Among the diseases transmitted by vectors, there are those caused by viruses named arboviruses (arthropod-borne viruses). In past years, viruses transmitted by mosquitoes have been of relevance in global health, such as Chikungunya (CHIKV), Dengue (DENV), and Zika (ZIKV), which have Aedes aegypti as a common vector, thus raising the possibility of multi-infection. Previous reports have described the general structure of RNA-dependent RNA polymerases termed right-hand fold, which is conserved in positive single-stranded RNA viruses. Here, we report a comparison between sequences and the computational structure of RNA-dependent RNA polymerases from CHIKV, DENV, and ZIKV and the conserved sites to be considered for the design of an antiviral drug against the three viruses. We show that the sequential identity between consensus sequences from CHIKV and DENV is 8.1% and the similarity is 15.1%; the identity between CHIKV and ZIKV is 9.3%, and the similarity is 16.6%; and the identity between DENV and ZIKV is 68.6%, and the similarity is 79.2%. Nevertheless, the structural alignment shows that the root-mean-square deviation (RMSD) measurement value in general structure comparison between CHIKV RdRp and ZIKV RdRp was 1.248 Å, RMSD between CHIKV RdRp and DENV RdRp was 1.070 Å, and RMSD between ZIKV RdRp and DENV RdRp was 1.106 Å. Despite the low identity and similarity of CHIKV sequence with DENV and ZIKV, we show that A, B, C, and E motifs are structurally well conserved. These structural similarities offer a window into drug design against these arboviruses giving clues about critical target sites.

Insights on Dengue and Zika NS5 RNA-dependent RNA polymerase (RdRp) inhibitors.

Over recent years, many outbreaks caused by (re)emerging RNA viruses have been reported worldwide, including life-threatening Flaviviruses, such as Dengue (DENV) and Zika (ZIKV). Currently, there is only one licensed vaccine against Dengue, Dengvaxia®. However, its administration is not recommended for children under nine years. Still, there are no specific inhibitors available to treat these infectious diseases. Among the flaviviral proteins, NS5 RNA-dependent RNA polymerase (RdRp) is a metalloenzyme essential for viral replication, suggesting that it is a promising macromolecular target since it has no human homolog. Nowadays, several NS5 RdRp inhibitors have been reported, while none inhibitors are currently in clinical development. In this context, this review constitutes a comprehensive work focused on RdRp inhibitors from natural, synthetic, and even repurposing sources. Furthermore, their main aspects associated with the structure-activity relationship (SAR), proposed mechanisms of action, computational studies, and other topics will be discussed in detail.

In vitro study of Hesperetin and Hesperidin as inhibitors of zika and chikungunya virus proteases.

The potential outcome of flavivirus and alphavirus co-infections is worrisome due to the development of severe diseases. Hundreds of millions of people worldwide live under the risk of infections caused by viruses like chikungunya virus (CHIKV, genus Alphavirus), dengue virus (DENV, genus Flavivirus), and zika virus (ZIKV, genus Flavivirus). So far, neither any drug exists against the infection by a single virus, nor against co-infection. The results described in our study demonstrate the inhibitory potential of two flavonoids derived from citrus plants: Hesperetin (HST) against NS2B/NS3pro of ZIKV and nsP2pro of CHIKV and, Hesperidin (HSD) against nsP2pro of CHIKV. The flavonoids are noncompetitive inhibitors and the determined IC50 values are in low µM range for HST against ZIKV NS2B/NS3pro (12.6 ± 1.3 µM) and against CHIKV nsP2pro (2.5 ± 0.4 µM). The IC50 for HSD against CHIKV nsP2pro was 7.1 ± 1.1 µM. The calculated ligand efficiencies for HST were > 0.3, which reflect its potential to be used as a lead compound. Docking and molecular dynamics simulations display the effect of HST and HSD on the protease 3D models of CHIKV and ZIKV. Conformational changes after ligand binding and their effect on the substrate-binding pocket of the proteases were investigated. Additionally, MTT assays demonstrated a very low cytotoxicity of both the molecules. Based on our results, we assume that HST comprise a chemical structure that serves as a starting point molecule to develop a potent inhibitor to combat CHIKV and ZIKV co-infections by inhibiting the virus proteases.

Identification of a Potential Zika Virus Inhibitor Targeting NS5 Methyltransferase Using Virtual Screening and Molecular Dynamics Simulations.

The NS5 methyltransferase (MTase) has been reported as an attractive molecular target for antivirals discovery against the Zika virus (ZIKV). Here, we report structure-based virtual screening of 42 390 structures from the Development Therapeutics Program (DTP) AIDS Antiviral Screen Database. Among the docked compounds, ZINC1652386 stood out due to its high affinity for MTase in comparison to the cocrystallized ligand MS2042, which interacts with the Asp146 residue in the MTase binding site by hydrogen bonding. Subsequent molecular dynamics simulations predicted that this compound forms a stable complex with MTase within 50 ns. Thus, ZINC1652386 may represent a promising ZIKV methyltransferase inhibitor.

Identification of Zika Virus NS2B-NS3 Protease Inhibitors by Structure-Based Virtual Screening and Drug Repurposing Approaches.

The NS2B-NS3 protease has been identified as an attractive target for drug development against Zika virus (ZIKV) and combined drug repurposing and structure-based virtual screening has improved the development of antiviral drugs. In this study, we performed a structure-based virtual screening of 1861 Food and Administration (FDA) approved drugs available in DrugBank by the selection and docking validation of crystal structure of ZIKV NS2B-NS3 protease (PDB ID 5H4I ) using Glide and DOCK 6 software. The antihistaminic chlorcyclizine (Grid score -24.8 kcal/mol) exhibited the most promising interaction with NS2B-NS3 protease in comparison to crystallography ligand (Grid score -15.6 kcal/mol) by interaction to Tyr161 by hydrophobic interactions in the binding site of NS2B-NS3 which is recognized as an important amino acid in substrate molecular recognition. Cytotoxicity and global antiviral activity assay in Vero cells by MTT method showed that chlorcyclizine reduced the ZIKV induced cytopathic effect (EC50 of 69.0 ± 7.3 µM and SI = 1.9), and explicit molecular dynamics simulations implemented on a NAMD program indicated great stability of chlorcyclizine in protease binding site, suggesting the repurposing of chlorcyclizine as a promising finding in anti-ZIKV drug development.

Proteolytic cleavage of host proteins by the Group IV viral proteases of Venezuelan equine encephalitis virus and Zika virus.

The alphaviral nonstructural protein 2 (nsP2) cysteine proteases (EC 3.4.22.-) are essential for the proteolytic processing of the nonstructural (ns) polyprotein and are validated drug targets. A common secondary role of these proteases is to antagonize the effects of interferon (IFN). After delineating the cleavage site motif of the Venezuelan equine encephalitis virus (VEEV) nsP2 cysteine protease, we searched the human genome to identify host protein substrates. Here we identify a new host substrate of the VEEV nsP2 protease, human TRIM14, a component of the mitochondrial antiviral-signaling protein (MAVS) signalosome. Short stretches of homologous host-pathogen protein sequences (SSHHPS) are present in the nonstructural polyprotein and TRIM14. A 25-residue cyan-yellow fluorescent protein TRIM14 substrate was cleaved in vitro by the VEEV nsP2 protease and the cleavage site was confirmed by tandem mass spectrometry. A TRIM14 cleavage product also was found in VEEV-infected cell lysates. At least ten other Group IV (+)ssRNA viral proteases have been shown to cleave host proteins involved in generating the innate immune responses against viruses, suggesting that the integration of these short host protein sequences into the viral protease cleavage sites may represent an embedded mechanism of IFN antagonism. This interference mechanism shows several parallels with those of CRISPR/Cas9 and RNAi/RISC, but with a protease recognizing a protein sequence common to both the host and pathogen. The short host sequences embedded within the viral genome appear to be analogous to the short phage sequences found in a host's CRISPR spacer sequences. To test this algorithm, we applied it to another Group IV virus, Zika virus (ZIKV), and identified cleavage sites within human SFRP1 (secreted frizzled related protein 1), a retinal Gs alpha subunit, NT5M, and Forkhead box protein G1 (FOXG1) in vitro. Proteolytic cleavage of these proteins suggests a possible link between the protease and the virus-induced phenotype of ZIKV. The algorithm may have value for selecting cell lines and animal models that recapitulate virus-induced phenotypes, predicting host-range and susceptibility, selecting oncolytic viruses, identifying biomarkers, and de-risking live virus vaccines. Inhibitors of the proteases that utilize this mechanism may both inhibit viral replication and alleviate suppression of the innate immune responses.

Molecular dynamics simulation revealed binding of nucleotide inhibitors to ZIKV polymerase over 444 nanoseconds.

In the year 2015, new Zika virus (ZIKV) broke out in Brazil and spread away in more than 80 countries. Scientists directed their efforts toward viral polymerase in attempt to find inhibitors that might interfere with its function. In this study, molecular dynamics simulation (MDS) was performed over 444 ns for a ZIKV polymerase model. Molecular docking (MD) was then performed every 10 ns during the MDS course to ensure the binding of small molecules to the polymerase over the entire time of the simulation. MD revealed the binding ability of four suggested guanosine inhibitors (GIs); (Guanosine substituted with OH and SH (phenyl) oxidanyl in the 2' carbon of the ribose ring). The GIs were compared to guanosine triphosphate (GTP) and five anti-hepatitis C virus drugs (either approved or under clinical trials). The mode of binding and the binding performance of GIs to ZIKV polymerase were found to be the same as GTP. Hence, these compounds were capable of competing GTP for the active site. Moreover, GIs bound to ZIKV active site more tightly compared to ribavirin, the wide-range antiviral drug.

Crystal structure of Zika virus NS5 RNA-dependent RNA polymerase.

The current Zika virus (ZIKV) outbreak became a global health threat of complex epidemiology and devastating neurological impacts, therefore requiring urgent efforts towards the development of novel efficacious and safe antiviral drugs. Due to its central role in RNA viral replication, the non-structural protein 5 (NS5) RNA-dependent RNA-polymerase (RdRp) is a prime target for drug discovery. Here we describe the crystal structure of the recombinant ZIKV NS5 RdRp domain at 1.9 Å resolution as a platform for structure-based drug design strategy. The overall structure is similar to other flaviviral homologues. However, the priming loop target site, which is suitable for non-nucleoside polymerase inhibitor design, shows significant differences in comparison with the dengue virus structures, including a tighter pocket and a modified local charge distribution.

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

America is still suffering with the outbreak of Zika virus (ZIKV) infection. Congenital ZIKV syndrome has already caused a public health emergency of international concern. However, there are still no vaccines to prevent or drugs to treat the infection caused by ZIKV. The ZIKV NS3 helicase (NS3h) protein is a promising target for drug discovery due to its essential role in viral genome replication. NS3h unwinds the viral RNA to enable the replication of the viral genome by the NS5 protein. NS3h contains two important binding sites: the NTPase binding site and the RNA binding site. Here, we used molecular dynamics (MD) simulations to study the molecular behavior of ZIKV NS3h in the presence and absence of ssRNA and the potential implications for NS3h activity and inhibition. Although there is conformational variability and poor electron densities of the RNA binding loop in various apo flaviviruses NS3h crystallographic structures, the MD trajectories of NS3h-ssRNA demonstrated that the RNA binding loop becomes more stable when NS3h is occupied by RNA. Our results suggest that the presence of RNA generates important interactions with the RNA binding loop, and these interactions stabilize the loop sufficiently that it remains in a closed conformation. This closed conformation likely keeps the ssRNA bound to the protein for a sufficient duration to enable the unwinding/replication activities of NS3h to occur. In addition, conformational changes of this RNA binding loop can change the nature and location of the optimal ligand binding site, according to ligand binding site prediction results. These are important findings to help guide the design and discovery of new inhibitors of NS3h as promising compounds to treat the ZIKV infection.

Structure and sequence based functional annotation of Zika virus NS2b protein: Computational insights.

While Zika virus (ZIKV) outbreaks are a growing concern for global health, a deep understanding about the virus is lacking. Here we report a contribution to the basic science on the virus- a detailed computational analysis of the non structural protein NS2b. This protein acts as a cofactor for the NS3 protease (NS3Pro) domain that is important on the viral life cycle, and is an interesting target for drug development. We found that ZIKV NS2b cofactor is highly similar to other virus within the Flavivirus genus, especially to West Nile Virus, suggesting that it is completely necessary for the protease complex activity. Furthermore, the ZIKV NS2b has an important role to the function and stability of the ZIKV NS3 protease domain even when presents a low conservation score. In addition, ZIKV NS2b is mostly rigid, which could imply a non dynamic nature in substrate recognition. Finally, by performing a computational alanine scanning mutagenesis, we found that residues Gly 52 and Asp 83 in the NS2b could be important in substrate recognition.