Inhibitor binding influences the protonation states of histidines in SARS-CoV-2 main protease

Prediction of Binding Pose and Affinity of Nelfinavir, a SARS-CoV-2 Main Protease Repositioned Drug, by Combining Docking, Molecular Dynamics, and Fragment Molecular Orbital Calculations

A novel in silico drug design procedure is described targeting the Main protease (Mpro) of the SARS-CoV-2 virus. The procedure combines molecular docking, molecular dynamics (MD), and fragment molecular orbital (FMO) calculations. The binding structure and properties of Mpro were predicted for Nelfinavir (NFV), which had been identified as a candidate compound through drug repositioning, targeting Mpro. Several poses of the Mpro and NFV complexes were generated by docking, from which four docking poses were selected by scoring with FMO energy. Then, each pose was subjected to MD simulation, 100 snapshot structures were sampled from each of the generated MD trajectories, and the structures were evaluated by FMO calculations to rank the pose based on binding energy. Several residues were found to be important in ligand recognition, including Glu47, Asp48, Glu166, Asp187, and Gln189, all of which interacted strongly with NFV. Asn142 is presumably regarded to form hydrogen bonds or CH/π interaction with NFV; however, in the present calculation, their interactions were transient. Moreover, the tert-butyl group of NFV had no interaction with Mpro. Identifying such strong and weak interactions provides candidates for maintaining and substituting ligand functional groups and important suggestions for drug discovery using drug repositioning. Besides the interaction between NFV and the amino acid residues of Mpro, the desolvation effect of the binding pocket also affected the ranking order. A similar procedure of drug design was applied to Lopinavir, and the calculated interaction energy and experimental inhibitory activity value trends were consistent. Our approach provides a new guideline for structure-based drug design starting from a candidate compound whose complex crystal structure has not been obtained.

Structural insights into the modulation of coronavirus spike tilting and infectivity by hinge glycans

Coronavirus spike glycoproteins presented on the virion surface mediate receptor binding, and membrane fusion during virus entry and constitute the primary target for vaccine and drug development. How the structure dynamics of the full-length spikes incorporated in viral lipid envelope correlates with the virus infectivity remains poorly understood. Here we present structures and distributions of native spike conformations on vitrified human coronavirus NL63 (HCoV-NL63) virions without chemical fixation by cryogenic electron tomography (cryoET) and subtomogram averaging, along with site-specific glycan composition and occupancy determined by mass spectrometry. The higher oligomannose glycan shield on HCoV-NL63 spikes than on SARS-CoV-2 spikes correlates with stronger immune evasion of HCoV-NL63. Incorporation of cryoET-derived native spike conformations into all-atom molecular dynamic simulations elucidate the conformational landscape of the glycosylated, full-length spike that reveals a role of hinge glycans in modulating spike bending. We show that glycosylation at N1242 at the upper portion of the stalk is responsible for the extensive orientational freedom of the spike crown. Subsequent infectivity assays implicated involvement of N1242-glyan in virus entry. Our results suggest a potential therapeutic target site for HCoV-NL63.

Discovery of the covalent SARS-CoV-2 Mpro inhibitors from antiviral herbs via integrating target-based high-throughput screening and chemoproteomic approaches

The main proteases (Mpro ) are highly conserved cysteine-rich proteins that can be covalently modified by numerous natural and synthetic compounds. Herein, we constructed an integrative approach to efficiently discover covalent inhibitors of Mpro from complex herbal matrices. This work begins with biological screening of 60 clinically used antiviral herbal medicines, among which Lonicera japonica Flos (LJF) demonstrated the strongest anti-Mpro effect (IC50  = 37.82 μg/mL). Mass spectrometry (MS)-based chemical analysis and chemoproteomic profiling revealed that LJF extract contains at least 50 constituents, of which 22 exhibited the capability to covalently modify Mpro . We subsequently verified the anti-Mpro effects of these covalent binders. Gallic acid and quercetin were found to potently inhibit severe acute respiratory syndrome coronavirus 2 Mpro in dose- and time- dependent manners, with the IC50 values below 10 µM. The inactivation kinetics, binding affinity and binding mode of gallic acid and quercetin were further characterized by fluorescence resonance energy transfer, surface plasmon resonance, and covalent docking simulations. Overall, this study established a practical approach for efficiently discovering the covalent inhibitors of Mpro from herbal medicines by integrating target-based high-throughput screening and MS-based assays, which would greatly facilitate the discovery of key antiviral constituents from medicinal plants.

The molecular mechanism of non-covalent inhibitor WU-04 targeting SARS-CoV-2 3CLpro and computational evaluation of its effectiveness against mainstream coronaviruses

There is an urgent need for highly effective therapeutic agents to interrupt the continued spread of SARS-CoV-2. As a pivotal protease in the replication process of coronaviruses, the 3CLpro protein is considered as a potential target of drug development to stop the spread and infection of the virus. In this work, molecular dynamics (MD) simulations were used to elucidate the molecular mechanism of a novel and highly effective non-covalent inhibitor, WU-04, targeting the SARS-CoV-2 3CLpro protein. The difference in dynamic behavior between the apo-3CLpro and the holo-3CLpro systems suggests that the presence of WU-04 inhibits the motion amplitude of the 3CLpro protein relative to the apo-3CLpro system, thus maintaining a stable conformational binding state. The energy calculations and interaction analysis show that the hot-spot residues Q189, M165, M49, E166, and H41 and the warm-spot residues H163 and C145 have a strong binding capacity to WU-04 by forming multiple hydrogen bonds and hydrophobic interactions, which stabilizes the binding of the inhibitor. After that, the resistance of WU-04 to the six SARS-CoV-2 variants (Alpha, Beta, Gamma, Delta, Lambda, and Omicron) and two other mainstream coronavirus (SARS-CoV and MERS-CoV) 3CLpro proteins was further investigated. Excitingly, the slight difference in energy values relative to the SARS-CoV-2 system indicates that WU-04 is still highly effective against the coronaviruses, which becomes crucial evidence that WU-04 is a pan-inhibitor of the 3CLpro protein in various SARS-CoV-2 variants and other mainstream coronaviruses. The study will hopefully provide theoretical insights for the future rational design and improvement of novel non-covalent inhibitors targeting the 3CLpro protein.

