Hydride Transfer Step Research Articles

Coupling of protein motions and hydrogen transfer during catalysis byEscherichia colidihydrofolate reductase

The enzyme DHFR (dihydrofolate reductase) catalyses hydride transfer from NADPH to, and protonation of, dihydrofolate. The physical basis of the hydride transfer step catalysed by DHFR from Escherichia coli has been studied through the measurement of the temperature dependence of the reaction rates and the kinetic isotope effects. Single turnover experiments at pH 7.0 revealed a strong dependence of the reaction rates on temperature. The observed relatively large difference in the activation energies for hydrogen and deuterium transfer led to a temperature dependence of the primary kinetic isotope effects from 3.0+/-0.2 at 5 degrees C to 2.2+/-0.2 at 40 degrees C and an inverse ratio of the pre-exponential factors of 0.108+/-0.04. These results are consistent with theoretical models for hydrogen transfer that include contributions from quantum mechanical tunnelling coupled with protein motions that actively modulate the tunnelling distance. Previous work had suggested a coupling of a remote residue,Gly121, with the kinetic events at the active site. However, pre-steady-state experiments at pH 7.0 with the mutant G121V-DHFR, in which Gly121 was replaced with valine, revealed that the chemical mechanism of DHFR catalysis was robust to this replacement. The reduced catalytic efficiency of G121V-DHFR was mainly a consequence of the significantly reduced pre-exponential factors, indicating the requirement for significant molecular reorganization during G121V-DHFR catalysis. In contrast, steady-state measurements at pH 9.5, where hydride transfer is rate limiting, revealed temperature-independent kinetic isotope effects between 15 and 35 degrees C and a ratio of the pre-exponential factors above the semi-classical limit, suggesting a rigid active site configuration from which hydrogen tunnelling occurs. The mechanism by which hydrogen tunnelling in DHFR is coupled with the environment appears therefore to be sensitive to pH.

Kinetic and Chemical Mechanisms of the fabG-Encoded Streptococcus pneumoniae β-Ketoacyl-ACP Reductase

Beta-ketoacyl-acyl carrier protein reductase (KACPR) catalyzes the NADPH-dependent reduction of beta-ketoacyl-acyl carrier protein (AcAc-ACP) to generate (3S)-beta-hydroxyacyl-ACP during the chain-elongation reaction of bacterial fatty acid biosynthesis. We report the evaluation of the kinetic and chemical mechanisms of KACPR using acetoacetyl-CoA (AcAc-CoA) as a substrate. Initial velocity, product inhibition, and deuterium kinetic isotope effect studies were consistent with a random bi-bi rapid-equilibrium kinetic mechanism of KACPR with formation of an enzyme-NADP(+)-AcAc-CoA dead-end complex. Plots of log V/K(NADPH) and log V/K(AcAc)(-)(CoA) indicated the presence of a single basic group (pK = 5.0-5.8) and a single acidic group (pK = 8.0-8.8) involved in catalysis, while the plot of log V vs pH indicated that at high pH an unprotonated form of the ternary enzyme complex was able to undergo catalysis. Significant and identical primary deuterium kinetic isotope effects were observed for V (2.6 +/- 0.4), V/K(NADPH) (2.6 +/- 0.1), and V/K(AcAc)(-)(CoA) (2.6 +/- 0.1) at pH 7.6, but all three values attenuated to values of near unity (1.1 +/- 0.03 or 0.91 +/- 0.02) at pH 10. Similarly, the large alpha-secondary deuterium kinetic isotope effect of 1.15 +/- 0.02 observed for [4R-(2)H]NADPH on V/K(AcAc)(-)(CoA) at pH 7.6 was reduced to a value of unity (1.00 +/- 0.04) at high pH. The complete analysis of the pH profiles and the solvent, primary, secondary, and multiple deuterium isotope effects were most consistent with a chemical mechanism of KACPR that is stepwise, wherein the hydride-transfer step is followed by protonation of the enolate intermediate. Estimations of the intrinsic primary and secondary deuterium isotope effects ((D)k = 2.7, (alpha)(-D)k = 1.16) and the correspondingly negligible commitment factors suggest a nearly full expression of the intrinsic isotope effects on (D)V/K and (alpha)(-D)V/K, and are consistent with a late transition state for the hydride transfer step. Conversely, the estimated intrinsic solvent effect ((D)2(O)k) of 5.3 was poorly expressed in the experimentally derived parameters (D)2(O)V/K and (D)2(O)V (both = 1.2 +/- 0.1), in agreement with the estimation that the catalytic commitment factor for proton transfer to the enolate intermediate is large. Such detailed knowledge of the chemical mechanism of KAPCR may now help guide the rational design of, or inform screening assay-design strategies for, potent inhibitors of this and related enzymes of the short chain dehydrogenase enzyme class.

