Structural basis for replicase polyprotein cleavage and substrate specificity of main protease from SARS-CoV-2

Crystal structures of coronaviral main proteases in complex with the non-covalent inhibitor X77

Main proteases (Mpros) are a class of conserved cysteine hydrolases among coronaviruses and play a crucial role in viral replication. Therefore, Mpros are ideal targets for the development of pan-coronavirus drugs. X77, previously developed against SARS-CoV Mpro, was repurposed as a non-covalent tight binder inhibitor against SARS-CoV-2 Mpro during COVID-19 pandemic. Many novel inhibitors with favorable efficacy have been discovered using X77 as a reference, suggesting that X77 could be a valuable scaffold for drug design. However, the broad-spectrum performance of X77 and underlying mechanism remain less understood. Here, we reported the crystal structures of Mpros from SARS-CoV-2, SARS-CoV, and MERS-CoV, and several Mpro mutants from SARS-CoV-2 variants bound to X77. A detailed analysis of these structures revealed key structural determinants essential for interaction and elucidated the binding modes of X77 with different coronaviral Mpros. The potencies of X77 against these investigated Mpros were further evaluated through molecular dynamic simulation and binding free energy calculation. These data provide molecular insights into broad-spectrum inhibition against coronaviral Mpros by X77 and the similarities and differences of X77 when bound to various Mpros, which will promote X77-based design of novel antivirals with broad-spectrum efficacy against different coronaviruses and SARS-CoV-2 variants.

3-chymotrypsin-like protease in SARS-CoV-2

Coronaviruses constitute a significant threat to the human population. Severe acute respiratory syndrome coronavirus-2, SARS-CoV-2, is a highly pathogenic human coronavirus that has caused the coronavirus disease 2019 (COVID-19) pandemic. It has led to a global viral outbreak with an exceptional spread and a high death toll, highlighting the need for effective antiviral strategies. 3-Chymotrypsin-like protease (3CLpro), the main protease in SARS-CoV-2, plays an indispensable role in the SARS-CoV-2 viral life cycle by cleaving the viral polyprotein to produce 11 individual non-structural proteins necessary for viral replication. 3CLpro is one of two proteases that function to produce new viral particles. It is a highly conserved cysteine protease with identical structural folds in all known human coronaviruses. Inhibitors binding with high affinity to 3CLpro will prevent the cleavage of viral polyproteins, thus impeding viral replication. Multiple strategies have been implemented to screen for inhibitors against 3CLpro, including peptide-like and small molecule inhibitors that covalently and non-covalently bind the active site, respectively. In addition, allosteric sites of 3CLpro have been identified to screen for small molecules that could make non-competitive inhibitors of 3CLpro. In essence, this review serves as a comprehensive guide to understanding the structural intricacies and functional dynamics of 3CLpro, emphasizing key findings that elucidate its role as the main protease of SARS-CoV-2. Notably, the review is a critical resource in recognizing the advancements in identifying and developing 3CLpro inhibitors as effective antiviral strategies against COVID-19, some of which are already approved for clinical use in COVID-19 patients.

Structure-Guided Design of Potent Coronavirus Inhibitors with a 2-Pyrrolidone Scaffold: Biochemical, Crystallographic, and Virological Studies

Zoonotic coronaviruses are known to produce severe infections in humans and have been the cause of significant morbidity and mortality worldwide. SARS-CoV-2 was the largest and latest contributor of fatal cases, even though MERS-CoV has the highest case-fatality ratio among zoonotic coronaviruses. These infections pose a high risk to public health worldwide warranting efforts for the expeditious discovery of antivirals. Hence, we hereby describe a novel series of inhibitors of coronavirus 3CLpro embodying an N-substituted 2-pyrrolidone scaffold envisaged to exploit favorable interactions with the S3-S4 subsites and connected to an invariant Leu-Gln P2-P1 recognition element. Several inhibitors showed nanomolar antiviral activity in enzyme and cell-based assays, with no significant cytotoxicity. High-resolution crystal structures of inhibitors bound to the 3CLpro were determined to probe and identify the molecular determinants associated with binding, to inform the structure-guided optimization of the inhibitors, and to confirm the mechanism of action of the inhibitors.

Computational Insights into SARS-CoV-2 Main Protease Mutations and Nirmatrelvir Efficacy: The Effects of P132H and P132H-A173V

Nirmatrelvir, a pivotal component of the oral antiviral Paxlovid for COVID-19, targets the SARS-CoV-2 main protease (Mpro) as a covalent inhibitor. Here, we employed combined computational methods to explore how the prevalent Omicron variant mutation P132H, alone and in combination with A173V (P132H-A173V), affects nirmatrelvir's efficacy. Our findings suggest that P132H enhances the noncovalent binding affinity of Mpro for nirmatrelvir, whereas P132H-A173V diminishes it. Although both mutants catalyze the rate-limiting step more efficiently than the wild-type (WT) Mpro, P132H slows the overall rate of covalent bond formation, whereas P132H-A173V accelerates it. Comprehensive analysis of noncovalent and covalent contributions to the overall binding free energy of the covalent complex suggests that P132H likely enhances Mpro sensitivity to nirmatrelvir, while P132H-A173V may confer resistance. Per-residue decompositions of the binding and activation free energies pinpoint key residues that significantly affect the binding affinity and reaction rates, revealing how the mutations modulate these effects. The mutation-induced conformational perturbations alter drug-protein local contact intensities and the electrostatic preorganization of the protein, affecting noncovalent binding affinity and the stability of key reaction states, respectively. Our findings inform the mechanisms of nirmatrelvir resistance and sensitivity, facilitating improved drug design and the detection of resistant strains.

Mechanistic insights into bismuth(iii) inhibition of SARS-CoV-2 helicase

The COVID-19 pandemic caused by SARS-CoV-2 resulted in a global public health crisis. In addition to vaccines, the development of effective therapy is highly desirable. Targeting a protein that plays a critical role in virus replication may allow pan-spectrum antiviral drugs to be developed. Among SARS-CoV-2 proteins, helicase (i.e., non-structural protein 13) is considered as a promising antiviral drug target due to its highly conserved sequence, unique structure and function. Herein, we demonstrate SARS-CoV-2 helicase as a target of bismuth-based antivirals in virus-infected mammalian cells by a metal-tagged antibody approach. To search for more potent bismuth-based antivirals, we further screened a panel of bismuth compounds towards inhibition of ATPase and DNA unwinding activity of nsp13 and identified a highly potent bismuth compound Bi(5-aminotropolonate)3, namely Bi(Tro-NH2)3 with an IC50 of 30 nM for ATPase. We show that bismuth-based compounds inhibited nsp13 unwinding activity via disrupting the binding of ATP and the DNA substrate to viral helicase. Binding of Bi(iii) to nsp13 also abolished the interaction between nsp12 and nsp13 as evidenced by immunofluorescence and co-immunoprecipitation assays. Finally, we validate our in vitro data in SARS-CoV-2 infected mammalian cells. Notably, Bi(6-TG)3 exhibited an EC50 of 1.18 ± 0.09 μM with a selective index of 847 in VeroE6-TMPRSS2 infected cells. This study highlights the important role of helicase for the development of more effective antiviral drugs to combat SARS-CoV-2 infection.

A SARS-CoV-2 Mpro fluorescent sensor for exploring pharmacodynamic substances from traditional Chinese medicine

