Ensemble‐Averaged QM/MM Kinetic Isotope Effects for the SN2 Reaction of Cyanide Anions with Chloroethane in DMSO Solution

The existence of solvent fluctuations leads to populations of reactant‐state (RS) and transition‐state (TS) configurations and implies that property calculations must include appropriate averaging over distributions of values for individual configurations. Average kinetic isotope effects 〈KIE〉 for NC^?+EtCl→NCEt+Cl^? in DMSO solution at 30 °C are best obtained as the ratio 〈f_RS〉/〈f_TS〉 of isotopic partition function ratios separately averaged over all RS and TS configurations. In this way the hybrid AM1/OPLS‐AA potential yields 〈KIE〉 values for all six isotopic substitutions (2° α‐²H₂, 2° β‐²H₃, α‐¹¹C/¹⁴C, leaving group ³⁷Cl, and nucleophile ¹³C and ¹⁵N) for this reaction in the correct direction as measured experimentally. These thermally‐averaged calculated KIEs may be compared meaningfully with experiment, and only one of them differs in magnitude from the experimental value by more than one standard deviation from the mean. This success contrasts with previous KIE calculations based upon traditional methods without averaging. The isotopic partition function ratios are best evaluated using all (internal) vibrational and (external) librational frequencies obtained from Hessians determined for subsets of atoms, relaxed to local minima or saddle points, within frozen solvent environments of structures sampled along molecular dynamics trajectories for RS and TS. The current method may perfectly well be implemented with other QM or QM/MM methods, and thus provides a useful tool for investigating KIEs in relation to studies of chemical reaction mechanisms in solution or catalyzed by enzymes.     

A stepwise solvent-promoted SNi reaction of α-D-glucopyranosyl fluoride: mechanistic implications for retaining glycosyltransferases

The solvolysis of α-d-glucopyranosyl fluoride in hexafluoro-2-propanol gives two products, 1,1,1,3,3,3-hexafluoropropan-2-yl α-d-glucopyranoside and 1,6-anhydro-β-D-glucopyranose. The ratio of these two products is essentially unchanged for reactions that are performed between 56 and 100 °C. The activation parameters for the solvolysis reaction are as follows: ΔH(++) = 81.4 ± 1.7 kJ mol(-1), and ΔS(++) = -90.3 ± 4.6 J mol(-1) K(-1). To characterize, by use of multiple kinetic isotope effect (KIE) measurements, the TS for the solvolysis reaction in hexafluoro-2-propanol, we synthesized a series of isotopically labeled α-d-glucopyranosyl fluorides. The measured KIEs for the C1 deuterium, C2 deuterium, C5 deuterium, anomeric carbon, ring oxygen, O6, and solvent deuterium are 1.185 ± 0.006, 1.080 ± 0.010, 0.987 ± 0.007, 1.008 ± 0.007, 0.997 ± 0.006, 1.003 ± 0.007, and 1.68 ± 0.07, respectively. The transition state for the solvolysis reaction was modeled computationally using the experimental KIE values as constraints. Taken together, the reported data are consistent with the retained solvolysis product being formed in an S(N)i (D(N)(++)*A(Nss)) reaction with a late transition state in which cleavage of the glycosidic bond is coupled to the transfer of a proton from a solvating hexafluoro-2-propanol molecule. In comparison, the inverted product, 1,6-anhydro-β-D-glucopyranose, is formed by intramolecular capture of a solvent-equilibrated glucopyranosylium ion, which results from dissociation of the solvent-separated ion pair formed in the rate-limiting ionization reaction (D(N)(++) + A(N)). The implications that this model reaction have for the mode of action of retaining glycosyltransferases are discussed.     

Deuterium kinetic isotope effects on the dissociation of a protein-fatty acid complex in the gas phase