The Ferredoxin-NADP(H) Reductase fromRhodobacter capsulatus: Molecular Structure and Catalytic Mechanism,

The photosynthetic bacterium Rhodobacter capsulatus contains a ferredoxin (flavodoxin)-NADP(H) oxidoreductase (FPR) that catalyzes electron transfer between NADP(H) and ferredoxin or flavodoxin. The structure of the enzyme, determined by X-ray crystallography, contains two domains harboring the FAD and NADP(H) binding sites, as is typical of the FPR structural family. The FAD molecule is in a hairpin conformation in which stacking interactions can be established between the dimethylisoalloxazine and adenine moieties. The midpoint redox potentials of the various transitions undergone by R. capsulatus FPR were similar to those reported for their counterparts involved in oxygenic photosynthesis, but its catalytic activity is orders of magnitude lower (1-2 s(-)(1) versus 200-500 s(-)(1)) as is true for most of its prokaryotic homologues. To identify the mechanistic basis for the slow turnover in the bacterial enzymes, we dissected the R. capsulatus FPR reaction into hydride transfer and electron transfer steps, and determined their rates using stopped-flow methods. Hydride exchange between the enzyme and NADP(H) occurred at 30-150 s(-)(1), indicating that this half-reaction does not limit FPR activity. In contrast, electron transfer to flavodoxin proceeds at 2.7 s(-)(1), in the range of steady-state catalysis. Flavodoxin semiquinone was a better electron acceptor for FPR than oxidized flavodoxin under both single turnover and steady-state conditions. The results indicate that one-electron reduction of oxidized flavodoxin limits the enzyme activity in vitro, and support the notion that flavodoxin oscillates between the semiquinone and fully reduced states when FPR operates in vivo.

Is Arg418 the Catalytic Base Required for the Hydrolysis Step of the IMP Dehydrogenase Reaction?

The first committed step of guanine nucleotide biosynthesis is the oxidation of inosine 5'-monophosphate (IMP) to xanthosine 5'-monophosphate (XMP) catalyzed by IMP dehydrogenase. The reaction involves the reduction of NAD(+) with the formation of a covalent enzyme intermediate (E-XMP). Hydrolysis of E-XMP requires the enzyme to adopt a closed conformation and is rate-limiting. Thr321, Arg418, and Tyr419 are candidates for the residue that activates water. The substitution of Thr321 has similar, but small, effects on both the hydride transfer and hydrolysis steps. This result suggests that Thr321 influences the reactivity of Cys319, either through a direct interaction or by stabilizing the structure of the active site loop. The hydrolysis of E-XMP is accelerated by the deprotonation of a residue with a pK(a) of approximately 8. A similar deprotonation stabilizes the closed conformation; this residue has a pK(a) of >or=6 in the closed conformation. The substitution of Tyr419 with Phe does not change the pH dependence of either the hydrolysis of E-XMP or the conformational change, which suggests that Tyr419 is not the residue that activates water. In contrast, the conformational change becomes pH-independent when Arg418 is substituted with Gln. Lys can replace the function of Arg418 in the hydrolysis reaction but does not stabilize the closed conformation. The simplest explanation for these observations is that Arg418 serves as the base that activates water in the IMPDH reaction.

Structure and Hydride Transfer Mechanism of a Moderate Thermophilic Dihydrofolate Reductase from Bacillus stearothermophilus and Comparison to Its Mesophilic and Hyperthermophilic Homologues,

