Structural and functional parameters of the flaviviral protease: a promising antiviral drug target

An Update on the Entomology, Virology, Pathogenesis, and Epidemiology Status of West Nile and Dengue Viruses in Europe (2018-2023)

In recent decades, increases in temperature and tropical rainfall have facilitated the spread of mosquito species into temperate zones. Mosquitoes are vectors for many viruses, including West Nile virus (WNV) and dengue virus (DENV), and pose a serious threat to public health. This review covers most of the current knowledge on the mosquito species associated with the transmission of WNV and DENV and their geographical distribution and discusses the main vertebrate hosts involved in the cycles of WNV or DENV. It also describes virological and pathogenic aspects of WNV or DENV infection, including emerging concepts linking WNV and DENV to the reproductive system. Furthermore, it provides an epidemiological analysis of the human cases of WNV and DENV reported in Europe, from 1 January 2018 to 31 December 2023, with a particular focus on Italy. The first autochthonous cases of DENV infection, with the most likely vector being Aedes albopictus, have been observed in several European countries in recent years, with a high incidence in Italy in 2023. The lack of treatments and effective vaccines is a serious challenge. Currently, the primary strategy to prevent the spread of WNV and DENV infections in humans remains to limit the spread of mosquitoes.

Characterizing the West Nile Virus's polyprotein from nucleotide sequence to protein structure - Computational tools

Objectives

West Nile virus (WNV) belongs to the Flaviviridae family and causes West Nile fever. The mechanism of transmission involves the culex mosquito species. Infected individuals are primarily asymptomatic, and few exhibit common symptoms. Moreover, 10 % of neuronal infection caused by this virus cause death. The proteins encoded by these genes had been uncharacterized, although understanding their function and structure is important for formulating antiviral drugs.

Methods

Herein, we used in silico approaches, including various bioinformatic tools and databases, to analyse the proteins from the WNV polyprotein individually. The characterization included GC content, physicochemical properties, conserved domains, soluble and transmembrane regions, signal localization, protein disorder, and secondary structure features and their respective 3D protein structures.

Results

Among 11 proteins, eight had >50 % GC content, eight proteins had basic pI values, three proteins were unstable under in vitro conditions, four were thermostable according to >100 AI values and some had negative GRAVY values in physicochemical analyses. All protein-conserved domains were shared among Flaviviridae family members. Five proteins were soluble and lacked transmembrane regions. Two proteins had signals for localization in the host endoplasmic reticulum. Non-structural (NS) 2A showed low protein disorder. The secondary structural features and tertiary structure models provide a valuable biochemical resource for designing selective substrates and synthetic inhibitors.

Conclusions

WNV proteins NS2A, NS2B, PM, NS3 and NS5 can be used as drug targets for the pharmacological design of lead antiviral compounds.

Investigating the folding dynamics of NS2B protein of Zika virus

NS2B protein of the Zika virus acts as a co-factor for NS3 protease and also involves in remodeling NS3 protease structure. Therefore, we investigated the overall dynamics of NS2B protein. We find surprising similarities between selected flavivirus NS2B model structures predicted from Alphafold2. Further, the simulated ZIKV NS2B protein structure shows a disordered cytosolic domain (residues 45-95) as a part of a full-length protein. Since only the cytosolic domain of NS2B is sufficient for the protease activity, we also investigated the conformational dynamics of only ZIKV NS2B cytosolic domain (residues 49-95) in the presence of TFE, SDS, Ficoll, and PEG using simulation and spectroscopy. The presence of TFE induces α-helix in NS2B cytosolic domain (residues 49-95). On the other hand, the presence of SDS, ficoll, and PEG does not induce secondary structural change. This dynamics study could have implications for some unknown folds of the NS2B protein.

Allosteric Inhibitors of Zika Virus NS2B-NS3 Protease Targeting Protease in "Super-Open" Conformation

The Zika virus (ZIKV), a member of the Flaviviridae family, is considered a major health threat causing multiple cases of microcephaly in newborns and Guillain-Barré syndrome in adults. In this study, we targeted a transient, deep, and hydrophobic pocket of the "super-open" conformation of ZIKV NS2B-NS3 protease to overcome the limitations of the active site pocket. After virtual docking screening of approximately seven million compounds against the novel allosteric site, we selected the top six candidates and assessed them in enzymatic assays. Six candidates inhibited ZIKV NS2B-NS3 protease proteolytic activity at low micromolar concentrations. These six compounds, targeting the selected protease pocket conserved in ZIKV, serve as unique drug candidates and open new opportunities for possible treatment against several flavivirus infections.

Flavivirus Persistence in Wildlife Populations

