Conformational selection in the flaviviral NS2B-NS3 protease

Alternate recognition by dengue protease: Proteolytic and binding assays provide functional evidence beyond an induced-fit

Proteases are key enzymes in viral replication, and interfering with these targets is the basis for therapeutic interventions. We previously introduced a hypothesis about conformational selection in the protease of dengue virus and related flaviviruses, based on conformational plasticity noted in X-ray structures. The present work presents the first functional evidence for alternate recognition by the dengue protease, in a mechanism based primarily on conformational selection rather than induced-fit. Recognition of distinct substrates and inhibitors in proteolytic and binding assays varies to a different extent, depending on factors reported to influence the protease structure. The pH, salinity, buffer type, and temperature cause a change in binding, proteolysis, or inhibition behavior. Using representative inhibitors with distinct structural scaffolds, we identify two contrasting binding profiles to dengue protease. Noticeable effects are observed in the binding assay upon inclusion of a non-ionic detergent in comparison to the proteolytic assay. The findings highlight the impact of the selection of testing conditions on the observed ligand affinity or inhibitory potency. From a broader scope, the dengue protease presents an example, where the induced-fit paradigm appears insufficient to explain binding events with the biological target. Furthermore, this protein reveals the complexity of comparing or combining biochemical assay data obtained under different conditions. This can be particularly critical for artificial intelligence (AI) approaches in drug discovery that rely on large datasets of compounds activity, compiled from different sources using non-identical testing procedures. In such cases, mismatched results will compromise the model quality and its predictive power.

Combined NMR and molecular dynamics conformational filter identifies unambiguously dynamic ensembles of Dengue protease NS2B/NS3pro

The dengue protease NS2B/NS3pro has been reported to adopt either an 'open' or a 'closed' conformation. We have developed a conformational filter that combines NMR with MD simulations to identify conformational ensembles that dominate in solution. Experimental values derived from relaxation parameters for the backbone and methyl side chains were compared with the corresponding back-calculated relaxation parameters of different conformational ensembles obtained from free MD simulations. Our results demonstrate a high prevalence for the 'closed' conformational ensemble while the 'open' conformation is absent, indicating that the latter conformation is most probably due to crystal contacts. Conversely, conformational ensembles in which the positioning of the co-factor NS2B results in a 'partially' open conformation, previously described in both MD simulations and X-ray studies, were identified by our conformational filter. Altogether, we believe that our approach allows for unambiguous identification of true conformational ensembles, an essential step for reliable drug discovery.

Designing photoaffinity tool compounds for the investigation of the DENV NS2B-NS3 protease allosteric binding pocket

Dengue virus (DENV) infection still lacks specific antiviral therapy, making the NS2B-NS3 protease an attractive target for drug development. However, allosteric inhibitors that bind to a site other than the active site still need to be better understood. In this study, we designed and synthesised tool compounds for photoaffinity labelling (PAL) to investigate the binding site of allosteric inhibitors on the DENV protease. These tool compounds contained an affinity moiety, a photoreactive group, and a reporter tag for detection. Upon irradiation, the photoreactive group formed a covalent bond with the protease, allowing for binding site identification. SDS-PAGE-based assays confirmed the qualitative binding of the designed inhibitors to the allosteric pocket, and pull-down experiments validated the interaction. Tryptic protein digestion following liquid chromatography/mass spectrometry analysis further supported the binding of the inhibitor to the proposed pocket revealing photo-attachment to an NS3 loop close to the C-terminus. These results enhance our understanding of allosteric inhibitors and their mechanism of action against the DENV protease. The developed tool compounds and PAL are potent tools for future drug discovery efforts and investigations targeting the DENV protease.

Boophilin D1, a Kunitz type protease inhibitor, as a source of inhibitors for the ZIKA virus NS2B-NS3 protease

