Crystallographic structure of wild-type SARS-CoV-2 main protease acyl-enzyme intermediate with physiological C-terminal autoprocessing site

Proteolytic cleavage and inactivation of the TRMT1 tRNA modification enzyme by SARS-CoV-2 main protease

Nonstructural protein 5 (Nsp5) is the main protease of SARS-CoV-2 that cleaves viral polyproteins into individual polypeptides necessary for viral replication. Here, we show that Nsp5 binds and cleaves human tRNA methyltransferase 1 (TRMT1), a host enzyme required for a prevalent post-transcriptional modification in tRNAs. Human cells infected with SARS-CoV-2 exhibit a decrease in TRMT1 protein levels and TRMT1-catalyzed tRNA modifications, consistent with TRMT1 cleavage and inactivation by Nsp5. Nsp5 cleaves TRMT1 at a specific position that matches the consensus sequence of SARS-CoV-2 polyprotein cleavage sites, and a single mutation within the sequence inhibits Nsp5-dependent proteolysis of TRMT1. The TRMT1 cleavage fragments exhibit altered RNA binding activity and are unable to rescue tRNA modification in TRMT1-deficient human cells. Compared to wild-type human cells, TRMT1-deficient human cells infected with SARS-CoV-2 exhibit reduced levels of intracellular viral RNA. These findings provide evidence that Nsp5-dependent cleavage of TRMT1 and perturbation of tRNA modification patterns contribute to the cellular pathogenesis of SARS-CoV-2 infection.

Glycyrrhizic acid conjugates with amino acid methyl esters target the main protease, exhibiting antiviral activity against wild-type and nirmatrelvir-resistant SARS-CoV-2 variants

COVID-19 pandemic is predominantly caused by SARS-CoV-2, with its main protease, Mpro, playing a pivotal role in viral replication and serving as a potential target for inhibiting different variants. In this study, potent Mpro inhibitors were identified from glycyrrhizic acid (GL) derivatives with amino acid methyl/ethyl esters. Out of the 17 derivatives semisynthesized, Compounds 2, 6, 9, and 15, with methionine methyl esters, D-tyrosine methyl esters, glutamic acid methyl esters, and methionines in the carbohydrate moiety, respectively, significantly inhibited wild-type SARS-CoV-2 Mpro-mediated proteolysis, with IC50 values ranging from 0.06 μM to 0.84 μM. They also demonstrated efficacy in inhibiting trans-cleavage by mutant Mpro variants (Mpro_P132H, Mpro_E166V, Mpro_P168A, Mpro_Q189I), with IC50 values ranging from 0.05 to 0.92 μM, surpassing nirmatrelvir (IC50: 1.17∼152.9 μM). Molecular modeling revealed stronger interactions with Valine166 in the structural complex of Mpro_E166V with the compounds compared to nirmatrelvir. Moreover, these compounds efficiently inhibited the post-entry viral processes of wild-type SARS-CoV-2 single-round infectious particles (SRIPs), mitigating viral cytopathic effects and reducing replicon-driven GFP reporter signals, as well as in vitro infectivity of wild-type, Mpro_E166V, and Mpro_Q189I SRIPs, with EC50 values ranging from 0.02 to 0.53 μM. However, nirmatrelvir showed a significant decrease in inhibiting the replication of mutant SARS-CoV-2 SRIPs carrying Mpro_E166V (EC50: > 20 μM) and Mpro_Q189I (EC50: 13.2 μM) compared to wild-type SRIPs (EC50: 0.06 μM). Overall, this study identifies four GL derivatives as promising lead compounds for developing treatments against various SARS-CoV-2 strains, including Omicron, and nirmatrelvir-resistant variants.

SARS-CoV-2 3CLPro Dihedral Angles Reveal Allosteric Signaling

In allosteric proteins, identifying the pathways that signals take from allosteric ligand-binding sites to enzyme active sites or binding pockets and interfaces remains challenging. This avenue of research is motivated by the goals of understanding particular macromolecular systems of interest and creating general methods for their study. An especially important protein that is the subject of many investigations in allostery is the SARS-CoV-2 main protease (Mpro), which is necessary for coronaviral replication. It is both an attractive drug target and, due to intense interest in it for the development of pharmaceutical compounds, a gauge of the state-of-the-art approaches in studying protein inhibition. Here we develop a computational method for characterizing protein allostery and use it to study Mpro. We propose a role of the protein's C-terminal tail in allosteric modulation and warn of unintuitive traps that can plague studies of the role of protein dihedrals angles in transmitting allosteric signals.

Macrocyclic Azapeptide Nitriles: Structure-Based Discovery of Potent SARS-CoV-2 Main Protease Inhibitors as Antiviral Drugs

Given the crucial role of the main protease (Mpro) in the replication cycle of SARS-CoV-2, this viral cysteine protease constitutes a high-profile drug target. We investigated peptidomimetic azapeptide nitriles as auspicious, irreversibly acting inhibitors of Mpro. Our systematic approach combined an Mpro active-site scanning by combinatorially assembled azanitriles with structure-based design. Encouraged by the bioactive conformation of open-chain inhibitors, we conceptualized the novel chemotype of macrocyclic azanitriles whose binding mode was elucidated by cocrystallization. This strategy provided a favorable entropic contribution to target binding and resulted in the development of the extraordinarily potent Mpro inhibitor 84 with an IC50 value of 3.23 nM and a second-order rate constant of inactivation, kinac/Ki, of 448,000 M-1s-1. The open-chain Mpro inhibitor 58, along with the macrocyclic compounds 83 and 84, a broad-spectrum anticoronaviral agent, demonstrated the highest antiviral activity with EC50 values in the single-digit micromolar range. Our findings are expected to promote the future development of peptidomimetic Mpro inhibitors as anti-SARS-CoV-2 agents.

SARS-CoV-2 NSP5 antagonizes MHC II expression by subverting histone deacetylase 2

SARS-CoV-2 interferes with antigen presentation by downregulating major histocompatibility complex (MHC) II on antigen-presenting cells, but the mechanism mediating this process is unelucidated. Herein, analysis of protein and gene expression in human antigen-presenting cells reveals that MHC II is downregulated by the SARS-CoV-2 main protease, NSP5. This suppression of MHC II expression occurs via decreased expression of the MHC II regulatory protein CIITA. CIITA downregulation is independent of the proteolytic activity of NSP5, and rather, NSP5 delivers HDAC2 to the transcription factor IRF3 at an IRF-binding site within the CIITA promoter. Here, HDAC2 deacetylates and inactivates the CIITA promoter. This loss of CIITA expression prevents further expression of MHC II, with this suppression alleviated by ectopic expression of CIITA or knockdown of HDAC2. These results identify a mechanism by which SARS-CoV-2 limits MHC II expression, thereby delaying or weakening the subsequent adaptive immune response.

De novo design of potential peptide analogs against the main protease of Omicron variant using in silico studies

SARS-CoV-2 and its variants are crossing the immunity barrier induced through vaccination. Recent Omicron sub-variants are highly transmissible and have a low mortality rate. Despite the low severity of Omicron variants, these new variants are known to cause acute post-infectious syndromes. Nowadays, novel strategies to develop new potential inhibitors for SARS-CoV-2 and other Omicron variants have gained prominence. For viral replication and survival the main protease of SARS-CoV-2 plays a vital role. Peptide-like inhibitors that mimic the substrate peptide have already proved to be effective in inhibiting the Mpro of SARS-CoV-2 variants. Our systematic canonical amino acid point mutation analysis on the native peptide has revealed various ways to improve the native peptide of the main protease. Multi mutation analysis has led us to identify and design potent peptide-analog inhibitors that act against the Mpro of the Omicron sub-variants. Our in-depth analysis of all-atom molecular dynamics studies has paved the way to characterize the atomistic behavior of Mpro in Omicron variants. Our goal is to develop potent peptide-analogs that could be therapeutically effective against Omicron and its sub-variants.

