pH-dependent conformational flexibility of the SARS-CoV main proteinase (M(pro)) dimer: molecular dynamics simulations and multiple X-ray structure analyses

Competitive Interaction of the SGFRKMAF Peptide with 3CLpro Dimerization Intermediates: A Brownian Dynamics Investigation

The SGFRKMAF peptide is known to inhibit the dimerization of 3CLpro monomers, which is essential for SARS-CoV-2 replication. The mechanism behind this, however, is largely unknown. In this work, we used Brownian dynamics simulations to compare and contrast 3CLpro monomer-monomer interactions and 3CLpro monomer-SGFRKMAF peptide interactions. We found that formation of the 3CLpro wild-type dimer could potentially involve formation of three intermediates that are primarily stabilized by G11-G124, S1-S301, and T118-G278 interactions. Analysis of 3CLpro monomer interaction with the SGFRKMAF peptide, however, revealed the presence of eight basins of interactions where the peptide assumes the highest local densities at the 3CLPro monomer surface. The second highest-density basin was found to coincide with the interface region of the wild-type 3CLpro dimer, thereby directly blocking the 3CLpro dimer-dimer interactions. The other basins, however, were found to lie far from the interface region. Notably, we found that only 6% of the BD trajectories end up directly into the basin at the interface region and ∼39% of the trajectories end up into those basins lying away from the interface region, indicating a greater role for peptide binding at sites away from the dimer interface region. Importantly, the locations of the basins lying away from the interface were found to coincide with the 3CLpro residues involved in stabilization of the 3CLpro monomer-monomer intermediates. Given that the rate constant of the peptide reaching the monomer surface was found to be almost an order of magnitude higher than the rate constant of monomer-monomer association, the SGFRKMAF peptide has the potential to inhibit dimerization of 3CLpro monomers not only through blocking the interface region but also through blocking the formation of the intermediates involved in the dimerization process. This could potentially open new avenues for 3CLpro dimerization inhibitors that transcend traditional X-ray-based discovery approaches.

On Inactivation of the Coronavirus Main Protease

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

Improved expression and purification of highly-active 3 chymotrypsin-like protease from SARS-CoV-2

