Drug repurposing: identification of SARS-CoV-2 potential inhibitors by virtual screening and pharmacokinetics strategies.

INTRODUCTION: The coronavirus disease 2019 (COVID-19) pandemic caused global health, economic, and population loss. Variants of the coronavirus contributed to the severity of the disease and persistent rise in infections. This study aimed to identify potential drug candidates from fifteen approved antiviral drugs against SARS-CoV-2 (6LU7), SARS-CoV (5B6O), and SARS-CoV-2 spike protein (6M0J) using virtual screening and pharmacokinetics to gain insights into COVID-19 therapeutics. METHODOLOGY: We employed drug repurposing approach to analyze binding performance of fifteen clinically approved antiviral drugs against the main protease of SARS-CoV-2 (6LU7), SARS-CoV (5B6O), and SARS-CoV-2 spike proteins bound to ACE-2 receptor (6M0J), to provide an insight into the therapeutics of COVID-19. AutoDock Vina was used for docking studies. The binding affinities were calculated, and 2-3D structures of protein-ligand interactions were drawn. RESULTS: Rutin, hesperidin, and nelfinavir are clinically approved antiviral drugs with high binding affinity to proteins 6LU7, 5B6O, and 6M0J. These ligands have excellent pharmacokinetics, ensuring efficient absorption, metabolism, excretion, and digestibility. Hesperidin showed the most potent interaction with spike protein 6M0J, forming four H-bonds. Nelfinavir had a high human intestinal absorption (HIA) score of 0.93, indicating maximum absorption in the body and promising interactions with 6LU7. CONCLUSIONS: Our results indicated that rutin, hesperidin, and nelfinavir had the highest binding results against the proposed drug targets. The computational approach effectively identified SARS-CoV-2 inhibitors. COVID-19 is still a recurrent threat globally and predictive analysis using natural compounds might serve as a starting point for new drug development against SARS-CoV-2 and related viruses.

Computational Estimation of Residues Involving Resistance to the SARS-CoV-2 Main Protease Inhibitor Ensitrelvir Based on Virtual Alanine Scan of the Active Site.

Ensitrelvir is a noncovalent inhibitor of the main protease (Mpro) of severe acute respiratory syndrome coronavirus 2. Acquisition of drug resistance in virus-derived proteins is a serious therapeutic concern, and drug resistance occurs due to amino acid mutations. In this study, we computationally constructed 24 mutants, in which one residue around the active site was replaced with alanine and performed molecular dynamics simulations to the complex of Mpro and ensitrelvir to predict the residues involved in drug resistance. We evaluated the changes in the entire protein structure and ligand configuration in each of these mutants and estimated which residues were involved in ensitrelvir recognition. This method is called a virtual alanine scan. In nine mutants (S1A, T26A, H41A, M49A, L141A, H163A, E166A, V186A, and R188A), although the entire protein structure and catalytic dyad (cysteine (Cys)145 and histidine (His)41) were not significantly moved, the ensitrelvir configuration changed. Thus, it is considered that these mutants did not recognize ensitrelvir while maintaining Mpro enzymatic activities, and Ser1, Thr26, His41, Met49, Leu141, His163, Glu166, Val186, and Arg188 may be related to ensitrelvir resistance. The ligand shift noted in M49A was similar to that observed in M49I, which has been shown to be experimentally ensitrelvir resistant. These findings suggest that our research approach can predict mutations that incite drug resistance.

Isolation, Virtual Screening, and Evaluation of Hazelnut-Derived Immunoactive Peptides for the Inhibition of SARS-CoV-2 Main Protease.

