Small temperature dependence of the kinetic isotope effect for the hydride transfer reaction catalyzed by Escherichia coli dihydrofolate reductase

The H/D primary kinetic isotope effect (KIE) for the hydride transfer reaction catalyzed by Escherichia coli dihydrofolate reductase (ecDHFR) is calculated as a function of temperature employing ensemble-averaged variational transition-state theory with multidimensional tunneling. The calculated KIEs display only a small temperature dependence over the temperature range of 5 to 45 degrees C. We identify two key features that contribute to canceling most of the temperature dependence of the KIE that would be expected on the basis of simpler models. Related issues such as the isotope effects on Arrhenius preexponential factors, large differences between free energies of activation and Arrhenius activation energy, and fluctuations of effective barriers are also discussed.     

Studies on the temperature dependence of isotope effects arising in low-temperature acetaldehyde photolysis

The temperature dependence of hydrogen isotope effects arising in the photolysis of mixtures of acetaldehyde-h₄ and -d₄ has been studied experimentally. The temperature maximum of these effects (ξ = 40 at T 27 K) has been found. Good agreement of theory and experiment has been obtained in the approximation of strong exciton-phonon coupling.     

Magnesium isotope effects in enzymatic phosphorylation

Recent discovery of magnesium isotope effect in the rate of enzymatic synthesis of adenosine triphosphate (ATP) offers a new insight into the mechanochemistry of enzymes as the molecular machines. The activity of phosphorylating enzymes (ATP-synthase, phosphocreatine, and phosphoglycerate kinases) in which Mg(2+) ion has a magnetic isotopic nucleus 25Mg was found to be 2-3 times higher than that of enzymes in which Mg(2+) ion has spinless, nonmagnetic isotopic nuclei 24Mg or 26Mg. This isotope effect demonstrates unambiguously that the ATP synthesis is a spin-dependent ion-radical process. The reaction schemes, suggested to explain the effect, imply a reversible electron transfer from the terminal phosphate anion of ADP to Mg(2+) ion as a first step, generating ion-radical pair with singlet and triplet spin states. The yields of ATP along the singlet and triplet channels are controlled by hyperfine coupling of unpaired electron in 25Mg+ ion with magnetic nucleus 25Mg. There is no difference in the ATP yield for enzymes with 24Mg and 26Mg; it gives evidence that in this reaction magnetic isotope effect (MIE) operates rather than classical, mass-dependent one. Similar effects have been also found for the pyruvate kinase. Magnetic field dependence of enzymatic phosphorylation is in agreement with suggested ion-radical mechanism.     

Proton-coupled electron transfer in soybean lipoxygenase: dynamical behavior and temperature dependence of kinetic isotope effects

The dynamical behavior and the temperature dependence of the kinetic isotope effects (KIEs) are examined for the proton-coupled electron transfer reaction catalyzed by the enzyme soybean lipoxygenase. The calculations are based on a vibronically nonadiabatic formulation that includes the quantum mechanical effects of the active electrons and the transferring proton, as well as the motions of all atoms in the complete solvated enzyme system. The rate constant is represented by the time integral of a probability flux correlation function that depends on the vibronic coupling and on time correlation functions of the energy gap and the proton donor-acceptor mode, which can be calculated from classical molecular dynamics simulations of the entire system. The dynamical behavior of the probability flux correlation function is dominated by the equilibrium protein and solvent motions and is not significantly influenced by the proton donor-acceptor motion. The magnitude of the overall rate is strongly influenced by the proton donor-acceptor frequency, the vibronic coupling, and the protein/solvent reorganization energy. The calculations reproduce the experimentally observed magnitude and temperature dependence of the KIE for the soybean lipoxygenase reaction without fitting any parameters directly to the experimental kinetic data. The temperature dependence of the KIE is determined predominantly by the proton donor-acceptor frequency and the distance dependence of the vibronic couplings for hydrogen and deuterium. The ratio of the overlaps of the hydrogen and deuterium vibrational wavefunctions strongly impacts the magnitude of the KIE but does not significantly influence its temperature dependence. For this enzyme reaction, the large magnitude of the KIE arises mainly from the dominance of tunneling between the ground vibronic states and the relatively large ratio of the overlaps between the corresponding hydrogen and deuterium vibrational wavefunctions. The weak temperature dependence of the KIE is due in part to the dominance of the local component of the proton donor-acceptor motion.     