Arboviruses are a global concern for a multitude of reasons, including their increased incidence and human mortality. Vectors associated with arboviruses include the mosquito Aedes sp., which is responsible for transmitting the Zika virus. Flaviviruses, like the Zika virus, present only one chymotrypsin-like serine protease (NS3) in their genome. Together with host enzymes, the NS2B co-factor NS3 protease complex are essential for the viral replication cycle by virus polyprotein processing. To search for Zika virus NS2B-NS3 protease (ZIKVPro) inhibitors, a phage display library was constructed using the Boophilin domain 1 (BoophD1), a thrombin inhibitor from the Kunitz family. A BoophilinD1 library mutated at positions P1-P4' was constructed, presenting a titer of 2.9x106 (cfu), and screened utilizing purified ZIKVPro. The results demonstrated at the P1-P4' positions the occurrence of 47% RALHA sequence (mut 12) and 11.8% RASWA sequence (mut14), SMRPT, or KALIP (wt) sequence. BoophD1-wt and mutants 12 and 14 were expressed and purified. The purified BoophD1 wt, mut 12 and 14, presented Ki values for ZIKVPro of 0.103, 0.116, and 0.101 μM, respectively. The BoophD1 mutant inhibitors inhibit the Dengue virus 2 protease (DENV2) with Ki values of 0.298, 0.271, and 0.379 μM, respectively. In conclusion, BoophD1 mut 12 and 14 selected for ZIKVPro demonstrated inhibitory activity like BoophD1-wt, suggesting that these are the strongest Zika inhibitors present in the BoophD1 mutated phage display library. Furthermore, BoophD1 mutants selected for ZIKVPro inhibit both Zika and Dengue 2 proteases making them potential pan-flavivirus inhibitors.

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.

Therapeutics for flaviviral infections

Flaviviruses are vector-borne pathogens capable of causing devastating human diseases. The re-emergence of Zika in 2016 notoriously led to a widescale epidemic in the Americas. New daunting evidence suggests that a single mutation in Zika virus genome may increase transmission and pathogenesis, further highlighting the need to be prepared for flavivirus outbreaks. Dengue, in particular infects about 400 million people each year, leading to reoccurring local outbreaks. Public health efforts to mitigate flavivirus transmission is largely dependent on vector control strategies, as only a limited number of flavivirus vaccines have been developed thus far. There are currently no commercially available antivirals for flaviviruses, leaving supportive care as the primary treatment option. In this review, we will briefly paint a broad picture of the flavivirus landscape in terms of therapeutics, with particular focus on viral targets, promising novel compounds entering the drug discovery pipeline, as well as model systems for evaluating drug efficacy.

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.

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.

Discovery of potent benzoxaborole inhibitors against SARS-CoV-2 main and dengue virus proteases

The RNA viruses SARS-CoV-2 and dengue pose a major threat to human health worldwide and their proteases (M^pro; NS2B/NS3) are considered as promising targets for drug development. We present the synthesis and biological evaluation of novel benzoxaborole inhibitors of these two proteases. The most active compound achieves single-digit micromolar activity against SARS-CoV-2 M^pro in a biochemical assay. The most active substance against dengue NS2B/NS3 protease has submicromolar activity in cells (EC₅₀ 0.54 μM) and inhibits DENV-2 replication in cell culture. Most benzoxaboroles had no relevant cytotoxicity or significant off-target inhibition. Furthermore, the class demonstrated passive membrane penetration and stability against the evaluated proteases. This compound class may contribute to the development of antiviral agents with activity against DENV or SARS-CoV-2.

1H, 13C and 15N resonance assignment of backbone and IVL-methyl side chain of the S135A mutant NS3pro/NS2B protein of Dengue II virus reveals unique secondary structure features in solution

The serotype II Dengue (DENV 2) virus is the most prevalent of all four known serotypes. Herein, we present nearly complete ¹H, ¹⁵N, and ¹³C backbone and ¹H, ¹³C isoleucine, valine, and leucine methyl resonance assignment of the apo S135A catalytically inactive variant of the DENV 2 protease enzyme folded as a tandem formed between the serine protease domain NS3pro and the cofactor NS2B, as well as the secondary structure prediction of this complex based on the assigned chemical shifts using the TALOS-N software. Our results provide a solid ground for future elucidation of the structure and dynamic of the apo NS3pro/NS2B complex, key for adequate development of inhibitors, and a thorough molecular understanding of their function(s).