Severe acute respiratory syndrome-coronavirus 2 (SARS-CoV-2) is the causative pathogen of coronavirus disease-19 (COVID-19). The COVID-19 pandemic has resulted in millions of deaths and widespread socio-economic damage worldwide. Therefore, numerous studies have been conducted to identify effective measures to control the spreading of the virus. Among various potential targets, the 3 chymotrypsin-like protease (3CLpro), also known as Mpro, stands out as the key protease of SARS-CoV-2, playing an essential role in virus replication and assembly, is the most prospective. In this study, we modified the commercial vector, pETM33-Nsp5-Mpro (plasmid # 156475, Addgene, USA), by inserting an autocleavage site (AVLQ) of 3CLpro and 6 × His-tag encoding sequences before and after the Nsp5-Mpro sequence, respectively. This modification enabled the expression of 3CLpro as an authentic N terminal protease (au3CLpro), which was purified to electrophoretic homogeneity by a single-step chromatography using two tandem Glutathione- and Ni-Sepharose columns. The enzyme au3CLpro demonstrated significantly higher activity (3169 RFU/min/μg protein) and catalytic efficiency (Kcat/Km of 0.007 μM-1.s-1) than that of the 3CLpro (com3CLpro) expressed from the commercial vector (pETM33-Nsp5-Mpro) with specific activity 889 RFU/min/μg and Kcat/Km of 0.0015 μM-1.s-1, respectively. Optimal conditions for au3CLpro activity included a 50 mM Tris-HCl buffer at pH 7, containing 150 mM NaCl and 0.1 mg/ml BSA at 37 °C.

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

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

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

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

Analyzing 3D structures of the SARS-CoV-2 main protease reveals structural features of ligand binding for COVID-19 drug discovery

The severe acute respiratory syndrome-coronavirus 2 (SARS-CoV-2) main protease has an essential role in viral replication and has become a major target for coronavirus 2019 (COVID-19) drug development. Various inhibitors have been discovered or designed to bind to the main protease. The availability of more than 550 3D structures of the main protease provides a wealth of structural details on the main protease in both ligand-free and ligand-bound states. Therefore, we examined these structures to ascertain the structural features for the role of the main protease in the cleavage of polyproteins, the alternative conformations during main protease maturation, and ligand interactions in the main protease. The structural features unearthed could promote the development of COVID-19 drugs targeting the SARS-CoV-2 main protease.

Antiviral activity of Cenostigma pluviosum var. peltophoroides extract and fractions against SARS-CoV-2

Few extracts of plant species from the Brazilian flora have been validated from a pharmacological and clinical point of view, and it is important to determine whether their traditional use is proven by pharmacological effects. Cenostigma pluviosum var. peltophoroides is one of those plants, which belongs to the Fabaceae family that is widely used in traditional medicine and is very rich in tannins. Due to the lack of effective drugs to treat severe cases of Covid-19, the main protease of SARS-CoV-2 (Mpro) becomes an attractive target in the research for new antivirals since this enzyme is crucial for virus replication and does not have homologs in humans. This study aimed to prospect inhibitor candidates among the compounds from C. pluviosum extract, by virtual screening simulations using SARS-CoV-2 Mpro as target. Experimental validation was made by inhibitory proteolytic assays of recombinant Mpro and by antiviral activity with infected Vero cells. Docking simulations identify four compounds with potential inhibitory activity of Mpro present in the extract. The compound pentagalloylglucose showed the best result in proteolytic kinetics experiments, with suppression of recombinant Mpro activity by approximately 60%. However, in experiments with infected cells ethyl acetate fraction and sub-fractions, F2 and F4 of C. pluviosum extract performed better than pentagalloylglucose, reaching close to 100% of antiviral activity. The prominent activity of the extract fractions in infected cells may be a result of a synergistic effect from the different hydrolyzable tannins present, performing simultaneous action on Mpro and other targets from SARS-CoV-2 and host.Communicated by Ramaswamy H. Sarma.

P1 Glutamine isosteres in the design of inhibitors of 3C/3CL protease of human viruses of the Pisoniviricetes class

Viral infections are one of the leading causes of acute morbidity in humans and much endeavour has been made by the synthetic community for the development of drugs to treat associated diseases. Peptide-based enzyme inhibitors, usually short sequences of three or four residues, are one of the classes of compounds currently under development for enhancement of their activity and pharmaceutical properties. This review reports the advances made in the design of inhibitors targeting the family of highly conserved viral proteases 3C/3CLpro, which play a key role in viral replication and present minimal homology with mammalian proteases. Particular focus is put on the reported development of P1 glutamine isosteres to generate potent inhibitors mimicking the natural substrate sequence at the site of recognition.'

Systematic Analyses of the Resistance Potential of Drugs Targeting SARS-CoV-2 Main Protease

Drugs that target the main protease (Mpro) of SARS-CoV-2 are effective therapeutics that have entered clinical use. Wide-scale use of these drugs will apply selection pressure for the evolution of resistance mutations. To understand resistance potential in Mpro, we performed comprehensive surveys of amino acid changes that can cause resistance to nirmatrelvir (Pfizer), and ensitrelvir (Xocova) in a yeast screen. We identified 142 resistance mutations for nirmatrelvir and 177 for ensitrelvir, many of which have not been previously reported. Ninety-nine mutations caused apparent resistance to both inhibitors, suggesting likelihood for the evolution of cross-resistance. The mutation with the strongest drug resistance score against nirmatrelvir in our study (E166V) was the most impactful resistance mutation recently reported in multiple viral passaging studies. Many mutations that exhibited inhibitor-specific resistance were consistent with the distinct interactions of each inhibitor in the substrate binding site. In addition, mutants with strong drug resistance scores tended to have reduced function. Our results indicate that strong pressure from nirmatrelvir or ensitrelvir will select for multiple distinct-resistant lineages that will include both primary resistance mutations that weaken interactions with drug while decreasing enzyme function and compensatory mutations that increase enzyme activity. The comprehensive identification of resistance mutations enables the design of inhibitors with reduced potential of developing resistance and aids in the surveillance of drug resistance in circulating viral populations.

Integrative Approach to Dissect the Drug Resistance Mechanism of the H172Y Mutation of SARS-CoV-2 Main Protease

Nirmatrelvir is an orally available inhibitor of SARS-CoV-2 main protease (Mpro) and the main ingredient of Paxlovid, a drug approved by the U.S. Food and Drug Administration for high-risk COVID-19 patients. Recently, a rare natural mutation, H172Y, was found to significantly reduce nirmatrelvir's inhibitory activity. As the COVID-19 cases skyrocket in China and the selective pressure of antiviral therapy builds in the US, there is an urgent need to characterize and understand how the H172Y mutation confers drug resistance. Here, we investigated the H172Y Mpro's conformational dynamics, folding stability, catalytic efficiency, and inhibitory activity using all-atom constant pH and fixed-charge molecular dynamics simulations, alchemical and empirical free energy calculations, artificial neural networks, and biochemical experiments. Our data suggest that the mutation significantly weakens the S1 pocket interactions with the N-terminus and perturbs the conformation of the oxyanion loop, leading to a decrease in the thermal stability and catalytic efficiency. Importantly, the perturbed S1 pocket dynamics weaken the nirmatrelvir binding in the P1 position, which explains the decreased inhibitory activity of nirmatrelvir. Our work demonstrates the predictive power of the combined simulation and artificial intelligence approaches, and together with biochemical experiments, they can be used to actively surveil continually emerging mutations of SARS-CoV-2 Mpro and assist the optimization of antiviral drugs. The presented approach, in general, can be applied to characterize mutation effects on any protein drug targets.

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

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

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

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

Natural Product-Based Screening for Lead Compounds Targeting SARS CoV-2 Mpro

Drugs that cure COVID-19 have been marketed; however, this disease continues to ravage the world without becoming extinct, and thus, drug discoveries are still relevant. Since Mpro has known advantages as a drug target, such as the conserved nature of the active site and the absence of homologous proteins in the body, it receives the attention of many researchers. Meanwhile, the role of traditional Chinese medicine (TCM) in the control of epidemics in China has also led to a focus on natural products, with the hope of finding some promising lead molecules through screening. In this study, we selected a commercial library of 2526 natural products from plants, animals and microorganisms with known biological activity for drug discovery, which had previously been reported for compound screening of the SARS CoV-2 S protein, but had not been tested on Mpro. This library contains compounds from a variety of Chinese herbs, including Lonicerae Japonicae Flos, Forsythiae Fructus and Scutellariae Radix, which are derived from traditional Chinese medicine prescriptions that have been shown to be effective against COVID-19. We used the conventional FRET method for the initial screening. After two rounds of selection, the remaining 86 compounds were divided into flavonoids, lipids, phenylpropanoids, phenols, quinones, alkaloids, terpenoids and steroids according to the skeleton structures, with inhibition rates greater than 70%. The top compounds in each group were selected to test the effective concentration ranges; the IC50 values were as follows: (-)-gallocatechin gallate (1.522 ± 0.126 μM), ginkgolic acid C15:1 (9.352 ± 0.531 μM), hematoxylin (1.025 ± 0.042 μM), fraxetin (2.486 ± 0.178 μM), wedelolactone (1.003 ± 0.238 μM), hydroxytyrosol acetate (3.850 ± 0.576 μM), vanitiolide (2.837 ± 0.225 μM), β,β-dimethylacrylalkannin (2.731 ± 0.308 μM), melanin (7.373 ± 0.368 μM) and cholesteryl sodium sulfate (2.741 ± 0.234μM). In the next step, we employed two biophysical techniques, SPR and nanoDSF, to obtain KD/Kobs values: hematoxylin (0.7 μM), (-)-gallocatechin gallate (126 μM), ginkgolic acid C15:1 (227 μM), wedelolactone (0.9770 μM), β,β-dimethylacrylalkannin (1.9004 μM,), cholesteryl sodium sulfate (7.5950 μM) and melanin (11.5667 μM), which allowed better assessments of the binding levels. Here, seven compounds were the winners. Then, molecular docking experiments were specially performed by AutoDock Vina to analyze the mode of interactions within Mpro and ligands. We finally formulated the present in silico study to predict pharmacokinetic parameters as well as drug-like properties, which is presumably the step that tells humans whether the compounds are drug-like or not. Moreover, hematoxylin, melanin, wedelolactone, β,β-dimethylacrylalkannin and cholesteryl sodium sulfate are in full compliance with the "Lipinski" principle and possess reasonable ADME/T properties, they have a greater potential of being lead compounds. The proposed five compounds are also the first to be found to have potential inhibitory effects on SARS CoV-2 Mpro. We hope that the results in this manuscript may serve as benchmarks for the above potentials.

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.

SARS-CoV-2 proteases Mpro and PLpro: Design of inhibitors with predicted high potency and low mammalian toxicity using artificial neural networks, ligand-protein docking, molecular dynamics simulations, and ADMET calculations

The main (Mpro) and papain-like (PLpro) proteases are highly conserved viral proteins essential for replication of the COVID-19 virus, SARS-COV-2. Therefore, a logical plan for producing new drugs against this pathogen is to discover inhibitors of these enzymes. Accordingly, the goal of the present work was to devise a computational approach to design, characterize, and select compounds predicted to be potent dual inhibitors - effective against both Mpro and PLpro. The first step employed LigDream, an artificial neural network, to create a virtual ligand library. Ligands with computed ADMET profiles indicating drug-like properties and low mammalian toxicity were selected for further study. Initial docking of these ligands into the active sites of Mpro and PLpro was done with GOLD, and the highest-scoring ligands were redocked with AutoDock Vina to determine binding free energies (ΔG). Compounds 89-00, 89-07, 89-32, and 89-38 exhibited favorable ΔG values for Mpro (-7.6 to -8.7 kcal/mol) and PLpro (-9.1 to -9.7 kcal/mol). Global docking of selected compounds with the Mpro dimer identified prospective allosteric inhibitors 89-00, 89-27, and 89-40 (ΔG -8.2 to -8.9 kcal/mol). Molecular dynamics simulations performed on Mpro and PLpro active site complexes with the four top-scoring ligands from Vina demonstrated that the most stable complexes were formed with compounds 89-32 and 89-38. Overall, the present computational strategy generated new compounds with predicted drug-like characteristics, low mammalian toxicity, and high inhibitory potencies against both target proteases to form stable complexes. Further preclinical studies will be required to validate the in silico findings before the lead compounds could be considered for clinical trials.

Disruption of 3CLpro protease self-association by short peptides as a potential route to broad spectrum coronavirus inhibitors

Coronaviruses have posed a persistent threat to human health over the last two decades. Despite the accumulated knowledge about coronavirus-related pathogens, development of an effective treatment for its new variant COVID-19 is highly challenging. For the highly-conserved and main coronavirus protease 3CLpro, dimerization is known to be essential for its catalytic activity and thereby for virus proliferation. Here, we assess the potential of short peptide segments to disrupt dimerization of the 3CLpro protease as a route to block COVID-19 proliferation. Based on the X-ray structure of the 3CLpro dimer, we identified the SPSGVY₁₂₆QCAMRP dodecapeptide segment as overlapping the hotspot regions on the 3CLpro dimer interface. Using computational blind docking of the peptide to the 3CLpro monomer, we found that the SPSGVY₁₂₆QCAMRP peptide has favourable thermodynamic binding (ΔG= -5.93 kcal/mol) to the hotspot regions at the 3CLpro dimer interface. Importantly, the peptide was also found to preferentially bind to the hotspot regions compared to other potential binding sites lying away from the dimer interface (ΔΔG=-1.31 kcal/mol). Docking of peptides corresponding to systematic mutation of the V125 and Y126 residues led to the identification of seven peptides, SPSGHAQCAMRP, SPSGVTQCAMRP, SPSGKPQCAMRP, SPSGATQCAMRP, SPSGWLQCAMRP, SPSGAPQCAMRP and SPSGHPQCAMRP, that outperform the wild-type SPSGVY₁₂₆QCAMRP peptide in terms of preferential binding to the 3CLpro dimer interface. These peptides have the potential to disrupt 3CLpro dimerization and therefore could provide lead structures for the development of broad spectrum COVID-19 inhibitors.Communicated by Ramaswamy H. Sarma.

Nirmatrelvir-resistant SARS-CoV-2 variants with high fitness in an infectious cell culture system

The oral protease inhibitor nirmatrelvir is of key importance for prevention of severe coronavirus disease 2019 (COVID-19). To facilitate resistance monitoring, we studied severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) escape from nirmatrelvir in cell culture. Resistant variants harbored combinations of substitutions in the SARS-CoV-2 main protease (Mpro). Reverse genetics revealed that E166V and L50F + E166V conferred high resistance in infectious culture, replicon, and Mpro systems. While L50F, E166V, and L50F + E166V decreased replication and Mpro activity, L50F and L50F + E166V variants had high fitness in the infectious system. Naturally occurring L50F compensated for fitness cost of E166V and promoted viral escape. Molecular dynamics simulations revealed that E166V and L50F + E166V weakened nirmatrelvir-Mpro binding. Polymerase inhibitor remdesivir and monoclonal antibody bebtelovimab retained activity against nirmatrelvir-resistant variants, and combination with nirmatrelvir enhanced treatment efficacy compared to individual compounds. These findings have implications for monitoring and ensuring treatments with efficacy against SARS-CoV-2 and emerging sarbecoviruses.