Statine-based peptidomimetic compounds as inhibitors for SARS-CoV-2 main protease (SARS-CoV‑2 Mpro)

COVID-19 is a multisystemic disease caused by the SARS-CoV-2 airborne virus, a member of the Coronaviridae family. It has a positive sense single-stranded RNA genome and encodes two non-structural proteins through viral cysteine-proteases processing. Blocking this step is crucial to control virus replication. In this work, we reported the synthesis of 23 statine-based peptidomimetics to determine their ability to inhibit the main protease (Mpro) activity of SARS-CoV-2. Among the 23 peptidomimetics, 15 compounds effectively inhibited Mpro activity by 50% or more, while three compounds (7d, 8e, and 9g) exhibited maximum inhibition above 70% and IC50 < 1 µM. Compounds 7d, 8e, and 9g inhibited roughly 80% of SARS-CoV-2 replication and proved no cytotoxicity. Molecular docking simulations show putative hydrogen bond and hydrophobic interactions between specific amino acids and these inhibitors. Molecular dynamics simulations further confirmed the stability and persisting interactions in Mpro's subsites, exhibiting favorable free energy binding (ΔGbind) values. These findings suggest the statine-based peptidomimetics as potential therapeutic agents against SARS-CoV-2 by targeting Mpro.

In silico study of alkaloids with quercetin nucleus for inhibition of SARS-CoV-2 protease and receptor cell protease

Covid-19 disease caused by the deadly SARS-CoV-2 virus is a serious and threatening global health issue declared by the WHO as an epidemic. Researchers are studying the design and discovery of drugs to inhibit the SARS-CoV-2 virus due to its high mortality rate. The main Covid-19 virus protease (Mpro) and human transmembrane protease, serine 2 (TMPRSS2) are attractive targets for the study of antiviral drugs against SARS-2 coronavirus. Increasing consumption of herbal medicines in the community and a serious approach to these drugs have increased the demand for effective herbal substances. Alkaloids are one of the most important active ingredients in medicinal plants that have wide applications in the pharmaceutical industry. In this study, seven alkaloid ligands with Quercetin nucleus for the inhibition of Mpro and TMPRSS2 were studied using computational drug design including molecular docking and molecular dynamics simulation (MD). Auto Dock software was used to evaluate molecular binding energy. Three ligands with the most negative docking score were selected to be entered into the MD simulation procedure. To evaluate the protein conformational changes induced by tested ligands and calculate the binding energy between the ligands and target proteins, GROMACS software based on AMBER03 force field was used. The MD results showed that Phyllospadine and Dracocephin-A form stable complexes with Mpro and TMPRSS2. Prolinalin-A indicated an acceptable inhibitory effect on Mpro, whereas it resulted in some structural instability of TMPRSS2. The total binding energies between three ligands, Prolinalin-A, Phyllospadine and Dracocephin-A and two proteins MPro and TMRPSS2 are (-111.235 ± 15.877, - 75.422 ± 11.140), (-107.033 ± 9.072, -84.939 ± 10.155) and (-102.941 ± 9.477, - 92.451 ± 10.539), respectively. Since the binding energies are at a minimum, this indicates confirmation of the proper binding of the ligands to the proteins. Regardless of some Prolinalin-A-induced TMPRSS2 conformational changes, it may properly bind to TMPRSS2 binding site due to its acceptable binding energy. Therefore, these three ligands can be promising candidates for the development of drugs to treat infections caused by the SARS-CoV-2 virus.

Protocol to identify flavonoid antagonists of the SARS-CoV-2 main protease

Flavonoids are naturally occurring metabolites of plants that can inhibit the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) main protease (Mpro), which is required for viral replication. Here, we present a protocol to identify flavonoid antagonists of the SARS-CoV-2 Mpro. We describe steps for the expression and purification of Mpro and a kinetic enzymatic assay for Mpro activity using a dequenching fluorescence resonance energy transfer peptide substrate. We then detail procedures for using this enzymatic assay to test flavonoid antagonism and reversible inhibition. For complete details on the use and execution of this protocol, please refer to Lin et al.1.

Distal Protein-Protein Interactions Contribute to SARS-CoV-2 Main Protease Substrate Binding and Nirmatrelvir Resistance

SARS-CoV-2 main protease, M pro , is responsible for the processing of the viral polyproteins into individual proteins, including the protease itself. M pro is a key target of anti-COVID-19 therapeutics such as nirmatrelvir (the active component of Paxlovid). Resistance mutants identified clinically and in viral passage assays contain a combination of active site mutations (e.g. E166V, E166A, L167F), which reduce inhibitor binding and enzymatic activity, and non-active site mutations (e.g. P252L, T21I, L50F), which restore the fitness of viral replication. Although the mechanism of resistance for the active site mutations is apparent, the role of the non-active site mutations in fitness rescue remains elusive. In this study, we use the model system of a M pro triple mutant (L50F/E166A/L167F) that confers not only nirmatrelvir drug resistance but also a similar fitness of replication compared to the wild-type both in vitro and in vivo. By comparing peptide and full-length M pro protein as substrates, we demonstrate that the binding of M pro substrate involves more than residues in the active site. In particular, L50F and other non-active site mutations can enhance the M pro dimer-dimer interactions and help place the nsp5-6 substrate at the enzyme catalytic center. The structural and enzymatic activity data of M pro L50F, L50F/E166A/L167F, and others underscore the importance of considering the whole substrate protein in studying M pro and substrate interactions, and offers important insights into M pro function, resistance development, and inhibitor design.