The H163A mutation unravels an oxidized conformation of the SARS-CoV-2 main protease

The main protease of SARS-CoV-2 (Mpro) is an important target for developing COVID-19 therapeutics. Recent work has highlighted Mpro's susceptibility to undergo redox-associated conformational changes in response to cellular and immune-system-induced oxidation. Despite structural evidence indicating large-scale rearrangements upon oxidation, the mechanisms of conformational change and its functional consequences are poorly understood. Here, we present the crystal structure of an Mpro point mutant (H163A) that shows an oxidized conformation with the catalytic cysteine in a disulfide bond. We hypothesize that Mpro adopts this conformation under oxidative stress to protect against over-oxidation. Our metadynamics simulations illustrate a potential mechanism by which H163 modulates this transition and suggest that this equilibrium exists in the wild type enzyme. We show that other point mutations also significantly shift the equilibrium towards this state by altering conformational free energies. Unique avenues of SARS-CoV-2 research can be explored by understanding how H163 modulates this equilibrium.

Comparative Performance of High-Throughput Methods for Protein pKa Predictions

The medically relevant field of protein-based therapeutics has triggered a demand for protein engineering in different pH environments of biological relevance. In silico engineering workflows typically employ high-throughput screening campaigns that require evaluating large sets of protein residues and point mutations by fast yet accurate computational algorithms. While several high-throughput pKa prediction methods exist, their accuracies are unclear due to the lack of a current comprehensive benchmarking. Here, seven fast, efficient, and accessible approaches including PROPKA3, DeepKa, PKAI, PKAI+, DelPhiPKa, MCCE2, and H++ were systematically tested on a nonredundant subset of 408 measured protein residue pKa shifts from the pKa database (PKAD). While no method outperformed the null hypotheses with confidence, as illustrated by statistical bootstrapping, DeepKa, PKAI+, PROPKA3, and H++ had utility. More specifically, DeepKa consistently performed well in tests across multiple and individual amino acid residue types, as reflected by lower errors, higher correlations, and improved classifications. Arithmetic averaging of the best empirical predictors into simple consensuses improved overall transferability and accuracy up to a root-mean-square error of 0.76 pKa units and a correlation coefficient (R2) of 0.45 to experimental pKa shifts. This analysis should provide a basis for further methodological developments and guide future applications, which require embedding of computationally inexpensive pKa prediction methods, such as the optimization of antibodies for pH-dependent antigen binding.

Theoretical insights into the binding interaction of Nirmatrelvir with SARS-CoV-2 Mpro mutants (C145A and C145S): MD simulations and binding free-energy calculation to understand drug resistance

Mpro, the main protease and a crucial enzyme in SARS-CoV-2 is the most fascinating molecular target for pharmacological treatment and is also liable for viral protein maturation. For antiviral therapy, no drugs have been approved clinically to date. Targeting the Mpro with a compound having inhibitory properties against it can hinder viral replication. The therapeutic potential of the antiviral compound Nirmatrelvir (NMV) against SARS-CoV-2 Mpro was investigated using a systematic approach of molecular docking, MD simulations, and binding free energy calculation based on the MM-GBSA method. NMV, a covalent inhibitor with a recently revealed chemical structure, is a promising oral antiviral clinical candidate with significant in vitro anti-SARS-CoV-2 action in third-phase clinical trials. To explore the therapeutic ability and possible drug resistance, the Mpro system was studied for WT and two of its primary mutants (C145A & C145S). The protein-ligand (Mpro/NMV) complexes were further examined through long MD simulations to check the possible drug resistance in the mutants. To understand the binding affinity, the MM-GBSA method was applied to the Mpro/NMV complexes. Moreover, PCA analysis confirms the detachment of the linker region from the major domains in C145S and C145A mutants allowing for conformational alterations in the active-site region. Based on the predicted biological activities and binding affinities of NMV to WT and mutant (C145A & C145S) Mpro, it can be stipulated that NMV may have conventional potency to act as an anti-viral agent against WT Mpro, while the catalytic-dyad mutations may show substantial mutation-induced drug resistance.Communicated by Ramaswamy H. Sarma.

Targeting SARS-CoV-2 Main Protease (MPro) with Kinase Inhibitors: A Promising Approach for Discovering Antiviral and Anti-inflammatory Molecules against SARS-CoV-2

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virus infected over 688 million people worldwide, causing public health concern and approximately 6.8 million deaths due to COVID-19. COVID-19, especially severe cases, is characterized by exacerbated lung inflammation with an increase of pro-inflammatory cytokines. In addition to antiviral drugs, there is a need for anti-inflammatory therapies to treat all phases of COVID-19. One of the most attractive drug targets for COVID-19 is the SARS-CoV-2 main protease (MPro), an enzyme responsible for cleaving polyproteins formed after the translation of viral RNA, which is essential for viral replication. MPro inhibitors, therefore, have the potential to stop viral replication and act as antiviral drugs. Considering that several kinase inhibitors are known for their action in inflammatory pathways, this could also be investigated toward a potential anti-inflammatory treatment for COVID-19. Therefore, the use of kinase inhibitors against SARS-CoV-2 MPro may be a promising strategy to find molecules with dual activity─antiviral and anti-inflammatory. Considering this, the potential of six kinase inhibitors against SARS-CoV-2 MPro were evaluated in silico and in vitro, including Baricitinib, Tofacitinib, Ruxolitinib, BIRB-796, Skepinone-L, and Sorafenib. To assess the inhibitory potential of the kinase inhibitors, a continuous fluorescent-based enzyme activity assay was optimized with SARS-CoV-2 MPro and MCA-AVLQSGFR-K(Dnp)-K-NH2 (substrate). BIRB-796 and Baricitinib were identified as inhibitors of SARS-CoV-2 MPro, presenting IC₅₀ values of 7.99 and 25.31 μM, respectively. As they are also known for their anti-inflammatory action, both are prototype compounds with the potential to present antiviral and anti-inflammatory activity against SARS-CoV-2 infection.

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.