Nirmatrelvir-resistant SARS-CoV-2 variants with high fitness in an infectious cell culture system

The oral protease inhibitor nirmatrelvir is of key importance for prevention of severe coronavirus disease 2019 (COVID-19). To facilitate resistance monitoring, we studied severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) escape from nirmatrelvir in cell culture. Resistant variants harbored combinations of substitutions in the SARS-CoV-2 main protease (Mpro). Reverse genetics revealed that E166V and L50F + E166V conferred high resistance in infectious culture, replicon, and Mpro systems. While L50F, E166V, and L50F + E166V decreased replication and Mpro activity, L50F and L50F + E166V variants had high fitness in the infectious system. Naturally occurring L50F compensated for fitness cost of E166V and promoted viral escape. Molecular dynamics simulations revealed that E166V and L50F + E166V weakened nirmatrelvir-Mpro binding. Polymerase inhibitor remdesivir and monoclonal antibody bebtelovimab retained activity against nirmatrelvir-resistant variants, and combination with nirmatrelvir enhanced treatment efficacy compared to individual compounds. These findings have implications for monitoring and ensuring treatments with efficacy against SARS-CoV-2 and emerging sarbecoviruses.

GPU-Accelerated All-Atom Particle-Mesh Ewald Continuous Constant pH Molecular Dynamics in Amber

Constant pH molecular dynamics (MD) simulations sample protonation states on the fly according to the conformational environment and user specified pH conditions; however, the current accuracy is limited due to the use of implicit-solvent models or a hybrid solvent scheme. Here, we report the first GPU-accelerated implementation, parametrization, and validation of the all-atom continuous constant pH MD (CpHMD) method with particle-mesh Ewald (PME) electrostatics in the Amber22 pmemd.cuda engine. The titration parameters for Asp, Glu, His, Cys, and Lys were derived for the CHARMM c22 and Amber ff14sb and ff19sb force fields. We then evaluated the PME-CpHMD method using the asynchronous pH replica-exchange titration simulations with the c22 force field for six benchmark proteins, including BBL, hen egg white lysozyme (HEWL), staphylococcal nuclease (SNase), thioredoxin, ribonuclease A (RNaseA), and human muscle creatine kinase (HMCK). The root-mean-square deviation from the experimental pK_a's of Asp, Glu, His, and Cys is 0.76 pH units, and the Pearson's correlation coefficient for the pK_a shifts with respect to model values is 0.80. We demonstrated that a finite-size correction or much enlarged simulation box size can remove a systematic error of the calculated pK_a's and improve agreement with experiment. Importantly, the simulations captured the relevant biology in several challenging cases, e.g., the titration order of the catalytic dyad Glu35/Asp52 in HEWL and the coupled residues Asp19/Asp21 in SNase, the large pK_a upshift of the deeply buried catalytic Asp26 in thioredoxin, and the large pK_a downshift of the deeply buried catalytic Cys283 in HMCK. We anticipate that PME-CpHMD will offer proper pH control to improve the accuracies of MD simulations and enable mechanistic studies of proton-coupled dynamical processes that are ubiquitous in biology but remain poorly understood due to the lack of experimental tools and limitation of current MD simulations.

pH profiles of 3-chymotrypsin-like protease (3CLpro) from SARS-CoV-2 elucidate its catalytic mechanism and a histidine residue critical for activity

3-Chymotrypsin-like protease (3CLpro) is a promising drug target for coronavirus disease 2019 and related coronavirus diseases because of the essential role of this protease in processing viral polyproteins after infection. Understanding the detailed catalytic mechanism of 3CLpro is essential for designing effective inhibitors of infection by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Molecular dynamics studies have suggested pH-dependent conformational changes of 3CLpro, but experimental pH profiles of SARS-CoV-2 3CLpro and analyses of the conserved active-site histidine residues have not been reported. In this work, pH-dependence studies of the kinetic parameters of SARS-CoV-2 3CLpro revealed a bell-shaped pH profile with 2 pK_a values (6.9 ± 0.1 and 9.4 ± 0.1) attributable to ionization of the catalytic dyad His41 and Cys145, respectively. Our investigation of the roles of conserved active-site histidines showed that different amino acid substitutions of His163 produced inactive enzymes, indicating a key role of His163 in maintaining catalytically active SARS-CoV-2 3CLpro. By contrast, the H164A and H172A mutants retained 75% and 26% of the activity of WT, respectively. The alternative amino acid substitutions H172K and H172R did not recover the enzymatic activity, whereas H172Y restored activity to a level similar to that of the WT enzyme. The pH profiles of H164A, H172A, and H172Y were similar to those of the WT enzyme, with comparable pK_a values for the catalytic dyad. Taken together, the experimental data support a general base mechanism of SARS-CoV-2 3CLpro and indicate that the neutral states of the catalytic dyad and active-site histidine residues are required for maximum enzyme activity.

