Nonperfect synchronization of reaction center rehybridization in the transition state of the hydride transfer catalyzed by dihydrofolate reductase

It has been suggested that the magnitudes of secondary kinetic isotope effects (2 degrees KIEs) of enzyme-catalyzed reactions are an indicator of the extent of reaction-center rehybridization at the transition state. A 2 degrees KIE value close to the corresponding secondary equilibrium isotope effects (2 degrees EIE) is conventionally interpreted as indicating a late transition state that resembles the final product. The reliability of using this criterion to infer the structure of the transition state is examined by carrying out a theoretical investigation of the hybridization states of the hydride donor and acceptor in the Escherichia coli dihydrofolate reductase (ecDHFR)-catalyzed reaction for which a 2 degrees KIE close to the 2 degrees EIE was reported. Our results show that the donor carbon at the hydride transfer transition state resembles the reactant state more than the product state, whereas the acceptor carbon is more productlike, which is a symptom of transition state imbalance. The conclusion that the isotopically substituted carbon is reactant-like disagrees with the conclusion that would have been derived from the criterion of 2 degrees KIEs and 2 degrees EIEs, but the breakdown of the correlation with the equilibrium isotope effect can be explained by considering the effect of tunneling.     

Am1 and pm3 transition structure for the hydride transfer. A model of reaction catalyzed by dihydrofolate reductase

The transition state structure for the hydride transfer in dihydrofolate reductase, DHFR, enzyme has been calculated with analytical gradients at semiempirical levels: AM1 and PM3. The geometry, electronic structure and transition vector components are qualitatively semiempirical level independent. Comparing the transition structures for the hydride transfer step in models of liver alcohol dehydrogenase, formate dehydrogenase, lactate dehydrogenase, and glutathione reductase, the geometries of these stationary points are transferable and invariant. The topology of the transition structures in these enzymes resembles the one calculated in this paper.     

Effects of the donor-acceptor distance and dynamics on hydride tunneling in the dihydrofolate reductase catalyzed reaction

A significant contemporary question in enzymology involves the role of protein dynamics and hydrogen tunneling in enhancing enzyme catalyzed reactions. Here, we report a correlation between the donor-acceptor distance (DAD) distribution and intrinsic kinetic isotope effects (KIEs) for the dihydrofolate reductase (DHFR) catalyzed reaction. This study compares the nature of the hydride-transfer step for a series of active-site mutants, where the size of a side chain that modulates the DAD (I14 in E. coli DHFR) is systematically reduced (I14V, I14A, and I14G). The contributions of the DAD and its dynamics to the hydride-transfer step were examined by the temperature dependence of intrinsic KIEs, hydride-transfer rates, activation parameters, and classical molecular dynamics (MD) simulations. Results are interpreted within the framework of the Marcus-like model where the increase in the temperature dependence of KIEs arises as a direct consequence of the deviation of the DAD from its distribution in the wild type enzyme. Classical MD simulations suggest new populations with larger average DADs, as well as broader distributions, and a reduction in the population of the reactive conformers correlated with the decrease in the size of the hydrophobic residue. The more flexible active site in the mutants required more substantial thermally activated motions for effective H-tunneling, consistent with the hypothesis that the role of the hydrophobic side chain of I14 is to restrict the distribution and dynamics of the DAD and thus assist the hydride-transfer. These studies establish relationships between the distribution of DADs, the hydride-transfer rates, and the DAD's rearrangement toward tunneling-ready states. This structure-function correlation shall assist in the interpretation of the temperature dependence of KIEs caused by mutants far from the active site in this and other enzymes, and may apply generally to C-H→C transfer reactions.     

Measurement of the alpha-secondary kinetic isotope effect for the reaction catalyzed by mammalian protein farnesyltransferase

Protein farnesytransferase (FTase) catalyzes the transfer of a 15-carbon prenyl group from farnesyl diphosphate (FPP) to the cysteine residue of target proteins and is a member of the newest class of zinc metalloenzymes that catalyze sulfur alkylation. Common substrates of FTase include oncogenic Ras proteins, and therefore inhibitors are under development for the treatment of various cancers. An increased understanding of the salient features of the chemical transition state of FTase may aid in the design of potent inhibitors and enhance our understanding of the mechanism of this class of zinc enzymes. To investigate the transition state of FTase we have used transient kinetics to measure the alpha-secondary 3H kinetic isotope effect at the sensitive C1 position of FPP. The isotope effect for the FTase single turnover reaction using a peptide substrate that is farnesylated rapidly is near unity, indicating that a conformational change, rather than farnesylation, is the rate-limiting step. To look at the chemical step, the kinetic isotope effect was measured as 1.154 +/- 0.006 for a peptide that is farnesylated slowly, and these data suggest that FTase proceeds via a concerted mechanism with dissociative character.     