Deuterium kinetic isotope effects (KIEs) are reported for the first time for the dissociation of a protein-ligand complex in the gas phase. Temperature-dependent rate constants were measured for the loss of neutral ligand from the deprotonated ions of the 1:1 complex of bovine β-lactoglobulin (Lg) and palmitic acid (PA), (Lg + PA)(n-) → Lg(n-) + PA, at the 6- and 7- charge states. At 25 °C, partial or complete deuteration of the acyl chain of PA results in a measurable inverse KIE for both charge states. The magnitude of the KIEs is temperature dependent, and Arrhenius analysis of the rate constants reveals that deuteration of PA results in a decrease in activation energy. In contrast, there is no measurable deuterium KIE for the dissociation of the (Lg + PA) complex in aqueous solution at pH 8. Deuterium KIEs were calculated using conventional transition-state theory with an assumption of a late dissociative transition state (TS), in which the ligand is free of the binding pocket. The vibrational frequencies of deuterated and non-deuterated PA in the gas phase and in various solvents (n-hexane, 1-chlorohexane, acetone, and water) were established computationally. The KIEs calculated from the corresponding differences in zero-point energies account qualitatively for the observation of an inverse KIE but do not account for the magnitude of the KIEs nor their temperature dependence. It is proposed that the dissociation of the (Lg + PA) complex in aqueous solution also proceeds through a late TS in which the acyl chain is extensively hydrated such that there is no significant differential change in the vibrational frequencies along the reaction coordinate and, consequently, no significant KIE.     

Measurements of the 12C/13C kinetic isotope effects in the gas-phase reactions of light alkanes with chlorine atoms

The carbon kinetic isotope effects (KIEs) of the reactions of several light non-methane hydrocarbons (NMHC) with Cl atoms were determined at room temperature and ambient pressure. All measured KIEs, defined as the ratio of the Cl reaction rate constants of the light isotopologue over that of the heavy isotopologue (Clk12/Clk13) are greater than unity or normal KIEs. For simplicity, measured KIEs are reported in per mil according to Clepsilon=(Clk12/Clk13 -1)x1000 per thousand unless noted otherwise. The following average KIEs were obtained (all in per thousand): 10.73+/-0.20 (ethane), 6.44+/-0.14 (propane), 6.18+/-0.18 (methylpropane), 3.94+/-0.01 (n-butane), 1.79+/-0.42 (methylbutane), 3.22+/-0.17 (n-pentane), 2.02+/-0.40 (n-hexane), 2.06+/-0.19 (n-heptane), 1.54+/-0.15 (n-octane), 3.04+/-0.09 (cyclopentane), 2.30+/-0.09 (cyclohexane), and 2.56+/-0.25 (methylcyclopentane). Measurements of the 12C/13C KIEs for the Cl atom reactions of the C2-C8 n-alkanes were also made at 348 K, and no significant temperature dependence was observed. To our knowledge, these 12C/13C KIE measurements for alkanes+Cl reactions are the first of their kind. Simultaneous to the KIE measurement, the rate constant for the reaction of each alkane with Cl atoms was measured using a relative rate method. Our measurements agree with published values within+/-20%. The measured rate constant for methylcyclopentane, for which no literature value is available, is (2.83+/-0.11)x10-10 cm3 molecule-1 s-1, 1sigma standard error. The Clepsilon values presented here for the C2-C8 alkanes are an order of magnitude smaller than reported methane Clepsilon values (Geophys. Res. Lett., 2000, 27, 1715), in contrast to reported OHepsilon values for methane (J. Geophys. Res. (Atmos.), 2001, 106, 23, 127) and C2-C8 alkanes (J. Phys. Chem. A, 2004, 108, 11537), which are all smaller than 10 per thousand. This has important implications for atmospheric modeling of saturated NMHC stable carbon isotope ratios. 13C-structure reactivity relationship values (13C-SRR) for alkane-Cl reactions have been determined and are similar to previously reported values for alkane-OH reactions.     

Linear free energy relationship and kinetic isotope effects as measures for the transition state variation. A computational study

