Theoretical Kinetic Isotope Effects in Establishing the Precise Biodegradation Mechanisms of Organic Pollutants.

Testosterone 6beta-hydroxylation is a prototypic reaction of cytochrome P450 (P450) 3A4, the major human P450. Biomimetic reactions produced a variety of testosterone oxidation products with 6beta-hydroxylation being only a minor reaction, indicating that P450 3A4 has considerable control over the course of steroid hydroxylation because 6beta-hydroxylation is not dominant in a thermodynamically controlled oxidation of the substrate. Several isotopically labeled testosterone substrates were prepared and used to probe the catalytic mechanism of P450 3A4: (i) 2,2,4,6,6-(2)H(5); (ii) 6,6-(2)H(2); (iii) 6alpha-(2)H; (iv) 6beta-(2)H; and (v) 6beta-(3)H testosterone. Only the 6beta-hydrogen was removed by P450 3A4 and not the 6alpha, indicating that P450 3A4 abstracts hydrogen and rebounds oxygen only at the beta face. Analysis of the rates of hydroxylation of 6beta-(1)H-, 6beta-(2)H-, and 6beta-(3)H-labeled testosterone and application of the Northrop method yielded an apparent intrinsic kinetic deuterium isotope effect ((D)k) of 15. The deuterium isotope effects on k(cat) and k(cat)/K(m) in non-competitive reactions were only 2-3. Some "switching" to other hydroxylations occurred because of 6beta-(2)H substitution. The high (D)k value is consistent with an initial hydrogen atom abstraction reaction. Attenuation of the high (D)k in the non-competitive experiments implies that C-H bond breaking is not a dominant rate-limiting step. Considerable attenuation of a high (D)k value was also seen with a slower P450 3A4 reaction, the O-dealkylation of 7-benzyloxyquinoline. Thus P450 3A4 is an enzyme with regioselective flexibility but also considerable regioselectivity and stereoselectivity in product formation, not necessarily dominated by the ease of C-H bond breaking.

Compound-specific isotope analysis was used to monitor the delta(13)C signature of chloroform produced upon the chlorination of model compounds representing natural organic matter functional groups (resorcinol, acetylacetone, acetophenone, phenol, and 2,4,6-trichlorophenol) and a natural water sample. For each model compound, a different apparent kinetic isotope effect was found for chloroform formation. Normal isotope effects were found for resorcinol, acetylacetone, and acetophenone, and ranged from 1.009 +/- 0.002 to 1.024 +/- 0.004. For the two phenols, an inverse effect was found (0.980 +/- 0.004). Lake Zürich water also had a inverse effect (0.997 +/- < 0.001) indicating that phenols are likely chloroform precursors in NOM, but that other functional groups may also participate. The apparent 13C kinetic isotope effect for the addition/ elimination reaction of 1,1,1-trichloropropanone mediated by OH- to yield chloroform is 1.014 +/- 0.002. A comparison of this value to those found for the chlorination of the model precursors and an evaluation of the differences in chloroform production kinetics for the different model precursors argue against a mechanism in which all NOM precursors react via a common intermediate. Compound specific isotope analysis may give additional insights into chloroform formation mechanisms beyond those allowed by current techniques.

