Structurally- and dynamically-driven allostery of the chymotrypsin-like proteases of SARS, Dengue and Zika viruses

Duck Tembusu virus NS3 protein induces apoptosis by activating the PERK/PKR pathway and mitochondrial pathway

IMPORTANCE

Duck Tembusu virus (DTMUV) is an emerging pathogenic flavivirus that replicates well in mosquito, bird, and mammalian cells. An in vivo study revealed that BALB/c mice and Kunming mice were susceptible to DTMUV after intracerebral inoculation. Moreover, there are no reports about DTMUV-related human disease, but antibodies against DTMUV and viral RNA were detected in the serum samples of duck industry workers. This information implies that DTMUV has expanded its host range and poses a threat to mammalian health. Thus, understanding the pathogenic mechanism of DTMUV is crucial for identifying potential antiviral targets. In this study, we discovered that NS3 can induce the mitochondria-mediated apoptotic pathway through the PERK/PKR pathway; it can also interact with voltage-dependent anion channel 2 to induce apoptosis. Our findings provide a theoretical basis for understanding the pathogenic mechanism of DTMUV infection and identifying potential antiviral targets and may also serve as a reference for exploring the pathogenesis of other flaviviruses.

Structural Modifications Introduced by NS2B Cofactor Binding to the NS3 Protease of the Kyasanur Forest Disease Virus

Kyasanur Forest Disease virus (KFDV), a neglected human pathogenic virus, is a Flavivirus that causes severe hemorrhagic fever in humans. KFDV is transmitted to humans by the bite of the hard tick (Haemaphysalis spinigera), which acts as a reservoir of KFDV. The recent expansion of the endemic area of KFDV is of concern and requires the development of new preventive measures against KFDV. Currently, there is no antiviral therapy against KFDV, and the existing vaccine has limited efficacy. To develop a new antiviral therapy against KFDV, we focused on the nonstructural proteins NS2B and NS3 of KFDV, which are responsible for serine protease activity. Viral proteases have shown to be suitable therapeutic targets in the development of antiviral drugs against many diseases. However, success has been limited in flaviviruses, mainly because of the important features of the active site, which is flat and highly charged. In this context, the present study focuses on the dynamics of NS2B and NS3 to identify potential allosteric sites in the NS2B/NS3 protease of KDFV. To our knowledge, there are no reports on the dynamics of NS2B and NS3 in KFDV, and the crystal structure of the NS2B/NS3 protease of KFDV has not yet been solved. Overall, we created the structure of the NS2B/NS3 protease of KFDV using AlphaFold and performed molecular dynamics simulations with and without NS2B cofactor to investigate structural rearrangements due to cofactor binding and to identify alternative allosteric sites. The identified allosteric site is promising due to its geometric and physicochemical properties and druggability and can be used for new drug development. The applicability of the proposed allosteric binding sites was verified for the best-hit molecules from the virtual screening and MD simulations.

Investigating the folding dynamics of NS2B protein of Zika virus

NS2B protein of the Zika virus acts as a co-factor for NS3 protease and also involves in remodeling NS3 protease structure. Therefore, we investigated the overall dynamics of NS2B protein. We find surprising similarities between selected flavivirus NS2B model structures predicted from Alphafold2. Further, the simulated ZIKV NS2B protein structure shows a disordered cytosolic domain (residues 45-95) as a part of a full-length protein. Since only the cytosolic domain of NS2B is sufficient for the protease activity, we also investigated the conformational dynamics of only ZIKV NS2B cytosolic domain (residues 49-95) in the presence of TFE, SDS, Ficoll, and PEG using simulation and spectroscopy. The presence of TFE induces α-helix in NS2B cytosolic domain (residues 49-95). On the other hand, the presence of SDS, ficoll, and PEG does not induce secondary structural change. This dynamics study could have implications for some unknown folds of the NS2B protein.

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.

Myricetin Allosterically Inhibits the Dengue NS2B-NS3 Protease by Disrupting the Active and Locking the Inactive Conformations

