Homology modeling and molecular dynamics simulations of Dengue virus NS2B/NS3 protease: insight into molecular interaction

The pathogenic West Nile virus (WNV) and Dengue virus (DV) are growing global threats for which there are no specific treatments. Both viruses possess a two component NS2B/NS3 protease which cleaves viral precursor proteins. Whereas for the WNV protease two crystal structures in complex with an inhibitor have been solved recently, no such information is available for the DV protease. Here, we report the generation of a homology model of DV NS2B/NS3 protease. Since it is known from the related WNV protease that it adopts a distinct conformation in free and in inhibitor‐complexed form, a special emphasis was given to the analysis of the protease flexibility. Therefore, several models of DV NS2B/NS3 protease complexed with the peptidic inhibitor (Bz‐Nle(P4)‐Lys(P3)‐Arg(P2)‐Arg(P1)‐H) were generated. The first DV protease model (DV‐1) was constructed using the available crystal structure of the apo DV NS2B/NS3 protease. The second model (DV‐2) was built taking the WNV NS3/NS2B protease in the inhibitor‐complexed form as the template structure. Molecular dynamics simulations which were carried out for the WNV crystal structures as well as for the DV models provided an understanding of the role of NS2B for maintaining the protease in the active conformation. It was also demonstrated that NS2B is not only important for maintaining NS3 in the active form, but is also essential for establishing the interaction between residues from the S2 pocket and the peptidic inhibitor. The DV NS2B/NS3 model in the productive conformation can now be used for structure‐based design purposes. Copyright © 2009 John Wiley & Sons, Ltd.     

Identification of novel thiadiazoloacrylamide analogues as inhibitors of dengue-2 virus NS2B/NS3 protease

Dengue virus is endemic throughout tropical and subtropical regions, and cause severe epidemic diseases. The NS2B/NS3 protease is a promising drug target for dengue virus. Herein, we report the discovery and modification of a novel class of thiadiazoloacrylamide derivatives with potent inhibitory activity against the NS2B/NS3 protease. Thiadiazolopyrimidinone 1 was firstly determined as a new chemical structure against NS2B/NS3 from a commercial compound library. Then, we sought to identify similar compounds with the thiadiazoloacrylamide core that would exhibit better activity. A series of analogues were synthesized and fourteen of them were identified with strong inhibitory activities, in which the nitrile group in the linker part was discovered as an essential group for the inhibitory activity. The best of these (8b) demonstrated an IC₅₀ at 2.24 μM based on in vitro DENV2 NS2B-NS3pro assays.     

An inhibition model of BPTI to unlinked dengue virus NS2B-NS3 protease

One approach to treating the dengue virus infection is to inhibit its NS2B-NS3 protease that plays a vital role in virus maturation. However, the lack of structural information on the active conformation of the protease hindered related drug design. With a co-expression system, we obtained the active two-component protease in its unlinked form. BPTI shows strong competitive inhibitory activity (K_i = 6.5 nM) against this unlinked protease, which adopts a closed conformation. Based on the biochemical and NMR perturbation information, an inhibition model of BPTI to NS2B-NS3 protease is proposed.     

Discovery of novel cyclic peptide inhibitors of dengue virus NS2B-NS3 protease with antiviral activity

NS2B-NS3 protease is an essential enzyme for the replication of dengue virus (DENV), which continues to be a serious threat to worldwide public health. We designed and synthesized a series of cyclic peptides mimicking the substrates of this enzyme, and assayed their activity against the DENV-2 NS2B-NS3 protease. The introduction of aromatic residues at the appropriate positions and conformational restriction generated the most promising cyclic peptide with an IC₅₀ of 0.95 μM against NS2B-NS3 protease. Cyclic peptides with proper positioning of additional arginines and aromatic residues exhibited antiviral activity against DENV. Furthermore, replacing the C-terminal amide bond of the polybasic amino acid sequence with an amino methylene moiety stabilized the cyclic peptides against hydrolysis by NS2B-NS3 protease, while maintaining their enzyme inhibitory activity and antiviral activity.     

Phosphonate inhibitors of West Nile virus NS2B/NS3 protease

West Nile virus (WNV) is a member of the flavivirus genus belonging to the Flaviviridae family. The viral serine protease NS2B/NS3 has been considered an attractive target for the development of anti-WNV agents. Although several NS2B/NS3 protease inhibitors have been described so far, most of them are reversible inhibitors. Herein, we present a series of α-aminoalkylphosphonate diphenyl esters and their peptidyl derivatives as potent inhibitors of the NS2B/NS3 protease. The most potent inhibitor identified was Cbz-Lys-Arg-(4-GuPhe)^P(OPh)₂ displaying K_i and k₂/K_i values of 0.4 µM and 28 265 M^?1s^?1, respectively, with no significant inhibition of trypsin, cathepsin G, and HAT protease.     