Kinetic isotope effects (KIEs) provide a powerful tool to interrogate transition states of both enzymic and non-enzymic reactions, provided that one can measure the intrinsic KIE on the chemical step of interest, that is, the isotope effect undiminished by other isotope-insensitive steps that may contribute to the overall rate of reaction. For small molecules reacting in solution the KIE measured usually represents the intrinsic value; however, for enzymes this is seldom the case. Here we report the first measurement of the intrinsic KIE in an adenosylcobalamin enzyme (AdoCbl, coenzyme B12) for hydrogen atom transfer from substrate to coenzyme, which is a key step in the mechanism of this class of enzymes. For the B12 enzyme glutamate mutase the intrinsic deuterium KIE for hydrogen transfer from the substrate, (2S,3S)-3-methylaspartate, to 5’-deoxyadensosine is 4.1. This value is well within the semiclassical limit for a deuterium isotope effect and is much smaller than the anomalously large KIEs previously measured in other B12 enzymes and non-enzymatic model reactions, which were attributed to extensive hydrogen tunneling. Glutamate mutase one is of a group of AdoCbl-dependent enzymes that catalyze unusual isomerization reactions that formally involve a 1,2 hydrogen atom migration and proceed through a mechanism involving carbon-based free radical intermediates (Scheme 1). Radicals are generated by homolysis of the reactive cobalt–carbon bond of the coenzyme to form cob(II)alamin, a cobalt(II) intermediate, and the 5’-deoxyadenosyl radical. The adenosyl radical then abstracts the migrating hydrogen from the substrate to form 5’-deoxyadenosine and the substrate radical. The substrate radical next undergoes rearrangement to give the product radical, which is then quenched by hydrogen transfer from 5’deoxyadenosine to give the product and regenerate the 5’deoxyadenosyl radical. Finally, recombination of the adenosyl radical and cob(II)alamin to reform the coenzyme completes the catalytic cycle. Our interest in the mechanisms by which enzymes generate free radicals, as exemplified by dependent glutamate mutase, led us to undertake an extensive set of KIE measurements to examine how hydrogen abstraction from the substrate and coenzyme homolysis are coupled together. KIE measurements using deuteriumand tritium-labeled substrates and coenzyme have proved especially informative probes of the key steps of Co C bond homolysis and hydrogen atom abstraction from substrate. Pre-steady-state measurements on a number of enzymes have shown that hydrogen abstraction is kinetically coupled to Co C bond homolysis, as evidenced by the appearance of a kinetic isotope effect on cobalt–carbon bond homolysis when the enzymes are reacted with deuterated substrates. This observation implies that the 5’-dA radical is a highenergy intermediate that only has a fleeting existence. Furthermore, the KIEs reported for several AdoCbl enzymes are extremely large (ranging from 10 to 50), which has generally been attributed to hydrogen tunneling. In particular, extensive hydrogen tunneling in methylmalonyl-CoA mutase has been deduced from the temperature dependence of the deuterium isotope effect on hydrogen transfer and Co C bond homolysis. The KIEs discussed above were all measured indirectly by UV/Vis stopped-flow spectroscopy, using the spectroscopic changes associated with Co C bond homolysis as a convenient reporter of the kinetics. Although we originally reported a very large primary deuterium isotope effect, suggestive of tunneling, for the formation of cob(II)alamin and 5’-dA upon reaction of hologlutamate mutase with deuterated substrates (kH/kD= 28 with glutamate and 35 with methylaspartate), [11]

Bacterial protein toxins are the most powerful human poisons known, exhibiting an LD(50) of 0.1-1 ng kg(-)(1). A major subset of such toxins is the NAD(+)-dependent ADP-ribosylating exotoxins, which include pertussis, cholera, and diphtheria toxin. Diphtheria toxin catalyzes the ADP ribosylation of the diphthamide residue of eukaryotic elongation factor 2 (eEF-2). The transition state of ADP ribosylation catalyzed by diphtheria toxin has been characterized by measuring a family of kinetic isotope effects using (3)H-, (14)C-, and (15)N-labeled NAD(+) with purified yeast eEF-2. Isotope trapping experiments yield a commitment to catalysis of 0.24 at saturating eEF-2 concentrations, resulting in suppression of the intrinsic isotope effects. Following correction for the commitment factor, intrinsic primary kinetic isotope effects of 1.055 +/- 0.003 and 1.022 +/- 0.004 were observed for [1(N)'-(14)C]- and [1(N)-(15)N]NAD(+), respectively; the double primary isotope effect was 1.066 +/- 0.004 for [1(N)'-(14)C, 1(N)-(15)N]NAD(+). Secondary kinetic isotope effects of 1.194 +/- 0.002, 1.101 +/- 0.003, 1.013 +/- 0.005, and 0.988 +/- 0.002 were determined for [1(N)'-(3)H]-, [2(N)'-(3)H]-, [4(N)'-(3)H]-, and [5(N)'-(3)H]NAD(+), respectively. The transition state structure was modeled using density functional theory (B1LYP/6-31+G) as implemented in Gaussian 98, and theoretical kinetic isotope effects were subsequently calculated using Isoeff 98. Constraints were varied in a systematic manner until the calculated kinetic isotope effects matched the intrinsic isotope effects. The transition state model most consistent with the intrinsic isotope effects is characterized by the substantial loss in bond order of the nicotinamide leaving group (bond order = 0.18, 1.99 A) and weak participation of the attacking imidazole nucleophile (bond order = 0.03, 2.58 A). The transition state structure imparts strong oxacarbenium ion character to the ribose ring even though significant bond order remains to the nicotinamide leaving group. The transition state model presented here is asymmetric and consistent with a dissociative S(N)1 type mechanism in which attack of the diphthamide nucleophile lags behind departure of the nicotinamide.

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.

