Intrinsic isotope effects on benzylic hydroxylation by the aromatic amino acid hydroxylases: evidence for hydrogen tunneling, coupled motion, and similar reactivities

Biochemical and biophysical approaches to characterization of the aromatic amino acid hydroxylases

The aromatic amino acid hydroxylases phenylalanine hydroxylase, tyrosine hydroxylase, and tryptophan hydroxylase utilize a non-heme iron to catalyze the hydroxylation of the aromatic rings of their amino acid substrates, with a tetrahydropterin serving as the source of the electrons necessary for the monooxygenation reaction. These enzymes have been subjected to a variety of biochemical and biophysical approaches, resulting in a detailed understanding of their structures and mechanism. We summarize here the experimental approaches that have led to this understanding.

Intermolecular Interactions and Spectroscopic Signatures of the Hydrogen-Bonded System-n-Octanol in Experimental and Theoretical Studies

n-Octanol is the object of experimental and theoretical study of spectroscopic signatures and intermolecular interactions. The FTIR measurements were carried out at 293 K for n-octanol and its deuterated form. Special attention was paid to the vibrational features associated with the O-H stretching and the isotope effect. Density Functional Theory (DFT) in its classical formulations was applied to develop static models describing intermolecular hydrogen bond (HB) and isotope effect in the gas phase and using solvent reaction field reproduced by Polarizable Continuum Model (PCM). The Atoms in Molecules (AIM) theory enabled electronic structure and molecular topology study. The Symmetry-Adapted Perturbation Theory (SAPT) was used for energy decomposition in the dimers of n-octanol. Finally, time-evolution methods, namely classical molecular dynamics (MD) and Car-Parrinello Molecular Dynamics (CPMD) were employed to shed light onto dynamical nature of liquid n-octanol with emphasis put on metric and vibrational features. As a reference, CPMD gas phase results were applied. Nuclear quantum effects were included using Path Integral Molecular Dynamics (PIMD) and a posteriori method by solving vibrational Schrödinger equation. The latter applied procedure allowed to study the deuterium isotope effect.

Dioxygen activation by nonheme iron enzymes with the 2-His-1-carboxylate facial triad that generate high-valent oxoiron oxidants

The 2-His-1-carboxylate facial triad is a widely used scaffold to bind the iron center in mononuclear nonheme iron enzymes for activating dioxygen in a variety of oxidative transformations of metabolic significance. Since the 1990s, over a hundred different iron enzymes have been identified to use this platform. This structural motif consists of two histidines and the side chain carboxylate of an aspartate or a glutamate arranged in a facial array that binds iron(II) at the active site. This triad occupies one face of an iron-centered octahedron and makes the opposite face available for the coordination of O₂ and, in many cases, substrate, allowing the tailoring of the iron-dioxygen chemistry to carry out a plethora of diverse reactions. Activated dioxygen-derived species involved in the enzyme mechanisms include iron(III)-superoxo, iron(III)-peroxo, and high-valent iron(IV)-oxo intermediates. In this article, we highlight the major crystallographic, spectroscopic, and mechanistic advances of the past 20 years that have significantly enhanced our understanding of the mechanisms of O₂ activation and the key roles played by iron-based oxidants.

A common late-stage intermediate in catalysis by 2-hydroxyethyl-phosphonate dioxygenase and methylphosphonate synthase

2-Hydroxyethylphosphonate dioxygenase (HEPD) and methylphosphonate synthase (MPnS) are nonheme iron oxygenases that both catalyze the carbon-carbon bond cleavage of 2-hydroxyethylphosphonate but generate different products. Substrate labeling experiments led to a mechanistic hypothesis in which the fate of a common intermediate determined product identity. We report here the generation of a bifunctional mutant of HEPD (E176H) that exhibits the activity of both HEPD and MPnS. The product distribution of the mutant is sensitive to a substrate isotope effect, consistent with an isotope-sensitive branching mechanism involving a common intermediate. The X-ray structure of the mutant was determined and suggested that the introduced histidine does not coordinate the active site metal, unlike the iron-binding glutamate it replaced.

