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

ZIKV Inhibitors Based on Pyrazolo[3,4-d]pyridazine-7-one Core: Rational Design, In Vitro Evaluation, and Theoretical Studies

The Zika virus (ZIKV) is believed to cause birth defects, and no anti-ZIKV drugs have been approved by medical organizations to date. Starting from antimicrobial lead compounds with a pyrazolo[3,4-d]pyridazine-7-one scaffold, we synthesized 16 derivatives and screened their ability to interfere with ZIKV infection utilizing a cell-based phenotypic assay. Of these, five compounds showed significant inhibition of ZIKV with a selective index value greater than 4.6. In particular, compound 9b showed the best anti-ZIKV activity with a selectivity index of 22.4 (half-maximal effective concentration = 25.6 μM and 50% cytotoxic concentration = 572.4 μM). Through the brine shrimp lethality bioassay, 9b, 10b, 12, 17a, and 19a showed median lethal dose values in a range of 87.2-100.3 μg/mL. Compound 9b was also targeted to the NS2B-NS3 protease of ZIKV using molecular docking protocols, in which it acted as a noncompetitive inhibitor and strongly bound to five key amino acids (His51, Asp75, Ser135, Ala132, Tyr161). Utilizing the pharmacophore model of 9b, the top 20 hits were identified as prospective inhibitors of NS2B-NS3 protease, and six of them were confirmed for their stability with the protease via redocking and molecular dynamics simulations.

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.

Quantum Biochemistry and MM-PBSA Description of the ZIKV NS2B-NS3 Protease: Insights into the Binding Interactions beyond the Catalytic Triad Pocket

The Zika virus protease NS2B-NS3 has a binding site formed with the participation of a H51-D75-S135 triad presenting two forms, active and inactive. Studies suggest that the inactive conformation is a good target for the design of inhibitors. In this paper, we evaluated the co-crystallized structures of the protease with the inhibitors benzoic acid (5YOD) and benzimidazole-1-ylmethanol (5H4I). We applied a protocol consisting of two steps: first, classical molecular mechanics energy minimization followed by classical molecular dynamics were performed, obtaining stabilized molecular geometries; second, the optimized/relaxed geometries were used in quantum biochemistry and molecular mechanics/Poisson-Boltzmann surface area (MM-PBSA) calculations to estimate the ligand interactions with each amino acid residue of the binding pocket. We show that the quantum-level results identified essential residues for the stabilization of the 5YOD and 5H4I complexes after classical energy minimization, matching previously published experimental data. The same success, however, was not observed for the MM-PBSA simulations. The application of quantum biochemistry methods seems to be more promising for the design of novel inhibitors acting on NS2B-NS3.

C-Terminal Extended Hexapeptides as Potent Inhibitors of the NS2B-NS3 Protease of the ZIKA Virus

The Zika virus (ZIKV) protease is an attractive drug target for the design of novel inhibitors to control the ZIKV infection. As the protease substrate-binding site contains acidic residues, inhibitors with basic residues can be beneficial for the inhibition of protease activities. Molecular dynamics (MD) simulation and molecular mechanics with generalized Born and surface area solvation (MM/GBSA) techniques are employed herein to design potent peptide inhibitors and to understand the nature of the basic residues that can potentially stabilize the acidic residues of the protease substrate-binding site. It is found that the inclusion of K, R, and K at P1, P2, and P3 positions, respectively, and Y at the P4 position (YKRK) would generate a highly stable tetrapeptide-protease complex with a ΔG_bind of ~ −80 kcal/mol. We have also shown that the C-terminal extension of this and the second most stable tetrapeptide (YRRR) with small polar residues, such as S and T would generate even more stable hexapeptide-protease complexes. The modes of interactions of these inhibitors are discussed in detail, which are in agreement with earlier experimental studies. Thus, this study is expected to aid in the design of novel antiviral drugs against the ZIKV.

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.

Intermolecular interactions of cn-716 and acyl-KR-aldehyde dipeptide inhibitors against Zika virus

Structural representation and graphic panel showing the most relevant residues that contribute to the ZIKV NS2B–NS3–ligand complexes.