Multiple hydrogen kinetic isotope effects for enzymes catalyzing exchange with solvent: application to alanine racemase

Alanine racemase catalyzes the pyridoxal phosphate-dependent interconversion of the D- and L-isomers of alanine. Previous studies have shown that the enzyme employs a two-base mechanism in which Lys39 and Tyr265 are the acid/base catalysts. It is thus possible that stereoisomerization of the external aldimine intermediates occurs through a concerted double proton transfer without the existence of a distinct carbanionic intermediate. This possibility was tested by the application of multiple kinetic isotope effect (KIE) methodology to alanine racemase. The mutual dependence of primary substrate and solvent deuterium KIEs has been measured using equilibrium perturbation-type experiments. The conceptually straightforward measurement of the substrate KIE in H(2)O is complemented with a less intuitive protium washout perturbation-type measurement in D(2)O. The primary substrate KIE in the D --> L direction at 25 degrees C is reduced from 1.297 in H(2)O to 1.176 in D(2)O, while in the L --> D direction it is reduced from 1.877 in H(2)O to 1.824 in D(2)O. Similar reductions are also observed at 65 degrees C, the temperature to which the Bacillus stearothermophilus enzyme is adapted. These data strongly support a stepwise racemization of stereoisomeric aldimine intermediates in which a substrate-based carbanion is an obligatory intermediate. The ionizations observed in k(cat)/K(M) pH profiles have been definitively assigned based on the DeltaH(ion) values of the observed pK(a)'s with alanine and on the pH dependence of k(cat)/K(M) for the alternative substrate serine. The acidic pK(a) in the bell-shaped curve is due to the phenolic hydroxyl of Tyr265, which must be unprotonated for reaction with either isomer of alanine. The basic pK(a) is due to the substrate amino group, which must be protonated to react with Tyr265-unprotonated enzyme. A detailed reaction mechanism incorporating these results is proposed.     

Role of Y94 in proton and hydride transfers catalyzed by thymidylate synthase

Thymidylate synthase (TS) catalyzes the substitution of a carbon-bound proton in a uracil base by a methyl group to yield thymine in the de novo biosynthesis of this DNA base. The enzymatic mechanism involves making and breaking several covalent bonds. Traditionally, a conserved tyrosine (Y94 in Escherichia coli, Y146 in Lactobacillus casei, and Y135 in humans) was assumed to serve as the general base catalyzing the proton abstraction. That assumption was examined here by comparing the nature of the proton abstraction using wild-type (wt) E. coli TS (ecTS) and its Y94F mutant (with a turnover rate reduced by 2 orders of magnitude). A subsequent hydride transfer was also studied using the wt and Y94F. The physical nature of both H-transfer steps was examined by determining intrinsic kinetic isotope effects (KIEs). Surprisingly, the findings did not suggest a direct role for Y94 in the proton abstraction step. The effect of this mutation on the subsequent hydride transfer was examined by a comparison of the temperature dependency of the intrinsic KIE on both the wt and the mutant. The intrinsic KIEs for Y94F at physiological temperatures were slightly smaller than those for wt but, otherwise, were as temperature-independent, suggesting a perfectly preorganized reaction coordinate for both enzymes. At reduced temperatures, however, the KIE for the mutant increased with a decrease in temperature, indicating a poorly preorganized reaction coordinate. Other kinetic and structural properties were also compared, and the findings suggested that Y94 is part of a H-bond network that plays a critical role at a step between the proton and the hydride transfers, presumably the dissociation of H4folate from the covalently bound intermediate. The possibility that no single residue serves as the general base in question but, rather, that the whole network of H-bonds at the active site catalyzes proton abstraction is discussed.     

Importance of barrier shape in enzyme-catalyzed reactions. Vibrationally assisted hydrogen tunneling in tryptophan tryptophylquinone-dependent amine dehydrogenases