Factors affecting hydrogen-tunneling contribution in hydroxylation reactions promoted by oxoiron(IV) porphyrin π-cation radical complexes

Hydrogen atom transfer with a tunneling effect (H-tunneling) has been proposed to be involved in aliphatic hydroxylation reactions catalyzed by cytochrome P450 and synthetic heme complexes as a result of the observation of large hydrogen/deuterium kinetic isotope effects (KIEs). In the present work, we investigate the factors controlling the H-tunneling contribution to the H-transfer process in hydroxylation reaction by examining the kinetics of hydroxylation reactions at the benzylic positions of xanthene and 1,2,3,4-tetrahydronaphthalene by oxoiron(IV) 5,10,15,20-tetramesitylporphyrin π-cation radical complexes ((TMP(+•))Fe(IV)O(L)) under single-turnover conditions. The Arrhenius plots for these hydroxylation reactions of H-isotopomers have upwardly concave profiles. The Arrhenius plots of D-isotopomers, clear isosbestic points, and product analysis rule out the participation of thermally dependent other reaction processes in the concave profiles. These results provide evidence for the involvement of H-tunneling in the rate-limiting H-transfer process. These profiles are simulated using an equation derived from Bell's tunneling model. The temperature dependence of the KIE values (k(H)/k(D)) determined for these reactions indicates that the KIE value increases as the reaction temperature becomes lower, the bond dissociation energy (BDE) of the C-H bond of a substrate becomes higher, and the reactivity of (TMP(+•))Fe(IV)O(L) decreases. In addition, we found correlation of the slope of the ln(k(H)/k(D)) - 1/T plot and the bond strengths of the Fe═O bond of (TMP(+•))Fe(IV)O(L) estimated from resonance Raman spectroscopy. These observations indicate that these factors modulate the extent of the H-tunneling contribution by modulating the ratio of the height and thickness of the reaction barrier.

Pulsed EPR study of amino acid and tetrahydropterin binding in a tyrosine hydroxylase nitric oxide complex: evidence for substrate rearrangements in the formation of the oxygen-reactive complex

Tyrosine hydroxylase is a nonheme iron enzyme found in the nervous system that catalyzes the hydroxylation of tyrosine to form l-3,4-dihydroxyphenylalanine, the rate-limiting step in the biosynthesis of the catecholamine neurotransmitters. Catalysis requires the binding of three substrates: tyrosine, tetrahydrobiopterin, and molecular oxygen. We have used nitric oxide as an O₂ surrogate to poise Fe(II) at the catalytic site in an S = 3/2, {FeNO}⁷ form amenable to EPR spectroscopy. ²H-electron spin echo envelope modulation was then used to measure the distance and orientation of specifically deuterated substrate tyrosine and cofactor 6-methyltetrahydropterin with respect to the magnetic axes of the {FeNO}⁷ paramagnetic center. Our results show that the addition of tyrosine triggers a conformational change in the enzyme that reduces the distance from the {FeNO}⁷ center to the closest deuteron on 6,7-²H-6-methyltetrahydropterin from >5.9 Å to 4.4 ± 0.2 Å. Conversely, the addition of 6-methyltetrahydropterin to enzyme samples treated with 3,5-²H-tyrosine resulted in reorientation of the magnetic axes of the S = 3/2, {FeNO}⁷ center with respect to the deuterated substrate. Taken together, these results show that the coordination of both substrate and cofactor direct the coordination of NO to Fe(II) at the active site. Parallel studies of a quaternary complex of an uncoupled tyrosine hydroxylase variant, E332A, show no change in the hyperfine coupling to substrate tyrosine and cofactor 6-methyltetrahydropterin. Our results are discussed in the context of previous spectroscopic and X-ray crystallographic studies done on tyrosine hydroxylase and phenylalanine hydroxylase.

