0000000000486186
AUTHOR
Ian H. Williams
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.
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…
The transition state and cognate concepts
Abstract This review aims firstly to clarify the meanings of key terms and concepts associated with the idea of the transition state, as developed by theoreticians and applied by experimentalist, and secondly to provide an update to the meaning and significance of the transition state in an era when computational simulation, in which complexity is being increasingly incorporated, is commonly employed as a means by which to bridge the realms of theory and experiment. The relationship between the transition state and the potential-energy surface for an elementary reaction is explored, with discussion of the following terms: saddle point, minimum-energy reaction path, reaction coordinate, acti…
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…
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
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…
Computational simulation of the lifetime of the methoxymethyl cation in water. A simple model for a glycosyl cation: when is an intermediate an intermediate?
A two-dimensional free-energy surface is constructed for transfer of the methoxymethyl cation between two water molecules. These atoms are treated quantum mechanically within a box of >1000 classical solvent water molecules, and the molecular dynamics of the whole system is considered at 300 K. This provides a simple model for glycosyl transfer in water. The best surface obtained (MPWB1K/6-31+G(d,p) corrected AMI/TIP3P) contains a shallow free-energy well corresponding to an oxacarbenium ion intermediate in a stepwise mechanism. Molecular dynamics analysis at three temperatures leads to a classical estimate of the lifetime of the methoxymethyl cation in water; when quantum corrections fo…
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 …
Does glycosyl transfer involve an oxacarbenium intermediate? Computational simulation of the lifetime of the methoxymethyl cation in water
2D free-energy surfaces for transfer of the methoxymethyl cation between two water molecules are constructed from molecular dynamics (MD) simulations in which these atoms are treated quantum-mechanically within a box of 1030 classical solvent water molecules at 300 K. This provides a simple model for glycosyl transfer in water. The AM1/TIP3P surfaces with 2D-spline corrections at either MPWB1K/6-31+G(d,p) or MP2/6-31+G(d,p) contain a shallow free-energy well corresponding to an oxacarbenium ion intermediate in a DN*AN mechanism. MD analysis at three temperatures leads to a classical estimate of the lifetime of the methoxymethyl cation in water; when quantum corrections for vibrational zero-…