C-H bond breakage by tryptophan tryptophylquinone (TTQ)-dependent methylamine dehydrogenase (MADH) occurs by vibrationally assisted tunneling (Basran, J., Sutcliffe, M. J., and Scrutton, N. S. (1999) Biochemistry 38, 3218--3222). We show here a similar mechanism in TTQ-dependent aromatic amine dehydrogenase (AADH). The rate of TTQ reduction by dopamine in AADH has a large, temperature independent kinetic isotope effect (KIE = 12.9 +/- 0.2), which is highly suggestive of vibrationally assisted tunneling. H-transfer is compromised with benzylamine as substrate and the KIE is deflated (4.8 +/- 0.2). The KIE is temperature-independent, but reaction rates are strongly dependent on temperature. With tryptamine as substrate reaction rates can be determined only at low temperature as C-H bond cleavage is rapid, and an exceptionally large KIE (54.7 +/- 1.0) is observed. Studies with deuterated tryptamine suggest vibrationally assisted tunneling is the mechanism of deuterium and, by inference, hydrogen transfer. Bond cleavage by MADH using a slow substrate (ethanolamine) occurs with an inflated KIE (14.7 +/- 0.2 at 25 degrees C). The KIE is temperature-dependent, consistent with differential tunneling of protium and deuterium. Our observations illustrate the different modes of H-transfer in MADH and AADH with fast and slow substrates and highlight the importance of barrier shape in determining reaction rate.     

The interpretation of multiple-step transient-state kinetic isotope effects

In contrast to steady-state kinetic isotope effects (KIEs), transient-state tKIEs are both time and signal dependent and therefore require a very different form of theory for their interpretation. We have previously provided such a theory for the case of single-step isotopic substitutions. No such properly derived theory applicable to the analysis of multiple-step isotopic substitutions required by transient-state solvent isotope effect studies has been available up to this time. Here, we set forth a more general form of that theory which is applicable to multiple-step substituted cases. We prove three theorems: 1. the observed transient-state KIE for any given reactive component in the reaction sequence evaluated at zero time (tKIE(0)) is in fact the arithmetic product of the intrinsic KIEs of all the steps that precede the formation of that component. 2. The observed tKIE(0) is completely independent of the intrinsic KIEs of any reverse step in the reaction. 3. The intrinsic KIE of any step may be obtained by dividing the value of the tKIE(0) for that step by the value of the tKIE(0) of the immediately preceding step in the reaction sequence.     

Triple isotopic labeling and kinetic isotope effects: exposing H-transfer steps in enzymatic systems

Kinetic isotope effect (KIE) studies can provide insight into the mechanism and kinetics of specific chemical steps in complex catalytic cascades. Recent results from hydrogen KIE measurements have examined correlations between enzyme dynamics and catalytic function, leading to a surge of studies in this area. Unfortunately, most enzymatic H-transfer reactions are not rate limiting, and the observed KIEs do not reliably reflect the intrinsic KIEs on the chemical step of interest. Given their importance to understanding the chemical step under study, accurate determination of the intrinsic KIE from the observed data is essential. In 1975, Northrop developed an elegant method to assess intrinsic KIEs from their observed values [Northrop, D. B. (1975) Steady-state analysis of kinetic isotope effects in enzymic reactions. Biochemistry 14, 2644-2651]. The Northrop method involves KIE measurements using all three hydrogen isotopes, where one of them serves as the reference isotope. This method has been successfully used with different combinations of observed KIEs over the years, but criteria for a rational choice of reference isotope have never before been experimentally determined. Here we compare different reference isotopes (and hence distinct experimental designs) using the reduction of dihydrofolate and dihydrobiopterin by two dissimilar enzymes as model reactions. A number of isotopic labeling patterns have been applied to facilitate the comparative study of reference isotopes. The results demonstrate the versatility of the Northrop method and that such experiments are limited only by synthetic techniques, availability of starting materials, and the experimental error associated with the use of distinct combinations of isotopologues.     

Strength of the Calpha H..O hydrogen bond of amino acid residues

Although the peptide C(alpha)H group has historically not been thought to form hydrogen bonds within proteins, ab initio quantum calculations show it to be a potent proton donor. Its binding energy to a water molecule lies in the range between 1.9 and 2.5 kcal/mol for nonpolar and polar amino acids; the hydrogen bond (H-bond) involving the charged lysine residue is even stronger than a conventional OH..O interaction. The preferred H-bond lengths are quite uniform, about 3.32 A. Formation of each interaction results in a downfield shift of the bridging hydrogen's chemical shift and a blue shift in the C(alpha)H stretching frequency, potential diagnostics of the presence of such an H-bond within a protein.     

The proton inventory technique

The proton inventory technique uses the dependence of enzymic reaction rate on the atom fraction of deuterium present in mixtures of protium oxide and deuterium oxide to deduce for simple cases the number of exchangeable hydrogenic sites that produce isotope effects, and the magnitude of the isotope effect generated at each site. For more complex cases, other information, such as the participation of more than a single step in limiting the rate, may be obtained. The background of the method, the conduct of the experiments and the interpretation of the results are briefly reviewed. The method is then illustrated in its application to various enzyme systems by a series of case histories.     