The dengue NS2B-NS3 protease existing in equilibrium between the active and inactive forms is essential for virus replication, thus representing a key drug target. Here, myricetin, a plant flavonoid, was characterized to noncompetitively inhibit the dengue protease. Further NMR study identified the protease residues perturbed by binding to myricetin, which were utilized to construct the myricetin-protease complexes. Strikingly, in the active form, myricetin binds to a new allosteric site (AS2) far away from the active site pocket and the allosteric site (AS1) for binding curcumin, while in the inactive form, it binds to both AS1 and AS2. To decipher the mechanism for the allosteric inhibition by myricetin, we conducted molecular dynamics simulations on different forms of dengue NS2B-NS3 proteases. Unexpectedly, the binding of myricetin to AS2 is sufficient to disrupt the active conformation by displacing the characteristic NS2B C-terminal β-hairpin from the active site pocket. By contrast, the binding of myricetin to AS1 and AS2 results in locking the inactive conformation. Therefore, myricetin represents the first small molecule, which allosterically inhibits the dengue protease by both disrupting the active conformation and locking the inactive conformation. The results enforce the notion that a global allosteric network exists in the dengue NS2B-NS3 protease, which is susceptible to allosteric inhibition by small molecules such as myricetin and curcumin. As myricetin has been extensively used as a food additive, it might be directly utilized to fight the dengue infections and as a promising starting material for further design of potent allosteric inhibitors.

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.

Impact of dimerization and N3 binding on molecular dynamics of SARS-CoV and SARS-CoV-2 main proteases

SARS-CoV-2 main protease is one of the major targets in drug development efforts against Covid-19. Even though several structures were reported to date, its dynamics is not understood well. In particular, impact of dimerization and ligand binding on the dynamics is an important issue to investigate. In this study, we performed molecular dynamics simulations of SARS-CoV and SARS-CoV-2 main proteases to investigate influence of dimerization on the dynamics by modeling monomeric and dimeric apo and holo forms. The dimerization causes an organization of the interdomain dynamics as well as some local structural changes. Moreover, we investigated impact of a peptide mimetic (N3) on the dynamics of SARS-CoV and SARS-CoV-2 Mpro. The ligand binding to the dimeric forms causes some key local changes at the dimer interface and it causes an allosteric interaction between the active sites of two protomers. Our results support the idea that only one protomer is active on SARS-CoV-2 due to this allosteric interaction. Additionally, we analyzed the molecular dynamics trajectories from graph theoretical perspective and found that the most influential residues – as measured by eigenvector centrality – are a group of residues in active site and dimeric interface of the protease. This study may form a bridge between what we know about the dynamics of SARS-CoV and SARS-CoV-2 Mpro. We think that enlightening allosteric communication of the active sites and the role of dimerization in SARS-CoV-2 Mpro can contribute to development of novel drugs against this global health problem as well as other similar proteases.

Communicated by Ramaswamy H. Sarma

Overview of antiviral drug candidates targeting coronaviral 3C-like main proteases

Coronaviruses (CoVs) are positive single-stranded RNA viruses that cause severe respiratory syndromes in humans, including severe acute respiratory syndrome (SARS) and Middle East respiratory syndrome (MERS). Coronavirus disease 2019 (COVID-19) caused by a novel severe acute respiratory syndrome CoV (SARS-CoV-2) at the end of 2019 became a global pandemic. The 3C-like cysteine protease (3CLpro) processes viral polyproteins to yield mature non-structural proteins, thus playing an important role in the CoV life cycle, and therefore is considered as a prominent target for antiviral drugs. To date, many 3CLpro inhibitors have been reported, and their molecular mechanisms have been illustrated. Here, we briefly introduce the structural features of 3CLpro of the human-related SARS-CoV, MERS-CoV and SARS-CoV-2, and explore the potency and mechanism of their cognate inhibitors. This information will shed light on the development and optimization of CoV 3CLpro inhibitors, which may benefit the further designation of therapeutic strategies for treating CoV diseases.

Experiments and simulation on ZIKV NS2B-NS3 protease reveal its complex folding

