Effect of deuterium substitution on electron transfer at cytochrome c/SAM interfaces

The involvement of protons in the heterogeneous electron transfer between cytochrome c and a gold electrode to which it is attached was studied by comparing the electron transfer rate constants for H2O and D2O solutions. Rate constants were measured as a function of the electrochemical cell solution and the protein incubant solution, i.e., k (0)(incubant, cell). Two separate isotope effects exist: a cell "isotope effect", KIE cell = k (0)(H2O, H2O): k (0)(H2O, D2O), which is manifest at short time scales (<30 s) and arises from the viscosity difference between H2O and D2O, and an incubant isotope effect, KIE inc= k (0)(H2O, H2O): k (0)(D2O, H2O), which is manifest at longer times (>2 h) and results from H/D exchange. The two isotope effects are approximately equal ( approximately 1.2) and a total isotope effect KIE total = k (0)(H2O, H2O): k (0)(D2O, D2O) can be constructed that is the product of KIE cell and KIE inc. The nature of the electron transfer process, possible coupling to a proton transfer process, and the involvement of specific hydrogens in the transfer mechanism are discussed.     

Calculations find that tunneling plays a major role in the reductive elimination of methane from hydridomethylbis(trimethylphosphine)platinum: how to confirm this computational prediction experimentally

DFT calculations, including the effects of small curvature tunneling, have been performed on the reductive elimination of methane from hydridomethylbis(trimethylphosphine)platinum (1d). The calculations find that at 250 K tunneling results in an increase in the rate constant for reductive elimination by a factor of 4, a lowering of Ea by 1.7 kcal/mol, and a decrease in A by a factor of nearly 10. Tunneling is also calculated to increase the primary H/D kinetic isotope effect (KIE) from k(1d)/k(1f) = 2.26 to k(1d)/k(1f) = 4.12 and to result in a large secondary KIE of k(1d)/k(1e) = 1.35. In addition, tunneling is predicted to result in a violation of the rule of the geometric mean, so that the secondary KIE for reductive elimination of methane-d1 from 1f is calculated to be k(1f)/k(1g) = 1.06, which is much smaller than the secondary KIE of k(1d)/k(1e) = 1.35 for reductive elimination of methane from 1d. Comparison of the measured values of k(1d)/k(1e) and k(1f)/k(1g) is therefore proposed as an experimental test of the prediction that tunneling plays an important role in the reductive elimination of methane from 1d.     

Hydride reduction of NAD+ analogues by isopropyl alcohol: kinetics, deuterium isotope effects and mechanism

Observed pseudo-first-order rate constants (k(obs)) of the hydride-transfer reactions from isopropyl alcohol (i-PrOH) to two NAD(+) analogues, 9-phenylxanthylium ion (PhXn(+)) and 10-methylacridinium ion (MA(+)), were determined at temperatures ranging from 49 to 82 degrees C in i-PrOH containing various amounts of AN or water. Formations of the alcohol-cation ether adducts (ROPr-i) were observed as side equilibria. The equilibrium constants for the conversion of PhXn(+) to PhXnOPr-i in i-PrOH/AN (v/v = 1) were determined, and the equilibrium isotope effect (EIE = K(i-PrOH)/K(i-PrOD)) at 62 degrees C was calculated to be 2.67. The k(H) of the hydride-transfer step for both reactions were calculated on the basis of the k(obs) and K. The corresponding deuterium kinetic isotope effects (e.g., KIE(OD)(H) = k(H)(i-PrOH)/k(H)(i-PrOD) and KIE(beta-D6)(H) = k(obs)(i-PrOH)/k(obs)((CD3)2CHOH)), as well as the activation parameters, were derived. For the reaction of PhXn(+) (62 degrees C) and MA(+) (67 degrees C), primary KIE(alpha-D)(H) (4.4 and 2.1, respectively) as well as secondary KIE(OD)(H) (1.07 and 1.18) and KIE(beta-D6)(H) (1.1 and 1.5) were observed. The observed EIE and KIE(OD)(H) were explained in terms of the fractionation factors for deuterium between OH and OH(+)(OH(delta+)) sites. The observed inverse kinetic solvent isotope effect for the reaction of PhXn(+) (k(obs)(i-PrOH)/k(obs)(i-PrOD) = 0.39) is consistent with the intermolecular hydride-transfer mechanism. The dramatic reduction of the reaction rate for MA(+), when the water or i-PrOH cosolvent was replaced by AN, suggests that the hydride-transfer T.S. is stabilized by H-bonding between O of the solvent OH and the substrate alcohol OH(delta+). This result suggests an H-bonding stabilization effect on the T.S. of the alcohol dehydrogenase reactions.     