&lt;p&gt;The nucleotide 2'-deoxythymidine 5'-monophosphate (thymidylate, dTMP) is phosphorylated twice to become a substrate for DNA polymerases, which copy a cell’s genetic information in advance of cell division. The main route to dTMP is mediated by the enzyme thymidylate synthase (TSase) and goes through 2'-deoxyuridine 5'-monophosphate (dUMP); dUMP’s heterocyclic aromatic pyrimidine ring loses a proton from its C5 position and gains a methylene and a hydride from the other reactant, methylene tetrahydrofolate (MTHF). In general, intricate knowledge of an enzyme’s mechanism can yield insight that leads to the development of precision-targeted inhibitors tailored exactly to thymidylate synthase. In fact, even more careful targeting could be achievable: Although E. coli TSase has served as a model system, investigators have increasingly been directing their lines of inquiry toward human TSase.&lt;/p&gt; &lt;p&gt;A general enzymatic catalytic cascade is complex, comprising substrate binding, the chemical steps and product release; typically, the product release step is rate-limiting. TSase, however, is partially rate-limited by the chemistry portion of the process. The enzymatic mechanism has been considered for decades, yet recently has undergone a reassessment. After substrate binding – for which there is strong evidence for preference to dUMP as the first ligand in the wild-type E. coli enzyme – the important events are methylene transfer from MTHF to dUMP, proton abstraction and hydride transfer. The last of these – hydride transfer – is irreversible and rate-limiting (to a large degree without Mg2+, and to a small but noticeable degree with Mg2+). The studies described here are aimed at three therapeutically relevant questions: (a) determining the extent of negative charge accumulation at the O4 position of the hydride transfer acceptor; (b) expanding knowledge of the differential properties of E. coli and human TSase; and (c) gaining insight into the molecular origin of the drug resistance seen in a clinically relevant human TSase mutant.&lt;/p&gt; &lt;p&gt;The properties touched on in this work include steady-state kinetics; inhibition constants toward 5-fluoro dUMP, substrate binding sequence and the temperature dependency of intrinsic hydride transfer kinetic isotope effects (KIEs). Intrinsic KIEs are a specialized measurement that permits the investigator to examine a particular hydrogen transfer step in isolation; it is achieved by labeling the bond to hydrogen broken in the reaction with protium (1H, also written as H), deuterium (2H, also written as D) or tritium (3H, also written as T). The latter is radioactive. The reaction is conducted with a mixture of two hydrogen isotopes at a time, and the extent to which the heavier isotope is disfavored against reaction is assessed; this covers multiple steps. Heavier isotopes directly participating in a chemical step react slower both because of zero-point vibrational energies if a semi-classical view is taken and because of the mass-dependence of tunneling probabilities if a quantum-mechanical view is taken. Each of the two-way isotopic comparisons mentioned above furnishes an observed KIE for that competition between two isotopes. Mathematical combination of two isotopic comparisons cancels out the effect of isotopically insensitive steps and provides rich insight into the hydride transfer alone. The ultimate result is the ratio of rate constants for the isotopologues; this ratio’s magnitude and variation with temperature report on the compactness of the active site and its resistance to thermal fluctuation, respectively.&lt;/p&gt; &lt;p&gt;Our results reveal a possible role for E. coli asparagine 177 (N177) in the hydride transfer transition state (TS) stabilization, as revealed by its disruption in the aspartate mutant, N177D. This disruption was found to be alleviated to a high extent when the substrate was changed to dCMP, consistent with the N177 stabilizing partial negative charge at the TS for hydride transfer. This has drug design implications. Our work on human TSase underscores slightly weaker substrate binding preference, insensitivity to Mg2+ and mild alteration of hydride transfer TS when compared with E. coli TSase. Finally, analysis of the Y33H mutant of human TSase – the affected residue being remote from the active site – indicated the drug resistance was because of a higher inhibition constant for 5F-dUMP and that the hydride transfer step is disrupted, with a wider variation among donor-acceptor distances (between the two carbons involved in the hydride transfer at the TS for that step). Other researchers’ crystallographic evidence reveals greater positional uncertainty for a set of active-site side chains in the E. coli equivalent mutant. In totality, the data available implicate enzyme motions as relevant to drug binding and to catalysis for human TSase.&lt;/p&gt; &lt;p&gt;In summary, the research described herein enriches the understanding of several aspects of the behavior of multiple TSase variants – the overall performance as seen via steady-state kinetics; the pattern of substrate binding as seen with observed KIEs for the proton abstraction step; and the efficiency of active site preparation for hydride transfer as evidenced in the temperature dependency of intrinsic hydride transfer KIEs.&lt;/p&gt;