The main protease of SARS-CoV-2 (SARS-CoV-2 Mpro) plays a critical role in the replication and life cycle of the virus. Currently, how to screen SARS-CoV-2 Mpro inhibitors from complex traditional Chinese medicine (TCM) is the bottleneck for exploring the pharmacodynamic substances of TCM against SARS-CoV-2. In this study, a simple, cost-effective, rapid, and selective fluorescent sensor (TPE-S-TLG sensor) was designed with an AIE (aggregation-induced emission) probe (TPE-Ph-In) and the SARS-CoV-2 Mpro substrate (S-TLG). The TPE-S-TLG sensor was characterized using UV-Vis absorption spectroscopy, fluorescence spectroscopy, dynamic light scattering (DLS), transmission electron microscopy (TEM), zeta potential, and Fourier transform infrared (FTIR) spectroscopy techniques. The limit of detection of this method to detect SARS-CoV-2 Mpro was measured to be 5 ng mL-1. Furthermore, the TPE-S-TLG sensor was also successfully applied to screen Mpro inhibitors from Xuebijing injection using the separation and collection of the HPLC-fully automatic partial fraction collector (HPLC-FC). Six active compounds, including protocatechualdehyde, chlorogenic acid, hydroxysafflower yellow A, caffeic acid, isoquercetin, and pentagalloylglucose, were identified using UHPLC-Q-TOF/MS that could achieve 90% of the Mpro inhibition rate for the Xuebijing injection. Accordingly, the strategy can be broadly applied in the detection of disease-related proteases as well as screening active substances from TCM.

Computational evaluation and benchmark study of 342 crystallographic holo-structures of SARS-CoV-2 Mpro enzyme

The coronavirus disease 19 pandemic, caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has led to a global health crisis with millions of confirmed cases and related deaths. The main protease (Mpro) of SARS-CoV-2 is crucial for viral replication and presents an attractive target for drug development. Despite the approval of some drugs, the search for effective treatments continues. In this study, we systematically evaluated 342 holo-crystal structures of Mpro to identify optimal conformations for structure-based virtual screening (SBVS). Our analysis revealed limited structural flexibility among the structures. Three docking programs, AutoDock Vina, rDock, and Glide were employed to assess the efficiency of virtual screening, revealing diverse performances across selected Mpro structures. We found that the structures 5RHE, 7DDC, and 7DPU (PDB Ids) consistently displayed the lowest EF, AUC, and BEDROCK scores. Furthermore, these structures demonstrated the worst pose prediction results in all docking programs. Two structural differences contribute to variations in docking performance: the absence of the S1 subsite in 7DDC and 7DPU, and the presence of a subpocket in the S2 subsite of 7DDC, 7DPU, and 5RHE. These findings underscore the importance of selecting appropriate Mpro conformations for SBVS, providing valuable insights for advancing drug discovery efforts.

Optimization of quenched fluorescent peptide substrates of SARS-CoV-2 3CLpro main protease (Mpro) from proteomic identification of P6-P6' active site specificity

SARS-CoV-2 3C-like main protease (3CLpro) is essential for protein excision from the viral polyprotein. 3CLpro inhibitor drug development to block SARS-CoV-2 replication focuses on the catalytic non-prime (P) side for specificity and potency, but the importance of the prime (P') side in substrate specificity and for drug development remains underappreciated. We determined the P6-P6' specificity for 3CLpro from >800 cleavage sites that we identified using Proteomic Identification of Cleavage site Specificity (PICS). Cleavage occurred after the canonical P1-Gln and non-canonical P1-His and P1-Met residues. Moreover, P3 showed a preference for Arg/Lys and P3' for His. Essential H-bonds between the N-terminal Ser1 of protomer-B in 3CLpro dimers form with P1-His, but not with P1-Met. Nonetheless, cleavage occurs at P1-Met456 in native MAP4K5. Elevated reactive oxygen species in SARS-CoV-2 infection oxidize methionines. Molecular simulations revealed P1-MetOX forms an H-bond with Ser1 and notably, strong positive cooperativity between P1-Met with P3'-His was revealed, which enhanced peptide-cleavage rates. The highly plastic S3' subsite accommodates P3'-His that displays stabilizing backbone H-bonds with Thr25 lying central in a "'threonine trio" (Thr24-Thr25-Thr26) in the P'-binding domain I. Molecular docking simulations unveiled structure-activity relationships impacting 3CLpro-substrate interactions, and the role of these structural determinants was confirmed by MALDI-TOF-MS cleavage assays of P1'- and P3'-positional scanning peptide libraries carrying a 2nd optimal cut-site as an internal positive control. These data informed the design of two new and highly soluble 3CLproquenched-fluorescent peptide substrates for improved FRET monitoring of 3CLpro activity with 15× improved sensitivity over current assays.IMPORTANCEFrom global proteomics identification of >800 cleavage sites, we characterized the P6-P6' active site specificity of SARS-CoV-2 3CLpro using proteome-derived peptide library screens, molecular modeling simulations, and focussed positional peptide libraries. In P1', we show that alanine and serine are cleaved 3× faster than glycine and the hydrophobic small amino acids Leu, Ile, or Val prevent cleavage of otherwise optimal non-prime sequences. In characterizing non-canonical non-prime P1 specificity, we explored the unusual P1-Met specificity, discovering enhanced cleavage when in the oxidized state (P1-MetOX). We unveiled unexpected amino acid cooperativity at P1-Met with P3'-His and noncanonical P1-His with P2-Phe, and the importance of the threonine trio (Thr24-Thr25-Thr26) in the prime side binding domain I in defining prime side binding in SARS-CoV-2 3CLpro. From these analyses, we rationally designed quenched-fluorescence natural amino acid peptide substrates with >15× improved sensitivity and high peptide solubility, facilitating handling and application for screening of new antiviral drugs.

Subpocket-Based Analysis Approach for the Protein Pocket Dynamics

Structural and dynamic characteristics of protein pockets remarkably influence their biological functions and are also important for enzyme engineering and new drug research and development. To date, several softwares have been developed to analyze the dynamic properties of protein pockets. However, due to the complexity and diversity of the pocket information during the kinetic relaxation, further improvement and capacity expansion of current tools are required. Here, we developed a platform software AlphaTraj in which a computational strategy that divides the whole protein pocket into subpockets and examines various properties of the subpockets such as survival time, stability, and correlation was proposed and implemented. We also proposed a scoring function for the subpockets as well as the whole pocket to visualize the quality of the pocket. Furthermore, we implemented automated conformational search functions for ligand docking and ligand optimization. These functions may help us to gain a deep understanding of the dynamic properties of protein pockets and accelerate the protein engineering and the design of inhibitors and small-molecule drugs. The software is freely available at https://github.com/dooo12332/AlphaTraj.git under the GNU GPL license.

3-Chymotrypsin-like Protease (3CLpro) of SARS-CoV-2: Validation as a Molecular Target, Proposal of a Novel Catalytic Mechanism, and Inhibitors in Preclinical and Clinical Trials

Proteases represent common targets in combating infectious diseases, including COVID-19. The 3-chymotrypsin-like protease (3CLpro) is a validated molecular target for COVID-19, and it is key for developing potent and selective inhibitors for inhibiting viral replication of SARS-CoV-2. In this review, we discuss structural relationships and diverse subsites of 3CLpro, shedding light on the pivotal role of dimerization and active site architecture in substrate recognition and catalysis. Our analysis of bioinformatics and other published studies motivated us to investigate a novel catalytic mechanism for the SARS-CoV-2 polyprotein cleavage by 3CLpro, centering on the triad mechanism involving His41-Cys145-Asp187 and its indispensable role in viral replication. Our hypothesis is that Asp187 may participate in modulating the pKa of the His41, in which catalytic histidine may act as an acid and/or a base in the catalytic mechanism. Recognizing Asp187 as a crucial component in the catalytic process underscores its significance as a fundamental pharmacophoric element in drug design. Next, we provide an overview of both covalent and non-covalent inhibitors, elucidating advancements in drug development observed in preclinical and clinical trials. By highlighting various chemical classes and their pharmacokinetic profiles, our review aims to guide future research directions toward the development of highly selective inhibitors, underscore the significance of 3CLpro as a validated therapeutic target, and propel the progression of drug candidates through preclinical and clinical phases.

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.

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.

Efficient synthesis of novel colchicine-magnolol hybrids and evaluation of their inhibitory activity on key proteases of 2019-nCoV replication and acute lung injury

Coronavirus disease 2019 (COVID-19), caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2 or 2019-nCoV), is a life-threatening infectious condition. Acute lung injury is a common complication in patients with COVID-19. 3-chymotrypsin-like protease (3CLpro) of 2019-nCoV and neutrophil elastase are critical targets of COVID-19 and acute lung injury, respectively. Colchicine and magnolol are reported to exert inhibitory effects on inflammatory response, the severe comorbidity in both COVID-19 and acute lung injury. We thus designed and synthesized a series of novel colchicine-magnolol hybrids based on a two-step synthetic sequence. It was found that these novel hybrids provided unexpected inhibition on 3CLpro and neutrophil elastase, a bioactivity that colchicine and magnolol did not possess. These findings not only provide perquisites for further in vitro and in vivo investigation to confirm the therapeutic potentiality of novel colchicine-magnolol hybrids, but also suggest that the concurrent inhibition of 3CLpro and neutrophil elastase may enable novel colchicine-magnolol hybrids as effective multi-target drug compounds.