Conversion and origin of normal and abnormal temperature dependences of kinetic isotope effect in hydride transfer reactions

The effects of substituents on the temperature dependences of kinetic isotope effect (KIE) for the reactions of the hydride transfer from the substituted 5-methyl-6-phenyl-5,6-dihydrophenanthridine (G-PDH) to thioxanthylium (TX(+)) in acetonitrile were examined, and the results show that the temperature dependences of KIE for the hydride transfer reactions can be converted by adjusting the nature of the substituents in the molecule of the hydride donor. In general, electron-withdrawing groups can make the KIE to have normal temperature dependence, but electron-donating groups can make the KIE to have abnormal temperature dependence. Thermodynamic analysis on the possible pathways of the hydride transfer from G-PDH to TX(+) in acetonitrile suggests that the transfers of the hydride anion in the reactions are all carried out by the concerted one-step mechanism whether the substituent is an electron-withdrawing group or an electron-donating group. But the examination of Hammett-type free energy analysis on the hydride transfer reactions supports that the concerted one-step hydride transfer is not due to an elementary chemical reaction. The experimental values of KIE at different temperatures for the hydride transfer reactions were modeled by using a kinetic equation formed according to a multistage mechanism of the hydride transfer including a returnable charge-transfer complex as the reaction intermediate; the real mechanism of the hydride transfer and the root that why the temperature dependences of KIE can be converted as the nature of the substituents are changed were discovered.     

Origin of the temperature dependence of isotope effects in enzymatic reactions: the case of dihydrofolate reductase

The origin of the temperature dependence of kinetic isotope effects (KIEs) in enzyme reactions is a problem of general interest and a major challenge for computational chemistry. The present work simulates the nuclear quantum mechanical (NQM) effects and the corresponding KIE in dihydrofolate reductase (DHFR) and two of its mutants by using the empirical valence bond (EVB) and the quantum classical path (QCP) centroid path integral approach. Our simulations reproduce the overall observed trend while using a fully microscopic rather than a phenomenological picture and provide an interesting insight. It appears that the KIE increases when the distance between the donor and acceptor increases, in a somewhat counter intuitive way. The temperature dependence of the KIE appears to reflect mainly the temperature dependence of the distance between the donor and acceptor. This trend is also obtained from a simplified vibronic treatment, but as demonstrated here, the vibronic treatment is not valid at short and medium distances, where it is essential to use the path integral or other approaches capable of moving seamlessly from the adiabatic to the diabatic limits. It is pointed out that although the NQM effects do not contribute to catalysis in DHFR, the observed temperature dependence can be used to refine the potential of mean force for the donor and acceptor distance and its change due to distanced mutations.     

Simulations of the large kinetic isotope effect and the temperature dependence of the hydrogen atom transfer in lipoxygenase

Elucidating the role of nuclear quantum mechanical (NQM) effects in enzyme catalysis is a topic of significant current interest. Despite the great experimental progress in this field it is important to have theoretical approaches capable of evaluating and analyzing nuclear quantum mechanical contributions to catalysis. In this study, we use the catalytic reaction of lipoxygenase, which is characterized by an extremely large kinetic isotope effect, as a challenging test case for our simulation approach. This is done by applying the quantum classical path (QCP) method with an empirical valence bond potential energy surface. Our computational strategy evaluates the relevant NQM corrections and reproduces the large observed kinetic isotope effect and the temperature dependence of the H atom transfer reaction while being less successful with the D atom transfer reaction. However, the main point of our study is not so much to explore the temperature dependence of the isotope effect but rather to develop and validate an approach for calculations of nuclear quantum mechanical contributions to activation free energies. Here, we find that the deviation between the calculated and observed activation free energies is small for both H and D at all investigated temperatures. The present study also explores the nature of the reorganization energy in the enzyme and solution reactions. It is found that the outer-sphere reorganization energy is extremely small. This reflects the fact that the considered reaction involves a very small charge transfer. The implication of this finding is discussed in the framework of the qualitative vibronic model. The main point of the present study is, however, that the rigorous QCP approach provides a reliable computational tool for evaluating NQM contributions to catalysis even when the given reaction includes large tunneling contributions. Interestingly, our results indicate that the NQM effects in the lipoxygenase reaction are similar in the enzyme and in the reference solution reactions, and thus do not contribute to catalysis. We also reached similar conclusions in studies of other enzymes.     

Tunneling and coupled motion in the Escherichia coli dihydrofolate reductase catalysis