7-Ethoxy (OEt) coumarin has been used as a model substrate in many cytochrome P450 (P450) studies, including the use of kinetic isotope effects to probe facets of P450 kinetics. P450s 1A2 and 2E1 are known to be the major catalysts of 7-OEt coumarin O-deethylation in human liver microsomes. Human P450 1A2 also catalyzed 3-hydroxylation of 7-methoxy (OMe) coumarin at appreciable rates but P450 2E1 did not. Intramolecular kinetic isotope effects were used as estimates of the intrinsic kinetic deuterium isotope effects for both 7-OMe and 7-OEt coumarin dealkylation reactions. The apparent intrinsic isotope effect for P450 1A2 (9.4 for O-demethylation, 6.1 for O-deethylation) showed little attenuation in other competitive and noncompetitive experiments. With P450 2E1, the intrinsic isotope effect (9.6 for O-demethylation, 6.1 for O-deethylation) was attenuated in the noncompetitive intermolecular experiments. High noncompetitive intermolecular kinetic isotope effects were seen for 7-OEt coumarin O-deethylation in a baculovirus-based microsomal system and five samples of human liver microsomes (7.3-8.1 for O-deethylation), consistent with the view that P450 1A2 is the most efficient P450 catalyzing this reaction in human liver microsomes and indicating that the C-H bond-breaking step makes a major contribution to the rate of this P450 (1A2) reaction. Thus, the rate-limiting step appears to be the chemistry of the breaking of this bond by the activated iron-oxygen complex, as opposed to steps involved in the generation of the reactive complex. The conclusion about the rate-limiting step applies to all of the systems studied with this model P450 1A2 reaction including human liver microsomes, the most physiologically relevant.

Measuring stable isotope fractionation of carbon, hydrogen, and other elements by Compound Specific Isotope Analysis (CSIA) is a new, innovative approach to assess organic pollutant degradation in the environment. Central to this concept is the Rayleigh equation which relates degradation-induced decreases in concentrations directly to concomitant changes in bulk (= average over the whole compound) isotope ratios. The extent of in situ transformation may therefore be inferred from measured isotope ratios in field samples, provided that an appropriate enrichment factor (epsilonbulk) is known. This epsilonbulk value, however, is usually only valid for a specific compound and for specific degradation conditions. Therefore, a direct comparison of epsilonbulk values for different compounds and for different types of reactions has in general not been feasible. In addition, it is often uncertain how robust and reproducible epsilonbulk values are and how confidently they can be used to quantify contaminant degradation in the field. To improve this situation and to achieve a more in-depth understanding, this critical review aims to relate fundamental insight about kinetic isotope effects (KIE) found in the physico(bio)chemical literature to apparent kinetic isotope effects (AKIE) derived from epsilonbulk values reported in environmentally oriented studies. Starting from basic rate laws, a quite general derivation of the Rayleigh equation is given, resulting in a novel set of simple equations that take into account the effects of (1) nonreacting positions and (2) intramolecular competition and that lead to position-specific AKIE values rather than bulk enrichment factors. Reevaluation of existing epsilonbulk literature values result in consistent ranges of AKIE values that generally are in good agreement with previously published data in the (bio)-chemical literature and are typical of certain degradation reactions (subscripts C and H indicate values for carbon and hydrogen): AKIEc = 1.01-1.03 and AKIEH = 2-23 for oxidation of C-H bonds; AKIEc = 1.03-1.07 for SN2-reactions; AKIEc = 1.02-1.03 for reductive cleavage of C-Cl bonds; AKIEc = 1.00-1.01 for C=C bond epoxidation; AKIEc = 1.02-1.03 for C=C bond oxidation by permanganate. Hence, the evaluation scheme presented bridges a gap between basic and environmental (bio)chemistry and provides insight into factors that control the magnitude of bulk isotope fractionation factors. It also serves as a basis to identify degradation pathways using isotope data. It is shown how such an analysis may be even possible in complex field situations and/or in cases where AKIE values are smaller than intrinsic KIE values, provided that isotope fractionation is measured for two elements simultaneously ("two-dimensional isotope analysis"). Finally, the procedure is used (1) to point outthe possibility of estimating approximate epsilonbulk values for new compounds and (2) to discuss the moderate, but non-negligible variability that may quite generally be associated with epsilonbulk values. Future research is suggested to better understand and take into account the various factors that may cause such variability.