A substantial number of humans are at risk for infection by vector-borne flaviviruses, resulting in considerable morbidity and mortality worldwide. These viruses also infect wildlife at a considerable rate, persistently cycling between ticks/mosquitoes and small mammals and reptiles and non-human primates and humans. Substantially increasing evidence of viral persistence in wildlife continues to be reported. In addition to in humans, viral persistence has been shown to establish in mammalian, reptile, arachnid, and mosquito systems, as well as insect cell lines. Although a considerable amount of research has centered on the potential roles of defective virus particles, autophagy and/or apoptosis-induced evasion of the immune response, and the precise mechanism of these features in flavivirus persistence have yet to be elucidated. In this review, we present findings that aid in understanding how vector-borne flavivirus persistence is established in wildlife. Research studies to be discussed include determining the critical roles universal flavivirus non-structural proteins played in flaviviral persistence, the advancement of animal models of viral persistence, and studying host factors that allow vector-borne flavivirus replication without destructive effects on infected cells. These findings underscore the viral–host relationships in wildlife animals and could be used to elucidate the underlying mechanisms responsible for the establishment of viral persistence in these animals.

Antiviral phytochemicals as potent inhibitors against NS3 protease of dengue virus

Dengue, a mosquito-borne disease, has appeared as a major infectious disease globally. The virus requires its proteins to replicate and reproduce in the host cell. The NS3 protease converts the polyprotein to functional proteins with the help of the NS2B cofactor. Thus, NS3 protease is a promising target to develop antiviral inhibitors against the dengue virus. A systematic screening including ADMET properties, molecular docking, molecular dynamics (MD) simulation, binding free energy calculation, and QSAR studies is carried out to predict potent inhibitors against the NS3 protease. From the screening of 40 antiviral phytochemicals, ADMET properties analysis was used to screen out ligands that violate ADME rules and have probable toxicity. Cyanidin 3-Glucoside, Dithymoquinone, and Glabridin were predicted to be potent inhibitors against the NS3 protease according to their binding affinity. These ligands showed several noncovalent interactions, including hydrogen bond, hydrophobic interaction, electrostatic interaction, pi-sulfur interactions. The ligand-protein complexes were further scrutinized using 250 ns molecular dynamics simulation. The MM-PBSA binding free energy calculation was conducted to investigate their binding stability in dynamic conditions. The calculated pIC50(mM) value was predicted using the QSAR model with 89.91% goodness of fit. The predicted biologocal activity value for the ligands indicates they might have good potency.

Experiments and simulation on ZIKV NS2B-NS3 protease reveal its complex folding

Zika virus has been identified in various body fluids such as semen, urine, saliva, cerebrospinal fluid, and vaginal secretion of an infected individual. The pH of these fluids varies from mildly acidic to mildly alkaline. So it is imperative to understand the impact of these conditions on viral protein functioning. We investigated the NS2B-NS3 protease stability and its activity in different denaturing environments. Finding indicates that NS2B-NS3 protease maintains stability at pH 4.8-8.7. Thus it suggests that the complex remains functionally active to hydrolyze the polyprotein within a diverse environmental condition such as variable pH. Despite a stable structure at a broad pH range, a change in environmental conditions dramatically influence its protease activity. Moreover, it is susceptible to structural transformation leading to increased β-strand or helix content in the presence of alcohol. This study may help further to understand the folding-function relationship of the general flaviviral protease complex.

Intermolecular interactions of cn-716 and acyl-KR-aldehyde dipeptide inhibitors against Zika virus

Structural representation and graphic panel showing the most relevant residues that contribute to the ZIKV NS2B–NS3–ligand complexes.

Hydroxychloroquine Inhibits Zika Virus NS2B-NS3 Protease

Zika virus is a mosquito-transmitted flavivirus that causes devastating fetal outcomes in the context of maternal infection during pregnancy. An important target for drugs combatting Zika virus pathogenicity is NS2B-NS3 protease, which plays an essential role in hydrolysis and maturation of the flavivirus polyprotein. We identify hydroxychloroquine, a drug that already has approved uses in pregnancy, as a possible inhibitor of NS2B-NS3 protease by using a Food and Drug Administration-approved drug library, molecular docking, and molecular dynamics simulations. Further, to gain insight into its inhibitory potential toward NS2B-NS3 protease, we performed enzyme kinetic studies, which revealed that hydroxychloroquine inhibits protease activity with an inhibition constant (K i) of 92.34 ± 11.91 μM. Additionally, hydroxychloroquine significantly decreases Zika virus infection in placental cells.

Pharmacophore anchor models of flaviviral NS3 proteases lead to drug repurposing for DENV infection

Background

Viruses of the flaviviridae family are responsible for some of the major infectious viral diseases around the world and there is an urgent need for drug development for these diseases. Most of the virtual screening methods in flaviviral drug discovery suffer from a low hit rate, strain-specific efficacy differences, and susceptibility to resistance. It is because they often fail to capture the key pharmacological features of the target active site critical for protein function inhibition. So in our current work, for the flaviviral NS3 protease, we summarized the pharmacophore features at the protease active site as anchors (subsite-moiety interactions).

