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.

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.

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.

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.

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.

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.

A QM/MM study on the origin of retro-aldolase activity of a catalytic antibody.

The retro-aldolase mechanism of methodol catalysed by the catalytic antibody 33F12 is described based on the exploration of the free energy landscape obtained with QM/MM methods. The amino acids involved in the reaction have been identified, as well as their specific role played in the active site and in the flexibility of the loops. Finally, the comparison with a de novo enzyme RA95.5-8F provides a deeper understanding of catalytic differences between such different protein scaffolds.

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.

Theoretical Studies on Mechanism of Inactivation of Kanamycin A by 4'-O-Nucleotidyltransferase.

This work is focused on mechanistic studies of the transfer of an adenylyl group (Adenoside-5'-monophosfate) from adenosine 5'-triphosphate (ATP) to a OH-4' hydroxyl group of an antibiotic. Using hybrid quantum mechanics/molecular mechanics (QM/MM) techniques, we study the substrate and base-assisted mechanisms of the inactivation process of kanamycin A (KAN) catalyzed by 4'-O-Nucleotidyltransferase [ANT(4')], an active enzyme against almost all aminoglycoside antibiotics. Free energy surfaces, obtained with Free Energy Perturbation methods at the M06-2X/MM level of theory, show that the most favorable reaction path presents a barrier of 12.2 kcal·mol-1 that corresponds to the concerted activation of O4' from KAN by Glu145. In addition, the primary and secondary 18O kinetic isotope effects (KIEs) have been computed for bridge O3α, and non-bridge O1α, O2α, and O5' atoms of ATP. The observed normal 1°-KIE of 1.2% and 2°-KIE of 0.07% for the Glu145-assisted mechanism are in very good agreement with experimentally measured data. Additionally, based on the obtained results, the role of electrostatic and compression effects in enzymatic catalysis is discussed.

Theoretical Study of the Mechanism of Exemestane Hydroxylation Catalyzed by Human Aromatase Enzyme.

Human aromatase (CYP19A1) aromatizes the androgens to form estrogens via a three-step oxidative process. The estrogens are necessary in humans, mainly in women, because of the role they play in sexual and reproductive development. However, these also are involved in the development and growth of hormone-dependent breast cancer. Therefore, inhibition of the enzyme aromatase, by means of drugs known as aromatase inhibitors, is the frontline therapy for these types of cancers. Exemestane is a suicidal third-generation inhibitor of aromatase, currently used in breast cancer treatment. In this study, the hydroxylation of exemestane catalyzed by aromatase has been studied by means of hybrid QM/MM methods. The Free Energy Perturbation calculations provided a free energy of activation for the hydrogen abstraction step (rate-limiting step) of 17 kcal/mol. The results reveal that the hydroxylation of exemestane is not the inhibition stage, suggesting a possible competitive mechanism between the inhibitor and the natural substrate androstenedione in the first catalytic subcycle of the enzyme. Furthermore, the analysis of the interaction energy for the substrate and the cofactor in the active site shows that the role of the enzymatic environment during this reaction consists of a transition state stabilization by means of electrostatic effects.

Joint Use of Bonding Evolution Theory and QM/MM Hybrid Method for Understanding the Hydrogen Abstraction Mechanism via Cytochrome P450 Aromatase.

Bonding evolution theory (BET), as a combination of the electron localization function (ELF) and Thom's catastrophe theory (CT), has been coupled with quantum mechanics/molecular mechanics (QM/MM) method in order to study biochemical reaction paths. The evolution of the bond breaking/forming processes and electron pair rearrangements in an inhomogeneous dynamic environment provided by the enzyme has been elucidated. The proposed methodology is applied in an enzymatic system in order to clarify the reaction mechanism for the hydrogen abstraction of the androstenedione (ASD) substrate catalyzed by the cytochrome P450 aromatase enzyme. The use of a QM/MM Hamiltonian allows inclusion of the polarization of the charges derived from the amino acid residues in the wave function, providing a more accurate and realistic description of the chemical process. The hydrogen abstraction step is found to have five different ELF structural stability domains, whereas the C-H breaking and O-H forming bond process rearrangements are taking place in an asynchronous way.

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.

Protein Conformational Landscapes and Catalysis. Influence of Active Site Conformations in the Reaction Catalyzed by L-Lactate Dehydrogenase.

In the last decade L-Lactate Dehydrogenase (LDH) has become an extremely useful marker in both clinical diagnosis and in monitoring the course of many human diseases. It has been assumed from the 80s that the full catalytic process of LDH starts with the binding of the cofactor and the substrate followed by the enclosure of the active site by a mobile loop of the protein before the reaction to take place. In this paper we show that the chemical step of the LDH catalyzed reaction can proceed within the open loop conformation, and the different reactivity of the different protein conformations would be in agreement with the broad range of rate constants measured in single molecule spectrometry studies. Starting from a recently solved X-ray diffraction structure that presented an open loop conformation in two of the four chains of the tetramer, QM/MM free energy surfaces have been obtained at different levels of theory. Depending on the level of theory used to describe the electronic structure, the free energy barrier for the transformation of pyruvate into lactate with the open conformation of the protein varies between 12.9 and 16.3 kcal/mol, after quantizing the vibrations and adding the contributions of recrossing and tunneling effects. These values are very close to the experimentally deduced one (14.2 kcal·mol-1) and ~2 kcal·mol-1 smaller than the ones obtained with the closed loop conformer. Calculation of primary KIEs and IR spectra in both protein conformations are also consistent with our hypothesis and in agreement with experimental data. Our calculations suggest that the closure of the active site is mainly required for the inverse process; the oxidation of lactate to pyruvate. According to this hypothesis H4 type LDH enzyme molecules, where it has been propose that lactate is transformed into pyruvate, should have a better ability to close the mobile loop than the M4 type LDH molecules.

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.

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.

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.

Studying the role of protein dynamics in an SN2 enzyme reaction using free-energy surfaces and solvent coordinates.

Conformational changes are known to be able to drive an enzyme through its catalytic cycle, allowing, for example, substrate binding or product release. However, the influence of protein motions on the chemical step is a controversial issue. One proposal is that the simple equilibrium fluctuations incorporated into transition-state theory are insufficient to account for the catalytic effect of enzymes and that protein motions should be treated dynamically. Here, we propose the use of free-energy surfaces, obtained as a function of both a chemical coordinate and an environmental coordinate, as an efficient way to elucidate the role of protein structure and motions during the reaction. We show that the structure of the protein provides an adequate environment for the progress of the reaction, although a certain degree of flexibility is needed to attain the full catalytic effect. However, these motions do not introduce significant dynamical corrections to the rate constant and can be described as equilibrium fluctuations.

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.

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.