Mechanisms of tryptophan and tyrosine hydroxylase

The aromatic amino acid hydroxylases tryptophan hydroxylase and tyrosine hydroxylase are responsible for the initial steps in the formation of serotonin and the catecholamine neurotransmitters, respectively. Both enzymes are nonheme iron-dependent monooxygenases that catalyze the insertion of one atom of molecular oxygen onto the aromatic ring of their amino acid substrates, using a tetrahydropterin as a two electron donor to reduce the second oxygen atom to water. This review discusses the current understanding of the catalytic mechanism of these two enzymes. The reaction occurs as two sequential half reactions: a reaction between the active site iron, oxygen, and the tetrahydropterin to form a reactive Fe(IV) O intermediate and hydroxylation of the amino acid by the Fe(IV) O. The mechanism of formation of the Fe(IV) O is unclear; however, considerable evidence suggests the formation of an Fe(II) -peroxypterin intermediate. The amino acid is hydroxylated by the Fe(IV) O intermediate in an electrophilic aromatic substitution mechanism.

Mutagenesis of a specificity-determining residue in tyrosine hydroxylase establishes that the enzyme is a robust phenylalanine hydroxylase but a fragile tyrosine hydroxylase

The aromatic amino acid hydroxylases tyrosine hydroxylase (TyrH) and phenylalanine hydroxylase (PheH) have essentially identical active sites; however, PheH is nearly incapable of hydroxylating tyrosine, while TyrH can readily hydroxylate both tyrosine and phenylalanine. Previous studies have indicated that Asp425 of TyrH is important in determining the substrate specificity of that enzyme [Daubner, S. C., Melendez, J., and Fitzpatrick, P. F. (2000) Biochemistry 39, 9652-9661]. Alanine-scanning mutagenesis of amino acids 423-427, a mobile loop containing Asp425, shows that only mutagenesis of Asp425 alters the activity of the enzyme significantly. Saturation mutagenesis of Asp425 results in large (up to 10(4)) decreases in the V(max) and V(max)/K(tyr) values for tyrosine hydroxylation, but only small decreases or even increases in the V(max) and V(max)/K(phe) values for phenylalanine hydroxylation. The decrease in the tyrosine hydroxylation activity of the mutant proteins is due to an uncoupling of tetrahydropterin oxidation from amino acid hydroxylation with tyrosine as the amino acid substrate. In contrast, with the exception of the D425W mutant, the extent of coupling of tetrahydropterin oxidation and amino acid hydroxylation is unaffected or increases with phenylalanine as the amino acid substrate. The decrease in the V(max) value with tyrosine as the substrate shows a negative correlation with the hydrophobicity of the amino acid residue at position 425. The results are consistent with a critical role of Asp425 being to prevent a hydrophobic interaction that results in a restricted active site in which hydroxylation of tyrosine does not occur.

Evidence for a high-spin Fe(IV) species in the catalytic cycle of a bacterial phenylalanine hydroxylase

Phenylalanine hydroxylase is a mononuclear non-heme iron protein that uses tetrahydropterin as the source of the two electrons needed to activate dioxygen for the hydroxylation of phenylalanine to tyrosine. Rapid-quench methods have been used to analyze the mechanism of a bacterial phenylalanine hydroxylase from Chromobacterium violaceum. Mössbauer spectra of samples prepared by freeze-quenching the reaction of the enzyme-(57)Fe(II)-phenylalanine-6-methyltetrahydropterin complex with O(2) reveal the accumulation of an intermediate at short reaction times (20-100 ms). The Mössbauer parameters of the intermediate (δ = 0.28 mm/s, and |ΔE(Q)| = 1.26 mm/s) suggest that it is a high-spin Fe(IV) complex similar to those that have previously been detected in the reactions of other mononuclear Fe(II) hydroxylases, including a tetrahydropterin-dependent tyrosine hydroxylase. Analysis of the tyrosine content of acid-quenched samples from similar reactions establishes that the Fe(IV) intermediate is kinetically competent to be the hydroxylating intermediate. Similar chemical-quench analysis of a reaction allowed to proceed for several turnovers shows a burst of tyrosine formation, consistent with rate-limiting product release. All three data sets can be modeled with a mechanism in which the enzyme-substrate complex reacts with oxygen to form an Fe(IV)═O intermediate with a rate constant of 19 mM(-1) s(-1), the Fe(IV)═O intermediate hydroxylates phenylalanine with a rate constant of 42 s(-1), and rate-limiting product release occurs with a rate constant of 6 s(-1) at 5 °C.