Results

For each of the four flaviviral NS3 proteases (i.e., HCV, DENV, WNV, and JEV), the anchors were obtained and summarized into ‘Pharmacophore anchor (PA) models’. To capture the conserved pharmacophore anchors across these proteases, were merged the four PA models. We identified five consensus core anchors (CEH1, CH3, CH7, CV1, CV3) in all PA models, represented as the “Core pharmacophore anchor (CPA) model” and also identified specific anchors unique to the PA models. Our PA/CPA models complied with 89 known NS3 protease inhibitors. Furthermore, we proposed an integrated anchor-based screening method using the anchors from our models for discovering inhibitors. This method was applied on the DENV NS3 protease to screen FDA drugs discovering boceprevir, telaprevir and asunaprevir as promising anti-DENV candidates. Experimental testing against DV2-NGC virus by in-vitro plaque assays showed that asunaprevir and telaprevir inhibited viral replication with EC₅₀ values of 10.4 μM & 24.5 μM respectively. The structure-anchor-activity relationships (SAAR) showed that our PA/CPA model anchors explained the observed in-vitro activities of the candidates. Also, we observed that the CEH1 anchor engagement was critical for the activities of telaprevir and asunaprevir while the extent of inhibitor anchor occupation guided their efficacies.

Conclusion

These results validate our NS3 protease PA/CPA models, anchors and the integrated anchor-based screening method to be useful in inhibitor discovery and lead optimization, thus accelerating flaviviral drug discovery.

Electronic supplementary material

The online version of this article (10.1186/s12859-017-1957-5) contains supplementary material, which is available to authorized users.

STD-NMR experiments identify a structural motif with novel second-site activity against West Nile virus NS2B-NS3 protease

West Nile virus (WNV) belongs to the genus Flavivirus of the family Flaviviridae. This mosquito-borne virus that is highly pathogenic to humans has been evolving into a global threat during the past two decades. Despite many efforts, neither antiviral drugs nor vaccines are available. The viral protease NS2B-NS3^pro is essential for viral replication, and therefore it is considered a prime drug target. However, success in the development of specific NS2B-NS3^pro inhibitors had been moderate so far. In the search for new structural motifs with binding affinity for NS2B-NS3^pro, we have screened a fragment library, the Maybridge Ro5 library, employing saturation transfer difference (STD) NMR experiments as readout. About 30% of 429 fragments showed binding to NS2B-NS3^pro. Subsequent STD-NMR competition experiments using the known active site fragment A as reporter ligand yielded 14 competitively binding fragments, and 22 fragments not competing with A. In a fluorophore-based protease assay, all of these fragments showed inhibition in the micromolar range. Interestingly, 10 of these 22 fragments showed a notable increase of STD intensities in the presence of compound A suggesting cooperative binding. The most promising non-competitive inhibitors 1 and 2 (IC₅₀ ∼ 500 μM) share a structural motif that may guide the development of novel second-site (potentially allosteric) inhibitors of NS2B-NS3^pro. To identify the matching protein binding site, chemical shift perturbation studies employing ¹H,¹⁵N-TROSY-HSQC experiments with uniformly ²H,¹⁵N-labeled protease were performed in the presence of 1, and in the concomitant absence or presence of A. The data suggest that 1 interacts with Met 52* of NS2B, identifying a secondary site adjacent to the binding site of A. Therefore, our study paves the way for the synthesis of novel bidentate NS2B-NS3^pro inhibitors.

Characterization of the Zika virus two-component NS2B-NS3 protease and structure-assisted identification of allosteric small-molecule antagonists

The recent re-emergence of Zika virus (ZIKV)¹, a member of the Flaviviridae family, has become a global emergency. Currently, there are no effective methods of preventing or treating ZIKV infection, which causes severe neuroimmunopathology and is particularly harmful to the developing fetuses of infected pregnant women. However, the pathology induced by ZIKV is unique among flaviviruses, and knowledge of the biology of other family members cannot easily be extrapolated to ZIKV. Thus, structure-function studies of ZIKV proteins are urgently needed to facilitate the development of effective preventative and therapeutic agents. Like other flaviviruses, ZIKV expresses an NS2B-NS3 protease, which consists of the NS2B cofactor and the NS3 protease domain and is essential for cleavage of the ZIKV polyprotein precursor and generation of fully functional viral proteins. Here, we report the enzymatic characterization of ZIKV protease, and we identify structural scaffolds for allosteric small-molecule inhibitors of this protease. Molecular modeling of the protease-inhibitor complexes suggests that these compounds bind to the druggable cavity in the NS2B-NS3 protease interface and affect productive interactions of the protease domain with its cofactor. The most potent compound demonstrated efficient inhibition of ZIKV propagation in vitro in human fetal neural progenitor cells and in vivo in SJL mice. The inhibitory scaffolds could be further developed into valuable research reagents and, ultimately, provide a roadmap for the selection of efficient inhibitors of ZIKV infection.

