Potent Noncovalent Inhibitors of the Main Protease of SARS-CoV-2 from Molecular Sculpting of the Drug Perampanel Guided by Free Energy Perturbation Calculations

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.

Surely you are joking, Mr Docking!

In the wake of recent COVID-19 pandemics scientists around the world rushed to deliver numerous CADD (Computer-Aided Drug Discovery) methods and tools that could be reliably used to discover novel drug candidates against the SARS-CoV-2 virus. With that, there emerged a trend of a significant democratization of CADD that contributed to the rapid development of various COVID-19 drug candidates currently undergoing different stages of validation. On the other hand, this democratization also inadvertently led to the surge rapidly performed molecular docking studies to nominate multiple scores of novel drug candidates supported by computational arguments only. Albeit driven by best intentions, most of such studies also did not follow best practices in the field that require experience and expertise learned through multiple rigorously designed benchmarking studies and rigorous experimental validation. In this Viewpoint we reflect on recent disbalance between small number of rigorous and comprehensive studies and the proliferation of purely computational studies enabled by the ease of docking software availability. We further elaborate on the hyped oversale of CADD methods' ability to rapidly yield viable drug candidates and reiterate the critical importance of rigor and adherence to the best practices of CADD in view of recent emergence of AI and Big Data in the field.

A Systematic Survey of Reversibly Covalent Dipeptidyl Inhibitors of the SARS-CoV-2 Main Protease

SARS-CoV-2 is the coronavirus pathogen of the currently prevailing COVID-19 pandemic. It relies on its main protease (M Pro ) for replication and pathogenesis. M Pro is a demonstrated target for the development of antivirals for SARS-CoV-2. Past studies have systematically explored tripeptidyl inhibitors such as nirmatrelvir as M Pro inhibitors. However, dipeptidyl inhibitors especially those with a spiro residue at their P2 position have not been systematically investigated. In this work, we synthesized about 30 reversibly covalent dipeptidyl M Pro inhibitors and characterized them on in vitro enzymatic inhibition potency, structures of their complexes with M Pro , cellular M Pro inhibition potency, antiviral potency, cytotoxicity, and in vitro metabolic stability. Our results indicated that M Pro has a flexible S2 pocket that accommodates dipeptidyl inhibitors with a large P2 residue and revealed that dipeptidyl inhibitors with a large P2 spiro residue such as ( S )-2-azaspiro[4,4]nonane-3-carboxylate and ( S )-2-azaspiro[4,5]decane-3-carboxylate have optimal characteristics. One compound MPI60 containing a P2 ( S )-2-azaspiro[4,4]nonane-3-carboxylate displayed high antiviral potency, low cellular cytotoxicity, and high in vitro metabolic stability and can be potentially advanced to further preclinical tests.

Chetomin, a SARS-CoV-2 3C-like Protease (3CLpro) Inhibitor: In Silico Screening, Enzyme Docking, Molecular Dynamics and Pharmacokinetics Analysis

The emergence of the Coronavirus Disease 2019 (COVID-19) pandemic caused by severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) has led to over 6 million deaths. The 3C-like protease (3CLpro) enzyme of the SARS-CoV-2 virus is an attractive druggable target for exploring therapeutic drug candidates to combat COVID-19 due to its key function in viral replication. Marine natural products (MNPs) have attracted considerable attention as alternative sources of antiviral drug candidates. In looking for potential 3CLpro inhibitors, the MNP database (>14,000 molecules) was virtually screened against 3CLpro with the assistance of molecular docking computations. The performance of AutoDock and OEDocking software in anticipating the ligand-3CLpro binding mode was first validated according to the available experimental data. Based on the docking scores, the most potent MNPs were further subjected to molecular dynamics (MD) simulations, and the binding affinities of those molecules were computed using the MM-GBSA approach. According to MM-GBSA//200 ns MD simulations, chetomin (UMHMNP1403367) exhibited a higher binding affinity against 3CLpro than XF7, with ΔGbinding values of −55.5 and −43.7 kcal/mol, respectively. The steadiness and tightness of chetomin with 3CLpro were evaluated, revealing the high stabilization of chetomin (UMHMNP1403367) inside the binding pocket of 3CLpro throughout 200 ns MD simulations. The physicochemical and pharmacokinetic features of chetomin were also predicted, and the oral bioavailability of chetomin was demonstrated. Furthermore, the potentiality of chetomin analogues −namely, chetomin A-D− as 3CLpro inhibitors was investigated. These results warrant further in vivo and in vitro assays of chetomin (UMHMNP1403367) as a promising anti-COVID-19 drug candidate.

Synthesis of conjugates of (aR,7S)-colchicine with monoterpenoids and investigation of their biological activity

Conjugates of the natural alkaloid (aR,7S)-colchicine with bicyclic monoterpenoids and their derivatives were synthesized for the first time. Molecular docking of the synthesized agents in the active site of the main viral protease of the SARS-CoV-2 virus was carried out. The cytotoxic properties of the agents against different cell lines and the ability to inhibit the main viral protease 3CLPro were studied.

Searching for potential inhibitors of SARS-COV-2 main protease using supervised learning and perturbation calculations

Inhibiting the biological activity of SARS-CoV-2 Mpro can prevent viral replication. In this context, a hybrid approach using knowledge- and physics-based methods was proposed to characterize potential inhibitors for SARS-CoV-2 Mpro. Initially, supervised machine learning (ML) models were trained to predict a ligand-binding affinity of ca. 2 million compounds with the correlation on a test set of R=0.748±0.044. Atomistic simulations were then used to refine the outcome of the ML model. Using LIE/FEP calculations, nine compounds from the top 100 ML inhibitors were suggested to bind well to the protease with the domination of van der Waals interactions. Furthermore, the binding affinity of these compounds is also higher than that of nirmatrelvir, which was recently approved by the US FDA to treat COVID-19. In addition, the ligands altered the catalytic triad Cys145 - His41 - Asp187, possibly disturbing the biological activity of SARS-CoV-2.

Exploring the dual effect of novel 1,4-diarylpyranopyrazoles as antiviral and anti-inflammatory for the management of SARS-CoV-2 and associated inflammatory symptoms