Use of an altered sugar-nucleotide to unmask the transition state for alpha(2-->6) sialyltransferase

Rat liver alpha(2-->6) sialyltransferase catalyzes the formation of a glycosidic bond between N-acetylneuraminic acid and the 6-hydroxyl group of a galactose residue at the nonreducing terminus of an oligosaccharide. This reaction has been investigated through the use of the novel sugar-nucleotide donor substrate UMP-NeuAc. A series of UMP-NeuAc radioisotopomers were prepared by chemical deamination of the corresponding CMP-NeuAc precursors. Kinetic isotope effects (KIEs) on V/K were measured using mixtures of radiolabeled UMP-NeuAc's as the donor substrate and N-acetyllactosamine as the acceptor. The secondary beta-(2)H KIE was 1.218 +/- 0.010, and the primary (14)C KIE was 1.030 +/- 0.010. A large inverse (3)H binding isotope effect of 0.944 +/- 0.010 was measured at the terminal carbon of the NeuAc glycerol side chain. These KIEs observed using UMP-NeuAc are much larger than those previously measured with CMP-NeuAc [Bruner, M., and Horenstein, B. A. (1998) Biochemistry 37, 289-297]. Solvent deuterium isotope effects of 1.3 and 2.6 on V/K and V(max) were observed with CMP-NeuAc as the donor, and it is revealing that these isotope effects vanished with use of the slow donor substrate UMP-NeuAc. Bell-shaped pH versus rate profiles were observed for V(max) (pK(a) values = 5.5, 9.0) and V/K(UMP)(-)(NeuAc) (pK(a)values = 6.2, 9.0). The results are considered in terms of a mechanism involving an isotopically sensitive conformational change which is independent of the glycosyl transfer step. The isotope effects reveal that the enzyme-bound transition state bears considerable charge on the N-acetylneuraminic acid residue, and this and other features of this mechanism provide new directions for sialyltransferase inhibitor design.     

Origins of deuterium kinetic isotope effects on the proton transfers of the bacteriorhodopsin photocycle

Deuterium kinetic isotope effects (KIE) were measured, and proton inventory plots were constructed, for the rates of reactions in the photocycles of wild-type bacteriorhodopsin and several site-specific mutants. Consistent with earlier reports from many groups, very large KIEs were observed for the third (and largest) rise component for the M state and for the decay of the O state, processes both linked to proton transfers in the extracellular region. The proton inventory plots (ratio of reaction rates in mixtures of H(2)O and D(2)O to that in H(2)O vs mole fraction of D(2)O) were approximately linear for the first and second M rise components and for M decay, as well as for O decay, indicating that the rates of these reactions are limited by simple proton transfer. Uniquely, the third rise component of M (and in the D96N mutant also a fourth rise component) exhibited a strongly curved proton inventory plot, suggesting that its rate, which largely accounts for the rate of deprotonation of the retinal Schiff base, depends on a complex multiproton process. This curvature is observed also in the E194Q, E204Q, and Y57F mutants but not in the R82A mutant. From these findings, and from the locations of bound water in the extracellular region in the crystal structure of the protein [Luecke, Schobert, Richter, Cartailler, and Lanyi (1999) J. Mol. Biol. 291, 899-911], we suspect that the effects of deuterium substitution on the formation of the M state originate from cooperative rearrangements of the extensively hydrogen-bonded water molecules 401, 402, and 406 near Asp-85 and Arg-82.     

Kinetics of enzyme acylation and deacylation in the penicillin acylase-catalyzed synthesis of beta-lactam antibiotics.

