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

Elucidating the inhibitory mechanism of Zika virus NS2B-NS3 protease with dipeptide inhibitors: Insights from molecular docking and molecular dynamics simulations

Microcephaly, Guillain-Barré syndrome, and potential sexual transmission stand as prominent complications associated with Zika virus (ZIKV) infection. The absence of FDA-approved drugs or vaccines presents a substantial obstacle in combatting the virus. Furthermore, the inclusion of pregnancy in the pharmacological screening process complicates and extends the endeavor to ensure molecular safety and minimal toxicity. Given its pivotal role in viral assembly and maturation, the NS2B-NS3 viral protease emerges as a promising therapeutic target against ZIKV. In this context, a dipeptide inhibitor was specifically chosen as a control against 200 compounds for docking analysis. Subsequent molecular dynamics simulations extending over 200 ns were conducted to ascertain the stability of the docked complex and confirm the binding of the inhibitor at the protein's active site. The simulation outcomes exhibited conformity to acceptable thresholds, encompassing parameters such as root mean square deviation (RMSD), root mean square fluctuation (RMSF), ligand-protein interaction analysis, ligand characterization, and surface area analysis. Notably, analysis of ligand angles bolstered the identification of prospective ligands capable of inhibiting viral protein activity and impeding virus dissemination. In this study, the integration of molecular docking and dynamics simulations has pinpointed the dipeptide inhibitor as a potential candidate ligand against ZIKV protease, thereby offering promise for therapeutic intervention against the virus.

Investigating the antiviral therapeutic potentialities of marine polycyclic lamellarin pyrrole alkaloids as promising inhibitors for SARS-CoV-2 and Zika main proteases (Mpro)

The new coronavirus variant (SARS-CoV-2) and Zika virus are two world-wide health pandemics. Along history, natural products-based drugs have always crucially recognized as a main source of valuable medications. Considering the SARS-CoV-2 and Zika main proteases (Mpro) as the re-production key element of the viral cycle and its main target, herein we report an intensive computer-aided virtual screening for a focused list of 39 marine lamellarins pyrrole alkaloids, against SARS-CoV-2 and Zika main proteases (Mpro) using a set of combined modern computational methodologies including molecular docking (MDock), molecule dynamic simulations (MDS) and structure-activity relationships (SARs) as well. Indeed, the molecular docking studies had revealed four promising marine alkaloids including [lamellarin H (14)/K (17)] and [lamellarin S (26)/Z (39)], according to their notable ligand-protein energy scores and relevant binding affinities with the SARS-CoV-2 and Zika (Mpro) pocket residues, respectively. Consequentially, these four chemical hits were further examined thermodynamically though investigating their MD simulations at 100 ns, where they showed prominent stability within the accommodated (Mpro) pockets. Moreover, in-deep SARs studies suggested the crucial roles of the rigid fused polycyclic ring system, particularly aromatic A- and F- rings, position of the phenolic -OH and δ-lactone functionalities as essential structural and pharmacophoric features. Finally, these four promising lamellarins alkaloids were investigated for their in-silico ADME using the SWISS ADME platform, where they displayed appropriated drug-likeness properties. Such motivating outcomes are greatly recommending further in vitro/vivo examinations regarding those lamellarins pyrrole alkaloids (LPAs).Communicated by Ramaswamy H. Sarma.

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.

Repurposing of antiparasitic drugs against the NS2B-NS3 protease of the Zika virus

To date, no approved drugs are available to treat the Zika virus (ZIKV) infection. Therefore, it is necessary to urgently identify potential drugs against the ZIKV infection. Here, the repurposing of 30 antiparasitic drugs against the NS2B-NS3 protease of the ZIKV has been carried out by using combined docking and molecular dynamics- (MD) simulations. Based on the docking results, 5 drugs, such as Amodiaquine, Primaquine, Paromomycin, Dichlorophene, and Ivermectin were screened for further analysis by MD simulations and free energy calculations. Among these drugs, Amodiaquine and Dichlorophen are found to produce the most stable complexes and possess relative binding free energies of about -44.3 ± 3.7 kcal/mol and -41.1 ± 5.3 kcal/mol respectively. Therefore, they would act as potent small-molecule inhibitors of the ZIKV protease.However, evaluations of biological and safety activities of these drugs against the ZIKV protease are required before their clinical use.Communicated by Ramaswamy H. Sarma.

A low-background, fluorescent assay to evaluate inhibitors of diverse viral proteases