H-transfer was studied in the complex kinetic cascade of dihydrofolate reductase. Intrinsic kinetic isotope effects, their temperature dependence, and other temperature-dependent parameters indicated H-tunneling, but no 1 degrees to 2 degrees coupled motion. The data also suggested environmentally coupled tunneling and commitment to catalysis on pre-steady-state isotope effects.     

Hydride transfer reaction catalyzed by hyperthermophilic dihydrofolate reductase is dominated by quantum mechanical tunneling and is promoted by both inter- and intramonomeric correlated motions

Simulations of hydride and deuteride transfer catalyzed by dihydrofolate reductase from the hyperthermophile Thermotoga maritima (TmDHFR) are presented. TmDHFR was modeled with its active homodimeric quaternary structure, where each monomer has three subdomains. The potential energy function was a combined quantum mechanical and molecular mechanical potential (69 atoms were treated quantum mechanically, and 35 287, by molecular mechanics). The calculations of the rate constants by ensemble-averaged variational transition state theory with multidimensional tunneling predicted that hydride and deuteride transfer at 278 K proceeded with 81 and 80% by tunneling. These percentages decreased to 50 and 49% at 338 K. The kinetic isotope effect was dominated by contributions of bound vibrations and decreased from 3.0 to 2.2 over the temperature range. The calculated rates for hydride and deuteride transfer catalyzed by the hypothetical monomer were smaller by approximately 2 orders of magnitude. At 298 K tunneling contributed 73 and 66% to hydride and deuteride transfer in the monomer. The decreased catalytic efficiency of the monomer was therefore not the result of a decrease of the tunneling contribution but an increase in the quasi-classical activation free energy. The catalytic effect was associated in the dimer with correlated motions between domains as well as within and between subunits. The intrasubunit correlated motions were decreased in the monomer when compared to both native dimeric TmDHFR and monomeric E. coli enzyme. TmDHFR and its E. coli homologue involve similar patterns of correlated interactions that affect the free energy barrier of hydride transfer despite only 27% sequence identity and different quaternary structures.     

Light activates reduction of methotrexate by NADPH in the ternary complex with Escherichia coli dihydrofolate reductase

Methotrexate (MTX), a strong inhibitor of dihydrofolate reductase (DHFR), has been widely used for chemotherapy for many types of cancer as well as for juvenile rheumatoid arthritis. It mimics folate substrates and binds tightly to the active site of DHFR, perhaps in a conformation close to the transition state of the folate catalyzed reaction. Absorption, fluorescence and ultrasensitive Raman difference spectroscopies show that light-activated MTX reacts with NADPH in the enzyme active site, producing 5,8-dihydromethotrexate (5,8-dihydro-MTX) and NADP+. The reaction, which proceeds with a hydride transfer between C4 (pro-R side) of the nicotinamide ring and N5 of the pteridine ring, is similar to that between folate and NADPH except that the hydride is transferred to C6 in this case. Hence, MTX is catalytically competent in its excited state. Most experiments were performed on the Escherichia coli enzyme, but preliminary studies show that the reaction also occurs with human DHFR.     

Comparison of coupled motions in Escherichia coli and Bacillus subtilis dihydrofolate reductase

Hybrid quantum/classical molecular dynamics simulations are used to compare the role of protein motion in the hydride transfer reaction catalyzed by Escherichia coli and Bacillus subtilis dihydrofolate reductase (DHFR). These two enzymes have 44% sequence identity, and the experimentally determined structures and hydride transfer rates are similar. The simulations indicate that the tertiary structures of both enzymes evolve in a similar manner during the hydride transfer reaction. In both enzymes, the donor-acceptor distance decreases to approximately 2.7 Angstroms at the transition state configurations to enable hydride transfer. Zero point energy and hydrogen tunneling effects are found to be significant for both enzymes. Covariance and rank correlation analyses of motions throughout the protein and ligands illustrate that E. coli and B. subtilis DHFR exhibit both similarities and differences in the equilibrium fluctuations and the conformational changes along the collective reaction coordinate for hydride transfer. A common set of residues that play a significant role in the network of coupled motions leading to configurations conducive to hydride transfer for both E. coli and B. subtilis DHFR was identified. These results suggest a balance between conservation and flexibility in the thermal motions and conformational changes during hydride transfer. Homologous protein structures, in conjunction with conformational sampling, enable enzymes with different sequences to catalyze the same hydride transfer reaction with similar efficiency.     

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.     

Simulations of remote mutants of dihydrofolate reductase reveal the nature of a network of residues coupled to hydride transfer

