Structure-guided design of a perampanel-derived pharmacophore targeting the SARS-CoV-2 main protease

Investigating the Inhibition of Diindolylmethane Derivatives on SARS-CoV-2 Main Protease

The SARS-CoV-2 main protease (Mpro) is an essential enzyme that promotes viral transcription and replication. Mpro conserved nature in different variants and its nonoverlapping nature with human proteases make it an attractive target for therapeutic intervention against SARS-CoV-2. In this work, the interaction mechanism between Mpro and diindolylmethane derivatives was investigated by molecular docking, enzymatic inhibition assay, UV-vis, fluorescence spectroscopy, and circular dichroism spectroscopy. Results of IC50 values show that 1p (9.87 μM) was the strongest inhibitor for Mpro in this work, which significantly inhibited the activity of Mpro. The binding constant (4.07 × 105 Lmol-1), the quenching constant (5.41 × 105 Lmol-1), and thermodynamic parameters indicated that the quenching mode of 1p was static quenching, and the main driving forces between 1p and Mpro are hydrogen bond and van der Waals force. The influence of molecular structure on the binding is investigated. Chlorine atoms and methoxy groups are favorable for the diindolylmethane derivative inhibitors of Mpro. This work confirms the changes in the microenvironment of Mpro by 1p, and provides clues for the design of potential inhibitors.

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.

Leak Proof PDBBind: A Reorganized Dataset of Protein-Ligand Complexes for More Generalizable Binding Affinity Prediction

Many physics-based and machine-learned scoring functions (SFs) used to predict protein-ligand binding free energies have been trained on the PDBBind dataset. However, it is controversial as to whether new SFs are actually improving since the general, refined, and core datasets of PDBBind are cross-contaminated with proteins and ligands with high similarity, and hence they may not perform comparably well in binding prediction of new protein-ligand complexes. In this work we have carefully prepared a cleaned PDBBind data set of non-covalent binders that are split into training, validation, and test datasets to control for data leakage, defined as proteins and ligands with high sequence and structural similarity. The resulting leak-proof (LP)-PDBBind data is used to retrain four popular SFs: AutoDock Vina, Random Forest (RF)-Score, InteractionGraphNet (IGN), and DeepDTA, to better test their capabilities when applied to new protein-ligand complexes. In particular we have formulated a new independent data set, BDB2020+, by matching high quality binding free energies from BindingDB with co-crystalized ligand-protein complexes from the PDB that have been deposited since 2020. Based on all the benchmark results, the retrained models using LP-PDBBind consistently perform better, with IGN especially being recommended for scoring and ranking applications for new protein-ligand systems.

Proof-of-concept studies with a computationally designed Mpro inhibitor as a synergistic combination regimen alternative to Paxlovid

As the SARS-CoV-2 virus continues to spread and mutate, it remains important to focus not only on preventing spread through vaccination but also on treating infection with direct-acting antivirals (DAA). The approval of Paxlovid, a SARS-CoV-2 main protease (Mpro) DAA, has been significant for treatment of patients. A limitation of this DAA, however, is that the antiviral component, nirmatrelvir, is rapidly metabolized and requires inclusion of a CYP450 3A4 metabolic inhibitor, ritonavir, to boost levels of the active drug. Serious drug-drug interactions can occur with Paxlovid for patients who are also taking other medications metabolized by CYP4503A4, particularly transplant or otherwise immunocompromised patients who are most at risk for SARS-CoV-2 infection and the development of severe symptoms. Developing an alternative antiviral with improved pharmacological properties is critical for treatment of these patients. By using a computational and structure-guided approach, we were able to optimize a 100 to 250 μM screening hit to a potent nanomolar inhibitor and lead compound, Mpro61. In this study, we further evaluate Mpro61 as a lead compound, starting with examination of its mode of binding to SARS-CoV-2 Mpro. In vitro pharmacological profiling established a lack of off-target effects, particularly CYP450 3A4 inhibition, as well as potential for synergy with the currently approved alternate antiviral, molnupiravir. Development and subsequent testing of a capsule formulation for oral dosing of Mpro61 in B6-K18-hACE2 mice demonstrated favorable pharmacological properties, efficacy, and synergy with molnupiravir, and complete recovery from subsequent challenge by SARS-CoV-2, establishing Mpro61 as a promising potential preclinical candidate.

On the origins of SARS-CoV-2 main protease inhibitors

In order to address the world-wide health challenge caused by the COVID-19 pandemic, the 3CL protease/SARS-CoV-2 main protease (SARS-CoV-2-Mpro) coded by its nsp5 gene became one of the biochemical targets for the design of antiviral drugs. In less than 3 years of research, 4 inhibitors of SARS-CoV-2-Mpro have actually been authorized for COVID-19 treatment (nirmatrelvir, ensitrelvir, leritrelvir and simnotrelvir) and more such as EDP-235, FB-2001 and STI-1558/Olgotrelvir or five undisclosed compounds (CDI-988, ASC11, ALG-097558, QLS1128 and H-10517) are undergoing clinical trials. This review is an attempt to picture this quite unprecedented medicinal chemistry feat and provide insights on how these cysteine protease inhibitors were discovered. Since many series of covalent SARS-CoV-2-Mpro inhibitors owe some of their origins to previous work on other proteases, we first provided a description of various inhibitors of cysteine-bearing human caspase-1 or cathepsin K, as well as inhibitors of serine proteases such as human dipeptidyl peptidase-4 or the hepatitis C protein complex NS3/4A. This is then followed by a description of the results of the approaches adopted (repurposing, structure-based and high throughput screening) to discover coronavirus main protease inhibitors.

Development of de-novo coronavirus 3-chymotrypsin-like protease (3CLpro) inhibitors since COVID-19 outbreak: A strategy to tackle challenges of persistent virus infection

