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.     

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.     

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.     

13C kinetic isotope effects in the copper(I)-mediated living radical polymerization of methyl methacrylate

Carbon-13 kinetic isotope effects (KIEs) have been determined for free-radical and copper-mediated living radical polymerizations of methyl methacrylate at 60 degrees C. While free-radical polymerization shows only one primary 13C KIE, on the least-substituted double bond carbon (k12/k13 = 1.045), two significant KIEs are observed, one on each double bond carbon, for copper-mediated polymerization (k12/k13(H2C=) = 1.050, k12/k13(=C <) = 1.010), showing that copper-mediated living radical polymerization does not propagate via a simple free radical process.     

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.     

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.     

A theoretical investigation of alpha-carbon kinetic isotope effects and their relationship to the transition-state structure of S(N)2 reactions

[reaction: see text] The transition structures and alpha-carbon 12C/13C kinetic isotope effects for 22 S(N)2 reactions between methyl chloride and a wide variety of nucleophiles have been calculated using the B1LYP/aug-cc-pVDZ level of theory. Anionic, neutral, and radical anion nucleophiles were used to give a wide range of S(N)2 transition states so the relationship between the magnitude of the alpha-carbon kinetic isotope effect and transition-state structure could be determined. The results suggest that the alpha-carbon 12C/13C kinetic isotope effects for S(N)2 reactions will be large (near the experimental maximum) and that the curve relating the magnitude of the KIE to the percent transfer of the alpha-carbon from the nucleophile to the leaving group in the transition state has a broad maximum. This means very similar KIEs will be found for early, symmetric, and late transition states and that one cannot use the magnitude of these KIEs to estimate transition-state structure.     

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.     

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.     

A kinetic isotope effect study on the hydrolysis reactions of methyl xylopyranosides and methyl 5-thioxylopyranosides: oxygen versus sulfur stabilization of carbenium ions.

The following kinetic isotope effects, KIEs (k(light)/k(heavy)), have been measured for the hydrolyses of methyl alpha- and beta-xylopyranosides, respectively, in aqueous HClO(4) (mu = 1.0 M, NaClO(4)) at 80 degrees C: alpha-D, 1.128 +/- 0.004, 1.098 +/- 0.005; beta-D, 1.088 +/- 0.008, 1.042 +/- 0.004; gamma-D(2), (C5) 0.986 +/- 0.001, 0.967 +/- 0.003; leaving-group (18)O, 1.023 +/- 0.002, 1.023 +/- 0.003; ring (18)O, 0.983 +/- 0.001, 0.978 +/- 0.001; anomeric (13)C, 1.006 +/- 0.001, 1.006 +/- 0.003; and solvent, 0.434 +/- 0.017, 0.446 +/- 0.012. In conjunction with the reported (J. Am. Chem. Soc. 1986, 108, 7287-7294) KIEs for the acid-catalyzed hydrolysis of methyl alpha- and beta-glucopyranosides, it is possible to conclude that at the transition state for xylopyranoside hydrolysis resonance stabilization of the developing carbenium ion by the ring oxygen atom is coupled to exocyclic C-O bond cleavage, and the corresponding methyl glucopyranosides hydrolyze via transition states in which charge delocalization lags behind aglycon departure. In the analogous hydrolysis reactions of methyl 5-thioxylopyranosides, the measured KIEs in aqueous HClO(4) (mu = 1.0 M, NaClO(4)) at 80 degrees C for the alpha- and beta-anomers were, respectively, alpha-D, 1.142 +/- 0.010, 1.094 +/- 0.002; beta-D 1.061 +/- 0.003, 1.018(5) +/- 0.001; gamma-D(2), (C5) 0.999 +/- 0.001, 0.986 +/- 0.002; leaving-group (18)O, 1.027 +/- 0.001, 1.035 +/- 0.001; anomeric (13)C, 1.031 +/- 0.002, 1.028 +/- 0.002; solvent, 0.423 +/- 0.015, 0.380 +/- 0.014. The acid-catalyzed hydrolyses of methyl 5-thio-alpha- and beta-xylopyranosides, which occur faster than methyl alpha- and beta-xylopyranosides by factors of 13.6 and 18.5, respectively, proceed via reversibly formed O-protonated conjugate acids that undergo slow, rate-determining exocyclic C-O bond cleavage. These hydrolysis reactions do not have a nucleophilic solvent component as a feature of the thiacarbenium ion-like transition states.     