The aim of this study is to validate the activity of hazelnut (Corylus avellana L.)-derived immunoactive peptides inhibiting the main protease (Mpro) of SARS-CoV-2 and further unveil their interaction mechanism using in vitro assays, molecular dynamics (MD) simulations, and binding free energy calculations. In general, the enzymatic hydrolysis components, especially molecular weight < 3 kDa, possess good immune activity as measured by the proliferation ability of mouse splenic lymphocytes and phagocytic activity of mouse peritoneal macrophages. Over 866 unique peptide sequences were isolated, purified, and then identified by nanohigh-performance liquid chromatography/tandem mass spectrometry (NANO-HPLC-MS/MS) from hazelnut protein hydrolysates, but Trp-Trp-Asn-Leu-Asn (WWNLN) and Trp-Ala-Val-Leu-Lys (WAVLK) in particular are found to increase the cell viability and phagocytic capacity of RAW264.7 macrophages as well as promote the secretion of the cytokines nitric oxide (NO), tumor necrosis factor-α (TNF-α), and interleukin-1ß (IL-1ß). Fluorescence resonance energy transfer assay elucidated that WWNLN and WAVLK exhibit excellent inhibitory potency against Mpro, with IC50 values of 6.695 and 16.750 µM, respectively. Classical all-atom MD simulations show that hydrogen bonds play a pivotal role in stabilizing the complex conformation and protein-peptide interaction. Molecular Mechanics/Generalized Born Surface Area (MM/GBSA) calculation indicates that WWNLN has a lower binding free energy with Mpro than WAVLK. Furthermore, adsorption, distribution, metabolism, excretion, and toxicity (ADMET) predictions illustrate favorable drug-likeness and pharmacokinetic properties of WWNLN compared to WAVLK. This study provides a new understanding of the immunomodulatory activity of hazelnut hydrolysates and sheds light on peptide inhibitors targeting Mpro.

Phytoconstituents of Artemisia Annua as potential inhibitors of SARS CoV2 main protease: an in silico study.

BACKGROUND: In November 2019, the world faced a pandemic called SARS-CoV-2, which became a major threat to humans and continues to be. To overcome this, many plants were explored to find a cure. METHODS: Therefore, this research was planned to screen out the active constituents from Artemisia annua that can work against the viral main protease Mpro as this non-structural protein is responsible for the cleavage of replicating enzymes of the virus. Twenty-five biocompounds belonging to different classes namely alpha-pinene, beta-pinene, carvone, myrtenol, quinic acid, caffeic acid, quercetin, rutin, apigenin, chrysoplenetin, arteannunin b, artemisinin, scopoletin, scoparone, artemisinic acid, deoxyartemisnin, artemetin, casticin, sitogluside, beta-sitosterol, dihydroartemisinin, scopolin, artemether, artemotil, artesunate were selected. Virtual screening of these ligands was carried out against drug target Mpro by CB dock. RESULTS: Quercetin, rutin, casticin, chrysoplenetin, apigenin, artemetin, artesunate, sopolin and sito-gluside were found as hit compounds. Further, ADMET screening was conducted which represented Chrysoplenetin as a lead compound. Azithromycin was used as a standard drug. The interactions were studied by PyMol and visualized in LigPlot. Furthermore, the RMSD graph shows fluctuations at various points at the start of simulation in Top1 (Azithromycin) complex system due to structural changes in the helix-coil-helix and beta-turn-beta changes at specific points resulting in increased RMSD with a time frame of 50 ns. But this change remains stable after the extension of simulation time intervals till 100 ns. On other side, the Top2 complex system remains highly stable throughout the time scale. No such structural dynamics were observed bu the ligand attached to the active site residues binds strongly. CONCLUSION: This study facilitates researchers to develop and discover more effective and specific therapeutic agents against SARS-CoV-2 and other viral infections. Finally, chrysoplenetin was identified as a more potent drug candidate to act against the viral main protease, which in the future can be helpful.

Synthesis and in silico inhibitory action studies of azo-anchored imidazo[4,5-b]indole scaffolds against the COVID-19 main protease (Mpro).

