Theoretical insights in enzyme catalysis
In this tutorial review we show how the methods and techniques of computational chemistry have been applied to the understanding of the physical basis of the rate enhancement of chemical reactions by enzymes. This is to answer the question: Why is the activation free energy in enzyme catalysed reactions smaller than the activation free energy observed in solution? Two important points of view are presented: Transition State (TS) theories and Michaelis Complex (MC) theories. After reviewing some of the most popular computational methods employed, we analyse two particular enzymatic reactions: the conversion of chorismate to prephenate catalysed by Bacillus subtilis chorismate mutase, and a m…
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.
Calculation of binding energy using BLYP/MM for the HIV-1 integrase complexed with the S-1360 and two analogues.
Abstract Integrase (IN) is one of the three human immunodeficiency virus type 1 (HIV-1) enzymes essential for effective viral replication. S-1360 is a potent and selective inhibitor of HIV-1 IN. In this work, we have carried out molecular dynamics (MD) simulations using a hybrid Quantum Mechanics/Molecular Mechanics (QM/MM) approach, to determine the protein–ligand interaction energy for S-1360 and two analogues. Analysis of the MD trajectories reveals that the strongest protein–inhibitor interactions, observed in the three studied complexes, are established with Lys-159 residue and Mg 2+ cation. Calculations of binding energy using BLYP/MM level of theory reveal that there is a direct rela…
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–871…
Theoretical Modeling on the Reaction Mechanism of p-Nitrophenylmethylphosphate Alkaline Hydrolysis and its Kinetic Isotope Effects
We have studied the alkaline hydrolysis of p-nitrophenylmethylphosphate (p-NPmP) in aqueous solution by means of polarizable continuum models and by hybrid quantum-mechanical/molecular-mechanical (QM/MM) methods. The theoretical predictions of kinetic isotope effects (KIEs) are in very good agreement with the experimental data, confirming a concerted asynchronous molecular mechanism. In addition, comparison of high level DFT theory with semiempirical AM1/d Hamiltonian has allowed checking the reliability of the later to be used in modeling very large molecular models containing phosphorus atoms.
Promiscuity in alkaline phosphatase superfamily. Unraveling evolution through molecular simulations.
We here present a theoretical study of the alkaline hydrolysis of a phosphodiester (methyl p-nitrophenyl phosphate or MpNPP) in the active site of Escherichia coli alkaline phosphatase (AP), a monoesterase that also presents promiscuous activity as a diesterase. The analysis of our simulations, carried out by means of molecular dynamics (MD) simulations with hybrid quantum mechanics/molecular mechanics (QM/MM) potentials, shows that the reaction takes place through a D(N)A(N) or dissociative mechanism, the same mechanism employed by AP in the hydrolysis of monoesters. The promiscuous activity observed in this superfamily can be then explained on the basis of a conserved reaction mechanism. …
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.
Computational design of biological catalysts
The purpose of this tutorial review is to illustrate the way to design new and powerful catalysts. The first possibility to get a biological catalyst for a given chemical process is to use existing enzymes that catalyze related reactions. The second possibility is the use of immune systems that recognize stable molecules resembling the transition structure of the target reaction. We finally show how computational techniques are able to provide an enormous quantity of information, providing clues to guide the development of new biological catalysts
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…
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 …
Preorganization and reorganization as related factors in enzyme catalysis: the chorismate mutase case.
In this paper a deeper insight into the chorismate-to prephenate-rearrangement, catalyzed by Bacillus subtilis chorismate mutase, is provided by means of a combination of statistical quantum mechanics/molecular mechanics simulation methods and hybrid potential energy surface exploration techniques. The main aim of this work is to present an estimation of the preorganization and reorganization terms of the enzyme catalytic rate enhancement. To analyze the first of these, we have studied different conformational equilibria of chorismate in aqueous solution and in the enzyme active site. Our conclusion is that chorismate mutase preferentially binds the reactive conformer of the substrate--that…
Do dynamic effects play a significant role in enzymatic catalysis? A theoretical analysis of formate dehydrogenase.
A theoretical study of the protein dynamic effects on the hydride transfer between the formate anion and nicotinamide adenine dinucleotide (NAD + ), catalyzed by formate dehydrogenase (FDH), is presented in this paper. The analysis of free downhill molecular dynamic trajectories, performed in the enzyme and compared with the reaction in aqueous solution, has allowed the study of the dynamic coupling between the reacting fragments and the protein or the solvent water molecules, as well as an estimation of the dynamic effect contribution to the catalytic effect from calculation of the transmission coefficient in the enzyme and in solution. The obtained transmission coefficients for the enzyme…
On Transition Structures for Hydride Transfer Step: A Theoretical Study of the Reaction Catalyzed by Dihydrofolate Reductase Enzyme
Abstract A theoretical study is presented of the catalytic mechanism of dihydrofolate reductase (DHFR) enzyme based upon the characterization of the transition structure (TS) for the hydride transfer step. Analytical gradients at AM1 and PM3 semiempirical levels have been used to characterize the saddle point of index one (SPi-1) on global energy hypersurface for the hydride transfer in the active site of DHFR enzyme. The geometry, stereochemistry, electronic structure, and transition vector (TV) components associated to SPi-1 are qualitatively computational level independent. The TV amplitudes show primary and secondary isotope effects to be strongly coupled. The geometrical arrangement of…
Theoretical study of the temperature dependence of dynamic effects in thymidylate synthase.
A theoretical study of the temperature dependence of dynamic effects in the rate limiting step of the reaction catalyzed by thymidylate synthase is presented in this paper. From hybrid Quantum Mechanics/Molecular Mechanics (QM/MM) optimizations of transition state structures within a fully flexible molecular model, free downhill molecular dynamics trajectories have been performed at four different temperatures. The analysis of the reactive and non-reactive trajectories in the enzyme environment has allowed us to study the geometric and electronic coupling between the substrate, the cofactor and the protein. The results show how the contribution of dynamic effects to the rate enhancement mea…
Transition structures for hydride transfer reactions in vacuo and their role in enzyme catalysis
A general discussion as to the role of in vacuo transition structure in enzyme catalysis is presented. Quantum mechanical aspects are emphasized. The transition structures defined as saddle points ...
A theoretical study of carbon-carbon bond formation by a Michael-type addition.
A theoretical study of the Michael-type addition of 1,3-dicarbonyl compounds to α,β-unsaturated carbonyl compounds has been performed in the gas phase by means of the AM1 semiempirical method and by density functional theory (DFT) calculations within the B3LYP and M06-2X hybrid functionals. A molecular model has been selected to mimic the role of a base, which is traditionally used as a catalyst in Michael reactions, an acetate moiety to modulate its basicity, and point charges to imitate the stabilization of the negative charge developed in the substrate during the reaction when taking place in enzymatic environments. Results of the study of six different reactions obtained at the three di…
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 …
QM/MM Determination of Kinetic Isotope Effects for COMT-Catalyzed Methyl Transfer Does Not Support Compression Hypothesis
Secondary alpha-D3 kinetic isotope effects calculated by the hybrid AM1/TIP3P/CHARMM method for the reaction of S-adenosylmethionine with catecholate anion in aqueous solution and catalyzed by rat liver catechol O-methyltransferase at 298 K are 0.94 and 0.85, respectively, in good accord with experiment. The large inverse effect for the enzymatic reaction is not due to compression but arises from significant increases in the stretching and bending force constants involving the isotopically substituted atoms of the transferring methyl group as between the reactant complex and the transition structure, larger than for the reaction in water.
A quantum mechanics/molecular mechanics study of the protein-ligand interaction for inhibitors of HIV-1 integrase.
Human immunodeficiency virus type-1 integrase (HIV-1 IN) is an essential enzyme for effective viral replication. Diketo acids such as L-731,988 and S-1360 are potent and selective inhibitors of HIV-1 IN. In this study, we used molecular dynamics simulations, within the hybrid quantum mechanics/molecular mechanics (QM/MM) approach, to determine the protein-ligand interaction energy between HIV-1 IN and L-731,988 and 10 of its derivatives and analogues. This hybrid methodology has the advantage that it includes quantum effects such as ligand polarisation upon binding, which can be very important when highly polarisable groups are embedded in anisotropic environments, as for example in metal-c…
QM/MM calculations of kinetic isotope effects in the chorismate mutase active site.
Kinetic isotope effects have been computed for the Claisen rearrangement of chorismate to prephenate in aqueous solution and in the active site of chorismate mutase from B. subtilus. These included primary 13C and 18O and secondary 3H effects for substitutions at the bond-making and bond-breaking positions. The initial structures of the putative stationary points on the potential energy surface, required for the calculations of isotope effects using the CAMVIB/CAMISO programs, have been selected from hybrid QM/MM molecular dynamical simulations using the DYNAMO program. Refinement of the reactant complex and transition-state structures has been carried out by means of AM1/CHARMM24/TIP3P cal…
Application of Grote-Hynes theory to the reaction catalyzed by thymidylate synthase.
A theoretical study of dynamic effects on the rate-limiting step of the thymidylate synthase catalyzed reaction has been carried out by means of Grote-Hynes theory, successfully predicting the values of the recrossing effects for a chemical reaction that involves the transfer of a classical light particle. The transmission coefficients, obtained at 278, 293, 303, and 313 K, are almost invariant and in all cases far from unity, revealing a significant coupling of the environment motions and the reaction coordinate. Nevertheless, their energetic contribution to the activation free energy represents less than 0.50 kcal/mol for each of the four tested temperatures. Calculation of the transmissi…
Do zwitterionic species exist in the non-enzymatic peptide bond formation?
The use of proper computational methods and models has allowed answering the controversial question of whether zwitterionic species exist in the mechanism of peptide bond synthesis in aqueous solution. In fact, the different conformations of zwitterionic species open the door to different mechanistic paths.
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…
Modeling methods for studying post-translational and transcriptional modifying enzymes.
Biological catalysis is a complex chemical process that involves not only electronic reorganization in the substrate but also the reorganization of the catalyst. This complexity is even larger in the case of post-transcriptional and post-translational modifications which may involve the interaction between two biomacromolecules. However, the development over the past decades of new computational methods and strategies is offering a detailed molecular picture of the catalytic event and a deep understanding of the mechanisms of chemical reactions in biological environments. Here we review the efforts made in the last years to model catalysis in post-transcriptional and post-translational proc…
Protein Conformational Landscapes and Catalysis. Influence of Active Site Conformations in the Reaction Catalyzed by L-Lactate Dehydrogenase
In the past 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 since the 1980s 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 takes 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 spectrom…
Conformational equilibrium of chorismate. A QM/MM theoretical study combining statistical simulations and geometry optimisations in gas phase and in aqueous solution
We report a theoretical study on the conformational equilibrium of chorismate that precedes its rearrangement to prephenate, an important enzyme-catalyzed reaction. In first place we show that the usual classification of chorismate conformers based on the relative position of the hydroxyl and ether bridge, pseudo-diaxial and pseudo-diequatorial, is not the only relevant factor from the point of view of the a posteriori rearrangement. Here we also analyse another complementary geometrical classification based on the interatomic distance between the carbon atoms to be bounded. Using the umbrella sampling approach and this distance as distinguished internal reaction coordinate, the gas phase A…
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…
A Theoretical Study of the Favorskii Rearrangement. Calculation of Gas-Phase Reaction Paths and Solvation Effects on the Molecular Mechanism for the Transposition of the α-Chlorocyclobutanone
The molecular mechanism of the α-chlorocyclobutanone transposition to yield cyclopropanecarboxylic acid, as a model of the Favorskii rearrangement, has been theoretically characterized in vacuo by means of ab initio molecular orbital procedures at the Hartree−Fock (HF) level of theory with the 6-31G* and 6-31+G* basis sets. The electron correlation has been estimated at the MP2/6-31G* level and calculations based on density functional theory, BLYP/6-31G*. The solvent effects are included at HF/6-31G* level by means of a polarizable continuum model. The questions related to the two accepted molecular mechanisms, the semibenzilic acid and the cyclopropanone transpositions, as well as the comp…
Origin of Enzymatic Kinetic Isotope Effects in Human Purine Nucleoside Phosphorylase
Here we report a study of the effect of heavy isotope labeling on the reaction catalyzed by human purine nucleoside phosphorylase (hPNP) to elucidate the origin of its catalytic effect and of the enzymatic kinetic isotope effect (EKIE). Using quantum mechanical and molecular mechanical (QM/MM) molecular dynamics (MD) simulations, we study the mechanism of the hPNP enzyme and the dynamic effects by means of the calculation of the recrossing transmission coefficient. A free energy surface (FES), as a function of both a chemical and an environmental coordinate, is obtained to show the role of the environment on the chemical reaction. Analysis of reactive and nonreactive trajectories allows us …
An AM1 theoretical study on the effect of Zn2+ Lewis acid catalysis on the mechanism of the cycloaddition between 3-phenyl-1-(2-pyridyl)-2-propen-1-one and cyclopentadiene
Abstract The mechanism of the Diels–Alder reaction between 3-phenyl-1-(2-pyridyl)-2-propen-1-one and cyclopentadiene has been investigated with the AM1 semiempirical method. Stationary points for two reactive channels, endo - cis and exo - cis , have been characterized. The role of the Lewis acid catalyst has been modeled taking into account the formation of a complex between Zn 2+ and the carbonyl oxygen atom and the pyridyl nitrogen atom of the 3-phenyl-1-(2-pyridyl)-2-propen-1-one system with and without the presence of two molecules of water around the cation. The mechanism of the uncatalyzed reaction corresponds to a concerted process, but in the presence of Lewis acid catalyst the mec…
On Transition Structures for Hydride Transfer Step in Enzyme Catalysis. A Comparative Study on Models of Glutathione Reductase Derived from Semiempirical, HF, and DFT Methods
As a model of the chemical reactions that take place in the active site of gluthatione reductase, the nature of the molecular mechanism for the hydride transfer step has been characterized by means of accurate quantum chemical characterizations of transition structures. The calculations have been carried out with analytical gradients at AM1 and PM3 semiempirical procedures, ab initio at HF level with 3-21G, 4-31G, 6-31G, and 6-31G basis sets and BP86 and BLYP as density functional methods. The results of this study suggest that the endo relative orientation on the substrate imposed by the active site is optimal in polarizing the C4-Ht bond and situating the system in the neighborhood of the…
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 …
Computing kinetic isotope effects for chorismate mutase with high accuracy. A new DFT/MM strategy.
A novel procedure has been applied to compute experimentally unobserved intrinsic kinetic isotope effects upon the rearrangement of chorismate to prephenate catalyzed by B. subtilis chorismate mutase. In this modified QM/MM approach, the "low-level" QM description of the quantum region is corrected during the optimization procedure by means of a "high-level" calculation in vacuo, keeping the QM-MM interaction contribution at a quantum "low-level". This allows computation of energies, gradients, and Hessians including the polarization of the QM subsystem and its interaction with the MM environment, both terms calculated using the low-level method at a reasonable computational cost. New infor…
Hybrid Schemes Based on Quantum Mechanics/Molecular Mechanics Simulations
The development of characterization techniques, advanced synthesis methods, as well as molecular modeling has transformed the study of systems in a well-established research field. The current research challenges in biocatalysis and biotransformation evolve around enzyme discovery, design, and optimization. How can we find or create enzymes that catalyze important synthetic reactions, even reactions that may not exist in nature? What is the source of enzyme catalytic power? To answer these and other related questions, the standard strategies have evolved from trial-and-error methodologies based on chemical knowledge, accumulated experience, and common sense into a clearly multidisciplinary …
Dependence of enzyme reaction mechanism on protonation state of titratable residues and QM level description: lactate dehydrogenase
We have studied the dependence of the chemical reaction mechanism of L-lactate dehydrogenase (LDH) on the protonation state of titratable residues and on the level of the quantum mechanical (QM) description by means of hybrid quantum-mechanical/molecular-mechanical (QM/MM) methods; this methodology has allowed clarification of the timing of the hydride transfer and proton transfer components that hitherto had not been possible to state definitively. Ferrer Castillo, Silvia, Silvia.Ferrer@uv.es, Silla Santos, Estanislao, Estanislao.Silla@uv.es ; Tuñon Garcia de Vicuña, Ignacio Nilo, Ignacio.Tunon@uv.es
Towards a Rational Design of Antibody Catalysts through Computational Chemistry
Can hydrophobic interactions be correctly reproduced by the continuum models?
The ability of the continuum models to describe hydrophobic interactions is investigated. In this work we have studied the interactions between two methane molecules in aqueous solution by means of a continuum model. The resulting potential of mean force is in good agreement with those obtained using Monte Carlo and molecular dynamics techniques. The three energy contributions appearing in the continuum energy partition (electrostatic, dispersion−repulsion, and cavitation) have been analyzed. The cavitation free energy plays the most important role of the three, determining the existence of an energy barrier between the contact minimum and the separated methane monomers. This barrier, which…
Theoretical site-directed mutagenesis: Asp168Ala mutant of lactate dehydrogenase
Molecular simulations based on the use of hybrid quantum mechanics/molecular mechanics methods are able to provide detailed information about the complex enzymatic reactions and the consequences of specific mutations on the activity of the enzyme. In this work, the reduction of pyruvate to lactate catalysed by wild-type and Asp168Ala mutant lactate dehydrogenase (LDH) has been studied by means of simulations using a very flexible molecular model consisting of the full tetramer of the enzyme, together with the cofactor NADH, the substrate and solvent water molecules. Our results indicate that the Asp168Ala mutation provokes a shift in the p K a value of Glu199 that becomes unprotonated at n…
A density functional study of flavonoid compounds with anti-HIV activity.
Abstract Quantum chemical calculations at the DFT/B3LYP theory level, with the 6-31G* basis set, was employed to calculate a set of molecular properties of 26 flavonoid compounds with anti-HIV activity. The correlation between biological activity and structural properties was obtained by using the multiple linear regression method. The model obtained showed not only statistical significance but also predictive ability. We demonstrate in this paper that the anti-HIV activity of compounds can be related with the molecular hydrophobicity (ClogP), the electronegativity ( χ ) and the charges on some key atoms, while that the toxicity can be related with the electronic affinities (EA), ClogP and …
Computational strategies for the design of new enzymatic functions
In this contribution, recent developments in the design of biocatalysts are reviewed with particular emphasis in the de novo strategy. Studies based on three different reactions, Kemp elimination, Diels–Alder and Retro-Aldolase, are used to illustrate different success achieved during the last years. Finally, a section is devoted to the particular case of designed metalloenzymes. As a general conclusion, the interplay between new and more sophisticated engineering protocols and computational methods, based on molecular dynamics simulations with Quantum Mechanics/Molecular Mechanics potentials and fully flexible models, seems to constitute the bed rock for present and future successful desig…
Theoretical study of the molecular mechanism of the domino pathways for squarate ester sequential reactions
Are Heme-Dependent Enzymes Always Using a Redox Mechanism? A Theoretical Study of the Kemp Elimination Catalyzed by a Promiscuous Aldoxime Dehydratase
The design of biocatalysts is a goal to improve the rate, selectivity and environmental friendship of chemical processes in biotechnology. In this regard, the use of computational techniques has provided valuable assistance in the design of enzymes with remarkable catalytic activity. In this paper, hybrid QM/MM simulations have allowed getting an insight into the mechanism of a promiscuous aldoxime dehydratase (OxdA) for the Kemp elimination. We first demonstrate that, based on the use of linear response approximation (LRA) methods, the lowest energy electronic state of the benzisoxazole placed in the active sit of OxdA corresponds to a singlet state, being the triplet and the quintet state…
Heavy enzymes and the rational redesign of protein catalysts
Abstract An unsolved mystery in biology concerns the link between enzyme catalysis and protein motions. Comparison between isotopically labelled “heavy” dihydrofolate reductases and their natural‐abundance counterparts has suggested that the coupling of protein motions to the chemistry of the catalysed reaction is minimised in the case of hydride transfer. In alcohol dehydrogenases, unnatural, bulky substrates that induce additional electrostatic rearrangements of the active site enhance coupled motions. This finding could provide a new route to engineering enzymes with altered substrate specificity, because amino acid residues responsible for dynamic coupling with a given substrate present…
ChemInform Abstract: Computational Design of Biological Catalysts
The purpose of this tutorial review is to illustrate the way to design new and powerful catalysts. The first possibility to get a biological catalyst for a given chemical process is to use existing enzymes that catalyze related reactions. The second possibility is the use of immune systems that recognize stable molecules resembling the transition structure of the target reaction. We finally show how computational techniques are able to provide an enormous quantity of information, providing clues to guide the development of new biological catalysts.
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 salicylat…
Theoretical study of phosphodiester hydrolysis in nucleotide pyrophosphatase/phosphodiesterase. Environmental effects on the reaction mechanism.
We here present a theoretical study of the alkaline hydrolysis of methyl p-nitrophenyl phosphate (MpNPP(-)) in aqueous solution and in the active site of nucleotide pyrophosphatase/phosphodiesterase (NPP). The analysis of our simulations, carried out by means of hybrid quantum mechanics/molecular mechanics (QM/MM) methods, shows that the reaction takes place through different reaction mechanisms depending on the environment. Thus, while in aqueous solution the reaction occurs by means of an A(N)D(N) mechanism, the enzymatic process takes place through a D(N)A(N) mechanism. In the first case, we found associative transition-state (TS) structures, while in the enzyme TS structures have dissoc…
Protein Flexibility and Preorganization in the Design of Enzymes. The Kemp Elimination Catalyzed by HG3.17
A recently designed enzyme, HG3.17, obtained by directed evolution, has shown a catalytic activity close to natural enzymes. Hybrid QM/MM molecular dynamics simulations for the Kemp elimination in this new enzyme have provided a deep insight into the origin of its catalytic efficiency. In this case, we have first demonstrated the presence of different conformations with significantly different reactivity. The larger reactivity is related with a better electrostatic preorganization of the environment that creates a more favorable electrostatic potential for the reaction to proceed. In HG3.17, efforts to improve the catalytic properties must be focused in possible mutations increasing the pre…
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: …
Molecular Mechanism of the site-specific self-cleavage of the RNA phosphodiester backbone by a Twister Ribozyme
Published as part of the special collection of articles derived from the 10th Congress on Electronic Structure: Principles and Applications (ESPA-2016). The catalytic activity of some classes of natural RNA, named as ribozymes, has been discovered just in the past decades. In this paper, the cleavage of the RNA phosphodiester backbone has been studied in aqueous solution and in a twister ribozyme from Oryza sativa. The free energy profiles associated with a baseline substrate-assisted mechanism for the reaction in the enzyme and in solution were computed by means of free energy perturbation methods within hybrid QM/MM potentials, describing the chemical system by the M06-2× functional and t…
Binding isotope effects as a tool for distinguishing hydrophobic and hydrophilic binding sites of HIV-1 RT.
The current treatment for HIV-1 infected patients consists of a cocktail of inhibitors, in an attempt to improve the potency of the drugs by adding the possible effects of each supplied compound. In this contribution, nine different inhibitors of HIV-1 RT, one of the three key proteins responsible for the virus replication, have been selected to develop and test a computational protocol that allows getting a deep insight into the inhibitors’ binding mechanism. The interaction between the inhibitors and the protein have been quantified by computing binding free energies through FEP calculations, while a more detailed characterization of the kind of inhibitor–protein interactions is based on …
Theoretical Study of Catalytic Efficiency of a Diels–Alderase Catalytic Antibody: An Indirect Effect Produced During the Maturation Process
The Diels–Alder reaction is one of the most important and versatile transformations available to organic chemists for the construction of complex natural products, therapeutics agents, and synthetic materials. Given the lack of efficient enzymes capable of catalyzing this kind of reaction, it is of interest to ask whether a biological catalyst could be designed from an antibody-combining site. In the present work, a theoretical study of the different behavior of a germline catalytic antibody (CA) and its matured form, 39 A-11, that catalyze a Diels–Alder reaction has been carried out. A free-energy perturbation technique based on a hybrid quantum-mechanics/molecular-mechanics scheme, togeth…
Theoretical Studies of the Self Cleavage Pistol Ribozyme Mechanism
AbstractRibozymes are huge complex biological catalysts composed of a combination of RNA and proteins. Nevertheless, there is a reduced number of small ribozymes, the self-cleavage ribozymes, that are formed just by RNA and, apparently, they existed in cells of primitive biological systems. Unveiling the details of these “fossils” enzymes can contribute not only to the understanding of the origins of life but also to the development of new simplified artificial enzymes. A computational study of the reactivity of the pistol ribozyme carried out by means of classical MD simulations and QM/MM hybrid calculations is herein presented to clarify its catalytic mechanism. Analysis of the geometries…
Transition structure selectivity in enzyme catalysis: a QM/MM study of chorismate mutase
Two different transition structures (TSs) have been located and characterized for the chorismate conversion to prephenate in Bacillus subtilis chorismate mutase by means of hybrid quantum-mechanical/molecular-mechanical (QM/MM) calculations. GRACE software, combined with an AM1/CHARMM24/TIP3P potential, has been used involving full gradient relaxation of the position of ca. 3300 atoms. These TSs have been connected with their respective reactants and products by the intrinsic reaction coordinate (IRC) procedure carried out in the presence of the protein environment, thus obtaining for the first time a realistic enzymatic reaction path for this reaction. Similar QM/MM computational schemes h…
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…
A Quantum Mechanic/Molecular Mechanic Study of the Wild-Type and N155S Mutant HIV-1 Integrase Complexed with Diketo Acid
Integrase (IN) is one of the three human immunodeficiency virus type 1 (HIV-1) enzymes essential for effective viral replication. Recently, mutation studies have been reported that have shown that a certain degree of viral resistance to diketo acids (DKAs) appears when some amino acid residues of the IN active site are mutated. Mutations represent a fascinating experimental challenge, and we invite theoretical simulations for the disclosure of still unexplored features of enzyme reactions. The aim of this work is to understand the molecular mechanisms of HIV-1 IN drug resistance, which will be useful for designing anti-HIV inhibitors with unique resistance profiles. In this study, we use mo…
Enzyme molecular mechanism as a starting point to design new inhibitors: a theoretical study of O-GlcNAcase.
O-Glycoprotein 2-acetamino-2-deoxy-β-d-glucopyranosidase (O-GlcNAcase) hydrolyzes O-linked 2-acetamido-2-deoxy-β-d-glucopyranoside (O-GlcNAc) residues from post-translationally modified serine/threonine residues of nucleocytoplasmic protein. The chemical process involves substrate-assisted catalysis, where two aspartate residues have been identified as the two key catalytic residues of O-GlcNAcase. In this report, the first step of the catalytic mechanism used by O-GlcNAcase involving substrate-assisted catalysis has been studied using a hybrid quantum mechanical/molecular mechanical (QM/MM) Molecular Dynamics (MD) calculations. The free energy profile shows that the formation of the oxazol…
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…
Insights on the origin of catalysis on glycine N-methyltransferase from computational modeling.
The origin of enzyme catalysis remains a question of debate despite much intense study. We report a QM/MM theoretical study of the SN2 methyl transfer reaction catalyzed by a glycine N-methyltransferase (GNMT) and three mutants to test whether recent experimental observations of rate-constant reductions and variations in inverse secondary α-3H kinetic isotope effects (KIEs) should be attributed to changes in the methyl donor−acceptor distance (DAD): is catalysis due to a compression effect? Semiempirical (AM1) and DFT (M06-2X) methods were used to describe the QM subset of atoms, while OPLS-AA and TIP3P classical force fields were used for the protein and water molecules, respectively. The …
Computational study on hydrolysis of cefotaxime in gas phase and in aqueous solution
We are presenting a theoretical study of the hydrolysis of a β-lactam antibiotic in gas phase and in aqueous solution by means of hybrid quantum mechanics/molecular mechanics potentials. After exploring the potential energy surfaces at semiempirical and density functional theory (DFT) level, potentials of mean force have been computed for the reaction in solution with hybrid PM3/TIP3P calculations and corrections with the B3LYP and M06-2X functionals. Inclusion of the full molecule of the antibiotic, Cefotaxime, in the gas phase molecular model has been demonstrated to be crucial since its carboxylate group can activate a nucleophilic water molecule. Moreover, the flexibility of the substra…
Theoretical Study of the Gas Phase Decomposition of Glycolic, Lactic, and 2-Hydroxyisobutyric Acids
The reaction mechanism associated with the decomposition of three α-hydroxycarboxylic acids (glycolic, lactic, and 2-hydroxyisobutyric) in the gas phase to form carbon monoxide, water, and the corresponding carbonyl compounds has been theoretically characterized by using ab initio analytical gradients at the MP2 level of theory with the 6-31G** and 6-31++G** basis sets. A detailed characterization of the potential energy surface points out the existence of three competitive reaction pathways for the decomposition process. The first pathway describes a two-step mechanism, with water elimination and formation of an α-lactone intermediate, achieved by the nucleophilic attack of the carbonylic …
A theoretical study of the addition of CH3MgCl to chiral α-alkoxy carbonyl compounds
Abstract A theoretical characterisation of the addition of CH 3 MgCl to the carbonyl group of 2-hydroxypropanal and 3,4-di- O -methyl-1- O -trimethylsilyl- l -erythrulose has been carried out. The PM3 semiempirical procedure and the ab initio method at HF level of theory with the 3–21G and 6–31G ∗ basis sets have been applied to identify the stationary points on the potential energy surface. The correlation effects have been included by using the MP2/6–31G ∗ method, as well as by means of MP3/6–31G ∗ // MP2/6–31G ∗ and MP4/6–31G ∗ //MP2/6–31G ∗ single point calculations. The geometry, harmonic vibrational frequencies, transition vector and electronic structure of the transition structures h…
Coupling of the guanosine glycosidic bond conformation and the ribonucleotide cleavage reaction: implications for barnase catalysis.
To examine the possible relationship of guanine-dependent GpA conformations with ribonucleotide cleavage, two potential of mean force (PMF) calculations were performed in aqueous solution. In the first calculation, the guanosine glycosidic (Gχ) angle was used as the reaction coordinate, and computations were performed on two GpA ionic species: protonated (neutral) or deprotonated (negatively charged) guanosine ribose O2 ′. Similar energetic profiles featuring two minima corresponding to the anti and syn Gχ regions were obtained for both ionic forms. For both forms the anti conformation was more stable than the syn, and barriers of ∼4 kcal/mol were obtained for the anti → syn transition. Str…
Convergence of Theory and Experiment on the Role of Preorganization, Quantum Tunneling, and Enzyme Motions into Flavoenzyme-Catalyzed Hydride Transfer
Hydride transfer is one of the most common reactions catalyzed by enzymatic systems, and it has become an object of study because of possible significant quantum tunneling effects. In the present work, we provide a combination of theoretical QM/MM simulations and experimental measurements of the rate constants and kinetic isotopic effects (KIEs) for the hydride transfer reaction catalyzed by morphinone reductase, MR. Quantum mechanical tunneling coefficients, computed in the framework of variational transition-state theory, play a significant role in this reaction, reaching values of 23.8 ± 5.5 for the lightest isotopologue—one of the largest values reported for enzymatic systems. This pred…
Catalysis in glycine N-methyltransferase: testing the electrostatic stabilization and compression hypothesis.
Glycine N-methyltransferase (GNMT) is an S-adenosyl-l-methionine dependent enzyme that catalyzes glycine transformation to sarcosine. Here, we present a hybrid quantum mechanics/molecular mechanics (QM/MM) computational study of the reaction compared to the counterpart process in water. The process takes place through an SN2 mechanism in both media with a transition state in which the transferring methyl group is placed in between the donor (SAM) and the acceptor (the amine group of glycine). Comparative analysis of structural, electrostatic, and electronic characteristics of the in-solution and enzymatic transition states allows us to get a deeper insight into the origins of the enzyme's c…
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…
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 dif…
Computational Modeling of Biological Systems: The LDH Story
Lactate dehydrogenases, LDH, catalyzed reaction has been used in this chapter as a conductor wire to present the evolution and difficulties on computing methods to model chemical reactions in enzymes, since the early calculations based at semiempirical level carried out in gas phase to the recent sophisticated simulations based on hybrid Quantum Mechanical/Molecular Mechanics Dynamics (QM/MM MD) schemes. LDH catalyzes the reversible transformation of pyruvate into lactate. The chemical step consists in a hydride and a proton transfer from the cofactor (NADH) and a protonated histidine (His195), respectively. This fact has generated a lot of controversy about the timing of both transfers in …
Revealing the Origin of the Efficiency of the De Novo Designed Kemp Eliminase HG-3.17 by Comparison with the Former Developed HG-3
The design of new biocatalysts is a goal in biotechnology to improve the rate, selectivity and environmental impact of industrial chemical processes. In this regard, the use of computational techniques has provided valuable assistance in the design of new enzymes with remarkable catalytic activity. In this paper, hybrid QM/MM molecular dynamics simulations have allowed insights to be gained on the origin of the limited efficiency of a computationally designed enzyme for the Kemp elimination; the HG-3. Comparison of results derived from this enzyme with those of a more evolved protein containing additional point mutations, HG-3.17, rendered important information that should be taken into acc…
Theoretical Study of Primary Reaction of Pseudozyma antarctica Lipase B as the Starting Point To Understand Its Promiscuity
Pseudozyma antarctica lipase B (PALB) is a serine hydrolase that catalyzes the hydrolysis of carboxylic acid esters in aqueous medium but it has also shown catalytic activity for a plethora of reactions. This promiscuous activity has found widespread applications. In the present paper, the primary reaction of PALB, its native hydrolytic activity, has been studied using hybrid quantum mechanical/molecular mechanical (QM/MM) potentials. Free energy surfaces, obtained from QM/MM Molecular Dynamics (MD) simulations, show that the reaction takes place by means of a multi-step mechanism where the first step, the activation of the carbonyl group of the substrate and the nucleophilic attack of Ser1…
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 com…
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…