The pyridoxal phosphate dependent alanine racemase catalyzes the interconversion of L- and D-alanine. The latter is an essential component of peptidoglycan in cell walls of Gram-negative and -positive bacteria, making alanine racemase an attractive target for antibacterials. Global analysis of protiated and deuterated progress curves simultaneously enables determination of intrinsic kinetic and equilibrium isotope effects for alanine racemase. The intrinsic primary kinetic isotope effects for Calpha hydron abstraction are 1.57 +/- 0.05 in the D --> L direction and 1.66 +/- 0.09 in the L --> D direction. Secondary kinetic isotope effects were found for the external aldimine formation steps in both the L --> D (1.13 +/- 0.05, forward; 0.90 +/- 0.03, reverse) and D --> L (1.13 +/- 0.06, forward; 0.89 +/- 0.03, reverse) directions. The secondary equilibrium isotope effects calculated from these are 1.26 +/- 0.07 and 1.27 +/- 0.07 for the L --> D and D --> L directions, respectively. These equilibrium isotope effects imply substantial ground-state destabilization of the C-H bond via hyperconjugation with the conjugated Schiff base/pyridine ring pi system. The magnitudes of the intrinsic primary kinetic isotope effects, the lower boundary on the energy of the quinonoid intermediate, and the protonation states of the active site catalytic acids/bases (K39-epsilonNH2 and Y265-OH) suggest that the pKa of the substrate Calpha-H bond in the external aldimine lies between those of the two catalytic bases, such that the proton abstraction transition state is early in the D --> L direction and late in the L --> D direction.

Kinetic isotope effects (KIEs) and computer modeling using density functional theory were used to approximate the transition state of human 5'-methylthioadenosine phosphorylase (MTAP). KIEs were measured on the arsenolysis of 5'-methylthioadenosine (MTA) catalyzed by MTAP and were corrected for the forward commitment to catalysis. Intrinsic KIEs were obtained for [1'-(3)H], [1'-(14)C], [2'-(3)H], [4'-(3)H], [5'-(3)H(2)], [9-(15)N], and [Me-(3)H(3)] MTAs. The primary intrinsic KIEs (1'-(14)C and 9-(15)N) suggest that MTAP has a dissociative S(N)1 transition state with its cationic center at the anomeric carbon and insignificant bond order to the leaving group. The 9-(15)N intrinsic KIE of 1.039 also establishes an anionic character for the adenine leaving group, whereas the alpha-primary 1'-(14)C KIE of 1.031 indicates significant nucleophilic participation at the transition state. Computational matching of the calculated EIEs to the intrinsic isotope effects places the oxygen nucleophile 2.0 Angstrom from the anomeric carbon. The 4'-(3)H KIE is sensitive to the polarization of the 3'-OH group. Calculations suggest that a 4'-(3)H KIE of 1.047 is consistent with ionization of the 3'-OH group, indicating formation of a zwitterion at the transition state. The transition state has cationic character at the anomeric carbon and is anionic at the 3'-OH oxygen, with an anionic leaving group. The isotope effects predicted a 3'-endo conformation for the ribosyl zwitterion, corresponding to a H1'-C1'-C2'-H2' torsional angle of 33 degrees. The [Me-(3)H(3)] and [5'-(3)H(2)] KIEs arise predominantly from the negative hyperconjugation of the lone pairs of sulfur with the sigma (C-H) antibonding orbitals. Human MTAP is characterized by a late S(N)1 transition state with significant participation of the phosphate nucleophile.