Quantitative Assessment of Energetic Contributions of Residues in a SARS-CoV-2 Viral Enzyme/Nanobody Interface

The highly conserved protease enzyme from SARS-CoV-2 (MPro) is crucial for viral replication and is an attractive target for the design of novel inhibitory compounds. MPro is known to be conformationally flexible and has been stabilized in an extended conformation in a complex with a novel nanobody (NB2B4), which inhibits the dimerization of the enzyme via binding to an allosteric site. However, the energetic contributions of the nanobody residues stabilizing the MPro/nanobody interface remain unresolved. We probed these residues using all-atom MD simulations in combination with alchemical free energy calculations by studying the physical residue-residue interactions and discovered the role of hydrophobic and electrostatic interactions in stabilizing the complex. Specifically, we found via mutational analysis that three interfacial nanobody residues (Y59, R106, and L109) contributed significantly, two residues (L107 and P110) contributed moderately, and two residues (H112 and T113) contributed minimally to the overall binding affinity of the nanobody. We also discovered that the nanobody affinity could be enhanced via a charge-reversal mutation (D62R) that alters the local interfacial electrostatic environment of this residue in the complex. These findings are potentially useful in designing novel synthetic nanobodies as allosteric inhibitors of MPro.

Bioengineered amyloid peptide for rapid screening of inhibitors against main protease of SARS-CoV-2

The coronavirus disease 2019 (COVID-19) caused by severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) has evoked a worldwide pandemic. As the emergence of variants has hampered the neutralization capacity of currently available vaccines, developing effective antiviral therapeutics against SARS-CoV-2 and its variants becomes a significant challenge. The main protease (Mpro) of SARS-CoV-2 has received increased attention as an attractive pharmaceutical target because of its pivotal role in viral replication and proliferation. Here, we generated a de novo Mpro-inhibitor screening platform to evaluate the efficacies of Mpro inhibitors based on Mpro cleavage site-embedded amyloid peptide (MCAP)-coated gold nanoparticles (MCAP-AuNPs). We fabricated MCAPs comprising an amyloid-forming sequence and Mpro-cleavage sequence, mimicking in vivo viral replication process mediated by Mpro. By measuring the proteolytic activity of Mpro and the inhibitory efficacies of various drugs, we confirmed that the MCAP-AuNP-based platform was suitable for rapid screening potential of Mpro inhibitors. These results demonstrated that our MCAP-AuNP-based platform has great potential for discovering Mpro inhibitors and may accelerate the development of therapeutics against COVID-19.

Structure-based design of pan-coronavirus inhibitors targeting host cathepsin L and calpain-1

Respiratory disease caused by coronavirus infection remains a global health crisis. Although several SARS-CoV-2-specific vaccines and direct-acting antivirals are available, their efficacy on emerging coronaviruses in the future, including SARS-CoV-2 variants, might be compromised. Host-targeting antivirals provide preventive and therapeutic strategies to overcome resistance and manage future outbreak of emerging coronaviruses. Cathepsin L (CTSL) and calpain-1 (CAPN1) are host cysteine proteases which play crucial roles in coronaviral entrance into cells and infection-related immune response. Here, two peptidomimetic α-ketoamide compounds, 14a and 14b, were identified as potent dual target inhibitors against CTSL and CAPN1. The X-ray crystal structures of human CTSL and CAPN1 in complex with 14a and 14b revealed the covalent binding of α-ketoamide groups of 14a and 14b to C25 of CTSL and C115 of CAPN1. Both showed potent and broad-spectrum anticoronaviral activities in vitro, and it is worth noting that they exhibited low nanomolar potency against SARS-CoV-2 and its variants of concern (VOCs) with EC50 values ranging from 0.80 to 161.7 nM in various cells. Preliminary mechanistic exploration indicated that they exhibited anticoronaviral activity through blocking viral entrance. Moreover, 14a and 14b exhibited good oral pharmacokinetic properties in mice, rats and dogs, and favorable safety in mice. In addition, both 14a and 14b treatments demonstrated potent antiviral potency against SARS-CoV-2 XBB 1.16 variant infection in a K18-hACE2 transgenic mouse model. And 14b also showed effective antiviral activity against HCoV-OC43 infection in a mouse model with a final survival rate of 60%. Further evaluation showed that 14a and 14b exhibited excellent anti-inflammatory effects in Raw 264.7 mouse macrophages and in mice with acute pneumonia. Taken together, these results suggested that 14a and 14b are promising drug candidates, providing novel insight into developing pan-coronavirus inhibitors with antiviral and anti-inflammatory properties.

A Broad-Spectrum Multi-Antigen mRNA/LNP-Based Pan-Coronavirus Vaccine Induced Potent Cross-Protective Immunity Against Infection and Disease Caused by Highly Pathogenic and Heavily Spike-Mutated SARS-CoV-2 Variants of Concern in the Syrian Hamster Model

