0000000000873791
AUTHOR
J. Javier Ruiz-pernía
Minimization of dynamic effects in the evolution of dihydrofolate reductase
Protein isotope labeling is a powerful technique to probe functionally important motions in enzyme catalysis and can be applied to investigate the conformational dynamics of proteins.
Critical Role of Substrate Conformational Change in the Proton Transfer Process Catalyzed by 4-Oxalocrotonate Tautomerase
4-Oxalocrotonate tautomerase enzyme (4-OT) catalyzes the isomerization of 2-oxo-4-hexenedioate to 2-oxo-3-hexenedioate. The chemical process involves two proton transfers, one from a carbon of the substrate to the nitrogen of Pro1 and another from this nitrogen atom to a different carbon of the substrate. In this paper the isomerization has been studied using the combined quantum mechanical and molecular mechanical (QM/MM) method with a dual-level treatment of the quantum subsystem employing the MPW1BK density functional as the higher level. Exploration of the potential energy surface shows that the process is stepwise, with a stable intermediate state corresponding to the deprotonated subs…
Heavy enzymes—experimental and computational insights in enzyme dynamics
The role of protein motions in the chemical step of enzyme-catalyzed reactions is the subject of an open debate in the scientific literature. The systematic use of isotopically substituted enzymes has been revealed as a useful tool to quantify the role of these motions. According to the Born-Oppenheimer approximation, changing the mass of the protein does not change the forces acting on the system but alters the frequencies of the protein motions, which in turn can affect the rate constant. Experimental and theoretical studies carried out in this field are presented in this article and discussed in the framework of Transition State Theory.
Temperature dependence of dynamic, tunnelling and kinetic isotope effects in formate dehydrogenase
The origin of the catalytic power of enzymes has been a question of debate for a long time. In this regard, the possible contribution of protein dynamics in enzymatic catalysis has become one of the most controversial topics. In the present work, the hydride transfer step in the formate dehydrogenase (FDH EC 1.2.1.2) enzyme is studied by means of molecular dynamic (MD) simulations with quantum mechanics/molecular mechanics (QM/MM) potentials in order to explore any correlation between dynamics, tunnelling effects and the rate constant. The temperature dependence of the kinetic isotope effects (KIEs), which is one of the few tests that can be studied by experiments and simulations to shed li…
Multiscale Simulations of SARS-CoV-2 3CL Protease Inhibition with Aldehyde Derivatives. Role of Protein and Inhibitor Conformational Dynamics in the Reaction Mechanism
<p>We here investigate the mechanism of SARS-CoV-2 3CL protease inhibition by one of the most promising families of inhibitors, those containing an aldehyde group as warhead. These compounds are covalent inhibitors that inactivate the protease forming a stable hemithioacetal complex. Inhibitor 11a is a potent inhibitor that has been already tested in vitro and in animals. Using a combination of classical and QM/MM simulations we determined the binding mode of the inhibitor into the active site and the preferred rotameric state of the catalytic histidine. In the noncovalent complex the aldehyde group is accommodated into the oxyanion hole formed by the NH main chain groups of residues …
Unraveling the SARS-CoV-2 Main Protease Mechanism Using Multiscale DFT/MM Methods
<p>We present a detailed theoretical analysis of the reaction mechanism of proteolysis catalyzed by the main protease of SARS-CoV-2. Using multiscale simulation methods, we have characterized the interactions stablished by a peptidic substrate in the active site and then we have explored the free energy landscape associated to the acylation and de-acylation steps of the proteolysis reaction, characterizing the transition states of the process. Our mechanistic proposals can explain most of the experimental observations made on the highly similar ortholog protease of SARS-CoV. We point out to some key interactions that may facilitate the acylation process and thus can be crucial in the …
Enzymatic effects on reactant and transition states. The case of chalcone isomerase.
Chalcone isomerase catalyzes the transformation of chalcone to naringerin as a part of flavonoid biosynthetic pathways. The global reaction takes place through a conformational change of the substrate followed by chemical reaction, being thus an excellent example to analyze current theories about enzyme catalysis. We here present a detailed theoretical study of the enzymatic action on the conformational pre-equilibria and on the chemical steps for two different substrates of this enzyme. Free-energy profiles are obtained in terms of potentials of mean force using hybrid quantum mechanics/molecular mechanics potentials. The role of the enzyme becomes clear when compared to the counterpart eq…
A molecular dynamics study on the role of the protonation state in the biosynthesis of R-PAC by AHAS
Abstract The effect of the protonation state of the hydroxyl-ethylthiamin diphosphate intermediate, HEThDP, on the enzyme-substrate interactions and their consequences on the biosynthesis of R-phenylacetylcarbinol, R-PAC, by the acetohydroxy acid synthase, AHAS, is addressed by molecular dynamics simulations. It is found that the form of HEThDP, which favors the formation of R-PAC, is that having the 4-aminopyrimidine ring with the N1′ atom protonated and the N4′ atom as aminopyrimidinium ion. Under this form both active sites of AHAS have the ability to perform the catalysis, unlike that observed for the other possible protonation states of N1′ and N4′ atoms.
A QM/MM study of the reaction mechanism for the 3′-processing step catalyzed by HIV-1 integrase
Integrase (IN) is one of the three human immunodeficiency virus type 1 enzymes (HIV-1) essential for effective viral replication. This viral enzyme is involved in the integration of HIV DNA into host chromosomal DNA. In this work we have carried out molecular dynamics simulations using a hybrid Quantum Mechanics/Molecular Mechanics (QM/MM) approach to study the reaction mechanism for the 3′-processing step of DNA integration using a model substrate. The results obtained by QM(AM1)/MM and QM(PM3)/MM simulations have been improved by single-point corrections using an ab initio method to describe the quantum subsystem. The results obtained within this computational model can be used to obtain …
Hybrid Quantum Mechanics/Molecular Mechanics Simulations with Two-Dimensional Interpolated Corrections: Application to Enzymatic Processes
Hybrid quantum mechanics/molecular mechanics (QM/MM) techniques are widely used to study chemical reactions in large systems. Because of the computational cost associated with the high dimensionality of these systems, the quantum description is usually restricted to low-level methods, such as semiempirical Hamiltonians. In some cases, the description obtained at this computational level is quite poor and corrections must be considered. We here propose a simple but efficient way to include higher-level corrections to be used in potential energy surface explorations and in the calculation of potentials of mean force. We evaluate a correction energy term as the difference between a high-level …
Exploring Chemical Reactivity in Enzyme Catalyzed Processes Using QM/MM Methods: An Application to Dihydrofolate Reductase
Enzymes are the catalysts used by living organisms to accelerate chemical processes under physiological conditions. In this chapter, we illustrate the current view about the origin of their extraordinary rate enhancement based on molecular simulations and, in particular, on methods based on the combination of Quantum Mechanics and Molecular Mechanics potentials which provide a solution to treat the chemical reactivity of these large and complex molecular systems. Computational studies on Dihydrofolate Reductase have been selected as a conductor wire to present the evolution and difficulties to model chemical reactivity in enzymes. The results discussed here show that experimental observatio…
A Microscopic Description of SARS-CoV-2 Main Protease Inhibition with Michael Acceptors. Strategies for Improving Inhibitors 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 …
Increased dynamic effects in a catalytically compromised variant of Escherichia coli dihydrofolate reductase
Isotopic substitution (15N, 13C, 2H) 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 motion…
Hydrolysis of Phosphotriesters: A Theoretical Analysis of the Enzymatic and Solution Mechanisms
A theoretical study on the alkaline hydrolysis of paraoxon, one of the most popular organophosphorus pesticides, in aqueous solution and in the active site of Pseudomonas diminuta phosphotriesterase (PTE) is presented. Simulations by means of hybrid quantum mechanics/molecular mechanics (QM/MM) potentials show that the hydrolysis of paraoxon takes place through an A(N)D(N) or associative mechanism both in solution and in the active site of PTE. The results correctly reproduce the magnitude of the activation free energies and can be used to rationalize the observed kinetic isotope effects (KIEs) for the hydrolysis of paraoxon in both media. Enzymatic hydrolysis of O,O-diethyl p-chlorophenyl …
Computational simulations on the binding and reactivity of a nitrile inhibitor of the SARS-CoV-2 main protease.
We present a detailed computational analysis of the binding mode and reactivity of the novel oral inhibitor PF-07321332 developed against the SARS-CoV-2 3CL protease. Alchemical free energy calculations suggest that positions P3 and P4 could be susceptible to improvement in order to get a larger binding strength. QM/MM simulations unveil the reaction mechanism for covalent inhibition, showing that the nitrile warhead facilitates the recruitment of a water molecule for the proton transfer step.
Unraveling the SARS-CoV-2 Main Protease Mechanism Using Multiscale Methods
We present a detailed theoretical analysis of the reaction mechanism of proteolysis catalyzed by the main protease of SARS-CoV-2. Using multiscale simulation methods, we have characterized the interactions established by a peptidic substrate in the active site, and then we have explored the free energy landscape associated with the acylation and deacylation steps of the proteolysis reaction, characterizing the transition states of the process. Our mechanistic proposals can explain most of the experimental observations made on the highly similar ortholog protease of SARS-CoV. We point to some key interactions that may facilitate the acylation process and thus can be crucial in the design of …
Translocation of enzymes into a mesoporous MOF for enhanced catalytic activity under extreme conditions
Translocation of protease into mesoporous MIL-101-NH2 results in enhanced catalytic activity, excellent recyclability and tolerance to competing enzymes.
Multiscale Simulations of SARS-CoV-2 3CL Protease Inhibition with Aldehyde Derivatives. Role of Protein and Inhibitor Conformational Changes in the Reaction Mechanism
We here investigate the mechanism of SARS-CoV-2 3CL protease inhibition by one of the most promising families of inhibitors, those containing an aldehyde group as a warhead. These compounds are covalent inhibitors that inactivate the protease, forming a stable hemithioacetal complex. Inhibitor 11a is a potent inhibitor that has been already tested in vitro and in animals. Using a combination of classical and QM/MM simulations, we determined the binding mode of the inhibitor into the active site and the preferred rotameric state of the catalytic histidine. In the noncovalent complex, the aldehyde group is accommodated into the oxyanion hole formed by the NH main-chain groups of residues 143 …
Loss of Hyperconjugative Effects Drives Hydride Transfer during Dihydrofolate Reductase Catalysis
Hydride transfer is widespread in nature and has an essential role in applied research. However, the mechanisms of how this transformation occurs in living organisms remain a matter of vigorous debate. Here, we examined dihydrofolate reductase (DHFR), an enzyme that catalyzes hydride from C4′ of NADPH to C6 of 7,8-dihydrofolate (H2F). Despite many investigations of the mechanism of this reaction, the contribution of polarization of the π-bond of H2F in driving hydride transfer remains unclear. H2F was stereospecifically labeled with deuterium β to the reacting center, and β-deuterium kinetic isotope effects were measured. Our experimental results combined with analysis derived from QM/MM si…
A Novel Strategy to Study Electrostatic Effects in Chemical Reactions: Differences between the Role of Solvent and the Active Site of Chalcone Isomerase in a Michael Addition
The electrostatic behavior of active site residues in enzyme catalysis is quite different from that of water molecules in solution. To highlight the electrostatic differences between both environments, we propose a QM/MM strategy to study the role of the environment in chemical reactions. The novelty of the present communication is that free energy surfaces are generated by means of two distinguished reaction coordinates: a solute coordinate and the electrostatic potential created by the environment. This is applied to analyze the origin of catalysis in the transformation of a chalcone into a flavanone, a Michael addition that requires the desolvation of the nucleophile.
Why Are Some Enzymes Dimers? Flexibility and Catalysis in Thermotoga maritima Dihydrofolate Reductase
[Image: see text] Dihydrofolate reductase from Thermotoga maritima (TmDFHFR) is a dimeric thermophilic enzyme that catalyzes the hydride transfer from the cofactor NADPH to dihydrofolate less efficiently than other DHFR enzymes, such as the mesophilic analogue Escherichia coli DHFR (EcDHFR). Using QM/MM potentials, we show that the reduced catalytic efficiency of TmDHFR is most likely due to differences in the amino acid sequence that stabilize the M20 loop in an open conformation, which prevents the formation of some interactions in the transition state and increases the number of water molecules in the active site. However, dimerization provides two advantages to the thermophilic enzyme: …
On the Nature of the Enzyme–Substrate Complex and the Reaction Mechanism in Human Arginase I. A Combined Molecular Dynamics and QM/MM Study
We present here a detailed theoretical analysis of L-arginine hydrolysis catalyzed by Human Arginase I (HARGI). Our study combines classical molecular dynamic simulations of different model for the...
Toward an Automatic Determination of Enzymatic Reaction Mechanisms and Their Activation Free Energies.
We present a combination of the string method and a path collective variable for the exploration of the free energy surface associated to a chemical reaction in condensed environments. The on-the-fly string method is employed to find the minimum free energy paths on a multidimensional free energy surface defined in terms of interatomic distances, which is a convenient selection to study bond forming/breaking processes. Once the paths have been determined, a reaction coordinate is defined as a measure of the advance of the system along these paths. This reaction coordinate can be then used to trace the reaction Potential of Mean Force from which the activation free energy can be obtained. Th…
Unraveling the role of protein dynamics in dihydrofolate reductase catalysis
Protein dynamics have controversially been proposed to be at the heart of enzyme catalysis, but identification and analysis of dynamical effects in enzyme-catalyzed reactions have proved very challenging. Here, we tackle this question by comparing an enzyme with its heavy ((15)N, (13)C, (2)H substituted) counterpart, providing a subtle probe of dynamics. The crucial hydride transfer step of the reaction (the chemical step) occurs more slowly in the heavy enzyme. A combination of experimental results, quantum mechanics/molecular mechanics simulations, and theoretical analyses identify the origins of the observed differences in reactivity. The generally slightly slower reaction in the heavy e…
Comparative computational analysis of different active site conformations and substrates in a chalcone isomerase catalyzed reaction.
Chalcone isomerase catalyzes the transformation of chalcones to flavanones. We present a computational study of the rate-limiting chemical step, an intramolecular Michael addition of a 2'-oxyanion to the alpha,beta-double bound. By using quantum mechanical/molecular mechanical hybrid methods we traced the free-energy profiles associated with the reaction of two different substrates (chalcone and 6'-deoxychalcone) in two different conformations of the active site that are described in the different crystallographic structures available. We have obtained significant differences (about 4 kcal/mol) in the free-energy barriers calculated for the two active sites. According to our results, the ac…
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 sho…
Mechanistic study of the biosynthesis of R-phenylcarbinol by acetohydroxyacid synthase enzyme using hybrid quantum mechanics/molecular mechanics simulations
Abstract The biosynthesis of R-phenylacetylcarbinol (R-PAC) by the acetohydroxy acid synthase, (AHAS) is addressed by molecular dynamics simulations (MD), hybrid quantum mechanics/molecular mechanics (QM/MM), and QM/MM free energy calculations. The results show the reaction starts with the nucleophilic attack of the C2α atom of the HEThDP intermediate on the Cβ atom of the carbonyl group of benzaldehyde substrate via the formation of a transition state (TS1) with the HEThDP intermediate under 4′-aminopyrimidium (APH+) form. The calculated activation free energy for this step is 17.4 kcal mol−1 at 27 °C. From this point, the reaction continues with the abstraction of Hβ atom of the HEThDP in…
Dynamic Effects on Reaction Rates in a Michael Addition Catalyzed by Chalcone Isomerase. Beyond the Frozen Environment Approach
We present a detailed microscopic study of the dynamics of the Michael addition reaction leading from 6'-deoxychalcone to the corresponding flavanone. The reaction dynamics are analyzed for both the uncatalyzed reaction in aqueous solution and the reaction catalyzed by Chalcone Isomerase. By means of rare event simulations of trajectories started at the transition state, we have computed the transmission coefficients, obtaining 0.76 +/- 0.04 and 0.87 +/- 0.03, in water and in the enzyme, respectively. According to these simulations, the Michael addition can be seen as a formation of a new intramolecular carbon-oxygen bond accompanied by a charge transfer essentially taking place from the nu…
Inhibition Mechanism of SARS‐CoV‐2 Main Protease with Ketone‐Based Inhibitors Unveiled by Multiscale Simulations: Insights for Improved Designs**
Abstract We present the results of classical and QM/MM simulations for the inhibition of SARS‐CoV‐2 3CL protease by a hydroxymethylketone inhibitor, PF‐00835231. In the noncovalent complex the carbonyl oxygen atom of the warhead is placed in the oxyanion hole formed by residues 143 to 145, while P1–P3 groups are accommodated in the active site with interactions similar to those observed for the peptide substrate. According to alchemical free energy calculations, the P1′ hydroxymethyl group also contributes to the binding free energy. Covalent inhibition of the enzyme is triggered by the proton transfer from Cys145 to His41. This step is followed by the nucleophilic attack of the Sγ atom on …
Protein isotope effects in dihydrofolate reductase from Geobacillus stearothermophilus show entropic-enthalpic compensatory effects on the rate constant.
Catalysis by dihydrofolate reductase from the moderately thermophilic bacterium Geobacillus stearothermophilus (BsDHFR) was investigated by isotope substitution of the enzyme. The enzyme kinetic isotope effect for hydride transfer was close to unity at physiological temperatures but increased with decreasing temperatures to a value of 1.65 at 5 °C. This behavior is opposite to that observed for DHFR from Escherichia coli (EcDHFR), where the enzyme kinetic isotope effect increased slightly with increasing temperature. These experimental results were reproduced in the framework of variational transition-state theory that includes a dynamical recrossing coefficient that varies with the mass of…