Zika virus has been identified in various body fluids such as semen, urine, saliva, cerebrospinal fluid, and vaginal secretion of an infected individual. The pH of these fluids varies from mildly acidic to mildly alkaline. So it is imperative to understand the impact of these conditions on viral protein functioning. We investigated the NS2B-NS3 protease stability and its activity in different denaturing environments. Finding indicates that NS2B-NS3 protease maintains stability at pH 4.8-8.7. Thus it suggests that the complex remains functionally active to hydrolyze the polyprotein within a diverse environmental condition such as variable pH. Despite a stable structure at a broad pH range, a change in environmental conditions dramatically influence its protease activity. Moreover, it is susceptible to structural transformation leading to increased β-strand or helix content in the presence of alcohol. This study may help further to understand the folding-function relationship of the general flaviviral protease complex.

Protein structural heterogeneity: A hypothesis for the basis of proteolytic recognition by the main protease of SARS-CoV and SARS-CoV-2

The main protease (M^pro) of SARS-CoV and SARS-CoV-2 is a key enzyme in viral replication and a promising target for the development of antiviral therapeutics. The understanding of this protein is based on a number of observations derived from earlier x-ray structures, which mostly consider substrates or ligands as the main reason behind modulation of the active site. This lead to the concept of substrate-induced subsite cooperativity as an initial attempt to explain the dual binding specificity of this enzyme in recognizing the cleavage sequences at its N- and C-termini, which are important processing steps in obtaining the mature protease. The presented hypothesis proposes that structural heterogeneity is a property of the enzyme, independent of the presence of a substrate or ligand. Indeed, the analysis of M^pro structures of SARS-CoV and SARS-CoV-2 reveals a conformational diversity for the catalytically competent state in ligand-free structures. Variation in the binding site appears to result from flexibility at residues lining the S₁ subpocket and segments incorporating methionine 49 and glutamine 189. The structural evidence introduces “structure-based recognition” as a new paradigm in substrate proteolysis by M^pro. In this concept, the binding space in subpockets of the enzyme varies in a non-cooperative manner, causing distinct conformations, which recognize and process different cleavage sites, as the N- and C-termini. Insights into the recognition basis of the protease provide explanation to the ordered processing of cleavage sites. The hypothesis expands the conformational space of the enzyme and consequently opportunities for drug development and repurposing efforts.

Tethering-induced destabilization and ATP-binding for tandem RRM domains of ALS-causing TDP-43 and hnRNPA1

TDP-43 and hnRNPA1 contain tandemly-tethered RNA-recognition-motif (RRM) domains, which not only functionally bind an array of nucleic acids, but also participate in aggregation/fibrillation, a pathological hallmark of various human diseases including amyotrophic lateral sclerosis (ALS), frontotemporal dementia (FTD), alzheimer's disease (AD) and Multisystem proteinopathy (MSP). Here, by DSF, NMR and MD simulations we systematically characterized stability, ATP-binding and conformational dynamics of TDP-43 and hnRNPA1 RRM domains in both tethered and isolated forms. The results reveal three key findings: (1) upon tethering TDP-43 RRM domains become dramatically coupled and destabilized with Tm reduced to only 49 °C. (2) ATP specifically binds TDP-43 and hnRNPA1 RRM domains, in which ATP occupies the similar pockets within the conserved nucleic-acid-binding surfaces, with the affinity slightly higher to the tethered than isolated forms. (3) MD simulations indicate that the tethered RRM domains of TDP-43 and hnRNPA1 have higher conformational dynamics than the isolated forms. Two RRM domains become coupled as shown by NMR characterization and analysis of inter-domain correlation motions. The study explains the long-standing puzzle that the tethered TDP-43 RRM1–RRM2 is particularly prone to aggregation/fibrillation, and underscores the general role of ATP in inhibiting aggregation/fibrillation of RRM-containing proteins. The results also rationalize the observation that the risk of aggregation-causing diseases increases with aging.

Distinct conformational dynamics and allosteric networks in alpha tryptophan synthase during active catalysis