Dihydrofolate reductase (DHFR) from a moderate thermophilic organism, Bacillus stearothermophilus, has been cloned and expressed. Physical characterization of the protein (BsDHFR) indicates that it is a monomeric protein with a molecular mass of 18,694.6 Da (0.8), coincident with the mass of 18 694.67 Da calculated from the primary sequence. Determination of the X-ray structure of BsDHFR provides the first structure for a monomeric DHFR from a thermophilic organism, indicating a high degree of conservation of structure in relation to all chromosomal DHFRs. Structurally based sequence alignment of DHFRs indicates the following levels of sequence identity and similarity for BsDHFR: 38 and 58% with Escherichia coli, 35 and 56% with Lactobacillus casei, and 23 and 40% with Thermotoga maritima, respectively. Steady state kinetic isotope effect studies indicate an ordered kinetic mechanism at elevated temperatures, with NADPH binding first to the enzyme. This converts to a more random mechanism at reduced temperatures, reflected in a greatly reduced K(m) for dihydrofolate at 20 degrees C in relation to that at 60 degrees C. A reduction in either temperature or pH reduces the degree to which the hydride transfer step is rate-determining for the second-order reaction of DHF with the enzyme-NADPH binary complex. Transient state kinetics have been used to study the temperature dependence of the isotope effect on hydride transfer at pH 9 between 10 and 50 degrees C. The data support rate-limiting hydride transfer with a moderate enthalpy of activation (E(a) = 5.5 kcal/mol) and a somewhat greater temperature dependence for the kinetic isotope effect than predicted from classical behavior [A(H)/A(D) = 0.57 (0.15)]. Comparison of kinetic parameters for BsDHFR to published data for DHFR from E. coli and T. maritima shows a decreasing trend in efficiency of hydride transfer with increasing thermophilicity of the protein. These results are discussed in the context of the capacity of each enzyme to optimize H-tunneling from donor (NADPH) to acceptor (DHF) substrates.

Mechanistic investigation of a highly active phosphite dehydrogenase mutant and its application for NADPH regeneration

NAD(P)H regeneration is important for biocatalytic reactions that require these costly cofactors. A mutant phosphite dehydrogenase (PTDH-E175A/A176R) that utilizes both NAD and NADP efficiently is a very promising system for NAD(P)H regeneration. In this work, both the kinetic mechanism and practical application of PTDH-E175A/A176R were investigated for better understanding of the enzyme and to provide a basis for future optimization. Kinetic isotope effect studies with PTDH-E175A/A176R showed that the hydride transfer step is (partially) rate determining with both NAD and NADP giving (D)V values of 2.2 and 1.7, respectively, and (D)V/K(m,phosphite) values of 1.9 and 1.7, respectively. To better comprehend the relaxed cofactor specificity, the cofactor dissociation constants were determined utilizing tryptophan intrinsic fluorescence quenching. The dissociation constants of NAD and NADP with PTDH-E175A/A176R were 53 and 1.9 microm, respectively, while those of the products NADH and NADPH were 17.4 and 1.22 microm, respectively. Using sulfite as a substrate mimic, the binding order was established, with the cofactor binding first and sulfite binding second. The low dissociation constant for the cofactor product NADPH combined with the reduced values for (D)V and k(cat) implies that product release may become partially rate determining. However, product inhibition does not prevent efficient in situ NADPH regeneration by PTDH-E175A/A176R in a model system in which xylose was converted into xylitol by NADP-dependent xylose reductase. The in situ regeneration proceeded at a rate approximately fourfold faster with PTDH-E175A/A176R than with either WT PTDH or a NADP-specific Pseudomonas sp.101 formate dehydrogenase mutant with a total turnover number for NADPH of 2500.

A Theoretical Comparison of Lewis Acid vs Bronsted Acid Catalysis for n-Hexane → Propane + Propene