The first-generation Spike-alone-based COVID-19 vaccines have successfully contributed to reducing the risk of hospitalization, serious illness, and death caused by SARS-CoV-2 infections. However, waning immunity induced by these vaccines failed to prevent immune escape by many variants of concern (VOCs) that emerged from 2020 to 2024, resulting in a prolonged COVID-19 pandemic. We hypothesize that a next-generation Coronavirus (CoV) vaccine incorporating highly conserved non-Spike SARS-CoV-2 antigens would confer stronger and broader cross-protective immunity against multiple VOCs. In the present study, we identified ten non-Spike antigens that are highly conserved in 8.7 million SARS-CoV-2 strains, twenty-one VOCs, SARS-CoV, MERS-CoV, Common Cold CoVs, and animal CoVs. Seven of the 10 antigens were preferentially recognized by CD8+ and CD4+ T-cells from unvaccinated asymptomatic COVID-19 patients, irrespective of VOC infection. Three out of the seven conserved non-Spike T cell antigens belong to the early expressed Replication and Transcription Complex (RTC) region, when administered to the golden Syrian hamsters, in combination with Spike, as nucleoside-modified mRNA encapsulated in lipid nanoparticles (LNP) (i.e., combined mRNA/LNP-based pan-CoV vaccine): (i) Induced high frequencies of lung-resident antigen-specific CXCR5+CD4+ T follicular helper (TFH) cells, GzmB+CD4+ and GzmB+CD8+ cytotoxic T cells (TCYT), and CD69+IFN-γ+TNFα+CD4+ and CD69+IFN-γ+TNFα+CD8+ effector T cells (TEFF); and (ii) Reduced viral load and COVID-19-like symptoms caused by various VOCs, including the highly pathogenic B.1.617.2 Delta variant and the highly transmittable heavily Spike-mutated XBB1.5 Omicron sub-variant. The combined mRNA/LNP-based pan-CoV vaccine could be rapidly adapted for clinical use to confer broader cross-protective immunity against emerging highly mutated and pathogenic VOCs.

Identification of SARS-CoV-2 Main Protease Inhibitors Using Chemical Similarity Analysis Combined with Machine Learning

SARS-CoV-2 Main Protease (Mpro) is an enzyme that cleaves viral polyproteins translated from the viral genome, which is critical for viral replication. Mpro is a target for anti-SARS-CoV-2 drug development. Herein, we performed a large-scale virtual screening by comparing multiple structural descriptors of reference molecules with reported anti-coronavirus activity against a library with >17 million compounds. Further filtering, performed by applying two machine learning algorithms, identified eighteen computational hits as anti-SARS-CoV-2 compounds with high structural diversity and drug-like properties. The activities of twelve compounds on Mpro's enzymatic activity were evaluated by fluorescence resonance energy transfer (FRET) assays. Compound 13 (ZINC13878776) significantly inhibited SARS-CoV-2 Mpro activity and was employed as a reference for an experimentally hit expansion. The structural analogues 13a (ZINC4248385), 13b (ZNC13523222), and 13c (ZINC4248365) were tested as Mpro inhibitors, reducing the enzymatic activity of recombinant Mpro with potency as follows: 13c > 13 > 13b > 13a. Then, their anti-SARS-CoV-2 activities were evaluated in plaque reduction assays using Vero CCL81 cells. Subtoxic concentrations of compounds 13a, 13c, and 13b displayed in vitro antiviral activity with IC50 in the mid micromolar range. Compounds 13a-c could become lead compounds for the development of new Mpro inhibitors with improved activity against anti-SARS-CoV-2.

Studies on the selectivity of the SARS-CoV-2 papain-like protease reveal the importance of the P2' proline of the viral polyprotein

The SARS-CoV-2 papain-like protease (PLpro) is an antiviral drug target that catalyzes the hydrolysis of the viral polyproteins pp1a/1ab, so releasing the non-structural proteins (nsps) 1-3 that are essential for the coronavirus lifecycle. The LXGG↓X motif in pp1a/1ab is crucial for recognition and cleavage by PLpro. We describe molecular dynamics, docking, and quantum mechanics/molecular mechanics (QM/MM) calculations to investigate how oligopeptide substrates derived from the viral polyprotein bind to PLpro. The results reveal how the substrate sequence affects the efficiency of PLpro-catalyzed hydrolysis. In particular, a proline at the P2' position promotes catalysis, as validated by residue substitutions and mass spectrometry-based analyses. Analysis of PLpro catalyzed hydrolysis of LXGG motif-containing oligopeptides derived from human proteins suggests that factors beyond the LXGG motif and the presence of a proline residue at P2' contribute to catalytic efficiency, possibly reflecting the promiscuity of PLpro. The results will help in identifying PLpro substrates and guiding inhibitor design.

Virtual Screening of Peptide Libraries: The Search for Peptide-Based Therapeutics Using Computational Tools

Over the last few decades, we have witnessed growing interest from both academic and industrial laboratories in peptides as possible therapeutics. Bioactive peptides have a high potential to treat various diseases with specificity and biological safety. Compared to small molecules, peptides represent better candidates as inhibitors (or general modulators) of key protein-protein interactions. In fact, undruggable proteins containing large and smooth surfaces can be more easily targeted with the conformational plasticity of peptides. The discovery of bioactive peptides, working against disease-relevant protein targets, generally requires the high-throughput screening of large libraries, and in silico approaches are highly exploited for their low-cost incidence and efficiency. The present review reports on the potential challenges linked to the employment of peptides as therapeutics and describes computational approaches, mainly structure-based virtual screening (SBVS), to support the identification of novel peptides for therapeutic implementations. Cutting-edge SBVS strategies are reviewed along with examples of applications focused on diverse classes of bioactive peptides (i.e., anticancer, antimicrobial/antiviral peptides, peptides blocking amyloid fiber formation).

Danshensu inhibits SARS-CoV-2 by targeting its main protease as a specific covalent inhibitor and discovery of bifunctional compounds eliciting antiviral and anti-inflammatory activity

The coronavirus disease 2019 (COVID-19) caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has posed a serious threat to human. Since there are still no effective treatment options against the new emerging variants of SARS-CoV-2, it is necessary to devote a continuous endeavor for more targeted drugs and the preparation for the next pandemic. Salvia miltiorrhiza and its active ingredients possess wide antiviral activities, including against SARS-CoV-2. Danshensu, as one of the most important active ingredients in Salvia miltiorrhiza, has been reported to inhibit the entry of SARS-CoV-2 into ACE2 (angiotensin-converting enzyme 2)-overexpressed HEK-293T cells and Vero-E6 cells. However, there is a paucity of information regarding its detailed target and mechanism against SARS-CoV-2. Here, we present Danshensu as a covalent inhibitor of 3-chymotrypsin-like protease (3CLpro) against SARS-CoV-2 by the time-dependent inhibition assay (TDI) and mass spectrometry analysis. Further molecular docking, site-directed mutagenesis, circular dichroism (CD) and fluorescence spectra revealed that Danshensu covalently binds to C145 of SARS-CoV-2 3CLpro, meanwhile forming the hydrogen bonds with S144, H163 and E166 in the S1 site. Structure-based optimization of Danshensu led to the discovery of the promising compounds with good inhibitory activity and microsomal stability in vitro. Due to Danshensu inhibiting lung inflammation in the mouse model, we found that Danshensu derivatives also showed better anti-inflammatory activity than Danshensu in lipopolysaccharide (LPS)-stimulated RAW264.7 macrophage cells. Thus, our study provides not only the clue of the efficacy of Salvia miltiorrhiza against SARS-CoV-2, but also a detailed mechanistic insight into the covalent mode of action of Danshensu for design of covalent inhibitors against SARS-CoV-2 3CLpro, highlighting its potential as a bifunctional molecule with antivirus and anti-inflammation.