NMR study of C-kinetic isotope effects at C natural abundance to characterize oxidations and an enzyme-catalyzed reduction

¹³C-kinetic isotope effects (KIEs) of four cinnamyl alcohol oxidations and a xylose reductase-catalyzed cinnamyl aldehyde reduction have been determined by ¹³C NMR using competition reactions with reactants at natural ¹³C-abundance. Differences in KIEs among oxidations indicate dissimilarities between the respective hydrogen transfers. Their mechanistic implications are discussed. A low primary KIE of the enzymatic reduction is consistent with a kinetically complex mechanism in which steps other than the chemical step of hydride transfer from NADH are slow.     

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.     

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.     

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.     

Mechanistic Insights into RNA Transphosphorylation from Kinetic Isotope Effects and Linear Free Energy Relationships of Model Reactions

Phosphoryl transfer reactions are ubiquitous in biology and the understanding of the mechanisms whereby these reactions are catalyzed by protein and RNA enzymes is central to reveal design principles for new therapeutics. Two of the most powerful experimental probes of chemical mechanism involve the analysis of linear free energy relations (LFERs) and the measurement of kinetic isotope effects (KIEs). These experimental data report directly on differences in bonding between the ground state and the rate‐controlling transition state, which is the most critical point along the reaction free energy pathway. However, interpretation of LFER and KIE data in terms of transition‐state structure and bonding optimally requires the use of theoretical models. In this work, we apply density‐functional calculations to determine KIEs for a series of phosphoryl transfer reactions of direct relevance to the 2′‐O‐transphosphorylation that leads to cleavage of the phosphodiester backbone of RNA. We first examine a well‐studied series of phosphate and phosphorothioate mono‐, di‐ and triesters that are useful as mechanistic probes and for which KIEs have been measured. Close agreement is demonstrated between the calculated and measured KIEs, establishing the reliability of our quantum model calculations. Next, we examine a series of RNA transesterification model reactions with a wide range of leaving groups in order to provide a direct connection between observed Brønsted coefficients and KIEs with the structure and bonding in the transition state. These relations can be used for prediction or to aid in the interpretation of experimental data for similar non‐enzymatic and enzymatic reactions. Finally, we apply these relations to RNA phosphoryl transfer catalyzed by ribonuclease A, and demonstrate the reaction coordinate–KIE correlation is reasonably preserved. A prediction of the secondary deuterium KIE in this reaction is also provided. These results demonstrate the utility of building up knowledge of mechanism through the systematic study of model systems to provide insight into more complex biological systems such as phosphoryl transfer enzymes and ribozymes.     

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.     

Intramolecular kinetic isotope effect in hydride transfer from dihydroacridine to a quinolinium ion. Rejection of a proposed two-step mechanism with a kinetically significant intermediate

The intramolecular kinetic isotope effect (KIE) for hydride transfer from 10-methyl-9,10-dihydroacridine to 1-benzyl-3-cyanoquinolinium ion has been found to be 5-6 by both (1)H NMR and mass spectrometry. This KIE is consistent with other hydride transfers. It is inconsistent with the high intermolecular KIEs derived by fitting to a two-step mechanism with a kinetically significant intermediate complex, and it is inconsistent with the strong temperature dependence of those KIEs. We therefore reject the two-step mechanism for this reaction, and we suggest that other cases proposed to follow this mechanism are in error.     

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.