Precursors of Viral Proteases as Distinct Drug Targets

Viral proteases are indispensable for successful virion maturation, thus making them a prominent drug target. Their enzyme activity is tightly spatiotemporally regulated by expression in the precursor form with little or no activity, followed by activation via autoprocessing. These cleavage events are frequently triggered upon transportation to a specific compartment inside the host cell. Typically, precursor oligomerization or the presence of a co-factor is needed for activation. A detailed understanding of these mechanisms will allow ligands with non-canonical mechanisms of action to be designed, which would specifically modulate the initial irreversible steps of viral protease autoactivation. Binding sites exclusive to the precursor, including binding sites beyond the protease domain, can be exploited. Both inhibition and up-regulation of the proteolytic activity of viral proteases can be detrimental for the virus. All these possibilities are discussed using examples of medically relevant viruses including herpesviruses, adenoviruses, retroviruses, picornaviruses, caliciviruses, togaviruses, flaviviruses, and coronaviruses.

Structure and Dynamics of Zika Virus Protease and Its Insights into Inhibitor Design

Zika virus (ZIKV)—a member of the Flaviviridae family—is an important human pathogen. Its genome encodes a polyprotein that can be further processed into structural and non-structural proteins. ZIKV protease is an important target for antiviral development due to its role in cleaving the polyprotein to release functional viral proteins. The viral protease is a two-component protein complex formed by NS2B and NS3. Structural studies using different approaches demonstrate that conformational changes exist in the protease. The structures and dynamics of this protease in the absence and presence of inhibitors were explored to provide insights into the inhibitor design. The dynamic nature of residues binding to the enzyme cleavage site might be important for the function of the protease. Due to the charges at the protease cleavage site, it is challenging to develop small-molecule compounds acting as substrate competitors. Developing small-molecule compounds to inhibit protease activity through an allosteric mechanism is a feasible strategy because conformational changes are observed in the protease. Herein, structures and dynamics of ZIKV protease are summarized. The conformational changes of ZIKV protease and other proteases in the same family are discussed. The progress in developing allosteric inhibitors is also described. Understanding the structures and dynamics of the proteases are important for designing potent inhibitors.

Ligand Accessibility Insights to the Dengue Virus NS3-NS2B Protease Assessed by Long-Timescale Molecular Dynamics Simulations

Dengue is a tropical disease caused by the dengue virus (DENV), with an estimate of 300 million new cases every year. Due to the limited vaccine efficiency and absence of effective antiviral treatment, new drug candidates are urgently needed. DENV NS3-NS2B protease complex is essential for viral post-translational processing and maturation, and this enzyme has been extensively studied as a relevant drug target. Crystal structures often underestimate NS3-NS2B flexibility, whereas they can adopt different conformational states depending on the bound substrate. We conducted molecular dynamics simulations (∼30 μs) with a non- and covalently bound inhibitor to understand the conformational changes in the DENV-3 NS3-NS2B complex. Our results show that the open-closing movement of the protease exposes multiple druggable subpockets that can be investigated in later drug discovery efforts.

Targeting the protease of West Nile virus

West Nile virus infections can cause severe neurological symptoms. During the last 25 years, cases have been reported in Asia, North America, Africa, Europe and Australia (Kunjin). No West Nile virus vaccines or specific antiviral therapies are available to date. Various viral proteins and host-cell factors have been evaluated as potential drug targets. The viral protease NS2B-NS3 is among the most promising viral targets. It releases viral proteins from a non-functional polyprotein precursor, making it a critical factor of viral replication. Despite strong efforts, no protease inhibitors have reached clinical trials yet. Substrate-derived peptidomimetics have facilitated structural elucidations of the active protease state, while alternative compounds with increased drug-likeness have recently expanded drug discovery efforts beyond the active site.