Linear free energy relationship (LFER) and kinetic isotope effects (KIEs) are frequently used experimental means to study reaction mechanisms, in particular the nature of transition states (TSs). Density functional theory (B3LYP/6-311+G**) calculations were carried out on a model reaction, acid-catalyzed ionization of phenylethyl alcohol, to analyze how experimentally observable properties, such as nonlinearity in the Hammett and Br?nsted relations and variation in KIE, are related to a variation of the transition state structure and the mechanism. Several conclusions and insights were obtained: (1) Linear Hammett plots with a dual parameter treatment may not be evidence for an invariable TS structure for a series of reactions. (2) Variations of KIEs indeed reflect the variations of TS structures. (3) Nonlinear Br?nsted plots cannot always be taken as evidence for a stepwise mechanism. (4) A TS structure in the gas phase may change much more easily than a TS structure in solution.     

Transition-state analysis of S. pneumoniae 5'-methylthioadenosine nucleosidase

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.     

Probing the transition states of four glucoside hydrolyses with 13C kinetic isotope effects measured at natural abundance by NMR spectroscopy

Kinetic isotope effects (KIEs) were measured for methyl glucoside (4) hydrolysis on unlabeled material by NMR. Twenty-eight (13)C KIEs were measured on the acid-catalyzed hydrolysis of alpha-4 and beta-4, as well as enzymatic hydrolyses with yeast alpha-glucosidase and almond beta-glucosidase. The 1-(13)C KIEs on the acid-catalyzed reactions of alpha-4 and beta-4, 1.007(2) and 1.010(6), respectively, were in excellent agreement with the previously reported values (1.007(1), 1.011(2): Bennet and Sinnott, J. Am. Chem. Soc. 1986, 108, 7287). Transition state analysis of the acid-catalyzed reactions using the (13)C KIEs, along with the previously reported (2)H KIEs, confirmed that both reactions proceed with a stepwise D(N)A(N) mechanism and showed that the glucosyl oxocarbenium ion intermediate exists in an E(3) sofa or (4)H(3) half-chair conformation. (13)C KIEs showed that the alpha-glucosidase reaction also proceeded through a D(N)*A(N) mechanism, with a 1-(13)C KIE of 1.010(4). The secondary (13)C KIEs showed evidence of distortions in the glucosyl ring at the transition state. For the beta-glucosidase-catalyzed reaction, the 1-(13)C KIE of 1.032(1) demonstrated a concerted A(N)D(N) mechanism. The pattern of secondary (13)C KIEs was similar to the acid-catalyzed reaction, showing no signs of distortion. KIE measurement at natural abundance makes it possible to determine KIEs much more quickly than previously, both by increasing the speed of KIE measurement and by obviating the need for synthesis of isotopically labeled compounds.     

Turnover is rate-limited by deglycosylation for Micromonospora viridifaciens sialidase-catalyzed hydrolyses: conformational implications for the Michaelis complex

A panel of seven isotopically substituted sialoside natural substrate analogues based on the core structure 7-(5-acetamido-3,5-dideoxy-d-glycero-α-d-galacto-non-2-ulopyranosylonic acid)-(2→6)-β-D-galactopyranosyloxy)-8-fluoro-4-methylcoumarin (1, Neu5Acα2,6GalβFMU) have been synthesized and used to probe the rate-limiting step for turnover by the M. viridifaciens sialidase. The derived kinetic isotope effects (KIEs) on k(cat) for the ring oxygen ((18)V), leaving group oxygen ((18)V), anomeric carbon ((13)V), C3-carbon ((13)V), C3-R deuterium ((D)V(R)), C3-S deuterium ((D)V(S)), and C3-dideuterium ((D)(2)V) are 0.986 ± 0.003, 1.003 ± 0.005, 1.021 ± 0.006, 1.001 ± 0.008, 1.029 ± 0.007, 0.891 ± 0.008, and 0.890 ± 0.006, respectively. The solvent deuterium KIE ((D(2)O)V) for the sialidase-catalyzed hydrolysis of 1 is 1.585 ± 0.004. In addition, a linear proton inventory was measured for the rate of hydrolysis, under saturating condition, as a function of n, the fraction of deuterium in the solvent. These KIEs are compatible with rate-determining cleavage of the enzymatic tyrosinyl β-sialoside intermediate. Moreover, the secondary deuterium KIEs are consistent with the accumulating Michaelis complex in which the sialosyl ring of the carbohydrate substrate is in a (6)S(2) skew boat conformation. These KIE measurements are also consistent with the rate-determining deglycosylation reaction occurring via an exploded transition state in which synchronous charge delocalization is occurring onto the ring oxygen atom. Finally, the proton inventory and the magnitude of the solvent KIE are consistent with deglycosylation involving general acid-catalyzed protonation of the departing tyrosine residue rather than general base-assisted attack of the nucleophilic water.     