Experimental observations of enzymes under active turnover conditions have brought new insight into the role of protein motions and allosteric networks in catalysis. Many of these studies characterize enzymes under dynamic chemical equilibrium conditions, in which the enzyme is actively catalyzing both the forward and reverse reactions during data acquisition. We have previously analyzed conformational dynamics and allosteric networks of the alpha subunit of tryptophan synthase under such conditions using NMR. We have proposed that this working state represents a four to one ratio of the enzyme bound with the indole-3-glycerol phosphate substrate (E:IGP) to the enzyme bound with the products indole and glyceraldehyde-3-phosphate (E:indole:G3P). Here, we analyze the inactive D60N variant to deconvolute the contributions of the substrate- and products-bound states to the working state. While the D60N substitution itself induces small structural and dynamic changes, the D60N E:IGP and E:indole:G3P states cannot entirely account for the conformational dynamics and allosteric networks present in the working state. The act of chemical bond breakage and/or formation, or possibly the generation of an intermediate, may alter the structure and dynamics present in the working state. As the enzyme transitions from the substrate-bound to the products-bound state, millisecond conformational exchange processes are quenched and new allosteric connections are made between the alpha active site and the surface which interfaces with the beta subunit. The structural ordering of the enzyme and these new allosteric connections may be important in coordinating the channeling of the indole product into the beta subunit.

The role of cysteine peptidases in coronavirus cell entry and replication: The therapeutic potential of cathepsin inhibitors

Over the last 2 decades, several coronaviruses (CoVs) have crossed the species barrier into humans, causing highly prevalent and severe respiratory diseases, often with fatal outcomes. CoVs are a large group of enveloped, single-stranded, positive-sense RNA viruses, which encode large replicase polyproteins that are processed by viral peptidases to generate the nonstructural proteins (Nsps) that mediate viral RNA synthesis. Papain-like peptidases (PLPs) and chymotrypsin-like cysteine 3C-like peptidase are essential for coronaviral replication and represent attractive antiviral drug targets. Furthermore, CoVs utilize the activation of their envelope spike glycoproteins by host cell peptidases to gain entry into cells. CoVs have evolved multiple strategies for spike protein activation, including the utilization of lysosomal cysteine cathepsins. In this review, viral and host peptidases involved in CoV cell entry and replication are discussed in depth, with an emphasis on papain-like cysteine cathepsins. Furthermore, important findings on cysteine peptidase inhibitors with regard to virus attenuation are highlighted as well as the potential of such inhibitors for future treatment strategies for CoV-related diseases.

Curcumin Allosterically Inhibits the Dengue NS2B-NS3 Protease by Disrupting Its Active Conformation

Flaviviruses including dengue virus and Zika virus encode a unique two-component NS2B-NS3 protease essential for maturation/infectivity, thus representing a key target for designing antiflavivirus drugs. Here, for the first time, by NMR and molecular docking, we reveal that curcumin allosterically inhibits the dengue protease by binding to a cavity with no overlap with the active site. Further molecular dynamics simulations decode that the binding of curcumin leads to unfolding/displacing the characteristic β-hairpin of the C-terminal NS2B and consequently disrupting the closed (active) conformation of the protease. Our study identified a cavity most likely conserved in all flaviviral NS2B-NS3 proteases, which could thus serve as a therapeutic target for the discovery/design of small-molecule allosteric inhibitors. Moreover, as curcumin has been used as a food additive for thousands of years in many counties, it can be directly utilized to fight the flaviviral infections and as a promising starting for further design of potent allosteric inhibitors.

Potential of coronavirus 3C-like protease inhibitors for the development of new anti-SARS-CoV-2 drugs: Insights from structures of protease and inhibitors

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), similar to SARS-CoV and Middle East respiratory syndrome coronavirus (MERS-CoV), which belong to the same Betacoronavirus genus, induces severe acute respiratory disease that is a threat to human health. Since the outbreak of infection by SARS-CoV-2 began, which causes coronavirus disease 2019 (COVID-19), the disease has rapidly spread worldwide. Thus, a search for effective drugs able to inhibit SARS-CoV-2 has become a global pursuit. The 3C-like protease (3CL^pro), which hydrolyses viral polyproteins to produce functional proteins, is essential for coronavirus replication and is considered an important therapeutic target for diseases caused by coronaviruses, including COVID-19. Many 3CL^pro inhibitors have been proposed and some new drug candidates have achieved success in preclinical studies. In this review, we briefly describe recent developments in determining the structure of 3CL^pro and its function in coronavirus replication and summarise new insights into 3CL^pro inhibitors and their mechanisms of action. The clinical application prospects and limitations of 3CL^pro inhibitors for COVID-19 treatment are also discussed.