The present work elicits a novel approach to combating COVID-19 by synthesizing a series of azo-anchored 3,4-dihydroimidazo[4,5-b]indole derivatives. The envisaged methodology involves the L-proline-catalyzed condensation of para-amino-functionalized azo benzene, indoline-2,3-dione, and ammonium acetate precursors with pertinent aryl aldehyde derivatives under ultrasonic conditions. The structures of synthesized compounds were corroborated through FT-IR, 1H NMR, 13C NMR, and mass analysis data. Molecular docking studies assessed the inhibitory potential of these compounds against the main protease (Mpro) of SARS-CoV-2. Remarkably, in silico investigations revealed significant inhibitory action surpassing standard drugs such as Remdesivir, Paxlovid, Molnupiravir, Chloroquine, Hydroxychloroquine (HCQ), and (N3), an irreversible Michael acceptor inhibitor. Furthermore, the highly active compound was also screened for cytotoxicity activity against HEK-293 cells and exhibited minimal toxicity across a range of concentrations, affirming its favorable safety profile and potential suitability. The pharmacokinetic properties (ADME) of the synthesized compounds have also been deliberated. This study paves the way for in vitro and in vivo testing of these scaffolds in the ongoing battle against SARS-CoV-2.

Accelerating reliable multiscale quantum refinement of protein-drug systems enabled by machine learning.

Biomacromolecule structures are essential for drug development and biocatalysis. Quantum refinement (QR) methods, which employ reliable quantum mechanics (QM) methods in crystallographic refinement, showed promise in improving the structural quality or even correcting the structure of biomacromolecules. However, vast computational costs and complex quantum mechanics/molecular mechanics (QM/MM) setups limit QR applications. Here we incorporate robust machine learning potentials (MLPs) in multiscale ONIOM(QM:MM) schemes to describe the core parts (e.g., drugs/inhibitors), replacing the expensive QM method. Additionally, two levels of MLPs are combined for the first time to overcome MLP limitations. Our unique MLPs+ONIOM-based QR methods achieve QM-level accuracy with significantly higher efficiency. Furthermore, our refinements provide computational evidence for the existence of bonded and nonbonded forms of the Food and Drug Administration (FDA)-approved drug nirmatrelvir in one SARS-CoV-2 main protease structure. This study highlights that powerful MLPs accelerate QRs for reliable protein-drug complexes, promote broader QR applications and provide more atomistic insights into drug development.

SARS-CoV-2 Mpro responds to oxidation by forming disulfide and NOS/SONOS bonds.

The main protease (Mpro) of SARS-CoV-2 is critical for viral function and a key drug target. Mpro is only active when reduced; turnover ceases upon oxidation but is restored by re-reduction. This suggests the system has evolved to survive periods in an oxidative environment, but the mechanism of this protection has not been confirmed. Here, we report a crystal structure of oxidized Mpro showing a disulfide bond between the active site cysteine, C145, and a distal cysteine, C117. Previous work proposed this disulfide provides the mechanism of protection from irreversible oxidation. Mpro forms an obligate homodimer, and the C117-C145 structure shows disruption of interactions bridging the dimer interface, implying a correlation between oxidation and dimerization. We confirm dimer stability is weakened in solution upon oxidation. Finally, we observe the protein's crystallization behavior is linked to its redox state. Oxidized Mpro spontaneously forms a distinct, more loosely packed lattice. Seeding with crystals of this lattice yields a structure with an oxidation pattern incorporating one cysteine-lysine-cysteine (SONOS) and two lysine-cysteine (NOS) bridges. These structures further our understanding of the oxidative regulation of Mpro and the crystallization conditions necessary to study this structurally.

SARS-CoV-2 Mpro oligomerization as a potential target for therapy.

The main protease (Mpro) of SARS-CoV-2 is critical in the virus's replication cycle, facilitating the maturation of polyproteins into functional units. Due to its conservation across taxa, Mpro is a promising target for broad-spectrum antiviral drugs. Targeting Mpro with small molecule inhibitors, such as nirmatrelvir combined with ritonavir (Paxlovid™), which the FDA has approved for post-exposure treatment and prophylaxis, can effectively interrupt the replication process of the virus. A key aspect of Mpro's function is its ability to form a functional dimer. However, the mechanics of dimerization and its influence on proteolytic activity remain less understood. In this study, we utilized biochemical, structural, and molecular modelling approaches to explore Mpro dimerization. We evaluated critical residues, specifically Arg4 and Arg298, that are essential for dimerization. Our results show that changes in the oligomerization state of Mpro directly affect its enzymatic activity and dimerization propensity. We discovered a synergistic relationship influencing dimer formation, involving both intra- and intermolecular interactions. These findings highlight the potential for developing allosteric inhibitors targeting Mpro, offering promising new directions for therapeutic strategies.