Determining the transition-state structure for different SN2 reactions using experimental nucleophile carbon and secondary alpha-deuterium kinetic isotope effects and theory

Nucleophile (11)C/ (14)C [ k (11)/ k (14)] and secondary alpha-deuterium [( k H/ k D) alpha] kinetic isotope effects (KIEs) were measured for the S N2 reactions between tetrabutylammonium cyanide and ethyl iodide, bromide, chloride, and tosylate in anhydrous DMSO at 20 degrees C to determine whether these isotope effects can be used to determine the structure of S N2 transition states. Interpreting the experimental KIEs in the usual fashion (i.e., that a smaller nucleophile KIE indicates the Nu-C alpha transition state bond is shorter and a smaller ( k H/ k D) alpha is found when the Nu-LG distance in the transition state is shorter) suggests that the transition state is tighter with a slightly shorter NC-C alpha bond and a much shorter C alpha-LG bond when the substrate has a poorer halogen leaving group. Theoretical calculations at the B3LYP/aug-cc-pVDZ level of theory support this conclusion. The results show that the experimental nucleophile (11)C/ (14)C KIEs can be used to determine transition-state structure in different reactions and that the usual method of interpreting these KIEs is correct. The magnitude of the experimental secondary alpha-deuterium KIE is related to the nucleophile-leaving group distance in the S N2 transition state ( R TS) for reactions with a halogen leaving group. Unfortunately, the calculated and experimental ( k H/ k D) alpha's change oppositely with leaving group ability. However, the calculated ( k H/ k D) alpha's duplicate both the trend in the KIE with leaving group ability and the magnitude of the ( k H/ k D) alpha's for the ethyl halide reactions when different scale factors are used for the high and the low energy vibrations. This suggests it is critical that different scaling factors for the low and high energy vibrations be used if one wishes to duplicate experimental ( k H/ k D) alpha's. Finally, neither the experimental nor the theoretical secondary alpha-deuterium KIEs for the ethyl tosylate reaction fit the trend found for the reactions with a halogen leaving group. This presumably is found because of the bulky (sterically hindered) leaving group in the tosylate reaction. From every prospective, the tosylate reaction is too different from the halogen reactions to be compared.     

A new interpretation of chlorine leaving group kinetic isotope effects; a theoretical approach

The chlorine leaving group kinetic isotope effects (KIEs) for the S(N)2 reactions between methyl chloride and a wide range of anionic, neutral, and radical anion nucleophiles were calculated in the gas phase and, in several cases, using a continuum solvent model. In contrast to the expected linear dependence of the chlorine KIEs on the C(alpha)-Cl bond order in the transition state, the KIEs fell in a very small range (1.0056-1.0091), even though the C(alpha)-Cl transition state bond orders varied widely from approximately 0.32 to 0.78, a range from reactant-like to very product-like. This renders chlorine KIEs, and possibly other leaving-group KIEs, less useful for studies of reaction mechanisms than commonly assumed. A partial explanation for this unexpected relationship between the C(alpha)-Cl transition state bond order and the magnitude of the chlorine KIE is presented.     

Transition state analysis of acid-catalyzed dAMP hydrolysis

