Computational study of the mechanism of a polyurethane esterase A (PueA) from Pseudomonas chlororaphis.

The effective management of plastic waste has become a global imperative, given our reliance on a linear model in which plastics are manufactured, used once, and then discarded. This has led to the pervasive accumulation of plastic debris in landfills and environmental contamination. Recognizing this issue, numerous initiatives are underway to address the environmental repercussions associated with plastic disposal. In this study, we investigate the possible molecular mechanism of polyurethane esterase A (PueA), which has been previously identified as responsible for the degradation of a polyester polyurethane (PU) sample in Pseudomonas chlororaphis, as an effort to develop enzymatic biodegradation solutions. After generating the unsolved 3D structure of the protein by AlphaFold2 from its known genome, the enzymatic hydrolysis of the same model PU compound previously used in experiments has been explored employing QM/MM molecular dynamics simulations. This required a preliminary analysis of the 3D structure of the apo-enzyme, identifying the putative active site, and the search for the optimal protein-substrate binding site. Finally, the resulting free energy landscape indicates that wild-type PueA can degrade PU chains, although with low-level activity. The reaction takes place by a characteristic four-step path of the serine hydrolases, involving an acylation followed by a diacylation step. Energetics and structural analysis of the evolution of the active site along the reaction suggests that PueA can be considered a promising protein scaffold for further development to achieve efficient biodegradation of PU.

Computational Study of the Inhibition of RgpB Gingipain, a Promising Target for the Treatment of Alzheimer's Disease.

Alzheimer's disease represents one of the most ambitious challenges for biomedical sciences due to the growing number of cases worldwide in the elderly population and the lack of efficient treatments. One of the recent attempts to develop a treatment points to the cysteine protease RgpB as a promising drug target. In this attempt, several small-molecule covalent inhibitors of this enzyme have been proposed. Here, we report a computational study at the atomic level of the inhibition mechanism of the most promising reported compounds. Molecular dynamics simulations were performed on six of them, and their binding energies in the active site of the protein were computed. Contact maps and interaction energies were decomposed by residues to disclose those key interactions with the enzyme. Finally, quantum mechanics/molecular mechanics (QM/MM) molecular dynamics (MD) simulations were performed to evaluate the reaction mechanism by which these drug candidates lead to covalently bound complexes, inhibiting the RgpB protease. The results provide a guide for future re-design of prospective and efficient inhibitors for the treatment of Alzheimer's disease.

Mechanistic Modeling of Lys745 Sulfonylation in EGFR C797S Reveals Chemical Determinants for Inhibitor Activity and Discriminates Reversible from Irreversible Agents.

Targeted covalent inhibitors hold promise for drug discovery, particularly for kinases. Targeting the catalytic lysine of epidermal growth factor receptor (EGFR) has attracted attention as a new strategy to overcome resistance due to the emergence of C797S mutation. Sulfonyl fluoride derivatives able to inhibit EGFRL858R/T790M/C797S by sulfonylation of Lys745 have been reported. However, atomistic details of this process are still poorly understood. Here, we describe the mechanism of inhibition of an innovative class of compounds that covalently engage the catalytic lysine of EGFR, through a sulfur(VI) fluoride exchange (SuFEx) process, with the help of hybrid quantum mechanics/molecular mechanics (QM/MM) and path collective variables (PCVs) approaches. Our simulations identify the chemical determinants accounting for the irreversible activity of agents targeting Lys745 and provide hints for the further optimization of sulfonyl fluoride agents.

Assessment of the PETase conformational changes induced by poly(ethylene terephthalate) binding.

