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

Synthesis and biological activity of imidazole phenazine derivatives as potential inhibitors for NS2B-NS3 dengue protease

A series of new imidazole-phenazine derivatives were synthesized via a two-step process. The condensation of 2,3-diaminophenazine and benzaldehyde derivatives proceeds with intermediate formation of an aniline Schiff base, which undergoes subsequent cyclodehydrogenation in situ. The structures of the synthesized compounds were characterized by 1D and 2D NMR, FTIR and HRMS. A total of thirteen imidazole phenazine derivatives were synthesized and validated for their inhibitory activity as anti-dengue agents by an in vitro DENV2 NS2B-NS3 protease assay using a fluorogenic Boc-Gly-Arg-Arg-AMC substrate. Two para-substituted imidazole phenazines, 3e and 3k, were found to be promising lead molecules for novel NS2B-NS3 protease inhibitors with IC50 of 54.8 μM and 71.9 μM, respectively, compared to quercetin as a control (IC50 104.8 μM). The in silico study was performed using AutoDock Vina to identify the binding energy and conformation of 3e and 3k with the active site of the DENV2 NS2B-NS3 protease Wichapong model. The results indicate better binding properties of 3e and 3k with calculated binding energies of -8.5 and -8.4 kcal mol-1, respectively, compared to the binding energy of quercetin of -7.2 kcal mol-1, which corroborates well with the experimental observations.

Computational study of nitro-benzylidene phenazine as dengue virus-2 NS2B-NS3 protease inhibitor

According to the World Health Organisation (WHO), as of week 23 of 2022, there were more than 1,311 cases of dengue in Malaysia, with 13 deaths reported. Furthermore, there was an increase of 65.7% during the same period in 2021. Despite the increase in cumulative dengue incidence, there is no effective antiviral drug available for dengue treatment. This work aimed to evaluate several nitro-benzylidene phenazine compounds, especially those that contain 4-hydroxy-3,5-bis((2-(4-nitrophenyl)hydrazinylidene)-methyl)benzoate through pharmacophore queries selection method as potential dengue virus 2 (DENV2) NS2B-NS3 protease inhibitors. Herein, molecular docking was employed to correlate the energies of selected hits' free binding and their binding affinities. Pan assay interference compounds (PAINS) filter was also adopted to identify and assess the drug-likeness, toxicity, mutagenicity potentials, and pharmacokinetic profiles to select hit compounds that can be considered as lead DENV2 NS2B-NS3 protease inhibitors. Molecular dynamics assessment of two nitro-benzylidene phenazine derivatives bearing dinitro and hydroxy groups at the benzylidene ring showed their stability at the main binding pocket of DENV2 protease, where their MM-PBSA binding energies were between -22.53 and -17.01 kcal/mol. This work reports those two nitro-benzylidene phenazine derivatives as hits with 52-55% efficiency as antiviral candidates. Therefore, further optimisation is required to minimise the lead compounds' toxicity and mutagenicity.

Potential Role of Flavivirus NS2B-NS3 Proteases in Viral Pathogenesis and Anti-flavivirus Drug Discovery Employing Animal Cells and Models: A Review

Flaviviruses are known to cause a variety of diseases in humans in different parts of the world. There are very limited numbers of antivirals to combat flavivirus infection, and therefore new drug targets must be explored. The flavivirus NS2B-NS3 proteases are responsible for the cleavage of the flavivirus polyprotein, which is necessary for productive viral infection and for causing clinical infections; therefore, they are a promising drug target for devising novel drugs against different flaviviruses. This review highlights the structural details of the NS2B-NS3 proteases of different flaviviruses, and also describes potential antiviral drugs that can interfere with the viral protease activity, as determined by various studies. Moreover, optimized in vitro reaction conditions for studying the NS2B-NS3 proteases of different flaviviruses may vary and have been incorporated in this review. The increasing availability of the in silico and crystallographic/structural details of flavivirus NS2B-NS3 proteases in free and drug-bound states can pave the path for the development of promising antiflavivirus drugs to be used in clinics. However, there is a paucity of information available on using animal cells and models for studying flavivirus NS2B-NS3 proteases, as well as on the testing of the antiviral drug efficacy against NS2B-NS3 proteases. Therefore, on the basis of recent studies, an effort has also been made to propose potential cellular and animal models for the study of flavivirus NS2B-NS3 proteases for the purposes of exploring flavivirus pathogenesis and for testing the efficacy of possible drugs targets, in vitro and in vivo.