Penicillin acylase catalyses the hydrolysis and synthesis of semisynthetic beta-lactam antibiotics via formation of a covalent acyl-enzyme intermediate. The kinetic and mechanistic aspects of these reactions were studied. Stopped-flow experiments with the penicillin and ampicillin analogues 2-nitro-5-phenylacetoxy-benzoic acid (NIPAOB) and d-2-nitro-5-[(phenylglycyl)amino]-benzoic acid (NIPGB) showed that the rate-limiting step in the conversion of penicillin G and ampicillin is the formation of the acyl-enzyme. The phenylacetyl- and phenylglycyl-enzymes are hydrolysed with rate constants of at least 1000 s-1 and 75 s-1, respectively. A normal solvent deuterium kinetic isotope effect (KIE) of 2 on the hydrolysis of 2-nitro-5-[(phenylacetyl)amino]-benzoic acid (NIPAB), NIPGB and NIPAOB indicated that the formation of the acyl-enzyme proceeds via a general acid-base mechanism. In agreement with such a mechanism, the proton inventory of the kcat for NIPAB showed that one proton, with a fractionation factor of 0.5, is transferred in the transition state of the rate-limiting step. The overall KIE of 2 for the kcat of NIPAOB resulted from an inverse isotope effect at low concentrations of D2O, which is overridden by a large normal isotope effect at large molar fractions of D2O. Rate measurements in the presence of glycerol indicated that the inverse isotope effect originated from the higher viscosity of D2O compared to H2O. Deacylation of the acyl-enzyme was studied by nucleophile competition and inhibition experiments. The beta-lactam compound 7-aminodesacetoxycephalosporanic acid (7-ADCA) was a better nucleophile than 6-aminopenicillanic acid, caused by a higher affinity of the enzyme for 7-ADCA and complete suppression of hydrolysis of the acyl-enzyme upon binding of 7-ADCA. By combining the results of the steady-state, presteady state and nucleophile binding experiments, values for the relevant kinetic constants for the synthesis and hydrolysis of beta-lactam antibiotics were obtained.     

Environmentally coupled hydrogen tunneling. Linking catalysis to dynamics.

Many biological C-H activation reactions exhibit nonclassical kinetic isotope effects (KIEs). These nonclassical KIEs are too large (kH/kD > 7) and/or exhibit unusual temperature dependence such that the Arrhenius prefactor KIEs (AH/AD) fall outside of the semiclassical range near unity. The focus of this minireview is to discuss such KIEs within the context of the environmentally coupled hydrogen tunneling model. Full tunneling models of hydrogen transfer assume that protein or solvent fluctuations generate a reactive configuration along the classical, heavy-atom coordinate, from which the hydrogen transfers via nuclear tunneling. Environmentally coupled tunneling also invokes an environmental vibration (gating) that modulates the tunneling barrier, leading to a temperature-dependent KIE. These properties directly link enzyme fluctuations to the reaction coordinate for hydrogen transfer, making a quantum view of hydrogen transfer necessarily a dynamic view of catalysis. The environmentally coupled hydrogen tunneling model leads to a range of magnitudes of KIEs, which reflect the tunneling barrier, and a range of AH/AD values, which reflect the extent to which gating modulates hydrogen transfer. Gating is the primary determinant of the temperature dependence of the KIE within this model, providing insight into the importance of this motion in modulating the reaction coordinate. The potential use of variable temperature KIEs as a direct probe of coupling between environmental dynamics and the reaction coordinate is described. The evolution from application of a tunneling correction to a full tunneling model in enzymatic H transfer reactions is discussed in the context of a thermophilic alcohol dehydrogenase and soybean lipoxygenase-1.     

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.     

Kinetic isotope effect studies of the reaction catalyzed by uracil DNA glycosylase: evidence for an oxocarbenium ion-uracil anion intermediate

The DNA repair enzyme uracil DNA glycosylase catalyzes the first step in the uracil base excision repair pathway, the hydrolytic cleavage of the N-glycosidic bond of deoxyuridine in DNA. Here we report kinetic isotope effect (KIE) measurements that have allowed the determination of the transition-state structure for this important reaction. The small primary (13)C KIE (=1.010 +/- 0.009) and the large secondary alpha-deuterium KIE (=1.201 +/- 0.021) indicate that (i) the glycosidic bond is essentially completely broken in the transition state and (ii) there is significant sp(2) character at the anomeric carbon. Large secondary beta-deuterium KIEs were observed when [2'R-(2)H] = 1.102 +/- 0.011 and [2'S-(2)H] = 1.106 +/- 0.010. The nearly equal and large magnitudes of the two stereospecific beta-deuterium KIEs indicate strong hyperconjugation between the elongated glycosidic bond and both of the C2'-H2' bonds. Geometric interpretation of these beta-deuterium KIEs indicates that the furanose ring adopts a mild 3'-exo sugar pucker in the transition state, as would be expected for maximal stabilization of an oxocarbenium ion. Taken together, these results strongly indicate that the reaction proceeds through a dissociative transition state, with complete dissociation of the uracil anion followed by addition of water. To our knowledge, this is the first transition-state structure determined for enzymatic cleavage of the glycosidic linkage in a pyrimidine deoxyribonucleotide.     