Recently, a bacterium strain of Ideonella sakaiensis was identified with the uncommon ability to degrade the poly(ethylene terephthalate) (PET). The PETase from I. sakaiensis strain 201-F6 (IsPETase) catalyzes the hydrolysis of PET converting it to mono(2-hydroxyethyl) terephthalic acid (MHET), bis(2-hydroxyethyl)-TPA (BHET), and terephthalic acid (TPA). Despite the potential of this enzyme for mitigation or elimination of environmental contaminants, one of the limitations of the use of IsPETase for PET degradation is the fact that it acts only at moderate temperature due to its low thermal stability. Besides, molecular details of the main interactions of PET in the active site of IsPETase remain unclear. Herein, molecular docking and molecular dynamics (MD) simulations were applied to analyze structural changes of IsPETase induced by PET binding. Results from the essential dynamics revealed that the ß1-ß2 connecting loop is very flexible. This loop is located far from the active site of IsPETase and we suggest that it can be considered for mutagenesis to increase the thermal stability of IsPETase. The free energy landscape (FEL) demonstrates that the main change in the transition between the unbound to the bound state is associated with the ß7-α5 connecting loop, where the catalytic residue Asp206 is located. Overall, the present study provides insights into the molecular binding mechanism of PET into the IsPETase structure and a computational strategy for mapping flexible regions of this enzyme, which can be useful for the engineering of more efficient enzymes for recycling plastic polymers using biological systems.

QMCube (QM3 ): An all-purpose suite for multiscale QM/MM calculations.

QMCube (QM3 ) is a suite written in the Python programming language, initially focused on multiscale QM/MM simulations of biological systems, but open enough to address other kinds of problems. It allows the user to combine highly efficient QM and MM programs, providing unified access to a wide range of computational methods. The suite also supplies additional modules with extra functionalities. These modules facilitate common tasks such as performing the setup of the models or process the data generated during the simulations. The design of QM3 has been carried out considering the least number of external dependencies (only an algebra library, already included in the distribution), which makes it extremely portable. Also, the modular structure of the suite should help to expand and develop new computational methods.

Unrevealing the Proteolytic Activity of RgpB Gingipain from Computational Simulations.

Alzheimer's disease represents one of the greatest medical concerns for today's population and health services. Its multifactorial inherent nature represents a challenge for its treatment and requires the development of a broad spectrum of drugs. Recently, the cysteine protease gingipain RgpB has been related to neurodegenerative diseases, including Alzheimer's disease, and its inhibition appears to be a promising neuroprotective strategy. Given these features, a computational study that integrates molecular dynamics (MD) simulations with classical and hybrid quantum mechanics/molecular mechanics (QM/MM) potentials was carried out to unravel the atomistic details of RgpB activity. First, a preliminary study based on principal component analysis (PCA), determined the protonation state of the Cys/His catalytic dyad, as well as the crucial role of a flexible loop that favors reactive interactions of the catalytic residues and the peptide in the precatalytic state in its closed conformation. Then, different mechanisms were explored by means of QM/MM MD simulations. The most favorable mechanism consists of two stages. First is an acylation stage that takes place in two steps where, initially, the sulfur atom of the C244 residue attacks the carbonylic carbon of the peptide and the proton of the C244 residue is transferred to the amino group of the peptide in a concerted manner. Subsequently, the peptide bond is broken, and a fragment of the peptide is released. After that, the deacylation stage takes place in a single step where a water molecule attacks the carbonylic carbon of the peptide and a proton of the water is transferred to the C244 residue. The free energy barrier of the rate limiting step is in very good agreement with available experimental data. The mechanism exhibits an unusual role of H211 residue compared with other cysteine proteases but a crucial role of the peptide in triggering the catalysis. Notably, the atomic and energetic particularities found represent a significant contribution to the comprehension of the reaction mechanism and a great opportunity for the design of efficient inhibitors of gingipain RgpB.

Peptide Bond Formation Mechanism Catalyzed by Ribosome.

In this paper we present a study of the peptide bond formation reaction catalyzed by ribosome. Different mechanistic proposals have been explored by means of Free Energy Perturbation methods within hybrid QM/MM potentials, where the chemical system has been described by the M06-2X functional and the environment by means of the AMBER force field. According to our results, the most favorable mechanism in the ribosome would proceed through an eight-membered ring transition state, involving a proton shuttle mechanism through the hydroxyl group of the sugar and a water molecule. This transition state is similar to that described for the reaction in solution (J. Am. Chem. Soc. 2013, 135, 8708-8719), but the reaction mechanisms are noticeably different. Our simulations reproduce the experimentally determined catalytic effect of ribosome that can be explained by the different behavior of the two environments. While the solvent reorganizes during the chemical process involving an entropic penalty, the ribosome is preorganized in the formation of the Michaelis complex and does not suffer important changes along the reaction, dampening the charge redistribution of the chemical system.