Investigating inner-sphere reorganization via secondary kinetic isotope effects in the C-H cleavage reaction catalyzed by soybean lipoxygenase: tunneling in the substrate backbone as well as the transferred hydrogen

This work describes the application of NMR to the measurement of secondary deuterium (2° (2)H) and carbon-13 ((13)C) kinetic isotope effects (KIEs) at positions 9-13 within the substrate linoleic acid (LA) of soybean lipoxygenase-1. The KIEs have been measured using LA labeled with either protium (11,11-h2-LA) or deuterium (11,11-d2-LA) at the reactive C11 position, which has been previously shown to yield a primary deuterium isotope effect of ca. 80. The conditions of measurement yield the intrinsic 2° (2)H and (13)C KIEs on k(cat)/K(m) directly for 11,11-d2-LA, whereas the values for the 2° (2)H KIEs for 11,11-h2-LA are obtained after correction for a kinetic commitment. The pattern of the resulting 2° (2)H and (13)C isotope effects reveals values that lie far above those predicted from changes in local force constants. Additionally, many of the experimental values cannot be modeled by electronic effects, torsional strain, or the simple inclusion of a tunneling correction to the rate. Although previous studies have shown the importance of extensive tunneling for cleavage of the primary hydrogen at C11 of LA, the present findings can only be interpreted by extending the conclusion of nonclassical behavior to the secondary hydrogens and carbons that flank the position undergoing C-H bond cleavage. A quantum mechanical method introduced by Buhks et al. [J. Phys. Chem. 1981, 85, 3763] to model the inner-sphere reorganization that accompanies electron transfer has been shown to be able to reproduce the scale of the 2° (2)H KIEs.

Single turnover kinetics of tryptophan hydroxylase: evidence for a new intermediate in the reaction of the aromatic amino acid hydroxylases

Tryptophan hydroxylase (TrpH) uses a non-heme mononuclear iron center to catalyze the tetrahydropterin-dependent hydroxylation of tryptophan to 5-hydroxytryptophan. The reactions of the TrpH.Fe(II), TrpH.Fe(II).tryptophan, TrpH.Fe(II).6MePH(4).tryptophan, and TrpH.Fe(II).6MePH(4).phenylalanine complexes with O(2) were monitored by stopped-flow absorbance spectroscopy and rapid quench methods. The second-order rate constant for the oxidation of TrpH.Fe(II) has a value of 104 M(-1) s(-1) irrespective of the presence of tryptophan. Stopped-flow absorbance analyses of the reaction of the TrpH.Fe(II).6MePH(4).tryptophan complex with oxygen are consistent with the initial step being reversible binding of oxygen, followed by the formation with a rate constant of 65 s(-1) of an intermediate I that has maximal absorbance at 420 nm. The rate constant for decay of I, 4.4 s(-1), matches that for formation of the 4a-hydroxypterin product monitored at 248 nm. Chemical-quench analyses show that 5-hydroxytryptophan forms with a rate constant of 1.3 s(-1) and that overall turnover is limited by a subsequent slow step, presumably product release, with a rate constant of 0.2 s(-1). All of the data with tryptophan as substrate can be described by a five-step mechanism. In contrast, with phenylalanine as substrate, the reaction can be described by three steps: a second-order reaction with oxygen to form I, decay of I as tyrosine forms, and slow product release.