Korupensamine A, but not its atropisomer, korupensamine B, inhibits SARS-CoV-2 in vitro by targeting its main protease (Mpro)

By combining docking and molecular dynamics simulations, we explored a library of 65 mostly axially chiral naphthylisoquinoline alkaloids and their analogues, with most different molecular architectures and structural analogues, for their activity against SARS-CoV-2. Although natural biaryls are often regarded without consideration of their axial chirality, they can bind to protein targets in an atroposelective manner. By combining docking results with steered molecular dynamics simulations, we identified one alkaloid, korupensamine A, that atropisomer-specifically inhibited the main protease (M^pro) activity of SARS-CoV-2 significantly in comparison to the reference covalent inhibitor GC376 (IC₅₀ = 2.52 ± 0.14 and 0.88 ± 0.15 μM, respectively) and reduced viral growth by five orders of magnitude in vitro (EC₅₀ = 4.23 ± 1.31 μM). To investigate the binding pathway and mode of interaction of korupensamine A within the active site of the protease, we utilized Gaussian accelerated molecular dynamics simulations, which reproduced the docking pose of korupensamine A inside the active site of the enzyme. The study presents naphthylisoquinoline alkaloids as a new class of potential anti-COVID-19 agents.

Integrated analyses reveal a hinge glycan regulates coronavirus spike tilting and virus infectivity

Coronavirus spike glycoproteins presented on the virion surface mediate receptor binding, and membrane fusion during virus entry and constitute the primary target for vaccine and drug development. How the structure dynamics of the full-length spikes incorporated in viral lipid envelope correlates with the virus infectivity remains poorly understood. Here we present structures and distributions of native spike conformations on vitrified human coronavirus NL63 (HCoV-NL63) virions without chemical fixation by cryogenic electron tomography (cryoET) and subtomogram averaging, along with site-specific glycan composition and occupancy determined by mass spectroscopy. The higher oligomannose glycan shield on HCoV-NL63 spikes than on SARS-CoV-2 spikes correlates with stronger immune evasion of HCoV-NL63. Incorporation of cryoET-derived native spike conformations into all-atom molecular dynamic simulations elucidate the conformational landscape of the glycosylated, full-length spike that reveals a novel role of stalk glycans in modulating spike bending. We show that glycosylation at N1242 at the upper portion of the stalk is responsible for the extensive orientational freedom of the spike crown. Subsequent infectivity assays support the hypothesis that this glycan-dependent motion impacts virus entry. Our results suggest a potential therapeutic target site for HCoV-NL63.

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.

Methodology-Centered Review of Molecular Modeling, Simulation, and Prediction of SARS-CoV-2

Despite tremendous efforts in the past two years, our understanding of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), virus-host interactions, immune response, virulence, transmission, and evolution is still very limited. This limitation calls for further in-depth investigation. Computational studies have become an indispensable component in combating coronavirus disease 2019 (COVID-19) due to their low cost, their efficiency, and the fact that they are free from safety and ethical constraints. Additionally, the mechanism that governs the global evolution and transmission of SARS-CoV-2 cannot be revealed from individual experiments and was discovered by integrating genotyping of massive viral sequences, biophysical modeling of protein-protein interactions, deep mutational data, deep learning, and advanced mathematics. There exists a tsunami of literature on the molecular modeling, simulations, and predictions of SARS-CoV-2 and related developments of drugs, vaccines, antibodies, and diagnostics. To provide readers with a quick update about this literature, we present a comprehensive and systematic methodology-centered review. Aspects such as molecular biophysics, bioinformatics, cheminformatics, machine learning, and mathematics are discussed. This review will be beneficial to researchers who are looking for ways to contribute to SARS-CoV-2 studies and those who are interested in the status of the field.

EWALD: A macromolecular diffractometer for the second target station

Revealing the positions of all the atoms in large macromolecules is powerful but only possible with neutron macromolecular crystallography (NMC). Neutrons provide a sensitive and gentle probe for the direct detection of protonation states at near-physiological temperatures and clean of artifacts caused by x rays or electrons. Currently, NMC use is restricted by the requirement for large crystal volumes even at state-of-the-art instruments such as the macromolecular neutron diffractometer at the Spallation Neutron Source. EWALD's design will break the crystal volume barrier and, thus, open the door for new types of experiments, the study of grand challenge systems, and the more routine use of NMC in biology. EWALD is a single crystal diffractometer capable of collecting data from macromolecular crystals on orders of magnitude smaller than what is currently feasible. The construction of EWALD at the Second Target Station will cause a revolution in NMC by enabling key discoveries in the biological, biomedical, and bioenergy sciences.

Crystal structure of SARS-CoV 3C-like protease with baicalein