COVID-19 and associated substantial inflammations continue to threaten humankind triggering death worldwide. So, the development of new effective antiviral and anti-inflammatory medications is a major scientific goal. Pyranopyrazoles have occupied a crucial position in medicinal chemistry because of their biological importance. Here, we report the design and synthesis of a series of sixteen pyranopyrazole derivatives substituted with two aryl groups at N-1 and C-4. The designed compounds are suggested to show dual activity to combat the emerging Coronaviruses and associated substantial inflammations. All compounds were evaluated for their in vitro antiviral activity and cytotoxicity against SARS-CoV infected Vero cells. As well, the in vitro assay of all derivatives against the SARS-CoV M^pro target was performed. Results revealed the potential of three pyranopyrazoles (22, 27, and 31) to potently inhibit the viral main protease with IC₅₀ values of 2.01, 1.83, and 4.60 μM respectively compared with 12.85 and 82.17 μM for GC-376 and lopinavir. Additionally, in vivo anti-inflammatory testing for the most active compound 27 proved its ability to reduce levels of two cytokines (TNF-α and IL-6). Molecular docking and dynamics simulation revealed consistent results with the in vitro enzymatic assay and indicated the stability of the putative complex of 27 with SARS-CoV-2 M^pro. The assessment of metabolic stability and physicochemical properties of 27 have also been conducted. This investigation identified a set of metabolically stable pyranopyrazoles as effective anti-SARS-CoV-2 M^pro and suppressors of host cell cytokine release. We believe that the new compounds deserve further chemical optimization and evaluation for COVID-19 treatment.

Synthesis of non-symmetric N-benzylbispidinol amides and study of their inhibitory activity against the main protease of the SARS-CoV-2 virus

Based on the data obtained by molecular modeling of the non-covalent interaction of non-symmetric N-benzylbispidin-9-ol amides with the active site of the main protease 3CLpro of the SARS-CoV-2 virus, a series of compounds was synthesized, and their inhibitory activity against 3CLpro was studied and compared with that of the known inhibitor ML188 (IC₅₀ = 1.56±0.55 µmol L^-1). It was found that only compound 1g containing the 1,4-dihydroindeno[1,2-c]pyrazole fragment showed moderate activity (IC₅₀ = 100±5.7µmol L^-1) and was characterized by the highest calculated binding energy among the studied bispidine derivatives according to molecular docking data.

Gene set proximity analysis: expanding gene set enrichment analysis through learned geometric embeddings, with drug-repurposing applications in COVID-19

MOTIVATION

Gene set analysis methods rely on knowledge-based representations of genetic interactions in the form of both gene set collections and protein-protein interaction (PPI) networks. However, explicit representations of genetic interactions often fail to capture complex interdependencies among genes, limiting the analytic power of such methods.

RESULTS

We propose an extension of gene set enrichment analysis to a latent embedding space reflecting PPI network topology, called gene set proximity analysis (GSPA). Compared with existing methods, GSPA provides improved ability to identify disease-associated pathways in disease-matched gene expression datasets, while improving reproducibility of enrichment statistics for similar gene sets. GSPA is statistically straightforward, reducing to a version of traditional gene set enrichment analysis through a single user-defined parameter. We apply our method to identify novel drug associations with SARS-CoV-2 viral entry. Finally, we validate our drug association predictions through retrospective clinical analysis of claims data from 8 million patients, supporting a role for gabapentin as a risk factor and metformin as a protective factor for severe COVID-19.

AVAILABILITY AND IMPLEMENTATION

GSPA is available for download as a command-line Python package at https://github.com/henrycousins/gspa.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Virtual screening and structure-based 3D pharmacophore approach to identify small-molecule inhibitors of SARS-CoV-2 Mpro

The outbreak caused by a coronavirus 2 has required quick and potential treatment strategies. The main protease enzyme M^pro plays an important role in the viral replication which renders it an important target for discovering SARS-CoV-2 inhibitors. In this study, 3D pharmacophore structure-based virtual screening and molecular docking were conducted using MOE and Bristol University Docking Engine (BUDE). Around 400,000 molecules of ZINC15 database were docked against the crystal structure of main protease, followed by 3D pharmacophore filtration. Six top-ranked hits (ZINC58717986, ZINC60399606, ZINC58662884, ZINC45988635, ZINC54706757 and ZINC17320595) were identified based on their strong spatial affinity and forming H-bonds with key residues H41, E166, Q189 and T190 of the binding pocket of M^pro SARS-CoV-2. The 6 hits subjected to molecular dynamics simulations for 100 ns followed by binding free energy calculations using MM-PBSA technique. Interestingly, three hits showed free binding energy (ΔG_binding) lower than tert-butyl N-[1-[(2S)-1-[[(2S)-4-(benzylamino)-3,4-dioxo-1-[(3S)-2-oxopyrrolidin-3-yl]butan-2-yl]amino]-3-cyclopropyl-1-oxopropan-2-yl]-2-oxopyridin-3-yl]carbamate (α-ketoamide 13 b) (ΔG_binding) -76.67 ± 0.5 kJ/mol which suggested their potential against SARS-CoV-2. The best binding free energy candidates, ZINC58717986 (ΔG_binding) -98.41 ± 0.7 kJ/mol. The second-best hit candidate, ZINC54706757 (ΔG_binding) -83.4 ± 0.6 kJ/mol and the third one ZINC17320595 (ΔG_binding) -78.85 ± 0.5 kJ/mol. Per residue decomposition free energy indicates H41, S46, H164, E166, D187, Q189 and Q192 are hot spot residues while residues M49, M165, L167 and P168 contribute to the hydrophobic interactions. The pharmacokinetic study suggests that the selected 6 hits possess drug-like properties. The 3D pharmacophore virtual screening, molecular dynamics and MM-PBSA approaches facilitated identification 3 promising hits with low free binding energy as SARS-CoV-2 inhibitors.Communicated by Ramaswamy H. Sarma.

Discovery of novel inhibitors of SARS-CoV-2 main protease