A small compound targeting the interaction between nonstructural proteins 2B and 3 inhibits dengue virus replication

The non-structural protein NS2B/NS3 serine-protease complex of the dengue virus (DENV) is required for the maturation of the viral polyprotein. Dissociation of the NS2B cofactor from NS3 diminishes the enzymatic activity of the complex. In this study, we identified a small molecule inhibitor that interferes with the interaction between NS2B and NS3 using structure-based screening and a cell-based viral replication assay. A library containing 661,417 small compounds derived from the Molecular Operating Environment lead-like database was docked to the NS2B/NS3 structural model. Thirty-nine compounds with high scores were tested in a secondary screening using a cell-based viral replication assay. SK-12 was found to inhibit replication of all DENV serotypes (EC₅₀ = 0.74–4.92 μM). In silico studies predicted that SK-12 pre-occupies the NS2B-binding site of NS3. Steady-state kinetics using a fluorogenic short peptide substrate demonstrated that SK-12 is a noncompetitive inhibitor against the NS2B/NS3 protease. Inhibition to Japanese encephalitis virus by SK-12 was relatively weak (EC₅₀ = 29.81 μM), and this lower sensitivity was due to difference in amino acid at position 27 of NS3. SK-12 is the promising small-molecule inhibitor that targets the interaction between NS2B and NS3.     

Effects of NS2B-NS3 protease and furin inhibition on West Nile and Dengue virus replication

West Nile virus (WNV) and Dengue virus (DENV) replication depends on the viral NS2B-NS3 protease and the host enzyme furin, which emerged as potential drug targets. Modification of our previously described WNV protease inhibitors by basic phenylalanine analogs provided compounds with reduced potency against the WNV and DENV protease. In a second series, their decarboxylated P1-trans-(4-guanidino)cyclohexylamide was replaced by an arginyl-amide moiety. Compound 4-(guanidinomethyl)-phenylacetyl-Lys-Lys-Arg-NH₂ inhibits the NS2B-NS3 protease of WNV with an inhibition constant of 0.11?µM. Due to the similarity in substrate specificity, we have also tested the potency of our previously described multibasic furin inhibitors. Their further modification provided chimeric inhibitors with additional potency against the WNV and DENV proteases. A strong inhibition of WNV and DENV replication in cell culture was observed for the specific furin inhibitors, which reduced virus titers up to 10,000-fold. These studies reveal that potent inhibitors of furin can block the replication of DENV and WNV.     

Homology modeling and molecular dynamics study of West Nile virus NS3 protease: a molecular basis for the catalytic activity increased by the NS2B cofactor

The West Nile virus (WNV) NS3 serine protease, which plays an important role in assembly of infective virion, is an attractive target for anti-WNV drug development. Cofactors NS2B and NS4A increase the catalytic activity of NS3 in dengue virus and Hepatitis C virus, respectively. Recent studies on the WNV-NS3 characterize the catalytically active form of NS3 by tethering the 40-residue cofactor NS2B. It is suggested that NS2B is essential for the NS3 activity in WNV, while there is no information of the WNV-NS3-related crystal structure. To understand the role of NS2B/substrate in the NS3 catalytic activity, we built a series of models: WNV-NS3 and WNV-NS3-NS2B and WNV-NS3-NS2B-substrate using homology modeling and molecular modeling techniques. Molecular dynamics (MD) simulations were performed for 2.75 ns on each model, to investigate the structural stabilization and catalytic triad motion of the WNV NS3 protease with and without NS2B/substrate. The simulations show that the NS3 rearrangement occurs upon the NS2B binding, resulting in the stable D75-OD1...H51-NH hydrogen bonding. After the substrate binds to the NS3-NS2B active site, the NS3 protease becomes more stable, and the catalytic triad is formed. These results provide a structural basis for the activation and stabilization of the enzyme by its cofactor and substrate.     

Competitive inhibition of the dengue virus NS3 serine protease by synthetic peptides representing polyprotein cleavage sites