Although no longer a public health emergency of international concern, COVID-19 remains a persistent and critical health concern. The development of effective antiviral drugs could serve as the ultimate piece of the puzzle to curbing this global crisis. 3-chymotrypsin-like protease (3CLpro), with its substrate specificity mirroring that of the main picornavirus 3C protease and conserved across various coronaviruses, emerges as an ideal candidate for broad-spectrum antiviral drug development. Moreover, it holds the potential as a reliable contingency option to combat emerging SARS-CoV-2 variants. In this light, the approved drugs, promising candidates, and de-novo small molecule therapeutics targeting 3CLpro since the COVID-19 outbreak in 2020 are discussed. Emphasizing the significance of diverse structural characteristics in inhibitors, be they peptidomimetic or nonpeptidic, with a shared mission to minimize the risk of cross-resistance. Moreover, the authors propose an innovative optimization strategy for 3CLpro reversible covalent PROTACs, optimizing pharmacodynamics and pharmacokinetics to better prepare for potential future viral outbreaks.

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.

Structure and function of SARS-CoV and SARS-CoV-2 main proteases and their inhibition: A comprehensive review

Severe acute respiratory syndrome-associated coronavirus (SARS-CoV) identified in 2003 infected ∼8000 people in 26 countries with 800 deaths, which was soon contained and eradicated by syndromic surveillance and enhanced quarantine. A closely related coronavirus SARS-CoV-2, the causative agent of COVID-19 identified in 2019, has been dramatically more contagious and catastrophic. It has infected and caused various flu-like symptoms of billions of people in >200 countries, including >6 million people died of or with the virus. Despite the availability of several vaccines and antiviral drugs against SARS-CoV-2, finding new therapeutics is needed because of viral evolution and a possible emerging coronavirus in the future. The main protease (Mpro) of these coronaviruses plays important roles in their life cycle and is essential for the viral replication. This article represents a comprehensive review of the function, structure and inhibition of SARS-CoV and -CoV-2 Mpro, including structure-activity relationships, protein-inhibitor interactions and clinical trial status.

Identification of novel 1,2,3-triazole isatin derivatives as potent SARS-CoV-2 3CLpro inhibitors via click-chemistry-based rapid screening

SARS-CoV-2 3-chymotrypsin-like protease (3CLpro) is considered an attractive target for the development of anti-COVID-19 agents due to its vital function. The N-substituted isatin derivative L-26 is a potential SARS-CoV-2 3CLpro inhibitor, but it has poor cell-based antiviral activity and high cytotoxicity. With L-26 as the lead compound, 58 isatin derivatives were prepared using click-chemistry-based miniaturized synthesis and their 3CLpro inhibitory activities were determined by a fluorescence resonance energy transfer-based enzymatic assay. Compounds D1N8 (IC50 = 0.44 ± 0.12 μM) and D1N52 (IC50 = 0.53 ± 0.21 μM) displayed excellent inhibitory potency against SARS-CoV-2 3CLpro, being equivalent to that of L-26 (IC50 = 0.30 ± 0.14 μM). In addition, the cytotoxicity of D1N8 (CC50 >20 μM) and D1N52 (CC50 >20 μM) decreased significantly compared with L-26 (CC50 <2.6 μM). Further molecular dynamics simulations revealed the potential binding interactions between D1N52 and SARS-CoV-2 3CLpro. These efforts lay a solid foundation for the research of novel anti-SARS-CoV-2 agents targeting 3CLpro.

Design and Optimization of Novel Competitive, Non-peptidic, SARS-CoV-2 Mpro Inhibitors

The SARS-CoV-2 main protease (Mpro) has been proven to be a highly effective target for therapeutic intervention, yet only one drug currently holds FDA approval status for this target. We were inspired by a series of publications emanating from the Jorgensen and Anderson groups describing the design of potent, non-peptidic, competitive SARS-CoV-2 Mpro inhibitors, and we saw an opportunity to make several design modifications to improve the overall pharmacokinetic profile of these compounds without losing potency. To this end, we created a focused virtual library using reaction-based enumeration tools in the Schrödinger suite. These compounds were docked into the Mpro active site and subsequently prioritized for synthesis based upon relative binding affinity values calculated by FEP+. Fourteen compounds were selected, synthesized, and evaluated both biochemically and in cell culture. Several of the synthesized compounds proved to be potent, competitive Mpro inhibitors with improved metabolic stability profiles.

Computation of Absolute Binding Free Energies for Noncovalent Inhibitors with SARS-CoV-2 Main Protease

Accurate, routine calculation of absolute binding free energies (ABFEs) for protein-ligand complexes remains a key goal of computer-aided drug design since it can enable screening and optimization of drug candidates. For development and testing of related methods, it is important to have high-quality datasets. To this end, from our own experimental studies, we have selected a set of 16 inhibitors of the SARS-CoV-2 main protease (Mpro) with structural diversity and well-distributed BFEs covering a 5 kcal/mol range. There is also minimal structural uncertainty since X-ray crystal structures have been deposited for 12 of the compounds. For methods testing, we report ABFE results from 2 μs molecular dynamics (MD) simulations using free energy perturbation (FEP) theory. The correlation of experimental and computed results is encouraging, with a Pearson's r2 of 0.58 and a Kendall τ of 0.24. The results indicate that current FEP-based ABFE calculations can be used for identification of active compounds (hits). While their accuracy for lead optimization is not yet sufficient, this activity remains addressable in separate lead series by relative BFE calculations.

Drug Potency Prediction of SARS-CoV-2 Main Protease Inhibitors Based on a Graph Generative Model

The prediction of a ligand potency to inhibit SARS-CoV-2 main protease (M-pro) would be a highly helpful addition to a virtual screening process. The most potent compounds might then be the focus of further efforts to experimentally validate their potency and improve them. A computational method to predict drug potency, which is based on three main steps, is defined: (1) defining the drug and protein in only one 3D structure; (2) applying graph autoencoder techniques with the aim of generating a latent vector; and (3) using a classical fitting model to the latent vector to predict the potency of the drug. Experiments in a database of 160 drug-M-pro pairs, from which the pIC50 is known, show the ability of our method to predict their drug potency with high accuracy. Moreover, the time spent to compute the pIC50 of the whole database is only some seconds, using a current personal computer. Thus, it can be concluded that a computational tool that predicts, with high reliability, the pIC50 in a cheap and fast way is achieved. This tool, which can be used to prioritize which virtual screening hits, will be further examined in vitro.

SARS-CoV-2 3CL-protease inhibitors derived from ML300: investigation of P1 and replacements of the 1,2,3-benzotriazole

