Structure-guided mutagenesis of active site residues in the dengue virus two-component protease NS2B-NS3

Exploring the antiviral inhibitory activity of Niloticin against the NS2B/NS3 protease of Dengue virus (DENV2)

Dengue virus (DENV2) is the cause of dengue disease and a worldwide health problem. DENV2 replicates in the host cell using polyproteins such as NS3 protease in conjugation with NS2B cofactor, making NS3 protease a promising antiviral drug-target. This study investigated the efficacy of 'Niloticin' against NS2B/NS3-protease. In silico and in vitro analyses were performed which included interaction of niloticin with NS2B/NS3-protease, protein stability and flexibility, mutation effect, betweenness centrality of residues and analysis of cytotoxicity, protein expression and WNV NS3-protease activity. Similar like acyclovir, niloticin forms strong H-bonds and hydrophobic interactions with residues LEU149, ASN152, LYS74, GLY148 and ALA164. The stability of the niloticin-NS2B/NS3-protease complex was found to be stable compared to the apo NS2B/NS3-protease in structural deviation, PCA, compactness and FEL analysis. The IC50 value of niloticin was 0.14 μM in BHK cells based on in vitro cytotoxicity analysis and showed significant activity at 2.5 μM in a concentration-dependent manner. Western blotting and qRT-PCR analyses showed that niloticin reduced DENV2 protein transcription in a dose-dependent manner. Besides, niloticin confirmed the inhibition of NS3-protease by the SensoLyte 440 WNV protease detection kit. These promising results suggest that niloticin could be an effective antiviral drug against DENV2 and other flaviviruses.

Discovery of azaleatin as a potential allosteric inhibitor for dengue NS2B-NS3 protease using in vitro and in silico studies

The rise in dengue cases in tropical and sub-tropical areas has become a significant health concern. At present, there is no definitive cure for dengue fever, which underscores the importance of identifying potent inhibitors. Dengue NS2B-NS3 protease is the prime drug target due to its vital function for replication. Quercetin, a flavone, has anti-dengue virus properties but is limited by low bioavailability. Previous studies have shown that methoxy substitution in flavones improves bioavailability and metabolic stability. Azaleatin is a derivative of quercetin with a methoxy substitution at the C5 position, however its ability to inhibit dengue is unknown. In this study, azaleatin was investigated for its inhibition against dengue NS2B-NS3 protease using in vitro and in silico techniques. The fluorescence assay was used to determine the IC50 value and inhibition kinetics. The molecular interaction between azaleatin and NS2B-NS3 was studied using CB-Dock2 and AutoDock Vina. The complex's stability was then analysed using GROMACS. Besides, the ADMETlab 2.0 was utilized to predict pharmacokinetic of the azaleatin. Results showed that azaleatin inhibits dengue NS2B-NS3 protease non-competitively with a Ki of 26.82 µg/ml and an IC50 of 38 µg/ml. Molecular docking indicated binding of the azaleatin to the allosteric pocket of NS2B-NS3 with a docking score of -8.2 kcal/mol. Azaleatin was found stable in the pocket along 100 ns, supporting its inhibitory mode. The compound has favourable pharmacokinetic profiles and conformed to Lipinski's Rule of Five. Taken together, azaleatin inhibits NS2B-NS3 protease in a non-competitive mode, suggesting its potential as safer anti-dengue compound.Communicated by Ramaswamy H. Sarma.

Loop Dynamics and Conformational Flexibility in Dengue Serine Protease Activity: Noninvasive Perturbation by Solvent Exchange