The NS3 serine protease of dengue virus is required for the maturation of the viral polyprotein and consequently represents a promising target for the development of antiviral inhibitors. However, the substrate specificity of this enzyme has been characterized only to a limited extent. In this study, we have investigated product inhibition of the NS3 protease by synthetic peptides derived from the P6-P1 and the P1'-P5' regions of the natural polyprotein substrate. N-terminal cleavage site peptides corresponding to the P6-P1 region of the polyprotein were found to act as competitive inhibitors of the enzyme with K(i) values ranging from 67 to 12 microM. The lowest K(i) value was found for the peptide representing the NS2A/NS2B cleavage site, RTSKKR. Inhibition by this cleavage site sequence was analyzed by using shorter peptides, SKKR, KKR, KR, AGRR, and GKR. With the exception of the peptide AGRR which did not inhibit the protease at a concentration of 1mM, all other peptides displayed K(i) values in the range from 188 to 22 microM. Peptides corresponding to the P1'-P5' region of the polyprotein cleavage sites had no effect on enzymatic activity at a concentration of 1mM. Molecular docking data of peptide inhibitors to a homology-based model of the dengue virus type 2 NS2B(H)-NS3p co-complex indicate that binding of the non-prime site product inhibitors is similar to ground-state binding of the corresponding substrates.     

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.     

Complex of NS3 protease and NS4A peptide of BK strain hepatitis C virus: a 2.2 A resolution structure in a hexagonal crystal form.

The crystal structure of the NS3 protease of the hepatitis C virus (BK strain) has been determined in the space group P6(3)22 to a resolution of 2.2 A. This protease is bound with a 14-mer peptide representing the central region of the NS4A protein. There are two molecules of the NS3(1-180)-NS4A(21'-34') complex per asymmetric unit. Each displays a familiar chymotrypsin-like fold that includes two beta-barrel domains and four short alpha-helices. The catalytic triad (Ser-139, His-57, and Asp-81) is located in the crevice between the beta-barrel domains. The NS4A peptide forms an almost completely enclosed peptide surface association with the protease. In contrast to the reported H strain complex of NS3 protease-NS4A peptide in a trigonal crystal form (Kim JL et al., 1996, Cell 87:343-355), the N-terminus of the NS3 protease is well-ordered in both molecules in the asymmetric unit of our hexagonal crystal form. The folding of the N-terminal region of the NS3 protease is due to the formation of a three-helix bundle as a result of crystal packing. When compared with the unbound structure (Love RA et al., 1996, Cell 87:331-342), the binding of the NS4A peptide leads to the ordering of the N-terminal 28 residues of the NS3 protease into a beta-strand and an alpha-helix and also causes local rearrangements important for a catalytically favorable conformation at the active site. Our analysis provides experimental support for the proposal that binding of an NS4A-mimicking peptide, which increases catalytic rates, is necessary but not sufficient for formation of a well-ordered, compact and, hence, highly active protease molecule.     

Characterisation of the role of zinc in the hepatitis C virus NS2/3 auto-cleavage and NS3 protease activities

Cleavage of the hepatitis C virus polyprotein between the non-structural NS2 and NS3 proteins is mediated by a poorly characterised auto-proteolytic activity that maps to the C terminus of NS2 and the N terminus of NS3, but is distinct from the NS3 protease activity responsible for downstream cleavages in the polyprotein. We have exploited the fact that the minimal precursor (residues 904-1206 of the HCV polyprotein) can be expressed as an insoluble protein in Escherichia coli and subsequently refolded into a form active for both auto-cleavage and NS3 protease activity, to further characterise the NS2/3 auto-cleavage activity. We show that both activities are zinc-dependent and show an absolute requirement for cysteine residues 1123, 1125 and 1171 within NS3. In contrast cysteine 922 (within NS2) is only required for NS2/3 auto-cleavage activity and histidine 1175 is only required for NS3 activity. Although the complete NS3 protease domain (including the C-terminal alpha-helix) is required for NS2/3 auto-cleavage, the activity of the NS3 protease is not essential. Lastly we show that the NS2/3 auto-cleavage activity is more sensitive to zinc chelation by 1,10-phenanthroline than the NS3 protease activity. This observation is consistent with different conformations of the precursor competent for either NS2/3 auto-cleavage or NS3 protease activity; these two conformations can be distinguished by their relative strength and geometry of zinc coordination.     

NS3 protease from flavivirus as a target for designing antiviral inhibitors against dengue virus

The development of novel therapeutic agents is essential for combating the increasing number of cases of dengue fever in endemic countries and among a large number of travelers from non-endemic countries. The dengue virus has three structural proteins and seven non-structural (NS) proteins. NS3 is a multifunctional protein with an N-terminal protease domain (NS3pro) that is responsible for proteolytic processing of the viral polyprotein, and a C-terminal region that contains an RNA triphosphatase, RNA helicase and RNA-stimulated NTPase domain that are essential for RNA replication. The serine protease domain of NS3 plays a central role in the replicative cycle of dengue virus. This review discusses the recent structural and biological studies on the NS2B-NS3 protease-helicase and considers the prospects for the development of small molecules as antiviral drugs to target this fascinating, multifunctional protein.     

