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

Alternate recognition by dengue protease: Proteolytic and binding assays provide functional evidence beyond an induced-fit

Proteases are key enzymes in viral replication, and interfering with these targets is the basis for therapeutic interventions. We previously introduced a hypothesis about conformational selection in the protease of dengue virus and related flaviviruses, based on conformational plasticity noted in X-ray structures. The present work presents the first functional evidence for alternate recognition by the dengue protease, in a mechanism based primarily on conformational selection rather than induced-fit. Recognition of distinct substrates and inhibitors in proteolytic and binding assays varies to a different extent, depending on factors reported to influence the protease structure. The pH, salinity, buffer type, and temperature cause a change in binding, proteolysis, or inhibition behavior. Using representative inhibitors with distinct structural scaffolds, we identify two contrasting binding profiles to dengue protease. Noticeable effects are observed in the binding assay upon inclusion of a non-ionic detergent in comparison to the proteolytic assay. The findings highlight the impact of the selection of testing conditions on the observed ligand affinity or inhibitory potency. From a broader scope, the dengue protease presents an example, where the induced-fit paradigm appears insufficient to explain binding events with the biological target. Furthermore, this protein reveals the complexity of comparing or combining biochemical assay data obtained under different conditions. This can be particularly critical for artificial intelligence (AI) approaches in drug discovery that rely on large datasets of compounds activity, compiled from different sources using non-identical testing procedures. In such cases, mismatched results will compromise the model quality and its predictive power.

Lessons Learnt from COVID-19: Computational Strategies for Facing Present and Future Pandemics

Since its outbreak in December 2019, the COVID-19 pandemic has caused the death of more than 6.5 million people around the world. The high transmissibility of its causative agent, the SARS-CoV-2 virus, coupled with its potentially lethal outcome, provoked a profound global economic and social crisis. The urgency of finding suitable pharmacological tools to tame the pandemic shed light on the ever-increasing importance of computer simulations in rationalizing and speeding up the design of new drugs, further stressing the need for developing quick and reliable methods to identify novel active molecules and characterize their mechanism of action. In the present work, we aim at providing the reader with a general overview of the COVID-19 pandemic, discussing the hallmarks in its management, from the initial attempts at drug repurposing to the commercialization of Paxlovid, the first orally available COVID-19 drug. Furthermore, we analyze and discuss the role of computer-aided drug discovery (CADD) techniques, especially those that fall in the structure-based drug design (SBDD) category, in facing present and future pandemics, by showcasing several successful examples of drug discovery campaigns where commonly used methods such as docking and molecular dynamics have been employed in the rational design of effective therapeutic entities against COVID-19.

A new inactive conformation of SARS-CoV-2 main protease

The SARS‐CoV‐2 main protease (M^pro) has a pivotal role in mediating viral genome replication and transcription of the coronavirus, making it a promising target for drugs against the COVID‐19 pandemic. Here, a crystal structure is presented in which M^pro adopts an inactive state that has never been observed before, called new‐inactive. It is shown that the oxyanion loop, which is involved in substrate recognition and enzymatic activity, adopts a new catalytically incompetent conformation and that many of the key interactions of the active conformation of the enzyme around the active site are lost. Solvation/desolvation energetic contributions play an important role in the transition from the inactive to the active state, with Phe140 moving from an exposed to a buried environment and Asn142 moving from a buried environment to an exposed environment. In new‐inactive M^pro a new cavity is present near the S2′ subsite, and the N‐terminal and C‐terminal tails, as well as the dimeric interface, are perturbed, with partial destabilization of the dimeric assembly. This novel conformation is relevant both for comprehension of the mechanism of action of M^pro within the catalytic cycle and for the successful structure‐based drug design of antiviral drugs.

Binary-QSAR guided virtual screening of FDA approved drugs and compounds in clinical investigation against SARS-CoV-2 main protease

With the emergence of the new SARS-CoV-2 virus, drug repurposing studies have gained substantial importance. Combined with the efficacy of recent improvements in ligand- and target-based virtual screening approaches, virtual screening has become faster and more productive than ever. In the current study, an FDA library of approved drugs and compounds under clinical investigation were screened for their antiviral activity using the antiviral therapeutic activity binary QSAR model of the MetaCore/MetaDrug platform. Among 6733-compound collection, we found 370 compounds with a normalized therapeutic activity value greater than a cutoff of 0.75. Only these selected compounds were used for molecular docking studies against the SARS-CoV-2 main protease (M^pro). After initial short (10 ns) molecular dynamics (MD) simulations with the top-50 docking scored compounds and following molecular mechanics generalized born surface area (MM/GBSA) calculations, top-10 compounds were subjected to longer (100 ns) MD simulations and end-point MM/GBSA estimations. Our virtual screening protocol yielded Cefuroxime pivoxetil, an ester prodrug of second-generation cephalosporin antibiotic Cefuroxime, as being a considerable molecule for drug repurposing against the SARS-CoV-2 M^pro.

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.

Scaffold Hopping of α-Rubromycin Enables Direct Access to FDA-Approved Cromoglicic Acid as a SARS-CoV-2 MPro Inhibitor