Measurement of the intramolecular isotope effect on aliphatic hydroxylation by Chromobacterium violaceum phenylalanine hydroxylase

The non-heme iron enzyme phenylalanine hydroxylase from Chromobacterium violaceum has previously been shown to catalyze the hydroxylation of benzylic and aliphatic carbons in addition to the normal aromatic hydroxylation reaction. The intrinsic isotope effect for hydroxylation of 3-cyclochexylalanine by the enzyme was determined in order to gain insight into the reactivity of the iron center. With 3-[(2)H(11)-cyclohexyl]alanine as the substrate, the isotope effect on the k(cat) value was 1, consistent with an additional step in the overall reaction being significantly slower than hydroxylation. Consequently, the isotope effect was determined as an intramolecular effect by measuring the amount of deuterium lost in the hydroxylation of 3-[1,2,3,4,5,6-(2)H(6)-cyclohexyl]alanine. The ratio of 4-HO-cyclohexylalanine that retained deuterium to that which lost one deuterium atom was 2.8. This gave a calculated value of 12.6 for the ratio of the primary deuterium kinetic isotope effect to the secondary isotope effect. This value is consistent with hydrogen atom abstraction by an electrophilic Fe(O) center and a contribution of quantum-mechanical tunneling to the reaction.

Redox biocatalysis and metabolism: molecular mechanisms and metabolic network analysis

Whole-cell biocatalysis utilizes native or recombinant enzymes produced by cellular metabolism to perform synthetically interesting reactions. Besides hydrolases, oxidoreductases represent the most applied enzyme class in industry. Oxidoreductases are attributed a high future potential, especially for applications in the chemical and pharmaceutical industries, as they enable highly interesting chemistry (e.g., the selective oxyfunctionalization of unactivated C-H bonds). Redox reactions are characterized by electron transfer steps that often depend on redox cofactors as additional substrates. Their regeneration typically is accomplished via the metabolism of whole-cell catalysts. Traditionally, studies towards productive redox biocatalysis focused on the biocatalytic enzyme, its activity, selectivity, and specificity, and several successful examples of such processes are running commercially. However, redox cofactor regeneration by host metabolism was hardly considered for the optimization of biocatalytic rate, yield, and/or titer. This article reviews molecular mechanisms of oxidoreductases with synthetic potential and the host redox metabolism that fuels biocatalytic reactions with redox equivalents. The tools discussed in this review for investigating redox metabolism provide the basis for studies aiming at a deeper understanding of the interplay between synthetically active enzymes and metabolic networks. The ultimate goal of rational whole-cell biocatalyst engineering and use for fine chemical production is discussed.

Regulation of phenylalanine hydroxylase: conformational changes upon phenylalanine binding detected by hydrogen/deuterium exchange and mass spectrometry

Phenylalanine acts as an allosteric activator of the tetrahydropterin-dependent enzyme phenylalanine hydroxylase. Hydrogen/deuterium exchange monitored by mass spectrometry has been used to gain insight into local conformational changes accompanying activation of rat phenylalanine hydroxylase by phenylalanine. Peptides in the regulatory and catalytic domains that lie in the interface between these two domains show large increases in the extent of deuterium incorporation from solvent in the presence of phenylalanine. In contrast, the effects of phenylalanine on the exchange kinetics of a mutant enzyme lacking the regulatory domain are limited to peptides surrounding the binding site for the amino acid substrate. These results support a model in which the N-terminus of the protein acts as an inhibitory peptide, with phenylalanine binding causing a conformational change in the regulatory domain that alters the interaction between the catalytic and regulatory domains.

Measurement of intrinsic rate constants in the tyrosine hydroxylase reaction