Molecular mechanics play an important role in enzyme action and understanding the dynamics of loop motion is key for designing inhibitors of an enzyme, particularly targeting the allosteric sites. For the successful creation of new protease inhibitors targeting the dengue serine protease, our current investigation detailed the intricate structural dynamics of NS2B/NS3 dengue protease. This enzyme is one of the most essential enzymes in the life cycle of the dengue virus, which is responsible for the activation/processing of viral polyprotein, thus making it a potential target for drug discovery. We showed that the internal dynamics of two regions, fingers 1 and 2 (R24-G39 and L149-A164, respectively) adjacent to the active site triad of this protease, control the enzyme action. Each of these regions is composed of two antiparallel β-strands connected by β-turn/hairpin loops. The correlated bending and rocking motions in the two β-turns on either side of the active site were found to modulate the activity of the enzyme to a large extent. With increasing concentration of cosolvent dimethyl sulfoxide, correlated motions in the finger 2 region get diminished and bending of finger 1 increases, which are also reflected in the loss of enzyme activity. Decreasing temperature and mutations in neighboring nonsubstrate binding residues show similar effects on loop motion and enzyme kinetics. Therefore, in vitro noninvasive perturbation of these motions by the solvent exchange as well as cold stress in combination with in silico molecular dynamics simulations established the importance of the two β-turns in the functioning of dengue virus serotype 2 NS2B/NS3 serine protease.

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.

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.

MD simulations reveal alternate conformations of the oxyanion hole in the Zika virus NS2B/NS3 protease

Recent crystallography studies have shown that the binding site oxyanion hole plays an important role in inhibitor binding, but can exist in two conformations (active/inactive). We have undertaken molecular dynamics (MD) calculations to better understand oxyanion hole dynamics and thermodynamics. We find that the Zika virus (ZIKV) NS2B/NS3 protease maintains a stable closed conformation over multiple 100-ns conventional MD simulations in both the presence and absence of inhibitors. The S1, S2, and S3 pockets are stable as well. However, in two of eight simulations, the A132-G133 peptide bond in the binding pocket of S1' spontaneously flips to form a 3₁₀ -helix that corresponds to the inactive conformation of the oxyanion hole, and then maintains this conformation until the end of the 100-ns conventional MD simulations without inversion of the flip. This conformational change affects the S1' pocket in ZIKV NS2B/NS3 protease active site, which is important for small molecule binding. The simulation results provide evidence at the atomic level that the inactive conformation of the oxyanion hole is more favored energetically when no specific interactions are formed between substrate/inhibitor and oxyanion hole residues. Interestingly, however, transition between the active and inactive conformation of the oxyanion hole can be observed by boosting the valley potential in accelerated MD simulations. This supports a proposed induced-fit mechanism of ZIKV NS2B/NS3 protease from computational methods and provides useful direction to enhance inhibitor binding predictions in structure-based drug design.

Discovering key residues of dengue virus NS2b-NS3-protease: New binding sites for antiviral inhibitors design

The NS2B-NS3 protease is essential for the Dengue Virus (DENV) replication process. This complex constitutes a target for efficient antiviral discovery because a drug could inhibit the viral polyprotein processing. Furthermore, since the protease is highly conserved between the four Dengue virus serotypes, it is probable that a drug would be equally effective against all of them. In this article, a strategy is reported that allowed us to identify influential residues on the function of the Dengue NS2b-NS3 Protease. Moreover, this is a strategy that could be applied to virtually any protein for the search of alternative influential residues, and for non-competitive inhibitor development. First, we incorporated several features derived from computational alanine scanning mutagenesis, sequence, structure conservation, and other structure-based characteristics. Second, these features were used as variables to obtain a multilayer perceptron model to identify defined groups (clusters) of key residues as possible candidate pockets for binding sites of new leads on the DENV protease. The identified residues included: i) amino acids close to the beta sheet-loop-beta sheet known to be important in its closed conformation for NS2b ii) residues close to the active site, iii) several residues evenly spread on the NS2b-NS3 contact surface, and iv) some inner residues most likely related to the overall stability of the protease. In addition, we found concordance on our list of residues with previously identified amino acids part of a highly conserved peptide studied for vaccine development.

Sequence, structure and function relationships in flaviviruses as assessed by evolutive aspects of its conserved non-structural protein domains

Flaviviruses are responsible for serious diseases such as dengue, yellow fever, and zika fever. Their genomes encode a polyprotein which, after cleavage, results in three structural and seven non-structural proteins. Homologous proteins can be studied by conservation and coevolution analysis as detected in multiple sequence alignments, usually reporting positions which are strictly necessary for the structure and/or function of all members in a protein family or which are involved in a specific sub-class feature requiring the coevolution of residue sets. This study provides a complete conservation and coevolution analysis on all flaviviruses non-structural proteins, with results mapped on all well-annotated available sequences. A literature review on the residues found in the analysis enabled us to compile available information on their roles and distribution among different flaviviruses. Also, we provide the mapping of conserved and coevolved residues for all sequences currently in SwissProt as a supplementary material, so that particularities in different viruses can be easily analyzed.