Kinetic isotope effects (KIEs) and computer modeling are used to approximate the transition state of S. pneumoniae 5'-methylthioadenosine/S-adenosylhomocysteine nucleosidase (MTAN). Experimental KIEs were measured and corrected for a small forward commitment factor. Intrinsic KIEs were obtained for [1'-3H], [1'-14C], [2'-3H], [4'-3H], [5'-3H(2)], [9-15N] and [Me-3H(3)] MTAs. The intrinsic KIEs suggest an SN1 transition state with no covalent participation of the adenine or the water nucleophile. The transition state was modeled as a stable ribooxacarbenium ion intermediate and was constrained to fit the intrinsic KIEs. The isotope effects predicted a 3-endo conformation for the ribosyl oxacarbenium-ion corresponding to H1'-C1'-C2'-H2' dihedral angle of 70 degrees. Ab initio Hartree-Fock and DFT calculations were performed to study the effect of polarization of ribosyl hydroxyls, torsional angles, and the effect of base orientation on isotope effects. Calculations suggest that the 4'-3H KIE arises from hyperconjugation between the lonepair (n(p)) of O4' and the sigma* (C4'-H4') antibonding orbital owing to polarization of the 3'-hydroxyl by Glu174. A [methyl-3H(3)] KIE is due to hyperconjugation between np of sulfur and sigma* of methyl C-H bonds. The van der Waal contacts increase the 1'-3H KIE because of induced dipole-dipole interactions. The 1'-3H KIE is also influenced by the torsion angles of adjacent atoms and by polarization of the 2'-hydroxyl. Changing the virtual solvent (dielectric constant) does not influence the isotope effects. Unlike most N-ribosyltransferases, N7 of the leaving group adenine is not protonated at the transition state of S. pneumoniae MTAN. This feature differentiates the S. pneumoniae and E. coli transition states and explains the 10(3)-fold decrease in the catalytic efficiency of S. pneumoniae MTAN relative to that from E. coli.

Methylthioadenosine/S-adenosylhomocysteine nucleosidase (MTAN) catalyzes reactions linked to polyamine metabolism, quorum sensing pathways, methylation reactions, and adenine salvage. It is a candidate target for antimicrobial drug design. Kinetic isotope effects (KIEs) were measured on the MTAN-catalyzed hydrolysis of 5'-methylthioadenosine (MTA) to determine the transition state structure. KIEs measured at pH 7.5 were near unity due to the large forward commitment to catalysis. Intrinsic KIEs were expressed by increasing the pH to 8.5. Intrinsic KIEs from MTAs labeled at 1'-(3)H, 1'-(14)C, 2'-(3)H, 4'-(3)H, 5'-(3)H, 9-(15)N, and Me-(3)H(3) were 1.160 +/- 0.004, 1.004 +/- 0.003, 1.044 +/- 0.004, 1.015 +/- 0.002, 1.010 +/- 0.002, 1.018 +/- 0.006, and 1.051 +/- 0.002, respectively. The large 1'-(3)H and small 1'-(14)C KIEs indicate that the Escherichia coli MTAN reaction undergoes a dissociative (D(N)A(N)) (S(N)1) mechanism with little involvement of the leaving group or participation of the attacking nucleophile at the transition state, causing the transition state to have significant ribooxacarbenium ion character. A transition state constrained to match the intrinsic KIEs was located with density functional theory [B3LYP/6-31G(d,p)]. The leaving group (N9) is predicted to be 3.0 A from the anomeric carbon. The small beta-secondary 2'-(3)H KIE of 1.044 corresponds to a modest 3'-endo conformation for ribose and a H1'-C1'-C2'-H2' dihedral angle of 53 degrees at the transition state. Natural bond orbital analysis of the substrate and the transition state suggests that the 4'-(3)H KIE is due to hyperconjugation between the lone pair (n(p)) of O3' and the antibonding (sigma) orbital of the C4'-H4' group, and the methyl-(3)H(3) KIE is due to hyperconjugation between the n(p) of sulfur and the sigma of methyl C-H bonds. Transition state analogues that resemble this transition state structure are powerful inhibitors, and their molecular electrostatic potential maps closely resemble that of the transition state.

The apparent kinetic isotope effect (KIE) for a multistepsteady-statereaction can be expressed simply as a sum of terms, one for each transitionstate (TS) in the serial sequence, each of which is the product ofthe KIE for an individual TS (with respect to a common reference state)and a weighting factor, which is the degree of kinetic significancefor that TS. This requires knowledge of the relative Gibbs energiesof the sequential TSs but not of any intermediates, and it involvesa much simpler expression than the conventional method for analysisof KIEs for enzyme reactions. A numerical example is presented toillustrate how the same apparent KIE may result from numerous combinationsof individual KIEs and weighting factors. It is proposed that computedapparent KIEs should be compared directly with experimentally observedKIEs rather than with derived intrinsic KIEs of possibly dubious validity.The results of DFT calculations for an SN1 nucleophilicdisplacement are presented to show how the apparent KIE varies, asthe relative concentration of the nucleophilic species ranges from0.1 to 10, between limiting values corresponding to either the firstor second step being completely rate limiting.

Natural abundance kinetic isotope effects: expt. vs theory.

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.