The 3C-like protease (M^pro, 3CL^pro) plays a key role in the replication process in coronaviruses (CoVs). The M^pro is an essential enzyme mediates CoVs replication and is a promising target for development of antiviral drugs. Until now, baicalein has been shown the specific activity for SARS-CoV M^pro in vitro experiments. In this study, we resolved the SARS-CoV M^pro with baicalein by X-ray diffraction at 2.25 Å (PDB code 7XAX), which provided a structural basis for the research and development of baicalein as an anti-CoVs drug.

Evidence for a high pKa of an aspartic acid residue in the active site of CALB by a fully atomistic multiscale approach

Candida antarctica Lipase B (CALB) is a paradigm for the family of lipases. At pH 7, the optimal pH for catalysis, the protonation state of an aspartic acid of the active site (Asp134) could not be conclusively assigned. In fact, the pK_a estimate provided by a widely used computational tool, namely PropKa, that predicts pKa values of ionizable groups in proteins based on the crystallographic structure, is only slightly above 7 (pKa = 7.25). This, along with the lack of an experimental evaluation, makes the assignment of its protonation state at neutral pH challenging. Here, we calculate the pK_a of Asp134 by means of a fully atomistic multiscale computational approach based on classical molecular dynamics (MD) simulation and the perturbed matrix method (PMM), namely the MD-PMM approach. MD-PMM is able to take into account the dynamics of the system and, at the same time, to treat the deprotonation step at the quantum level. The calculations provide a pK_a value of 8.9 ± 1.1, hence suggesting that Asp134 in CALB should be protonated at neutral, and even at slightly basic, pH.

Multiple protonation states in ligand-free SARS-CoV-2 main protease revealed by large-scale quantum molecular dynamics simulations

The main protease (M^pro) in severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) catalyzes the cleavage of polyproteins for viral replication. Here, large-scale quantum molecular dynamics and metadynamics simulations for ligand-free M^pro were performed, where all the atoms were treated quantum-mechanically, focusing on elucidation of the controversial active-site protonation state. The simulations clarified that the interconverting multiple protonation states exist in unliganded M^pro, and the catalytically relevant ion-pair state is more stable than the neutral state, which is consistent with neutron crystallography. The results highlight the importance of the ion-pair state for repurposing or discovering antiviral drugs that target M^pro.

Hit Expansion of a Noncovalent SARS-CoV-2 Main Protease Inhibitor

Inhibition of the SARS-CoV-2 main protease (M^pro) is a major focus of drug discovery efforts against COVID-19. Here we report a hit expansion of non-covalent inhibitors of M^pro. Starting from a recently discovered scaffold (The COVID Moonshot Consortium. Open Science Discovery of Oral Non-Covalent SARS-CoV-2 Main Protease Inhibitor Therapeutics. bioRxiv 2020.10.29.339317) represented by an isoquinoline series, we searched a database of over a billion compounds using a cheminformatics molecular fingerprinting approach. We identified and tested 48 compounds in enzyme inhibition assays, of which 21 exhibited inhibitory activity above 50% at 20 μM. Among these, four compounds with IC₅₀ values around 1 μM were found. Interestingly, despite the large search space, the isoquinolone motif was conserved in each of these four strongest binders. Room-temperature X-ray structures of co-crystallized protein–inhibitor complexes were determined up to 1.9 Å resolution for two of these compounds as well as one of the stronger inhibitors in the original isoquinoline series, revealing essential interactions with the binding site and water molecules. Molecular dynamics simulations and quantum chemical calculations further elucidate the binding interactions as well as electrostatic effects on ligand binding. The results help explain the strength of this new non-covalent scaffold for M^pro inhibition and inform lead optimization efforts for this series, while demonstrating the effectiveness of a high-throughput computational approach to expanding a pharmacophore library.

Fullerenes against COVID-19: Repurposing C60 and C70 to Clog the Active Site of SARS-CoV-2 Protease

The persistency of COVID-19 in the world and the continuous rise of its variants demand new treatments to complement vaccines. Computational chemistry can assist in the identification of moieties able to lead to new drugs to fight the disease. Fullerenes and carbon nanomaterials can interact with proteins and are considered promising antiviral agents. Here, we propose the possibility to repurpose fullerenes to clog the active site of the SARS-CoV-2 protease, M^pro. Through the use of docking, molecular dynamics, and energy decomposition techniques, it is shown that C₆₀ has a substantial binding energy to the main protease of the SARS-CoV-2 virus, M^pro, higher than masitinib, a known inhibitor of the protein. Furthermore, we suggest the use of C₇₀ as an innovative scaffold for the inhibition of SARS-CoV-2 M^pro. At odds with masitinib, both C₆₀ and C₇₀ interact more strongly with SARS-CoV-2 M^pro when different protonation states of the catalytic dyad are considered. The binding of fullerenes to M^pro is due to shape complementarity, i.e., vdW interactions, and is aspecific. As such, it is not sensitive to mutations that can eliminate or invert the charges of the amino acids composing the binding pocket. Fullerenic cages should therefore be more effective against the SARS-CoV-2 virus than the available inhibitors such as masinitib, where the electrostatic term plays a crucial role in the binding.

Energy landscape of the SARS-CoV-2 reveals extensive conformational heterogeneity

Cryo-electron microscopy (cryo-EM) has produced a number of structural models of the SARS-CoV-2 spike, already prompting biomedical outcomes. However, these reported models and their associated electrostatic potential maps represent an unknown admixture of conformations stemming from the underlying energy landscape of the spike protein. As with any protein, some of the spike's conformational motions are expected to be biophysically relevant, but cannot be interpreted only by static models. Using experimental cryo-EM images, we present the energy landscape of the glycosylated spike protein, and identify the diversity of low-energy conformations in the vicinity of its open (so called 1RBD-up) state. The resulting atomic refinement reveal global and local molecular rearrangements that cannot be inferred from an average 1RBD-up cryo-EM model. Here we report varied degrees of "openness" in global conformations of the 1RBD-up state, not revealed in the single-model interpretations of the density maps, together with conformations that overlap with the reported models. We discover how the glycan shield contributes to the stability of these low-energy conformations. Five out of six binding sites we analyzed, including those for engaging ACE2, therapeutic mini-proteins, linoleic acid, two different kinds of antibodies, switch conformations between their known apo- and holo-conformations, even when the global spike conformation is 1RBD-up. This apo-to-holo switching is reminiscent of a conformational preequilibrium. We found only one binding site, namely that of AB-C135 remains in apo state within all the sampled free energy-minimizing models, suggesting an induced fit mechanism for the docking of this antibody to the spike.