Cracking of an all-trans n-alkane, via idealized Lewis acid and Bronsted acid catalysis, was examined using density functional theory. Optimized geometries and transitions states were determined for catalyst-reactant complexes, using AlCl3 and HCl.AlCl3 as the Lewis and Bronsted acids. For the Lewis acid cycle, hydride-transfer steps are seen to have large barriers in both forward and reverse directions, and an unstable physisorbed carbenium ion (lying 20 kcal mol(-1) above the chemisorbed intermediate) is the launching point for the beta-scission that leads to products. For the Bronsted acid cycle, proton-transfer steps have smaller barriers in both forward and reverse directions, and a semistable physisorbed alkanium ion is the launching point for the alkanium alpha-scission that leads to products. In the idealized Lewis cycle, formation of HCl units (and hence Bronsted acids) was found to be a common side reaction. A recent ionic-liquid catalysis study is mentioned as motivation, although our study is not a computational modeling study; we are more interested in the fundamental differences between Brosnted and Lewis mechanisms rather than merely mimicking a particular system. However, results of exploratory optimizations of various intermediates with Al2Cl7- as the catalyst are presented to provide the first step for future modeling studies on the ionic liquid system.

Reaction Profiles and Substituent Effects for the Reduction of Carbonyl Compounds with Monomer and Dimer Simple Metal Hydride Reagents

Reaction profiles and energetics for the reactions of substituted benzaldehydes with a series of different simple metal hydrides (BH(3), BMeH(2), BMe(2)H, AlH(3), and AlMe(2)H) are examined computationally. B3LYP/6-31G optimizations and MP2/6-311G single point energy calculations revealed that the Al reagents are more reactive than B reagents. Replacement of H with Me on BH(3) or AlH(3) makes the reduction transition state (TS) less stable. Accordingly the overall reactivity is in the order AlH(3) > AlMe(2)H > BH(3) > BMe(2)H. The Hammett rho value for substituted benzaldehydes (BAs) is negative for the initial complex formation and positive for the hydride-transfer step. The size and the sign of the apparent rho value depend on the relative stabilities of the separated reactants and the complex. The TS structures vary according to the Hammond postulate for substituted BAs and the variation is reflected in carbonyl-carbon and aldehyde-deuterium isotope effects. Comparison of the reaction profiles of the monomer and dimer reagents reveals that the real reacting species is the monomer in the gas phase for BH(3) but the dimer for BMe(2)H.

A Catalytic Triad Is Responsible for Acid−Base Chemistry in the Ascaris suum NAD−Malic Enzyme

The pH dependence of kinetic parameters of several active site mutants of the Ascaris suum NAD-malic enzyme was investigated to determine the role of amino acid residues likely involved in catalysis on the basis of three-dimensional structures of malic enzyme. Lysine 199 is positioned to act as the general base that accepts a proton from the 2-hydroxyl of malate during the hydride transfer step. The pH dependence of V/K(malate) for the K199R mutant enzyme reveals a pK of 5.3 for an enzymatic group required to be unprotonated for activity and a second pK of 6.3 that leads to a 10-fold loss in activity above the pK of 6.3 to a new constant value up to pH 10. The V profile for K199R is pH independent from pH 5.5 to pH 10 and decreases below a pK of 4.9. Tyrosine 126 is positioned to act as the general acid that donates a proton to the enolpyruvate intermediate to form pyruvate. The pH dependence of V/K(malate) for the Y126F mutant is qualitatively similar to K199R, with a requirement for a group to be unprotonated for activity with a pK of 5.6 and a partial activity loss of about 3-fold above a pK of 6.7 to a new constant value. The Y126F mutant enzyme is about 60000-fold less active than the wild-type enzyme. In contrast to K199R, the V rate profile for Y126F also shows a partial activity loss above pH 6.6. The wild-type pH profiles were reinvestigated in light of the discovery of the partial activity change for the mutant enzymes. The wild-type V/K(malate) pH-rate profile exhibits the requirement for a group to be unprotonated for catalysis with a pK of 5.6 and also shows the partial activity loss above a pK of 6.4. The wild-type V pH-rate profile decreases below a pK of 5.2 and is pH independent from pH 5.5 to pH 10. Aspartate 294 is within hydrogen-bonding distance to K199 in the open and closed forms of malic enzyme. D294A is about 13000-fold less active than the wild-type enzyme, and the pH-rate profile for V/K(malate) indicates the mutant is only active above pH 9. The data suggest that the pK present at about pH 5.6 in all of the pH profiles represents D294, and during catalysis D294 accepts a proton from K199 to allow K199 to act as a general base in the reaction. The pK for the general acid in the reaction is not observed, consistent with rapid tautomerization of enolpyruvate. No other ionizable group in the active site is likely responsible for the partial activity change observed in the pH profiles, and thus the group responsible is probably remote from the active site and the effect on activity is transmitted through the protein by a conformational change.