On the mechanism of isomerization of α-hydroxypropionaldehyde to hydroxyacetone

Dimeric lactaldehyde and lactaldehyde-2-d in anhydrous or aqueous pyridine isomerize to hydroxyacetone. The isomerization is catalyzed by protic species that includes self-catalysis by the starting material and product, and catalysis by water and other Brönsted acids. The deuterium isotope effect which has values between 4 and 7, the isotopic exchange data, and the proton nmr assignment of the open dimer are consonant with the enediol mechanism with a rate-determining step involving the enolization of the open dimer.     

The monovalent cation-induced association of formyltetrahydrofolate synthetase subunits: a solvent isotope effect.

In the presence of specific monovalent cations (K+, Cs+, NH4+), inactive monomers of formyltetrahydrofolate synthetase associate to a catalytically active tetramer. The rate and extent of association of enzyme monomers prepared from C. cylindrosporum are enhanced 3.3-and about 50-fold, respectively, by the substitution of D2O for H2O. Both rate and equilibrium solvent isotope effects are due to a decrease in D2O of the dissociation constant of the monomer-cation complex. Analysis of rate and equilibria data obtained in solvent mixtures of varying deuterium/protium ratios indicates that the isotope effect may be due to the change in bonding of a single monomer proton during the association process. The data are most consistent with a model in which this proton is in a very weak potential in the cation-free monomer and is converted to a "normal" water-like proton in the monomer-cation complex.     

Hydrogen donor-acceptor fluctuations from kinetic isotope effects: a phenomenological model

Kinetic isotope effects (KIEs) and their temperature dependence can probe the structural and dynamic nature of enzyme-catalyzed proton or hydride transfers. The molecular interpretation of their temperature dependence requires expensive and specialized quantum mechanics/molecular mechanics (QM/MM) calculations to provide a quantitative molecular understanding. Currently available phenomenological models use a nonadiabatic assumption that is not appropriate for most hydride and proton-transfer reactions, while others require more parameters than the experimental data justify. Here we propose a phenomenological interpretation of KIEs based on a simple method to quantitatively link the size and temperature dependence of KIEs to a conformational distribution of the catalyzed reaction. This model assumes adiabatic hydrogen tunneling, and by fitting experimental KIE data, the model yields a population distribution for fluctuations of the distance between donor and acceptor atoms. Fits to data from a variety of proton and hydride transfers catalyzed by enzymes and their mutants, as well as nonenzymatic reactions, reveal that steeply temperature-dependent KIEs indicate the presence of at least two distinct conformational populations, each with different kinetic behaviors. We present the results of these calculations for several published cases and discuss how the predictions of the calculations might be experimentally tested. This analysis does not replace molecular QM/MM investigations, but it provides a fast and accessible way to quantitatively interpret KIEs in the context of a Marcus-like model.     

Kinetic isotope effect studies on the de novo rate of chromophore formation in fast- and slow-maturing GFP variants

The maturation process of green fluorescent protein (GFP) entails a protein oxidation reaction triggered by spontaneous backbone condensation. The chromophore is generated by full conjugation of the Tyr66 phenolic group with the heterocycle, a process that requires C-H bond scission at the benzylic carbon. We have prepared isotope-enriched protein bearing tyrosine residues deuterated at the beta carbon, and have determined kinetic isotope effects (KIEs) on the GFP self-processing reaction. Progress curves for the production of H 2O 2 and the mature chromophore were analyzed by global curve fitting to a three-step mechanism describing preoxidation, oxidation and postoxidation events. Although a KIE for protein oxidation could not be discerned ( k H/ k D = 1.1 +/- 0.2), a full primary KIE of 5.9 (+/-2.8) was extracted for the postoxidation step. Therefore, the exocyclic carbon is not involved in the reduction of molecular oxygen. Rather, C-H bond cleavage proceeds from the oxidized cyclic imine form, and is the rate-limiting event of the final step. Substantial pH-dependence of maturation was observed upon substitution of the catalytic glutamate (E222Q), indicating an apparent p K a of 9.4 (+/-0.1) for the base catalyst. For this variant, a KIE of 5.8 (+/-0.4) was determined for the intrinsic time constant that is thought to describe the final step, as supported by ultra-high resolution mass spectrometric results. The data are consistent with general base catalysis of the postoxidation events yielding green color. Structural arguments suggest a mechanism in which the highly conserved Arg96 serves as catalytic base in proton abstraction from the Tyr66-derived beta carbon.     