Binding mechanism of inhibitors to SARS-CoV-2 main protease deciphered by multiple replica molecular dynamics simulations

The outbreak caused by SARS-CoV-2 has received extensive worldwide attention. As the main protease (M^pro) in SARS-CoV-2 has no human homologues, it is feasible to reduce the possibility of targeting the host protein by accidental drugs. Thus, M^pro has been an attractive target of efficient drug design for anti-SARS-CoV-2 treatment. In this work, multiple replica molecular dynamics (MRMD) simulations, principal component analysis (PCA), free energy landscapes (FELs), and the molecular mechanics-generalized Born surface area (MM-GBSA) method were integrated together to decipher the binding mechanism of four inhibitors masitinib, O6K, FJC and GQU to M^pro. The results indicate that the binding of four inhibitors clearly affects the structural flexibility and internal dynamics of M^pro along with dihedral angle changes of key residues. The analysis of FELs unveils that the stability in the relative orientation and geometric position of inhibitors to M^pro is favorable for inhibitor binding. Residue-based free energy decomposition reveals that the inhibitor-M^pro interaction networks involving hydrogen bonding interactions and hydrophobic interactions provide significant information for the design of potent inhibitors against M^pro. The hot spot residues including H41, M49, F140, N142, G143, C145, H163, H164, M165, E166 and Q189 identified by computational alanine scanning are considered as reliable targets of clinically available inhibitors inhibiting the activities of M^pro.

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

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

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

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

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.

Characterization of the non-covalent interaction between the PF-07321332 inhibitor and the SARS-CoV-2 main protease

We have studied the non-covalent interaction between PF-07321332 and SARS-CoV-2 main protease at the atomic level using a computational approach based on extensive molecular dynamics simulations with explicit solvent. PF-07321332, whose chemical structure has been recently disclosed, is a promising oral antiviral clinical candidate with well-established anti-SARS-CoV-2 activity in vitro. The drug, currently in phase III clinical trials in combination with ritonavir, relies on the electrophilic attack of a nitrile warhead to the catalytic cysteine of the protease. Nonbonded interaction between the inhibitor and the residues of the binding pocket, as well as with water molecules on the protein surface, have been characterized using two different force fields and the two possible protonation states of the main protease catalytic dyad HIS41-CYS145. When the catalytic dyad is in the neutral state, the non-covalent binding is likely to be stronger. Molecular dynamics simulations seems to lend support for an inhibitory mechanism in two steps: a first non-covalent addition with the dyad in neutral form and then the formation of the thiolate-imidazolium ion pair and the ligand relocation for finalising the electrophilic attack.

Crystal structures of human coronavirus NL63 main protease at different pH values

Human coronavirus NL63 (HCoV-NL63), which belongs to the genus Alphacoronavirus, mainly infects children and the immunocompromized and is responsible for a series of clinical manifestations, including cough, fever, rhinorrhoea, bronchiolitis and croup. HCoV-NL63, which was first isolated from a seven-month-old child in 2004, has led to infections worldwide and accounts for 10% of all respiratory illnesses caused by etiological agents. However, effective antivirals against HCoV-NL63 infection are currently unavailable. The HCoV-NL63 main protease (Mpro), also called 3C-like protease (3CLpro), plays a vital role in mediating viral replication and transcription by catalyzing the cleavage of replicase polyproteins (pp1a and pp1ab) into functional subunits. Moreover, Mpro is highly conserved among all coronaviruses, thus making it a prominent drug target for antiviral therapy. Here, four crystal structures of HCoV-NL63 Mpro in the apo form at different pH values are reported at resolutions of up to 1.78 Å. Comparison with Mpro from other human betacoronaviruses such as SARS-CoV-2 and SARS-CoV reveals common and distinct structural features in different genera and extends knowledge of the diversity, function and evolution of coronaviruses.

Characterization of Phytochemicals in Ulva intestinalis L. and Their Action Against SARS-CoV-2 Spike Glycoprotein Receptor-Binding Domain

Coronavirus disease-2019 (COVID-19) has caused a severe impact on almost all aspects of human life and economic development. Numerous studies are being conducted to find novel therapeutic strategies to overcome COVID-19 pandemic in a much effective way. Ulva intestinalis L. (Ui), a marine microalga, known for its antiviral property, was considered for this study to determine the antiviral efficacy against severe acute respiratory syndrome-associated Coronavirus-2 (SARS-CoV-2). The algal sample was dried and subjected to ethanolic extraction, followed by purification and analysis using gas chromatography-coupled mass spectrometry (GC-MS). Forty-three known compounds were identified and docked against the S1 receptor binding domain (RBD) of the spike (S) glycoprotein. The compounds that exhibited high binding affinity to the RBD of S1 protein were further analyzed for their chemical behaviour using conceptual density-functional theory (C-DFT). Finally, pharmacokinetic properties and drug-likeliness studies were carried out to test if the compounds qualified as potential leads. The results indicated that mainly phenols, polyenes, phytosteroids, and aliphatic compounds from the extract, such as 2,4-di-tert-butylphenol (2,4-DtBP), doconexent, 4,8,13-duvatriene-1,3-diol (DTD), retinoyl-β-glucuronide 6',3'-lactone (RBGUL), and retinal, showed better binding affinity to the target. Pharmacokinetic validation narrowed the list to 2,4-DtBP, retinal and RBGUL as the possible antiviral candidates that could inhibit the viral spike protein effectively.

