Characterization of an Allosteric Pocket in Zika Virus NS2B-NS3 Protease

Probing allosteric communication with combined molecular dynamics simulations and network analysis

Understanding the allosteric mechanisms within biomolecules involved in diseases is of paramount importance for drug discovery. Indeed, characterizing communication pathways and critical hotspots in signal transduction can guide a rational approach to leverage allosteric modulation for therapeutic purposes. While the atomistic signatures of allosteric processes are difficult to determine experimentally, computational methods can be a remarkable resource. Network analysis built on Molecular Dynamics simulation data is particularly suited in this respect and is gradually becoming of routine use. Herein, we collect the recent literature in the field, discussing different aspects and available options for network construction and analysis. We further highlight interesting refinements and extensions, eventually providing our perspective on this topic.

Crystallographic fragment screening delivers diverse chemical scaffolds for Zika virus NS2B-NS3 protease inhibitor development

Zika virus (ZIKV) infections cause microcephaly in new-borns and Guillain-Barre syndrome in adults raising a significant global public health concern, yet no vaccines or antiviral drugs have been developed to prevent or treat ZIKV infections. The viral protease NS3 and its co-factor NS2B are essential for the cleavage of the Zika polyprotein precursor into individual structural and non-structural proteins and is therefore an attractive drug target. Generation of a robust crystal system of co-expressed NS2B-NS3 protease has enabled us to perform a crystallographic fragment screening campaign with 1076 fragments. 48 binders with diverse chemical scaffolds were identified in the active site of the protease, with another 6 fragment hits observed in a potential allosteric binding site. Our work provides potential starting points for the development of potent NS2B-NS3 protease inhibitors. Furthermore, we have structurally characterized a potential allosteric binding pocket, identifying opportunities for allosteric inhibitor development.

Aminopyrimidine Derivatives as Multiflavivirus Antiviral Compounds Identified from a Consensus Virtual Screening Approach

Around three billion people are at risk of infection by the dengue virus (DENV) and potentially other flaviviruses. Worldwide outbreaks of DENV, Zika virus (ZIKV), and yellow fever virus (YFV), the lack of antiviral drugs, and limitations on vaccine usage emphasize the need for novel antiviral research. Here, we propose a consensus virtual screening approach to discover potential protease inhibitors (NS3pro) against different flavivirus. We employed an in silico combination of a hologram quantitative structure-activity relationship (HQSAR) model and molecular docking on characterized binding sites followed by molecular dynamics (MD) simulations, which filtered a data set of 7.6 million compounds to 2,775 hits. Lastly, docking and MD simulations selected six final potential NS3pro inhibitors with stable interactions along the simulations. Five compounds had their antiviral activity confirmed against ZIKV, YFV, DENV-2, and DENV-3 (ranging from 4.21 ± 0.14 to 37.51 ± 0.8 μM), displaying aggregator characteristics for enzymatic inhibition against ZIKV NS3pro (ranging from 28 ± 7 to 70 ± 7 μM). Taken together, the compounds identified in this approach may contribute to the design of promising candidates to treat different flavivirus infections.

Structural Insights into Plasticity and Discovery of Flavonoid Allosteric Inhibitors of Flavivirus NS2B–NS3 Protease

Flaviviruses are among the most critical pathogens in tropical regions; they cause various severe diseases in developing countries but are not restricted to these countries. The development of antiviral therapeutics is crucial for managing flavivirus outbreaks. Ten proteins are encoded in the flavivirus RNA. The N2B–NS3pro protein complex plays a fundamental role in flavivirus replication and is a promising drug target; however, no flavivirus protease inhibitors have progressed to the preclinical stage. This study analyzed the structural models and plasticity of the NS2B–NS3pro protein complex of five medically important non-dengue flaviviruses (West Nile, Rocio, Ilhéus, yellow fever, and Saint Louis encephalitis). The flavonoids amentoflavone, tetrahydrorobustaflavone, and quercetin were selected for their exceptional binding energies as potential inhibitors of the NS2B–NS3pro protein complex. AutoDock Vina results ranged from −7.0 kcal/mol to −11.5 kcal/mol and the compounds preferentially acted non-competitively. Additionally, the first structural model for the NS2B–NS3pro protein complex was proposed for Ilhéus and Rocio viruses. The NS2B–NS3pro protease is an attractive molecular target for drug development. The three identified natural flavonoids showed great inhibitory potential against the viral species. Nevertheless, further in silico and in vitro studies are required to obtain more information regarding NS2B–NS3pro inhibition by these flavonoids and their therapeutic potential.