New meroterpenoids and polyketides from the endophytic fungus Paraphaeosphaeria sp. C-XB-J-1 and their anti-inflammatory and SARS-CoV-2 Mpro inhibitory activities.

Seven new meroterpenoids, paraphaeones A-G (1-7), and two new polyketides, paraphaeones H-I (8-9), along with eight known compounds (10-17), were isolated from the endophytic fungus Paraphaeosphaeria sp. C-XB-J-1. The structures of 1-9 were identified through the analysis of 1H, 13C, and 2D NMR spectra, assisted by HR-ESI-MS data. Compounds 1 and 7 exhibited a dose-dependent decrease in lactate dehydrogenase levels, with IC50 values of 1.78 µM and 1.54 µM, respectively. Moreover, they inhibited the secretion of IL-1ß and CASP-1, resulting in a reduction in the activity levels of NLRP3 inflammasomes. Fluorescence microscopy results indicated that compound 7 concentration-dependently attenuated cell pyroptosis. Additionally, compounds 4 and 7 showed potential inhibitory effects on the severe acute respiratory syndrome coronavirus-2 main protease (SARS-CoV-2 Mpro), with IC50 values of 10.8 ± 0.9 µM and 12.9 ± 0.7 µM, respectively.

In vitro biological evaluation and in silico insights into the antiviral activity of standardized olive leaves extract against SARS-CoV-2.

There is still a great global need for efficient treatments for the management of SARS-CoV-2 illness notwithstanding the availability and efficacy of COVID-19 vaccinations. Olive leaf is an herbal remedy with a potential antiviral activity that could improve the recovery of COVID-19 patients. In this work, the olive leaves major metabolites were screened in silico for their activity against SARS-CoV-2 by molecular docking on several viral targets such as methyl transferase, helicase, Plpro, Mpro, and RdRp. The results of in silico docking study showed that olive leaves phytoconstituents exhibited strong potential antiviral activity against SARS-CoV-2 selected targets. Verbacoside demonstrated a strong inhibition against methyl transferase, helicase, Plpro, Mpro, and RdRp (docking scores = -17.2, -20, -18.2, -19.8, and -21.7 kcal/mol.) respectively. Oleuropein inhibited 5rmm, Mpro, and RdRp (docking scores = -15, -16.6 and -18.6 kcal/mol., respectively) respectively. Apigenin-7-O-glucoside exhibited activity against methyl transferase and RdRp (docking score = -16.1 and -19.4 kcal/mol., respectively) while Luteolin-7-O-glucoside inhibited Plpro and RdRp (docking score = -15.2 and -20 kcal/mol., respectively). The in vitro antiviral assay was carried out on standardized olive leaf extract (SOLE) containing 20% oleuropein and IC50 was calculated. The results revealed that 20% SOLE demonstrated a moderate antiviral activity against SARS-CoV-2 with IC50 of 118.3 µg /mL. Accordingly, olive leaf could be a potential herbal therapy against SARS-CoV-2 but more in vivo and clinical investigations are recommended.

De novo design of potential peptide analogs against the main protease of Omicron variant using in silico studies.

SARS-CoV-2 and its variants are crossing the immunity barrier induced through vaccination. Recent Omicron sub-variants are highly transmissible and have a low mortality rate. Despite the low severity of Omicron variants, these new variants are known to cause acute post-infectious syndromes. Nowadays, novel strategies to develop new potential inhibitors for SARS-CoV-2 and other Omicron variants have gained prominence. For viral replication and survival the main protease of SARS-CoV-2 plays a vital role. Peptide-like inhibitors that mimic the substrate peptide have already proved to be effective in inhibiting the Mpro of SARS-CoV-2 variants. Our systematic canonical amino acid point mutation analysis on the native peptide has revealed various ways to improve the native peptide of the main protease. Multi mutation analysis has led us to identify and design potent peptide-analog inhibitors that act against the Mpro of the Omicron sub-variants. Our in-depth analysis of all-atom molecular dynamics studies has paved the way to characterize the atomistic behavior of Mpro in Omicron variants. Our goal is to develop potent peptide-analogs that could be therapeutically effective against Omicron and its sub-variants.