Study of combining virtual screening and antiviral treatments of the Sars-CoV-2 (Covid-19)

The recent epidemic outbreak of a novel human coronavirus called SARS-CoV-2 and causing the respiratory tract disease COVID-19 has reached worldwide resonance and a global effort is being undertaken to characterize the molecular features and evolutionary origins of this virus. Therefore, rapid and accurate identification of pathogenic viruses plays a vital role in selecting appropriate treatments, saving people's lives and preventing epidemics. Additionally, general treatments, coronavirus-specific treatments, and antiviral treatments useful in fighting COVID-19 are addressed. This review sets out to shed light on the SARS-CoV-2 and host receptor recognition, a crucial factor for successful virus infection and taking immune-informatics approaches to identify B- and T-cell epitopes for surface glycoprotein of SARS-CoV-2. A variety of improved or new approaches also have been developed. It is anticipated that this will assist researchers and clinicians in developing better techniques for timely and effective detection of coronavirus infection. Moreover, the genomic sequence of the virus responsible for COVID-19, as well as the experimentally determined three-dimensional structure of the Main protease (Mpro) is available. The reported structure of the target Mpro was described in this review to identify potential drugs for COVID-19 using virtual high throughput screening.

Structural Genomics of SARS-CoV-2 Indicates Evolutionary Conserved Functional Regions of Viral Proteins

During its first two and a half months, the recently emerged 2019 novel coronavirus, SARS-CoV-2, has already infected over one-hundred thousand people worldwide and has taken more than four thousand lives. However, the swiftly spreading virus also caused an unprecedentedly rapid response from the research community facing the unknown health challenge of potentially enormous proportions. Unfortunately, the experimental research to understand the molecular mechanisms behind the viral infection and to design a vaccine or antivirals is costly and takes months to develop. To expedite the advancement of our knowledge, we leveraged data about the related coronaviruses that is readily available in public databases and integrated these data into a single computational pipeline. As a result, we provide comprehensive structural genomics and interactomics roadmaps of SARS-CoV-2 and use this information to infer the possible functional differences and similarities with the related SARS coronavirus. All data are made publicly available to the research community.

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.

Thioguanine-based DENV-2 NS2B/NS3 protease inhibitors: Virtual screening, synthesis, biological evaluation and molecular modelling

Dengue virus Type 2 (DENV-2) is predominant serotype causing major dengue epidemics. There are a number of studies carried out to find its effective antiviral, however to date, there is still no molecule either from peptide or small molecules released as a drug. The present study aims to identify small molecules inhibitor from National Cancer Institute database through virtual screening. One of the hits, D0713 (IC₅₀ = 62 μM) bearing thioguanine scaffold was derivatised into 21 compounds and evaluated for DENV-2 NS2B/NS3 protease inhibitory activity. Compounds 18 and 21 demonstrated the most potent activity with IC₅₀ of 0.38 μM and 16 μM, respectively. Molecular dynamics and MM/PBSA free energy of binding calculation were conducted to study the interaction mechanism of these compounds with the protease. The free energy of binding of 18 calculated by MM/PBSA is -16.10 kcal/mol compared to the known inhibitor, panduratin A (-11.27 kcal/mol), which corroborates well with the experimental observation. Results from molecular dynamics simulations also showed that both 18 and 21 bind in the active site and stabilised by the formation of hydrogen bonds with Asn174.

Development of a NS2B/NS3 protease inhibition assay using AlphaScreen® beads for screening of anti-dengue activities

Background

Dengue infection is an endemic infectious disease and it can lead to dengue fever, dengue hemorrhagic fever, and/or dengue shock syndromes. Dengue NS2B/NS3 protease complex is essential for viral replication and is a primary target for anti-dengue drug development. In this study, a NS2B/NS3 protease inhibition assay was developed using AlphaScreen^® beads and was used to screen compounds for their protease inhibition activities.

Methods

The assay system utilized a known NS2B/NS3 peptide substrate, a recombinant of NS2B/NS3 protease with proprietary StrepTactin^® donor and nickel chelate acceptor beads in 384-well format.

Results

