Serotype-specific interactions among functional domains of dengue virus 2 nonstructural proteins (NS) 5 and NS3 are crucial for viral RNA replication

Predicting repurposed drugs targeting the NS3 protease of dengue virus using machine learning-based QSAR, molecular docking, and molecular dynamics simulations

Dengue fever, prevalent in Southeast Asian countries, currently lacks effective pharmaceutical interventions for virus replication control. This study employs a strategy that combines machine learning (ML)-based quantitative-structure-activity relationship (QSAR), molecular docking, and molecular dynamics simulations to discover potential inhibitors of the NS3 protease of the dengue virus. We used nine molecular fingerprints from PaDEL to extract features from the NS3 protease dataset of dengue virus type 2 in the ChEMBL database. Feature selection was achieved through the low variance threshold, F-Score, and recursive feature elimination (RFE) methods. Our investigation employed three ML models - support vector machine (SVM), random forest (RF), and extreme gradient boosting (XGBoost) - for classifier development. Our SVM model, combined with SVM-RFE, had the best accuracy (0.866) and ROC_AUC (0.964) in the testing set. We identified potent inhibitors on the basis of the optimal classifier probabilities and docking binding affinities. SHAP and LIME analyses highlighted the significant molecular fingerprints (e.g. ExtFP69, ExtFP362, ExtFP576) involved in NS3 protease inhibitory activity. Molecular dynamics simulations indicated that amphotericin B exhibited the highest binding energy of -212 kJ/mol and formed a hydrogen bond with the critical residue Ser196. This approach enhances NS3 protease inhibitor identification and expedites the discovery of dengue therapeutics.

Dengue virus serotypic replacement of NS3 protease or helicase domain causes chimeric viral attenuation but can be recovered by a compensated mutation at helicase domain or NS2B, respectively

