A head-to-head comparison of the inhibitory activities of 15 peptidomimetic SARS-CoV-2 3CLpro inhibitors

A SAR and QSAR study on 3CLpro inhibitors of SARS-CoV-2 using machine learning methods

The 3C-like Proteinase (3CLpro) of novel coronaviruses is intricately linked to viral replication, making it a crucial target for antiviral agents. In this study, we employed two fingerprint descriptors (ECFP_4 and MACCS) to comprehensively characterize 889 compounds in our dataset. We constructed 24 classification models using machine learning algorithms, including Support Vector Machine (SVM), Random Forest (RF), extreme Gradient Boosting (XGBoost), and Deep Neural Networks (DNN). Among these models, the DNN- and ECFP_4-based Model 1D_2 achieved the most promising results, with a remarkable Matthews correlation coefficient (MCC) value of 0.796 in the 5-fold cross-validation and 0.722 on the test set. The application domains of the models were analysed using dSTD-PRO calculations. The collected 889 compounds were clustered by K-means algorithm, and the relationships between structural fragments and inhibitory activities of the highly active compounds were analysed for the 10 obtained subsets. In addition, based on 464 3CLpro inhibitors, 27 QSAR models were constructed using three machine learning algorithms with a minimum root mean square error (RMSE) of 0.509 on the test set. The applicability domains of the models and the structure-activity relationships responded from the descriptors were also analysed.

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.

In-silico study: docking simulation and molecular dynamics of peptidomimetic fullerene-based derivatives against SARS-CoV-2 Mpro

COVID-19 is caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2, has become a global pandemic resulting in significant morbidity and mortality. This study presents 12 new peptidomimetic fullerene-based derivatives in three groups that are investigated theoretically as SARS-CoV-2 M^pro inhibitors to increase the chance of treating COVID-19. Studied compounds are designed and optimized at B88-LYP/DZVP method. Molecular descriptors results show the stability and reactivity of the compounds with M^pro, especially in the 3rd group (Ser compounds). However, Lipinski's Rule of Five values indicates that the compounds are not suitable as oral drugs. Furthermore, molecular docking simulations are carried out to investigate the binding affinity and interaction modes of the top five compounds (compounds 1, 9, 11, 2, and 10) with the M^pro protein, which have the lowest binding energy. Molecular dynamics simulations are also performed to evaluate the stability of the protein-ligand complexes with compounds 1 and 9 and compare them with natural substrate interaction. The analysis of RMSD, H-bonds, Rg, and SASA indicates that both compounds 1 (Gly-α acid) and 9 (Ser-α acid) have good stability and strong binding affinity with the M^pro protein. However, compound 9 shows slightly better stability and binding affinity compared to compound 1.

Multidimensional virtual screening approaches combined with drug repurposing to identify potential covalent inhibitors of SARS-CoV-2 3CL protease

The outbreak of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) has caused an unprecedented global pandemic, and new cases are still on the rise due to the absence of effective medicines. However, developing new drugs within a short time is extremely difficult. Repurposing the existing drugs provides a fast and effective strategy to identify promising inhibitors. Here we focus on the screening of drugs database for discovering potential covalent inhibitors that target 3-chymotrypsin-like protease (3CLpro), an essential enzyme mediating viral replication and transcription. Firstly, we constructed a receptor-ligand pharmacophore model and verified it through decoy set. The importance of pharmacophore features was evaluated by combining molecular dynamics simulation with interaction analyses. Then, covalent docking was used to perform further screening. According to docking score and Prime/Molecular Mechanics Generalized Born Surface Area (MM-GBSA) score, total ten compounds obtained good scores and successfully established covalent bonds with the catalytic Cys145 residue. They also formed favorable interactions with key residues in active sites and closely integrated with 3CLpro with binding modes similar to known 3CLpro inhibitor. Finally, the top four hits DB08732, DB04653, DB01871 and DB07299 were further subjected to 100 ns molecular dynamics (MD) simulation and MM-GBSA binding free energy calculations. The results suggest that the four candidates show good binding affinities for 3CLpro, which warrants further evaluation for their in-vitro/in-vivo activities. Overall, our research methods provide a valuable reference for discovering promising inhibitors against SARS-CoV-2 and help to fight against the epidemic.Communicated by Ramaswamy H. Sarma.