Ensemble-averaged QM/MM kinetic isotope effects for the S(N)2 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° α-(2)H(2), 2° β-(2)H(3), α-(11)C/(14)C, leaving group (37)Cl, and nucleophile (13)C and (15)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.     

Co(+)-assisted decomposition of h6-acetone and d6-acetone: acquisition of reaction rate constants and dynamics of the dissociative mechanism

Reaction rate constants have been acquired for the gaseous unimolecular decomposition reaction of the Co(+)(OC(CH(3))(2)) cluster ion and its deuterium labeled analog. Each rate constant is measured at a well resolved cluster internal energy within the range 12,300-16,100 cm(-1). The weighted, averaged kinetic isotope effect (KIE), k(H)/k(D) = 1.54 ± 0.05, is about three times smaller than the KIE measured for the rate-determining rate constants in the similar Ni(+)(OC(CH(3))(2)) decomposition reaction. These reactions likely follow the same oxidative addition-reductive elimination mechanism. Thus, this unexpected change in the KIE magnitudes is not due to differences in the dissociative reaction coordinates. Rather, we propose that the unique dissociation dynamics of these two similar systems is due to differences in the low-lying electronic structure of each transition metal ion.     

Gauging the relative oxidative powers of compound I, ferric-hydroperoxide, and the ferric-hydrogen peroxide species of cytochrome P450 toward C-H hydroxylation of a radical clock substrate

Density functional calculations were performed in response to the controversies regarding the identity of the oxidant species in cytochrome P450. The calculations were used to gauge the relative C-H hydroxylation reactivity of three potential oxidant species of the enzyme, the high-valent oxo-iron species Compound I (Cpd I), the ferric hydroperoxide Compound 0 (Cpd 0), and the ferric-hydrogen peroxide complex Fe(H(2)O(2)). The results for the hydroxylation of a radical probe substrate, 1, show the following trends: (a) Cpd I is the most reactive species; in its presence the other two reagents will be silent. (b) In the absence of Cpd I, substrate oxidation by Cpd 0 and Fe(H(2)O(2)) will take place via a stepwise mechanism that involves initial O-O homolysis followed by H-abstraction from 1. (c) Cpd 0 will undergo mostly porphyrin hydroxylation and only approximately 15% of substrate oxidation producing mostly the rearranged alcohol, 3 (Scheme 2). (d) Fe(H(2)O(2)) will generate mostly free hydrogen peroxide (uncoupling). A small fraction will perform substrate oxidation and lead mostly to 3. Reactivity probes for these reagents are kinetic isotope effect (KIE) and the product ratio of unrearranged to rearranged alcohols, [2/3]. Thus, for substrate oxidation by Cpd 0 or Fe(H(2)O(2)) KIE will be small, approximately 2, while Cpd I will have large KIE values. Typically both Cpd 0 and Fe(H(2)O(2)) will lead to a [2/3] ratio < 1, while Cpd I will lead to ratios > 1. In addition, the product isotope effect (KIE(2)/KIE(3) not equal 1) is expected from the reactivity of Cpd I.     

Ultrafast excited-state intermolecular proton transfer of cyanine fluorochrome dyes

Steady-state and time-resolved emission spectroscopy techniques were employed to study the excited-state proton transfer (ESPT) to water and D(2)O from QCy7, a recently synthesized near-infrared (NIR)-emissive dye with a fluorescence band maximum at 700 nm. We found that the ESPT rate constant, k(PT), of QCy7 excited from its protonated form, ROH, is ~1.5 × 10(12) s(-1). This is the highest ever reported value in the literature thus far, and it is comparable to the reciprocal of the longest solvation dynamics time component in water, τ(S) = 0.8 ps. We found a kinetic isotope effect (KIE) on the ESPT rate of ~1.7. This value is lower than that of weaker photoacids, which usually have KIE value of ~3, but comparable to the KIE on proton diffusion in water of ~1.45, for which the average time of proton transfer between adjacent water molecules is similar to that of QCy7.     

Reactions of PhSCH(2)Li and NCCH(2)Li with benzaldehyde and benzophenone: when does the mechanism change from ET to polar?