Corona Virus Disease 2019 (COVID-19), referred to as 'New Coronary Pneumonia', is a type of acute infectious disease caused by the Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) infection. M^pro is one of the main targets for treating COVID-19. The current research on M^pro mainly focuses on the repurposing of old drugs, and there are only a few novel ligands that inhibit M^pro. In this research, we used computational free energy calculation to screen a compound library against M^pro, and discovered four novel compounds with the two best compounds (AG-690/13507628 and AG-690/13507724) having experimental measured IC₅₀ of just under 3 μM and low cell toxicity. Detailed decomposition of the interactions between the inhibitors and M^pro reveals key interacting residues and interactions that determine the activity. The results from this study should provide a basis for further development of anti-SARS-CoV-2 drugs.Communicated by Ramaswamy H. Sarma.

Structure-Based Identification of Naphthoquinones and Derivatives as Novel Inhibitors of Main Protease Mpro and Papain-like Protease PLpro of SARS-CoV-2

The worldwide COVID-19 pandemic caused by the coronavirus SARS-CoV-2 urgently demands novel direct antiviral treatments. The main protease (Mpro) and papain-like protease (PLpro) are attractive drug targets among coronaviruses due to their essential role in processing the polyproteins translated from the viral RNA. In this study, we virtually screened 688 naphthoquinoidal compounds and derivatives against Mpro of SARS-CoV-2. Twenty-four derivatives were selected and evaluated in biochemical assays against Mpro using a novel fluorogenic substrate. In parallel, these compounds were also assayed with SARS-CoV-2 PLpro. Four compounds inhibited Mpro with half-maximal inhibitory concentration (IC50) values between 0.41 μM and 9.0 μM. In addition, three compounds inhibited PLpro with IC50 ranging from 1.9 μM to 3.3 μM. To verify the specificity of Mpro and PLpro inhibitors, our experiments included an assessment of common causes of false positives such as aggregation, high compound fluorescence, and inhibition by enzyme oxidation. Altogether, we confirmed novel classes of specific Mpro and PLpro inhibitors. Molecular dynamics simulations suggest stable binding modes for Mpro inhibitors with frequent interactions with residues in the S1 and S2 pockets of the active site. For two PLpro inhibitors, interactions occur in the S3 and S4 pockets. In summary, our structure-based computational and biochemical approach identified novel naphthoquinonal scaffolds that can be further explored as SARS-CoV-2 antivirals.

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

Vaccines and drugs have helped reduce disease severity and blunt the spread of 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/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.

Optimization of potential non-covalent inhibitors for the SARS-CoV-2 main protease inspected by a descriptor of the subpocket occupancy

The main protease is regarded as an essential drug target for treating Coronavirus Disease 2019. In the present study, 13 marketed drugs were investigated to explore the possible binding mechanism, utilizing molecular docking, molecular dynamics simulation, and MM-PB(GB)SA binding energy calculations. Our results suggest that fusidic acid, polydatin, SEN-1269, AZD6482, and UNC-2327 have high binding affinities of more than 23 kcal mol^-1. A descriptor was defined for the energetic occupancy of the subpocket, and it was found that S4 had a low occupancy of less than 10% on average. The molecular optimization of ADZ6482 via reinforcement learning algorithms was carried out to screen out three lead compounds, in which slight structural changes give more considerable binding energies and an occupancy of the S4 subpocket of up to 43%. The energetic occupancy could be a useful descriptor for evaluating the local binding affinity for drug design.

Characterizing the ligand-binding affinity toward SARS-CoV-2 Mpro via physics- and knowledge-based approaches

Computational approaches, including physics- and knowledge-based methods, have commonly been used to determine the ligand-binding affinity toward SARS-CoV-2 main protease (Mpro or 3CLpro). Strong binding ligands can thus be suggested as potential inhibitors for blocking the biological activity of the protease. In this context, this paper aims to provide a short review of computational approaches that have recently been applied in the search for inhibitor candidates of Mpro. In particular, molecular docking and molecular dynamics (MD) simulations are usually combined to predict the binding affinity of thousands of compounds. Quantitative structure-activity relationship (QSAR) is the least computationally demanding and therefore can be used for large chemical collections of ligands. However, its accuracy may not be high. Moreover, the quantum mechanics/molecular mechanics (QM/MM) method is most commonly used for covalently binding inhibitors, which also play an important role in inhibiting the activity of SARS-CoV-2. Furthermore, machine learning (ML) models can significantly increase the searching space of ligands with high accuracy for binding affinity prediction. Physical insights into the binding process can then be confirmed via physics-based calculations. Integration of ML models into computational chemistry provides many more benefits and can lead to new therapies sooner.

3CLpro inhibitors: DEL-based molecular generation

Molecular generation (MG) via machine learning (ML) has speeded drug structural optimization, especially for targets with a large amount of reported bioactivity data. However, molecular generation for structural optimization is often powerless for new targets. DNA-encoded library (DEL) can generate systematic, target-specific activity data, including novel targets with few or unknown activity data. Therefore, this study aims to overcome the limitation of molecular generation in the structural optimization for the new target. Firstly, we generated molecules using the structure-affinity data (2.96 million samples) for 3C-like protease (3CLpro) from our own-built DEL platform to get rid of using public databases (e.g., CHEMBL and ZINC). Subsequently, to analyze the effect of transfer learning on the positive rate of the molecule generation model, molecular docking and affinity model based on DEL data were applied to explore the enhanced impact of transfer learning on molecule generation. In addition, the generated molecules are subjected to multiple filtering, including physicochemical properties, drug-like properties, and pharmacophore evaluation, molecular docking to determine the molecules for further study and verified by molecular dynamics simulation.