Consensus Pharmacophore Strategy For Identifying Novel SARS-Cov-2 Mpro Inhibitors from Large Chemical Libraries

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) main Protease (Mpro) is an enzyme that cleaves viral polyproteins translated from the viral genome and is critical for viral replication. Mpro is a target for anti-SARS-CoV-2 drug development, and multiple Mpro crystals complexed with competitive inhibitors have been reported. In this study, we aimed to develop an Mpro consensus pharmacophore as a tool to expand the search for inhibitors. We generated a consensus model by aligning and summarizing pharmacophoric points from 152 bioactive conformers of SARS-CoV-2 Mpro inhibitors. Validation against a library of conformers from a subset of ligands showed that our model retrieved poses that reproduced the crystal-binding mode in 77% of the cases. Using models derived from a consensus pharmacophore, we screened >340 million compounds. Pharmacophore-matching and chemoinformatics analyses identified new potential Mpro inhibitors. The candidate compounds were chemically dissimilar to the reference set, and among them, demonstrating the relevance of our model. We evaluated the effect of 16 candidates on Mpro enzymatic activity finding that seven have inhibitory activity. Three compounds (1, 4, and 5) had IC50 values in the midmicromolar range. The Mpro consensus pharmacophore reported herein can be used to identify compounds with improved activity and novel chemical scaffolds against Mpro. The method developed for its generation is provided as an open-access code (https://github.com/AngelRuizMoreno/ConcensusPharmacophore) and can be applied to other pharmacological targets.

On Inactivation of the Coronavirus Main Protease

A deeper understanding of the inactive conformations of the coronavirus main protease (MPro) could inform the design of allosteric drugs. Based on extensive molecular dynamics simulations, we built a Markov State Model to investigate structural changes that can inactivate the SARS-CoV-2 MPro. In a subset of structures, one subunit of the homodimer assumes an inactive conformation that resembles an inactive crystal structure. However, contradicting the widely held half-of-sites activity hypothesis, the most populated enzyme structures have two active subunits. We then used transition path theory (TPT) and the Jensen-Shannon Divergence (JSD) to pinpoint residues involved in the inactivation process. A π stack between Phe140 and His163 is a key feature that can distinguish active and inactive conformations of MPro. Each subunit has unique inactive conformations stabilized by π stacking interactions involving residues Phe140, Tyr118, His163, and His172, a hydrogen bonding network centered around His163 and His172, and a modified network of interactions in the dimer interface. The importance of these residues in maintaining an active structure explains the sensitivity of enzymatic activity to site-directed mutagenesis.

Exploration of the P1 residue in 3CL protease inhibitors leading to the discovery of a 2-tetrahydrofuran P1 replacement

The virally encoded 3C-like protease (3CLpro) is a well-validated drug target for the inhibition of coronaviruses including Severe Acute Respiratory Syndrome Coronavirus-2 (SARS-CoV-2). Most inhibitors of 3CLpro are peptidomimetic, with a γ-lactam in place of Gln at the P1 position of the pseudopeptide chain. An effort was pursued to identify a viable alternative to the γ-lactam P1 mimetic which would improve physicochemical properties while retaining affinity for the target. Discovery of a 2-tetrahydrofuran as a suitable P1 replacement that is a potent enzymatic inhibitor of 3CLpro in SARS-CoV-2 virus is described herein.

Unexpected Single-Ligand Occupancy and Negative Cooperativity in the SARS-CoV-2 Main Protease

Many homodimeric enzymes tune their functions by exploiting either negative or positive cooperativity between subunits. In the SARS-CoV-2 Main protease (Mpro) homodimer, the latter has been suggested by symmetry in most of the 500 reported protease/ligand complex structures solved by macromolecular crystallography (MX). Here we apply the latter to both covalent and noncovalent ligands in complex with Mpro. Strikingly, our experiments show that the occupation of both active sites of the dimer originates from an excess of ligands. Indeed, cocrystals obtained using a 1:1 ligand/protomer stoichiometry lead to single occupation only. The empty binding site exhibits a catalytically inactive geometry in solution, as suggested by molecular dynamics simulations. Thus, Mpro operates through negative cooperativity with the asymmetric activity of the catalytic sites. This allows it to function with a wide range of substrate concentrations, making it resistant to saturation and potentially difficult to shut down, all properties advantageous for the virus' adaptability and resistance.

Proteolytic cleavage and inactivation of the TRMT1 tRNA modification enzyme by SARS-CoV-2 main protease

Nonstructural protein 5 (Nsp5) is the main protease of SARS-CoV-2 that cleaves viral polyproteins into individual polypeptides necessary for viral replication. Here, we show that Nsp5 binds and cleaves human tRNA methyltransferase 1 (TRMT1), a host enzyme required for a prevalent post-transcriptional modification in tRNAs. Human cells infected with SARS-CoV-2 exhibit a decrease in TRMT1 protein levels and TRMT1-catalyzed tRNA modifications, consistent with TRMT1 cleavage and inactivation by Nsp5. Nsp5 cleaves TRMT1 at a specific position that matches the consensus sequence of SARS-CoV-2 polyprotein cleavage sites, and a single mutation within the sequence inhibits Nsp5-dependent proteolysis of TRMT1. The TRMT1 cleavage fragments exhibit altered RNA binding activity and are unable to rescue tRNA modification in TRMT1-deficient human cells. Compared to wildtype human cells, TRMT1-deficient human cells infected with SARS-CoV-2 exhibit reduced levels of intracellular viral RNA. These findings provide evidence that Nsp5-dependent cleavage of TRMT1 and perturbation of tRNA modification patterns contribute to the cellular pathogenesis of SARS-CoV-2 infection.

Multiple redox switches of the SARS-CoV-2 main protease in vitro provide opportunities for drug design

Besides vaccines, the development of antiviral drugs targeting SARS-CoV-2 is critical for preventing future COVID outbreaks. The SARS-CoV-2 main protease (Mpro), a cysteine protease with essential functions in viral replication, has been validated as an effective drug target. Here, we show that Mpro is subject to redox regulation in vitro and reversibly switches between the enzymatically active dimer and the functionally dormant monomer through redox modifications of cysteine residues. These include a disulfide-dithiol switch between the catalytic cysteine C145 and cysteine C117, and generation of an allosteric cysteine-lysine-cysteine SONOS bridge that is required for structural stability under oxidative stress conditions, such as those exerted by the innate immune system. We identify homo- and heterobifunctional reagents that mimic the redox switching and inhibit Mpro activity. The discovered redox switches are conserved in main proteases from other coronaviruses, e.g. MERS-CoV and SARS-CoV, indicating their potential as common druggable sites.

Computational Estimation of Residues Involving Resistance to the SARS-CoV-2 Main Protease Inhibitor Ensitrelvir Based on Virtual Alanine Scan of the Active Site

Ensitrelvir is a noncovalent inhibitor of the main protease (Mpro) of severe acute respiratory syndrome coronavirus 2. Acquisition of drug resistance in virus-derived proteins is a serious therapeutic concern, and drug resistance occurs due to amino acid mutations. In this study, we computationally constructed 24 mutants, in which one residue around the active site was replaced with alanine and performed molecular dynamics simulations to the complex of Mpro and ensitrelvir to predict the residues involved in drug resistance. We evaluated the changes in the entire protein structure and ligand configuration in each of these mutants and estimated which residues were involved in ensitrelvir recognition. This method is called a virtual alanine scan. In nine mutants (S1A, T26A, H41A, M49A, L141A, H163A, E166A, V186A, and R188A), although the entire protein structure and catalytic dyad (cysteine (Cys)145 and histidine (His)41) were not significantly moved, the ensitrelvir configuration changed. Thus, it is considered that these mutants did not recognize ensitrelvir while maintaining Mpro enzymatic activities, and Ser1, Thr26, His41, Met49, Leu141, His163, Glu166, Val186, and Arg188 may be related to ensitrelvir resistance. The ligand shift noted in M49A was similar to that observed in M49I, which has been shown to be experimentally ensitrelvir resistant. These findings suggest that our research approach can predict mutations that incite drug resistance.

Inhibitory effect of thymoquinone from Nigella sativa against SARS-CoV-2 main protease. An in-silico study

Abstract Nigella sativa is known for the safety profile, containing a wealth of useful antiviral compounds. The main protease (Mpro, 3CLpro) of severe acute respiratory syndrome 2 (SARS-CoV-2) is being considered as one of the most attractive viral target, processing the polyproteins during viral pathogenesis and replication. In the current investigation we analyzed the potency of active component, thymoquinone (TQ) of Nigella sativa against SARS-CoV-2 Mpro. The structures of TQ and Mpro was retrieved from PubChem (CID10281) and Protein Data Bank (PDB ID 6MO3) respectively. The Mpro and TQ were docked and the complex was subjected to molecular dynamic (MD) simulations for a period 50ns. Protein folding effect was analyzed using radius of gyration (Rg) while stability and flexibility was measured, using root means square deviations (RMSD) and root means square fluctuation (RMSF) respectively. The simulation results shows that TQ is exhibiting good binding activity against SARS-CoV-2 Mpro, interacting many residues, present in the active site (His41, Cys145) and also the Glu166, facilitating the pocket shape. Further, experimental approaches are needed to validate the role of TQ against virus infection. The TQ is interfering with pocket maintaining residues as well as active site of virus Mpro which may be used as a potential inhibitor against SARS-CoV-2 for better management of COVID-19.

Lessons learnt from the preclinical discovery and development of ensitrelvir as a COVID-19 therapeutic option

INTRODUCTION

The COVID-19 pandemic stimulated the development of several therapeutic tools with several degrees of success. Ensitrelvir, a protease inhibitor that blocks the replication of SARS-CoV-2, can reduce the viral load and the severity of symptoms in infected patients and become available for emergency use in Japan. Clinical trials showed a good tolerability profile although the potential for interactions with substrates, inhibitors, and inducers of CYP3A must be considered. The occurrence of resistance is also a matter of investigation.

AREAS COVERED

In this article, the authors describe the development of ensitrelvir starting from the identification of the molecule to the pre-clinical and clinical trials up to the post-authorization phase.

EXPERT OPINION

Ensitrelvir was developed in a late phase of the pandemic when the availability of patients that can be candidate to enter the clinical trial was limited with consequences for the possibility of assessing certain outcomes and for the robustness of results. Although the evidence about the benefits of ensitrelvir in COVID-19 is not questionable, the problems of interactions with other drugs, emerging resistant variants, the availability of alternative therapeutic options, costs, and accessibility will concur to its probable limited clinical use in the future.

Assessing Genomic Mutations in SARS-CoV-2: Potential Resistance to Antiviral Drugs in Viral Populations from Untreated COVID-19 Patients

Naturally occurring SARS-CoV-2 variants mutated in genomic regions targeted by antiviral drugs have not been extensively studied. This study investigated the potential of the RNA-dependent RNA polymerase (RdRp) complex subunits and non-structural protein (Nsp)5 of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) to accumulate natural mutations that could affect the efficacy of antiviral drugs. To this aim, SARS-CoV-2 genomic sequences isolated from 4155 drug-naive individuals from southern Italy were analyzed using the Illumina MiSeq platform. Sequencing of the 4155 samples showed the following viral variant distribution: 71.2% Delta, 22.2% Omicron, and 6.4% Alpha. In the Nsp12 sequences, we found 84 amino acid substitutions. The most common one was P323L, detected in 3777/4155 (91%) samples, with 2906/3777 (69.9%) also showing the G671S substitution in combination. Additionally, we identified 28, 14, and 24 different amino acid substitutions in the Nsp5, Nsp7, and Nsp8 genomic regions, respectively. Of note, the V186F and A191V substitutions, affecting residues adjacent to the active site of Nsp5 (the target of the antiviral drug Paxlovid), were found in 157/4155 (3.8%) and 3/4155 (0.07%) samples, respectively. In conclusion, the RdRp complex subunits and the Nsp5 genomic region exhibit susceptibility to accumulating natural mutations. This susceptibility poses a potential risk to the efficacy of antiviral drugs, as these mutations may compromise the drug ability to inhibit viral replication.