In Silico Screening, Alanine Mutation, and DFT Approaches for Identification of NS2B/NS3 Protease Inhibitors

To identify the ligand that binds to a target protein with high affinity is a nontrivial task in computer-assisted approaches. Antiviral drugs have been identified for NS2B/NS3 protease enzyme on the mechanism to cleave the viral protein using the computational tools. The consequence of the molecular docking, free energy calculations, and simulation protocols explores the better ligand. It provides in-depth structural insights with the catalytic triad of His51, Asp75, Ser135, and Gly133. The MD simulation was employed here to predict the stability of the complex. The alanine mutation has been performed and its stability was monitored by using the molecular dynamics simulation. The minimal RMSD value suggests that the derived complexes are close to equilibrium. The DFT outcome reveals that the HOMO-LUMO gap of Ligand19 is 2.86 kcal/mol. Among the considered ligands, Ligand19 shows the lowest gap and it is suggested that the HOMO of Ligand19 may transfer the electrons to the LUMO in the active regions. The calculated binding energy of Ligand19 using the DFT method is in good agreement with the docking studies. The pharmacological activity of ligand was performed and satisfies Lipinski rule of 5. Moreover, the computational results are compared with the available IC₅₀ values of experimental results.

Discovery of antiviral molecules for dengue: In silico search and biological evaluation

BACKGROUND

Dengue disease is a global disease that has no effective treatment. The dengue virus (DENV) NS2B/NS3 protease complex is a target for designing specific antivirals due to its importance in viral replication and its high degree of conservation.

METHODS

NS2B/NS3 protease complex structural information was employed to find small molecules that are capable of inhibiting the activity of the enzyme complex. This inhibitory activity was probed with in vitro assays using a fluorescent substrate and the complex NS2B/NS3 obtained by recombinant DNA techniques. HepG2 cells infected with dengue virus serotype 2 were used to test the activity against dengue virus replication.

RESULTS

A total of 210,903 small molecules from PubChem were docked in silico to the NS2B/NS3 structure (PDB: 2FOM) to find molecules that were capable of inhibiting this protein complex. Five of the best 500 leading compounds, according to their affinity values (-11.6 and -13.5 kcal/mol), were purchased. The inhibitory protease activity on the recombinant protein and antiviral assays was tested.

CONCLUSIONS

Chemicals CID 54681617, CID 54692801 and CID 54715399 were strong inhibitors of NS2B/NS3, with IC50 values (μM) and percentages of viral titer reductions of 19.9, 79.9%; 17.5, 69.8%; and 9.1, 73.9%, respectively. Multivariate methods applied to the molecular descriptors showed two compounds that were structurally different from other DENV inhibitors.

GENERAL SIGNIFICANCE

This discovery opens new possibilities for obtaining drug candidates against Dengue virus.

A straightforward experimental approach to expression, purification, refolding, and enzymatic analysis of recombinant dengue virus NS2B(H)-NS3pro protease

Dengue virus threatens around 2.5 billion people worldwide; about 50 million become infected every year, and yet no vaccine or drug is available for prevention and/or treatment. The flaviviral NS2B-NS3pro complex is indispensable for flaviviral replication and is considered to be an important drug target. The aim of this study was to develop a simple and generally applicable experimental strategy to construct, purify, and assay a highly active recombinant NS2B(H)-NS3pro complex that would be useful for high-throughput screening of potential inhibitors. The sequence of NS2B(H)-NS3pro was generated by overlap extension PCR (SOE-PCR) and cloned into the pTrcHisA vector. Hexahistidine-tagged NS2B(H)-NS3pro complex was expressed in E. coli predominantly as insoluble protein and purified to >95% purity by single-step immobilized metal affinity chromatography. SDS-PAGE followed by immunoblotting of the purified enzyme demonstrated the presence of the NS2B(H)-NS3pro precursor and its autocleavage products, NS3pro and NS2B(H), as 37, 21, and 10 kDa bands, respectively. Kinetic parameters, Km, kcat, and kcat/Km for the fluorophore-linked protease model substrate Ac-nKRR-amc were obtained using inner-filter effect correction. The kinetic parameters Km, kcat, and kcat/Km for Ac-nKRR-amc substrate were 100 µM, 0.112 s(-1), and 1120 M(-1)·s(-1), respectively. A simplified procedure for the cloning, overexpression, and purification of the NS2B(H)-NS3pro complex was applied, and a highly active recombinant NS2B(H)-NS3pro complex was obtained that could be useful for the design of high-throughput assays aimed at flaviviral inhibitor discovery.