Starting from compound 5 (CCF0058981), a structure-based optimization of the P1 subsite was performed against the severe acute respiratory syndrome coronavirus (SARS-CoV-2) main protease (3CLpro). Inhibitor 5 and the compounds disclosed bind to 3CLpro using a non-covalent mode of action that utilize a His163 H-bond interaction in the S1 subpocket. In an effort to examine more structurally diverse P1 groups a number of azoles and heterocycles were designed. Several azole ring systems and replacements, including C-linked azoles, with similar or enhanced potency relative to 5 were discovered (28, 29, and 30) with demonstrated IC50 values less than 100 nM. In addition, pyridyl and isoquinoline P1 groups were successful as P1 replacements leading to 3-methyl pyridyl 36 (IC50 = 85 nM) and isoquinoline 27 (IC50 = 26 nM). High resolution X-ray crystal structures of these inhibitors were utilized to confirm binding orientation and guide optimization. These findings have implications towards antiviral development and preparedness to combat SARS-like zoonotic coronavirus outbreaks.

Structure-guided design of direct-acting antivirals that exploit the gem-dimethyl effect and potently inhibit 3CL proteases of severe acute respiratory syndrome Coronavirus-2 (SARS-CoV-2) and middle east respiratory syndrome coronavirus (MERS-CoV)

The high morbidity and mortality associated with SARS-CoV-2 infection, the etiological agent of COVID-19, has had a major impact on global public health. Significant progress has been made in the development of an array of vaccines and biologics, however, the emergence of SARS-CoV-2 variants and breakthrough infections are an ongoing major concern. Furthermore, there is an existing paucity of small-molecule host and virus-directed therapeutics and prophylactics that can be used to counter the spread of SARS-CoV-2, and any emerging and re-emerging coronaviruses. We describe herein our efforts to address this urgent need by focusing on the structure-guided design of potent broad-spectrum inhibitors of SARS-CoV-2 3C-like protease (3CL^pro or Main protease), an enzyme essential for viral replication. The inhibitors exploit the directional effects associated with the presence of a gem-dimethyl group that allow the inhibitors to optimally interact with the S4 subsite of the enzyme. Several compounds were found to potently inhibit SARS-CoV-2 and MERS-CoV 3CL proteases in biochemical and cell-based assays. Specifically, the EC50 values of aldehyde 1c and its corresponding bisulfite adduct 1d against SARS-CoV-2 were found to be 12 and 10 nM, respectively, and their CC50 values were >50 μM. Furthermore, deuteration of these compounds yielded compounds 2c/2d with EC50 values 11 and 12 nM, respectively. Replacement of the aldehyde warhead with a nitrile (CN) or an α-ketoamide warhead or its corresponding bisulfite adduct yielded compounds 1g, 1eand1f with EC50 values 60, 50 and 70 nM, respectively. High-resolution cocrystal structures have identified the structural determinants associated with the binding of the inhibitors to the active site of the enzyme and, furthermore, have illuminated the mechanism of action of the inhibitors. Overall, the high Safety Index (SI) (SI=CC50/EC50) displayed by these compounds suggests that they are well-suited to conducting further preclinical studies.

Thermal Titration Molecular Dynamics (TTMD): Not Your Usual Post-Docking Refinement

Molecular docking is one of the most widely used computational approaches in the field of rational drug design, thanks to its favorable balance between the rapidity of execution and the accuracy of provided results. Although very efficient in exploring the conformational degrees of freedom available to the ligand, docking programs can sometimes suffer from inaccurate scoring and ranking of generated poses. To address this issue, several post-docking filters and refinement protocols have been proposed throughout the years, including pharmacophore models and molecular dynamics simulations. In this work, we present the first application of Thermal Titration Molecular Dynamics (TTMD), a recently developed method for the qualitative estimation of protein-ligand unbinding kinetics, to the refinement of docking results. TTMD evaluates the conservation of the native binding mode throughout a series of molecular dynamics simulations performed at progressively increasing temperatures with a scoring function based on protein-ligand interaction fingerprints. The protocol was successfully applied to retrieve the native-like binding pose among a set of decoy poses of drug-like ligands generated on four different pharmaceutically relevant biological targets, including casein kinase 1δ, casein kinase 2, pyruvate dehydrogenase kinase 2, and SARS-CoV-2 main protease.

Identification of new dihydrophenanthrene derivatives as promising anti-SARS-CoV-2 drugs through in silico investigations

To research, evaluate, and invent novel compounds that inhibit SARS-CoV-2 activity, a series of reported 39 substituted 9, 10-dihydrophenanthrene derivatives were subjected to a quantitative structure-activity relationship (QSAR) study. Gaussian 09 and ChemOffice programs were used to calculate the molecular descriptors employed to determine their impact on the studied activity. Then we reduced the number of descriptors by eliminating the redundant information using principal component analysis (PCA). The creation of molecular models was done by using multiple linear regression (MLR) according to the principles established by the Organization for Economic Co-operation and Development (OECD) and the validation by using external and internal validation, Y-randomization tests, and domain of applicability. Moreover, we evaluated the toxicity of developed compounds using ADMET and Molecular docking to determine their optimal position to form a stable complex. As a result, four molecules may be used to develop a novel drug that can inhibit SARS-CoV-2 without causing the side effect.

Phytonutrient Inhibitors of SARS-CoV-2/NSP5-Encoded Main Protease (Mpro) Autocleavage Enzyme Critical for COVID-19 Pathogenesis