Anti-Coronavirus Potential of Polyether Ionophores: The New Application of Veterinary Antibiotics in Livestock.

Coronaviruses have consistently posed a major global concern in the field of livestock industry and public health. However, there is currently a lack of efficient drugs with broad-spectrum antiviral activity to address the challenges presented by emerging mutated strains or drug resistance. Additionally, the method for identifying multitarget drugs is also insufficient. Aminopeptidase N (APN) and 3C-like proteinase (3CLpro) represent promising targets for host-directed and virus-directed strategies, respectively, in the development of effective drugs against various coronaviruses. In this study, maduramycin ammonium demonstrated a broad-spectrum antiviral effect by targeting both of the proteins. The binding domains 4 Å from the ligand of both target proteins shared a structural similarity, suggesting that screening and designing drugs based on these domains might exhibit broad-spectrum and highly effective antiviral activity. Furthermore, it was identified that the polyether ionophores' ability to carry zinc ion might be one of the reasons why they were able to target APN and exhibit antiviral effect. The findings of this experiment provide novel perspectives for future drug screening and design, while also offering valuable references for the utilization of polyether ionophores in the management of livestock health.

Exploring the Efficacy of Noncovalent SARS-CoV-2 Main Protease Inhibitors: A Computational Simulation Analysis Study.

The SARS-CoV-2 main protease, as a key target for antiviral therapeutics, is instrumental in maintaining virus stability, facilitating translation, and enabling the virus to evade innate immunity. Our research focused on designing non-covalent inhibitors to counteract the action of this protease. Utilizing a 3D-QSAR model and contour map, we successfully engineered eight novel non-covalent inhibitors. Further evaluation and comparison of these novel compounds through methodologies including molecular docking, ADMET analysis, frontier molecular orbital studies, molecular dynamics simulations, and binding free energy revealed that the inhibitors N02 and N03 demonstrated superior research performance (N02 ΔGbind=-206.648 kJ/mol, N03 ΔGbind=-185.602 kJ/mol). These findings offer insightful guidance for the further refinement of molecular structures and the development of more efficacious inhibitors. Consequently, future investigations can draw upon these findings to unearth more potent inhibitors, thereby amplifying their impact in the treatment and prevention of associated diseases.

Discovery of Covalent Lead Compounds Targeting 3CL Protease with a Lateral Interactions Spiking Neural Network.