The carbonyl-carbon kinetic isotope effect (KIE) and the substituent effect were measured for the reaction of phenylthiomethyllithium (PhSCH(2)Li, 1) with benzaldehyde and benzophenone, and cyanomethyllithium (NCCH(2)Li, 2) with benzaldehyde, and the results were compared with those for other lithium reagents such as MeLi, PhLi, CH(2)=CHCH(2)Li, and CH(2)=C(OLi)C(CH(3))(3). It was previously shown that the reactions of MeLi, PhLi, and CH(2)=CHCH(2)Li proceed via a rate-determining electron transfer (ET) process whereas the reaction of lithium pinacolone enolate goes through the polar (PL) mechanism. The reaction of 1 with benzaldehyde gave no carbonyl-carbon KIE ((12)k/(13)k = 0.999 +/- 0.004), similar to that measured previously for the MeLi reaction with benzophenone ((12)k/(14)k = 1.000). The effect of substituents of the aromatic ring of benzaldehyde and benzophenone on the reactivity gave very small Hammett rho values of 0.17 +/- 0.03 and 0.26 +/- 0.05, respectively. These small rho values are again similar to that observed for the reaction of MeLi. Likewise the reactions of 2 with benzaldehydes gave small KIE and the rho value ((12)k/(13)k = 0.996 +/- 0.004, rho = 0.14 +/- 0.02). Dehalogenation and enone-isomerization probe experiments for 2 showed no evidence for the presence of radical-ion pair of sufficient lifetime during the course of the reaction. It is concluded that the reactions of 1 and 2 with the aromatic carbonyl compounds proceed via the electron transfer-radical coupling mechanism with rate-determining ET as in the reactions of MeLi, PhLi, and CH(2)=CHCH(2)Li.     

A new insight into using chlorine leaving group and nucleophile carbon kinetic isotope effects to determine substituent effects on the structure of SN2 transition states

Chlorine leaving group k(35)/k(37), nucleophile carbon k(11)/k(14), and secondary alpha-deuterium [(kH/kD)alpha] kinetic isotope effects (KIEs) have been measured for the SN2 reactions between para-substituted benzyl chlorides and tetrabutylammonium cyanide in tetrahydrofuran at 20 degrees C to determine whether these isotope effects can be used to determine the substituent effect on the structure of the transition state. The secondary alpha-deuterium KIEs indicate that the transition states for these reactions are unsymmetric. The theoretical calculations at the B3LYP/aug-cc-pVDZ level of theory support this conclusion; i.e., they suggest that the transition states for these reactions are unsymmetric with a long NC-C(alpha) and reasonably short C(alpha)-Cl bonds. The chlorine isotope effects suggest that these KIEs can be used to determine the substituent effects on transition state structure with the KIE decreasing when a more electron-withdrawing para-substituent is present. This conclusion is supported by theoretical calculations. The nucleophile carbon k(11)/k(14) KIEs for these reactions, however, do not change significantly with substituent and, therefore, do not appear to be useful for determining how the NC-C(alpha) transition-state bond changes with substituent. The theoretical calculations indicate that the NC-C(alpha) bond also shortens as a more electron-withdrawing substituent is placed on the benzene ring of the substrate but that the changes in the NC-C(alpha) transition-state bond with substituent are very small and may not be measurable. The results also show that using leaving group and nucleophile carbon KIEs to determine the substituent effect on transition-state structure is more complicated than previously thought. The implication of using both chlorine leaving group and nucleophile carbon KIEs to determine the substituent effect on transition-state structure is discussed.     

Multidimensional tunneling in the [1,5] shift in (Z)-1,3-pentadiene: how useful are Swain-Schaad exponents at detecting tunneling?

Rates, kinetic isotope effects (KIE), and Swain-Schaad exponents (SSE) have been calculated for a variety of isotopologues for the [1,5] shift in (Z)-1,3-pentadiene using mPW1K/6-31+G(d,p). Quantum mechanical effects along the reaction coordinate were incorporated with the zero-curvature tunneling (ZCT) model and with the multidimensional small curvature tunneling (SCT) model, which allows for coupling of modes perpendicular to the reaction coordinate. The latter model gives the best agreement with experimental rates and primary KIEs. The small quasiclassical primary KIE (2.6) is rationalized in terms of a nonlinear transition state. For sp3 to sp2 rehybridization, the quasiclassical alpha-secondary KIE shows an unusual inverse effect due to compression of the nonbonding hydrogens in the suprafacial transition state. SCT transmission coefficients (kappa) increase the rates by as much as one order of magnitude. Tunneling allows the reactant to evade 1-2.5 kcal/mol of the barrier depending on the isotope. Inclusion of tunneling in the secondary KIE increases it beyond the equilibrium isotope effect and converts the inverse effect (0.95) into a normal KIE (1.12). Tunneling was found to deflate the primary y SSE but by an amount too small to distinguish it from the quasiclassical SSE. On the other hand, when a specific labeling pattern is used, the difference between the quasiclassical secondary SSE (4.1) and the tunneling secondary SSE (2.3) may be sufficiently large to detect tunneling. The mixed secondary SSE shows even larger differences.     