Global consensus sequence development and analysis of dengue NS3 conserved domains

The dengue virus (DENV) genome encodes 10 different genes including the NS3 gene, which has a protease and helicase domain used in virus replication. This domain is a potential target for antiviral agents against dengue. Due to a high mutation rate, DENV is classified into four major serotypes (DENV1–DENV4). This study was designed to perform conservancy analysis of all four serotypes by drawing a consensus sequence for each serotype and then drawing a global consensus sequence to study conserved residues in all four serotypes. A total of 127 NS3 sequences belonging to all four serotypes were retrieved and aligned using multiple alignment feature of CLC Workbench and were subjected to phylogenetic tree construction. Conservancy analysis of NS3 revealed conserved peptides with active site residues that can be important in developing antiviral agents against dengue virus. Among conserved residues, residues G142, Ser144, and G145 (catalytic pocket residues), A219, D220, and D221 (divalent cations binding residues), and His56, Asp79, Ser144, 146 were highly conserved among all the serotypes. Residues from L138 to L149 and from L226 to L245 were also considerably conserved in all serotypes, while lysine141 mutated to serine in serotype 3. A total of 14 peptides from the conserved regions of DENV NS3 protein were identified, which may be helpful to develop peptide inhibitors. The DENV NS3 phylogenetic tree showed the evolutionary relationship among all four serotypes, and all serotypes of dengue were found to have evolved from the dengue 4 serotype. Because of its high variability, DENV has become a global health concern. It is important to study residues that are present in protease, helicase, the catalytic pocket Mg²⁺ binding site, and the AAA domain. This study revealed peptides with active site residues that are highly conserved among all four serotypes. These regions of the NS3 sequence may be helpful in developing antiviral agents.

New binding site conformations of the dengue virus NS3 protease accessed by molecular dynamics simulation

Dengue fever is caused by four distinct serotypes of the dengue virus (DENV1-4), and is estimated to affect over 500 million people every year. Presently, there are no vaccines or antiviral treatments for this disease. Among the possible targets to fight dengue fever is the viral NS3 protease (NS3_PRO), which is in part responsible for viral processing and replication. It is now widely recognized that virtual screening campaigns should consider the flexibility of target protein by using multiple active conformational states. The flexibility of the DENV NS3_PRO could explain the relatively low success of previous virtual screening studies. In this first work, we explore the DENV NS3_PRO conformational states obtained from molecular dynamics (MD) simulations to take into account protease flexibility during the virtual screening/docking process. To do so, we built a full NS3_PRO model by multiple template homology modeling. The model comprised the NS2B cofactor (essential to the NS3_PRO activation), a glycine flexible link and the proteolytic domain. MD simulations had the purpose to sample, as closely as possible, the ligand binding site conformational landscape prior to inhibitor binding. The obtained conformational MD sample was clustered into four families that, together with principal component analysis of the trajectory, demonstrated protein flexibility. These results allowed the description of multiple binding modes for the Bz-Nle-Lys–Arg–Arg-H inhibitor, as verified by binding plots and pair interaction analysis. This study allowed us to tackle protein flexibility in our virtual screening campaign against the dengue virus NS3 protease.

Enzymatic analysis of recombinant Japanese encephalitis virus NS2B(H)-NS3pro protease with fluorogenic model peptide substrates

Background