Identification of amentoflavone as a potent SARS-CoV-2 Mpro inhibitor: a combination of computational studies and in vitro biological evaluation

Small-molecule inhibitors of SARS-CoV-2 Mpro that block the active site pocket of the viral main protease have been considered potential therapeutics for the development of drugs against SARS-CoV-2. Here, we report the identification of amentoflavone (a biflavonoid) through docking-based virtual screening of a library comprised of 231 compounds consisting of flavonoids and isoflavonoids. The docking results were further substantiated through extensive analysis of the data obtained from all-atom 150 ns MD simulation. End-state effective free energy calculations using MM-PBSA calculations further suggested that (Ra)-amentoflavone (C3'-C8''-atropisomer) may show a greater binding affinity towards the Mpro than (Sa)-amentoflavone. In vitro cytotoxicity assay established that amentoflavone showed a high CC50 value indicating much lower toxicity. Further, potent inhibition of the Mpro by amentoflavone was established by studying the effect on HEK293T cells treated with SARS-CoV-2 Mpro expressing plasmid.Communicated by Ramaswamy H. Sarma.

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.

Medicinal chemistry strategies towards the development of non-covalent SARS-CoV-2 Mpro inhibitors

The main protease (Mpro) of SARS-CoV-2 is an attractive target in anti-COVID-19 therapy for its high conservation and major role in the virus life cycle. The covalent Mpro inhibitor nirmatrelvir (in combination with ritonavir, a pharmacokinetic enhancer) and the non-covalent inhibitor ensitrelvir have shown efficacy in clinical trials and have been approved for therapeutic use. Effective antiviral drugs are needed to fight the pandemic, while non-covalent Mpro inhibitors could be promising alternatives due to their high selectivity and favorable druggability. Numerous non-covalent Mpro inhibitors with desirable properties have been developed based on available crystal structures of Mpro. In this article, we describe medicinal chemistry strategies applied for the discovery and optimization of non-covalent Mpro inhibitors, followed by a general overview and critical analysis of the available information. Prospective viewpoints and insights into current strategies for the development of non-covalent Mpro inhibitors are also discussed.

Targeted degrader technologies as prospective SARS-CoV-2 therapies

COVID-19 remains a severe public health threat despite the WHO declaring an end to the public health emergency in May 2023. Continual development of SARS-CoV-2 variants with resistance to vaccine-induced or natural immunity necessitates constant vigilance as well as new vaccines and therapeutics. Targeted protein degradation (TPD) remains relatively untapped in antiviral drug discovery and holds the promise of attenuating viral resistance development. From a unique structural design perspective, this review covers antiviral degrader merits and challenges by highlighting key coronavirus protein targets and their co-crystal structures, specifically illustrating how TPD strategies can refine existing SARS-CoV-2 3CL protease inhibitors to potentially produce superior protease-degrading agents.

Identification of Ebselen derivatives as novel SARS-CoV-2 main protease inhibitors: Design, synthesis, biological evaluation, and structure-activity relationships exploration

The main protease (Mpro) represents one of the most effective and attractive targets for designing anti-SARS-CoV-2 drugs. In this study, we designed and synthesized a novel series of Ebselen derivatives by incorporating privileged fragments from different pockets of the Mpro active site. Among these compounds, 11 compounds showed submicromolar activity in the FRET-based SARS-CoV-2 Mpro inhibition assay, with IC50 values ranging from 233 nM to 550 nM. Notably, compound 3a displayed submicromolar Mpro activity (IC50 = 364 nM) and low micromolar antiviral activity (EC50 = 8.01 µM), comparable to that of Ebselen (IC50 = 339 nM, EC50 = 3.78 µM). Time-dependent inhibition assay confirmed that these compounds acted as covalent inhibitors. Taken together, our optimization campaigns thoroughly explored the structural diversity of Ebselen and verified the impact of specific modifications on potency against Mpro.

Unsupervised deep learning for molecular dynamics simulations: a novel analysis of protein-ligand interactions in SARS-CoV-2 Mpro

Molecular dynamics (MD) simulations, which are central to drug discovery, offer detailed insights into protein-ligand interactions. However, analyzing large MD datasets remains a challenge. Current machine-learning solutions are predominantly supervised and have data labelling and standardisation issues. In this study, we adopted an unsupervised deep-learning framework, previously benchmarked for rigid proteins, to study the more flexible SARS-CoV-2 main protease (Mpro). We ran MD simulations of Mpro with various ligands and refined the data by focusing on binding-site residues and time frames in stable protein conformations. The optimal descriptor chosen was the distance between the residues and the center of the binding pocket. Using this approach, a local dynamic ensemble was generated and fed into our neural network to compute Wasserstein distances across system pairs, revealing ligand-induced conformational differences in Mpro. Dimensionality reduction yielded an embedding map that correlated ligand-induced dynamics and binding affinity. Notably, the high-affinity compounds showed pronounced effects on the protein's conformations. We also identified the key residues that contributed to these differences. Our findings emphasize the potential of combining unsupervised deep learning with MD simulations to extract valuable information and accelerate drug discovery.

The S1'-S3' Pocket of the SARS-CoV-2 Main Protease Is Critical for Substrate Selectivity and Can Be Targeted with Covalent Inhibitors

The main protease (Mpro ) of SARS-CoV-2 is a well-characterized target for antiviral drug discovery. To date, most antiviral drug discovery efforts have focused on the S4-S1' pocket of Mpro ; however, it is still unclear whether the S1'-S3' pocket per se can serve as a new site for drug discovery. In this study, the S1'-S3' pocket of Mpro was found to differentially recognize viral peptidyl substrates. For instance, S3' in Mpro strongly favors Phe or Trp, and S1' favors Ala. The peptidyl inhibitor D-4-77, which possesses an α-bromoacetamide warhead, was discovered to be a promising inhibitor of Mpro , with an IC50 of 0.95 μM and an antiviral EC50 of 0.49 μM. The Mpro /inhibitor co-crystal structure confirmed the binding mode of the inhibitor to the S1'-S3' pocket and revealed a covalent mechanism. In addition, D-4-77 functions as an immune protectant and suppresses SARS-CoV-2 Mpro -induced antagonism of the host NF-κB innate immune response. These findings indicate that the S1'-S3' pocket of SARS-CoV-2 Mpro is druggable, and that inhibiting SARS-CoV-2 Mpro can simultaneously protect human innate immunity and inhibit virion assembly.

Molecular mechanisms of SARS-CoV-2 resistance to nirmatrelvir

Nirmatrelvir is a specific antiviral drug that targets the main protease (Mpro) of SARS-CoV-2 and has been approved to treat COVID-191,2. As an RNA virus characterized by high mutation rates, whether SARS-CoV-2 will develop resistance to nirmatrelvir is a question of concern. Our previous studies have shown that several mutational pathways confer resistance to nirmatrelvir, but some result in a loss of viral replicative fitness, which is then compensated for by additional alterations3. The molecular mechanisms for this observed resistance are unknown. Here we combined biochemical and structural methods to demonstrate that alterations at the substrate-binding pocket of Mpro can allow SARS-CoV-2 to develop resistance to nirmatrelvir in two distinct ways. Comprehensive studies of the structures of 14 Mpro mutants in complex with drugs or substrate revealed that alterations at the S1 and S4 subsites substantially decreased the level of inhibitor binding, whereas alterations at the S2 and S4' subsites unexpectedly increased protease activity. Both mechanisms contributed to nirmatrelvir resistance, with the latter compensating for the loss in enzymatic activity of the former, which in turn accounted for the restoration of viral replicative fitness, as observed previously3. Such a profile was also observed for ensitrelvir, another clinically relevant Mpro inhibitor. These results shed light on the mechanisms by which SARS-CoV-2 evolves to develop resistance to the current generation of protease inhibitors and provide the basis for the design of next-generation Mpro inhibitors.

The molecular mechanism of non-covalent inhibitor WU-04 targeting SARS-CoV-2 3CLpro and computational evaluation of its effectiveness against mainstream coronaviruses

There is an urgent need for highly effective therapeutic agents to interrupt the continued spread of SARS-CoV-2. As a pivotal protease in the replication process of coronaviruses, the 3CLpro protein is considered as a potential target of drug development to stop the spread and infection of the virus. In this work, molecular dynamics (MD) simulations were used to elucidate the molecular mechanism of a novel and highly effective non-covalent inhibitor, WU-04, targeting the SARS-CoV-2 3CLpro protein. The difference in dynamic behavior between the apo-3CLpro and the holo-3CLpro systems suggests that the presence of WU-04 inhibits the motion amplitude of the 3CLpro protein relative to the apo-3CLpro system, thus maintaining a stable conformational binding state. The energy calculations and interaction analysis show that the hot-spot residues Q189, M165, M49, E166, and H41 and the warm-spot residues H163 and C145 have a strong binding capacity to WU-04 by forming multiple hydrogen bonds and hydrophobic interactions, which stabilizes the binding of the inhibitor. After that, the resistance of WU-04 to the six SARS-CoV-2 variants (Alpha, Beta, Gamma, Delta, Lambda, and Omicron) and two other mainstream coronavirus (SARS-CoV and MERS-CoV) 3CLpro proteins was further investigated. Excitingly, the slight difference in energy values relative to the SARS-CoV-2 system indicates that WU-04 is still highly effective against the coronaviruses, which becomes crucial evidence that WU-04 is a pan-inhibitor of the 3CLpro protein in various SARS-CoV-2 variants and other mainstream coronaviruses. The study will hopefully provide theoretical insights for the future rational design and improvement of novel non-covalent inhibitors targeting the 3CLpro protein.

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.

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.

Crystal structures of main protease (Mpro) mutants of SARS-CoV-2 variants bound to PF-07304814

There is an urgent need to develop effective antiviral drugs to prevent the viral infection caused by constantly circulating SARS-CoV-2 as well as its variants. The main protease (Mpro) of SARS-CoV-2 is a salient enzyme that plays a vital role in viral replication and serves as a fascinating therapeutic target. PF-07304814 is a covalent inhibitor targeting SARS-CoV-2 Mpro with favorable inhibition potency and drug-like properties, thus making it a promising drug candidate for the treatment of COVID-19. We previously solved the structure of PF-07304814 in complex with SARS-CoV-2 Mpro. However, the binding modes of PF-07304814 with Mpros from evolving SARS-CoV-2 variants is under-determined. In the current study, we expressed six Mpro mutants (G15S, K90R, M49I, S46F, V186F, and Y54C) that have been identified in Omicron variants including the recently emerged XBB.1.16 subvariant and solved the crystal structures of PF-07304814 bound to Mpro mutants. Structural analysis provided insight into the key molecular determinants responsible for the interaction between PF-07304814 and these mutant Mpros. Patterns for PF-07304814 to bind with these investigated Mpro mutants and the wild-type Mpro are generally similar but with some differences as revealed by detailed structural comparison. Structural insights presented in this study will inform the development of novel drugs against SARS-CoV-2 and the possible conformation changes of Mpro mutants when bound to an inhibitor.

Discovery and evaluation of active compounds from Xuanfei Baidu formula against COVID-19 via SARS-CoV-2 Mpro

BACKGROUND

The coronavirus disease 2019 (COVID-19) caused by the severe acute respiratory syndrome coronavirus (SARS-CoV-2) is still a widespread concern. As one of the effective traditional Chinese medicine (TCM) formulae, Xuanfei Baidu formula (XFBD) shows significant efficacy for treatment of COVID-19 patients. However, its antiviral active compounds and mechanism are still unclear.

PURPOSE

In this study, we explored the bioactive compounds of XFBD and its antiviral mechanism by integrating computational analysis and experimental testing.

METHODS

Focusing on the SARS-CoV-2 main protease (M^pro), as a key target in virus transcription and replication, the fluorescence resonance energy transfer (FRET) assay was built to screen out satisfactory natural inhibitors in XFBD. The surface plasmon resonance (SPR) and isothermal titration calorimetry (ITC) were undertaken to verify the binding affinity of ligand-M^pro. Omicron BA.1.1 and BA.2.3 variants were used to evaluate the antiviral activity of the focused compounds in non-cytotoxicity concentrations. For introducing the molecular mechanism, computational modeling and NMR spectra were employed to characterize the ligand-binding modes and identify the ligand-binding site on M^pro.

RESULTS

From a library of 83 natural compounds, acteoside, licochalcone B, licochalcone D, linoleic acid, and physcion showed the satisfactory inhibition effects on M^pro with IC₅₀ ranging from 1.93 to 42.96 µM, which were further verified by SPR. Showing the excellent binding affinity, acteoside was witnessed to gain valuable insights into the thermodynamic signatures by ITC and presented antiviral activity on Omicron BA.1.1 and BA.2.3 variants in vitro. The results revealed that acteoside inhibited M^pro via forming the hydrogen bond between 7-H of acteoside and M^pro.

CONCLUSION

Acteoside is regarded as a representative active natural compound in XFBD to inhibit replication of SARS-CoV-2, which provides the antiviral evidence and some insights into the identification of SARS-CoV-2 M^pro natural inhibitors.

Evaluation of the Inhibition Potency of Nirmatrelvir against Main Protease Mutants of SARS-CoV-2 Variants

SARS-CoV-2 continues to pose a threat to public health. Main protease (Mpro) is one of the most lucrative drug targets for developing specific antivirals against SARS-CoV-2 infection. By targeting Mpro, peptidomimetic nirmatrelvir is able to inhibit viral replication of SARS-CoV-2 and reduce the risk for progression to severe COVID-19. However, multiple mutations in the gene encoding Mpro of emerging SARS-CoV-2 variants raise a concern of drug resistance. In the present study, we expressed 16 previously reported SARS-CoV-2 Mpro mutants (G15S, T25I, T45I, S46F, S46P, D48N, M49I, L50F, L89F, K90R, P132H, N142S, V186F, R188K, T190I, and A191V). We evaluated the inhibition potency of nirmatrelvir against these Mpro mutants and solved the crystal structures of representative Mpro mutants of SARS-CoV-2 bound to nirmatrelvir. Enzymatic inhibition assays revealed that these Mpro variants remain susceptible to nirmatrelvir as the wildtype. Detailed analysis and structural comparison provided the inhibition mechanism of Mpro mutants by nirmatrelvir. These results informed the ongoing genomic surveillance of drug resistance of emerging SARS-CoV-2 variants to nirmatrelvir and facilitate the development of next-generation anticoronavirus drugs.

Comprehensive Understanding of the Kinetic Behaviors of Main Protease from SARS-CoV-2 and SARS-CoV: New Data and Comparison to Published Parameters

The main protease (M^pro) is a promising drug target for inhibiting the coronavirus due to its conserved properties and lack of homologous genes in humans. However, previous studies on M^pro's kinetic parameters have been confusing, hindering the selection of accurate inhibitors. Therefore, obtaining a clear view of M^pro's kinetic parameters is necessary. In our study, we investigated the kinetic behaviors of M^pro from SARS-CoV-2 and SARS-CoV using both FRET-based cleavage assay and the LC-MS method, respectively. Our findings indicate that the FRET-based cleavage assay could be used for preliminary screening of M^pro inhibitors, while the LC-MS method should be applied to select the effective inhibitors with higher reliability. Furthermore, we constructed the active site mutants (H41A and C145A) and measured the kinetic parameters to gain a deeper understanding of the atomic-level enzyme efficiency reduction compared to the wild type. Overall, our study provides valuable insights for inhibitor screening and design by offering a comprehensive understanding of M^pro's kinetic behaviors.

No Substrate Left behind-Mining of Shotgun Proteomics Datasets Rescues Evidence of Proteolysis by SARS-CoV-2 3CLpro Main Protease

Proteolytic processing is the most ubiquitous post-translational modification and regulator of protein function. To identify protease substrates, and hence the function of proteases, terminomics workflows have been developed to enrich and detect proteolytically generated protein termini from mass spectrometry data. The mining of shotgun proteomics datasets for such 'neo'-termini, to increase the understanding of proteolytic processing, is an underutilized opportunity. However, to date, this approach has been hindered by the lack of software with sufficient speed to make searching for the relatively low numbers of protease-generated semi-tryptic peptides present in non-enriched samples viable. We reanalyzed published shotgun proteomics datasets for evidence of proteolytic processing in COVID-19 using the recently upgraded MSFragger/FragPipe software, which searches data with a speed that is an order of magnitude greater than many equivalent tools. The number of protein termini identified was higher than expected and constituted around half the number of termini detected by two different N-terminomics methods. We identified neo-N- and C-termini generated during SARS-CoV-2 infection that were indicative of proteolysis and were mediated by both viral and host proteases-a number of which had been recently validated by in vitro assays. Thus, re-analyzing existing shotgun proteomics data is a valuable adjunct for terminomics research that can be readily tapped (for example, in the next pandemic where data would be scarce) to increase the understanding of protease function and virus-host interactions, or other diverse biological processes.

Design of Potent Inhibitors Targeting the Main Protease of SARS-CoV-2 Using QSAR Modeling, Molecular Docking, and Molecular Dynamics Simulations

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection is a serious global public health threat. The evolving strains of SARS-CoV-2 have reduced the effectiveness of vaccines. Therefore, antiviral drugs against SARS-CoV-2 are urgently needed. The main protease (Mpro) of SARS-CoV-2 is an extremely potent target due to its pivotal role in virus replication and low susceptibility to mutation. In the present study, a quantitative structure-activity relationship (QSAR) study was performed to design new molecules that might have higher inhibitory activity against SARS-CoV-2 Mpro. In this context, a set of 55 dihydrophenanthrene derivatives was used to build two 2D-QSAR models using the Monte Carlo optimization method and the Genetic Algorithm Multi-Linear Regression (GA-MLR) method. From the CORAL QSAR model outputs, the promoters responsible for the increase/decrease in inhibitory activity were extracted and interpreted. The promoters responsible for an increase in activity were added to the lead compound to design new molecules. The GA-MLR QSAR model was used to ensure the inhibitory activity of the designed molecules. For further validation, the designed molecules were subjected to molecular docking analysis and molecular dynamics simulations along with an absorption, distribution, metabolism, excretion, and toxicity (ADMET) analysis. The results of this study suggest that the newly designed molecules have the potential to be developed as effective drugs against SARS-CoV-2.

Transmissible SARS-CoV-2 variants with resistance to clinical protease inhibitors

Vaccines and drugs have helped reduce disease severity and blunt the spread of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). However, ongoing virus transmission, continuous evolution, and increasing selective pressures have the potential to yield viral variants capable of resisting these interventions. Here, we investigate the susceptibility of natural variants of the main protease [Mpro; 3C-like protease (3CLpro)] of SARS-CoV-2 to protease inhibitors. Multiple single amino acid changes in Mpro confer resistance to nirmatrelvir (the active component of Paxlovid). An additional clinical-stage inhibitor, ensitrelvir (Xocova), shows a different resistance mutation profile. Importantly, phylogenetic analyses indicate that several of these resistant variants have pre-existed the introduction of these drugs into the human population and are capable of spreading. These results encourage the monitoring of resistance variants and the development of additional protease inhibitors and other antiviral drugs with different mechanisms of action and resistance profiles for combinatorial therapy.

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.