A novel class of broad-spectrum active-site-directed 3C-like protease inhibitors with nanomolar antiviral activity against highly immune-evasive SARS-CoV-2 Omicron subvariants

Antivirals with broad coronavirus activity are important for treating high-risk individuals exposed to the constantly evolving SARS-CoV-2 variants of concern (VOCs) as well as emerging drug-resistant variants. We developed and characterized a novel class of active-site-directed 3-chymotrypsin-like protease (3CLpro) inhibitors (C2-C5a). Our lead direct-acting antiviral (DAA), C5a, is a non-covalent, non-peptide with a dissociation constant of 170 nM against recombinant SARS-CoV-2 3CLpro. The compounds C2-C5a exhibit broad-spectrum activity against Omicron subvariants (BA.5, BQ.1.1, and XBB.1.5) and seasonal human coronavirus-229E infection in human cells. Notably, C5a has median effective concentrations of 30-50 nM against BQ.1.1 and XBB.1.5 in two different human cell lines. X-ray crystallography has confirmed the unique binding modes of C2-C5a to the 3CLpro, which can limit virus cross-resistance to emerging Paxlovid-resistant variants. We tested the effect of C5a with two of our newly discovered host-directed antivirals (HDAs): N-0385, a TMPRSS2 inhibitor, and bafilomycin D (BafD), a human vacuolar H+-ATPase [V-ATPase] inhibitor. We demonstrated a synergistic action of C5a in combination with N-0385 and BafD against Omicron BA.5 infection in human Calu-3 lung cells. Our findings underscore that a SARS-CoV-2 multi-targeted treatment for circulating Omicron subvariants based on DAAs (C5a) and HDAs (N-0385 or BafD) can lead to therapeutic benefits by enhancing treatment efficacy. Furthermore, the high-resolution structures of SARS-CoV-2 3CLpro in complex with C2-C5a will facilitate future rational optimization of our novel broad-spectrum active-site-directed 3C-like protease inhibitors.

Discovery of antiviral SARS-CoV-2 main protease inhibitors by structure-guided hit-to-lead optimization of carmofur

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) main protease (Mpro) has been targeted for the development of anti-SARS-CoV-2 agents against COVID-19 infection because Mpro processes essential viral polyproteins and plays a key role in SARS-CoV-2 replication. In this study, we report the development of novel SARS-CoV-2 Mpro inhibitors derived from carmofur, a previously identified compound that has shown moderate potency as a covalent inhibitor of SARS-CoV-2 Mpro. To employ a structure-guided drug design strategy, a putative intact binding mode of carmofur at catalytic active site of Mpro was initially predicted by docking simulation. Based on the predicted binding mode, a series of carmofur derivatives aiming to occupy the Mpro substrate binding regions were investigated for structure-activity relationship analysis. As a result, an indole-based derivative, speculated to interact with the S4 binding pocket, 21b (IC50 = 1.5 ± 0.1 μM) was discovered. Its structure was further modified and evaluated in silico by combining docking simulation, free energy perturbation calculation and subpocket interaction analysis to optimize the interactions at the S2 and S4 binding pockets. Among the newly designed novel derivatives, 21h and 21i showed the best inhibitory potencies against Mpro with IC50 values of 0.35 and 0.37 μM, respectively. Moreover, their antiviral activities were confirmed with EC50 values of 20-30 μM in the SARS-CoV-2-infected cell-based assay, suggesting that these novel Mpro inhibitors could be applied as potential lead compounds for the development of substantial anti-SARS-CoV-2 agents.

Insights into the mechanism of SARS-CoV-2 main protease autocatalytic maturation from model precursors

A critical step for SARS-CoV-2 assembly and maturation involves the autoactivation of the main protease (MProWT) from precursor polyproteins. Upon expression, a model precursor of MProWT mediates its own release at its termini rapidly to yield a mature dimer. A construct with an E290A mutation within MPro exhibits time dependent autoprocessing of the accumulated precursor at the N-terminal nsp4/nsp5 site followed by the C-terminal nsp5/nsp6 cleavage. In contrast, a precursor containing E290A and R298A mutations (MProM) displays cleavage only at the nsp4/nsp5 site to yield an intermediate monomeric product, which is cleaved at the nsp5/nsp6 site only by MProWT. MProM and the catalytic domain (MPro1-199) fused to the truncated nsp4 region also show time-dependent conversion in vitro to produce MProM and MPro1-199, respectively. The reactions follow first-order kinetics indicating that the nsp4/nsp5 cleavage occurs via an intramolecular mechanism. These results support a mechanism involving an N-terminal intramolecular cleavage leading to an increase in the dimer population and followed by an intermolecular cleavage at the C-terminus. Thus, targeting the predominantly monomeric MPro precursor for inhibition may lead to the identification of potent drugs for treatment.

αβ,α'β'-Diepoxyketones are mechanism-based inhibitors of nucleophilic cysteine enzymes

Epoxides are an established class of electrophilic alkylating agents that react with nucleophilic protein residues. We report αβ,α'β'-diepoxyketones (DEKs) as a new type of mechanism-based inhibitors of nucleophilic cysteine enzymes. Studies with the L,D-transpeptidase LdtMt2 from Mycobacterium tuberculosis and the main protease from SARS-CoV-2 (Mpro) reveal that following epoxide ring opening by a nucleophilic cysteine, further reactions can occur, leading to irreversible alkylation.

Bayesian Inference Elucidates the Catalytic Competency of the SARS-CoV-2 Main Protease 3CLpro