The genomic reshuffling, mutagenicity, and high transmission rate of the SARS-CoV-2 pathogen highlights an urgent need for effective antiviral interventions for COVID-19 control. Targeting the highly conserved viral genes and/or gene-encoded viral proteins such as main proteinase (M^pro), RNA-dependent RNA polymerase (RdRp) and helicases are plausible antiviral approaches to prevent replication and propagation of the SARS-CoV-2 infection. Coronaviruses (CoVs) are prone to extensive mutagenesis; however, any genetic alteration to its highly conserved M^pro enzyme is often detrimental to the viral pathogen. Therefore, inhibitors that target the M^pro enzyme could reduce the risk of mutation-mediated drug resistance and provide effective antiviral protection. Several existing antiviral drugs and dietary bioactives are currently repurposed to treat COVID-19. Dietary bioactives from three ayurvedic medicinal herbs, 18 β-glycyrrhetinic acid (ΔG = 8.86 kcal/mol), Solanocapsine (ΔG = 8.59 kcal/mol), and Vasicoline (ΔG = 7.34 kcal/mol), showed high-affinity binding to M^pro enzyme than the native N3 inhibitor (ΔG = 5.41 kcal/mol). Flavonoids strongly inhibited SARS-CoV-2 M^pro with comparable or higher potency than the antiviral drug, remdesivir. Several tannin hydrolysates avidly bound to the receptor-binding domain and catalytic dyad (His₄₁ and Cys₁₄₅) of SARS-CoV-2 M^pro through H-bonding forces. Quercetin binding to M^pro altered the thermostability of the viral protein through redox-based mechanism and inhibited the viral enzymatic activity. Interaction of quercetin-derivatives with the M^pro seem to be influenced by the 7-OH group and the acetoxylation of sugar moiety on the ligand molecule. Based on pharmacokinetic and ADMET profiles, several phytonutrients could serve as a promising redox nutraceutical for COVID-19 management.

[The prediction of SARS-CoV-2 main protease inhibition with filtering by position of ligand]

The paper analyzes a set of equations that adequately predict the IC50 value for SARS-CoV-2 main protease inhibitors. The training set was obtained using filtering by criteria independent of prediction of target value. It included 76 compounds, and the test set included nine compounds. We used the values of energy contributions obtained in the calculation of the change of the free energy of complex by MMGBSA method and a number of characteristics of the physical and chemical properties of the inhibitors as independent variables. It is sufficient to use only seven independent variables without loss of prediction quality (Q² = 0.79; R²prediction = 0.89). The maximum error in this case does not exceed 0.92 lg(IC50) units with a full range of observed values from 1.26 to 4.95.

Qualitative Estimation of Protein-Ligand Complex Stability through Thermal Titration Molecular Dynamics Simulations

The prediction of ligand efficacy has long been linked to thermodynamic properties such as the equilibrium dissociation constant, which considers both the association and the dissociation rates of a defined protein-ligand complex. In the last 15 years, there has been a paradigm shift, with an increased interest in the determination of kinetic properties such as the drug-target residence time since they better correlate with ligand efficacy compared to other parameters. In this article, we present thermal titration molecular dynamics (TTMD), an alternative computational method that combines a series of molecular dynamics simulations performed at progressively increasing temperatures with a scoring function based on protein-ligand interaction fingerprints for the qualitative estimation of protein-ligand-binding stability. The protocol has been applied to four different pharmaceutically relevant test cases, including protein kinase CK1δ, protein kinase CK2, pyruvate dehydrogenase kinase 2, and SARS-CoV-2 main protease, on a variety of ligands with different sizes, structures, and experimentally determined affinity values. In all four cases, TTMD was successfully able to distinguish between high-affinity compounds (low nanomolar range) and low-affinity ones (micromolar), proving to be a useful screening tool for the prioritization of compounds in a drug discovery campaign.

Non-covalent cyclic peptides simultaneously targeting Mpro and NRP1 are highly effective against Omicron BA.2.75

Available vaccine-based immunity may at high risk of being evaded due to substantial mutations in the variant Omicron. The main protease (Mpro) of SARS-CoV-2 and human neuropilin-1 (NRP1), two less mutable proteins, have been reported to be crucial for SARS-CoV-2 replication and entry into host cells, respectively. Their dual blockade may avoid vaccine failure caused by continuous mutations of the SARS-CoV-2 genome and exert synergistic antiviral efficacy. Herein, four cyclic peptides non-covalently targeting both Mpro and NRP1 were identified using virtual screening. Among them, MN-2 showed highly potent affinity to Mpro (K_d = 18.2 ± 1.9 nM) and NRP1 (K_d = 12.3 ± 1.2 nM), which was about 3,478-fold and 74-fold stronger than that of the positive inhibitors Peptide-21 and EG3287. Furthermore, MN-2 exhibited significant inhibitory activity against Mpro and remarkable anti-infective activity against the pseudotyped variant Omicron BA.2.75 without obvious cytotoxicity. These data demonstrated that MN-2, a novel non-covalent cyclic peptide, is a promising agent against Omicron BA.2.75.

Discovery and Crystallographic Studies of Trisubstituted Piperazine Derivatives as Non-Covalent SARS-CoV-2 Main Protease Inhibitors with High Target Specificity and Low Toxicity

The continuous spread of SARS-CoV-2 calls for more direct-acting antiviral agents to combat the highly infectious variants. The main protease (Mpro) is an promising target for anti-SARS-CoV-2 drug design. Here, we report the discovery of potent non-covalent non-peptide Mpro inhibitors featuring a 1,2,4-trisubstituted piperazine scaffold. We systematically modified the non-covalent hit MCULE-5948770040 by structure-based rational design combined with multi-site binding and privileged structure assembly strategies. The optimized compound GC-14 inhibits Mpro with high potency (IC50 = 0.40 μM) and displays excellent antiviral activity (EC50 = 1.1 μM), being more potent than Remdesivir. Notably, GC-14 exhibits low cytotoxicity (CC50 > 100 μM) and excellent target selectivity for SARS-CoV-2 Mpro (IC50 > 50 μM for cathepsins B, F, K, L, and caspase 3). X-ray co-crystal structures prove that the inhibitors occupy multiple subpockets by critical non-covalent interactions. These studies may provide a basis for developing a more efficient and safer therapy for COVID-19.

X-ray crystallographic characterization of the SARS-CoV-2 main protease polyprotein cleavage sites essential for viral processing and maturation

Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2), the pathogen that causes COVID-19, produces polyproteins 1a and 1ab that contain, respectively, 11 or 16 non-structural proteins (nsp). Nsp5 is the main protease (M^pro) responsible for cleavage at eleven positions along these polyproteins, including at its own N- and C-terminal boundaries, representing essential processing events for viral assembly and maturation. Using C-terminally substituted M^pro chimeras, we have determined X-ray crystallographic structures of M^pro in complex with 10 of its 11 viral cleavage sites, bound at full occupancy intermolecularly in trans, within the active site of either the native enzyme and/or a catalytic mutant (C145A). Capture of both acyl-enzyme intermediate and product-like complex forms of a P2(Leu) substrate in the native active site provides direct comparative characterization of these mechanistic steps as well as further informs the basis for enhanced product release of M^pro’s own unique C-terminal P2(Phe) cleavage site to prevent autoinhibition. We characterize the underlying noncovalent interactions governing binding and specificity for this diverse set of substrates, showing remarkable plasticity for subsites beyond the anchoring P1(Gln)-P2(Leu/Val/Phe), representing together a near complete analysis of a multiprocessing viral protease. Collectively, these crystallographic snapshots provide valuable mechanistic and structural insights for antiviral therapeutic development.