[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.

Bat coronaviruses related to SARS-CoV-2: what about their 3CL proteases (MPro)?

Despite a huge effort by the scientific community to determine the animal reservoir of SARS-CoV-2, which led to the identification of several SARS-CoV-2-related viruses both in bats and in pangolins, the origin of SARS-CoV-2 is still not clear. Recently, Temmam et al. reported the discovery of bat coronaviruses with a high degree of genome similarity with SARS-CoV-2, especially concerning the RBDs of the S protein, which mediates the capability of such viruses to enter and therefore infect human cells through a hACE2-dependent pathway. These viruses, especially the one named BANAL-236, showed a higher affinity for the hACE2 compared to the original strain of SARS-CoV-2. In the present work, we analyse the similarities and differences between the 3CL protease (main protease, Mpro) of these newly reported viruses and SARS-CoV-2, discussing their relevance relative to the efficacy of existing therapeutic approaches against COVID-19, particularly concerning the recently approved orally available Paxlovid, and the development of future ones.

From the Wuhan-Hu-1 strain to the XD and XE variants: is targeting the SARS-CoV-2 spike protein still a pharmaceutically relevant option against COVID-19?

Since the outbreak of the COVID-19 pandemic in December 2019, the SARS-CoV-2 genome has undergone several mutations. The emergence of such variants has resulted in multiple pandemic waves, contributing to sustaining to date the number of infections, hospitalisations, and deaths despite the swift development of vaccines, since most of these mutations are concentrated on the Spike protein, a viral surface glycoprotein that is the main target for most vaccines. A milestone in the fight against the COVID-19 pandemic has been represented by the development of Paxlovid, the first orally available drug against COVID-19, which acts on the Main Protease (Mpro). In this article, we analyse the structural features of both the Spike protein and the Mpro of the recently reported SARS-CoV-2 variant XE, as well the closely related XD and XF ones, discussing their impact on the efficacy of existing treatments against COVID-19 and on the development of future ones.

Dual Inhibitors of Main Protease (MPro) and Cathepsin L as Potent Antivirals against SARS-CoV2

Given the current impact of SARS-CoV2 and COVID-19 on human health and the global economy, the development of direct acting antivirals is of paramount importance. Main protease (MPro), a cysteine protease that cleaves the viral polyprotein, is essential for viral replication. Therefore, MPro is a novel therapeutic target. We identified two novel MPro inhibitors, D-FFRCMKyne and D-FFCitCMKyne, that covalently modify the active site cysteine (C145) and determined cocrystal structures. Medicinal chemistry efforts led to SM141 and SM142, which adopt a unique binding mode within the MPro active site. Notably, these inhibitors do not inhibit the other cysteine protease, papain-like protease (PLPro), involved in the life cycle of SARS-CoV2. SM141 and SM142 block SARS-CoV2 replication in hACE2 expressing A549 cells with IC50 values of 8.2 and 14.7 nM. Detailed studies indicate that these compounds also inhibit cathepsin L (CatL), which cleaves the viral S protein to promote viral entry into host cells. Detailed biochemical, proteomic, and knockdown studies indicate that the antiviral activity of SM141 and SM142 results from the dual inhibition of MPro and CatL. Notably, intranasal and intraperitoneal administration of SM141 and SM142 lead to reduced viral replication, viral loads in the lung, and enhanced survival in SARS-CoV2 infected K18-ACE2 transgenic mice. In total, these data indicate that SM141 and SM142 represent promising scaffolds on which to develop antiviral drugs against SARS-CoV2.

Molecular Dynamics Simulation of Ligands from Anredera cordifolia (Binahong) to the Main Protease (M pro) of SARS-CoV-2

COVID-19 in Indonesia is considered to be entering the endemic phase, and the population is expected to live side by side with the SARS-CoV 2 viruses and their variants. In this study, procyanidin, oleic acid, methyl linoleic acid, and vitexin, four compounds from binahong leaves-tropical/subtropical plant, were examined for their interactions with the major protease (Mpro) of the SARS-CoV 2 virus. Molecular dynamics simulation shows that procyanidin and vitexin have the best docking scores of -9.132 and -8.433, respectively. Molecular dynamics simulation also shows that procyanidin and vitexin have the best Root Mean Square Displacement (RMSD) and Root Mean Square Fluctuation (RMSF) performance due to dominant hydrogen, hydrophobic, and water bridge interactions. However, further strain energy calculation obtained from ligand torsion analyses, procyanidin and vitexin do not conform as much as quercetin as ligand control even though these two ligands have good performance in terms of interaction with the target protein.

Gold-Based Metal Drugs as Inhibitors of Coronavirus Proteins: The Inhibition of SARS-CoV-2 Main Protease by Auranofin and Its Analogs

Gold compounds have a long tradition in medicine and offer many opportunities for new therapeutic applications. Herein, we evaluated the lead compound Auranofin and five related gold(I) complexes as possible inhibitors of SARS-CoV-2 Main Protease (SARS-CoV-2 Mpro), a validated drug target for the COVID-19 disease. The investigational panel of gold compounds included Auranofin; three halido analogues, i.e., Au(PEt3)Cl, Au(PEt3)Br, and Au(PEt3)I; and two gold carbene complexes, i.e., Au(NHC)Cl and [Au(NHC)2]PF6. Notably, all these gold compounds, with the only exception of [Au(NHC)2]PF6, turned out to be potent inhibitors of the catalytic activity of SARS-CoV-2 Mpro: the measured Ki values were in the range 2.1-0.4 μM. The reactions of the various gold compounds with SARS-CoV-2 Mpro were subsequently investigated through electrospray ionization (ESI) mass spectrometry (MS) upon a careful optimization of the experimental conditions; the ESI MS spectra provided clear evidence for the formation of tight metallodrug-protein adducts and for the coordination of well defined gold-containing fragments to the SARS-CoV-2 Mpro, again with the only exception of [Au(NHC)2]PF6, The metal-protein stoichiometry was unambiguously determined for the resulting species. The crystal structures of the metallodrug- Mpro adducts were solved in the case of Au(PEt3)Br and Au(NHC)Cl. These crystal structures show that gold coordination occurs at the level of catalytic Cys 145 in the case of Au(NHC)Cl and at the level of both Cys 145 and Cys 156 for Au(PEt3)Br. Tight coordination of gold atoms to functionally relevant cysteine residues is believed to represent the true molecular basis of strong enzyme inhibition.

A molecular evolution algorithm for ligand design in DOCK

As a complement to virtual screening, de novo design of small molecules is an alternative approach for identifying potential drug candidates. Here, we present a new 3D genetic algorithm to evolve molecules through breeding, mutation, fitness pressure, and selection. The method, termed DOCK_GA, builds upon and leverages powerful sampling, scoring, and searching routines previously implemented into DOCK6. Three primary experiments were used during development: Single-molecule evolution evaluated three selection methods (elitism, tournament, and roulette), in four clinically relevant systems, in terms of mutation type and crossover success, chemical properties, ensemble diversity, and fitness convergence, among others. Large scale benchmarking assessed performance across 651 different protein-ligand systems. Ensemble-based evolution demonstrated using multiple inhibitors simultaneously to seed growth in a SARS-CoV-2 target. Key takeaways include: (1) The algorithm is robust as demonstrated by the successful evolution of molecules across a large diverse dataset. (2) Users have flexibility with regards to parent input, selection method, fitness function, and molecular descriptors. (3) The program is straightforward to run and only requires a single executable and input file at run-time. (4) The elitism selection method yields more tightly clustered molecules in terms of 2D/3D similarity, with more favorable fitness, followed by tournament and roulette.

Novel ciprofloxacin and norfloxacin-tetrazole hybrids as potential antibacterial and antiviral agents: targeting S. aureus topoisomerase and SARS-CoV-2-MPro

This study was designed to synthesize hybridizing molecules from ciprofloxacin and norfloxacin by enhancing their biological activity with tetrazoles. The synthesized compounds were investigated in the interaction with the target enzyme of fluoroquinolones (DNA gyrase) and COVID-19 main protease using molecular similarity, molecular docking, and QSAR studies. A QSAR study was carried out to explore the antibacterial activity of our compounds over Staphylococcus aureus a QSAR study, using descriptors obtained from the docking with DNA gyrase, in combination with steric type descriptors, was done obtaining suitable statistical parameters (R2=87.00, QLMO2=71.67, and QEXT2=73.49) to support our results. The binding interaction of our compounds with CoV-2-Mpro was done by molecular docking and were compared with different covalent and non-covalent inhibitors of this enzyme. For the docking studies we used several crystallographic structures of the CoV-2-Mpro. The interaction energy values and binding mode with several key residues, by our compounds, support the capability of them to be CoV-2-Mpro inhibitors. The characterization of the compounds was completed using FT-IR, ¹H-NMR, ¹³C-NMR, ¹⁹F-NMR and HRMS spectroscopic methods. The results showed that compounds 1, 4, 5, 10 and 12 had the potential to be further studied as new antibacterial and antiviral compounds

Computational molecular interaction between SARS-CoV-2 main protease and theaflavin digallate using free energy perturbation and molecular dynamics

Our objective was to identify the molecule which can inhibit SARS-CoV-2 main protease and can be easily procured. Natural products may provide such molecules and can supplement the current custom chemical synthesis-based drug discovery for this objective. A combination of docking approaches, scoring functions, classical molecular dynamic simulation, binding pose metadynamics, and free energy perturbation calculations have been employed in this study. Theaflavin digallate has been observed in top-scoring compounds after the three independent virtual screening simulations of 598435 compounds (unique 27256 chemical entities). The main protease-theaflavin digallate complex interacts with critical active site residues of the main protease in molecular dynamics simulation independent of the explored computational framework, simulation time, initial structure, and force field used. Theaflavin digallate forms approximately three hydrogen bonds with Glutamate166 of main protease, primarily through hydroxyl groups in the benzene ring of benzo(7)annulen-6-one, along with other critical residues. Glu166 is the most critical amino acid for main protease dimerization, which is necessary for catalytic activity. The estimated binding free energy, calculated by Amber and Schrodinger MMGBSA module, reflects a high binding free energy between theaflavin digallate and main protease. Binding pose metadynamics simulation shows the highly persistent H-bond and a stable pose for the theaflavin digallate-main protease complex. Using method control, experimental controls, and test set, alchemical transformation studies confirm high relative binding free energy of theaflavin digallate with the main protease. Computational molecular interaction suggests that theaflavin digallate can inhibit the main protease of SARS-CoV-2.

An open-source molecular builder and free energy preparation workflow

Automated free energy calculations for the prediction of binding free energies of congeneric series of ligands to a protein target are growing in popularity, but building reliable initial binding poses for the ligands is challenging. Here, we introduce the open-source FEgrow workflow for building user-defined congeneric series of ligands in protein binding pockets for input to free energy calculations. For a given ligand core and receptor structure, FEgrow enumerates and optimises the bioactive conformations of the grown functional group(s), making use of hybrid machine learning/molecular mechanics potential energy functions where possible. Low energy structures are optionally scored using the gnina convolutional neural network scoring function, and output for more rigorous protein-ligand binding free energy predictions. We illustrate use of the workflow by building and scoring binding poses for ten congeneric series of ligands bound to targets from a standard, high quality dataset of protein-ligand complexes. Furthermore, we build a set of 13 inhibitors of the SARS-CoV-2 main protease from the literature, and use free energy calculations to retrospectively compute their relative binding free energies. FEgrow is freely available at https://github.com/cole-group/FEgrow, along with a tutorial.

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.

A multi-pronged evaluation of aldehyde-based tripeptidyl main protease inhibitors as SARS-CoV-2 antivirals

As an essential enzyme of SARS-CoV-2, the COVID-19 pathogen, main protease (M^Pro) is a viable target to develop antivirals for the treatment of COVID-19. By varying chemical compositions at both P2 and P3 positions and the N-terminal protection group, we synthesized 18 tripeptidyl M^Pro inhibitors that contained also an aldehyde warhead and β-(S-2-oxopyrrolidin-3-yl)-alaninal at the P1 position. Systematic characterizations of these inhibitors were conducted, including their in vitro enzymatic inhibition potency, X-ray crystal structures of their complexes with M^Pro, their inhibition of M^Pro transiently expressed in 293T cells, and cellular toxicity and SARS-CoV-2 antiviral potency of selected inhibitors. These inhibitors have a large variation of determined in vitro enzymatic inhibition IC₅₀ values that range from 4.8 to 650 nM. The determined in vitro enzymatic inhibition IC₅₀ values reveal that relatively small side chains at both P2 and P3 positions are favorable for achieving high in vitro M^Pro inhibition potency, the P3 position is tolerable toward unnatural amino acids with two alkyl substituents on the α-carbon, and the inhibition potency is sensitive toward the N-terminal protection group. X-ray crystal structures of M^Pro bound with 16 inhibitors were determined. In all structures, the M^Pro active site cysteine interacts covalently with the aldehyde warhead of the bound inhibitor to form a hemithioacetal that takes an S configuration. For all inhibitors, election density around the N-terminal protection group is weak indicating possible flexible binding of this group to M^Pro. In M^Pro, large structural variations were observed on residues N142 and Q189. Unlike their high in vitro enzymatic inhibition potency, most inhibitors showed low potency to inhibit M^Pro that was transiently expressed in 293T cells. Inhibitors that showed high potency to inhibit M^Pro transiently expressed in 293T cells all contain O-tert-butyl-threonine at the P3 position. These inhibitors also exhibited relatively low cytotoxicity and high antiviral potency. Overall, our current and previous studies indicate that O-tert-butyl-threonine at the P3 site is a key component to achieve high cellular and antiviral potency for tripeptidyl aldehyde inhibitors of M^Pro.

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.

Molecular Dynamics Simulations of Protein RNA Complexes by Using an Advanced Electrostatic Model

Protein-RNA interactions are integral to the biological functions of RNA. It is well recognized that molecular dynamics (MD) simulations of protein-RNA complexes are more challenging than those of each component. The difficulty arises from the strong electrostatic interactions and the delicate balance between various types of physical forces at the interface. Previously, MD simulations of protein-RNA complexes have predominantly employed fixed-charge force fields. Although force field modifications have been developed to address problems identified in the simulations, some protein-RNA structures are still hard to reproduce by simulations. Here, we present MD simulations of two representative protein-RNA complexes using the AMOEBA polarizable force field. The van der Waals parameters were refined to reproduce accurate quantum-mechanical data of base-base and base-amino acid interactions. It was found that the refined parameters produced a more stable hydrogen-bond network in the interface. One of the complexes remained stable during the short simulations, whereas it could quickly break down in previous simulations using fixed-charge force fields. There was reversible breaking and formation of hydrogen bonds that are observed in the crystal structure, which may indicate the difference in solution and crystal structures. While further improvement and validation of the force fields are still needed, this work demonstrates that polarizable force fields are promising for the study of protein-RNA complexes.

A reporter cell line for the automated quantification of SARS-CoV-2 infection in living cells

The SARS-CoV-2 pandemic and the urgent need for massive antiviral testing highlighted the lack of a good cell-based assay that allowed for a fast, automated screening of antivirals in high-throughput content with minimal handling requirements in a BSL-3 environment. The present paper describes the construction of a green fluorescent substrate that, upon cleavage by the SARS-CoV-2 main protease, re-localizes from the cytoplasm in non-infected cells to the nucleus in infected cells. The construction was stably expressed, together with a red fluorescent nuclear marker, in a highly susceptible clone derived from Vero-81 cells. With this fluorescent reporter cell line, named F1G-red, SARS-CoV-2 infection can be scored automatically in living cells by comparing the patterns of green and red fluorescence signals acquired by automated confocal microscopy in a 384-well plate format. We show the F1G-red system is sensitive to several SARS-CoV-2 variants of concern and that it can be used to assess antiviral activities of compounds in dose-response experiments. This high-throughput system will provide a reliable tool for antiviral screening against SARS-CoV-2.

Recent Drug Development and Medicinal Chemistry Approaches for the Treatment of SARS-CoV-2 Infection and COVID-19

COVID-19, caused by SARS-CoV-2 infection, continues to be a major public health crisis around the globe. Development of vaccines and the first cluster of antiviral drugs has brought promise and hope for prevention and treatment of severe coronavirus disease. However, continued development of newer, safer, and more effective antiviral drugs are critically important to combat COVID-19 and counter the looming pathogenic variants. Studies of the coronavirus life cycle revealed several important biochemical targets for drug development. In the present review, we focus on recent drug design and medicinal chemistry efforts in small molecule drug discovery, including the development of nirmatrelvir that targets viral protein synthesis and remdesivir and molnupiravir that target viral RdRp. These are recent FDA approved drugs for the treatment of COVID-19.

Structural Basis for the Inhibition of Coronaviral Main Proteases by a Benzothiazole-Based Inhibitor

The ongoing spread of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) has caused hundreds of millions of cases and millions of victims worldwide with serious consequences to global health and economies. Although many vaccines protecting against SARS-CoV-2 are currently available, constantly emerging new variants necessitate the development of alternative strategies for prevention and treatment of COVID-19. Inhibitors that target the main protease (M^pro) of SARS-CoV-2, an essential enzyme that promotes viral maturation, represent a key class of antivirals. Here, we showed that a peptidomimetic compound with benzothiazolyl ketone as warhead, YH-53, is an effective inhibitor of SARS-CoV-2, SARS-CoV, and MERS-CoV M^pros. Crystal structures of M^pros from SARS-CoV-2, SARS-CoV, and MERS-CoV bound to the inhibitor YH-53 revealed a unique ligand-binding site, which provides new insights into the mechanism of inhibition of viral replication. A detailed analysis of these crystal structures defined the key molecular determinants required for inhibition and illustrate the binding mode of M^pros from other coronaviruses. In consideration of the important role of M^pro in developing antivirals against coronaviruses, insights derived from this study should add to the design of pan-coronaviral M^pro inhibitors that are safer and more effective.