Covalent drugs exhibit advantages in that noncovalent drugs cannot match, and covalent docking is an important method for screening covalent lead compounds. However, it is difficult for covalent docking to screen covalent compounds on a large scale because covalent docking requires determination of the covalent reaction type of the compound. Here, we propose to use deep learning of a lateral interactions spiking neural network to construct a covalent lead compound screening model to quickly screen covalent lead compounds. We used the 3CL protease (3CL Pro) of SARS-CoV-2 as the screen target and constructed two classification models based on LISNN to predict the covalent binding and inhibitory activity of compounds. The two classification models were trained on the covalent complex data set targeting cysteine (Cys) and the compound inhibitory activity data set targeting 3CL Pro, respected, with good prediction accuracy (ACC > 0.9). We then screened the screening compound library with 6 covalent binding screening models and 12 inhibitory activity screening models. We tested the inhibitory activity of the 32 compounds, and the best compound inhibited SARS-CoV-2 3CL Pro with an IC50 value of 369.5 nM. Further assay implied that dithiothreitol can affect the inhibitory activity of the compound to 3CL Pro, indicating that the compound may covalently bind 3CL Pro. The selectivity test showed that the compound had good target selectivity to 3CL Pro over cathepsin L. These correlation assays can prove the rationality of the covalent lead compound screening model. Finally, covalent docking was performed to demonstrate the binding conformation of the compound with 3CL Pro. The source code can be obtained from the GitHub repository (https://github.com/guzh970630/Screen_Covalent_Compound_by_LISNN).

Primary Antioxidant Power and Mpro SARS-CoV-2 Non-Covalent Inhibition Capabilities of Miquelianin.

The antioxidant power of quercetin-3-O-glucuronide (miquelianin) has been studied, at the density functional level of theory, in both lipid-like and aqueous environments. In the aqueous phase, the computed pKa equilibria allowed the identification of the neutral and charged species present in solution that can react with the ⋅OOH radical. The Hydrogen Atom Transfer (HAT), Single Electron Transfer (SET) and Radical Adduct Formation (RAF) mechanisms were considered, and the individual, total and fraction corrected rate constants were obtained. Potential non-covalent inhibition of Mpro from SARS-CoV-2 by miquelianin has been also evaluated.

Structural Basis for Coronaviral Main Proteases Inhibition by the 3CLpro Inhibitor GC376.

The main protease (Mpro) of coronaviruses participates in viral replication, serving as a hot target for drug design. GC376 is able to effectively inhibit the activity of Mpro, which is due to nucleophilic addition of GC376 by binding covalently with Cys145 in Mpro active site. Here, we used fluorescence resonance energy transfer (FRET) assay to analyze the IC50 values of GC376 against Mpros from six different coronaviruses (SARS-CoV-2, HCoV-229E, HCoV-HUK1, MERS-CoV, SARS-CoV, HCoV-NL63) and five Mpro mutants (G15S, M49I, K90R, P132H, S46F) from SARS-CoV-2 variants. The results showed that GC376 displays effective inhibition to various coronaviral Mpros and SARS-CoV-2 Mpro mutants. In addition, the crystal structures of SARS-CoV-2 Mpro (wide type)-GC376, SARS-CoV Mpro-GC376, MERS-CoV Mpro-GC376, and SARS-CoV-2 Mpro mutants (G15S, M49I, S46F, K90R, and P132H)-GC376 complexes were solved. We found that GC376 is able to fit into the active site of Mpros from different coronaviruses and different SARS-CoV-2 variants properly. Detailed structural analysis revealed key molecular determinants necessary for inhibition and illustrated the binding patterns of GC376 to these different Mpros. In conclusion, we not only proved the inhibitory activity of GC376 against different Mpros including SARS-CoV-2 Mpro mutants, but also revealed the molecular mechanism of inhibition by GC376, which will provide scientific guidance for the development of broad-spectrum drugs against SARS-CoV-2 as well as other coronaviruses.

Structural and thermodynamic properties of conserved water molecules in Mpro native: A combined approach by MD simulation and Grid Inhomogeneous Solvation Theory.

The new viral strains of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) are continuously rising, becoming more virulent, and transmissible. Therefore, the development of new antiviral drugs is essential. Due to its significant role in the viral life cycle of SARS-CoV-2, the main protease (Mpro) enzyme is a leading target for antiviral drug design. The Mpro monomer consists of domain DI, DII, and DI-DII interface. Twenty-one conserved water molecules (W4-W24) are occupied at these domains according to multiple crystal structure analyses. The crystal and MD structures reveal the presence of eight conserved water sites in domain DI, DII and remaining in the DI-DII interface. Grid-based inhomogeneous fluid solvation theory (GIST) was employed on MD structures of Mpro native to predict structural and thermodynamic properties of each conserved water site for focusing to identify the specific conserved water molecules that can easily be displaced by proposed ligands. Finally, MD water W13 is emerged as a promising candidate for water mimic drug design due to its low mean interaction energy, loose binding character with the protein, and its involvement in a water-mediated H-bond with catalytic His41 via the interaction Thr25(OG)---W13---W---His41(NE2). In this context, water occupancy, relative interaction energy, entropy, and topologies of W13 are thermodynamically acceptable for the water displacement method. Therefore, the strategic use of W13's geometrical position in the DI domain may be implemented for drug discovery against COVID disease by designing new ligands with appropriately oriented chemical groups to mimic its structural, electronic, and thermodynamic properties.