The optimized assay to screen for NS2B/NS3 protease inhibitors was demonstrated to be potentially useful with reasonable zʹ factor, coefficient variance and signal to background ratio. However, screening of synthesized thioguanine derivatives using the optimized AlphaScreen^® assay revealed weak NS2B/NS3 inhibition activities.

Conclusion

The AlphaScreen^® assay to screen for NS2B/NS3 protease inhibitors is potentially applicable for high throughput screening.

Evaluation of the Leaf Essential Oil from Artemisia vulgaris and Its Larvicidal and Repellent Activity against Dengue Fever Vector Aedes aegypti-An Experimental and Molecular Docking Investigation

Aedes aegypti is a mosquito vector that spreads dengue fever and yellow fever worldwide in tropical and subtropical countries. Essential oil isolated from Artemisia vulgaris is found to have larvicidal and repellent action against this vector. The dried leaves were subjected to hydrodistillation using a clevenger-type apparatus for 4 h. The isolated essential oil was analyzed by using gas chromatography-mass spectrometry, and the major insecticidal compounds were identified as α-humulene (0.72%), β-caryophyllene (0.81%), and caryophyllene oxide (15.87%). Larvicidal activity results revealed that the essential oil exposure for 24 h period against the third stage larvae was LC50 = 6.87, LC90 = 59.197 ppm and for the fourth stage larvae LC50 = 4.269, LC90 = 50.363 ppm. Highest mortality rates were observed at 24 h exposure period of third and fourth stages, and the exposed A. aegypti larvae were subjected to histo chemical studies, and the studies revealed that larvae cells got totally damaged (midgut and cortex). The repellent activity results revealed that at 50% concentration of the essential oil showed the highest repellent activity at 60 min protection time against the A. aegypti female mosquitoes. To gain further insights into the insecticidal activity, density functional theory and molecular docking calculations were performed with the active components of this essential oil as the ligand and NS3 protease domain (PDB ID: 2FOM) as a receptor. Molecular docking calculation results show that (E)-β-caryophyllene strongly binds with NS3 protease domain than (Z)-β-caryophyllene, α-humulene, and β-caryophyllene oxide and is the major active component for the insecticidal action. It primarily interacts with the receptor through hydrophobic and ionic forces and using water bridges between the amino acid residues in the binding pocket and (E)-β-caryophyllene.

Molecular Dynamics Assisted Mechanistic Study of Isoniazid-Resistance against Mycobacterium tuberculosis InhA

Examination of InhA mutants I16T, I21V, I47T, S94A, and I95P showed that direct and water mediated H-bond interactions between NADH and binding site residues reduced drastically. It allowed conformational flexibility to NADH, particularly at the pyrophosphate region, leading to weakening of its binding at dinucleotide binding site. The highly scattered distribution of pyrophosphate dihedral angles and chi1 side chain dihedral angles of corresponding active site residues therein confirmed weak bonding between InhA and NADH. The average direct and water mediated bridged H-bond interactions between NADH and mutants were observed weaker as compared to wild type. Further, estimated NADH binding free energy in mutants supported the observed weakening of InhA-NADH interactions. Similarly, per residue contribution to NADH binding was also found little less as compared to corresponding residues in wild type. This investigation clearly depicted and supported the effect of mutations on NADH binding and can be accounted for isoniazid resistance as suggested by previous biochemical and mutagenic studies. Further, structural analysis of InhA provided the crucial points to enhance the NADH binding affinity towards InhA mutants in the presence of direct InhA inhibitors to combat isoniazid drug resistance. This combination could be a potential alternative for treatment of drug resistant tuberculosis.

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.

Inhibitors of Dengue virus and West Nile virus proteases based on the aminobenzamide scaffold

Dengue and West Nile viruses (WNV) are mosquito-borne members of flaviviruses that cause significant morbidity and mortality. There is no approved vaccine or antiviral drugs for human use to date. In this study, a series of functionalized meta and para aminobenzamide derivatives were synthesized and subsequently screened in vitro against Dengue virus and West Nile virus proteases. Four active compounds were identified which showed comparable activity toward the two proteases and shared in common a meta or para(phenoxy)phenyl group. The inhibition constants (K(i)) for the most potent compound 7n against Dengue and West Nile virus proteases were 8.77 and 5.55 μM, respectively. The kinetics data support a competitive mode of inhibition of both proteases by compound 7n. This conclusion is further supported by molecular modeling. This study reveals a new chemical scaffold which is amenable to further optimization to yield potent inhibitors of the viral proteases via the combined utilization of iterative medicinal chemistry/structure-activity relationship studies and in vitro screening.

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.