Multiple kinetic isotope effects (KIEs) on deoxyadenosine monophosphate (dAMP) hydrolysis in 0.1 M HCl were used to determine the transition state (TS) structure and probe its intrinsic reactivity. The experimental KIEs revealed a stepwise (SN1) mechanism, with a discrete oxacarbenium ion intermediate. This is the first direct evidence for the deoxyribosyl oxacarbenium ion in solution. In 50% methanol/0.1 M HCl the products were deoxyribose 5-phosphate (dRMP) and alpha- and beta-methyl dRMP. The alpha-Me-dRMP/beta-Me-dRMP ratio was 8.5:1. Assuming that a free oxacarbenium ion is equally susceptible to nucleophilic attack on either face, this indicated that approximately 20% proceeded through a solvent-separated ion pair complex, or free oxacarbenium ion, a DN+AN mechanism, while approximately 80% of the reaction proceeded through a contact ion pair complex. The oxacarbenium ion lifetime was estimated at 10(-11)-10(-10) s. Computational transition states were found for ANDN, DN*AN, DN*AN, and DN+AN mechanisms using hybrid density functional theory calculations. After taking into account 20% of DN+AN, there was an excellent match of calculated to experimental KIEs for 80% of the reaction having a DN*AN mechanism. That is, C-N bond cleavage is reversible, with dAMP and the {oxacarbenium ion*adenine} complex in equilibrium. The first irreversible step is water attack on the oxacarbenium ion. The calculated 1'-14C KIE for a stepwise mechanism with irreversible C-N bond cleavage (DN*AN) was 1.052, in the range previously associated only with ANDN transition states, and close to the calculated ANDN value, 1.059. The 1'-14C KIE was strongly dependent on the adenine protonation state.     

Alpha-secondary isotope effects as probes of "tunneling-ready" configurations in enzymatic H-tunneling: insight from environmentally coupled tunneling models

Using alpha-secondary kinetic isotope effects (2 degrees KIEs) in conjunction with primary (1 degrees ) KIEs, we have investigated the mechanism of environmentally coupled hydrogen tunneling in the reductive half-reactions of two homologous flavoenzymes, morphinone reductase (MR) and pentaerythritol tetranitrate reductase (PETNR). We find exalted 2 degrees KIEs (1.17-1.18) for both enzymes, consistent with hydrogen tunneling. These 2 degrees KIEs, unlike 1 degrees KIEs, are independent of promoting motions-a nonequilibrium pre-organization of cofactor and active site residues that is required to bring the reactants into a "tunneling-ready" configuration. That these 2 degrees KIEs are identical suggests the geometries of the "tunneling-ready" configurations in both enzymes are indistinguishable, despite the fact that MR, but not PETNR, has a clearly temperature-dependent 1 degrees KIE. The work emphasizes the benefit of combining studies of 1 degrees and 2 degrees KIEs to report on pre-organization and local geometries within the context of contemporary environmentally coupled frameworks for H-tunneling.     

The mechanism of base-promoted HF elimination from 4-fluoro-4-(4-nitrophenyl)butan-2-one is E1cB. evidence from double isotopic fractionation experiments