The anomalous temperature dependence of enzyme-catatlyzed reactions

Breaks or curvatures in the plot of log v_br vs. 1/T have been described in the literature for a number of reactions. If sufficiently small temperature intervals are used in the investigation of the temperature dependence (STI method) a trend line is obtained instead. It is then possible to distinguish regions of linear behavior and regions of anomalous behavior. The effects of e. g. the pH, the ionic strength, and added salts as well as of the size of the enzyme and substrate molecules on the anomalies were examined. The anomalous temperature dependence of the reaction rate was found to be a particularly sensitive indicator of the reaction event on the enzyme. The course of the reaction is evidently affected even by relatively small changes in conformation, particularly in the region of the active center.     

Mechanism of enzymatic Birch reduction: stereochemical course and exchange reactions of benzoyl-CoA reductase

Dearomatizing benzoyl-coenzyme A reductases (BCR) from facultatively anaerobic bacteria are key enzymes in the anaerobic degradation of aromatic compounds. They catalyze the ATP-dependent reduction of benzoyl-CoA (BCoA) to cyclohexa-1,5-diene-1-carboxyl-CoA (dienoyl-CoA). A Birch reduction mechanism involving alternate electron transfer and protonation steps has been proposed for BCR. In this work we reacted BCoA in H2O and D2O, and d5-BCoA in H2O with BCR and the second enzyme of the pathway, dienoyl-CoA hydratase (DCH). The 1,4 hydration product formed from the dienoyl-CoA, 6-hydroxycyclohex-1-ene-1-carbonyl-CoA, was analyzed by several NMR techniques. The results obtained indicate that BCR stereoselectively forms the trans-dienoyl-CoA product, and DCH stereoselectively catalyzes a trans-1,4 water addition. Moreover, unexpected proton exchanges at C-2 and C-6 were observed. They indicate that a free radical intermediate with an unusual low pKa is formed during BCR catalysis. This finding provides evidence for the proposed Birch reduction mechanism of BCR and is in agreement with the established radical mechanism of homologous alpha-hydroxyacyl-CoA dehydratases.     

Some symmetry-induced isotope effects in the kinetics of recombination reactions

Symmetry-induced isotope effects in recombination and collision-induced dissociation reactions are discussed. Progress on understanding the anomalous isotope effects in ozone is reviewed. Then, calculations are performed for the simpler reaction xNe+yNe+H<-->xNeyNe+H, where x and y label either identical or different isotopes. The atomic masses in the model are chosen so that symmetry is the only difference between the systems. Starting from a single potential energy surface, the properties of the bound, quasibound, and continuum states of the neon dimer are calculated. Then, the vibration rotation infinite order sudden approximation is used to calculate cross sections for all possible inelastic and dissociative processes. A rate constant matrix that exactly satisfies detailed balance is constructed. It allows recombination to occur both via direct three-body collisions and via tunneling into the quasibound states of the energy transfer mechanism. The eigenvalue rate coefficients are determined. Significant isotope effects are clearly found, and their behavior depends on the pressure, temperature, and mechanism of the reaction. Both spin statistics and symmetry breaking produce isotope effects. Under most conditions the breaking of symmetry enhances the rates, but a wide spectrum of effects is observed; they range from isotope effects with a normal mass dependence to huge, mass-independent isotope effects to cancellation and even to reversal of the isotope effects. This is the first calculation of symmetry-induced isotope effects in recombination rates from first principles. The relevance of the present effects to ozone recombination is discussed.     

On the temperature dependence of the rate in barrierless reactions in solution

We present a theoretical study of the temperature dependence of the relaxation rate in barrierless reactions. It is shown that a small change in the relative positions of the ground and excited potential surfaces can lead to a significant change in the relaxation rate. We suggest that a small shift of the potential surfaces, during the isoviscous temperature variation, may be responsible for the solvent-induced crossover from an apparent negative activation energy to a large positive activation energy observed in recent experiments.     

Isotope effects on the enzymatic and nonenzymatic reactions of chorismate