QM/MM modeling of the hydroxylation of the androstenedione substrate catalyzed by cytochrome P450 aromatase (CYP19A1).

CYP19A1 aromatase is a member of the Cytochrome P450 family of hemeproteins, and is the enzyme responsible for the final step of the androgens conversion into the corresponding estrogens, via a three-step oxidative process. For this reason, the inhibition of this enzyme plays an important role in the treatment of hormone-dependent breast cancer. The first catalytic subcycle, corresponding to the hydroxilation of androstenedione, has been proposed to occur through a first hydrogen abstraction and a subsequent oxygen rebound step. In present work, we have studied the mechanism of the first catalytic subcycle by means of hybrid quantum mechanics/molecular mechanics methods. The inclusion of the protein flexibility has been achieved by means of Free Energy Perturbation techniques, giving rise to a free energy of activation for the hydrogen abstraction step of 13.5 kcal/mol. The subsequent oxygen rebound step, characterized by a small free energy barrier (1.5 kcal/mol), leads to the hydroxylated products through a highly exergonic reaction. In addition, an analysis of the primary deuterium kinetic isotopic effects, calculated for the hydrogen abstraction step, reveals values (∼10) overpassing the semiclassical limit for the CH, indicating the presence of a substantial tunnel effect. Finally, a decomposition analysis of the interaction energy for the substrate and cofactor in the active site is also discussed. According to our results, the role of the enzymatic environment consists of a transition state stabilization by means of dispersive and polarization effects.

Quantum mechanics/molecular mechanics studies of the mechanism of falcipain-2 inhibition by the epoxysuccinate E64.

Because of the increasing resistance of malaria parasites to antimalarial drugs, the lack of highly effective vaccines, and an inadequate control of mosquito vectors, the problem is growing, especially in the developing world. New approaches to drug development are consequently required. One of the proteases involved in the degradation of human hemoglobin is named falcipain-2 (FP2), which has emerged as a promising target for the development of novel antimalarial drugs. However, very little is known about the inhibition of FP2. In this paper, the inhibition of FP2 by the epoxysuccinate E64 has been studied by molecular dynamics (MD) simulations using hybrid AM1d/MM and M06-2X/MM potentials to obtain a complete picture of the possible free energy reaction paths. A thorough analysis of the reaction mechanism has been conducted to understand the inhibition of FP2 by E64. According to our results, the irreversible attack of Cys42 on E64 can take place on both carbon atoms of the epoxy ring because both processes present similar barriers. While the attack on the C2 atom presents a slightly smaller barrier (12.3 vs 13.6 kcal mol(-1)), the inhibitor-protein complex derived from the attack on C3 appears to be much more stabilized. In contrast to previous hypotheses, our results suggest that residues such as Gln171, Asp170, Gln36, Trp43, Asn81, and even His174 would be anchoring the inhibitor in a proper orientation for the reaction to take place. These results may be useful for the rational design of new compounds with higher inhibitory activity.

Peptidyl nitroalkene inhibitors of main protease rationalized by computational and crystallographic investigations as antivirals against SARS-CoV-2.

The coronavirus disease 2019 (COVID-19) pandemic continues to represent a global public health issue. The viral main protease (Mpro) represents one of the most attractive targets for the development of antiviral drugs. Herein we report peptidyl nitroalkenes exhibiting enzyme inhibitory activity against Mpro (Ki: 1-10 µM) good anti-SARS-CoV-2 infection activity in the low micromolar range (EC50: 1-12 µM) without significant toxicity. Additional kinetic studies of compounds FGA145, FGA146 and FGA147 show that all three compounds inhibit cathepsin L, denoting a possible multitarget effect of these compounds in the antiviral activity. Structural analysis shows the binding mode of FGA146 and FGA147 to the active site of the protein. Furthermore, our results illustrate that peptidyl nitroalkenes are effective covalent reversible inhibitors of the Mpro and cathepsin L, and that inhibitors FGA145, FGA146 and FGA147 prevent infection against SARS-CoV-2.