Optimization of the expression of the main protease from SARS-CoV-2

The main protease (M^pro) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) plays a vital role in viral replication. To study the function of M^pro and screen inhibitors targeting M^pro, it is necessary to prepare high-purity and high-activity M^pro. In this study, four types of SARS-CoV-2 M^pros containing different termini were prepared, and their activities were determined successfully. The results showed that the activity of wild-type (WT) M^pro was the highest, and the additional residues at the N-terminus but not at the C-terminus had a major effect on the enzyme activity. To explain this, the alignment of structures of different forms of M^pro was determined, and the additional residues at the N-terminus were found to interfere with the formation of the substrate binding pocket. This study confirms the importance of the natural N-terminus to the activity of M^pro and suggests that WT-GPH₆ (M^pro with eight additional residues at the C-terminus) can be used as a substitute for authentic M^pro to screen inhibitors. In short, this study provides a reference for the expression and purification of new coronaviruses confronted in the future.

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 comprehensive survey of coronaviral main protease active site diversity in 3D: Identifying and analyzing drug discovery targets in search of broad specificity inhibitors for the next coronavirus pandemic

Although the rapid development of therapeutic responses to combat SARS-CoV-2 represents a great human achievement, it also demonstrates untapped potential for advanced pandemic preparedness. Cross-species efficacy against multiple human coronaviruses by the main protease (MPro) inhibitor nirmatrelvir raises the question of its breadth of inhibition and our preparedness against future coronaviral threats. Herein, we describe sequence and structural analyses of 346 unique MPro enzymes from all coronaviruses represented in the NCBI Virus database. Cognate substrates of these representative proteases were inferred from their polyprotein sequences. We clustered MPro sequences based on sequence identity and AlphaFold2-predicted structures, showing approximate correspondence with known viral subspecies. Predicted structures of five representative MPros bound to their inferred cognate substrates showed high conservation in protease:substrate interaction modes, with some notable differences. Yeast-based proteolysis assays of the five representatives were able to confirm activity of three on inferred cognate substrates, and demonstrated that of the three, only one was effectively inhibited by nirmatrelvir. Our findings suggest that comprehensive preparedness against future potential coronaviral threats will require continued inhibitor development. Our methods may be applied to candidate coronaviral MPro inhibitors to evaluate in advance the breadth of their inhibition and identify target coronaviruses potentially meriting advanced development of alternative countermeasures.