Inhibition of Dengue virus and West Nile virus proteases by click chemistry-derived benz[d]isothiazol-3(2H)-one derivatives

Two click chemistry-derived focused libraries based on the benz[d]isothiazol-3(2H)-one scaffold were synthesized and screened against Dengue virus and West Nile virus NS2B-NS3 proteases. Several compounds (4l, 7j-n) displayed noteworthy inhibitory activity toward Dengue virus NS2B-NS3 protease in the absence and presence of added detergent. These compounds could potentially serve as a launching pad for a hit-to-lead optimization campaign.

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.

Virtual ligand screening of the National Cancer Institute (NCI) compound library leads to the allosteric inhibitory scaffolds of the West Nile Virus NS3 proteinase

Viruses of the genus Flavivirus are responsible for significant human disease and mortality. The N-terminal domain of the flaviviral nonstructural (NS)3 protein codes for the serine, chymotrypsin-fold proteinase (NS3pro). The presence of the nonstructural (NS)2B cofactor, which is encoded by the upstream gene in the flaviviral genome, is necessary for NS3pro to exhibit its proteolytic activity. The two-component NS2B-NS3pro functional activity is essential for the viral polyprotein processing and replication. Both the structure and the function of NS2B-NS3pro are conserved in the Flavivirus family. Because of its essential function in the posttranslational processing of the viral polyprotein precursor, NS2B-NS3pro is a promising target for anti-flavivirus drugs. To identify selective inhibitors with the reduced cross-reactivity and off-target effects, we focused our strategy on the allosteric inhibitors capable of targeting the NS2B-NS3pro interface rather than the NS3pro active site. Using virtual ligand screening of the diverse, ∼275,000-compound library and the catalytic domain of the two-component West Nile virus (WNV) NS2B-NS3pro as a receptor, we identified a limited subset of the novel inhibitory scaffolds. Several of the discovered compounds performed as allosteric inhibitors and exhibited a nanomolar range potency in the in vitro cleavage assays. The inhibitors were also potent in cell-based assays employing the sub-genomic, luciferase-tagged WNV and Dengue viral replicons. The selectivity of the inhibitors was confirmed using the in vitro cleavage assays with furin, a human serine proteinase, the substrate preferences of which are similar to those of WNV NS2B-NS3pro. Conceptually, the similar in silico drug discovery strategy may be readily employed for the identification of inhibitors of other flaviviruses.

Novel dengue virus-specific NS2B/NS3 protease inhibitor, BP2109, discovered by a high-throughput screening assay

Dengue virus (DENV) causes disease globally, with an estimated 25 to 100 million new infections per year. At present, no effective vaccine is available, and treatment is supportive. In this study, we identified BP2109, a potent and selective small-molecule inhibitor of the DENV NS2B/NS3 protease, by a high-throughput screening assay using a recombinant protease complex consisting of the central hydrophilic portion of NS2B and the N terminus of the protease domain. BP2109 inhibited DENV (serotypes 1 to 4), but not Japanese encephalitis virus (JEV), replication and viral RNA synthesis without detectable cytotoxicity. The compound inhibited recombinant DENV-2 NS2B/NS3 protease with a 50% inhibitory concentration (IC(50)) of 15.43 ± 2.12 μM and reduced the reporter expression of the DENV-2 replicon with a 50% effective concentration (EC(50)) of 0.17 ± 0.01 μM. Sequencing analyses of several individual clones derived from BP2109-resistant DENV-2 RNAs revealed that two amino acid substitutions (R55K and E80K) are found in the region of NS2B, a cofactor of the NS2B/NS3 protease complex. The introduction of R55K and E80K double mutations into the dengue virus NS2B/NS3 protease and a dengue virus replicon construct conferred 10.3- and 73.8-fold resistance to BP2109, respectively. The E80K mutation was further determined to be the key mutation conferring dengue virus replicon resistance (61.3-fold) to BP2109, whereas the R55K mutation alone did not affect resistance to BP2109. Both the R55K and E80K mutations are located in the central hydrophilic portion of the NS2B cofactor, where extensive interactions with the NS3pro domain exist. Thus, our data provide evidence that BP2109 likely inhibits DENV by a novel mechanism.