Mechanisms of electron transfer in catalysis by copper zinc superoxide dismutase

Activated oxygen intermediates during copper zinc superoxide dismutase (SOD) catalysis were investigated using an isotope fractionation technique and natural abundance reagents. Competitive oxygen kinetic isotope effects (KIEs) are reported for the enzyme-catalyzed disproportionation of superoxide as well as the stoichiometric reaction of reduced SOD with molecular oxygen. Analysis within the context of quantum mechanical electron transfer theory provides evidence against an outer-sphere mechanism for O2*- oxidation. A CuII-O2-I intermediate is, therefore, proposed. The SOD-catalyzed oxidation of O2*- is characterized by an inverse (<1) KIE which is similar to those determined for the analogous reactions of synthetic copper compounds. An inverse kinetic isotope effect upon the enzymatic reduction of O2*- is also observed and proposed to arise from rate-determining proton transfer which leads to the formation of HO2* in the SOD active site.     

Measuring the kinetic isotope effect at natural isotopic abundances for discriminating between the homogeneous and heterogeneous catalytic mechanisms in the Heck and Suzuki reactions

The kinetic isotope effect (KIE) at the natural abundances of bromine and carbon isotopes in the substrates of the Heck and Suzuki reactions have been investigated to determine the true nature of catalyst in these reactions. Data processing has demonstrated that statistically significant differences between KIE values for the Suzuki reaction of nonactivated bromobenzene are observed upon the replacement of the soluble catalyst precursor with the insoluble one. This finding unambiguously indicates that the reaction takes place on heterogeneous palladium species. Similar experiments on the Heck reaction have demonstrated that the KIE values are insensitive to the nature of the catalyst precursor, which is consistent with the true homogeneous mechanism of catalysis.     

Secondary kinetic isotope effect on the photoenolization of triplet o-methylanthrones. A microcanonical transition state theory calculation

The energy profile for the tautomerization reaction of 1,4-dimethylanthrone in the first triplet electronic state obtained through electronic calculations (B3LYP/ 6-31G(d)) is used to calculate the rate constants for the process at a wide range of energies using a modified RRKM microcanonical statistical formalism that takes into account tunneling. Through partial or total substitution of the hydrogen atoms of the methyl groups by deuterium atoms, it is possible to evaluate different primary and secondary kinetic isotope effects (KIE). These results can be compared with experimental data for these processes taking place in solid matrix at extremely low temperatures (4-50 K). Such a comparison allows us to conclude that the reaction is taking place at energies just slightly below (around 0.5 kcal/mol) the adiabatic potential energy barrier, a result that was previously found for other related molecules so that this mechanism may be extended to the photoenolization of other o-aryl methyl ketones. Analysis of the different factors contributing to the primary and secondary KIEs discloses that at energies not far below the adiabatic barrier, the tunneling effect is not the only factor that accounts for the large KIE but the differences in the energy level distribution upon isotopic substitution may be the predominant factor at a certain range of negative energies (this is especially so for the case of primary KIE). At positive energies (above the barrier) the levels factor is always the dominant factor in the total KIE.     

Kinetic isotope effect and extended pH-dependent studies of the oxidation of hydroxylamine by hexachloroiridate(IV): ET and PCET

The oxidation of hydroxylamine by [IrCl6]2- has been studied spectrophotometrically in deoxygenated aqueous solutions in the range of pH 4-9 at 25 degrees C. The reaction is catalyzed by Cu2+, Fe2+, and impurities of aquochloroiridium complexes. Oxalate is a very effective inhibitor of catalysis by copper and iron ions. With excess hydroxylamine, the reaction follows pseudo-first-order kinetics, and the stoichiometric ratio (DeltanIr(IV)/Deltanhydroxylamine) is 1.05 at pH 5.9. Over the pH range 4.2-8.8, the empirical rate law is -d[IrCl(6)2-]/dt=k[IrCl6(2-)][NH2OH]tot, with k=k1Ka1/([H+]+Ka1)+k'Ka1/([H+]([H+]+Ka1)), where Ka1 is the dissociation constant of NH3OH+. Least-squares fitting yields k1=(17.05+/-0.47) M-1 s(-1) and k'=(2.59+/-0.09)x10(-6) s(-1) at ionic strength of 0.1 M (adjusted by NaClO4) and 25 degrees C. The kinetic isotope effects (KIE) (kH/kD) for k1 and k' are 4.4 and 9.8, correspondingly. A mechanism is inferred in which k1 corresponds to concerted proton-coupled electron transfer (PCET) and k' corresponds to electron transfer from NH2O-. In this mechanism, the large KIE for k' is due almost entirely to the equilibrium isotope effect for the pKa of NH2OH.     