Broad-Spectrum Cyclopropane-Based Inhibitors of Coronavirus 3C-like Proteases: Biochemical, Structural, and Virological Studies

The advent of SARS-CoV-2, the causative agent of COVID-19, and its worldwide impact on global health, have provided the impetus for the development of effective countermeasures that can be deployed against the virus, including vaccines, monoclonal antibodies, and direct-acting antivirals (DAAs). Despite these efforts, the current paucity of DAAs has created an urgent need for the creation of an enhanced and diversified portfolio of broadly acting agents with different mechanisms of action that can effectively abrogate viral infection. SARS-CoV-2 3C-like protease (3CLpro), an enzyme essential for viral replication, is a validated target for the discovery of SARS-CoV-2 therapeutics. In this report, we describe the structure-guided utilization of the cyclopropane moiety in the design of highly potent inhibitors of SARS-CoV-2 3CLpro, SARS-CoV-1 3CLpro, and MERS-CoV 3CLpro. High-resolution cocrystal structures were used to identify the structural determinants associated with the binding of the inhibitors to the active site of the enzyme and unravel the mechanism of action. Aldehydes 5c and 11c inhibited SARS-CoV-2 replication with EC50 values of 12 and 11 nM, respectively. Furthermore, the corresponding aldehyde bisulfite adducts 5d and 11d were equipotent with EC50 values of 13 and 12 nM, respectively. The safety index (SI) values for compounds 5c / 11c and 5d / 11d ranged between 7692 and 9090. Importantly, aldehydes 5c / 11c and bisulfite adducts 5d / 11d potently inhibited MERS-CoV 3CLpro with IC50 values of 80 and 120 nM, and 70 and 70 nM, respectively. Likewise, compounds 5c / 11c and 5d / 11d inhibited SARS-CoV-1 with IC50 values of 960 and 350 nM and 790 and 240 nM, respectively. Taken together, these studies suggest that the inhibitors described herein have low cytotoxicity and high potency and are promising candidates for further development as broad-spectrum direct-acting antivirals against highly pathogenic coronaviruses.