Designing cyclopentapeptide inhibitor as potential antiviral drug for dengue virus ns5 methyltransferase

NS5 methyltransferase (Mtase) has a crucial role in the replication of dengue virus. There are two active sites on NS5 Mtase i.e., SAM and RNA-cap binding sites. Inhibition of the NS5 Mtase activity is expected to prevent the propagation of dengue virus. This study was conducted to design cyclic peptide ligands as enzyme inhibitors of dengue virus NS5 Mtase through computational approach. Cyclopentapeptides were designed as ligand of SAM binding site as much as 1635 and 736 cyclopentpeptides were designed as ligand of RNA-cap binding site. Interaction between ligand and NS5 Mtase has been conducted on the Docking simulation. The result shows that cyclopentapeptide CTWYC was the best peptide candidate on SAM binding site, with estimated free binding energy -30.72 kca/mol. Cyclopentapeptide CYEFC was the best peptide on RNA-cap binding site with estimated free binding energy -22.89 kcal/mol. Both peptides did not have tendency toward toxicity properties. So it is expected that both CTWYC and CYEFC ligands could be used as a potential antiviral drug candidates, which can inhibit the SAM and RNA-cap binding sites of dengue virus NS5 Mtase.     

圣路易斯脑炎病毒NS2B-NS3蛋白酶的原核表达及其抑制剂的筛选

【目的】圣路易斯脑炎病毒(St. Louis encephalitis virus,SLEV)属于黄病毒科,是一种单股正链RNA病毒。黄病毒编码的非结构蛋白NS3在病毒复制以及多聚蛋白加工过程中起着重要作用,NS2B是其发挥作用的重要辅助因子。因此,NS2B-NS3蛋白酶复合物是抗病毒药物的重要靶标。本研究旨在构建SLEV NS2B-NS3蛋白酶的原核表达系统并建立其抑制剂的高通量筛选方法,从而发现其小分子抑制剂。【方法】通过PCR扩增SLEVNS2B-NS3蛋白的编码区,构建原核表达质粒;在大肠杆菌BL21(DE3)中,经异丙基硫代半乳糖苷(Isopropyl β-D-thiogalactoside)诱导得到可溶性的NS2B-NS3蛋白,并用镍亲和层析方法进行纯化;基于荧光共振能量转移(Fluorescence resonance energy transfer)技术检测NS2B-NS3蛋白酶活性,建立其抑制剂的高通量筛选平台。【结果】SLEV NS2B-NS3蛋白酶纯化程度高达95%以上,基于酶活测定的抑制剂筛选平台准确可行。对700多个上市药物进行筛选后,发现原花青素对SLEVNS2B-NS3蛋白酶具有明显的抑制活性。【结论】本研究为SLEVNS2B-NS3蛋白酶抑制剂提供了一种操作方便、高通量的筛选方法,并首次发现了原花青素具有抑制SLEV NS2B-NS3蛋白酶活性的功能,可以作为治疗SLEV感染的潜在靶向药物。     

An RNA ligand inhibits hepatitis C virus NS3 protease and helicase activities

The hepatitis C virus non-structural protein 3 (HCV NS3) possesses both protease and helicase activities that are essential for viral replication. In a previous study, we obtained RNA aptamers that specifically and efficiently inhibited NS3 protease activity (G9 aptamers). In order to add helicase-inhibition capability, we attached (U)14 to the 3'-terminal end of a minimized G9 aptamer, DeltaNEO-III. NEO-III-14U was shown to inhibit the NS3 protease activity more efficiently than the original aptamer and, furthermore, to efficiently inhibit the unwinding reaction by NS3 helicase. In addition, NEO-III-14U has the potential to diminish specific interactions between NS3 and the 3'-UTR of HCV-positive and -negative strands. NEO-III-14U showed effective inhibition against NS3 protease in living cells.     

Computer Aided Screening of Phytochemicals from Garcinia against the Dengue NS2B/NS3 Protease