Virus structure and structure-based antivirals

Structure-based antiviral developments in the past two years have been dominated by the structure determination and inhibition of SARS-CoV-2 proteins and new lead molecules for picornaviruses. The SARS-CoV-2 spike protein has been targeted successfully with antibodies, nanobodies, and receptor protein mimics effectively blocking receptor binding or fusion. The two most promising non-structural proteins sharing strong structural and functional conservation across virus families are the main protease and the RNA-dependent RNA polymerase, for which design and reuse of broad range inhibitors already approved for use has been an attractive avenue. For picornaviruses, the increasing recognition of the transient expansion of the capsid as a critical transition towards RNA release has been targeted through a newly identified, apparently widely conserved, druggable, interprotomer pocket preventing viral entry. We summarize some of the key papers in these areas and ponder the practical uses and contributions of molecular modeling alongside empirical structure determination.

Bioguided Isolation of Cyclopenin Analogues as Potential SARS-CoV-2 Mpro Inhibitors from Penicillium citrinum TDPEF34

SARS-CoV-2 virus mutations might increase its virulence, and thus the severity and duration of the ongoing pandemic. Global drug discovery campaigns have successfully developed several vaccines to reduce the number of infections by the virus. However, finding a small molecule pharmaceutical that is effective in inhibiting SARS-CoV-2 remains a challenge. Natural products are the origin of many currently used pharmaceuticals and, for this reason, a library of in-house fungal extracts were screened to assess their potential to inhibit the main viral protease Mpro in vitro. The extract of Penicillium citrinum, TDPEF34, showed potential inhibition and was further analysed to identify potential Mpro inhibitors. Following bio-guided isolation, a series of benzodiazepine alkaloids cyclopenins with good-to-moderate activity against SARS-CoV-2 Mpro were identified. The mode of enzyme inhibition of these compounds was predicted by docking and molecular dynamic simulation. Compounds 1 (isolated as two conformers of S- and R-isomers), 2, and 4 were found to have promising in vitro inhibitory activity towards Mpro, with an IC50 values range of 0.36-0.89 µM comparable to the positive control GC376. The in silico investigation revealed compounds to achieve stable binding with the enzyme active site through multiple H-bonding and hydrophobic interactions. Additionally, the isolated compounds showed very good drug-likeness and ADMET properties. Our findings could be utilized in further in vitro and in vivo investigations to produce anti-SARS-CoV-2 drug candidates. These findings also provide critical structural information that could be used in the future for designing potent Mpro inhibitors.

SARS-COV-2 Mpro conformational changes induced by covalently bound ligands

SARS-CoV-2's main protease (Mpro) interaction with ligands has been explored with a myriad of crystal structures, most of the monomers. Nonetheless, Mpro is known to be active as a dimer but the relevance of the dimerization in the ligand-induced conformational changes has not been fully elucidated. We systematically simulated different Mpro-ligand complexes aiming to study their conformational changes and interactions, through molecular dynamics (MD). We focused on covalently bound ligands (N1 and N3, ∼9 μs per system both monomers and dimers) and compared these trajectories against the apostructure. Our results suggest that the monomeric simulations led to an unrealistically flexible active site. In contrast, the Mpro dimer displayed a stable oxyanion-loop conformation along the trajectory. Also, ligand interactions with residues His41, Gly143, His163, Glu166 and Gln189 are postulated to impact the ligands' inhibitory activity significantly. In dimeric simulations, especially Gly143 and His163 have increased interaction frequencies. In conclusion, long-timescale MD is a more suitable tool for exploring in silico the activity of bioactive compounds that potentially inhibit the dimeric form of SARS-CoV-2 Mpro.Communicated by Ramaswamy H. Sarma.

The Se-S Bond Formation in the Covalent Inhibition Mechanism of SARS-CoV-2 Main Protease by Ebselen-like Inhibitors: A Computational Study

The inhibition mechanism of the main protease (Mpro) of SARS-CoV-2 by ebselen (EBS) and its analog with a hydroxyl group at position 2 of the benzisoselenazol-3(2H)-one ring (EBS-OH) was studied by using a density functional level of theory. Preliminary molecular dynamics simulations on the apo form of Mpro were performed taking into account both the hydrogen donor and acceptor natures of the Nδ and Nε of His41, a member of the catalytic dyad. The potential energy surfaces for the formation of the Se-S covalent bond mediated by EBS and EBS-OH on Mpro are discussed in detail. The EBS-OH shows a distinctive behavior with respect to EBS in the formation of the noncovalent complex. Due to the presence of canonical H-bonds and noncanonical ones involving less electronegative atoms, such as sulfur and selenium, the influence on the energy barriers and reaction energy of the Minnesota hybrid meta-GGA functionals M06, M06-2X and M08HX, and the more recent range-separated hybrid functional wB97X were also considered. The knowledge of the inhibition mechanism of Mpro by the small protease inhibitors EBS or EBS-OH can enlarge the possibilities for designing more potent and selective inhibitor-based drugs to be used in combination with other antiviral therapies.