Japanese encephalitis virus (JEV), a member of the Flaviviridae family, causes around 68,000 encephalitis cases annually, of which 20–30% are fatal, while 30–50% of the recovered cases develop severe neurological sequelae. Specific antivirals for JEV would be of great importance, particularly in those cases where the infection has become persistent. Being indispensable for flaviviral replication, the NS2B-NS3 protease is a promising target for design of anti-flaviviral inhibitors. Contrary to related flaviviral proteases, the JEV NS2B-NS3 protease is structurally and mechanistically much less characterized. Here we aimed at establishing a straightforward procedure for cloning, expression, purification and biochemical characterization of JEV NS2B(H)-NS3pro protease.

Methodology/Principal Findings

The full-length sequence of JEV NS2B-NS3 genotype III strain JaOArS 982 was obtained as a synthetic gene. The sequence of NS2B(H)-NS3pro was generated by splicing by overlap extension PCR (SOE-PCR) and cloned into the pTrcHisA vector. Hexahistidine-tagged NS2B(H)-NS3pro, expressed in E. coli as soluble protein, was purified to >95% purity by a single-step immobilized metal affinity chromatography. SDS-PAGE and immunoblotting of the purified enzyme demonstrated NS2B(H)-NS3pro precursor and its autocleavage products, NS3pro and NS2B(H), as 36, 21, and 10 kDa bands, respectively. Kinetic parameters, K_m and k_cat, for fluorogenic protease model substrates, Boc-GRR-amc, Boc-LRR-amc, Ac-nKRR-amc, Bz-nKRR-amc, Pyr-RTKR-amc and Abz-(R)₄SAG-nY-amide, were obtained using inner filter effect correction. The highest catalytic efficiency k_cat/K_m was found for Pyr-RTKR-amc (k_cat/K_m: 1962.96±85.0 M⁻¹ s⁻¹) and the lowest for Boc-LRR-amc (k_cat/K_m: 3.74±0.3 M⁻¹ s⁻¹). JEV NS3pro is inhibited by aprotinin but to a lesser extent than DEN and WNV NS3pro.

Conclusions/Significance

A simplified procedure for the cloning, overexpression and purification of the NS2B(H)-NS3pro was established which is generally applicable to other flaviviral proteases. Kinetic parameters obtained for a number of model substrates and inhibitors, are useful for the characterization of substrate specificity and eventually for the design of high-throughput assays aimed at antiviral inhibitor discovery.

Ligand-bound structures of the dengue virus protease reveal the active conformation

Dengue is a mosquito-borne viral hemorrhagic disease that is a major threat to human health in tropical and subtropical regions. Here we report crystal structures of a peptide covalently bound to dengue virus serotype 3 (DENV-3) protease as well as the serine-protease inhibitor aprotinin bound to the same enzyme. These structures reveal, for the first time, a catalytically active, closed conformation of the DENV protease. In the presence of the peptide, the DENV-3 protease forms the closed conformation in which the hydrophilic β-hairpin region of NS2B wraps around the NS3 protease core, in a manner analogous to the structure of West Nile virus (WNV) protease. Our results confirm that flavivirus proteases form the closed conformation during proteolysis, as previously proposed for WNV. The current DENV-3 protease structures reveal the detailed interactions at the P4' to P3 sites of the substrate. The new structural information explains the sequence preference, particularly for long basic residues in the nonprime side, as well as the difference in substrate specificity between the WNV and DENV proteases at the prime side. Structural analysis of the DENV-3 protease-peptide complex revealed a pocket that is formed by residues from NS2B and NS3; this pocket also exists in the WNV NS2B/NS3 protease structure and could be targeted for potential antivirus development. The structural information presented in the current study is invaluable for the design of specific inhibitors of DENV protease.

Important advances in the field of anti-dengue virus research

There are currently no licensed antivirals available for the treatment of dengue virus (DENV), which causes significant morbidity and mortality throughout tropical areas of the world and is now encroaching on the southern United States. Recent improvements in existing animal models and cell culture systems have been very important in elucidating the mechanisms of DENV pathogenesis in humans, including the identification of potential viral and host proteins that might be targeted for the treatment of DENV infection. The AG129 mouse model is a major advance in the development of antiviral and vaccine candidates for clinical use. It allows for testing of potential therapeutics in a relevant system that exhibits some aspects of disease that are similar to those observed in humans. This review focuses on recent developments in the AG129 mouse model and discusses compounds that have been found to be active in available cell and animal model systems within the past year.