Enhanced compound-protein binding affinity prediction by representing protein multimodal information via a coevolutionary strategy

Due to the lack of a method to efficiently represent the multimodal information of a protein, including its structure and sequence information, predicting compound-protein binding affinity (CPA) still suffers from low accuracy when applying machine-learning methods. To overcome this limitation, in a novel end-to-end architecture (named FeatNN), we develop a coevolutionary strategy to jointly represent the structure and sequence features of proteins and ultimately optimize the mathematical models for predicting CPA. Furthermore, from the perspective of data-driven approach, we proposed a rational method that can utilize both high- and low-quality databases to optimize the accuracy and generalization ability of FeatNN in CPA prediction tasks. Notably, we visually interpret the feature interaction process between sequence and structure in the rationally designed architecture. As a result, FeatNN considerably outperforms the state-of-the-art (SOTA) baseline in virtual drug evaluation tasks, indicating the feasibility of this approach for practical use. FeatNN provides an outstanding method for higher CPA prediction accuracy and better generalization ability by efficiently representing multimodal information of proteins via a coevolutionary strategy.

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.

Screening of Phytoconstituents from Traditional Plants against SARSCoV- 2 using Molecular Docking Approach

Background:

The emergence of COVID-19 as a fatal viral disease encourages researchers to develop effective and efficient therapeutic agents. The intervention of in silico studies has led to revolutionary changes in the conventional method of testing the bioactivity of plant constituents.

Objective:

The current study deals with the investigation of some traditional immunomodulators of plant origin to combat this ailment.

Materials and Methods:

A total of 151 phytomolecules of 12 immunomodulatory plants were evaluated for their inhibitory action against the main protease (PDB ID: 7D1M) and NSP15 endoribonuclease (PDB ID: 6WLC) by structure-based virtual screening. In addition, the promising molecules with ligand efficiency of more than -0.3(kcal/mol)/heavy atoms were further predicted for pharmacokinetic properties and druggability using the SwissADME web server, and their toxicity was also evaluated using Protox-II.

Result:

Myricetin-3-O-arabinofuranoside of cranberry plant was found to be the most potential candidate against both enzymes: main protease (–14.2 kcal/mol) and NSP15 endoribonuclease (–12.2 kcal/mol).

Conclusion:

The promising outcomes of the current study may be implemented in future drug development against coronavirus. The findings also help in the development of lead candidates of plant origin with a better ADMET profile in the future.

A systematic exploration of boceprevir-based main protease inhibitors as SARS-CoV-2 antivirals

Boceprevir is an HCV NSP3 inhibitor that was explored as a repurposed drug for COVID-19. It inhibits the SARS-CoV-2 main protease (MPro) and contains an α-ketoamide warhead, a P1 β-cyclobutylalanyl moiety, a P2 dimethylcyclopropylproline, a P3 tert-butylglycine, and a P4 N-terminal tert-butylcarbamide. By introducing modifications at all four positions, we synthesized 20 boceprevir-based MPro inhibitors including PF-07321332 and characterized their MPro inhibition potency in test tubes (in vitro) and 293T cells (in cellulo). Crystal structures of MPro bound with 10 inhibitors and cytotoxicity and antiviral potency of 4 inhibitors were characterized as well. Replacing the P1 site with a β-(S-2-oxopyrrolidin-3-yl)-alanyl (Opal) residue and the warhead with an aldehyde leads to high in vitro potency. The original moieties at P2, P3 and the P4 N-terminal cap positions in boceprevir are better than other tested chemical moieties for high in vitro potency. In crystal structures, all inhibitors form a covalent adduct with the MPro active site cysteine. The P1 Opal residue, P2 dimethylcyclopropylproline and P4 N-terminal tert-butylcarbamide make strong hydrophobic interactions with MPro, explaining high in vitro potency of inhibitors that contain these moieties. A unique observation was made with an inhibitor that contains a P4 N-terminal isovaleramide. In its MPro complex structure, the P4 N-terminal isovaleramide is tucked deep in a small pocket of MPro that originally recognizes a P4 alanine side chain in a substrate. Although all inhibitors show high in vitro potency, they have drastically different in cellulo potency to inhibit ectopically expressed MPro in human 293T cells. In general, inhibitors with a P4 N-terminal carbamide or amide have low in cellulo potency. This trend is reversed when the P4 N-terminal cap is changed to a carbamate. The installation of a P3 O-tert-butyl-threonine improves in cellulo potency. Three molecules that contain a P4 N-terminal carbamate were advanced to cytotoxicity tests on 293T cells and antiviral potency tests on three SARS-CoV-2 variants. They all have relatively low cytotoxicity and high antiviral potency with EC50 values around 1 μM. A control compound with a nitrile warhead and a P4 N-terminal amide has undetectable antiviral potency. Based on all observations, we conclude that a P4 N-terminal carbamate in a boceprevir derivative is key for high antiviral potency against SARS-CoV-2.