Relative binding free energy calculations with transformato: A molecular dynamics engine-independent tool

We present the software package transformato for the setup of large-scale relative binding free energy calculations. Transformato is written in Python as an open source project (https://github.com/wiederm/transformato); in contrast to comparable tools, it is not closely tied to a particular molecular dynamics engine to carry out the underlying simulations. Instead of alchemically transforming a ligand L₁ directly into another L₂, the two ligands are mutated to a common core. Thus, while dummy atoms are required at intermediate states, in particular at the common core state, none are present at the physical endstates. To validate the method, we calculated 76 relative binding free energy differences ΔΔGL1→L2bind for five protein–ligand systems. The overall root mean squared error to experimental binding free energies is 1.17 kcal/mol with a Pearson correlation coefficient of 0.73. For selected cases, we checked that the relative binding free energy differences between pairs of ligands do not depend on the choice of the intermediate common core structure. Additionally, we report results with and without hydrogen mass reweighting. The code currently supports OpenMM, CHARMM, and CHARMM/OpenMM directly. Since the program logic to choose and construct alchemical transformation paths is separated from the generation of input and topology/parameter files, extending transformato to support additional molecular dynamics engines is straightforward.

New Chemicals Suppressing SARS-CoV-2 Replication in Cell Culture

Candidates to being inhibitors of the main protease (Mpro) of SARS-CoV-2 were selected from the database of Voronezh State University using molecular modeling. The database contained approximately 19,000 compounds represented by more than 41,000 ligand conformers. These ligands were docked into Mpro using the SOL docking program. For one thousand ligands with best values of the SOL score, the protein–ligand binding enthalpy was calculated by the PM7 quantum-chemical method with the COSMO solvent model. Using the SOL score and the calculated protein–ligand binding enthalpies, eighteen compounds were selected for the experiments. Several of these inhibitors suppressed the replication of the coronavirus in cell culture, and we used the best three among them in the search for chemical analogs. Selection among analogs using the same procedure followed by experiments led to identification of seven inhibitors of the SARS-CoV-2 replication in cell culture with EC50 values at the micromolar level. The identified inhibitors belong to three chemical classes. The three inhibitors, 4,4-dimethyldithioquinoline derivatives, inhibit SARS-CoV-2 replication in Vero E6 cell culture just as effectively as the best published non-covalent inhibitors, and show low cytotoxicity. These results open up a possibility to develop antiviral drugs against the SARS-CoV-2 coronavirus.

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.

Structural Basis of Main Proteases of Coronavirus Bound to Drug Candidate PF-07304814

New variants of the severe acute respiratory syndrome Coronavirus 2 (SARS-CoV-2) emerged and spread rapidly all over the world, which strongly supports the need for pharmacological options to complement vaccine strategies. Main protease (M^pro or 3CL^pro) is a critical enzyme in the life cycle of SARS-CoV-2 and appears to be highly conserved among different genera of coronaviruses, making it an ideal target for the development of drugs with broad-spectrum property. PF-07304814 developed by Pfizer is an intravenously administered inhibitor targeting SARS-CoV-2 M^pro. Here we showed that PF-07304814 displays broad-spectrum inhibitory activity against M^pros from multiple coronaviruses. Crystal structures of M^pros of SARS-CoV-2, SARS-CoV, MERS-CoV, and HCoV-NL63 bound to the inhibitor PF-07304814 revealed a conserved ligand-binding site, providing new insights into the mechanism of inhibition of viral replication. A detailed analysis of these crystal structures complemented by comprehensive comparison defined the key structural determinants essential for inhibition and illustrated the binding mode of action of M^pros from different coronaviruses. In view of the importance of M^pro for the medications of SARS-CoV-2 infection, insights derived from the present study should accelerate the design of pan-coronaviral main protease inhibitors that are safer and more effective.

Discovery of 2-(furan-2-ylmethylene)hydrazine-1-carbothioamide derivatives as novel inhibitors of SARS-CoV-2 main protease

The COVID-19 posed a serious threat to human life and health, and SARS-CoV-2 M^pro has been considered as an attractive drug target for the treatment of COVID-19. Herein, we report 2-(furan-2-ylmethylene)hydrazine-1-carbothioamide derivatives as novel inhibitors of SARS-CoV-2 M^pro developed by in-house library screening and biological evaluation. Similarity search led to the identification of compound F8–S43 with the enzymatic IC₅₀ value of 10.76 μM. Further structure-based drug design and synthetic optimization uncovered compounds F8–B6 and F8–B22 as novel non-peptidomimetic inhibitors of M^pro with IC₅₀ values of 1.57 μM and 1.55 μM, respectively. Moreover, enzymatic kinetic assay and mass spectrometry demonstrated that F8–B6 was a reversible covalent inhibitor of M^pro. Besides, F8–B6 showed low cytotoxicity with CC₅₀ values of more than 100 μM in Vero and MDCK cells. Overall, these novel SARS-CoV-2 M^pro non-peptidomimetic inhibitors provide a useful starting point for further structural optimization.

Binding of Dual-Function Hybridized MetalOrganic Capsules to Enzymes for Cascade Catalysis

The combination of chemo- and biocatalysis for multistep syntheses provides attractive advantages in terms of evolvability, promiscuity, and sustainability striving for desirable catalytic performance. Through the encapsulation of flavin analogues by both NADH and heme mimics codecorated heteroleptic metal-organic capsules, herein, we report a progressive host-guest strategy to imitate cytochrome P450s catalysis for cascade oxidative coupling catalysis. Besides the construction of stable dual-function metal-organic capsules and the modification of cofactor-decorated capsules at the domain of enzymes, this supramolecular strategy involves multistage directional electron flow, affording reactive ferric peroxide species for inducing oxygenation. Under light irradiation, the metal-organic capsule selectively converts stilbene to oxidative coupling products (including 2-oxo-1,2-diphenylethyl formate, 2-alkoxy-1,2-diphenylethanone) in tandem with enzymatic reactions respectively, at the domain of natural enzymes. The ingenious combination of capsules and enzymes with the in situ-regenerated capsule-loaded NADH cofactor promises non-native coupling reactions by forming regional cooperation and division. This abiotic-biotic conjugated host-guest strategy is conducive to the de novo creation of multifunctional components approaching active enzymatic sites for reinforced matter and energy transporting, demonstrating a key role of multicomponent supramolecular catalysts for one-pot integrated catalytic conversions.

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.

Exploration of chalcones as 3-chymotrypsin-like protease (3CLpro) inhibitors of SARS-CoV-2 using computational approaches

The main protease 3CL^pro is one of the potential targets against coronavirus. Inhibiting this enzyme leads to the interruption of viral replication. Chalcone and its derivatives were reported to possess the ability to bind to 3CL^pro protease in the binding pocket. This study explored an in-house database of 269 chalcones as 3CL^pro inhibitors using in silico screening models, including molecular docking, molecular dynamics simulation, binding free energy calculation, and ADME prediction. C264 and C235 stand out as the two most potential structures. The top hit compound C264 was with the Jamda score of −2.8329 and the MM/GBSA binding energy mean value of −28.23 ± 3.53 kcal/mol, which was lower than the reference ligand. Despite the lower mean binding energy (−22.07 ± 3.39 kcal/mol), in-depth analysis of binding interaction suggested C235 could be another potential candidate. Further, in vitro and in vivo experiments are required to confirm the inhibitory ability.

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.

Primer for Designing Main Protease (Mpro) Inhibitors of SARS-CoV-2

The COVID-19 outbreak has been devastating, with hundreds of millions of infections and millions of deaths reported worldwide. In response, the application of structure-activity relationships (SAR) upon experimentally validated inhibitors of SARS-CoV-2 main protease (Mpro) may provide an avenue for the identification of new lead compounds active against COVID-19. Upon the basis of information gleaned from a combination of reported crystal structures and the docking of experimentally validated inhibitors, four "rules" for designing potent Mpro inhibitors have been proposed. The aim here is to guide medicinal chemists toward the most probable hits and to provide guidance on repurposing available structures as Mpro inhibitors. Experimental examination of our own previously reported inhibitors using the four "rules" identified a potential lead compound, the cathepsin inhibitor GB111-NH2, that was 2.3 times more potent than SARS-CoV-2 Mpro inhibitor N3.

Discovery of 2-Thiobenzimidazoles as Noncovalent Inhibitors of SARS-CoV-2 Main Protease

The discovery of antiviral agents against SARS-CoV-2 is an important step toward ending the COVID-19 pandemic and to tackle future outbreaks. In this context, the main protease (M^pro) represents an ideal target for developing coronavirus antivirals, being conserved among different strains and essential for survival. In this work, using in silico tools, we created and validated a docking protocol able to predict binders to the catalytic site of M^pro. The following structure-based virtual screening of a subset of the ZINC library (over 4.3 million unique structures), led to the identification of a hit compound having a 2-thiobenzimidazole scaffold. The inhibitory activity was confirmed using a FRET-based proteolytic assay against recombinant M^pro. Structure-activity relationships were obtained with the synthesis of a small library of analogs, guided by the analysis of the docking pose. Our efforts led to the identification of a micromolar M^pro inhibitor (IC₅₀ = 14.9 µM) with an original scaffold possessing ideal drug-like properties (predicted using the QikProp function) and representing a promising lead for the development of a novel class of coronavirus antivirals.

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.

Identifying Possible AChE Inhibitors from Drug-like Molecules via Machine Learning and Experimental Studies

Acetylcholinesterase (AChE) is one of the most important drug targets for Alzheimer's disease (AD) treatment. In this work, a machine learning model was trained to rapidly and accurately screen large chemical databases for the potential inhibitors of AChE. The obtained results were then validated via in vitro enzyme assay. Moreover, atomistic simulations including molecular docking and molecular dynamics simulations were then used to understand molecular insights into the binding process of ligands to AChE. In particular, two compounds including benzyl trifluoromethyl ketone and trifluoromethylstyryl ketone were indicated as highly potent inhibitors of AChE because they established IC₅₀ values of 0.51 and 0.33 μM, respectively. The obtained IC₅₀ of two compounds is significantly lower than that of galantamine (2.10 μM). The predicted log(BB) suggests that the compounds may be able to traverse the blood-brain barrier. A good agreement between computational and experimental studies was observed, indicating that the hybrid approach can enhance AD therapy.

Defining the substrate envelope of SARS-CoV-2 main protease to predict and avoid drug resistance

Coronaviruses can evolve and spread rapidly to cause severe disease morbidity and mortality, as exemplified by SARS-CoV-2 variants of the COVID-19 pandemic. Although currently available vaccines remain mostly effective against SARS-CoV-2 variants, additional treatment strategies are needed. Inhibitors that target essential viral enzymes, such as proteases and polymerases, represent key classes of antivirals. However, clinical use of antiviral therapies inevitably leads to emergence of drug resistance. In this study we implemented a strategy to pre-emptively address drug resistance to protease inhibitors targeting the main protease (Mpro) of SARS-CoV-2, an essential enzyme that promotes viral maturation. We solved nine high-resolution cocrystal structures of SARS-CoV-2 Mpro bound to substrate peptides and six structures with cleavage products. These structures enabled us to define the substrate envelope of Mpro, map the critical recognition elements, and identify evolutionarily vulnerable sites that may be susceptible to resistance mutations that would compromise binding of the newly developed Mpro inhibitors. Our results suggest strategies for developing robust inhibitors against SARS-CoV-2 that will retain longer-lasting efficacy against this evolving viral pathogen.

Production of a functionally active recombinant SARS-CoV-2 (COVID-19) 3C-Like protease and a soluble inactive 3C-like protease-RBD chimeric in a prokaryotic expression system

During the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) intracellular life-cycle, two large polyproteins, pp1a and pp1ab, are produced. Processing of these by viral cysteine proteases, the papain-like protease (PLpro) and the chymotrypsin-like 3C-like protease (3CL-pro) release non-structural proteins necessary for the establishment of the viral replication and transcription complex (RTC), crucial for viral replication. Hence, these proteases are considered prime targets against which anti-coronavirus disease 2019 (COVID-19) drugs could be developed. Here, we describe the expression of a highly soluble and functionally active recombinant 3CL-pro using Escherichia coli BL21 cells. We show that the enzyme functions in a dimeric form and exhibits an unexpected inhibitory profile because its activity is potently blocked by serine rather than cysteine protease inhibitors. In addition, we assessed the ability of our 3CL-pro to function as a carrier for the receptor binding domain (RBD) of the Spike protein. The co-expressed chimeric protein, 3CLpro-RBD, did not exhibit 3CL-pro activity, but its enhanced solubility made purification easier and improved RBD antigenicity when tested against serum from vaccinated individuals in ELISAs. Chimeric proteins containing the 3CL-pro could represent an innovative approach to developing new COVID-19 vaccines.

A Concise Synthetic Method for Constructing 3-Substituted Piperazine-2-Acetic Acid Esters from 1,2-Diamines

We report a short synthetic route for synthesizing 2,3-substituted piperazine acetic acid esters. Optically pure amino acids were efficiently converted into 1,2-diamines that could be utilized to deliver the title 2,3-substituted piperazines in five steps with a high enantiomeric purity. The novel route facilitated, for the first time, the synthesis of 3-phenyl substituted-2-piperazine acetic acid esters that were difficult to achieve using other methods; however, in this case, the products underwent racemization.