The main protease of SARS-CoV-2, 3CLpro, is a dimeric enzyme that is indispensable to viral replication and presents an attractive opportunity for therapeutic intervention. Previous reports regarding the key properties of 3CLpro and its highly similar SARS-CoV homologue conflict dramatically. Values of the dimeric Kd and enzymic kcat/KM differ by 106- and 103-fold, respectively. Establishing a confident benchmark of the intrinsic capabilities of this enzyme is essential for combating the current pandemic as well as potential future outbreaks. Here, we use enzymatic methods to characterize the dimerization and catalytic efficiency of the authentic protease from SARS-CoV-2. Specifically, we use the rigor of Bayesian inference in a Markov Chain Monte Carlo analysis of progress curves to circumvent the limitations of traditional Michaelis-Menten initial rate analysis. We report that SARS-CoV-2 3CLpro forms a dimer at pH 7.5 that has Kd = 16 ± 4 nM and is capable of catalysis with kcat = 9.9 ± 1.5 s-1, KM = 0.23 ± 0.01 mM, and kcat/KM = (4.3 ± 0.7) × 104 M-1 s-1. We also find that enzymatic activity decreases substantially in solutions of high ionic strength, largely as a consequence of impaired dimerization. We conclude that 3CLpro is a more capable catalyst than appreciated previously, which has important implications for the design of antiviral therapeutic agents that target 3CLpro.

Structural biology of SARS-CoV-2 Mpro and drug discovery

Since its outbreak in late 2019, the COVID-19 pandemic has drawn enormous attention worldwide as a consequence of being the most disastrous infectious disease in the past century. As one of the most immediately druggable targets of SARS-CoV-2, the main protease (Mpro) has been studied thoroughly. In this review, we provide a comprehensive summary of recent advances in structural studies of Mpro, which provide new knowledge about Mpro in terms of its biological function, structural characteristics, substrate specificity, and autocleavage process. We examine the remarkable strides made in targeting Mpro for drug discovery during the pandemic. We summarize insights into the current understanding of the structural features of Mpro and the discovery of existing Mpro-targeting drugs, illuminating pathways for the future development of anti-SARS-CoV-2 therapeutics.

Plant flavonoid inhibition of SARS-CoV-2 main protease and viral replication

Plant-based flavonoids have been evaluated as inhibitors of β-coronavirus replication and as therapies for COVID-19 on the basis of their safety profile and widespread availability. The SARS-CoV-2 main protease (Mpro) has been implicated as a target for flavonoids in silico. Yet no comprehensive in vitro testing of flavonoid activity against SARS-CoV-2 Mpro has heretofore been performed. We screened 1,019 diverse flavonoids for their ability to inhibit SARS-CoV-2 Mpro. Multiple structure-activity relationships were identified among active compounds such as enrichment of galloylated flavonoids and biflavones, including multiple biflavone analogs of apigenin. In a cell-based SARS-CoV-2 replication assay, the most potent inhibitors were apigenin and the galloylated pinocembrin analog, pinocembrin 7-O-(3''-galloyl-4'',6''-(S)-hexahydroxydiphenoyl)-beta-D-glucose (PGHG). Molecular dynamic simulations predicted that PGHG occludes the S1 binding site via a galloyl group and induces a conformational change in Mpro. These studies will advance the development of plant-based flavonoids-including widely available natural products-to target β-coronaviruses.

Alkyne as a Latent Warhead to Covalently Target SARS-CoV-2 Main Protease

There is an urgent need for improved therapy to better control the ongoing COVID-19 pandemic. The main protease Mpro plays a pivotal role in SARS-CoV-2 replications, thereby representing an attractive target for antiviral development. We seek to identify novel electrophilic warheads for efficient, covalent inhibition of Mpro. By comparing the efficacy of a panel of warheads installed on a common scaffold against Mpro, we discovered that the terminal alkyne could covalently modify Mpro as a latent warhead. Our biochemical and X-ray structural analyses revealed the irreversible formation of the vinyl-sulfide linkage between the alkyne and the catalytic cysteine of Mpro. Clickable probes based on the alkyne inhibitors were developed to measure target engagement, drug residence time, and off-target effects. The best alkyne-containing inhibitors potently inhibited SARS-CoV-2 infection in cell infection models. Our findings highlight great potentials of alkyne as a latent warhead to target cystine proteases in viruses and beyond.

Identification of motif-based interactions between SARS-CoV-2 protein domains and human peptide ligands pinpoint antiviral targets

The virus life cycle depends on host-virus protein-protein interactions, which often involve a disordered protein region binding to a folded protein domain. Here, we used proteomic peptide phage display (ProP-PD) to identify peptides from the intrinsically disordered regions of the human proteome that bind to folded protein domains encoded by the SARS-CoV-2 genome. Eleven folded domains of SARS-CoV-2 proteins were found to bind 281 peptides from human proteins, and affinities of 31 interactions involving eight SARS-CoV-2 protein domains were determined (KD ∼ 7-300 μM). Key specificity residues of the peptides were established for six of the interactions. Two of the peptides, binding Nsp9 and Nsp16, respectively, inhibited viral replication. Our findings demonstrate how high-throughput peptide binding screens simultaneously identify potential host-virus interactions and peptides with antiviral properties. Furthermore, the high number of low-affinity interactions suggest that overexpression of viral proteins during infection may perturb multiple cellular pathways.

The H163A mutation unravels an oxidized conformation of the SARS-CoV-2 main protease

The main protease of SARS-CoV-2 (Mpro) is an important target for developing COVID-19 therapeutics. Recent work has highlighted Mpro's susceptibility to undergo redox-associated conformational changes in response to cellular and immune-system-induced oxidation. Despite structural evidence indicating large-scale rearrangements upon oxidation, the mechanisms of conformational change and its functional consequences are poorly understood. Here, we present the crystal structure of an Mpro point mutant (H163A) that shows an oxidized conformation with the catalytic cysteine in a disulfide bond. We hypothesize that Mpro adopts this conformation under oxidative stress to protect against over-oxidation. Our metadynamics simulations illustrate a potential mechanism by which H163 modulates this transition and suggest that this equilibrium exists in the wild type enzyme. We show that other point mutations also significantly shift the equilibrium towards this state by altering conformational free energies. Unique avenues of SARS-CoV-2 research can be explored by understanding how H163 modulates this equilibrium.

Recognition and Cleavage of Human tRNA Methyltransferase TRMT1 by the SARS-CoV-2 Main Protease

The SARS-CoV-2 main protease (Mpro) is critical for the production of functional viral proteins during infection and, like many viral proteases, can also target host proteins to subvert their cellular functions. Here, we show that the human tRNA methyltransferase TRMT1 can be recognized and cleaved by SARS-CoV-2 Mpro. TRMT1 installs the N2,N2-dimethylguanosine (m2,2G) modification on mammalian tRNAs, which promotes global protein synthesis and cellular redox homeostasis. We find that Mpro can cleave endogenous TRMT1 in human cell lysate, resulting in removal of the TRMT1 zinc finger domain required for tRNA modification activity in cells. Evolutionary analysis shows that the TRMT1 cleavage site is highly conserved in mammals, except in Muroidea, where TRMT1 may be resistant to cleavage. In primates, regions outside the cleavage site with rapid evolution could indicate adaptation to ancient viral pathogens. We determined the structure of a TRMT1 peptide in complex with Mpro, revealing a substrate binding conformation distinct from the majority of available Mpro-peptide complexes. Kinetic parameters for peptide cleavage showed that the TRMT1(526-536) sequence is cleaved with comparable efficiency to the Mpro-targeted nsp8/9 viral cleavage site. Mutagenesis studies and molecular dynamics simulations together indicate that kinetic discrimination occurs during a later step of Mpro-mediated proteolysis that follows substrate binding. Our results provide new information about the structural basis for Mpro substrate recognition and cleavage that could help inform future therapeutic design and raise the possibility that proteolysis of human TRMT1 during SARS-CoV-2 infection suppresses protein translation and oxidative stress response to impact viral pathogenesis.

Investigating theobromine as a potential anti-human coronaviral agent