Recent experimental and theoretical studies have proposed that enzymes involve networks of coupled residues throughout the protein that participate in motions accompanying chemical barrier crossing. Here, we have examined portions of a proposed network in dihydrofolate reductase (DHFR) using quantum mechanics/molecular mechanics simulations. The simulations use a hybrid quantum mechanics‐molecular mechanics approach with a recently developed semiempirical AM1‐SRP Hamiltonian that provides accurate results for this reaction. The simulations reproduce experimentally determined catalytic rates for the wild type and distant mutants of E. coli DHFR, underscoring the accuracy of the simulation protocol. Additionally, the simulations provide detailed insight into how residues remote from the active site affect the catalyzed chemistry, through changes in the thermally averaged properties along the reaction coordinate. The mutations do not greatly affect the structure of the transition state near the bond activation, but we observe differences somewhat removed from the point of C? H cleavage that affect the rate. The mutations have global effects on the thermally averaged structure that propagate throughout the enzyme and the current simulations highlight several interactions that appear to be particularly important. © 2014 Wiley Periodicals, Inc.     

The temperature dependence of the kinetic isotope effect in the reaction of F atoms with CH4 and CD4

The kinetic isotope effect k_F+CH4/k_F+CD4 has been determined by reacting F atoms with mixtures of CH₄ and CD₄, using a discharge-flow-mass spectrometric technique. Experiments were carried out at four temperatures in the temperature range 183–298 K. The Arrhenius expression corresponding to the results is k_F+CH4/k_F+CD4=(0.99±0.02)×exp[(100±5)/T]. The present results are compared with previous published experimental and theoretical results.     

Negative kinetic temperature effect on the hydride transfer from NADH analogue BNAH to the radical cation of N-benzylphenothiazine in acetonitrile

The reaction rates of 1-(p-substituted benzyl)-1,4-dihydronicotinamide (G-BNAH) with N-benzylphenothiazine radical cation (PTZ(*+)) in acetonitrile were determined. The results show that the reaction rates (k(obs)) decreased from 2.80 x 10(7) to 2.16 x 10(7) M(-1) s(-1) for G = H as the reaction temperature increased from 298 to 318 K. The activation enthalpies of the reactions were estimated according to Eyring equation to give negative values (-3.4 to -2.9 kcal/mol). Investigation of the reaction intermediate shows that the charge-transfer complex (CT-complex) between G-BNAH and PTZ(*+) was formed in front of the hydride transfer from G-BNAH to PTZ(*+). The formation enthalpy of the CT-complex was estimated by using the Benesi-Hildebrand equation to give the values from -6.4 to -6.0 kcal/mol when the substituent G in G-BNAH changes from CH(3)O to Br. Detailed thermodynamic analyses on each elementary step in the possible reaction pathways suggest that the hydride transfer from G-BNAH to PTZ(*+) occurs by a concerted hydride transfer via a CT-complex. The effective charge distribution on the pyridine ring in G-BNAH at the various stages-the reactant G-BNAH, the charge-transfer complex, the transition-state, and the product G-BNA(+)-was estimated by using the method of Hammett-type linear free energy analysis, and the results show that the pyridine ring carries relative effective positive charges of 0.35 in the CT-complex and 0.45 in the transition state, respectively, which indicates that the concerted hydride transfer from G-BNAH to PTZ(*+) was practically performed by the initial charge (-0.35) transfer from G-BNAH to PTZ(*+) and then followed by the transfer of hydrogen atom with partial negative charge (-0.65). It is evident that the present work would be helpful in understanding the nature of the negative temperature effect, especially on the reaction of NADH coenzyme with the drug phenothiazine in vivo.     

The effect of solvent on a Lewis acid catalyzed Diels-Alder reaction,using computed and experimental kinetic isotope effects

A new transition structure for the Diels-Alder reaction between isoprene and acrolein catalyzed by Et(2)AlCl is found to reconcile reported discrepancies between computed and observed secondary kinetic isotope effects (KIEs). Including the effect of solvent realigns the computed results with experiment demonstrating the importance of nonbond interactions at transition structures. Comparison of experimental and newly predicted KIE data reaffirms the ability of theory and experiment to probe the mechanism and transition structure geometry of organic reactions.     

Biosynthesis of isoprenoids in Escherichia coli: stereochemistry of the reaction catalyzed by farnesyl diphosphate synthase

[formula: see text] Farnesyl diphosphate (FPP) synthase from Escherichia coli catalyzes the condensation of isopentenyl diphosphate (IPP) and geranyl diphosphate (GPP) with selective removal of the pro-R hydrogen at C2 of IPP, the same stereochemistry observed for the pig liver, yeast, and avian enzymes.