The drug target genes show higher evolutionary conservation than non-target genes

Although evidence indicates that drug target genes share some common evolutionary features, there have been few studies analyzing evolutionary features of drug targets from an overall level. Therefore, we conducted an analysis which aimed to investigate the evolutionary characteristics of drug target genes. We compared the evolutionary conservation between human drug target genes and non-target genes by combining both the evolutionary features and network topological properties in human protein-protein interaction network. The evolution rate, conservation score and the percentage of orthologous genes of 21 species were included in our study. Meanwhile, four topological features including the average shortest path length, betweenness centrality, clustering coefficient and degree were considered for comparison analysis. Then we got four results as following: compared with non-drug target genes, 1) drug target genes had lower evolutionary rates; 2) drug target genes had higher conservation scores; 3) drug target genes had higher percentages of orthologous genes and 4) drug target genes had a tighter network structure including higher degrees, betweenness centrality, clustering coefficients and lower average shortest path lengths. These results demonstrate that drug target genes are more evolutionarily conserved than non-drug target genes. We hope that our study will provide valuable information for other researchers who are interested in evolutionary conservation of drug targets.

Discovery of Nanomolar Dengue and West Nile Virus Protease Inhibitors Containing a 4-Benzyloxyphenylglycine Residue

The dengue virus (DENV) and West Nile Virus (WNV) NS2B-NS3 proteases are attractive targets for the development of dual-acting therapeutics against these arboviral pathogens. We present the synthesis and extensive biological evaluation of inhibitors that contain benzyl ethers of 4-hydroxyphenylglycine as non-natural peptidic building blocks synthesized via a copper-complex intermediate. A three-step optimization strategy, beginning with fragment growth of the C-terminal 4-hydroxyphenylglycine to the benzyloxy ether, followed by C- and N-terminal optimization, and finally fragment merging generated compounds with in vitro affinities in the low nanomolar range. The most promising derivative reached Ki values of 12 nM at the DENV-2 and 39 nM at the WNV proteases. Several of the newly discovered protease inhibitors yielded a significant reduction of dengue and West Nile virus titers in cell-based assays of virus replication, with an EC50 value of 3.4 μM at DENV-2 and 15.5 μM at WNV for the most active analogue.

Exploring Genomic, Geographic and Virulence Interactions among Epidemic and Non-Epidemic St. Louis Encephalitis Virus (Flavivirus) Strains

St. Louis encephalitis virus (SLEV) is a re-emerging arbovirus in South America. In 2005, an encephalitis outbreak caused by SLEV was reported in Argentina. The reason for the outbreak remains unknown, but may have been related to virological factors, changes in vectors populations, avian amplifying hosts, and/or environmental conditions. The main goal of this study was to characterize the complete genome of epidemic and non-epidemic SLEV strains from Argentina. Seventeen amino acid changes were detected; ten were non-conservative and located in proteins E, NS1, NS3 and NS5. Phylogenetic analysis showed two major clades based on geography: the North America and northern Central America (NAnCA) clade and the South America and southern Central America (SAsCA) clade. Interestingly, the presence of SAsCA genotype V SLEV strains in the NAnCA clade was reported in California, Florida and Texas, overlapping with known bird migration flyways. This work represents the first step in understanding the molecular mechanisms underlying virulence and biological variation among SLEV strains.

West Nile virus drug discovery

The outbreak of West Nile virus (WNV) in 1999 in the USA, and its continued spread throughout the Americas, parts of Europe, the Middle East and Africa, underscored the need for WNV antiviral development. Here, we review the current status of WNV drug discovery. A number of approaches have been used to search for inhibitors of WNV, including viral infection-based screening, enzyme-based screening, structure-based virtual screening, structure-based rationale design, and antibody-based therapy. These efforts have yielded inhibitors of viral or cellular factors that are critical for viral replication. For small molecule inhibitors, no promising preclinical candidate has been developed; most of the inhibitors could not even be advanced to the stage of hit-to-lead optimization due to their poor drug-like properties. However, several inhibitors developed for related members of the family Flaviviridae, such as dengue virus and hepatitis C virus, exhibited cross-inhibition of WNV, suggesting the possibility to re-purpose these antivirals for WNV treatment. Most promisingly, therapeutic antibodies have shown excellent efficacy in mouse model; one of such antibodies has been advanced into clinical trial. The knowledge accumulated during the past fifteen years has provided better rationale for the ongoing WNV and other flavivirus antiviral development.