Functional dynamics of SARS-CoV-2 3C-like protease as a member of clan PA

SARS-CoV-2 3C-like protease (3CL^pro), a potential therapeutic target for COVID-19, consists of a chymotrypsin fold and a C-terminal α-helical domain (domain III), the latter of which mediates dimerization required for catalytic activation. To gain further understanding of the functional dynamics of SARS-CoV-2 3CL^pro, this review extends the scope to the comparative study of many crystal structures of proteases having the chymotrypsin fold (clan PA of the MEROPS database). First, the close correspondence between the zymogen-enzyme transformation in chymotrypsin and the allosteric dimerization activation in SARS-CoV-2 3CL^pro is illustrated. Then, it is shown that the 3C-like proteases of family Coronaviridae (the protease family C30), which are closely related to SARS-CoV-2 3CL^pro, have the same homodimeric structure and common activation mechanism via domain III mediated dimerization. The survey extended to order Nidovirales reveals that all 3C-like proteases belonging to Nidovirales have domain III, but with various chain lengths, and 3CL^pro of family Mesoniviridae (family C107) has the same homodimeric structure as that of C30, even though they have no sequence similarity. As a reference, monomeric 3C proteases belonging to the more distant family Picornaviridae (family C3) lacking domain III are compared with C30, and it is shown that the 3C proteases are rigid enough to maintain their structures in the active state.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12551-022-01020-x.

Functional map of SARS-CoV-2 3CL protease reveals tolerant and immutable sites

The SARS-CoV-2 3CL protease (3CLpro) is an attractive therapeutic target, as it is essential to the virus and highly conserved among coronaviruses. However, our current understanding of its tolerance to mutations is limited. Here, we develop a yeast-based deep mutational scanning approach to systematically profile the activity of all possible single mutants of the 3CLpro and validate a subset of our results within authentic viruses. We reveal that the 3CLpro is highly malleable and is capable of tolerating mutations throughout the protein. Yet, we also identify specific residues that appear immutable, suggesting that these may be targets for future 3CLpro inhibitors. Finally, we utilize our screening as a basis to identify E166V as a resistance-conferring mutation against the clinically used 3CLpro inhibitor, nirmatrelvir. Collectively, the functional map presented herein may serve as a guide to better understand the biological properties of the 3CLpro and for drug development against coronaviruses.

High-throughput screening of SARS-CoV-2 main and papain-like protease inhibitors

The global COVID-19 coronavirus pandemic has infected over 109 million people, leading to over 2 million deaths up to date and still lacking of effective drugs for patient treatment. Here, we screened about 1.8 million small molecules against the main protease (M^pro) and papain like protease (PL^pro), two major proteases in severe acute respiratory syndrome-coronavirus 2 genome, and identified 1851M^pro inhibitors and 205 PL^pro inhibitors with low nmol/l activity of the best hits. Among these inhibitors, eight small molecules showed dual inhibition effects on both M^pro and PL^pro, exhibiting potential as better candidates for COVID-19 treatment. The best inhibitors of each protease were tested in antiviral assay, with over 40% of M^pro inhibitors and over 20% of PL^pro inhibitors showing high potency in viral inhibition with low cytotoxicity. The X-ray crystal structure of SARS-CoV-2 M^pro in complex with its potent inhibitor 4a was determined at 1.8 Å resolution. Together with docking assays, our results provide a comprehensive resource for future research on anti-SARS-CoV-2 drug development.

Electrostatic Features for the Receptor Binding Domain of SARS-COV-2 Wildtype and Its Variants. Compass to the Severity of the Future Variants with the Charge-Rule

Electrostatic intermolecular interactions are important in many aspects of biology. We have studied the main electrostatic features involved in the interaction of the receptor-binding domain (RBD) of the SARS-CoV-2 spike protein with the human receptor Angiotensin-converting enzyme 2 (ACE2). As the principal computational tool, we have used the FORTE approach, capable to model proton fluctuations and computing free energies for a very large number of protein-protein systems under different physical-chemical conditions, here focusing on the RBD-ACE2 interactions. Both the wild-type and all critical variants are included in this study. From our large ensemble of extensive simulations, we obtain, as a function of pH, the binding affinities, charges of the proteins, their charge regulation capacities, and their dipole moments. In addition, we have calculated the pK_as for all ionizable residues and mapped the electrostatic coupling between them. We are able to present a simple predictor for the RBD-ACE2 binding based on the data obtained for Alpha, Beta, Gamma, Delta, and Omicron variants, as a linear correlation between the total charge of the RBD and the corresponding binding affinity. This "RBD charge rule" should work as a quick test of the degree of severity of the coming SARS-CoV-2 variants in the future.

A Guide to the Continuous Constant pH Molecular Dynamics Methods in Amber and CHARMM [Article v1.0]

Like temperature and pressure, solution pH is an important environmental variable in biomolecular simulations. Virtually all proteins depend on pH to maintain their structure and function. In conventional molecular dynamics (MD) simulations of proteins, pH is implicitly accounted for by assigning and fixing protonation states of titratable sidechains. This is a significant limitation, as the assigned protonation states may be wrong and they may change during dynamics. In this tutorial, we guide the reader in learning and using the various continuous constant pH MD methods in Amber and CHARMM packages, which have been applied to predict pK a values and elucidate proton-coupled conformational dynamics of a variety of proteins including enzymes and membrane transporters.

Novel covalent and non-covalent complex-based pharmacophore models of SARS-CoV-2 main protease (Mpro) elucidated by microsecond MD simulations

As the world enters its second year of the pandemic caused by SARS-CoV-2, intense efforts have been directed to develop an effective diagnosis, prevention, and treatment strategies. One promising drug target to design COVID-19 treatments is the SARS-CoV-2 M^pro. To date, a comparative understanding of M^pro dynamic stereoelectronic interactions with either covalent or non-covalent inhibitors (depending on their interaction with a pocket called S1' or oxyanion hole) has not been still achieved. In this study, we seek to fill this knowledge gap using a cascade in silico protocol of docking, molecular dynamics simulations, and MM/PBSA in order to elucidate pharmacophore models for both types of inhibitors. After docking and MD analysis, a set of complex-based pharmacophore models was elucidated for covalent and non-covalent categories making use of the residue bonding point feature. The highest ranked models exhibited ROC-AUC values of 0.93 and 0.73, respectively for each category. Interestingly, we observed that the active site region of M^pro protein-ligand complex undergoes large conformational changes, especially within the S2 and S4 subsites. The results reported in this article may be helpful in virtual screening (VS) campaigns to guide the design and discovery of novel small-molecule therapeutic agents against SARS-CoV-2 M^pro protein.

Insighting the Therapeutic Potential of Fifty (50) Shogaol Derivatives Against Mpro of SARS-CoV-2

SARS-CoV-2, which causes COVID-19 disease, has proven to be a disastrous pandemic due to its contagious nature. This study has been planned to theoretically explore some antidotes against this virus from natural compounds. A total of 150 compounds from the shogaol class and shogaol derivatives (SDs) have been screened whereas 50 among those, which obeyed Lipinski’s Rule of Five (Ro5), have further been investigated using molecular docking techniques. Furthermore, reference antiviral drug chloroquine (ChQ) and Co-Crystallized inhibitor have also been studied against Mpro of SARS-CoV-2 for comparing the potential of our docked ligands. Surprisingly, 78% of our docked ligands have shown binding energies and inhibition constants lower than ChQ and all ligands showed these values lower than an inhibitor. We further visualized the nature of intermolecular interactions for the best docked six ligands, which have shown higher binding affinities. We have also assessed ADMET properties for three ligands that displayed visually the best intermolecular interactions. Quantum analysis of three selected ligands L4, L5, and L9 has proved their reactivity and kinetic stability. Moreover, molecular dynamic simulations over 60[Formula: see text]ns have been run for free Mpro and its selected three ligand-protein complexes for evaluating conformational stability and residual flexibility of docked complexes. Furthermore, 100[Formula: see text]ns the MD simulations have been performed for two ligand complexes L4, L5 (with negative binding free energy), and inhibitor. Available parameters suggest stable complexes for our ligands and could be active drugs against SARS-CoV-2 in near future.

Effect of pH on stability of dimer structure of the main protease of coronavirus-2

The viral main protease (M^pro) from a novel severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) is a key enzyme essential for viral replication and has become an attractive target for antiviral drug development. The M^pro forms a functional dimer and exhibits a pH-dependent enzyme activity and dimerization. Here, we report a molecular dynamics (MD) investigation to gain insights into the structural stability of the enzyme dimer at neutral and acidic pH. Our data shows larger changes in structure of the protein with the acidic pH than that with the neutral pH. Structural analysis of MD trajectories reveals a substantial increase in intersubunit separation, the loss of domain contacts, binding free energy and interaction energy of the dimer which implies the protein instability and tendency of dimer dissociation at acidic pH. The loss in the interaction energy is mainly driven by electrostatic interactions. We have identified the intersubunit hydrogen-bonding residues involved in the decreased dimer stability. These findings may be helpful for rational drug design and target evaluation against COVID-19.