The COVID-19 pandemic is still active around the globe despite the newly introduced vaccines. Hence, finding effective medications or repurposing available ones could offer great help during this serious situation. During our anti-COVID-19 investigation of microbial natural products (MNPs), we came across α-rubromycin, an antibiotic derived from Streptomyces collinus ATCC19743, which was able to suppress the catalytic activity (IC50 = 5.4 µM and Ki = 3.22 µM) of one of the viral key enzymes (i.e., MPro). However, it showed high cytotoxicity toward normal human fibroblasts (CC50 = 16.7 µM). To reduce the cytotoxicity of this microbial metabolite, we utilized a number of in silico tools (ensemble docking, molecular dynamics simulation, binding free energy calculation) to propose a novel scaffold having the main pharmacophoric features to inhibit MPro with better drug-like properties and reduced/minimal toxicity. Nevertheless, reaching this novel scaffold synthetically is a time-consuming process, particularly at this critical time. Instead, this scaffold was used as a template to explore similar molecules among the FDA-approved medications that share its main pharmacophoric features with the aid of pharmacophore-based virtual screening software. As a result, cromoglicic acid (aka cromolyn) was found to be the best hit, which, upon in vitro MPro testing, was 4.5 times more potent (IC50 = 1.1 µM and Ki = 0.68 µM) than α-rubromycin, with minimal cytotoxicity toward normal human fibroblasts (CC50 > 100 µM). This report highlights the potential of MNPs in providing unprecedented scaffolds with a wide range of therapeutic efficacy. It also revealed the importance of cheminformatics tools in speeding up the drug discovery process, which is extremely important in such a critical situation.

A higher flexibility at the SARS-CoV-2 main protease active site compared to SARS-CoV and its potentialities for new inhibitor virtual screening targeting multi-conformers

The main-protease (Mpro) catalyzes a crucial step for the SARS-CoV-2 life cycle. The recent SARS-CoV-2 presents the main protease (MCoV2pro) with 12 mutations compared to SARS-CoV (MCoV1pro). Recent studies point out that these subtle differences lead to mobility variances at the active site loops with functional implications. We use metadynamics simulations and a sort of computational analysis to probe the dynamic, pharmacophoric and catalytic environment differences between the monomers of both enzymes. So, we verify how much intrinsic distinctions are preserved in the functional dimer of MCoV2pro, as well as its implications for ligand accessibility and optimized drug screening. We find a significantly higher accessibility to open binding conformers in the MCoV2pro monomer compared to MCoV1pro. A higher hydration propensity for the MCoV2pro S2 loop with the A46S substitution seems to exercise a key role. Quantum calculations suggest that the wider conformations for MCoV2pro are less catalytically active in the monomer. However, the statistics for contacts involving the N-finger suggest higher maintenance of this activity at the dimer. Docking analyses suggest that the ability to vary the active site width can be important to improve the access of the ligand to the active site in different ways. So, we carry out a multiconformational virtual screening with different ligand bases. The results point to the importance of taking into account the protein conformational multiplicity for new promissors anti MCoV2pro ligands. We hope these results will be useful in prospecting, repurposing and/or designing new anti SARS-CoV-2 drugs.Communicated by Ramaswamy H. Sarma.

Olive-Derived Triterpenes Suppress SARS COV-2 Main Protease: A Promising Scaffold for Future Therapeutics

SARS CoV-2 pandemic is still considered a global health disaster, and newly emerged variants keep growing. A number of promising vaccines have been recently developed as a protective measure; however, cost-effective treatments are also of great importance to support this critical situation. Previously, betulinic acid has shown promising antiviral activity against SARS CoV via targeting its main protease. Herein, we investigated the inhibitory potential of this compound together with three other triterpene congeners (i.e., ursolic acid, maslinic acid, and betulin) derived from olive leaves against the viral main protease (Mpro) of the currently widespread SARS CoV-2. Interestingly, betulinic, ursolic, and maslinic acids showed significant inhibitory activity (IC50 = 3.22-14.55 µM), while betulin was far less active (IC50 = 89.67 µM). A comprehensive in-silico analysis (i.e., ensemble docking, molecular dynamic simulation, and binding-free energy calculation) was then performed to describe the binding mode of these compounds with the enzyme catalytic active site and determine the main essential structural features required for their inhibitory activity. Results presented in this communication indicated that this class of compounds could be considered as a promising lead scaffold for developing cost-effective anti-SARS CoV-2 therapeutics.

Asymmetric dynamics of dimeric SARS-CoV-2 and SARS-CoV main proteases in an apo form: Molecular dynamics study on fluctuations of active site, catalytic dyad, and hydration water

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has been widely spread around the world. It is necessary to examine the viral proteins that play a notorious role in the invasion of our body. The main protease (3CL^pro) facilitates the maturation of the coronavirus. It is thought that the dimerization of 3CL^pro leads to its catalytic activity; the detailed mechanism has, however, not been suggested. Furthermore, the structural differences between the predecessor SARS-CoV 3CL^pro and SARS-CoV-2 3CL^pro have not been fully understood. Here, we show the structural and dynamical differences between the two main proteases, and demonstrate the relationship between the dimerization and the activity via atomistic molecular dynamics simulations. Simulating monomeric and dimeric 3CL^pro systems for each protease, we show that (i) global dynamics between the two different proteases are not conserved, (ii) the dimerization stabilizes the catalytic dyad and hydration water molecules behind the dyad, and (iii) the substrate-binding site (active site) and hydration water molecules in each protomer fluctuate asymmetrically. We then speculate the roles of hydration water molecules in their catalytic activity.