Role of Methionine-13 in the Catalytic Mechanism of 6-Phosphogluconate Dehydrogenase from Sheep Liver

The crystal structure of sheep liver 6-phosphogluconate dehydrogenase (6PGDH) shows marked differences in the position of the nicotinamide mononucleotide (NMN) moiety of NADP(+) and NADPH (Adams, J. M., Grant, H. E., Gover, S., Naylor, C. E., and Phillips, C. (1994) Structure 2, 651-668). A methionine side chain (Met13) interacts with the si face of NADP(+) in the complex with the oxidized coenzyme, is likely to affect the binding mode of the nicotinamide ring of NADP(+), and may play a role in catalysis in the 6PGDH reaction. To check this possibility we performed site-directed mutagenesis, changing M13 to a number of residues including V, I, C, F, and Q. Mutant enzymes were characterized with respect to their kinetic parameters and primary deuterium isotope effects. All mutations resulted in a decrease in affinity of the enzyme for NADP(+), but not NADPH. In addition, the M13 to C (M13C), M13F, and M13Q mutant enzymes exhibited a decrease of at least an order of magnitude in V/E(t). The deuterium isotope effects on V and V/K(6PG) were decreased to about 1.2 for the M13F and M13C mutant enzymes, while they were increased to about 2.4 for the M13Q enzyme (a value of 1.8-1.9 is obtained for the wild-type enzyme). In at least three instances changes in the overall rate of the oxidative decarboxylation reaction relative to other steps along the reaction pathway were observed. Isotope effects indicate that the hydride transfer steps can become either more or less rate-determining dependent on the substitution. Data are consistent with a significant role of M13 in the orientation of the cofactor nicotinamide ring in the mechanism of 6PGDH, likely with respect to geometry and distance of the ring from C3 of 6PG.

Role of Rare Earth Cations in Y Zeolite for Hydrocarbon Cracking

Reaction kinetics data were collected for isobutane conversion over a series of ultra stable Y (USY) zeolite catalysts with and without rare earth cations and subjected to various extents of dealumination by steaming. We conducted these reaction studies at low temperatures (523-573 K) using isobutane feed streams containing known levels of isobutylene (100-400 ppm) so that the kinetics were controlled by bimolecular hydride transfer and oligomerization/beta-scission processes with little or no participation of monomolecular initiation reactions. These experimental conditions led to stable catalyst performance with the main products of isobutane conversion being propane, n-butane, and isopentane, with smaller amounts of propylene, trans-2-butene, and cis-2-butene. The rates of formation of these products per Brønsted acid site (as counted by pyridine adsorption) depended exponentially on Brønsted acid site density, regardless of whether the catalyst contained rare earth cations. Kinetic modeling showed an exponential dependence of hydride transfer and oligomerization/ beta-scission reaction rates on Brønsted acid site density which translated into composite activation energies for these reactions having a linear relationship with site density. Based on results in the literature from theoretical calculations, we suggest that increasing Brønsted acid site density in zeolite Y leads to larger zeolite elasticity, increased stabilization of cationic transition states, and lower composite activation barriers for hydride transfer and beta-scission steps. The role of rare earth cations, therefore, is to ensure the retention of high Brønsted acid site density under hydrothermal conditions, such as in fluid catalytic cracking (FCC) regenerators, where steam would dealuminate the Y zeolite framework and reduce this site density. It is for this reason that hydride transfer reaction rates are high in the presence of rare earth cations and lead to higher yields of less olefinic gasoline during FCC.

An appraisal of multiple NADPH binding-site models proposed for cytochrome P450 reductase, NO synthase, and related diflavin reductase systems.

The diflavin reductases exemplified by mammalian cytochrome P450 reductase catalyze NADPH dehydrogenation and electron transfer to an associated monooxygenase. It has recently been proposed that double occupancy of the NADPH dehydrogenation site inhibits the NADPH to FAD hydride transfer step in this series of enzymes. This has important implications for the mechanism of enzyme turnover. However, the conclusions are drawn from a series of pre-steady-state stopped-flow experiments in which the data analysis and interpretation are flawed. Recent data published for P450-BM3 reductase show a decrease in the rate constant for pre-steady-state flavin oxidation with increasing NADP(+) concentration. This is interpreted as evidence of inhibition by multiple substrate binding. A detailed reanalysis shows that the data are in fact consistent with a simple single-binding-site model in which reversible hydride transfer causes the observed effect. Data for the related systems are also discussed.