The SARS-CoV-2 protease Mpro is essential for viral replication. Here, the authors have determined the structures of Mpro in complex with 10 of the 11 viral cleavage sequences including a covalent acyl-enzyme intermediate, providing mechanistic and structural insights for antiviral development.

The SARS-CoV-2 main protease (Mpro): Structure, function, and emerging therapies for COVID-19

The main proteases (Mpro), also termed 3-chymotrypsin-like proteases (3CLpro), are a class of highly conserved cysteine hydrolases in β-coronaviruses. Increasing evidence has demonstrated that 3CLpros play an indispensable role in viral replication and have been recognized as key targets for preventing and treating coronavirus-caused infectious diseases, including COVID-19. This review is focused on the structural features and biological function of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) main protease Mpro (also known as 3CLpro), as well as recent advances in discovering and developing SARS-CoV-2 3CLpro inhibitors. To better understand the characteristics of SARS-CoV-2 3CLpro inhibitors, the inhibition activities, inhibitory mechanisms, and key structural features of various 3CLpro inhibitors (including marketed drugs, peptidomimetic, and non-peptidomimetic synthetic compounds, as well as natural compounds and their derivatives) are summarized comprehensively. Meanwhile, the challenges in this field are highlighted, while future directions for designing and developing efficacious 3CLpro inhibitors as novel anti-coronavirus therapies are also proposed. Collectively, all information and knowledge presented here are very helpful for understanding the structural features and inhibitory mechanisms of SARS-CoV-2 3CLpro inhibitors, which offers new insights or inspiration to medicinal chemists for designing and developing more efficacious 3CLpro inhibitors as novel anti-coronavirus agents.

Rational identification of small molecules derived from 9,10-dihydrophenanthrene as potential inhibitors of 3CLpro enzyme for COVID-19 therapy: a computer-aided drug design approach

Small molecules such as 9,10-dihydrophenanthrene derivatives have remarkable activity toward inhibition of SARS-CoV-2 3CLpro and COVID-19 proliferation, which show a strong correlation between their structures and bioactivity. Therefore, these small compounds could be suitable for clinical pharmaceutical use against COVID-19. The objective of this study was to remodel the structures of 9,10-dihydrophenanthrene derivatives to achieve a powerful biological activity against 3CLpro and favorable pharmacokinetic properties for drug design and discovery. Therefore, by the use of bioinformatics techniques, we developed robust 3D-QSAR models that are capable of describing the structure-activity relationship for 46 molecules based on 9,10-dihydrophenanthrene derivatives using CoMFA/SE (R 2 = 0.97, Q 2 = 0.81, R 2 pred = 0.95, c R 2 p = 0.71) and CoMSIA/SEHDA (R 2 = 0.94, Q 2 = 0.76, R 2 pred = 0.91, c R 2 p = 0.65) techniques. Accordingly, 96 lead compounds were generated based on a template molecule that showed the highest observed activity in vitro (T40, pIC50 = 5.81) and predicted their activities and bioavailability in silico. The rational screening outputs of 3D-QSAR, Molecular docking, ADMET, and MM-GBSA led to the identification of 9 novel modeled molecules as potent noncovalent drugs against SARS-CoV-2-3CLpro. Finally, by molecular dynamics simulations, the stability and structural dynamics of 3CLpro free and complex (PDB code: 6LU7) were discussed in the presence of samples of 9,10-dihydrophenanthrene derivative in an aqueous environment. Overall, the retrosynthesis of the proposed drug compounds in this study and the evaluation of their bioactivity in vitro and in vivo may be interesting for designing and discovering a new drug effective against COVID-19.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11224-022-02004-z.

Ligand-based and structure-based studies to develop predictive models for SARS-CoV-2 main protease inhibitors through the 3d-qsar.com portal

The main protease (Mpro) of SARS-Cov-2 is the essential enzyme for maturation of functional proteins implicated in viral replication and transcription. The peculiarity of its specific cleavage site joint with its high degree of conservation among all coronaviruses promote it as an attractive target to develop broad-spectrum inhibitors, with high selectivity and tolerable safety profile. Herein is reported a combination of three-dimensional quantitative structure-activity relationships (3-D QSAR) and comparative molecular binding energy (COMBINE) analysis to build robust and predictive ligand-based and structure-based statistical models, respectively. Models were trained on experimental binding poses of co-crystallized Mpro-inhibitors and validated on available literature data. By means of deep optimization both models' goodness and robustness reached final statistical values of r2/q2 values of 0.97/0.79 and 0.93/0.79 for the 3-D QSAR and COMBINE approaches respectively, and an overall predictiveness values of 0.68 and 0.57 for the SDEPPRED and AAEP metrics after application to a test set of 60 compounds covered by the training set applicability domain. Despite the different nature (ligand-based and structure-based) of the employed methods, their outcome fully converged. Furthermore, joint ligand- and structure-based structure-activity relationships were found in good agreement with nirmatrelvir chemical features properties, a novel oral Mpro-inhibitor that has recently received U.S. FDA emergency use authorization (EUA) for the oral treatment of mild-to-moderate COVID-19 infected patients. The obtained results will guide future rational design and/or virtual screening campaigns with the aim of discovering new potential anti-coronavirus lead candidates, minimizing both time and financial resources. Moreover, as most of calculation were performed through the well-established web portal 3d-qsar.com the results confirm the portal as a useful tool for drug design.