In Silico Characterization of Masitinib Interaction with SARS-CoV-2 Main Protease

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection continues to be a global health problem. Despite the current implementation of COVID-19 vaccination schedules, identifying effective antiviral drug treatments for this disease continues to be a priority. A recent study showed that masitinib (MST), a tyrosine kinase inhibitor, blocks the proteolytic activity of SARS-CoV-2 main protease (Mpro ). Although MST is a potential candidate for COVID-19 treatment, a comprehensive analysis of its interaction with Mpro has not been done. In this work, we performed molecular dynamics simulations of the MST-Mpro complex crystal structure. The effect of the protonation states of Mpro H163 residue and MST titratable groups were studied. Furthermore, we identified the MST substituents and Mpro mutations that affect the stability of the complex. Our results provide valuable insights into the design of new MST analogs as potential treatments for COVID-19.

Strengths and Weaknesses of Docking Simulations in the SARS-CoV-2 Era: the Main Protease (Mpro) Case Study

The scientific community is working against the clock to arrive at therapeutic interventions to treat patients with COVID-19. Among the strategies for drug discovery, virtual screening approaches have the capacity to search potential hits within millions of chemical structures in days, with the appropriate computing infrastructure. In this article, we first analyzed the published research targeting the inhibition of the main protease (Mpro), one of the most studied targets of SARS-CoV-2, by docking-based methods. An alarming finding was the lack of an adequate validation of the docking protocols (i.e., pose prediction and virtual screening accuracy) before applying them in virtual screening campaigns. The performance of the docking protocols was tested at some level in 57.7% of the 168 investigations analyzed. However, we found only three examples of a complete retrospective analysis of the scoring functions to quantify the virtual screening accuracy of the methods. Moreover, only two publications reported some experimental evaluation of the proposed hits until preparing this manuscript. All of these findings led us to carry out a retrospective performance validation of three different docking protocols, through the analysis of their pose prediction and screening accuracy. Surprisingly, we found that even though all tested docking protocols have a good pose prediction, their screening accuracy is quite limited as they fail to correctly rank a test set of compounds. These results highlight the importance of conducting an adequate validation of the docking protocols before carrying out virtual screening campaigns, and to experimentally confirm the predictions made by the models before drawing bold conclusions. Finally, successful structure-based drug discovery investigations published during the redaction of this manuscript allow us to propose the inclusion of target flexibility and consensus scoring as alternatives to improve the accuracy of the methods.

Neutron crystallography for the elucidation of enzyme catalysis

Hydrogen atoms and hydration water molecules in proteins are indispensable for many biochemical processes, especially enzymatic catalysis. The locations of hydrogen atoms in proteins are usually predicted based on X-ray structures, but it is still very difficult to know the ionization states of the catalytic residues, the hydration structure of the protein, and the characteristics of hydrogen-bonding interactions. Neutron crystallography allows the direct observation of hydrogen atoms that play crucial roles in molecular recognition and the catalytic reactions of enzymes. In this review, we present the current status of neutron crystallography in structural biology and recent neutron structural analyses of three enzymes: ascorbate peroxidase, the main protease of severe acute respiratory syndrome coronavirus 2, and copper-containing nitrite reductase.

Scaffold Hopping of α-Rubromycin Enables Direct Access to FDA-Approved Cromoglicic Acid as a SARS-CoV-2 MPro Inhibitor

The COVID-19 pandemic is still active around the globe despite the newly introduced vaccines. Hence, finding effective medications or repurposing available ones could offer great help during this serious situation. During our anti-COVID-19 investigation of microbial natural products (MNPs), we came across α-rubromycin, an antibiotic derived from Streptomyces collinus ATCC19743, which was able to suppress the catalytic activity (IC50 = 5.4 µM and Ki = 3.22 µM) of one of the viral key enzymes (i.e., MPro). However, it showed high cytotoxicity toward normal human fibroblasts (CC50 = 16.7 µM). To reduce the cytotoxicity of this microbial metabolite, we utilized a number of in silico tools (ensemble docking, molecular dynamics simulation, binding free energy calculation) to propose a novel scaffold having the main pharmacophoric features to inhibit MPro with better drug-like properties and reduced/minimal toxicity. Nevertheless, reaching this novel scaffold synthetically is a time-consuming process, particularly at this critical time. Instead, this scaffold was used as a template to explore similar molecules among the FDA-approved medications that share its main pharmacophoric features with the aid of pharmacophore-based virtual screening software. As a result, cromoglicic acid (aka cromolyn) was found to be the best hit, which, upon in vitro MPro testing, was 4.5 times more potent (IC50 = 1.1 µM and Ki = 0.68 µM) than α-rubromycin, with minimal cytotoxicity toward normal human fibroblasts (CC50 > 100 µM). This report highlights the potential of MNPs in providing unprecedented scaffolds with a wide range of therapeutic efficacy. It also revealed the importance of cheminformatics tools in speeding up the drug discovery process, which is extremely important in such a critical situation.

Influence of charge configuration on substrate binding to SARS-CoV-2 main protease

While the state-of-the-art computational simulations support the neutral state for the catalytic dyad of the SARS-CoV-2 main protease, the recently-reported neutron structure exhibits a zwitterionic form. To better compare the structural and dynamical features of the two charge configurations, we perform a Molecular Dynamics study of the dimeric enzyme in complex with a peptide substrate. The simulations show that the enzyme charge configuration from the neutron structure is not compatible with a catalytically-competent binding mode for peptide substrates.

Inhibition mechanism and hot-spot prediction of nine potential drugs for SARS-CoV-2 Mpro by large-scale molecular dynamic simulations combined with accurate binding free energy calculations