The important biosynthetic intermediate chorismate reacts thermally by two competitive pathways, one leading to 4-hydroxybenzoate via elimination of the enolpyruvyl side chain, and the other to prephenate by a facile Claisen rearrangement. Measurements with isotopically labeled chorismate derivatives indicate that both are concerted sigmatropic processes, controlled by the orientation of the enolpyruvyl group. In the elimination reaction of [4-2H]chorismate, roughly 60% of the label was found in pyruvate after 3 h at 60 degrees C. Moreover, a 1.846 +/- 0.057 2H isotope effect for the transferred hydrogen atom and a 1.0374 +/- 0.0005 18O isotope effect for the ether oxygen show that the transition state for this process is highly asymmetric, with hydrogen atom transfer from C4 to C9 significantly less advanced than C-O bond cleavage. In the competing Claisen rearrangement, a very large 18O isotope effect at the bond-breaking position (1.0482 +/- 0.0005) and a smaller 13C isotope effect at the bond-making position (1.0118 +/- 0.0004) were determined. Isotope effects of similar magnitude characterized the transformations catalyzed by evolutionarily unrelated chorismate mutases from Escherichia coli and Bacillus subtilis. The enzymatic reactions, like their solution counterpart, are thus concerted [3,3]-sigmatropic processes in which C-C bond formation lags behind C-O bond cleavage. However, as substantially larger 18O and smaller 13C isotope effects were observed for a mutant enzyme in which chemistry is fully rate determining, the ionic active site may favor a somewhat more polarized transition state than that seen in solution.     

Analysis of kinetic isotope effects for nonadiabatic reactions

Factors influencing the rates of quantum mechanical particle transfer reactions in many-body systems are discussed. The investigations are carried out on a simple model for a proton transfer reaction that captures generic features seen in more realistic models of condensed phase systems. The model involves a bistable quantum oscillator coupled to a one-dimensional double-well reaction coordinate, which is in turn coupled to a bath of harmonic oscillators. Reactive-flux correlation functions that involve quantum-classical Liouville dynamics for chemical species operators and quantum equilibrium sampling are used to estimate the reaction rates. Approximate analytical expressions for the quantum equilibrium structure are derived. Reaction rates are shown to be influenced significantly by both the quantum equilibrium structure and nonadiabatic dynamics. Nonadiabatic dynamical effects are found to play the major role in determining the magnitude of the kinetic isotope effect for the model transfer reaction.     

Technical aspects of quantum chemical modeling of enzymatic reactions: the case of phosphotriesterase

Quantum chemical methods are today a powerful tool in the study of enzymatic reaction mechanisms. In this paper we evaluate the adequacy of some of the technical approximations frequently used in the modeling of enzyme reactions with high level methods. These include the choice of basis set for geometry optimizations and energy evaluation, the choice of dielectric constant to model the enzyme surrounding, and the effects of locking the centers of truncation. As a test case, we choose the phosphotriesterase enzyme, which is a binuclear zinc enzyme that catalyzes the hydrolysis of organophosphate triesters.     

Ultrafast deactivation processes in aminopyridine clusters: excitation energy dependence and isotope effects

Fast excited-state relaxation in H-bonded aminopyridine clusters occurs via hydrogen transfer in the excited state. We used femtosecond pump-probe spectroscopy to characterize the excited-state reaction coordinate. Considerable isotope effects for partially deuterated clusters indicate that H-transfer is the rate-limiting step and validate ab initio calculations in the literature. A nonmonotonous dependence on the excitation energy, however, disagrees with the picture of a simple barrier along the reaction coordinate. An aminopyridine dimer serves as a model for Watson-Crick base pairs, where similar reactions have been predicted by theory.     

Intracrystalline oxygen isotope effects in CuSO4·5H2O and their dependence on crystallization temperature

In copper sulfate pentahydrate the water molecules occupy three different sites, connected with different oxygen isotope ratios. Results of measurements of the change of these isotope ratios with crystallisation temperature are reported. The temperature dependence found here provides the basis for the determination of crystallisation temperatures of hydrated crystals from such intracrystalline oxygen isotope fractionation. Suppositions necessary for the application of this method are discussed.     

Deuterium kinetic isotope effects in microsolvated gas-phase E2 reactions