Conformational Dynamics of the Dengue Virus Protease Revealed by Fluorescence Correlation and Single-Molecule FRET Studies

The dengue virus protease (DENV-PR) represents an attractive target for counteracting DENV infections. It is generally assumed that DENV-PR can exist in an open and a closed conformation and that active site directed ligands stabilize the closed state. While crystal structures of both the open and the closed conformation were successfully resolved, information about the prevalence of these conformations in solution remains elusive. Herein, we address the question of whether there is an equilibrium between different conformations in solution which can be influenced by addition of a competitive inhibitor. To this end, DENV-PR was statistically labeled by two dye molecules constituting a FRET (fluorescence resonance energy transfer) couple. Fluorescence correlation spectroscopy and photon-burst detection were employed to examine FRET pair labeled DENV-PRs freely diffusing in solution. The measurements were performed with two double mutants and with two dye couples. The data provide strong evidence that an equilibrium of at least two conformations of DENV-PR exists in solution. The competitive inhibitor stabilizes the closed state. Because the open and closed conformations appear to coexist in solution, our results support the picture of a conformational selection rather than that of an induced fit mechanism with respect to the inhibitor-induced formation of the closed state.

Protein structural heterogeneity: A hypothesis for the basis of proteolytic recognition by the main protease of SARS-CoV and SARS-CoV-2

The main protease (M^pro) of SARS-CoV and SARS-CoV-2 is a key enzyme in viral replication and a promising target for the development of antiviral therapeutics. The understanding of this protein is based on a number of observations derived from earlier x-ray structures, which mostly consider substrates or ligands as the main reason behind modulation of the active site. This lead to the concept of substrate-induced subsite cooperativity as an initial attempt to explain the dual binding specificity of this enzyme in recognizing the cleavage sequences at its N- and C-termini, which are important processing steps in obtaining the mature protease. The presented hypothesis proposes that structural heterogeneity is a property of the enzyme, independent of the presence of a substrate or ligand. Indeed, the analysis of M^pro structures of SARS-CoV and SARS-CoV-2 reveals a conformational diversity for the catalytically competent state in ligand-free structures. Variation in the binding site appears to result from flexibility at residues lining the S₁ subpocket and segments incorporating methionine 49 and glutamine 189. The structural evidence introduces “structure-based recognition” as a new paradigm in substrate proteolysis by M^pro. In this concept, the binding space in subpockets of the enzyme varies in a non-cooperative manner, causing distinct conformations, which recognize and process different cleavage sites, as the N- and C-termini. Insights into the recognition basis of the protease provide explanation to the ordered processing of cleavage sites. The hypothesis expands the conformational space of the enzyme and consequently opportunities for drug development and repurposing efforts.

Mechanisms of ligand binding

Many processes in chemistry and biology involve interactions of a ligand with its molecular target. Interest in the mechanism governing such interactions has dominated theoretical and experimental analysis for over a century. The interpretation of molecular recognition has evolved from a simple rigid body association of the ligand with its target to appreciation of the key role played by conformational transitions. Two conceptually distinct descriptions have had a profound impact on our understanding of mechanisms of ligand binding. The first description, referred to as induced fit, assumes that conformational changes follow the initial binding step to optimize the complex between the ligand and its target. The second description, referred to as conformational selection, assumes that the free target exists in multiple conformations in equilibrium and that the ligand selects the optimal one for binding. Both descriptions can be merged into more complex reaction schemes that better describe the functional repertoire of macromolecular systems. This review deals with basic mechanisms of ligand binding, with special emphasis on induced fit, conformational selection, and their mathematical foundations to provide rigorous context for the analysis and interpretation of experimental data. We show that conformational selection is a surprisingly versatile mechanism that includes induced fit as a mathematical special case and even captures kinetic properties of more complex reaction schemes. These features make conformational selection a dominant mechanism of molecular recognition in biology, consistent with the rich conformational landscape accessible to biological macromolecules being unraveled by structural biology.