Leaving-group fluorine and secondary deuterium multiple kinetic isotope effects (KIEs) have been determined for the base-promoted HF elimination from the 4-fluoro-4-(4'-nitrophenyl)-(1,1,1,3,3-(2)H(5))butan-2-one. The fluorine KIE was determined by using the accelerator-produced short-lived radionuclide (18)F in combination with the naturally abundant (19)F. The (19)F substrate was labeled with (14)C in a remote position to enable radioactivity measurements of both substrates. The size of the determined fluorine KIE is 1.0009 +/- 0.0010 when acetate is used as base. The secondary deuterium KIEs are 1.009 +/- 0.017, 1.000 +/- 0.018, and 1.010 +/- 0.023 for formate, acetate, and imidazole, respectively. The magnitudes of these KIEs are significantly smaller compared to the corresponding KIEs that we recently reported for the protic substrate. This new data clearly demonstrates that the elimination proceeds via an E1cB mechanism.     

Transition state structure of E. coli tRNA-specific adenosine deaminase

Bacterial tRNA-specific adenosine deaminase (TadA) catalyzes the essential deamination of adenosine to inosine at the wobble position of tRNAs and is necessary to permit a single tRNA species to recognize multiple codons. The transition state structure of Escherichia coli TadA was characterized by kinetic isotope effects (KIEs) and quantum chemical calculations. A stem loop of E. coli tRNA(Arg2) was used as a minimized TadA substrate, and its adenylate editing site was isotopically labeled as [1'-(3)H], [5'-(3)H2], [1'-(14)C], [6-(13)C], [6-(15)N], [6-(13)C, 6-(15)N] and [1-(15)N]. The intrinsic KIEs of 1.014, 1.022, 0.994, 1.014 and 0.963 were obtained for [6-(13)C]-, [6-(15)N]-, [1-(15)N]-, [1'-(3)H]-, [5'-(3)H2]-labeled substrates, respectively. The suite of KIEs are consistent with a late SNAr transition state with a complete, pro-S-face hydroxyl attack and nearly complete N1 protonation. A significant N6-C6 dissociation at the transition state of TadA is indicated by the large [6-(15)N] KIE of 1.022 and corresponds to an N6-C6 distance of 2.0 A in the transition state structure. Another remarkable feature of the E. coli TadA transition state structure is the Glu70-mediated, partial proton transfer from the hydroxyl nucleophile to the N6 leaving group. KIEs correspond to H-O and H-N distances of 2.02 and 1.60 A, respectively. The large inverse [5'-(3)H] KIE of -3.7% and modest normal [1'-(3)H] KIE of 1.4% indicate that significant ribosyl 5'-reconfiguration and purine rotation occur on the path to the transition state. The late SNAr transition-state established here for E. coli TadA is similar to the late transition state reported for cytidine deaminase. It differs from the early SNAr transition states described recently for the adenosine deaminases from human, bovine, and Plasmodium falciparum sources. The ecTadA transition state structure reveals the detailed architecture for enzymatic catalysis. This approach should be readily transferable for transition state characterization of other RNA editing enzymes.     

Heavy-atom kinetic isotope effects,cocatalysts, and the propagation transition state for polymerization of 1-hexene using the rac-(C(2)H(4)(1-indenyl)(2))ZrMe(2) catalyst precursor

Using 13C NMR techniques, the 12C/13C kinetic isotope effects (KIEs) for the polymerization of 1-hexene catalyzed by rac-(C2H4(1-indenyl)2)ZrMe2 in the presence of four different cocatalysts (tris(pentfluorophenyl)borane, tris(pentafluorophenyl)alane, anilinium tetrakis(pentafluorophenyl)borate, and methylalumoxane) have been determined. All cocatalysts yield similar KIEs within experimental uncertainty, with values of 1.009(1) and 1.017(2) at C1 and C2, respectively. Ab initio DFT computational modeling of the polymerization KIE indicates that alkene binding to the catalyst must be reversible, with the majority of the KIE developing in the subsequent migratory insertion reaction.     

Transition-state structure of human 5'-methylthioadenosine phosphorylase

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.     

Deep tunneling dominates the biologically important hydride transfer reaction from NADH to FMN in morphinone reductase

The temperature dependence of the primary kinetic isotope effect (KIE), combined temperature-pressure studies of the primary KIE, and studies of the alpha-secondary KIE previously led us to infer that hydride transfer from nicotinamide adenine dinucleotide to flavin mononucleotide in morphinone reductase proceeds via environmentally coupled hydride tunneling. We present here a computational analysis of this hydride transfer reaction using QM/MM molecular dynamics simulations and variational transition-state theory calculations. Our calculated primary and secondary KIEs are in good agreement with the corresponding experimental values. Although the experimentally observed KIE lies below the semiclassical limit, our calculations suggest that approximately 99% of the reaction proceeds via tunneling: this is the first "deep tunneling" reaction observed for hydride transfer. We also show that the dominant tunneling mechanism is controlled by the isotope at the primary rather than the secondary position: with protium in the primary position, large-curvature tunneling dominates, whereas with deuterium in this position, small-curvature tunneling dominates. Also, our study is consistent with tunneling being preceded by reorganization: in the reactant, the rings of the nicotinamide and isoalloxazine moieties are stacked roughly parallel to each other, and as the system moves toward a "tunneling-ready" configuration, the nicotinamide ring rotates to become almost perpendicular to the isoalloxazine ring.     

Transition-state variation in human, bovine, and Plasmodium falciparum adenosine deaminases

Adenosine deaminases (ADAs) from human, bovine, and Plasmodium falciparum sources were analyzed by kinetic isotope effects (KIEs) and shown to have distinct but related transition states. Human adenosine deaminase (HsADA) is present in most mammalian cells and is involved in B- and T-cell development. The ADA from Plasmodium falciparum (PfADA) is essential in this purine auxotroph, and its inhibition is expected to have therapeutic effects for malaria. Therefore, ADA is of continuing interest for inhibitor design. Stable structural mimics of ADA transition states are powerful inhibitors. Here we report the transition-state structures of PfADA, HsADA, and bovine ADA (BtADA) solved using competitive kinetic isotope effects (KIE) and density functional calculations. Adenines labeled at [6-13C], [6-15N], [6-13C, 6-15N], and [1-15N] were synthesized and enzymatically coupled with [1'-14C] ribose to give isotopically labeled adenosines as ADA substrates for KIE analysis. [6-13C], [6-15N], and [1-15N]adenosines reported intrinsic KIE values of (1.010, 1.011, 1.009), (1.005, 1.005, 1.002), and (1.004, 1.001, 0.995) for PfADA, HsADA, and BtADA, respectively. The differences in intrinsic KIEs reflect structural alterations in the transition states. The [1-15N] KIEs and computational modeling results indicate that PfADA, HsADA, and BtADA adopt early SNAr transition states, where N1 protonation is partial and the bond order to the attacking hydroxyl nucleophile is nearly complete. The key structural variation among PfADA, HsADA, and BtADA transition states lies in the degree of N1 protonation with the decreased bond lengths of 1.92, 1.55, and 1.28 A, respectively. Thus, PfADA has the earliest and BtADA has the most developed transition state. This conclusion is consistent with the 20-36-fold increase of kcat in comparing PfADA with HsADA and BtADA.     

Hydroxylation by the hydroperoxy-iron species in cytochrome P450 enzymes

Intramolecular and intermolecular kinetic isotope effects (KIEs) were determined for hydroxylation of the enantiomers of trans-2-(p-trifluoromethylphenyl)cyclopropylmethane (1) by hepatic cytochrome P450 enzymes, P450s 2B1, Delta2B4, Delta2B4 T302A, Delta2E1, and Delta2E1 T303A. Two products from oxidation of the methyl group were obtained, unrearranged trans-2-(p-trifluoromethylphenyl)cyclopropylmethanol (2) and rearranged 1-(p-trifluoromethylphenyl)but-3-en-1-ol (3). In intramolecular KIE studies with dideuteriomethyl substrates (1-d(2)) and in intermolecular KIE studies with mixtures of undeuterated (1-d(0)) and trideuteriomethyl (1-d(3)) substrates, the apparent KIE for product 2 was consistently larger than the apparent KIE for product 3 by a factor of ca. 1.2. Large intramolecular KIEs found with 1-d(2) (k(H)/k(D) = 9-11 at 10 degrees C) were shown not to be complicated by tunneling effects by variable temperature studies with two P450 enzymes. The results require two independent isotope-sensitive processes in the overall hydroxylation reactions that are either competitive or sequential. Intermolecular KIEs were partially masked in all cases and largely masked for some P450s. The intra- and intermolecular KIE results were combined to determine the relative rate constants for the unmasking and hydroxylation reactions, and a qualitative correlation was found for the unmasking reaction and release of hydrogen peroxide from four of the P450 enzymes in the absence of substrate. The results are consistent with the two-oxidants model for P450 (Vaz, A. D. N.; McGinnity, D. F.; Coon, M. J. Proc. Natl. Acad. Sci. U.S.A. 1998, 95, 3555), which postulates that a hydroperoxy-iron species (or a protonated analogue of this species) is a viable electrophilic oxidant in addition to the consensus oxidant, iron-oxo.