Advances in Computational Methods to Discover New NS2B-NS3 Inhibitors Useful Against Dengue and Zika Viruses

The Flaviviridae virus family consists of the genera Hepacivirus, Pestivirus, and Flavivirus, with approximately 70 viral types that use arthropods as vectors. Among these diseases, dengue (DENV) and zika virus (ZIKV) serotypes stand out, responsible for thousands of deaths worldwide. Due to the significant increase in cases, the World Health Organization (WHO) declared DENV a potential threat for 2019 due to being transmitted by infected travelers. Furthermore, ZIKV also has a high rate of transmissibility, highlighted in the outbreak in 2015, generating consequences such as Guillain-Barré syndrome and microcephaly. According to clinical outcomes, those infected with DENV can be asymptomatic, and in other cases, it can be lethal. On the other hand, ZIKV has severe neurological symptoms in newborn babies and adults. More serious symptoms include microcephaly, brain calcifications, intrauterine growth restriction, and fetal death. Despite these worrying data, no drug or vaccine is approved to treat these diseases. In the drug discovery process, one of the targets explored against these diseases is the NS2B-NS3 complex, which presents the catalytic triad His51, Asp75, and Ser135, with the function of cleaving polyproteins, with specificity for basic amino acid residues, Lys- Arg, Arg-Arg, Arg-Lys or Gln-Arg. Since NS3 is highly conserved in all DENV serotypes and plays a vital role in viral replication, this complex is an excellent drug target. In recent years, computer-aided drug discovery (CADD) is increasingly essential in drug discovery campaigns, making the process faster and more cost-effective, mainly explained by discovering new drugs against DENV and ZIKV. Finally, the main advances in computational methods applied to discover new compounds against these diseases will be presented here. In fact, molecular dynamics simulations and virtual screening is the most explored approach, providing several hit and lead compounds that can be used in further optimizations. In addition, fragment-based drug design and quantum chemistry/molecular mechanics (QM/MM) provides new insights for developing anti-DENV/ZIKV drugs. We hope that this review offers further helpful information for researchers worldwide and stimulates the use of computational methods to find a promising drug for treating DENV and ZIKV.

Functional Loop Dynamics and Characterization of the Inactive State of the NS2B-NS3 Dengue Protease due to Allosteric Inhibitor Binding

Dengue virus, a flavivirus that causes dengue shock syndrome and dengue hemorrhagic fever, is currently prevalent worldwide. A two-component protease (NS2B-NS3) is essential for maturation, representing an important target for designing anti-flavivirus drugs. Previously, consideration has been centered on developing active-site inhibitors of NS2B-NS3pro. However, the flat and charged nature of its active site renders difficulties in developing inhibitors, suggesting an alternative strategy for identifying allosteric inhibitors. The allosterically sensitive site of the dengue protease is located near Ala125, between the 120s loop and 150s loop. Using atomistic molecular dynamics simulations, we have explored the protease's conformational dynamics upon binding of an allosteric inhibitor. Furthermore, characterization of the inherent flexible loops (71-75s loop, 120s loop, and 150s loop) is carried out for allosteric-inhibitor-bound wild-type and mutant A125C variants and a comparison is performed with its unbound state to extract the structural changes describing the inactive state of the protease. Our study reveals that compared to the unliganded system, the inhibitor-bound system shows large structural changes in the 120s loop and 150s loop in contrast to the rigid 71-75s loop. The unliganded system shows a closed-state pocket in contrast to the open state for the wild-type complex that locks the protease into the open and inactive-state conformations. However, the mutant complex fluctuates between open and closed states. Also, we tried to see how mutation and binding of an allosteric inhibitor perturb the connectivity in a protein structure network (PSN) at contact levels. Altogether, our study reveals the mechanism of conformational rearrangements of loops at the molecular level, locking the protein in an inactive conformation, which may be useful for developing allosteric inhibitors.