The Functional Landscape of SARS-CoV-2 3CL Protease

SARS-CoV-2 (severe acute respiratory syndrome coronavirus 2) as the etiologic agent of COVID-19 (coronavirus disease 2019) has drastically altered life globally. Numerous efforts have been placed on the development of therapeutics to treat SARS-CoV-2 infection. One particular target is the 3CL protease (3CL ^pro ), which holds promise as it is essential to the virus and highly conserved among coronaviruses, suggesting that it may be possible to find broad inhibitors that treat not just SARS-CoV-2 but other coronavirus infections as well. While the 3CL protease has been studied by many groups for SARS-CoV-2 and other coronaviruses, our understanding of its tolerance to mutations is limited, knowledge which is particularly important as 3CL protease inhibitors become utilized clinically. Here, we develop a yeast-based deep mutational scanning approach to systematically profile the activity of all possible single mutants of the SARS-CoV-2 3CL ^pro , and validate our results both in yeast and in authentic viruses. We reveal that the 3CL ^pro is highly malleable and is capable of tolerating mutations throughout the protein, including within the substrate binding pocket. Yet, we also identify specific residues that appear immutable for function of the protease, suggesting that these interactions may be novel targets for the design of future 3CL ^pro inhibitors. Finally, we utilize our screening results as a basis to identify E166V as a resistance-conferring mutation against the therapeutic 3CL ^pro inhibitor, nirmatrelvir, in clinical use. Collectively, the functional map presented herein may serve as a guide for further understanding of the biological properties of the 3CL protease and for drug development for current and future coronavirus pandemics.

Comprehensive fitness landscape of SARS-CoV-2 Mpro reveals insights into viral resistance mechanisms

With the continual evolution of new strains of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) that are more virulent, transmissible, and able to evade current vaccines, there is an urgent need for effective anti-viral drugs. The SARS-CoV-2 main protease (Mpro) is a leading target for drug design due to its conserved and indispensable role in the viral life cycle. Drugs targeting Mpro appear promising but will elicit selection pressure for resistance. To understand resistance potential in Mpro, we performed a comprehensive mutational scan of the protease that analyzed the function of all possible single amino acid changes. We developed three separate high throughput assays of Mpro function in yeast, based on either the ability of Mpro variants to cleave at a defined cut-site or on the toxicity of their expression to yeast. We used deep sequencing to quantify the functional effects of each variant in each screen. The protein fitness landscapes from all three screens were strongly correlated, indicating that they captured the biophysical properties critical to Mpro function. The fitness landscapes revealed a non-active site location on the surface that is extremely sensitive to mutation, making it a favorable location to target with inhibitors. In addition, we found a network of critical amino acids that physically bridge the two active sites of the Mpro dimer. The clinical variants of Mpro were predominantly functional in our screens, indicating that Mpro is under strong selection pressure in the human population. Our results provide predictions of mutations that will be readily accessible to Mpro evolution and that are likely to contribute to drug resistance. This complete mutational guide of Mpro can be used in the design of inhibitors with reduced potential of evolving viral resistance.

Can DyeCycling break the photobleaching limit in single-molecule FRET?

Biomolecular systems, such as proteins, crucially rely on dynamic processes at the nanoscale. Detecting biomolecular nanodynamics is therefore key to obtaining a mechanistic understanding of the energies and molecular driving forces that control biomolecular systems. Single-molecule fluorescence resonance energy transfer (smFRET) is a powerful technique to observe in real-time how a single biomolecule proceeds through its functional cycle involving a sequence of distinct structural states. Currently, this technique is fundamentally limited by irreversible photobleaching, causing the untimely end of the experiment and thus, a narrow temporal bandwidth of ≤ 3 orders of magnitude. Here, we introduce "DyeCycling", a measurement scheme with which we aim to break the photobleaching limit in smFRET. We introduce the concept of spontaneous dye replacement by simulations, and as an experimental proof-of-concept, we demonstrate the intermittent observation of a single biomolecule for one hour with a time resolution of milliseconds. Theoretically, DyeCycling can provide > 100-fold more information per single molecule than conventional smFRET. We discuss the experimental implementation of DyeCycling, its current and fundamental limitations, and specific biological use cases. Given its general simplicity and versatility, DyeCycling has the potential to revolutionize the field of time-resolved smFRET, where it may serve to unravel a wealth of biomolecular dynamics by bridging from milliseconds to the hour range.

Electronic Supplementary Material

Supplementary material is available for this article at 10.1007/s12274-022-4420-5 and is accessible for authorized users.

Key dimer interface residues impact the catalytic activity of 3CLpro, the main protease of SARS-CoV-2

3C-like protease (3CLpro) is one of two proteases that process and liberate functional viral proteins essential for the maturation and infectivity of severe acute respiratory syndrome-coronavirus 2 (SARS-CoV-2), the virus responsible for COVID-19. It has been suggested that 3CLpro is catalytically active as a dimer, making the dimerization interface a target for antiviral development. Guided by structural analysis, here we introduced single amino acid substitutions at nine residues at three key sites of the dimer interface to assess their impact on dimerization and activity. We show that at site 1, alanine substitution of S1 or E166 increased by 2-fold or reduced relative activity, respectively. At site 2, alanine substitution of S10 or E14 eliminated activity, whereas K12A exhibited ∼60% relative activity. At site 3, alanine substitution of R4, E290, or Q299 eliminated activity, whereas S139A exhibited 46% relative activity. We further found the oligomerization states of the dimer interface mutants varied; the inactive mutants R4A, R4Q, S10A/C, E14A/D/Q/S, E290A, and Q299A/E were present as dimers, demonstrating that dimerization is not an indication of catalytically active 3CLpro. In addition, present mostly as monomers, K12A displayed residual activity, which could be attributed to the conspicuous amount of dimer present. Finally, differential scanning calorimetry did not reveal a direct relationship between the thermodynamic stability of mutants with oligomerization or catalytic activity. These results provide insights on two allosteric sites, R4/E290 and S10/E14, that may promote the design of antiviral compounds that target the dimer interface rather than the active site of SARS-CoV-2 3CLpro.