This work describes the first experimental studies of deuterium kinetic isotope effects (KIEs) for the gas-phase E2 reactions of microsolvated systems. The reactions of F(-)(H(2)O)(n) and OH(-)(H(2)O)(n), where n = 0, 1, with (CH(3))(3)CX (X = Cl, Br), as well as the deuterated analogs of the ionic and neutral reactants, were studied utilizing the flowing afterglow-selected ion flow tube technique. The E2 reactivity is found to decrease with solvation. Small, normal kinetic isotope effects are observed for the deuteration of the alkyl halide, while moderately inverse kinetic isotope effects are observed for the deuteration of the solvent. Minimal clustering of the product ions is observed, but there are intriguing differences in the nature and extent of the clustering process. Electronic structure calculations of the transition states provide qualitative insight into these microsolvated E2 reactions.     

Effect of mutation on enzyme motion in dihydrofolate reductase

Hybrid quantum-classical molecular dynamics simulations of a mutant Escherichia coli dihydrofolate reductase enzyme are presented. Although residue 121 is on the exterior of the enzyme, experimental studies have shown that the mutation of Gly-121 to valine reduces the rate of hydride transfer by a factor of 163. The simulations indicate that the decrease in the hydride transfer rate for the G121V mutant is due to an increase in the free energy barrier. The calculated free energy barrier is higher for the mutant than for the wild-type enzyme by an amount that is consistent with the experimentally observed rate reduction. The calculated transmission coefficients are comparable for the wild-type and mutant enzymes. The simulations suggest that this mutation may interrupt a network of coupled promoting motions proposed to play an important role in DHFR catalysis. This phenomenon has broad implications for protein engineering and drug design.     

Solvent isotope effects upon the kinetics of some simple electrode reactions

The effects of replacing H₂O with D₂O solvent upon the electrochemical kinetics of simple transition-metal redox couples containing aquo, ammine or ethylenediamine ligands have been investigated at mercury electrodes as a means of exploring the possible contribution of ligand-aqueous solvent interactions to the activation barrier to outer-sphere electron transfer. The general interpretation of solvent isotope effects upon electrode kinetics is discussed; it is concluded that double-layer corrected isotopic rate ratios (k^H/k^D)^E determined at a constant electrode potential vs. an aqueous reference electrode, as well as those determined at the respective standard potentials in H₂O and D₂O (k_S^H/k_S^D), have particular significance since the solvent liquid-junction potential can be arranged to be essentially zero. For aquo redox couples, values of (k_S^H/k_S^D) were observed that are substantially greater than unity and appear to be at least partly due to a greater solvent-reorganization barrier in D₂O arising from ligand-solvent hydrogen bonding. For ammine and ethylenediamine complexes values of (k^H/k^D)^E substantially greater than, and smaller than, unity were observed upon the separate deuteration of the ligands and the surrounding solvent respectively. Comparison of isotope rate ratios for corresponding electrochemical and homogeneous outer-sphere reactions involving cationic ammine and aquo complexes yields values of (k^H/k^D) for the former processes that are typically markedly larger than those predicted by the Marcus model from the homogeneous rate ratios. These discrepancies appear to arise from differences in the solvent environments in the transition states for electrochemical and homogeneous reactions.     

Quantum chemical modeling of enzymatic reactions: The case of histone lysine methyltransferase

Quantum chemical cluster models of enzyme active sites are today an important and powerful tool in the study of various aspects of enzymatic reactivity. This methodology has been applied to a wide spectrum of reactions and many important mechanistic problems have been solved. Herein, we report a systematic study of the reaction mechanism of the histone lysine methyltransferase (HKMT) SET7/9 enzyme, which catalyzes the methylation of the N‐terminal histone tail of the chromatin structure. In this study, HKMT SET7/9 serves as a representative case to examine the modeling approach for the important class of methyl transfer enzymes. Active site models of different sizes are used to evaluate the methodology. In particular, the dependence of the calculated energies on the model size, the influence of the dielectric medium, and the particular choice of the dielectric constant are discussed. In addition, we examine the validity of some technical aspects, such as geometry optimization in solvent or with a large basis set, and the use of different density functional methods. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2010