Open science discovery of potent noncovalent SARS-CoV-2 main protease inhibitors.

We report the results of the COVID Moonshot, a fully open-science, crowdsourced, and structure-enabled drug discovery campaign targeting the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) main protease. We discovered a noncovalent, nonpeptidic inhibitor scaffold with lead-like properties that is differentiated from current main protease inhibitors. Our approach leveraged crowdsourcing, machine learning, exascale molecular simulations, and high-throughput structural biology and chemistry. We generated a detailed map of the structural plasticity of the SARS-CoV-2 main protease, extensive structure-activity relationships for multiple chemotypes, and a wealth of biochemical activity data. All compound designs (>18,000 designs), crystallographic data (>490 ligand-bound x-ray structures), assay data (>10,000 measurements), and synthesized molecules (>2400 compounds) for this campaign were shared rapidly and openly, creating a rich, open, and intellectual property-free knowledge base for future anticoronavirus drug discovery.

Computational study on the mechanisms of inhibition of SARS-CoV-2 Mpro by aldehyde warheads based on DFT.

SARS-CoV-2 main protease, Mpro, plays a crucial role in the virus replication cycle, making it an important target for antiviral research. In this study, a simplified model obtained through truncation is used to explore the reaction mechanism of aldehyde warhead compounds inhibiting Mpro at the level of density functional theory. According to the calculation results, proton transfer (P_T)-nucleophilic attack (N_A) is the rate-determining step in the entire reaction pathway. The water molecule that plays a catalytic role occupies the oxyanion hole, which is unfavorable for the aldehyde warhead to approach the Cys145 SH. Through a hypothetical study of substituting the main chain NH with methylene, it is further confirmed that the P_T-N_A is a proton transfer-dominated process accompanied by a nucleophilic attack reaction. In this process, the oxyanion hole serves only to stabilize the aldehyde oxygen anion and therefore does not have a significant impact on the activation free energy barrier of the step. Our research results provide a unique perspective for understanding the covalent inhibition reaction of the Mpro active site. This study also offers theoretical guidance for the design of new Mpro covalent inhibitors.

Validation of X-ray Crystal Structure Ensemble Representations of SARS-CoV-2 Main Protease by Solution NMR Residual Dipolar Couplings.

Considerable debate has focused on whether sampling of molecular dynamics trajectories restrained by crystallographic data can be used to develop realistic ensemble models for proteins in their natural, solution state. For the SARS-CoV-2 main protease, Mpro, we evaluated agreement between solution residual dipolar couplings (RDCs) and various recently reported multi-conformer and dynamic-ensemble crystallographic models. Although Phenix-derived ensemble models showed only small improvements in crystallographic Rfree, substantially improved RDC agreement over fits to a conventionally refined 1.2-Å X-ray structure was observed, in particular for residues with above average disorder in the ensemble. For a set of six lower resolution (1.55-2.19 Å) Mpro X-ray ensembles, obtained at temperatures ranging from 100 to 310 K, no significant improvement over conventional two-conformer representations was found. At the residue level, large differences in motions were observed among these ensembles, suggesting high uncertainties in the X-ray derived dynamics. Indeed, combining the six ensembles from the temperature series with the two 1.2-Å X-ray ensembles into a single 381-member "super ensemble" averaged these uncertainties and substantially improved agreement with RDCs. However, all ensembles showed excursions that were too large for the most dynamic fraction of residues. Our results suggest that further improvements to X-ray ensemble refinement are feasible, and that RDCs provide a sensitive benchmark in such endeavors. Remarkably, a weighted ensemble of 350 PDB Mpro X-ray structures provided slightly better cross-validated agreement with RDCs than any individual ensemble refinement, implying that differences in lattice confinement also limit the fit of RDCs to X-ray coordinates.