Identification of bioactive molecules from Triphala (Ayurvedic herbal formulation) as potential inhibitors of SARS-CoV-2 main protease (Mpro) through computational investigations

Severe acute respiratory syndrome coronavirus disease (SARS-CoV-2) induced coronavirus disease 2019 (COVID-19) pandemic is the present worldwide health emergency. The global scientific community faces a significant challenge in developing targeted therapies to combat the SARS-CoV-2 infection. Computational approaches have been critical for identifying potential SARS-CoV-2 inhibitors in the face of limited resources and in this time of crisis. Main protease (M^pro) is an intriguing drug target because it processes the polyproteins required for SARS-CoV-2 replication. The application of Ayurvedic knowledge from traditional Indian systems of medicine may be a promising strategy to develop potential inhibitor for different target proteins of SARS-CoV-2. With this endeavor, we docked bioactive molecules from Triphala, an Ayurvedic formulation, against M^pro followed by molecular dynamics (MD) simulation (100 ns) to investigate their inhibitory potential against SARS-CoV-2. The top four best docked molecules (terflavin A, chebulagic acid, chebulinic acid, and corilagin) were selected for MD simulation study and the results obtained were compared to native ligand X77. From docking and MD simulation studies, the selected molecules showed promising binding affinity with the formation of stable complexes at the active binding pocket of M^pro and exhibited negative binding energy during MM-PBSA calculations, indication their strong binding affinity with the target protein. The identified bioactive molecules were further analyzed for drug-likeness by Lipinski's filter, ADMET and toxicity studies. Computational (in silico) investigations identified terflavin A, chebulagic acid, chebulinic acid, and corilagin from Triphala formulation as promising inhibitors of SARS-CoV-2 M^pro, suggesting experimental (in vitro/in vivo) studies to further explore their inhibitory mechanisms.

Targeting SARS-CoV-2 Proteases and Polymerase for COVID-19 Treatment: State of the Art and Future Opportunities

The newly emerged coronavirus, called SARS-CoV-2, is the causing pathogen of pandemic COVID-19. The identification of drugs to treat COVID-19 and other coronavirus diseases is an urgent global need, thus different strategies targeting either virus or host cell are still under investigation. Direct-acting agents, targeting protease and polymerase functionalities, represent a milestone in antiviral therapy. The 3C-like (or Main) protease (3CLpro) and the nsp12 RNA-dependent RNA-polymerase (RdRp) are the best characterized SARS-CoV-2 targets and show the highest degree of conservation across coronaviruses fostering the identification of broad-spectrum inhibitors. Coronaviruses also possess a papain-like protease, another essential enzyme, still poorly characterized and not equally conserved, limiting the identification of broad-spectrum agents. Herein, we provide an exhaustive comparative analysis of SARS-CoV-2 proteases and RdRp with respect to other coronavirus homologues. Moreover, we highlight the most promising inhibitors of these proteins reported so far, including the possible strategies for their further development.

Targeting SARS-CoV-2 Proteases and Polymerase for COVID-19 Treatment: State of the Art and Future Opportunities

The newly emerged coronavirus, called SARS-CoV-2, is the causing pathogen of pandemic COVID-19. The identification of drugs to treat COVID-19 and other coronavirus diseases is an urgent global need, thus different strategies targeting either virus or host cell are still under investigation. Direct-acting agents, targeting protease and polymerase functionalities, represent a milestone in antiviral therapy. The 3C-like (or Main) protease (3CL^pro) and the nsp12 RNA-dependent RNA-polymerase (RdRp) are the best characterized SARS-CoV-2 targets and show the highest degree of conservation across coronaviruses fostering the identification of broad-spectrum inhibitors. Coronaviruses also possess a papain-like protease, another essential enzyme, still poorly characterized and not equally conserved, limiting the identification of broad-spectrum agents. Herein, we provide an exhaustive comparative analysis of SARS-CoV-2 proteases and RdRp with respect to other coronavirus homologues. Moreover, we highlight the most promising inhibitors of these proteins reported so far, including the possible strategies for their further development.

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

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

Evaluation of RevX Solution Extract as a Potential Inhibitor of the Main Protease of SARS-CoV-2—In vitro Study and Molecular Docking

The main protease (M^pro) of SARS-CoV-2 is a protease necessary for viral polyprotein processing and maturation. M^pro cleaves the polypeptide sequence after the glutamine residues. There is no known cellular protease with this substrate specificity in humans; therefore, it is considered an attractive drug target. Previously, fermented sorghum extract RevX (trademark of Revolutrx INC.) solution significantly alleviated physical decline and complications in a patient with lung adenocarcinoma, suggesting the role of bioactive components in RevX solution. To further explore whether the bioactive components in RevX solution exhibit other biological activities, such as antiviral effects, we investigated its inhibitory effect on the M^pro of SARS-CoV-2 virus. We report herein that the solid extract of the RevX solution exhibits an efficacious M^pro inhibitory activity, with IC₅₀ of 2.07 ± 0.38 μg/mL. Molecular docking of sterol-like components in the RevX extracts identified by MS shows that the three sterol-like molecules can bind to the active region of the GC376-M^pro complex, supporting the structure-function relationship. Combined with its ability to significantly alleviate the body's immunity decline and to inhibit the activity of SARS-CoV-2 M^pro, RevX solution may provide a possible alternative supportive treatment for patients with COVID-19.

•

The unique substrate specificity of SARS-CoV-2 M^pro makes it a potential target for drug design.

•

Fermented sorghum extract RevX solution enhances adjuvant therapy of lung adenocarcinoma, suggesting the role of bioactive components in RevX solution.

•

The solid extract of RevX showed potent M^pro inhibitory activity with IC₅₀ of 2.07 ± 0.38 μg/mL.

•

The three sterol-like structures of RevX extract showed a similar binding cavity to Mpro-GC376, suggesting its putative inhibitory activity.

Applications of density functional theory in COVID-19 drug modeling

The rapidly evolving Coronavirus 2019 (COVID-19) pandemic has led to millions of deaths around the world, highlighting the pressing need to develop effective antiviral pharmaceuticals. Recent efforts with computer-aided rational drug discovery have allowed detailed examination of drug–macromolecule interactions primarily by molecular mechanics (MM) techniques. Less widely applied in COVID-19 drug modeling is density functional theory (DFT), a quantum mechanics (QM) method that enables electronic structure calculations and elucidations of reaction mechanisms. Here, we review recent advances in applying DFT in molecular modeling studies of COVID-19 pharmaceuticals. We start by providing an overview of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) drugs and targets, followed by a brief introduction to DFT. We then provide a discussion of different approaches by which DFT has been applied. Finally, we discuss essential factors to consider when incorporating DFT in future drug modeling research.

The dolabellane diterpenes as potential inhibitors of the SARS-CoV-2 main protease: molecular insight of the inhibitory mechanism through computational studies