Role of solvent on nonenzymatic peptide bond formation mechanisms and kinetic isotope effects.

Based on the hypothesis that similar mechanisms are involved in the peptide bond formation in aqueous solution and in the ribosome, the aminolysis of esters in aqueous solution has been the subject of numerous studies as the reference reaction for the catalyzed process. The mechanisms proposed in the literature have been explored in the present paper by hybrid QM/MM molecular dynamics simulations. The free energy profiles have been computed with the QM region of the system described at semiempirical AM1 level and by DFT within the M06-2X functional. According to the results, the formation of adduct zwitterion species is a preliminary step required for all possible mechanisms. Then, from different conformers of this species, four different paths were found: three of them taking place through concerted mechanisms of four-, six- and eight-membered ring transition states, and a stepwise mechanism through a neutral intermediate. Comparison of the free energy profiles indicates that the concerted mechanisms would be kinetically favored, with free energy barriers in very good agreement with experimental data. Calculations of kinetic isotope effects, when including the solute interactions with the first solvation shell, show that the 8-membered ring TS renders values in better agreement with available experimental data. Quantitative discrepancies can be attributed to different employed models in experiments and calculations.

Increased dynamic effects in a catalytically compromised variant of Escherichia coli dihydrofolate reductase.

Isotopic substitution ((15)N, (13)C, (2)H) of a catalytically compromised variant of Escherichia coli dihydrofolate reductase, EcDHFR-N23PP/S148A, has been used to investigate the effect of these mutations on catalysis. The reduction of the rate constant of the chemical step in the EcDHFR-N23PP/S148A catalyzed reaction is essentially a consequence of an increase of the quasi-classical free energy barrier and to a minor extent of an increased number of recrossing trajectories on the transition state dividing surface. Since the variant enzyme is less well set up to catalyze the reaction, a higher degree of active site reorganization is needed to reach the TS. Although millisecond active site motions are lost in the variant, there is greater flexibility on the femtosecond time scale. The "dynamic knockout" EcDHFR-N23PP/S148A is therefore a "dynamic knock-in" at the level of the chemical step, and the increased dynamic coupling to the chemical coordinate is in fact detrimental to catalysis. This finding is most likely applicable not just to hydrogen transfer in EcDHFR but also to other enzymatic systems.

New insight into the electronic structure of iron(IV)-oxo porphyrin compound I. A quantum chemical topological analysis.

The electronic structure of iron-oxo porphyrin π-cation radical complex Por(·+) Fe(IV)=O (S-H) has been studied for doublet and quartet electronic states by means of two methods of the quantum chemical topology analysis: electron localization function (ELF) η(r) and electron density ρ(r). The formation of this complex leads to essential perturbation of the topological structure of the carbon-carbon bonds in porphyrin moiety. The double C=C bonds in the pyrrole anion subunits, represented by pair of bonding disynaptic basins V(i=1,2)(C,C) in isolated porphyrin, are replaced by single attractor V(C,C)(i=1-20) after complexation with the Fe cation. The iron-nitrogen bonds are covalent dative bonds, N→Fe, described by the disynaptic bonding basins V(Fe,N)(i=1-4), where electron density is almost formed by the lone pairs of the N atoms. The nature of the iron-oxygen bond predicted by the ELF topological analysis, shows a main contribution of the electrostatic interaction, Fe(δ+)···O(δ-), as long as no attractors between the C(Fe) and C(O) core basins were found, although there are common surfaces between the iron and oxygen basines and coupling between iron and oxygen lone pairs, that could be interpreted as a charge-shift bond. The Fe-S bond, characterized by the disynaptic bonding basin V(Fe,S), is partially a dative bond with the lone pair donated from sulfur atom. The change of electronic state from the doublet (M = 2) to quartet (M = 4) leads to reorganization of spin polarization, which is observed only for the porphyrin skeleton (-0.43e to 0.50e) and S-H bond (-0.55e to 0.52e).