Coronaviruses (CoVs) have long been known to infect humans, mainly alpha-CoV and beta-CoV. The vaccines developed for SARS-CoV-2 are likely not effective against other coronavirus species, whereas the risk of the emergence of new strains that may cause the next epidemic/pandemic is high. The development of antiviral drugs that are effective across different CoVs represents a viable strategy for improving pandemic preparedness. In this study, we aim to identify pan-coronaviral agents by targeting the conserved main protease (Mpro). For drug screening, the catalytic dyad of four human CoVs (HCoVs: SARS-CoV-2, and seasonal CoV NL63, OC43, and 229E) was targeted by molecular docking. The identified leading candidate theobromine, a xanthine derivative, was further tested in cell culture models of coronavirus infection. Theobromine binds strongly with the catalytic dyad (His41 and Cys144/145) of SARS-CoV-2 and HCoV-NL63 Mpro, mildly with HCoV-OC43, but not with HCoV-229E. However, theobromine only shows dose-dependent inhibition in Calu3 cells inoculated with SARS-CoV-2, but not in cells inoculated with seasonal CoVs. Theobromine exerts antiviral activity against coronavirus infections potentially through targeting Mpro. However, the antiviral potency is distinct among different CoVs.

Theoretical insights into the binding interaction of Nirmatrelvir with SARS-CoV-2 Mpro mutants (C145A and C145S): MD simulations and binding free-energy calculation to understand drug resistance

Mpro, the main protease and a crucial enzyme in SARS-CoV-2 is the most fascinating molecular target for pharmacological treatment and is also liable for viral protein maturation. For antiviral therapy, no drugs have been approved clinically to date. Targeting the Mpro with a compound having inhibitory properties against it can hinder viral replication. The therapeutic potential of the antiviral compound Nirmatrelvir (NMV) against SARS-CoV-2 Mpro was investigated using a systematic approach of molecular docking, MD simulations, and binding free energy calculation based on the MM-GBSA method. NMV, a covalent inhibitor with a recently revealed chemical structure, is a promising oral antiviral clinical candidate with significant in vitro anti-SARS-CoV-2 action in third-phase clinical trials. To explore the therapeutic ability and possible drug resistance, the Mpro system was studied for WT and two of its primary mutants (C145A & C145S). The protein-ligand (Mpro/NMV) complexes were further examined through long MD simulations to check the possible drug resistance in the mutants. To understand the binding affinity, the MM-GBSA method was applied to the Mpro/NMV complexes. Moreover, PCA analysis confirms the detachment of the linker region from the major domains in C145S and C145A mutants allowing for conformational alterations in the active-site region. Based on the predicted biological activities and binding affinities of NMV to WT and mutant (C145A & C145S) Mpro, it can be stipulated that NMV may have conventional potency to act as an anti-viral agent against WT Mpro, while the catalytic-dyad mutations may show substantial mutation-induced drug resistance.Communicated by Ramaswamy H. Sarma.

Design of Tetra-Peptide Ligands of Antibody Fc Regions Using In Silico Combinatorial Library Screening

Peptides, or short chains of amino-acid residues, are becoming increasingly important as active ingredients of drugs and as crucial probes and/or tools in medical, biotechnological, and pharmaceutical research. Situated at the interface between small molecules and larger macromolecular systems, they pose a difficult challenge for computational methods. We report an in silico peptide library generation and prioritization workflow using CmDock for identifying tetrapeptide ligands that bind to Fc regions of antibodies that is analogous to known in vitro recombinant peptide libraries' display and expression systems. The results of our in silico study are in accordance with existing scientific literature on in vitro peptides that bind to antibody Fc regions. In addition, we postulate an evolving in silico library design workflow that will help circumvent the combinatorial problem of in vitro comprehensive peptide libraries by focusing on peptide subunits that exhibit favorable interaction profiles in initial in silico peptide generation and testing.

DNA-encoded chemical libraries yield non-covalent and non-peptidic SARS-CoV-2 main protease inhibitors

The development of SARS-CoV-2 main protease (Mpro) inhibitors for the treatment of COVID-19 has mostly benefitted from X-ray structures and preexisting knowledge of inhibitors; however, an efficient method to generate Mpro inhibitors, which circumvents such information would be advantageous. As an alternative approach, we show here that DNA-encoded chemistry technology (DEC-Tec) can be used to discover inhibitors of Mpro. An affinity selection of a 4-billion-membered DNA-encoded chemical library (DECL) using Mpro as bait produces novel non-covalent and non-peptide-based small molecule inhibitors of Mpro with low nanomolar Ki values. Furthermore, these compounds demonstrate efficacy against mutant forms of Mpro that have shown resistance to the standard-of-care drug nirmatrelvir. Overall, this work demonstrates that DEC-Tec can efficiently generate novel and potent inhibitors without preliminary chemical or structural information.

In Vitro Reconstitution and Analysis of SARS-CoV-2/Host Protein-Protein Interactions

The emergence of viral threats such as Ebola, ZIKA, and severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) requires a rapid and efficient approach for elucidating mechanisms of pathogenesis and development of therapeutics. In this context, cell-free protein synthesis (CFPS) holds a promise to resolve the bottlenecks of multiplexed protein production and interaction analysis among host and pathogen proteins. Here, we applied a eukaryotic CFPS system based on Leishmania tarentolae extract (LTE) protein expression in combination with AlphaLISA proximity-based protein interaction technology to identify intraviral and viral-human protein interactions of SARS-CoV-2 virus that can potentially be targeted by the existing or novel antiviral therapeutics. We produced and tested 54 putative human-viral protein pairs in vitro and identified 45 direct binary protein interactions. As a casing example of the assay's suitability for drug development applications, we analyzed the effect of a putative biologic on the human angiotensin-converting enzyme 2/receptor-binding domain (hACE2/RBD) interaction. This suggests that the presented pathogen characterization platform can facilitate the development of new therapeutic agents.

In Silico Substrate-Binding Profiling for SARS-CoV-2 Main Protease (Mpro) Using Hexapeptide Substrates

The SARS-CoV-2 main protease (M^pro) is essential for the life cycle of the COVID-19 virus. It cleaves the two polyproteins at 11 positions to generate mature proteins for virion formation. The cleavage site on these polyproteins is known to be Leu-Gln↓(Ser/Ala/Gly). A range of hexapeptides that follow the known sequence for recognition and cleavage was constructed using RDKit libraries and complexed with the crystal structure of M^pro (PDB ID 6XHM) through extensive molecular docking calculations. A subset of 131 of these complexes underwent 20 ns molecular dynamics simulations. The analyses of the trajectories from molecular dynamics included principal component analysis (PCA), and a method to compare PCA plots from separate trajectories was developed in terms of encoding PCA progression during the simulations. The hexapeptides formed stable complexes as expected, with reproducible molecular docking of the substrates given the extensiveness of the procedure. Only Lys-Leu-Gln*** (KLQ***) sequence complexes were studied for molecular dynamics. In this subset of complexes, the PCA analysis identified four classifications of protein motions across these sequences. KLQ*** complexes illustrated the effect of changes in substrate on the active site, with implications for understanding the substrate recognition of M^pro and informing the development of small molecule inhibitors.

Dynamical Nonequilibrium Molecular Dynamics Simulations Identify Allosteric Sites and Positions Associated with Drug Resistance in the SARS-CoV-2 Main Protease

The SARS-CoV-2 main protease (Mpro) plays an essential role in the coronavirus lifecycle by catalyzing hydrolysis of the viral polyproteins at specific sites. Mpro is the target of drugs, such as nirmatrelvir, though resistant mutants have emerged that threaten drug efficacy. Despite its importance, questions remain on the mechanism of how Mpro binds its substrates. Here, we apply dynamical nonequilibrium molecular dynamics (D-NEMD) simulations to evaluate structural and dynamical responses of Mpro to the presence and absence of a substrate. The results highlight communication between the Mpro dimer subunits and identify networks, including some far from the active site, that link the active site with a known allosteric inhibition site, or which are associated with nirmatrelvir resistance. They imply that some mutations enable resistance by altering the allosteric behavior of Mpro. More generally, the results show the utility of the D-NEMD technique for identifying functionally relevant allosteric sites and networks including those relevant to resistance.