An investigation has been carried out on natural products from dolabellane derivatives to understand their potential in inhibiting the SARS-CoV-2 main protease (3CLpro) using an in silico approach. Inhibition of the 3CLpro enzyme is a promising target in stopping the replication of the SARS-CoV-2 virus through inhibition of the subsite binding pocket. The redocking process aims to determine the 3CLpro active sites. The redocking requirement showed a good pose with an RMSD value of 1.39 Å. The combination of molecular docking and MD simulation shows the results of DD13 as a candidate which had a good binding affinity (kcal mol-1) to inhibit the 3CLpro enzyme activity. Prediction of binding free energy (kcal mol-1) of DD13 using the Molecular Mechanics-Poisson Boltzmann/Generalized Born Surface Area (MM-PB/GBSA) approach shows the results ΔG bind(MM-GBSA): -52.33 ± 0.34 and ΔG bind(MM-PBSA): -43.52 ± 0.42. The key residues responsible for the inhibition mechanism are Hie41, Ser46, Met49, Asn142, Cys145, Hie163, Met165, and Gln189. Additionally, pharmacokinetic prediction recommended that DD13 had promising criteria as a drug candidate. The results demonstrated in this study provide theoretical information to obtain a potential inhibitor against the SARS-CoV-2 main protease.

Allosteric Regulation of 3CL Protease of SARS-CoV-2 and SARS-CoV Observed in the Crystal Structure Ensemble

The 3C-like protease (3CL^pro) of SARS-CoV-2 is a potential therapeutic target for COVID-19. Importantly, it has an abundance of structural information solved as a complex with various drug candidate compounds. Collecting these crystal structures (83 Protein Data Bank (PDB) entries) together with those of the highly homologous 3CL^pro of SARS-CoV (101 PDB entries), we constructed the crystal structure ensemble of 3CL^pro to analyze the dynamic regulation of its catalytic function. The structural dynamics of the 3CL^pro dimer observed in the ensemble were characterized by the motions of four separate loops (the C-loop, E-loop, H-loop, and Linker) and the C-terminal domain III on the rigid core of the chymotrypsin fold. Among the four moving loops, the C-loop (also known as the oxyanion binding loop) causes the order (active)-disorder (collapsed) transition, which is regulated cooperatively by five hydrogen bonds made with the surrounding residues. The C-loop, E-loop, and Linker constitute the major ligand binding sites, which consist of a limited variety of binding residues including the substrate binding subsites. Ligand binding causes a ligand size dependent conformational change to the E-loop and Linker, which further stabilize the C-loop via the hydrogen bond between the C-loop and E-loop. The T285A mutation from SARS-CoV 3CL^pro to SARS-CoV-2 3CL^pro significantly closes the interface of the domain III dimer and allosterically stabilizes the active conformation of the C-loop via hydrogen bonds with Ser1 and Gly2; thus, SARS-CoV-2 3CL^pro seems to have increased activity relative to that of SARS-CoV 3CL^pro.

Unveiling the Effect of Low pH on the SARS-CoV-2 Main Protease by Molecular Dynamics Simulations

(1) Background: Main Protease (Mpro) is an attractive therapeutic target that acts in the replication and transcription of the SARS-CoV-2 coronavirus. Mpro is rich in residues exposed to protonation/deprotonation changes which could affect its enzymatic function. This work aimed to explore the effect of the protonation/deprotonation states of Mpro at different pHs using computational techniques. (2) Methods: The different distribution charges were obtained in all the evaluated pHs by the Semi-Grand Canonical Monte Carlo (SGCMC) method. A set of Molecular Dynamics (MD) simulations was performed to consider the different protonation/deprotonation during 250 ns, verifying the structural stability of Mpro at different pHs. (3) Results: The present findings demonstrate that active site residues and residues that allow Mpro dimerisation was not affected by pH changes. However, Mpro substrate-binding residues were altered at low pHs, allowing the increased pocket volume. Additionally, the results of the solvent distribution around Sγ, Hγ, Nδ1 and Hδ1 atoms of the catalytic residues Cys145 and His41 showed a low and high-water affinity at acidic pH, respectively. It which could be crucial in the catalytic mechanism of SARS-CoV-2 Mpro at low pHs. Moreover, we analysed the docking interactions of PF-00835231 from Pfizer in the preclinical phase, which shows excellent affinity with the Mpro at different pHs. (4) Conclusion: Overall, these findings indicate that SARS-CoV-2 Mpro is highly stable at acidic pH conditions, and this inhibitor could have a desirable function at this condition.

Natural Products-Based Drug Design against SARS-CoV-2 Mpro 3CLpro

Coronavirus disease 2019 (COVID-19), caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has received global attention due to the serious threat it poses to public health. Since the outbreak in December 2019, millions of people have been affected and its rapid global spread has led to an upsurge in the search for treatment. To discover hit compounds that can be used alone or in combination with repositioned drugs, we first analyzed the pharmacokinetic and toxicological properties of natural products from Brazil's semiarid region. After, we analyzed the site prediction and druggability of the SARS-CoV-2 main protease (Mpro), followed by docking and molecular dynamics simulation. The best SARS-CoV-2 Mpro complexes revealed that other sites were accessed, confirming that our approach could be employed as a suitable starting protocol for ligand prioritization, reinforcing the importance of catalytic cysteine-histidine residues and providing new structural data that could increase the antiviral development mainly against SARS-CoV-2. Here, we selected 10 molecules that could be in vitro assayed in response to COVID-19. Two compounds (b01 and b02) suggest a better potential for interaction with SARS-CoV-2 Mpro and could be further studied.

Discovery of SARS-CoV-2 Mpro peptide inhibitors from modelling substrate and ligand binding

The main protease (Mpro) of SARS-CoV-2 is central to viral maturation and is a promising drug target, but little is known about structural aspects of how it binds to its 11 natural cleavage sites. We used biophysical and crystallographic data and an array of biomolecular simulation techniques, including automated docking, molecular dynamics (MD) and interactive MD in virtual reality, QM/MM, and linear-scaling DFT, to investigate the molecular features underlying recognition of the natural Mpro substrates. We extensively analysed the subsite interactions of modelled 11-residue cleavage site peptides, crystallographic ligands, and docked COVID Moonshot-designed covalent inhibitors. Our modelling studies reveal remarkable consistency in the hydrogen bonding patterns of the natural Mpro substrates, particularly on the N-terminal side of the scissile bond. They highlight the critical role of interactions beyond the immediate active site in recognition and catalysis, in particular plasticity at the S2 site. Building on our initial Mpro-substrate models, we used predictive saturation variation scanning (PreSaVS) to design peptides with improved affinity. Non-denaturing mass spectrometry and other biophysical analyses confirm these new and effective 'peptibitors' inhibit Mpro competitively. Our combined results provide new insights and highlight opportunities for the development of Mpro inhibitors as anti-COVID-19 drugs.

Small molecule interactions with the SARS-CoV-2 main protease: In silico all-atom microsecond MD simulations, PELE Monte Carlo simulations, and determination of in vitro activity inhibition

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has caused the ongoing COVID-19 pandemic. With some notable exceptions, safe and effective vaccines, which are now being widely distributed globally, have largely begun to stabilise the situation. However, emerging variants of concern and vaccine hesitancy are apparent obstacles to eradication. Therefore, the need for the development of potent antivirals is still of importance. In this context, the SARS-CoV-2 main protease (M^pro) is a critical target and numerous clinical trials, predominantly in the private domain, are currently in progress. Here, our aim was to extend our previous studies, with hypericin and cyanidin-3-O-glucoside, as potential inhibitors of the SARS-CoV-2 M^pro. Firstly, we performed all-atom microsecond molecular dynamics simulations, which highlight the stability of the ligands in the M^pro active site over the duration of the trajectories. We also invoked PELE Monte Carlo simulations which indicate that both hypericin and cyanidin-3-O-glucoside preferentially interact with the M^pro active site and known allosteric sites. For further validation, we performed an in vitro enzymatic activity assay that demonstrated that hypericin and cyanidin-3-O-glucoside inhibit M^pro activity in a dose-dependent manner at biologically relevant (μM) concentrations. However, both ligands are much less potent than the well-known covalent antiviral GC376, which was used as a positive control in our experiments. Nevertheless, the biologically relevant activity of hypericin and cyanidin-3-O-glucoside is encouraging. In particular, a synthetic version of hypericin has FDA orphan drug designation, which could simplify potential clinical evaluation in the context of COVID-19.