Coronavirus disease 2019 (COVID-19), which is caused by a new coronavirus known as severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), is spreading around the world. However, a universally effective treatment regimen has not been developed to date. The main protease (M^pro), a key enzyme of SARS-CoV-2, plays a crucial role in the replication and transcription of this virus in cells and has become the ideal target for rational antiviral drug design. In this study, we performed molecular dynamics simulations three times for these complexes of M^pro (monomeric and dimeric) and nine potential drugs that have a certain effect on the treatment of COVID-19 to explore their binding mechanism. In addition, a total of 12 methods for calculating binding free energy were employed to determine the optimal drug. Ritonavir, Arbidol, and Chloroquine consistently showed an outstanding binding ability to monomeric M^pro under various methods. Ritonavir, Arbidol, and Saquinavir presented the best performance when binding to a dimer, which was independent of the protonated state of Hie41 (protonated at N_ε) and Hid41 (protonated at N_δ), and these findings suggest that Chloroquine may not effectively inhibit the activity of dimeric M^proin vivo. Furthermore, three common hot-spot residues of Met165, Hie41, and Gln189 of monomeric M^pro systems dominated the binding of Ritonavir, Arbidol, and Chloroquine. In dimeric M^pro, Gln189, Met165, and Met49 contributed significantly to binding with Ritonavir, Arbidol, and Saquinavir; therefore, Gln189 and Met165 might serve as the focus in the discovery and development of anti-COVID-19 drugs. In addition, the van der Waals interaction played a significant role in the binding process, and the benzene ring of the drugs showed an apparent inhibitory effect on the normal function of M^pro. The binding cavity had great flexibility to accommodate these different drugs. The results would be notably helpful for enabling a detailed understanding of the binding mechanisms for these important drug-M^pro interactions and provide valuable guidance for the design of potent inhibitors.

Tuning Proton Transfer Thermodynamics in SARS-CoV-2 Main Protease: Implications for Catalysis and Inhibitor Design

The catalytic reaction in SARS-CoV-2 main protease is activated by a proton transfer (PT) from Cys145 to His41. The same PT is likely also required for the covalent binding of some inhibitors. Here we use a multiscale computational approach to investigate the PT thermodynamics in the apo enzyme and in complex with two potent inhibitors, N3 and the α-ketoamide 13b. We show that with the inhibitors the free energy cost to reach the charge-separated state of the active-site dyad is lower, with N3 inducing the most significant reduction. We also show that a few key sites (including specific water molecules) significantly enhance or reduce the thermodynamic feasibility of the PT reaction, with selective desolvation of the active site playing a crucial role. The approach presented is a cost-effective procedure to identify the enzyme regions that control the activation of the catalytic reaction and is thus also useful to guide the design of inhibitors.

Olive-Derived Triterpenes Suppress SARS COV-2 Main Protease: A Promising Scaffold for Future Therapeutics

SARS CoV-2 pandemic is still considered a global health disaster, and newly emerged variants keep growing. A number of promising vaccines have been recently developed as a protective measure; however, cost-effective treatments are also of great importance to support this critical situation. Previously, betulinic acid has shown promising antiviral activity against SARS CoV via targeting its main protease. Herein, we investigated the inhibitory potential of this compound together with three other triterpene congeners (i.e., ursolic acid, maslinic acid, and betulin) derived from olive leaves against the viral main protease (Mpro) of the currently widespread SARS CoV-2. Interestingly, betulinic, ursolic, and maslinic acids showed significant inhibitory activity (IC50 = 3.22-14.55 µM), while betulin was far less active (IC50 = 89.67 µM). A comprehensive in-silico analysis (i.e., ensemble docking, molecular dynamic simulation, and binding-free energy calculation) was then performed to describe the binding mode of these compounds with the enzyme catalytic active site and determine the main essential structural features required for their inhibitory activity. Results presented in this communication indicated that this class of compounds could be considered as a promising lead scaffold for developing cost-effective anti-SARS CoV-2 therapeutics.

Direct Observation of Protonation State Modulation in SARS-CoV-2 Main Protease upon Inhibitor Binding with Neutron Crystallography

The main protease (3CL Mpro) from severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the virus that causes COVID-19, is an essential enzyme for viral replication with no human counterpart, making it an attractive drug target. To date, no small-molecule clinical drugs are available that specifically inhibit SARS-CoV-2 Mpro. To aid rational drug design, we determined a neutron structure of Mpro in complex with the α-ketoamide inhibitor telaprevir at near-physiological (22 °C) temperature. We directly observed protonation states in the inhibitor complex and compared them with those in the ligand-free Mpro, revealing modulation of the active-site protonation states upon telaprevir binding. We suggest that binding of other α-ketoamide covalent inhibitors can lead to the same protonation state changes in the Mpro active site. Thus, by studying the protonation state changes induced by inhibitors, we provide crucial insights to help guide rational drug design, allowing precise tailoring of inhibitors to manipulate the electrostatic environment of SARS-CoV-2 Mpro.

A microscopic description of SARS-CoV-2 main protease inhibition with Michael acceptors. Strategies for improving inhibitor design

The irreversible inhibition of the main protease of SARS-CoV-2 by a Michael acceptor known as N3 has been investigated using multiscale methods. The noncovalent enzyme-inhibitor complex was simulated using classical molecular dynamics techniques and the pose of the inhibitor in the active site was compared to that of the natural substrate, a peptide containing the Gln-Ser scissile bond. The formation of the covalent enzyme-inhibitor complex was then simulated using hybrid QM/MM free energy methods. After binding, the reaction mechanism was found to be composed of two steps: (i) the activation of the catalytic dyad (Cys145 and His41) to form an ion pair and (ii) a Michael addition where the attack of the Sγ atom of Cys145 to the Cβ atom of the inhibitor precedes the water-mediated proton transfer from His41 to the Cα atom. The microscopic description of protease inhibition by N3 obtained from our simulations is strongly supported by the excellent agreement between the estimated activation free energy and the value derived from kinetic experiments. Comparison with the acylation reaction of a peptide substrate suggests that N3-based inhibitors could be improved by adding chemical modifications that could facilitate the formation of the catalytic dyad ion pair.