A Warhead Substitution Study on the Coronavirus Main Protease Inhibitor Nirmatrelvir

The SARS-CoV-2 pandemic is currently causing an unprecedented global health emergency since its emergence in December 2019. In December 2021, the FDA granted emergency use authorization to nirmatrelvir, a SARS-CoV-2 main protease inhibitor, for treating infected patients. This peptidomimetic is designed with a nitrile warhead, which forms a covalent bond to the viral protease. Herein, we investigate nirmatrelvir analogs with different warheads and their inhibitory activities. In addition, antiviral activities against human alphacoronavirus 229E was also investigated along with a cell-based assay. We discovered that the hydroxymethylketone and ketobenzothiazole warheads were equipotent to the nitrile warhead, suggesting that these analogs can also be used for treating coronavirus infections.

Evolutionary and Structural Insights about Potential SARS-CoV-2 Evasion of Nirmatrelvir

The U.S. FDA approval of PAXLOVID, a combination therapy of nirmatrelvir and ritonavir has significantly boosted our morale in fighting the COVID-19 pandemic. Nirmatrelvir is an inhibitor of the main protease (MPro) of SARS-CoV-2. Since many SARS-CoV-2 variants that resist vaccines and antibodies have emerged, a concern of acquired viral resistance to nirmatrelvir naturally arises. Here, possible mutations in MPro to confer viral evasion of nirmatrelvir are analyzed and discussed from both evolutionary and structural standpoints. The analysis indicates that those mutations will likely reside in the whole aa45-51 helical region and residues including M165, L167, P168, R188, and Q189. Relevant mutations have also been observed in existing SARS-CoV-2 samples. Implications of this analysis to the fight against future drug-resistant viral variants and the development of broad-spectrum antivirals are discussed as well.

Efficient Identification of Anti-SARS-CoV-2 Compounds Using Chemical Structure- and Biological Activity-Based Modeling

Identification of anti-SARS-CoV-2 compounds through traditional high-throughput screening (HTS) assays is limited by high costs and low hit rates. To address these challenges, we developed machine learning models to identify compounds acting via inhibition of the entry of SARS-CoV-2 into human host cells or the SARS-CoV-2 3-chymotrypsin-like (3CL) protease. The optimal classification models achieved good performance with area under the receiver operating characteristic curve (AUC-ROC) values of >0.78. Experimental validation showed that the best performing models increased the assay hit rate by 2.1-fold for viral entry inhibitors and 10.4-fold for 3CL protease inhibitors compared to those of the original drug repurposing screens. Twenty-two compounds showed potent (<5 μM) antiviral activities in a SARS-CoV-2 live virus assay. In conclusion, machine learning models can be developed and used as a complementary approach to HTS to expand compound screening capacities and improve the speed and efficiency of anti-SARS-CoV-2 drug discovery.

Antiviral cyclic peptides targeting the main protease of SARS-CoV-2