Mechanistic Insights into the Inhibition of SARS-CoV-2 Main Protease by Clovamide and Its Derivatives: In Silico Studies

The novel coronavirus SARS-CoV-2 Main Protease (Mpro) is an internally encoded enzyme that hydrolyzes the translated polyproteins at designated sites. The protease directly mediates viral replication processes; hence, a promising target for drug design. Plant-based natural products, especially polyphenols and phenolic compounds, provide the scaffold for many effective antiviral medications, and have recently been shown to be able to inhibit Mpro of SARS-CoV-2. Specifically, polyphenolic compounds found in cacao and chocolate products have been shown by recent experimental studies to have strong inhibitory effects against Mpro activities. This work aims to uncover the inhibition processes of Mpro by a natural phenolic compound found in cacao and chocolate products, clovamide. Clovamide (caffeoyl-DOPA) is a naturally occurring caffeoyl conjugate that is found in the phenolic fraction of Theobroma Cacao L. and a potent radical-scavenging antioxidant as suggested by previous studies of our group. Here, we propose inhibitory mechanisms by which clovamide may act as a Mpro inhibitor as it becomes oxidized by scavenging reactive oxygen species (ROS) in the body, or becomes oxidized as a result of enzymatic browning. We use molecular docking, annealing-based molecular dynamics, and Density Functional Theory (DFT) calculations to study the interactions between clovamide with its derivatives and Mpro catalytic and allosteric sites. Our molecular modelling studies provide mechanistic insights of clovamide inhibition of Mpro, and indicate that clovamide may be a promising candidate as a drug lead molecule for COVID-19 treatments.

Pharmacoinformatics and molecular dynamics simulation studies reveal potential covalent and FDA-approved inhibitors of SARS-CoV-2 main protease 3CLpro

The SARS-CoV-2 was confirmed to cause the global pandemic of coronavirus disease 2019 (COVID-19). The 3-chymotrypsin-like protease (3CLpro), an essential enzyme for viral replication, is a valid target to combat SARS-CoV and MERS-CoV. In this work, we present a structure-based study to identify potential covalent inhibitors containing a variety of chemical warheads. The targeted Asinex Focused Covalent (AFCL) library was screened based on different reaction types and potential covalent inhibitors were identified. In addition, we screened FDA-approved protease inhibitors to find candidates to be repurposed against SARS-CoV-2 3CLpro. A number of compounds with significant covalent docking scores were identified. These compounds were able to establish a covalent bond (C-S) with the reactive thiol group of Cys145 and to form favorable interactions with residues lining the substrate-binding site. Moreover, paritaprevir and simeprevir from FDA-approved protease inhibitors were identified as potential inhibitors of SARS-CoV-2 3CLpro. The mechanism and dynamic stability of binding between the identified compounds and SARS-CoV-2 3CLpro were characterized by molecular dynamics (MD) simulations. The identified compounds are potential inhibitors worthy of further development as COVID-19 drugs. Importantly, the identified FDA-approved anti-hepatitis-C virus (HCV) drugs paritaprevir and simeprevir could be ready for clinical trials to treat infected patients and help curb COVID-19. Communicated by Ramaswamy H. Sarma.

Comprehensive virtual screening of 4.8 k flavonoids reveals novel insights into allosteric inhibition of SARS-CoV-2 MPRO

SARS-CoV-2 main protease is a common target for inhibition assays due to its high conservation among coronaviruses. Since flavonoids show antiviral activity, several in silico works have proposed them as potential SARS-CoV-2 main protease inhibitors. Nonetheless, there is reason to doubt certain results given the lack of consideration for flavonoid promiscuity or main protease plasticity, usage of short library sizes, absence of control molecules and/or the limitation of the methodology to a single target site. Here, we report a virtual screening study where dorsilurin E, euchrenone a11, sanggenol O and CHEMBL2171598 are proposed to inhibit main protease through different pathways. Remarkably, novel structural mechanisms were observed after sanggenol O and CHEMBL2171598 bound to experimentally proven allosteric sites. The former drastically affected the active site, while the latter triggered a hinge movement which has been previously reported for an inactive SARS-CoV main protease mutant. The use of a curated database of 4.8 k flavonoids, combining two well-known docking software (AutoDock Vina and AutoDock4.2), molecular dynamics and MMPBSA, guaranteed an adequate analysis and robust interpretation. These criteria can be considered for future screening campaigns against SARS-CoV-2 main protease.

S494 O-glycosylation site on the SARS-CoV-2 RBD affects the virus affinity to ACE2 and its infectivity; a molecular dynamics study

SARS-CoV-2 is a strain of Coronavirus family that caused the ongoing pandemic of COVID-19. Several studies showed that the glycosylation of virus spike (S) protein and the Angiotensin-Converting Enzyme 2 (ACE2) receptor on the host cell is critical for the virus infectivity. Molecular Dynamics (MD) simulations were used to explore the role of a novel mutated O-glycosylation site (D494S) on the Receptor Binding Domain (RBD) of S protein. This site was suggested as a key mediator of virus-host interaction. By exploring the dynamics of three O-glycosylated models and the control systems of unglcosylated S4944 and S494D complexes, it was shown that the decoration of S494 with elongated O-glycans results in stabilized interactions on the direct RBD-ACE2. Calculation of the distances between RBD and two major H1, H2 helices of ACE2 and the interacting pairs of amino acids in the interface showed that the elongated O-glycan maintains these interactions by forming several polar contacts with the neighbouring residues while it would not interfere in the direct binding interface. Relative binding free energy of RBD-ACE2 is also more favorable in the O-glycosylated models with longer glycans. The increase of RBD binding affinity to ACE2 depends on the size of attached O-glycan. By increasing the size of O-glycan, the RBD-ACE2 binding affinity will increase. Hence, this crucial factor must be taken into account for any further inhibitory approaches towards RBD-ACE2 interaction.

Interfacial Water Many-Body Effects Drive Structural Dynamics and Allosteric Interactions in SARS-CoV-2 Main Protease Dimerization Interface

Following our previous work ( Chem. Sci. 2021, 12, 4889-4907), we study the structural dynamics of the SARS-CoV-2 Main Protease dimerization interface (apo dimer) by means of microsecond adaptive sampling molecular dynamics simulations (50 μs) using the AMOEBA polarizable force field (PFF). This interface is structured by a complex H-bond network that is stable only at physiological pH. Structural correlations analysis between its residues and the catalytic site confirms the presence of a buried allosteric site. However, noticeable differences in allosteric connectivity are observed between PFFs and non-PFFs. Interfacial polarizable water molecules are shown to appear at the heart of this discrepancy because they are connected to the global interface H-bond network and able to adapt their dipole moment (and dynamics) to their diverse local physicochemical microenvironments. The water-interface many-body interactions appear to drive the interface volume fluctuations and to therefore mediate the allosteric interactions with the catalytic cavity.