Hydrogen tunneling in enzyme reactions

Structures of horse liver alcohol dehydrogenase complexed with NAD+ and substituted benzyl alcohols

Structures of the enzyme complexed with NAD+ and 2,3,4,5,6-pentafluorobenzyl alcohol were determined by X-ray crystallography at a resolution of 2.1 A and to a refinement R value of 18.3% for a monoclinic (P2(1)) form and to 2.4 A and an R value of 18.9% for a triclinic crystal form. The pentafluorobenzyl alcohol does not react, due to electron withdrawal by the fluorine atoms. A structure with NAD+ and p-bromobenzyl alcohol in the monoclinic form was also determined at 2.5 A and an R value of 16.7%. The conformations of the subunits in the monoclinic and triclinic crystal forms are very similar. The dimer is the asymmetric unit, and a rigid body rotation closes the cleft between the coenzyme and catalytic domains upon complex formation. In the monoclinic form, this conformational change is described by a rotation of 9 degrees in one subunit and 10 degrees in the other. The pentafluoro- and p-bromobenzyl alcohols bind in overlapping positions. The hydroxyl group of each alcohol is ligated to the catalytic zinc and participates in an extensive hydrogen-bonded network that includes the imidazole group of His-51, which can act as a base and shuttle a proton to solvent. The hydroxymethyl carbon of the pentafluorobenzyl alcohol is 3.4 A from C4 of the nicotinamide ring, and the pro-R hydrogen is in a good position for direct transfer to C4. The p-bromobenzyl alcohol may react after small rotations around single bonds of the alcohol. These structures should approximate the active Michaelis-Menten complexes.

Role of tyrosine 143 in lactate dehydrogenation by flavocytochrome b2. Primary kinetic isotope effect studies with a phenylalanine mutant

Flavocytochrome b2 catalyzes the oxidation of lactate at the expense of cytochrome c. After flavin (FMN) reduction by the substrate, reducing equivalents are transferred one by one to heme b2, and from there on to cytochrome c. The crystal structure of the enzyme is known at 2.4-A resolution, and specific roles in catalysis have been assigned to active side chains. Tyr143 in particular, located at the interface between the flavodehydrogenase moiety and the heme-binding domain, was thought to take part in substrate binding, as well as to orient the heme-binding domain for efficient electron transfer. A first study of the properties of a Tyr143Phe mutant showed that the major effect of the mutation was to decrease the rate of electron transfer from flavin to heme [Miles, C.S., Rouvière-Fourmy, N., Lederer, F., Mathews, F.S., Reid, G.A., Black, M.T., & Chapman, S.K. (1992) Biochem. J. 285, 187-192]. In the present paper, we focus on the effect of the mutation on catalysis of lactate dehydrogenation. We report the deuterium kinetic isotope effects on flavin reduction as measured with stopped-flow methods and on cytochrome c reduction in the steady-state using L-[2-2H]lactate. For the wild-type enzyme, isotope effects on FMN reduction, D(kredF) and D(kredF)/Km), were 7.2 +/- 0.9 and 4.2 +/- 1.3, respectively, and for the Y143F mutant values of 4.4 +/- 0.5 and 3.9 +/- 1.1 were obtained. Calculations, from deuterium isotope effects, of substrate Kd values, combined with knowledge of kcat/Km values, lead to the conclusion that Tyr143 does stabilize the Michaelis complex by hydrogen bonding to a substrate carboxylate, as was postulated; but the mutation does not destabilize the transition state more than the Michaelis complex.(ABSTRACT TRUNCATED AT 250 WORDS)

Molecular modeling methods. Basic techniques and challenging problems

An overview is presented of computer modeling and simulation methods that play an increasing role in drug design: quantum chemical methods, molecular mechanics, molecular dynamics and Brownian dynamics. The application of molecular dynamics for the prediction of thermodynamic properties like free energy differences and binding constants is discussed. The Brownian dynamics method is presented in connection with the calculation of effective electrostatic forces using the Poisson-Boltzmann equation, which allows one to sample ligand-binding geometries and to predict the kinetics of diffusion-limited enzyme reactions. New techniques that have recently been extensively developed, such as the global energy minimization and quantum-classical dynamics methods, are also introduced. The molecular modeling methods are illustrated with selected examples.

Vibrationally enhanced tunneling as a mechanism for enzymatic hydrogen transfer

We present a theory of enzymatic hydrogen transfer in which hydrogen tunneling is mediated by thermal fluctuations of the enzyme's active site. These fluctuations greatly increase the tunneling rate by shortening the distance the hydrogen must tunnel. The average tunneling distance is shown to decrease when heavier isotopes are substituted for the hydrogen or when the temperature is increased, leading to kinetic isotope effects (KIEs)--defined as the factor by which the reaction slows down when isotopically substituted substrates are used--that need be no larger than KIEs for nontunneling mechanisms. Within this theory we derive a simple KIE expression for vibrationally enhanced ground state tunneling that is able to fit the data for the bovine serum amine oxidase (BSAO) system, correctly predicting the large temperature dependence of the KIEs. Because the KIEs in this theory can resemble those for nontunneling dynamics, distinguishing the two possibilities requires careful measurements over a range of temperatures, as has been done for BSAO.

The consequences of the isotope effect on proline dehydrogenation rates estimated by the tritium loss method