Antivirals that specifically target SARS-CoV-2 are needed to control the COVID-19 pandemic. The main protease (M^pro) is essential for SARS-CoV-2 replication and is an attractive target for antiviral development. Here we report the use of the Random nonstandard Peptide Integrated Discovery (RaPID) mRNA display on a chemically cross-linked SARS-CoV-2 M^pro dimer, which yielded several high-affinity thioether-linked cyclic peptide inhibitors of the protease. Structural analysis of M^pro complexed with a selenoether analogue of the highest-affinity peptide revealed key binding interactions, including glutamine and leucine residues in sites S1 and S2, respectively, and a binding epitope straddling both protein chains in the physiological dimer. Several of these M^pro peptide inhibitors possessed antiviral activity against SARS-CoV-2 in vitro with EC₅₀ values in the low micromolar range. These cyclic peptides serve as a foundation for the development of much needed antivirals that specifically target SARS-CoV-2.

RaPID mRNA display was used for the discovery of antiviral cyclic peptides that potently and selectively inhibit SARS-CoV-2 M^pro. The most potent inhibitor exhibited a novel binding mode, interacting with residues across the homodimer interface.

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.

The N -Terminal Carbamate is Key to High Cellular and Antiviral Potency for Boceprevir-Based SARS-CoV-2 Main Protease Inhibitors

Boceprevir is an HCV NSP3 inhibitor that has been explored as a repurposed drug for COVID-19. It inhibits the SARS-CoV-2 main protease (M ^Pro ) and contains an α-ketoamide warhead, a P1 β-cyclobutylalanyl moiety, a P2 dimethylcyclopropylproline, a P3 tert -butyl-glycine, and a P4 N -terminal tert -butylcarbamide. By introducing modifications at all four positions, we synthesized 20 boceprevir-based M ^Pro inhibitors including PF-07321332 and characterized their M ^Pro inhibition potency in test tubes ( in vitro ) and human host cells ( in cellulo ). Crystal structures of M ^Pro bound with 10 inhibitors and antiviral potency of 4 inhibitors were characterized as well. Replacing the P1 site with a β-(S-2-oxopyrrolidin-3-yl)-alanyl (opal) residue and the warhead with an aldehyde leads to high in vitro potency. The original moieties at P2, P3 and the P4 N -terminal cap positions in boceprevir are better than other tested chemical moieties for high in vitro potency. In crystal structures, all inhibitors form a covalent adduct with the M ^Pro active site cysteine. The P1 opal residue, P2 dimethylcyclopropylproline and P4 N -terminal tert -butylcarbamide make strong hydrophobic interactions with M ^Pro , explaining high in vitro potency of inhibitors that contain these moieties. A unique observation was made with an inhibitor that contains an P4 N -terminal isovaleramide. In its M ^Pro complex structure, the P4 N -terminal isovaleramide is tucked deep in a small pocket of M ^Pro that originally recognizes a P4 alanine side chain in a substrate. Although all inhibitors show high in vitro potency, they have drastically different in cellulo potency in inhibiting ectopically expressed M ^Pro in human 293T cells. All inhibitors including PF-07321332 with a P4 N -terminal carbamide or amide have low in cellulo potency. This trend is reversed when the P4 N -terminal cap is changed to a carbamate. The installation of a P3 O-tert -butyl-threonine improves in cellulo potency. Three molecules that contain a P4 N -terminal carbamate were advanced to antiviral tests on three SARS-CoV-2 variants. They all have high potency with EC ₅₀ values around 1 μM. A control compound with a nitrile warhead and a P4 N -terminal amide has undetectable antiviral potency. Based on all observations, we conclude that a P4 N -terminal carbamate in a boceprevir derivative is key for high antiviral potency against SARS-CoV-2.

A Patent Review on SARS Coronavirus Main Protease (3CLpro ) Inhibitors

The severe acute respiratory syndrome coronavirus 2 (SARS‐CoV‐2) pandemic is an unprecedented global health emergency causing more than 4.2 million fatalities as of 30 July 2021. Only three antiviral therapies have been approved or granted emergency use authorization by the FDA. The SARS‐CoV‐2 3CL protease (3CL^pro) is deemed an attractive drug target as it plays an essential role in viral polyprotein processing and pathogenesis, although no inhibitors have been approved. This patent review discusses SARS coronavirus 3CL^pro inhibitors that have been filed up to 30 July 2021, giving an overview on the types of inhibitors that have generated commercial interest, especially amongst drug companies. Insights into the common structural motifs required for active site binding is also discussed.