Targeting SARS-CoV-2 Main Protease: A Successful Story Guided by an In Silico Drug Repurposing Approach

The SARS-CoV-2 main protease (Mpro) is a crucial enzyme for viral replication and has been considered an attractive drug target for the treatment of COVID-19. In this study, virtual screening techniques and in vitro assays were combined to identify novel Mpro inhibitors starting from around 8000 FDA-approved drugs. The docking analysis highlighted 17 promising best hits, biologically characterized in terms of their Mpro inhibitory activity. Among them, 7 cephalosporins and the oral anticoagulant betrixaban were able to block the enzyme activity in the micromolar range with no cytotoxic effect at the highest concentration tested. After the evaluation of the degree of conservation of Mpro residues involved in the binding with the studied ligands, the ligands' activity on SARS-CoV-2 replication was assessed. The ability of betrixaban to affect SARS-CoV-2 replication associated to its antithrombotic effect could pave the way for its possible use in the treatment of hospitalized COVID-19 patients.

An interaction-based drug discovery screen explains known SARS-CoV-2 inhibitors and predicts new compound scaffolds

The recent outbreak of the COVID-19 pandemic caused by severe acute respiratory syndrome-Coronavirus-2 (SARS-CoV-2) has shown the necessity for fast and broad drug discovery methods to enable us to react quickly to novel and highly infectious diseases. A well-known SARS-CoV-2 target is the viral main 3-chymotrypsin-like cysteine protease (M^pro), known to control coronavirus replication, which is essential for the viral life cycle. Here, we applied an interaction-based drug repositioning algorithm on all protein-compound complexes available in the protein database (PDB) to identify M^pro inhibitors and potential novel compound scaffolds against SARS-CoV-2. The screen revealed a heterogeneous set of 692 potential M^pro inhibitors containing known ones such as Dasatinib, Amodiaquine, and Flavin mononucleotide, as well as so far untested chemical scaffolds. In a follow-up evaluation, we used publicly available data published almost two years after the screen to validate our results. In total, we are able to validate 17% of the top 100 predictions with publicly available data and can furthermore show that predicted compounds do cover scaffolds that are yet not associated with M^pro. Finally, we detected a potentially important binding pattern consisting of 3 hydrogen bonds with hydrogen donors of an oxyanion hole within the active side of M^pro. Overall, these results give hope that we will be better prepared for future pandemics and that drug development will become more efficient in the upcoming years.

Contribution of the catalytic dyad of SARS-CoV-2 main protease to binding covalent and noncovalent inhibitors

The effect of mutations of the catalytic dyad residues of SARS-CoV-2 main protease (MPro^WT) on the thermodynamics of binding of covalent inhibitors comprising nitrile [nirmatrelvir (NMV), NBH2], aldehyde (GC373) and ketone (BBH1) warheads to MPro is examined together with room temperature X-ray crystallography. When lacking the nucleophilic C145, NMV binding is ∼400-fold weaker corresponding to 3.5 kcal/mol and 13.3 °C decreases in free energy (ΔG) and thermal stability (T_m), respectively, relative to MPro^WT. The H41A mutation results in a 20-fold increase in the dissociation constant (K_d), and 1.7 kcal/mol and 1.4 °C decreases in ΔG and T_m, respectively. Increasing the pH from 7.2 to 8.2 enhances NMV binding to MPro^H41A, whereas no significant change is observed in binding to MPro^WT. Structures of the four inhibitor complexes with MPro^1-304/C145A show that the active site geometries of the complexes are nearly identical to that of MPro^WT with the nucleophilic sulfur of C145 positioned to react with the nitrile or the carbonyl carbon. These results support a two-step mechanism for the formation of the covalent complex involving an initial non-covalent binding followed by a nucleophilic attack by the thiolate anion of C145 on the warhead carbon. Noncovalent inhibitor ensitrelvir (ESV) exhibits a binding affinity to MPro^WT that is similar to NMV but differs in its thermodynamic signature from NMV. The binding of ESV to MPro^C145A also results in a significant, but smaller, increase in K_d and decrease in ΔG and T_m, relative to NMV.

A Tale of Two Proteases: MPro and TMPRSS2 as Targets for COVID-19 Therapies

Considering the importance of the 2019 outbreak of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) resulting in the coronavirus disease 2019 (COVID-19) pandemic, an overview of two proteases that play an important role in the infection by SARS-CoV-2, the main protease of SARS-CoV-2 (M^Pro) and the host transmembrane protease serine 2 (TMPRSS2), is presented in this review. After summarising the viral replication cycle to identify the relevance of these proteases, the therapeutic agents already approved are presented. Then, this review discusses some of the most recently reported inhibitors first for the viral M^Pro and next for the host TMPRSS2 explaining the mechanism of action of each protease. Afterward, some computational approaches to design novel M^Pro and TMPRSS2 inhibitors are presented, also describing the corresponding crystallographic structures reported so far. Finally, a brief discussion on a few reports found some dual-action inhibitors for both proteases is given. This review provides an overview of two proteases of different origins (viral and human host) that have become important targets for the development of antiviral agents to treat COVID-19.

Unraveling the catalytic mechanism of SARS-CoV-2 papain-like protease with allosteric modulation of C270 mutation using multiscale computational approaches

Papain-like protease (PL^pro) is a promising therapeutic target against SARS-CoV-2, but its restricted S1/S2 subsites pose an obstacle in developing active site-directed inhibitors. We have recently identified C270 as a novel covalent allosteric site for SARS-CoV-2 PL^pro inhibitors. Here we present a theoretical investigation of the proteolysis reaction catalyzed by the wild-type SARS-CoV-2 PL^pro as well as the C270R mutant. Enhanced sampling MD simulations were first performed to explore the influence of C270R mutation on the protease dynamics, and sampled thermodynamically favorable conformations were then submitted to MM/PBSA and QM/MM MD simulations for thorough characterization of the protease-substrate binding and covalent reactions. The disclosed proteolysis mechanism of PL^pro, as characterized by the occurrence of proton transfer from the catalytic C111 to H272 prior to the substrate binding and with deacylation being the rate-determining step of the whole proteolysis process, is not completely identical to that of the 3C-like protease, another key cysteine protease of coronaviruses. The C270R mutation alters the structural dynamics of the BL2 loop that indirectly impairs the catalytic function of H272 and reduces the binding of the substrate with the protease, ultimately showing an inhibitory effect on PL^pro. Together, these results provide a comprehensive understanding at the atomic level of the key aspects of SARS-CoV-2 PL^pro proteolysis, including the catalytic activity allosterically regulated by C270 modification, which is crucial to the follow-up inhibitor design and development.

Dock-able linear and homodetic di, tri, tetra and pentapeptide library from canonical amino acids: SARS-CoV-2 Mpro as a case study

Peptide-based therapeutics are increasingly pushing to the forefront of biomedicine with their promise of high specificity and low toxicity. Although noncanonical residues can always be used, employing only the natural 20 residues restricts the chemical space to a finite dimension allowing for comprehensive in silico screening. Towards this goal, the dataset comprising all possible di-, tri-, and tetra-peptide combinations of the canonical residues has been previously reported. However, with increasing computational power, the comprehensive set of pentapeptides is now also feasible for screening as the comprehensive set of cyclic peptides comprising four or five residues. Here, we provide both the complete and prefiltered libraries of all di-, tri-, tetra-, and penta-peptide sequences from 20 canonical amino acids and their homodetic (N-to-C-terminal) cyclic homologues. The FASTA, simplified molecular-input line-entry system (SMILES), and structure-data file (SDF)-three dimension (3D) libraries can be readily used for screening against protein targets. We also provide a simple method and tool for conducting identity-based filtering. Access to this dataset will accelerate small peptide screening workflows and encourage their use in drug discovery campaigns. As a case study, the developed library was screened against severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) main protease to identify potential small peptide inhibitors.