Vibrationally enhanced hydrogen tunneling in the Escherichia coli thymidylate synthase catalyzed reaction.

The enzyme thymidylate synthase (TS) catalyzes a complex reaction that involves forming and breaking at least six covalent bonds. The physical nature of the hydride transfer step in this complex reaction cascade has been studied by means of isotope effects and their temperature dependence. Competitive kinetic isotope effects (KIEs) on the second-order rate constant (V/K) were measured over a temperature range of 5-45 degrees C. The observed H/T ((T)V/K(H)) and D/T ((T)V/K(D)) KIEs were used to calculate the intrinsic KIEs throughout the temperature range. The Swain-Schaad relationships between the H/T and D/T V/K KIEs revealed that the hydride transfer step is the rate-determining step at the physiological temperature of Escherichia coli (20-30 degrees C) but is only partly rate-determining at elevated and reduced temperatures. H/D KIE on the first-order rate constant k(cat) ((D)k = 3.72) has been previously reported [Spencer et al. (1997) Biochemistry 36, 4212-4222]. Additionally, the Swain-Schaad relationships between that (D)k and the V/K KIEs reported here suggested that at 20 degrees C the hydride transfer step is the rate-determining step for both rate constants. Intrinsic KIEs were calculated here and were found to be virtually temperature independent (DeltaE(a) = 0 within experimental error). The isotope effects on the preexponential Arrhenius factors for the intrinsic KIEs were A(H)/A(T) = 6.8 +/- 2.8 and A(D)/A(T) = 1.9 +/- 0.25. Both effects are significantly above the semiclassical (no-tunneling) predicted values and indicate a contribution of quantum mechanical tunneling to this hydride transfer reaction. Tunneling correction to transition state theory would predict that these isotope effects on activation parameters result from no energy of activation for all isotopes. Yet, initial velocity measurements over the same temperature range indicate cofactor inhibition and result in significant activation energy on k(cat) (4.0 +/- 0.1 kcal/mol). Taken together, the temperature-independent KIEs, the large isotope effects on the preexponential Arrhenius factors, and a significant energy of activation all suggest vibrationally enhanced hydride tunneling in the TS-catalyzed reaction.

Active site residues and mechanism of UDP-glucose dehydrogenase.

UDP-glucose dehydrogenase catalyzes the NAD+-dependent twofold oxidation of UDP-glucose to give UDP-glucuronic acid. A sequestered aldehyde intermediate is produced in the first oxidation step and a covalently bound thioester is produced in the second oxidation step. This work demonstrates that the Streptococcus pyogenes enzyme incorporates a single solvent-derived oxygen atom during catalysis and probably does not generate an imine intermediate. The reaction of UDP-[6",6"-di-2H]-d-glucose is not accompanied by a primary kinetic isotope effect, indicating that hydride transfer is not rate determining in this reaction. Studies with a mutant of the key active site nucleophile, Cys260Ala, show that it is capable of both reducing the aldehyde intermediate, and oxidizing the hydrated form of the aldehyde intermediate but is incapable of oxidizing UDP-glucose to UDP-glucuronic acid. In the latter case, a ternary Cys260Ala/aldehyde intermediate/NADH complex is presumably formed, but it does not proceed to product as both release and hydration of the bound aldehyde occur slowly. A washout experiment demonstrates that the NADH in this ternary complex is not exchangeable with external NADH, indicating that dissociation only occurs after the addition of a nucleophile to the aldehyde carbonyl. Studies on Thr118Ala show that the value of kcat is reduced 160-fold by this mutation, and that the reaction of UDP-D-[6",6"-di-2H]-glucose is now accompanied by a primary kinetic isotope effect. This indicates that the barriers for the hydride transfer steps have been selectively increased and supports a mechanism in which an ordered water molecule (H-bonded to Thr118) serves as the catalytic base in these steps.

Barriers to Hydride Transfer in Wild Type and Mutant Dihydrofolate Reductase from E. coli