Dengue virus NS2/NS3 protease because of its ability to cleave viral proteins is considered as an attractive target to screen antiviral agents. Medicinal plants contain a variety of phytochemicals that can be used as drug against different diseases and infections. Therefore, this study was designed to uncover possible phytochemical of different classes (Aromatic, Carbohydrates, Lignin, Saponins, Steroids, Tannins, Terpenoids, Xanthones) that could be used as inhibitors against the NS2B/NS3 protease of DENV. With the help of molecular docking, Garcinia phytochemicals found to be bound deeply inside the active site of DENV NS2B/NS3 protease among all tested phytochemicals and had interactions with catalytic triad (His51, Asp75, Ser135). Thus, it can be concluded from the study that these Gracinia phytochemicals could serve as important inhibitors to inhibit the viral replication inside the host cell. Further in-vitro investigations require confirming their efficacy.     

Discovery of novel P2 substituted 4-biaryl proline inhibitors of hepatitis C virus NS3 serine protease

Inhibitors of hepatitis C virus NS3 serine protease often incorporate a large P2 moiety to interact with the surface of the enzyme while shielding part of the catalytic triad. This feature is important in many inhibitors in order to have the necessary potency needed for efficacy. In this Letter we explore some new P2 motifs to further exploit this region of the enzyme. In a continuing effort to replace the often found 4-hydroxyproline P2 core found in the majority of inhibitors for this target, various directly attached aryl derivatives were evaluated. Of these, the 2,4-disubstituted thiazole core proved to be the most interesting. SAR around this motif has lead to compounds with K_i’s in the high picomolar range and provided cellular potencies in the single digit nM range.     

Structure‐guided screening of chemical database to identify NS3‐NS2B inhibitors for effective therapeutic application in dengue infection

Dengue infection is the most common arthropod‐borne disease caused by dengue viruses, predominantly affecting millions of human beings annually. To find out promising chemical entities for therapeutic application in Dengue, in the current research, a multi‐step virtual screening effort was conceived to screen out the entire “screening library” of the Asinex database. Initially, through “Lipinski rule of five” filtration criterion almost 0.6 million compounds were collected and docked with NS3‐NS2B protein. Thereby, the chemical space was reduced to about 3500 compounds through the analysis of binding affinity obtained from molecular docking study in AutoDock Vina. Further, the “Virtual Screening Workflow” (VSW) utility of Schrödinger suite was used, which follows a stepwise multiple docking programs such as ‐ high‐throughput virtual screening (HTVS), standard precision (SP), and extra precision (XP) docking, and in postprocessing analysis the MM‐GBSA based free binding energy calculation. Finally, five potent molecules were proposed as potential inhibitors for the dengue NS3‐NS2B protein based on the investigation of molecular interactions map and protein‐ligand fingerprint analyses. Different pharmacokinetics and drug‐likeness parameters were also checked, which favour the potentiality of selected molecules for being drug‐like candidates. The molecular dynamics (MD) simulation analyses of protein‐ligand complexes were explained that NS3‐NS2B bound with proposed molecules quite stable in dynamic states as observed from the root means square deviation (RMSD) and root means square fluctuation (RMSF) parameters. The binding free energy was calculated using MM‐GBSA method from the MD simulation trajectories revealed that all proposed molecules possess such a strong binding affinity towards the dengue NS3‐NS2B protein. Therefore, proposed molecules may be potential chemical components for effective inhibition of dengue NS3‐NS2B protein subjected to experimental validation.     

Flavonoid from Carica papaya inhibits NS2B-NS3 protease and prevents Dengue 2 viral assembly

Dengue virus belongs to the virus family Flaviviridae. Dengue hemorrhagic disease caused by dengue virus is a public health problem worldwide. The viral non structural 2B and 3 (NS2B-NS3) protease complex is crucial for virus replication and hence, it is considered to be a good anti-viral target. Leaf extracts from Carica papaya is generally prescribed for patients with dengue fever, but there are no scientific evidences for its anti-dengue activity; hence we intended to investigate the anti-viral activity of compounds present in the leaves of Carica papaya against dengue 2 virus (DENV-2). We analysed the anti-dengue activities of the extracts from Carica papaya by using bioinformatics tools. Interestingly, we find the flavonoid quercetin with highest binding energy against NS2B-NS3 protease which is evident by the formation of six hydrogen bonds with the amino acid residues at the binding site of the receptor. Our results suggest that the flavonoids from Carica papaya have significant anti-dengue activities.

Abbreviations

ADME - Absorption, distribution, metabolism and excretion, BBB - Blood brain barrier, CYP - Cytochrome P450, DENV - – Dengue virus, DHF - Dengue hemorrhagic fever, DSS - Dengue shock syndrome, GCMS - – Gas chromatography- Mass spectrometry, MOLCAD - Molecular Computer Aided Design, NS - Non structural, PDB - Protein data bank, PMF - Potential Mean Force.