SARS-CoV-2 Mpro Protease Variants of Concern Display Altered Viral Substrate and Cell Host Target Galectin-8 Processing but Retain Sensitivity toward Antivirals

The main protease of SARS-CoV-2 (M^pro) is the most promising drug target against coronaviruses due to its essential role in virus replication. With newly emerging variants there is a concern that mutations in M^pro may alter the structural and functional properties of protease and subsequently the potency of existing and potential antivirals. We explored the effect of 31 mutations belonging to 5 variants of concern (VOCs) on catalytic parameters and substrate specificity, which revealed changes in substrate binding and the rate of cleavage of a viral peptide. Crystal structures of 11 M^pro mutants provided structural insight into their altered functionality. Additionally, we show M^pro mutations influence proteolysis of an immunomodulatory host protein Galectin-8 (Gal-8) and a subsequent significant decrease in cytokine secretion, providing evidence for alterations in the escape of host-antiviral mechanisms. Accordingly, mutations associated with the Gamma VOC and highly virulent Delta VOC resulted in a significant increase in Gal-8 cleavage. Importantly, IC50s of nirmatrelvir (Pfizer) and our irreversible inhibitor AVI-8053 demonstrated no changes in potency for both drugs for all mutants, suggesting M^pro will remain a high-priority antiviral drug candidate as SARS-CoV-2 evolves.

SARS-CoV-2 polyprotein substrate regulates the stepwise Mpro cleavage reaction

The processing of the Coronavirus polyproteins pp1a and pp1ab by the main protease M^pro to produce mature proteins is a crucial event in virus replication and a promising target for antiviral drug development. M^pro cleaves polyproteins in a defined order, but how M^pro and/or the polyproteins determine the order of cleavage remains enigmatic due to a lack of structural information about polyprotein-bound M^pro. Here, we present the cryo-EM structures of SARS-CoV-2 M^pro in an apo form and in complex with the nsp7-10 region of the pp1a polyprotein. The complex structure shows that M^pro interacts with only the recognition site residues between nsp9 and nsp10, without any association with the rest of the polyprotein. Comparison between the apo form and polyprotein-bound structures of M^pro highlights the flexible nature of the active site region of M^pro, which allows it to accommodate 10 recognition sites found in the polyprotein. These observations suggest that the role of M^pro in selecting a preferred cleavage site is limited and underscore the roles of the structure, conformation and/or dynamics of the polyproteins in determining the sequence of polyprotein cleavage by M^pro.

Discovery of Potential Inhibitors of SARS-CoV-2 Main Protease by a Transfer Learning Method

The COVID-19 pandemic caused by SARS-CoV-2 remains a global public health threat and has prompted the development of antiviral therapies. Artificial intelligence may be one of the strategies to facilitate drug development for emerging and re-emerging diseases. The main protease (Mpro) of SARS-CoV-2 is an attractive drug target due to its essential role in the virus life cycle and high conservation among SARS-CoVs. In this study, we used a data augmentation method to boost transfer learning model performance in screening for potential inhibitors of SARS-CoV-2 Mpro. This method appeared to outperform graph convolution neural network, random forest and Chemprop on an external test set. The fine-tuned model was used to screen for a natural compound library and a de novo generated compound library. By combination with other in silico analysis methods, a total of 27 compounds were selected for experimental validation of anti-Mpro activities. Among all the selected hits, two compounds (gyssypol acetic acid and hyperoside) displayed inhibitory effects against Mpro with IC50 values of 67.6 μM and 235.8 μM, respectively. The results obtained in this study may suggest an effective strategy of discovering potential therapeutic leads for SARS-CoV-2 and other coronaviruses.

Diphenyl Diselenide and SARS-CoV-2: in silico Exploration of the Mechanisms of Inhibition of Main Protease (Mpro) and Papain-like Protease (PLpro)

The SARS-CoV-2 pandemic has prompted global efforts to develop therapeutics. The main protease of SARS-CoV-2 (M^pro) and the papain-like protease (PL^pro) are essential for viral replication and are key targets for therapeutic development. In this work, we investigate the mechanisms of SARS-CoV-2 inhibition by diphenyl diselenide (PhSe)₂ which is an archetypal model of diselenides and a renowned potential therapeutic agent. The in vitro inhibitory concentration of (PhSe)₂ against SARS-CoV-2 in Vero E6 cells falls in the low micromolar range. Molecular dynamics (MD) simulations and density functional theory (DFT) calculations [level of theory: SMD-B3LYP-D3(BJ)/6-311G(d,p), cc-pVTZ] are used to inspect non-covalent inhibition modes of both proteases via π-stacking and the mechanism of covalent (PhSe)₂ + M^pro product formation involving the catalytic residue C145, respectively. The in vitro CC₅₀ (24.61 μM) and EC₅₀ (2.39 μM) data indicate that (PhSe)₂ is a good inhibitor of the SARS-CoV-2 virus replication in a cell culture model. The in silico findings indicate potential mechanisms of proteases' inhibition by (PhSe)₂; in particular, the results of the covalent inhibition here discussed for M^pro, whose thermodynamics is approximatively isoergonic, prompt further investigation in the design of antiviral organodiselenides.

An in-solution snapshot of SARS-COV-2 main protease maturation process and inhibition

The main protease from SARS-CoV-2 (Mpro) is responsible for cleavage of the viral polyprotein. Mpro self-processing is called maturation, and it is crucial for enzyme dimerization and activity. Here we use C145S Mpro to study the structure and dynamics of N-terminal cleavage in solution. Native mass spectroscopy analysis shows that mixed oligomeric states are composed of cleaved and uncleaved particles, indicating that N-terminal processing is not critical for dimerization. A 3.5 Å cryo-EM structure provides details of Mpro N-terminal cleavage outside the constrains of crystal environment. We show that different classes of inhibitors shift the balance between oligomeric states. While non-covalent inhibitor MAT-POS-e194df51-1 prevents dimerization, the covalent inhibitor nirmatrelvir induces the conversion of monomers into dimers, even with intact N-termini. Our data indicates that the Mpro dimerization is triggered by induced fit due to covalent linkage during substrate processing rather than the N-terminal processing.

Virtual screening of natural products inspired in-house library to discover potential lead molecules against the SARS-CoV-2 main protease

SARS-CoV-2, a new coronavirus emerged in 2019, causing a global healthcare epidemic. Although a variety of drug targets have been identified as potential antiviral therapies, and effective candidate against SARS-CoV-2 remains elusive. One of the most promising targets for combating COVID-19 is SARS-CoV-2 Main protease (M^pro, a protein responsible for viral replication. In this work, an in-house curated library was thoroughly evaluated for druggability against M^pro. We identified four ligands (FG, Q5, P5, and PJ4) as potential inhibitors based on docking scores, predicted binding energies (MMGBSA), in silico ADME, and RMSD trajectory analysis. Among the selected ligands, FG, a natural product from Andrographis nallamalayana, exhibited the highest binding energy of -10.31 kcal/mol close to the docking score of clinical candidates Boceprevir and GC376. Other ligands (P5, natural product from cardiospermum halicacabum and two synthetic molecules Q5 and PJ4) have shown comparable docking scores ranging -7.65 kcal/mol to -7.18 kcal/mol. Interestingly, we found all four top ligands had Pi bond interaction with the main amino acid residues HIS41 and CYS145 (catalytic dyad), H-bonding interactions with GLU166, ARG188, and GLN189, and hydrophobic interactions with MET49 and MET165 in the binding site of M^pro. According to the ADME analysis, Q5 and P5 are within the acceptable range of drug likeliness, compared to FG and PJ4. The interaction stability of the lead molecules with viral protease was verified using replicated MD simulations. Thus, the present study opens up the opportunity of developing drug candidates targeting SARS-CoV-2 main protease (M^pro) to mitigate the disease.Communicated by Ramaswamy H. Sarma.