Computational study of the mechanism of half-reactions in class 1A dihydroorotate dehydrogenase from Trypanosoma cruzi.

Chagas' disease is considered to be a health problem affecting millions of people in Latin America. This disease is caused by the parasite Trypanosoma cruzi. Recently dihydroorotate dehydrogenase class 1A from Trypanosoma cruzi (TcDHODA) was shown to be essential for the survival and growth of T. cruzi and proposed as a drug target against Chagas' disease. This enzyme catalyzes the oxidation of (S)-dihydroorotate to orotate, with a proposed catalytic cycle consisting of two half-reactions. In the first half-reaction dihydroorotate is oxidized to orotate, with the consequent reduction of the flavin mononucleotide cofactor. In the second half-reaction fumarate is reduced to succinate. The first oxidation half-reaction may occur via a concerted or a stepwise mechanism. Herein, the catalytic mechanism of TcDHODA has been studied using hybrid Quantum Mechanical/Molecular Mechanical (QM/MM) Molecular Dynamics (MD) simulations. The free energy profiles derived from the bidimensional potential of mean force reveal more details for two half-reaction processes.

The catalytic mechanism of glyceraldehyde 3-phosphate dehydrogenase from Trypanosoma cruzi elucidated via the QM/MM approach.

Glyceraldehyde-3-phosphate dehydrogenase (GAPDH) has been identified as a key enzyme involved in glycolysis processes for energy production in the Trypanosoma cruzi parasite. This enzyme catalyses the oxidative phosphorylation of glyceraldehyde 3-phosphate (G3P) in the presence of inorganic phosphate (Pi) and nicotinamide adenosine dinucleotide (NAD+). The catalytic mechanism used by GAPDH has been intensively investigated. However, the individual roles of Pi and the C3 phosphate of G3P (Ps) sites, as well as some residues such as His194 in the catalytic mechanism, remain unclear. In this study, we have employed Molecular Dynamics (MD) simulations within hybrid quantum mechanical/molecular mechanical (QM/MM) potentials to obtain the Potential of Mean Force of the catalytic oxidative phosphorylation mechanism of the G3P substrate used by GAPDH. According to our results, the first stage of the reaction (oxidoreduction) takes place in the Pi site (energetically more favourable), with the formation of oxyanion thiohemiacetal and thioacylenzyme intermediates without acid-base assistance of His194. Analysis of the interaction energy by residues shows that Arg249 has an important role in the ability of the enzyme to bind the G3P substrate, which interacts with NAD+ and other important residues, such as Cys166, Glu109, Thr167, Ser247 and Thr226, in the GAPDH active site. Finally, the inhibition mechanism of the GAPDH enzyme by the 3-(p-nitrophenoxycarboxyl)-3-ethylene propyl dihydroxyphosphonate inhibitor was investigated in order to contribute to the design of new inhibitors of GAPDH from Trypanosoma cruzi.

Computational analysis of human OGA structure in complex with PUGNAc and NAG-thiazoline derivatives.

The substitution of serine and threonine residues in nucleocytoplasmic proteins with 2-acetamido-2-deoxy-ß-D-glucopyranose (O-GlcNAc) residues is an essential post-translational modification found in many multicellular eukaryotes. O-glycoprotein 2-acetamino-2-deoxy-ß-D-glucopyranosidase (O-GlcNAcase) hydrolyzes O-GlcNAc residues from post-translationally modified serine/threonine residues of nucleocytoplasmic protein. O-GlcNAc has been implicated in several disease states such as cancer, Alzheimer's disease, and type II diabetes. For this paper, a model of the human O-GlcNAcase (hOGA) enzyme based on the X-ray structures of bacterial Clostridium perfringens (CpNagJ) and Bacteroides thetaiotaomicrometer (BtOGA) homologues has been generated through molecular homology modeling. In addition, molecular docking, molecular dynamics (MD) simulations, and Linear Interaction Energy (LIE) were employed to determine the bind for derivatives of two potent inhibitors: O-(2-acetamido-2-deoxy-D-glucopyranosylidene)amino-N-phenylcarbamate (PUGNAc) and 1,2-dideoxy-2'-methyl-R-D-glucopyranoso-[2,1-d]-Δ2'-thiazoline (NAG-thiazoline), with hOGA. The results show that the binding free energy calculations using the Linear Interaction Energy (LIE) are correlated with inhibition constant values. Therefore, the model of the human O-GlcNAcase (hOGA) obtained here may be used as a target for rational design of new inhibitors.