Binding and Unbinding Pathways of Peptide Substrates on the SARS-CoV-2 3CL Protease

Based on many crystal structures of ligand complexes, much study has been devoted to understanding the molecular recognition of SARS-CoV-2 3C-like protease (3CL^pro), a potent drug target for COVID-19. In this research, to extend this present static view, we examined the kinetic process of binding/unbinding of an eight-residue substrate peptide to/from 3CL^pro by evaluating the path ensemble with the weighted ensemble simulation. The path ensemble showed the mechanism of how a highly flexible peptide folded into the bound form. At the early stage, the dominant motion was the diffusion on the protein surface showing a broad distribution, whose center was led into the cleft of the chymotrypsin fold. We observed a definite sequential formation of the hydrogen bonds at the later stage occurring in the cleft, initiated between Glu166 (3CL^pro) and P3_Val (peptide), followed by binding to the oxyanion hole and completed by the sequence-specific recognition at P1_Gln.

Innovative, rapid, high-throughput method for drug repurposing in a pandemic-A case study of SARS-CoV-2 and COVID-19

Several efforts to repurpose drugs for COVID-19 treatment have largely either failed to identify a suitable agent or agents identified did not translate to clinical use. Reasons that have been suggested to explain the failures include use of inappropriate doses, that are not clinically achievable, in the screening experiments, and the use of inappropriate pre-clinical laboratory surrogates to predict efficacy. In this study, we used an innovative algorithm, that incorporates dissemination and implementation considerations, to identify potential drugs for COVID-19 using iterative computational and wet laboratory methods. The drugs were screened at doses that are known to be achievable in humans. Furthermore, inhibition of viral induced cytopathic effect (CPE) was used as the laboratory surrogate to predict efficacy. Erythromycin, pyridoxine, folic acid and retapamulin were found to inhibit SARS-CoV-2 induced CPE in Vero cells at concentrations that are clinically achievable. Additional studies may be required to further characterize the inhibitions of CPE and the possible mechanisms.

Drug development targeting SARS-CoV-2 main protease

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and its variants are responsible for the devastating coronavirus disease 2019 (COVID-19) pandemic with more than 6.5 million deaths since 2019. Although a number of vaccines significantly reduced the mortality rate, a large number of the world population is yet being infected with highly contagious omicron variants/subvarints. Additional therapeutic interventions are needed to reduce hospitalization and curb the ongoing pandemic. The activity of the SARS-CoV-2 enzyme; chymotrypsin-like main protease (M^pro) is essential for the cleavage of viral nonstructural polypeptides into individual functional proteins and therefore M^pro is an attractive drug target. The aim of this review is to summarize recent progress toward the development of therapeutic drugs against M^pro protease.

Discovery and Crystallographic Studies of Nonpeptidic Piperazine Derivatives as Covalent SARS-CoV-2 Main Protease Inhibitors

The spread of SARS-CoV-2 keeps threatening human life and health, and small-molecule antivirals are in demand. The main protease (Mpro) is an effective and highly conserved target for anti-SARS-CoV-2 drug design. Herein, we report the discovery of potent covalent non-peptide-derived Mpro inhibitors. A series of covalent compounds with a piperazine scaffold containing different warheads were designed and synthesized. Among them, GD-9 was identified as the most potent compound with a significant enzymatic inhibition of Mpro (IC50 = 0.18 μM) and good antiviral potency against SARS-CoV-2 (EC50 = 2.64 μM), similar to that of remdesivir (EC50 = 2.27 μM). Additionally, GD-9 presented favorable target selectivity for SARS-CoV-2 Mpro versus human cysteine proteases. The X-ray co-crystal structure confirmed our original design concept showing that GD-9 covalently binds to the active site of Mpro. Our nonpeptidic covalent inhibitors provide a basis for the future development of more efficient COVID-19 therapeutics.

Biochemical and structural insights into SARS-CoV-2 polyprotein processing by Mpro

SARS-CoV-2, a human coronavirus, is the causative agent of the COVID-19 pandemic. Its genome is translated into two large polyproteins subsequently cleaved by viral papain-like protease and main protease (Mpro). Polyprotein processing is essential yet incompletely understood. We studied Mpro-mediated processing of the nsp7-11 polyprotein, whose mature products include cofactors of the viral replicase, and identified the order of cleavages. Integrative modeling based on mass spectrometry (including hydrogen-deuterium exchange and cross-linking) and x-ray scattering yielded a nsp7-11 structural ensemble, demonstrating shared secondary structural elements with individual nsps. The pattern of cross-links and HDX footprint of the C145A Mpro and nsp7-11 complex demonstrate preferential binding of the enzyme active site to the polyprotein junction sites and additional transient contacts to help orient the enzyme on its substrate for cleavage. Last, proteolysis assays were used to characterize the effect of inhibitors/binders on Mpro processing/inhibition using the nsp7-11 polyprotein as substrate.

Progress on COVID-19 Chemotherapeutics Discovery and Novel Technology

COVID-19 is an infectious disease caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), a novel highly contagious and pathogenic coronavirus that emerged in late 2019. SARS-CoV-2 spreads primarily through virus-containing droplets and small particles of air pollution, which greatly increases the risk of inhaling these virus particles when people are in close proximity. COVID-19 is spreading across the world, and the COVID-19 pandemic poses a threat to human health and public safety. To date, there are no specific vaccines or effective drugs against SARS-CoV-2. In this review, we focus on the enzyme targets of the virus and host that may be critical for the discovery of chemical compounds and natural products as antiviral drugs, and describe the development of potential antiviral drugs in the preclinical and clinical stages. At the same time, we summarize novel emerging technologies applied to the research on new drug development and the pathological mechanisms of COVID-19.