Quercetin and luteolin are single-digit micromolar inhibitors of the SARS-CoV-2 RNA-dependent RNA polymerase

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has rapidly become a global health pandemic. Among the viral proteins, RNA-dependent RNA polymerase (RdRp) is responsible for viral genome replication and has emerged as one of the most promising targets for pharmacological intervention against SARS-CoV-2. To this end, we experimentally tested luteolin and quercetin for their ability to inhibit the RdRp enzyme. These two compounds are ancestors of flavonoid natural compounds known for a variety of basal pharmacological activities. Luteolin and quercetin returned a single-digit IC50 of 4.6 µM and 6.9 µM, respectively. Then, through dynamic docking simulations, we identified possible binding modes of these compounds to a recently published cryo-EM structure of RdRp. Collectively, these data indicate that these two compounds are a valid starting point for further optimization and development of a new class of RdRp inhibitors to treat SARS-CoV-2 and potentially other viral infections.

Structure-Guided Design of Potent Spirocyclic Inhibitors of Severe Acute Respiratory Syndrome Coronavirus-2 3C-like Protease

The worldwide impact of the ongoing COVID-19 pandemic on public health has made imperative the discovery and development of direct-acting antivirals aimed at targeting viral and/or host targets. SARS-CoV-2 3C-like protease (3CLpro) has emerged as a validated target for the discovery of SARS-CoV-2 therapeutics because of the pivotal role it plays in viral replication. We describe herein the structure-guided design of highly potent inhibitors of SARS-CoV-2 3CLpro that incorporate in their structure novel spirocyclic design elements aimed at optimizing potency by accessing new chemical space. Inhibitors of both SARS-CoV-2 3CLpro and MERS-CoV 3CLpro that exhibit nM potency and high safety indices have been identified. The mechanism of action of the inhibitors and the structural determinants associated with binding were established using high-resolution cocrystal structures.

Penicillin Derivatives Inhibit the SARS-CoV-2 Main Protease by Reaction with Its Nucleophilic Cysteine

The SARS-CoV-2 main protease (Mpro) is a medicinal chemistry target for COVID-19 treatment. Given the clinical efficacy of β-lactams as inhibitors of bacterial nucleophilic enzymes, they are of interest as inhibitors of viral nucleophilic serine and cysteine proteases. We describe the synthesis of penicillin derivatives which are potent Mpro inhibitors and investigate their mechanism of inhibition using mass spectrometric and crystallographic analyses. The results suggest that β-lactams have considerable potential as Mpro inhibitors via a mechanism involving reaction with the nucleophilic cysteine to form a stable acyl-enzyme complex as shown by crystallographic analysis. The results highlight the potential for inhibition of viral proteases employing nucleophilic catalysis by β-lactams and related acylating agents.

Covalent narlaprevir- and boceprevir-derived hybrid inhibitors of SARS-CoV-2 main protease

Emerging SARS-CoV-2 variants continue to threaten the effectiveness of COVID-19 vaccines, and small-molecule antivirals can provide an important therapeutic treatment option. The viral main protease (M^pro) is critical for virus replication and thus is considered an attractive drug target. We performed the design and characterization of three covalent hybrid inhibitors BBH-1, BBH-2 and NBH-2 created by splicing components of hepatitis C protease inhibitors boceprevir and narlaprevir, and known SARS-CoV-1 protease inhibitors. A joint X-ray/neutron structure of the M^pro/BBH-1 complex demonstrates that a Cys145 thiolate reaction with the inhibitor’s keto-warhead creates a negatively charged oxyanion. Protonation states of the ionizable residues in the M^pro active site adapt to the inhibitor, which appears to be an intrinsic property of M^pro. Structural comparisons of the hybrid inhibitors with PF-07321332 reveal unconventional F···O interactions of PF-07321332 with M^pro which may explain its more favorable enthalpy of binding. BBH-1, BBH-2 and NBH-2 exhibit comparable antiviral properties in vitro relative to PF-07321332, making them good candidates for further design of improved antivirals.

Three covalent hybrid inhibitors of SARS-CoV-2 main protease (Mpro) have been designed and compared to Pfizer’s nirmatrelvir (PF-07321332), providing atomic and thermodynamic details of their binding to the enzyme, and antiviral potency.

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

COVID-19 is a deadly rampaging infectious disease with over 480 million cases worldwide. Unfortunately, effective therapies remain very limited. Novel antiviral agents are urgently needed to combat this global healthcare crisis. Here, we elucidate the structural basis for replicase polyprotein cleavage and substrate specificity of SARS-CoV-2 main protease (M^pro). Through analyzing a series of high-resolution structures of SARS-CoV-2 M^pro throughout the proteolytic process, we demonstrate the molecular mechanism of M^pro in proteolytic processing that confers substrate specificity. Substrate selectivity is revealed using structures of the H41A mutant in complex with six individual native cleavage substrates. Our study underscores the mechanistic function of M^pro in the viral life cycle, which provides structural insights to develop effective inhibitors against this essential target of SARS-CoV-2.

The main protease (M^pro) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a key enzyme, which extensively digests CoV replicase polyproteins essential for viral replication and transcription, making it an attractive target for antiviral drug development. However, the molecular mechanism of how M^pro of SARS-CoV-2 digests replicase polyproteins, releasing the nonstructural proteins (nsps), and its substrate specificity remain largely unknown. Here, we determine the high-resolution structures of SARS-CoV-2 M^pro in its resting state, precleavage state, and postcleavage state, constituting a full cycle of substrate cleavage. The structures show the delicate conformational changes that occur during polyprotein processing. Further, we solve the structures of the SARS-CoV-2 M^pro mutant (H41A) in complex with six native cleavage substrates from replicase polyproteins, and demonstrate that SARS-CoV-2 M^pro can recognize sequences as long as 10 residues but only have special selectivity for four subsites. These structural data provide a basis to develop potent new inhibitors against SARS-CoV-2.

How Ethanolic Disinfectants Disintegrate Coronavirus Model Membranes: A Dissipative Particle Dynamics Simulation Study