Theoretical studies of HIV-1 reverse transcriptase inhibition.

Computational methods for accurately calculating the binding affinity of a ligand for a protein play a pivotal role in rational drug design. We herein present a theoretical study of the binding of five different ligands to one of the proteins responsible for the human immunodeficiency virus type 1 (HIV-1) cycle replication; the HIV-1 reverse transcriptase (RT). Two types of approaches are used based on molecular dynamics (MD) simulations within hybrid QM/MM potentials: the alchemical free energy perturbation method, FEP, and the pathway method, in which the ligand is physically pulled away from the binding site, thus rendering a potential of mean force (PMF) for the binding process. Our comparative analysis stresses their advantages and disadvantages and, although the results are not in quantitative agreement, both methods are capable of distinguishing the most and the less potent inhibitors of HIV-1 RT activity on an RNase H site. The methods can then be used to select the proper scaffold to design new drugs. A deeper analysis of these inhibitors through molecular electrostatic potential (MEP) and calculation of the binding contribution of the individual residues shows that, in a rational design, apart from the strong interactions established with the two magnesium cations present in the RNase H site, it is important to take into account interactions with His539 and with those residues that are anchoring the metals; Asp443, Glu478, Asp498 and Asp549. The MEPs of the active site of the protein and the different ligands show a better complementarity in those inhibitors that present higher binding energies, but there are still possibilities of improving the favourable interactions and decreasing those that are repulsive in order to design compounds with higher inhibitory activity.

Understanding the different activities of highly promiscuous MbtI by computational methods.

Salicylate synthase from Mycobacterium tuberculosis, MbtI, is a highly promiscuous Mg(2+) dependent enzyme with up to four distinct activities detected in vitro: isochorismate synthase (IS), isochorismate pyruvate lyase (IPL), salicylate synthase (SS) and chorismate mutase (CM). In this paper, Molecular Dynamic (MD) simulations employing hybrid quantum mechanics/molecular mechanics (QM/MM) potentials have been carried out to get a detailed knowledge of the IS and the IPL activities at the molecular level. According to our simulations, the architecture of the MbtI active site allows catalyzing the two reactions: the isochorismate formation, by means of a stepwise mechanism, and the salicylate production from isochorismate, that appears to be pericyclic in nature. Findings also explain the role of the magnesium cation and the pH dependence activity experimentally observed in MbtI. Mg(2+) would be polarizing and pre-organizing the substrate and active site, as well as shifting the pK(a) values of key active site residues.

Molecular recognition and self-assembly special feature: Multipoint molecular recognition within a calix[6]arene funnel complex.

A multipoint recognition system based on a calix[6]arene is described. The calixarene core is decorated on alternating aromatic subunits by 3 imidazole arms at the small rim and 3 aniline groups at the large rim. This substitution pattern projects the aniline nitrogens toward each other when Zn(II) binds at the Tris-imidazole site or when a proton binds at an aniline. The XRD structure of the monoprotonated complex having an acetonitrile molecule bound to Zn(II) in the cavity revealed a constrained geometry at the metal center reminiscent of an entatic state. Computer modeling suggests that the aniline groups behave as a tritopic monobasic site in which only 1 aniline unit is protonated and interacts with the other 2 through strong hydrogen bonding. The metal complex selectively binds a monoprotonated diamine vs. a monoamine through multipoint recognition: coordination to the metal ion at the small rim, hydrogen bonding to the calix-oxygen core, CH/pi interaction within the cavity's aromatic walls, and H-bonding to the anilines at the large rim.

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.