Mosquito-borne dengue viruses (DENVs) have evolved to four serotypes with 69%-78% amino acid identities, resulting in incomplete immunity, where one serotype's infection does not cross-protect against secondary infections by other serotypes. Despite the amino acid differences, structural and nonstructural (NS) proteins among serotypes play similar functions. NS3 is an enzyme complex: NS3 has N-terminal protease (PRO) and C-terminal helicase (HEL) activities in addition to 5' RNA triphosphatase (5'RTP), which is involved in the RNA capping process. In this study, the effects of NS3 replacements among serotypes were tested. The replacement of NS3 full-length (FULL), PRO or HEL region suppressed viral replication in BHK-21 mammalian cells, while the single compensatory mutation improved the viral replications; P364S mutation in HEL revived PRO (DENV3)-replaced DENV1, while S68T alteration in NS2B recovered HEL (DENV1)-replaced DENV2. The results suggest that the interactions between PRO and HEL as well as HEL and NS2B are required for replication competence. Lower-frequency mutations also appeared at various locations in viral proteins, although after infecting C6/36 mosquito cells, the mutations' frequencies changed, and/or new mutations appeared. In contrast, the inter-domain region (INT, 12 amino acids)-replaced chimera quickly replicated without mutation in BHK-21 cells, although extended cell culture accumulated various mutations. These results suggest that NS3 variously interacts with DENV proteins, in which the chimeric NS3 domain replacements induced amino acid mutations, irrespective of replication efficiency. However, the viral sequences are further adjusted for replication efficiency, to fit in both mammalian cells and mosquito cells. IMPORTANCE Enzyme activities for replicating DENV 5' cap positive (+) sense RNA have been shown to reside in NS3 and NS5. However, it remains unknown how these enzymes coordinately synthesize negative (-) sense RNA, from which abundant 5' cap (+) sense RNA is produced. We previously revealed that NS5 dimerization and NS5 methyltransferase(MT)-NS3HEL interaction are important for DENV replication. Here, we found that replication incompetence due to NS3PRO or HEL replacement was compensated by a mutation at HEL or NS2B, respectively, suggesting that the interactions among NS2B, NS3PRO, and HEL are critical for DENV replication.

Flavivirus proteases: The viral Achilles heel to prevent future pandemics

Flaviviruses are important human pathogens and include dengue (DENV), West Nile (WNV), Yellow fever virus (YFV), Japanese encephalitis (JEV) and Zika virus (ZIKV). DENV, transmitted by mosquitoes, causes diseases ranging in severity from mild dengue fever with non-specific flu-like symptoms to fatal dengue hemorrhagic fever and dengue shock syndrome. DENV infections are caused by four serotypes, DENV1-4, which interact differently with antibodies in blood serum. The incidence of DENV infection has increased dramatically in recent decades and the CDC estimates 400 million dengue infections occur each year, resulting in ∼25,000 deaths mostly among children and elderly people. Similarly, ZIKV infections are caused by infected mosquito bites to humans, can be transmitted sexually and through blood transfusions. If a pregnant woman is infected, the virus can cross the placental barrier and can spread to her fetus, causing severe brain malformations in the child including microcephaly and other birth defects. It is noteworthy that the neurological manifestations of ZIKV were also observed in DENV endemic regions, suggesting that pre-existing antibody response to DENV could augment ZIKV infection. WNV, previously unknown in the US (and known to cause only mild disease in Middle East), first arrived in New York city in 1999 (NY99) and spread throughout the US and Canada by Culex mosquitoes and birds. WNV is now endemic in North America. Thus, emerging and re-emerging flaviviruses are significant threat to human health. However, vaccines are available for only a limited number of flaviviruses, and antiviral therapies are not available for any flavivirus. Hence, there is an urgent need to develop therapeutics that interfere with essential enzymatic steps, such as protease in the flavivirus lifecycle as these viruses possess significant threat to future pandemics. In this review, we focus on our E. coli expression of NS2B hydrophilic domain (NS2BH) covalently linked to NS3 protease domain (NS3Pro) in their natural context which is processed by the combined action of both subunits of the NS2B-NS3Pro precursor. Biochemical activities of the viral protease such as solubility and autoproteolysis of NS2BH-NS3Pro linkage depended on the C-terminal portion of NS2BH linked to the NS3Pro domain. Since 2008, we also focus on the use of the recombinant protease in high throughput screens and characterization of small molecular compounds identified in these screens.

In Silico Prediction of the Phosphorylation of NS3 as an Essential Mechanism for Dengue Virus Replication and the Antiviral Activity of Quercetin

Simple Summary

Dengue is a mosquito-borne virus that infects up to 400 million people worldwide annually. Dengue infection triggers high fever, severe body aches, rash, low platelet count, and could lead to Dengue hemorrhagic fever (DHF) in some cases. There is currently no cure, nor a broadly effective vaccine. The interaction of two viral proteins, nonstructural Proteins 3 and 5 (NS3 and NS5), is required for viral replication in the infected host’s cells. Our computational modeling of NS3 suggested that phosphorylation of a serine residue at position 137 of NS3 by a specific c-Jun N-terminal kinase (JNK) enhances viral replication by increasing the interaction of NS3 and NS5 through structural changes in amino acid residues 49–95. Experimental studies have shown that inhibition of JNK prevents viral replication and have suggested that the plants’ flavonoid Quercetin, Agathis flavone, and Myricetin inhibit Dengue infection. Our molecular simulations revealed that Quercetin binds NS3 and obstructs serine 137 phosphorylation, which may decrease viral replication. This work offers a molecular mechanism that can be used for anti-Dengue drug development.

Abstract

Dengue virus infection is a global health problem for which there have been challenges to obtaining a cure. Current vaccines and anti-viral drugs can only be narrowly applied in ongoing clinical trials. We employed computational methods based on structure-function relationships between human host kinases and viral nonstructural protein 3 (NS3) to understand viral replication inhibitors’ therapeutic effect. Phosphorylation at each of the two most evolutionarily conserved sites of NS3, serine 137 and threonine 189, compared to the unphosphorylated state were studied with molecular dynamics and docking simulations. The simulations suggested that phosphorylation at serine 137 caused a more remarkable structural change than phosphorylation at threonine 189, specifically located at amino acid residues 49–95. Docking studies supported the idea that phosphorylation at serine 137 increased the binding affinity between NS3 and nonstructural Protein 5 (NS5), whereas phosphorylation at threonine 189 decreased it. The interaction between NS3 and NS5 is essential for viral replication. Docking studies with the antiviral plant flavonoid Quercetin with NS3 indicated that Quercetin physically occluded the serine 137 phosphorylation site. Taken together, these findings suggested a specific site and mechanism by which Quercetin inhibits dengue and possible other flaviviruses.

Cross-serotypically conserved epitope recommendations for a universal T cell-based dengue vaccine

Dengue virus (DENV)-associated disease is a growing threat to public health across the globe. Co-circulating as four different serotypes, DENV poses a unique challenge for vaccine design as immunity to one serotype predisposes a person to severe and potentially lethal disease upon infection from other serotypes. Recent experimental studies suggest that an effective vaccine against DENV should elicit a strong T cell response against all serotypes, which could be achieved by directing T cell responses toward cross-serotypically conserved epitopes while avoiding serotype-specific ones. Here, we used experimentally-determined DENV T cell epitopes and patient-derived DENV sequences to assess the cross-serotypic variability of the epitopes. We reveal a distinct near-binary pattern of epitope conservation across serotypes for a large number of DENV epitopes. Based on the conservation profile, we identify a set of 55 epitopes that are highly conserved in at least 3 serotypes. Most of the highly conserved epitopes lie in functionally important regions of DENV non-structural proteins. By considering the global distribution of human leukocyte antigen (HLA) alleles associated with these DENV epitopes, we identify a potentially robust subset of HLA class I and class II restricted epitopes that can serve as targets for a universal T cell-based vaccine against DENV while covering ~99% of the global population.

Positive epistasis between viral polymerase and the 3' untranslated region of its genome reveals the epidemiologic fitness of dengue virus

Dengue virus (DENV) is a global health threat, causing repeated epidemics throughout the tropical world. While low herd immunity levels to any one of the four antigenic types of DENV predispose populations to outbreaks, viral genetic determinants that confer greater fitness for epidemic spread is an important but poorly understood contributor of dengue outbreaks. Here we report that positive epistasis between the coding and noncoding regions of the viral genome combined to elicit an epidemiologic fitness phenotype associated with the 1994 DENV2 outbreak in Puerto Rico. We found that five amino acid substitutions in the NS5 protein reduced viral genomic RNA (gRNA) replication rate to achieve a more favorable and relatively more abundant subgenomic flavivirus RNA (sfRNA), a byproduct of host 5'-3' exoribonuclease activity. The resulting increase in sfRNA relative to gRNA levels not only inhibited type I interferon (IFN) expression in infected cells through a previously described mechanism, but also enabled sfRNA to compete with gRNA for packaging into infectious particles. We suggest that delivery of sfRNA to new susceptible cells to inhibit type I IFN induction before gRNA replication and without the need for further de novo sfRNA synthesis could form a "preemptive strike" strategy against DENV.

Antiviral activity of the natural alkaloid anisomycin against dengue and Zika viruses

Flaviviruses constitute a public health concern because of their global burden and the lack of specific antiviral treatment. Here we investigated the antiviral activity of the alkaloid anisomycin against dengue (DENV) and Zika (ZIKV) viruses. At non-cytotoxic concentrations, anisomycin strongly inhibited the replication of reference strains and clinical isolates of all DENV serotypes and Asian and African strains of ZIKV in Vero cells. Anisomycin also prevented DENV and ZIKV multiplication in human cell lines. While initial steps of DENV and ZIKV replicative cycle were unaffected, a high inhibition of viral protein expression was demonstrated after treatment with anisomycin. DENV RNA synthesis was strongly reduced in anisomycin treated cultures, but the compound did not exert a direct inhibitory effect on 2' O-methyltransferase or RNA polymerase activities of DENV NS5 protein. Furthermore, anisomycin-mediated activation of p38 signaling was not related to the antiviral action of the compound. The evaluation of anisomycin efficacy in a mouse model of ZIKV morbidity and mortality revealed that animals treated with a low dose of anisomycin exhibited a significant reduction in viremia levels and died significantly later than the control group. This protective effect was lost at higher doses, though. In conclusion, anisomycin is a potent and selective in vitro inhibitor of DENV and ZIKV that impairs a post-entry step of viral replication; and a low-dose anisomycin treatment may provide some minimal benefit in a mouse model.

Flavivirus RNA-Dependent RNA Polymerase Interacts with Genome UTRs and Viral Proteins to Facilitate Flavivirus RNA Replication

Flaviviruses, most of which are emerging and re-emerging human pathogens and significant public health concerns worldwide, are positive-sense RNA viruses. Flavivirus replication occurs on the ER and is regulated by many mechanisms and factors. NS5, which consists of a C-terminal RdRp domain and an N-terminal methyltransferase domain, plays a pivotal role in genome replication and capping. The C-terminal RdRp domain acts as the polymerase for RNA synthesis and cooperates with diverse viral proteins to facilitate productive RNA proliferation within the replication complex. Here, we provide an overview of the current knowledge of the functions and characteristics of the RdRp, including the subcellular localization of NS5, as well as the network of interactions formed between the RdRp and genome UTRs, NS3, and the methyltransferase domain. We posit that a detailed understanding of RdRp functions may provide a target for antiviral drug discovery and therapeutics.

Zika virus NS3 is a canonical RNA helicase stimulated by NS5 RNA polymerase

Zika virus is a positive single-strand RNA virus whose replication involved RNA unwinding and synthesis. ZIKV NS3 contains a helicase domain, but its enzymatic activity is not fully characterized. Here, we established a dsRNA unwinding assay based on the FRET effect to study the helicase activity of ZIKV NS3, which provided kinetic information in real time. We found that ZIKV NS3 specifically unwound dsRNA/dsDNA with a 3' overhang in the 3' to 5' direction. The RNA unwinding ability of NS3 significantly decreased when the duplex was longer than 18 base pairs. The helicase activity of NS3 depends on ATP hydrolysis and binding to RNA. Mutations in the ATP binding region or the RNA binding region of NS3 impair its helicase activity, thus blocking viral replication in the cell. Furthermore, we showed that ZIKV NS5 interacted with NS3 and stimulated its helicase activity. Disrupting NS3-NS5 interaction resulted in a defect in viral replication, revealing the tight coupling of RNA unwinding and synthesis. We suggest that NS3 helicase activity is stimulated by NS5; thus, viral replication can be carried out efficiently. Our work provides a molecular mechanism of ZIKV NS3 unwinding and novel insights into ZIKV replication.

Infection of Aedes albopictus Mosquito C6/36 Cells with the wMelpop Strain of Wolbachia Modulates Dengue Virus-Induced Host Cellular Transcripts and Induces Critical Sequence Alterations in the Dengue Viral Genome

Dengue virus (DENV) causes frequent epidemics infecting ∼390 million people annually in over 100 countries. There are no approved vaccines or antiviral drugs for treatment of infected patients. However, there is a novel approach to control DENV transmission by the mosquito vectors, Aedes aegypti and Aedes albopictus, using the Wolbachia symbiont. The wMelPop strain of Wolbachia suppresses DENV transmission and shortens the mosquito life span. However, the underlying mechanism is poorly understood. To clarify this mechanism, either naive A. albopictus (C6/36) or wMelPop-C6/36 cells were infected with DENV serotype 2 (DENV2). Analysis of host transcript profiles by transcriptome sequencing (RNAseq) revealed that the presence of wMelPop dramatically altered the mosquito host cell transcription in response to DENV2 infection. The viral RNA evolved from wMelPop-C6/36 cells contained low-frequency mutations (∼25%) within the coding region of transmembrane domain 1 (TMD1) of E protein. Mutations with >97% frequencies were distributed within other regions of E, the NS5 RNA-dependent RNA polymerase (NS5POL) domain, and the TMDs of NS2A, NS2B, and NS4B. Moreover, while DENV2-infected naive C6/36 cells showed syncytium formation, DENV2-infected wMelPop-C6/36 cells did not. The Wolbachia-induced mutant DENV2 can readily infect and replicate in naive C6/36 cells, whereas in mutant DENV2-infected BHK-21 or Vero cells, virus replication was delayed. In LLC-MK2 cells, the mutant failed to produce plaques. Additionally, in BHK-21 cells, many mutations in the viral genome reverted to the wild type (WT) and compensatory mutations in NS3 gene appeared. Our results indicate that wMelPop impacts significantly the interactions of DENV2 with mosquito and mammalian host cells.IMPORTANCE Mosquito-borne diseases are of global significance causing considerable morbidity and mortality throughout the world. Dengue virus (DENV; serotypes 1 to 4), a member of the Flavivirus genus of the Flaviviridae family, causes millions of infections annually. Development of a safe vaccine is hampered due to absence of cross-protection and increased risk in secondary infections due to antibody-mediated immune enhancement. Infection of vector mosquitoes with Wolbachia bacteria offers a novel countermeasure to suppress DENV transmission, but the mechanisms are poorly understood. In this study, the host transcription profiles and viral RNA sequences were analyzed in naive A. albopictus (C6/36) and wMelPop-C6/36 cells by RNAseq. Our results showed that the wMelPop symbiont caused profound changes in host transcription profiles and morphology of DENV2-infected C6/36 cells. Accumulation of several mutations throughout DENV2 RNA resulted in loss of infectivity of progeny virions. Our findings offer new insights into the mechanism of Wolbachia-mediated suppression of DENV transmission.

Comparative Analysis of NS5 Protein for Tick Borne Encephalitis Virus Strains in three Virus Subtypes

Non-structural protein 5 (NS5) of tick-borne encephalitis virus is an enzyme which is responsible for a copying of viral RNA, and it has a strong structural similarity to RNA polymerases of another RNA virus families. The strains of the virus are separated into three subtypes, which differ by specific mutations in virus proteins, including NS5 protein. The methods of structural bioinformatics allow to construct a model of NS5 protein for several strains of the virus.

The paper presents the comparative analysis of sequences and structures of NS5 protein, for three subtypes of the tick-borne encephalitis virus. The segments of protein were identified where the highest difference between subtypes and within subtypes is observed. These segments, where most of the mutations are accumulated, are located in methyltransferase domain, in the inter-domain interface, and in the three subdomains of polymerase domain. The association between the locations of mutations in NS5 protein and the flexibility of a protein backbone was observed using normal mode analysis. Namely, the most important mutations are located in the parts of protein where the amplitude of synchronous oscillations estimated using normal mode analysis is the highest: in the second zinc binding pocket within polymerase domain, in the N-terminal extension within inter-domain interface, and around an active site of methyltransferase domain.

Genomic and structural features of the yellow fever virus from the 2016-2017 Brazilian outbreak

Southeastern Brazil has been suffering a rapid expansion of a severe sylvatic yellow fever virus (YFV) outbreak since late 2016, which has reached one of the most populated zones in Brazil and South America, heretofore a yellow fever-free zone for more than 70 years. In the current study, we describe the complete genome of 12 YFV samples from mosquitoes, humans and non-human primates from the Brazilian 2017 epidemic. All of the YFV sequences belong to the modern lineage (sub-lineage 1E) of South American genotype I, having been circulating for several months prior to the December 2016 detection. Our data confirm that viral strains associated with the most severe YF epidemic in South America in the last 70 years display unique amino acid substitutions that are mainly located in highly conserved positions in non-structural proteins. Our data also corroborate that YFV has spread southward into Rio de Janeiro state following two main sylvatic dispersion routes that converged at the border of the great metropolitan area comprising nearly 12 million unvaccinated inhabitants. Our original results can help public health authorities to guide the surveillance, prophylaxis and control measures required to face such a severe epidemiological problem. Finally, it will also inspire other workers to further investigate the epidemiological and biological significance of the amino acid polymorphisms detected in the Brazilian 2017 YFV strains.

The rearrangement of motif F in the flavivirus RNA-directed RNA polymerase

In the flavivirus genus, the non-structural protein NS5 plays a central role in RNA viral replication and constitutes a major target for drug discovery. One of the prime challenges in the study of NS5 protein is to investigate the interplay between the two protein domains, namely, the RNA-dependent RNA polymerase (RdRp) domain and the methyltransferase (MTase) domain. These investigations could clarify the multiple roles of NS5 protein in the virus life cycle. Here we present the results of sequence analyses and structural bioinformatics studies of NS5 protein, which suggest that the conserved motif F in the NS5 protein could act as a lock which controls the rearrangement of the domains and as a switch in the protein enzymatic activity.