We have calculated the distribution of activation energies for the key hydride transfer step in the reduction of dihydrofolate by the enzyme dihydrofolate reductase in wild-type DHFR as well as the G121S and G121V variants of the Escherichia coli enzyme using combined quantum mechanical/molecular mechanical methods. Comparison of the activation energy distributions present in the three protein systems demonstrates that the ensemble of energy barriers differs in each. The features of these energy barrier distributions are consistent with experimentally determined reaction rates for the three proteins and support experimental observations demonstrating time-dependent variation in the reaction rate of individual enzyme molecules. This occurrence suggests that several distinct conformational substates in the proteins modify the potential energy surface for the reaction. A select set of distances and dihedral angles was observed to correlate well with the progress of the reaction or with the presence of protei...

Correlated motion and the effect of distal mutations in dihydrofolate reductase.

Dihydrofolate reductase (DHFR) catalyzes the reduction of dihydrofolate to tetrahydrofolate. The catalytic rate in this system has been found to be significantly affected by mutations far from the site of chemical activity in the enzyme [Rajagopalan, P. T. R, Lutz, S., and Benkovic, S. J. (2002) Biochemistry 41, 12618-12628]. On the basis of extensive computer simulations for wild-type DHFR from Escherichia coli and four mutants (G121S, G121V, M42F, and M42F/G121S), we show that key parameters for catalysis are changed. The parameters we study are relative populations of different conformations sampled and hydrogen bonds. We find that the mutations result in long-range structural perturbations, rationalizing the effects that the mutations have on the kinetics of the enzyme. Such perturbations also provide a rationalization for the reported nonadditivity effect for double mutations. We finally examine the role a structural perturbation will have on the hydride transfer step. On the basis of our new findings, we discuss the role of coupled motions between distant regions in the enzyme, which previously was reported by Radkiewicz and Brooks.

Hydride transfer catalyzed by xylose isomerase: mechanism and quantum effects.

We have applied molecular dynamics umbrella-sampling simulation and ensemble-averaged variational transition state theory with multidimensional tunneling (EA-VTST/MT) to calculate the reaction rate of xylose-to- xylulose isomerization catalyzed by xylose isomerase in the presence of two Mg2+ ions. The calculations include determination of the free energy of activation profile and ensemble averaging in the transmission coefficient. The potential energy function is approximated by a combined QM/MM/SVB method involving PM3 for the quantum mechanical (QM) subsystem, CHARMM22 and TIP3P for the molecular mechanical (MM) environment, and a simple valence bond (SVB) local function of two bond distances for the hydride transfer reaction. The simulation confirms the essential features of a mechanism postulated on the basis of kinetics and X-ray data by Whitlow et al. (Whitlow, M.; Howard, A. J.; Finzel, B. C.; Poulos, T. L.; Winborne, E.; Gilliland, G. L. Proteins 1991, 9, 153) and Ringe, Petsko, and coworkers (Labie, A.; Allen, K.-N.; Petsko, G. A.; Ringe, D. Biochemistry 1994, 33, 5469). This mechanism involves a rate-determining 1,2-hydride shift with prior and post proton transfers. Inclusion of quantum mechanical vibrational energy is important for computing the free energy of activation, and quantum mechanical tunneling effects are essential for computing kinetic isotope effects (KIEs). It is found that 85% of the reaction proceeds by tunneling and 15% by overbarrier events. The computed KIE for the ratio of hydride to deuteride transfer is in good agreement with the experimental results. The molecular dynamics simulations reveal that proton and hydride transfer reactions are assisted by breathing motions of the mobile Mg2+ ion in the active site, providing evidence for concerted motion of Mg2+ during the hydride transfer step.

Hydride-exchange reactions between NADH and NAD+ model compounds under non-steady-state conditions. Apparent and real kinetic isotope effects.