Multiple coronaviruses (CoVs) can cause respiratory diseases in humans. While prophylactic vaccines designed to prevent infection are available for severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2), incomplete vaccine efficacy, vaccine hesitancy, and the threat of other pathogenic CoVs for which vaccines do not exist have highlighted the need for effective antiviral therapies. While antiviral compounds targeting the viral polymerase and protease are already in clinical use, their sensitivity to potential resistance mutations as well as their breadth against the full range of human and preemergent CoVs remain incompletely defined. To begin to fill that gap in knowledge, we report here the development of an improved, noninfectious, cell-based fluorescent assay with high sensitivity and low background that reports on the activity of viral proteases, which are key drug targets. We demonstrate that the assay is compatible with not only the SARS-CoV-2 Mpro protein but also orthologues from a range of human and nonhuman CoVs as well as clinically reported SARS-CoV-2 drug-resistant Mpro variants. We then use this assay to define the breadth of activity of two clinically used protease inhibitors, nirmatrelvir and ensitrelvir. Continued use of this assay will help define the strengths and limitations of current therapies and may also facilitate the development of next-generation protease inhibitors that are broadly active against both currently circulating and preemergent CoVs. IMPORTANCE Coronaviruses (CoVs) are important human pathogens with the ability to cause global pandemics. Working in concert with vaccines, antivirals specifically limit viral disease in people who are actively infected. Antiviral compounds that target CoV proteases are already in clinical use; their efficacy against variant proteases and preemergent zoonotic CoVs, however, remains incompletely defined. Here, we report an improved, noninfectious, and highly sensitive fluorescent method of defining the sensitivity of CoV proteases to small molecule inhibitors. We use this approach to assay the activity of current antiviral therapies against clinically reported SARS-CoV-2 protease mutants and a panel of highly diverse CoV proteases. Additionally, we show this system is adaptable to other structurally nonrelated viral proteases. In the future, this assay can be used to not only better define the strengths and limitations of current therapies but also help develop new, broadly acting inhibitors that more broadly target viral families.

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.

Discovery of Bispecific Lead Compounds from Azadirachta indica against ZIKA NS2B-NS3 Protease and NS5 RNA Dependent RNA Polymerase Using Molecular Simulations

Zika virus (ZIKV) has been characterized as one of many potential pathogens and placed under future epidemic outbreaks by the WHO. However, a lack of potential therapeutics can result in an uncontrolled pandemic as with other human pandemic viruses. Therefore, prioritized effective therapeutics development has been recommended against ZIKV. In this context, the present study adopted a strategy to explore the lead compounds from Azadirachta indica against ZIKV via concurrent inhibition of the NS2B-NS3 protease (ZIKV^pro) and NS5 RNA dependent RNA polymerase (ZIKV^RdRp) proteins using molecular simulations. Initially, structure-based virtual screening of 44 bioflavonoids reported in Azadirachta indica against the crystal structures of targeted ZIKV proteins resulted in the identification of the top four common bioflavonoids, viz. Rutin, Nicotiflorin, Isoquercitrin, and Hyperoside. These compounds showed substantial docking energy (−7.9 to −11.01 kcal/mol) and intermolecular interactions with essential residues of ZIKV^pro (B:His⁵¹, B:Asp⁷⁵, and B:Ser¹³⁵) and ZIKV^RdRp (Asp⁵⁴⁰, Ile⁷⁹⁹, and Asp⁶⁶⁵) by comparison to the reference compounds, O7N inhibitor (ZIKV^pro) and Sofosbuvir inhibitor (ZIKV^RdRp). Besides, long interval molecular dynamics simulation (500 ns) on the selected docked poses reveals stability of the respective docked poses contributed by intermolecular hydrogen bonds and hydrophobic interactions. The predicted complex stability was further supported by calculated end-point binding free energy using molecular mechanics generalized born surface area (MM/GBSA) method. Consequently, the identified common bioflavonoids are recommended as promising therapeutic inhibitors of ZIKV^pro and ZIKV^RdRp against ZIKV for further experimental assessment.

Structures and Dynamics of Dengue Virus Nonstructural Membrane Proteins

Dengue virus is an important human pathogen threating people, especially in tropical and sub-tropical regions. The viral genome has one open reading frame and encodes one polyprotein which can be processed into structural and nonstructural (NS) proteins. Four of the seven nonstructural proteins, NS2A, NS2B, NS4A and NS4B, are membrane proteins. Unlike NS3 or NS5, these proteins do not harbor any enzymatic activities, but they play important roles in viral replication through interactions with viral or host proteins to regulate important pathways and enzymatic activities. The location of these proteins on the cell membrane and the functional roles in viral replication make them important targets for antiviral development. Indeed, NS4B inhibitors exhibit antiviral activities in different assays. Structural studies of these proteins are hindered due to challenges in crystallization and the dynamic nature of these proteins. In this review, the function and membrane topologies of dengue nonstructural membrane proteins are presented. The roles of solution NMR spectroscopy in elucidating the structure and dynamics of these proteins are introduced. The success in the development of NS4B inhibitors proves that this class of proteins is an attractive target for antiviral development.

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.