Tyrosine hydroxylase (TyrH) is a pterin-dependent mononuclear non-heme aromatic amino acid hydroxylase that catalyzes the conversion of tyrosine to dihydroxyphenylalanine (DOPA). Chemical quench analyses of the enzymatic reaction show a burst of DOPA formation, followed by a linear rate equal to the k(cat) value at both 5 and 30 degrees C. The effects of increasing solvent viscosity confirm that k(cat) is approximately 84% limited by diffusion, most probably due to slow product release, and that tyrosine has a commitment to catalysis of 0.45. The effect of viscosity on the k(cat)/K(m) for 6-methyltetrahydropterin is greater than the theoretical limit, consistent with the coupling of pterin binding to the movement of a surface loop. The absorbance changes in the spectrum of the tetrahydropterin during the first turnover, the kinetics of DOPA formation during the first turnover, and the previously described kinetics for formation and decay of the Fe(IV)O intermediate [Eser, B. E., Barr, E. W., Frantom, P. A., Saleh, L., Bollinger, J. M., Jr., Krebs, C., and Fitzpatrick, P. F. (2007) J. Am. Chem. Soc. 129, 11334-11335] were analyzed globally, yielding a single set of rate constants for the TyrH reaction. Reversible binding of oxygen is followed by formation of Fe(IV)O and 4a-hydroxypterin with a rate constant of 13 s(-1) at 5 degrees C. Transfer of oxygen from Fe(IV)O to tyrosine to form DOPA follows with a rate constant of 22 s(-1). Release of DOPA and/or the 4a-hydroxypterin with a rate constant of 0.86 s(-1) completes the turnover.

Demonstration of a peroxide shunt in the tetrahydropterin-dependent aromatic amino acid monooxygenases

The nonheme iron enzyme phenylalanine hydroxylase, tyrosine hydroxylase, and tryptophan hydroxylase catalyze the hydroxylation of their aromatic amino acid substrates using a tetrahydropterin as the source of electrons. The hydroxylating intermediate is proposed to be an Fe(IV)O species. We report here that all three enzymes will catalyze hydroxylation reactions using H(2)O(2) in place of tetrahydropterin and oxygen, forming tyrosine and 3-hydroxyphenylalanine from phenylalanine, 4-HOCH(2)-phenylalanine from 4-CH(3)-phenylalanine, and hydroxycyclohexylalanine from 3-cyclohexylalanine. No peroxide-dependent reaction is seen with active site mutants of TyrH and PheH in which the stability or reactivity of the iron center is compromised. These results provide further support for an Fe(IV)O hydroxylating intermediate.

Kinetic isotope effects on aromatic and benzylic hydroxylation by Chromobacterium violaceum phenylalanine hydroxylase as probes of chemical mechanism and reactivity

Phenylalanine hydroxylase from Chromobacterium violaceum (CvPheH) is a non-heme iron monooxygenase that catalyzes the hydroxylation of phenylalanine to tyrosine. In this study, we used deuterium kinetic isotope effects to probe the chemical mechanisms of aromatic and benzylic hydroxylation to compare the reactivities of bacterial and eukaryotic aromatic amino acid hydroxylases. The (D) k cat value for the reaction of CvPheH with [(2)H 5]phenylalanine is 1.2 with 6-methyltetrahydropterin and 1.4 with 6,7-dimethyltetrahydropterin. With the mutant enzyme I234D, the (D) k cat value decreases to 0.9 with the latter pterin; this is likely to be the intrinsic effect for addition of oxygen to the amino acid. The isotope effect on the subsequent tautomerization of a dienone intermediate was determined to be 5.1 by measuring the retention of deuterium in tyrosine produced from partially deuterated phenylalanine; this large isotope effect is responsible for the normal effect on k cat. The isotope effect for hydroxylation of the methyl group of 4-CH 3-phenylalanine, obtained from the partitioning of benzylic and aromatic hydroxylation products, is 10. The temperature dependence of this isotope effect establishes the contribution of hydrogen tunneling to benzylic hydroxylation by this enzyme. The results presented here provide evidence that the reactivities of the prokaryotic and eukaryotic hydroxylases are similar and further define the reactivity of the iron center for the family of aromatic amino acid hydroxylases.

Characterization of metal ligand mutants of phenylalanine hydroxylase: Insights into the plasticity of a 2-histidine-1-carboxylate triad