We have developed dissipative particle dynamics models for pure dipalmitoylphosphatidylcholine (DPPC), dioleoylphosphatidylcholine (DOPC), and dimyristoylphosphatidylcholine (DMPC) as well as their binary and ternary mixed membranes, as coronavirus model membranes. The stabilities of pure and mixed membranes, surrounded by aqueous solutions containing up to 70 mol % ethanol (alcoholic disinfectants), have been investigated at room temperature. We found that aqueous solutions containing 5-10 mol % ethanol already have a significant weakening effect on the pure and mixed membranes. The magnitude of the effect depends on the membrane composition and the ethanol concentration. Ethanol permeabilizes the membrane, causing its lateral swelling and thickness shrinking and reducing the orientational order of the hydrocarbon tail of the bilayer. The free energy barrier for the permeation of ethanol in the bilayers is considerably reduced by the ethanol uptake. The rupture-critical ethanol concentrations causing the membrane failure are 20.7, 27.5, and 31.7 mol % in the aqueous phase surrounding pure DMPC, DOPC, and DPPC membranes, respectively. Characterizing the failure of lipid membranes by a machine-learning neural network framework, we found that all mixed binary and/or ternary membranes disrupt when immersed in an aqueous solution containing a rupture-critical ethanol concentration, ranging from 20.7 to 31.7 mol %, depending on the composition of the membrane; the DPPC-rich membranes are more intact, while the DMPC-rich membranes are least intact. Due to the tight packing of long, saturated hydrocarbon tails in DPPC, increasing the DPPC content of the mixed membrane increases its stability against the disinfectant. At high DPPC concentrations, where the DOPC and DMPC molecules are confined between the DPPC lipids, the ordered hydrocarbon tails of DPPC also induce order in the DOPC and DMPC molecules and, hence, stabilize the membrane more. Our simulations on pure and mixed membranes of a diversity of compositions reveal that a maximum ethanol concentration of 32 mol % (55 wt %) in the alcohol-based disinfectants is enough to disintegrate any membrane composed of these three lipids.

Inhibition of the main protease of SARS-CoV-2 (Mpro) by repurposing/designing drug-like substances and utilizing nature's toolbox of bioactive compounds

The emergence of the Severe Acute Respiratory Syndrome Coronavirus-2 (SARS-CoV-2) has resulted in a long pandemic, with numerous cases and victims worldwide and enormous consequences on social and economic life. Although vaccinations have proceeded and provide a valuable shield against the virus, the approved drugs are limited and it is crucial that further ways to combat infection are developed, that can also act against potential mutations. The main protease (Mpro) of the virus is an appealing target for the development of inhibitors, due to its importance in the viral life cycle and its high conservation among different coronaviruses. Several compounds have shown inhibitory potential against Mpro, both in silico and in vitro, with few of them also having entered clinical trials. These candidates include: known drugs that have been repurposed, molecules specifically designed based on the natural substrate of the protease or on structural moieties that have shown high binding affinity to the protease active site, as well as naturally derived compounds, either isolated or in plant extracts. The aim of this work is to collectively present the results of research regarding Mpro inhibitors to date, focusing on the function of the compounds founded by in silico simulations and further explored by in vitro and in vivo assays. Creating an extended portfolio of promising compounds that may block viral replication by inhibiting Mpro and by understanding involved structure-activity relationships, could provide a basis for the development of effective solutions against SARS-CoV-2 and future related outbreaks.

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

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

Covalent narlaprevir- and boceprevir-derived hybrid inhibitors of SARS-CoV-2 main protease: room-temperature X-ray and neutron crystallography, binding thermodynamics, and antiviral activity

The COVID-19 pandemic continues to disrupt everyday life, with constantly emerging SARS-CoV-2 variants threatening to render current vaccines ineffective. Small-molecule antivirals can provide an important therapeutic treatment option that is subject to challenges caused by the virus variants. The viral main protease (M pro ) is critical for the virus replication and thus is considered an attractive drug target for specific protease inhibitors. We performed the design and characterization of three reversible covalent hybrid inhibitors BBH-1, BBH-2 and NBH-2, whose structures were derived from those of hepatitis C protease inhibitors boceprevir and narlaprevir. A joint X-ray/neutron structure of the M pro /BBH-1 complex demonstrated that a Cys145 thiolate reaction with the inhibitor’s keto-warhead creates a negatively charged oxyanion, similar to that proposed for the M pro -catalyzed peptide bond hydrolysis. Protonation states of the ionizable residues in the M pro active site adapt to the inhibitor, which appears to be an intrinsic property of M pro . Structural comparisons of the hybrid inhibitors with PF-07321332 revealed unconventional interactions of PF-07321332 with M pro which may explain its more favorable enthalpy of binding and consequently higher potency. BBH-1, BBH-2 and NBH-2 demonstrated comparable antiviral properties in vitro relative to PF-07321332, making them good candidates for further design of improved antivirals.

Discovery of 9,10-dihydrophenanthrene derivatives as SARS-CoV-2 3CLpro inhibitors for treating COVID-19

The epidemic coronavirus disease 2019 (COVID-19) caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has now spread worldwide and efficacious therapeutics are urgently needed. 3-Chymotrypsin-like cysteine protease (3CLpro) is an indispensable protein in viral replication and represents an attractive drug target for fighting COVID-19. Herein, we report the discovery of 9,10-dihydrophenanthrene derivatives as non-peptidomimetic and non-covalent inhibitors of the SARS-CoV-2 3CLpro. The structure-activity relationships of 9,10-dihydrophenanthrenes as SARS-CoV-2 3CLpro inhibitors have carefully been investigated and discussed in this study. Among all tested 9,10-dihydrophenanthrene derivatives, C1 and C2 display the most potent SARS-CoV-2 3CLpro inhibition activity, with IC50 values of 1.55 ± 0.21 μM and 1.81 ± 0.17 μM, respectively. Further enzyme kinetics assays show that these two compounds dose-dependently inhibit SARS-CoV-2 3CLprovia a mixed-inhibition manner. Molecular docking simulations reveal the binding modes of C1 in the dimer interface and substrate-binding pocket of the target. In addition, C1 shows outstanding metabolic stability in the gastrointestinal tract, human plasma, and human liver microsome, suggesting that this agent has the potential to be developed as an orally administrated SARS-CoV-2 3CLpro inhibitor.

High-Throughput Virtual Screening and Validation of a SARS-CoV-2 Main Protease Noncovalent Inhibitor