Deuterium kinetic isotope effects in gas-phase S(N)2 and E2 reactions: comparison of experiment and theory

The competition between bimolecular nucleophilc substitution and base-induced elimination is investigated through kinetic isotope effect measurements for gas-phase reactions of RCl + ClO- (R = methyl, ethyl, isopropyl, and tert-butyl) utilizing a FA-SIFT instrument. The overall reaction rate constants and the kinetic isotope effect for the reaction of C2H5Cl + ClO- are compared to computational results. [Hu, W. P.; Truhlar, D. G. J. Am. Chem. Soc. 1996, 118, 860.] Experimental results show that as the degree of substitution in the neutral reactant increases the E2 channel becomes dominant. The systematic change in the overall kinetic isotope effects indicates that, for the reaction of ClO- with C2H5Cl, both the SN2 and E2 pathways do occur, as predicted by computation; however the experimental reaction rate constants and KIE deviate strongly from the computational result.     

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.     

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.     

Stable hydrogen and oxygen isotope ratios of bottled waters of the world

Bottled and packaged waters are an increasingly significant component of the human diet. These products are regulated at the regional, national, and international levels, and determining the authenticity of marketing and labeling claims represents a challenge to regulatory agencies. Here, we present a dataset of stable isotope ratios for bottled waters sampled worldwide, and consider potential applications of such data for regulatory, forensic and geochemical standardization applications. The hydrogen and oxygen isotope ratios of 234 samples of bottled water range from -147 per thousand to +15 per thousand and from -19.1 per thousand to +3.0 per thousand, respectively. These values fall within and span most of the normal range for meteoric waters, indicating that these commercially available products represent a source of waters for use as laboratory working standards in applications requiring standardization over a large range of isotope ratios. The measured values of bottled water samples cluster along the global meteoric water line, suggesting that bottled water isotope ratios preserve information about the water sources from which they were derived. Using the dataset, we demonstrate how bottled water isotope ratios provide evidence for substantial evaporative enrichment of water sources prior to bottling and for the marketing of waters derived from mountain and lowland sources under the same name. Comparison of bottled water isotope ratios with natural environmental water isotope ratios demonstrates that on average the isotopic composition of bottled water tends to be similar to the composition of naturally available local water sources, suggesting that in many cases bottled water need not be considered as an isotopically distinct component of the human diet. Our findings suggest that stable isotope ratios of bottled water have the power to distinguish ultimate (e.g., recharge) and proximal (e.g., reservoir) sources of bottled water and constitute a potential tool for use in the regulatory monitoring of water products.     

Modeling temperature dependent kinetic isotope effects for hydrogen transfer in a series of soybean lipoxygenase mutants: The effect of anharmonicity upon transfer distance

Soybean lipoxygenase-1 (SLO) catalyzes the oxidation of linoleic acid. The rate-limiting step in this transformation is the net abstraction of the pro-S hydrogen atom from the center of the 1,5-pentadienyl moiety in linoleic acid. The large deuterium kinetic isotope effect (KIE) for this step appears in the first order rate constant ((D)k(cat) = 81 ± 5 at T = 25 °C). Furthermore, the KIE and the rate for protium abstraction are weakly temperature dependent (E(A,D) - E(A,H) = 0.9 ± 0.2 kcal/mol and E(A,H) = 2.1 ± 0.2 kcal/mol, respectively). Mutations at a hydrophobic site about 13 ? from the active site Fe(III), Ile(553), induce a marked temperature dependence that varies roughly in accordance with the degree to which the residue is changed in bulk from the wild type Ile. While the temperature dependence for these mutants varies from the wild type enzyme, the magnitude of the KIE at 25 °C is on the same order of magnitude. A hydrogen tunneling model [Kuznetsov, A.M., Ulstrup, J. Can. J. Chem. 77 (1999) 1085-1096] is utilized to model the KIE temperature profiles for the wild type SLO and each Ile(553) mutant. Hydrogenic wavefunctions are modeled using harmonic oscillators and Morse oscillators in order to explore the effects of anharmonicity upon computed kinetic observables used to characterize this hydrogen transfer.