The iron atom in the nonheme iron monooxygenase phenylalanine hydroxylase is bound on one face by His285, His290, and Glu330. This arrangement of metal ligands is conserved in the other aromatic amino acid hydroxylases, tyrosine hydroxylase and tryptophan hydroxylase. A similar 2-His-1-carboxylate facial triad of two histidines and an acidic residue are the ligands to the iron in other nonheme iron enzymes, including the alpha-ketoglutarate-dependent hydroxylases and the extradiol dioxygenases. Previous studies of the effects of conservative mutations of the iron ligands in tyrosine hydroxylase established that there is some plasticity in the nature of the ligands and that the three ligands differ in their sensitivity to mutagenesis. To determine the generality of this finding for enzymes containing a 2-His-1-carboxylate facial triad, the His285, His290, and Glu330 in rat phenylalanine hydroxylase were mutated to glutamine, glutamate, and histidine. All of the mutant proteins had low but measurable activities for tyrosine formation. In general, mutation of Glu330 had the greatest effect on activity and mutation of His290 the least. All of the mutations resulted in an excess of tetrahydropterin oxidized relative to tyrosine formation, with mutation of His285 having the greatest effect on the coupling of the two partial reactions. The H285Q enzyme had the highest activity as tetrahydropterin oxidase at 20% the wild-type value. All of the mutations greatly decreased the affinity for iron, with mutation of Glu330 the most deleterious. The results complement previous results with tyrosine hydroxylase in establishing the plasticity of the individual iron ligands in this enzyme family.

Car−Parrinello Molecular Dynamics Study of Anharmonic Systems: A Mannich Base in Solution

A Car-Parrinello molecular dynamics study was performed for 4,5-dimethyl-2-(N,N-dimethylaminomethyl)phenol, a Mannich base, to investigate the vibrational properties in solution of its intramolecular hydrogen bond. The dynamic behavior of this hydrogen-bonded system was investigated using an explicit solvent model. Addition of a nonpolar solvent permitted inclusion of delicate environmental effects on the strongly anharmonic system which was studied from first principles. Molecular dynamics and a posteriori quantization of the O-H motion were applied to reproduce the vibrational features of the O-H stretching mode. Consistent application of Car-Parrinello dynamics based on the density functional theory with subsequent solution of the vibrational Schrödinger equation for the O-H stretching motion offers an effective method for strongly anharmonic systems, and this is supported by the comparison of the results with experimental spectra. As a further element of the intramolecular hydrogen bond study, the effects of deuteration were taken into account and a successful application of the O-H stretching mode quantization technique to the liquid phase is demonstrated. This provides a valuable computational methodology for investigations incorporating nuclear quantum effects in the liquid phase and enzyme active centers and can be used to investigate numerous systems that are not readily susceptible to experimental analysis.

Insights into the mechanisms of flavoprotein oxidases from kinetic isotope effects

Deuterium, solvent, and (15)N kinetic isotope effects have been used to probe the mechanisms by which flavoproteins oxidize carbon-oxygen and carbon-nitrogen bonds in amines, hydroxy acids, and alcohols. For the amine oxidases d-amino acid oxidase, N-methyltryptophan oxidase, and tryptophan monooxygenase, d-serine, sarcosine, and alanine are slow substrates for which CH bond cleavage is fully rate limiting. Inverse isotope effects for each of 0.992-0.996 are consistent with a common mechanism involving hydride transfer from the uncharged amine. Computational analyses of possible mechanisms support this conclusion. Deuterium and solvent isotope effects with wild-type and mutant variants of the lactate dehydrogenase flavocytochrome b(2) show that OH and CH bond cleavage are not concerted, but become so in the Y254F enzyme. This is consistent with a highly asynchronous reaction in which OH bond cleavage precedes hydride transfer. The results of Hammett analyses and solvent and deuterium isotope effects support a similar mechanism for alcohol oxidase.