Despite the recent availability of vaccines against the acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the search for inhibitory therapeutic agents has assumed importance especially in the context of emerging new viral variants. In this paper, we describe the discovery of a novel noncovalent small-molecule inhibitor, MCULE-5948770040, that binds to and inhibits the SARS-Cov-2 main protease (Mpro) by employing a scalable high-throughput virtual screening (HTVS) framework and a targeted compound library of over 6.5 million molecules that could be readily ordered and purchased. Our HTVS framework leverages the U.S. supercomputing infrastructure achieving nearly 91% resource utilization and nearly 126 million docking calculations per hour. Downstream biochemical assays validate this Mpro inhibitor with an inhibition constant (Ki) of 2.9 μM (95% CI 2.2, 4.0). Furthermore, using room-temperature X-ray crystallography, we show that MCULE-5948770040 binds to a cleft in the primary binding site of Mpro forming stable hydrogen bond and hydrophobic interactions. We then used multiple μs-time scale molecular dynamics (MD) simulations and machine learning (ML) techniques to elucidate how the bound ligand alters the conformational states accessed by Mpro, involving motions both proximal and distal to the binding site. Together, our results demonstrate how MCULE-5948770040 inhibits Mpro and offers a springboard for further therapeutic design.

Impact of Warhead Modulations on the Covalent Inhibition of SARS-CoV-2 Mpro Explored by QM/MM Simulations

The COVID-19 pandemic, caused by the severe acute respiratory syndrome coronavirus-2, SARS-CoV-2, shows the need for effective antiviral treatments. Here, we present a simulation study of the inhibition of the SARS-CoV-2 main protease (Mpro), a cysteine hydrolase essential for the life cycle of the virus. The free energy landscape for the mechanism of the inhibition process is explored by QM/MM umbrella sampling and free energy perturbation simulations at the M06-2X/MM level of theory for two proposed peptidyl covalent inhibitors that share the same recognition motif but feature distinct cysteine-targeting warheads. Regardless of the intrinsic reactivity of the modeled inhibitors, namely a Michael acceptor and a hydroxymethyl ketone activated carbonyl, our results confirm that the inhibitory process takes place by means of a two-step mechanism, in which the formation of an ion pair C145/H41 dyad precedes the protein-inhibitor covalent bond formation. The nature of this second step is strongly dependent on the functional groups in the warhead: while the nucleophilic attack of the C145 sulfur atom on the Cα of the double bond of the Michael acceptor takes place concertedly with the proton transfer from H41 to Cβ, in the compound with an activated carbonyl, the sulfur attacks the carbonyl carbon concomitant with a proton transfer from H41 to the carbonyl oxygen via the hydroxyl group. An analysis of the free energy profiles, structures along the reaction path, and interactions between the inhibitors and the different pockets of the active site on the protein shows a measurable effect of the warhead on the kinetics and thermodynamics of the process. These results and QM/MM methods can be used as a guide to select warheads to design efficient irreversible and reversible inhibitors of SARS-CoV-2 Mpro.

A Review of the Current Landscape of SARS-CoV-2 Main Protease Inhibitors: Have We Hit the Bullseye Yet?

In this review, we collected 1765 severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) M-pro inhibitors from the bibliography and other sources, such as the COVID Moonshot project and the ChEMBL database. This set of inhibitors includes only those compounds whose inhibitory capacity, mainly expressed as the half-maximal inhibitory concentration (IC50) value, against M-pro from SARS-CoV-2 has been determined. Several covalent warheads are used to treat covalent and non-covalent inhibitors separately. Chemical space, the variation of the IC50 inhibitory activity when measured by different methods or laboratories, and the influence of 1,4-dithiothreitol (DTT) are discussed. When available, we have collected the values of inhibition of viral replication measured with a cellular antiviral assay and expressed as half maximal effective concentration (EC50) values, and their possible relationship to inhibitory potency against M-pro is analyzed. Finally, the most potent covalent and non-covalent inhibitors that simultaneously inhibit the SARS-CoV-2 M-pro and the virus replication in vitro are discussed.

Structural, Electronic, and Electrostatic Determinants for Inhibitor Binding to Subsites S1 and S2 in SARS-CoV-2 Main Protease

Creating small-molecule antivirals specific for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) proteins is crucial to battle coronavirus disease 2019 (COVID-19). SARS-CoV-2 main protease (Mpro) is an established drug target for the design of protease inhibitors. We performed a structure-activity relationship (SAR) study of noncovalent compounds that bind in the enzyme's substrate-binding subsites S1 and S2, revealing structural, electronic, and electrostatic determinants of these sites. The study was guided by the X-ray/neutron structure of Mpro complexed with Mcule-5948770040 (compound 1), in which protonation states were directly visualized. Virtual reality-assisted structure analysis and small-molecule building were employed to generate analogues of 1. In vitro enzyme inhibition assays and room-temperature X-ray structures demonstrated the effect of chemical modifications on Mpro inhibition, showing that (1) maintaining correct geometry of an inhibitor's P1 group is essential to preserve the hydrogen bond with the protonated His163; (2) a positively charged linker is preferred; and (3) subsite S2 prefers nonbulky modestly electronegative groups.

An overview of potential inhibitors targeting non-structural proteins 3 (PLpro and Mac1) and 5 (3CLpro/Mpro) of SARS-CoV-2

There is an urgent need to develop effective treatments for coronavirus disease 2019 (COVID-19), which is caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). The rapid spread of SARS-CoV-2 has resulted in a global pandemic that has not only affected the daily lives of individuals but also had a significant impact on the global economy and public health. Although extensive research has been conducted to identify inhibitors targeting SARS-CoV-2, there are still no effective treatment strategies to combat COVID-19. SARS-CoV-2 comprises two important proteolytic enzymes, namely, the papain-like proteinase, located within non-structural protein 3 (nsp3), and nsp5, both of which cleave large replicase polypeptides into multiple fragments that are required for viral replication. Moreover, a domain within nsp3, known as the macrodomain (Mac1), also plays an important role in viral replication. Inhibition of their functions should be able to significantly interfere with the replication cycle of the virus, and therefore these key proteins may serve as potential therapeutic targets. The functions of the above viral targets and their corresponding inhibitors have been summarized in the current review. This review provides comprehensive updates of nsp3 and nsp5 inhibitor development and would help advance the discovery of novel anti-viral therapeutics against SARS-CoV-2.