Transition state analysis for human and Plasmodium falciparum purine nucleoside phosphorylases

Recent studies have shown that Plasmodium falciparum is sensitive to a purine salvage block at purine nucleoside phosphorylase (PNP) and that human PNP is a target for T-cell proliferative diseases. Specific tight-binding inhibitors might be designed on the basis of specific PNP transition state structures. Kinetic isotope effects (KIEs) were measured for arsenolysis of inosine catalyzed by P. falciparum and human purine nucleoside phosphorylases. Intrinsic KIEs from [1'-(3)H]-, [2'-(3)H]-, [1'-(14)C]-, [9-(15)N]-, and [5'-(3)H]inosines were 1.184 +/- 0.004, 1.031 +/- 0.004, 1.002 +/- 0.006, 1.029 +/- 0.006, and 1.062 +/- 0.002 for the human enzyme and 1.116 +/- 0.007, 1.036 +/- 0.003, 0.996 +/- 0.006, 1.019 +/- 0.005, and 1.064 +/- 0.003 for P. falciparum PNPs, respectively. Analysis of KIEs indicated a highly dissociative D(N)A(N) (S(N)1) stepwise mechanism with very little leaving group involvement. The near-unity 1'-(14)C KIEs for both human and P. falciparum PNP agree with the theoretical value for a 1'-(14)C equilibrium isotope effect for oxacarbenium ion formation when computed at the B1LYP/6-31G(d) level of theory. The 9-(15)N KIE for human PNP is also in agreement with theory for equilibrium formation of hypoxanthine and oxacarbenium ion at this level of theory. The 9-(15)N KIE for P. falciparum PNP shows a constrained vibrational environment around N9 at the transition state. A relatively small beta-secondary 2'-(3)H KIE for both enzymes indicates a 3'-endo conformation for ribose and relatively weak hyperconjugation at the transition state. The large 5'-(3)H KIE reveals substantial distortion at the 5'-hydroxymethyl group which causes loosening of the C5'-H5' bonds during the reaction coordinate.     

A method to determine 18 O kinetic isotope effects in the hydrolysis of nucleotide triphosphates

A method to determine 18 O kinetic isotope effects (KIEs) in the hydrolysis of GTP that is generally applicable to reactions involving other nucleotide triphosphates is described. Internal competition, where the substrate of the reaction is a mixture of 18 O-labeled and unlabeled nucleotides, is employed, and the change in relative abundance of the two species in the course of the reaction is used to calculate KIE. The nucleotide labeled with 18 O at sites of mechanistic interest also contains 13C at all carbon positions, whereas the 16 O-labeled nucleotide is depleted of 13C. The relative abundance of the labeled and unlabeled substrates or products is reflected in the carbon isotope ratio (13C/12C) in GTP or GDP, which is determined by the use of a liquid chromatography-coupled isotope ratio mass spectrometer (LC-coupled IRMS). The LC is coupled to the IRMS by an Isolink interface. Carbon isotope ratios can be determined with accuracy and precision greater than 0.04% and are consistent over an order of magnitude in sample amount. KIE values for Ras/NF1(333)-catalyzed hydrolysis of [beta18 O3,13C]GTP were determined by change in the isotope ratio of GTP or GDP or the ratio of the isotope ratio of GDP to that of GTP. KIE values computed in the three ways agree within 0.1%, although the method using the ratio of isotope ratios of GDP and GTP gives superior precision (<0.1%). A single KIE measurement can be conducted in 25 min with less than 5 microg nucleotide reaction product.     

Vibrationally enhanced tunneling as a mechanism for enzymatic hydrogen transfer.

We present a theory of enzymatic hydrogen transfer in which hydrogen tunneling is mediated by thermal fluctuations of the enzyme's active site. These fluctuations greatly increase the tunneling rate by shortening the distance the hydrogen must tunnel. The average tunneling distance is shown to decrease when heavier isotopes are substituted for the hydrogen or when the temperature is increased, leading to kinetic isotope effects (KIEs)--defined as the factor by which the reaction slows down when isotopically substituted substrates are used--that need be no larger than KIEs for nontunneling mechanisms. Within this theory we derive a simple KIE expression for vibrationally enhanced ground state tunneling that is able to fit the data for the bovine serum amine oxidase (BSAO) system, correctly predicting the large temperature dependence of the KIEs. Because the KIEs in this theory can resemble those for nontunneling dynamics, distinguishing the two possibilities requires careful measurements over a range of temperatures, as has been done for BSAO.