Loss of tritium from a substrate is often used to estimate the rate of dehydrogenation. However, loss of 3H may be much slower than loss of H because of the tritium isotope effect. In order to assess the impact of the tritium isotope effect, loss of 3H from the C-5 position of proline during dehydrogenation by rat liver mitochondria and bacteroids from soybean (Glycine max [L.] Merrill) nodules was compared with appearance of 14C in products of [14C]proline dehydrogenation. Incubations were carried out in the presence of o-aminobenzaldehyde (added to trap the initial product, delta 1-pyroline-5-carboxylate). The fraction of total 14C products trapped by o-aminobenzaldehyde varied from 0.07 to 0.75 depending upon experimental conditions. With rat liver mitochondria, dehydrogenation of [14C]proline was between 3.27 and 9.25 times faster than dehydrogenation of 3H proline, depending upon assay conditions. Soybean nodule bacteriods dehydrogenated [14C]proline about 5 times faster than [3H]proline. We conclude the following: (i) the rate of proline dehydrogenation may be greatly underestimated by the tritium assay because of the tritium isotope effect, and (ii) the 14C assay may underestimate the rate of proline dehydrogenation if it is assumed that o-aminobenzaldehyde quantitatively traps delta 1-pyrroline-5-carboxylate under all conditions. The simplicity of the tritium assay makes it attractive for routine use. However, its use requires determination of the tritium isotope effect, under the specific conditions of the assay, in order to correct the results. The considerations discussed here have broad applicability to any dehydrogenase assay employing tritium loss.

Human erythrocyte glutathione reductase: chemical mechanism and structure of the transition state for hydride transfer

Kinetic parameters and primary deuterium kinetic isotope effects for NADH and five pyridine nucleotide substrates have been determined at pH 8.1 for human erythrocyte glutathione reductase. DV/KNADH and DV are equal to 1.4 and are pH independent below pH 8.1, but DV decreases to 1.0 at high pH as a group exhibiting a pK of 8.6 is deprotonated. This result suggests that as His-467' is deprotonated, the rate of the isotopically insensitive oxidative half-reaction is specifically decreased and becomes rate-limiting. For all substrates, equivalent V and V/K primary deuterium kinetic isotope effects are observed at pH values below 8.1. The primary deuterium kinetic isotope effect on V, but not V/K, is sensitive to solvent isotopic composition. The primary tritium kinetic isotope effects agree well with the corresponding value calculated from the primary deuterium kinetic isotope effects by using the Swain-Schaad relationship. This suggests that the primary deuterium kinetic isotope effects observed in these steady-state experiments are the intrinsic primary deuterium kinetic isotope effects for hydride transfer. The magnitude of the primary deuterium kinetic isotope effect is dependent on the redox potential of the pyridine nucleotide substrate used, varying from approximately 1.4 for NADH and -320 mV reductants to 2.7 for thioNADH to 4.2-4.8 for 3-acetylpyridine adenine dinucleotide (3APADH). The alpha-secondary tritium kinetic isotope effects also increase as the redox potential of the pyridine nucleotide substrate becomes more positive. Together, these data indicate that the transition state for hydride transfer is very early for NADH and becomes later for thioNADH and 3APADH, as predicted by Hammond's postulate.(ABSTRACT TRUNCATED AT 250 WORDS)

Kinetic isotope effect studies on milk xanthine oxidase and on chicken liver xanthine dehydrogenase

The effect of isotopic substitution of the 8-H of xanthine (with 2H and 3H) on the rate of oxidation by bovine xanthine oxidase and by chicken xanthine dehydrogenase has been measured. V/K isotope effects were determined from competition experiments. No difference in H/T(V/K) values was observed between xanthine oxidase (3.59 +/- 0.1) and xanthine dehydrogenase (3.60 +/- 0.09). Xanthine dehydrogenase exhibited a larger T/D(V/K) value (0.616 +/- 0.028) than that observed for xanthine oxidase (0.551 +/- 0.016). Observed H/T(V/K) values for either enzyme are less than those H/T(V/K) values calculated with D/T(V/K) data. These discrepancies are suggested to arise from the presence of a rate-limiting step(s) prior to the irreversible C-H bond cleavage step in the mechanistic pathways of both enzymes. These kinetic complexities preclude examination of whether tunneling contributes to the reaction coordinate for the H-transfer step in each enzyme. No observable exchange of tritium with solvent is observed during the anaerobic incubation of [8-3H]xanthine with either enzyme, which suggests the reverse commitment to catalysis (Cr) is essentially zero. With the assumption of adherence to reduced mass relationships, the intrinsic deuterium isotope effect (Dk) for xanthine oxidation is calculated to be 7.4 +/- 0.7 for xanthine oxidase and 4.2 +/- 0.2 for xanthine dehydrogenase. By use of these values and steady-state kinetic data, the minimal rate for the hydrogen-transfer step is calculated to be approximately 75-fold faster than kcat for xanthine oxidase and approximately 10-fold faster than kcat for xanthine dehydrogenase. This calculated rate is consistent with data obtained by rapid-quench experiments with XO. A stoichiometry of 1.0 +/- 0.3 mol of uric acid/mol of functional enzyme is formed within the mixing time of the instrument (5-10 ms). The kinetic isotope effect data also permitted the calculation of the Kd values [Klinman, J. P., & Mathews, R. G. (1985) J. Am. Chem. Soc. 107, 1058-1060] for substrate dissociation, including all reversible steps prior to C-H bond cleavage. Values calculated for each enzyme (Kd = 120 microM) were found to be identical within experimental uncertainty.

Quantum mechanical effects in enzyme-catalysed hydrogen transfer reactions

Hydrogen tunneling at room temperature has recently been demonstrated in the reactions catalysed by yeast alcohol dehydrogenase and bovine serum amine oxidase. These results suggest that quantum mechanical effects may be a fundamental and pervasive feature of enzyme-catalysed hydrogen transfer reactions. Our ability to detect such behavior introduces a new probe for the investigation of reaction barrier shape and protein dynamics in enzyme catalysis.