The kinetics of the hydride exchange reaction between NADH model compound 10-methyl-9,10-dihydroacridine (MAH) and 1-benzyl-3-cyanoquinolinium (BQCN+) ion in acetonitrile were studied at temperatures ranging from 291 to 325 K. The extent of reaction-time profiles during the first half-lives are compared with theoretical data for the simple single-step mechanism and a 2-step mechanism involving initial donor/acceptor complex formation followed by unimolecular hydride transfer. The profiles for the reactions of MAH deviate significantly from those expected for the simple single-step mechanism with the deviation increasing with increasing temperature. The deviation from simple mechanism behavior is much less pronounced for the reactions of 10-methyl-9,10-dihydroacridine-10,10-d2 (MAD) which gives rise to extent of reaction dependent apparent kinetic isotope effects (KIEapp). Excellent fits of the experimental extent of reaction-time profiles with theoretical data for the 2-step mechanism, in the pre-steady-state time period, were observed in all cases. Resolution of the kinetics of the hydride exchange reaction into the microscopic rate constants over the entire temperature range resulted in real kinetic isotope effects for the hydride transfer step ranging from 40 (291 K) to 8.2 (325 K). That the reaction involves significant hydride tunnelling was verified by the magnitudes of the Arrhenius parameters; Ea D - EaH = 8.7 kcal mol-1 and AD/AH = 8 x 10(4). An electron donor acceptor complex (lambda max = 526 nm) was observed to be a reaction intermediate. Theoretical extent of reaction-time profile data are discussed for the case where a reaction intermediate is formed in a non-productive side equilibrium as compared to the case where it is a real intermediate on the reaction coordinate between reactants and products. The common assumption that the two cases are kinetically indistinguishable is shown to be incorrect.

Coupling interactions of distal residues enhance dihydrofolate reductase catalysis: mutational effects on hydride transfer rates.

Recently, the participation in catalysis of residues spatially removed from the enzyme's active site has received considerable attention. The influence of the distal Gly-121 residue on the chemical step of hydride transfer in dihydrofolate reductase (DHFR) catalysis had been demonstrated previously [Cameron, C. E., and Benkovic, S. J. (1997) Biochemistry 36, 15792-15800]. In our continuing effort to identify catalytically important residues that are distal from the active site, we used sequence conservation information, kinetic data on site-directed mutants, dynamic motion information from NMR methods, and correlated motions from MD simulations to identify a subset of residues. Among them, the region spanning positions 41-45 is distal to the active site and was chosen as the focus for the mutagenesis and kinetic studies reported here. Specifically, the highly conserved Met-42 was selected for site-directed mutagenesis. While the reaction kinetics for the M42F mutant enzyme did not deviate from wild-type behavior, a 41-fold reduction in the forward hydride transfer rate was found for the M42W mutant. Given the established role of Gly-121 in the hydride transfer process, double mutant enzymes involving positions 42 and 121 were constructed and characterized. These double mutant enzymes generally showed little changes in substrate and cofactor binding but synergistic decreases in forward hydride transfer rates, while the decreases in reverse rates were additive. Along with supporting information from mixed quantum/classical MD simulations [Agarwal, P. K., Billeter, S. R., Rajagopalan, P. T., Benkovic, S. J., and Hammes-Schiffer, S. (2002) Proc. Natl. Acad. Sci. U.S.A. 99, 2794-2799], the data suggest that a coupled dynamic motion of these distal residues enhances crossing of the chemical reaction barrier and imply an expanded nonstatic role for the protein fold in catalysis.

Computational Methods for the Study of Enzymic Reaction Mechanisms. 1. Application to the Hydride Transfer Step in the Catalysis of Dihydrofolate Reductase

We report results and assessment from the use of several “state-of-the-art” computational methods to investigate the effect of the Asp27 ionization state on the hydride-ion transfer step in the enzymic reduction of folate and dihydrofolate (DHF) by DHFR from Escherichia coli (E. coli). The active forms of the DHFR complex are assumed protonated on the pterin/dihydropterin ring of the folate/dihydrofolate molecule prior to transfer of the hydride ion from the nicotinamide adenine dinucleotide phosphate (NADPH) cofactor. The calculations have been carried out for both protonated (neutral) and unprotonated (negatively charged) states of the conserved active-site Asp27 (E. coli) residue in the DHFR complexes. First, geometry optimizations at the semiempirical (PM3), density functional (DFT) and ab initio levels were performed on reaction-analogue clusters to obtain reactant, transition state (TS), and product complexes. The DFT and ab initio reaction energies obtained for the unprotonated Asp complexes were e...