Enzyme mechanisms from molecular modeling and isotope effects

The application of kinetic isotope effects and molecular modeling to characterize three enzyme-catalyzed reactions is presented; the mechanism of the chloroacid dehalogenase catalyzed reaction is approached using chlorine kinetic isotope effects and solvent kinetic isotope effects. The pre-steady-state phase of the reaction catalyzed by methylmalonyl-CoA mutase is approached by different QM/MM schemes and the results are validated by comparison with the experimental value of the deuterium kinetic isotope effect. Finally, a procedure for improving QM/MM calculations is illustrated by analysis of the trihydroxynaphthalene reductase-catalyzed reaction.     

Tritium secondary kinetic isotope effect on phenylalanine ammonia-lyase-catalyzed reaction.

The mechanism by which phenylalanine ammonia-lyase (PAL, EC 4.3.1.5) catalyzes the reversible elimination of ammonia from phenylalanine yielding (E)-cinnamic acid has gained much attention in the recent years. Dehydroalanine is essential for the catalysis. It was assumed that this prostetic group acts as the electrophile, leading to a covalently bonded enzyme-intermediate complex with quarternary nitrogen of phenylalanine. Recently, an alternative mechanism has been suggested in which the enzyme-intermediate complex is formed in a Friedel-Crafts reaction between dehydroalanine and orthocarbon of the aromatic ring. Using semiempirical calculations we have shown that these two alternative mechanisms can be distinguished on the basis of the hydrogen secondary kinetic isotope effect when tritium label is placed in the orthopositions. Our calculations indicated also that the kinetic isotope effect measured using ring-labeled d(5)-phenylalanine could not be used to differentiate these alternative mechanisms. Measured secondary tritium kinetic isotope effect shows strong dependence on the reaction progress, starting at the inverse value of k(H)/k(T) = 0.85 for 5% conversion and reaching the normal value of about 1.15 as the conversion increases to 20%. This dependence has been interpreted in terms of a complex mechanism with initial formation of the Friedel-Crafts type intermediate.     

pH and kinetic isotope effects in d-amino acid oxidase catalysis.

The effects of pH, solvent isotope, and primary isotope replacement on substrate dehydrogenation by Rhodotorula gracilis d-amino acid oxidase were investigated. The rate constant for enzyme-FAD reduction by d-alanine increases approximately fourfold with pH, reflecting apparent pKa values of approximately 6 and approximately 8, and reaches plateaus at high and low pH. Such profiles are observed in all presteady-state and steady-state kinetic experiments, using both d-alanine and d-asparagine as substrates, and are inconsistent with the operation of a base essential to catalysis. A solvent deuterium isotope effect of 3.1 +/- 1.1 is observed on the reaction with d-alanine at pH 6; it decreases to 1.2 +/- 0.2 at pH 10. The primary substrate isotope effect on the reduction rate with [2-D]d-alanine is 9.1 +/- 1.5 at low and 2.3 +/- 0.3 at high pH. At pH 6.0, the solvent isotope effect is 2.9 +/- 0.8 with [2-D]d-alanine, and the primary isotope effect is 8.4 +/- 2.4 in D2O. Thus, primary and solvent kinetic isotope effects (KIEs) are independent of the presence of the other isotope, i.e. the 'double' kinetic isotope effect is the product of the individual KIEs, consistent with a transition state in which rupture of the two bonds of the substrate to hydrogen is concerted. These results support a hydride transfer mechanism for the dehydrogenation reaction in d-amino acid oxidase and argue against the occurrence of any intermediates in the process. A pKa,app of approximately 8 is interpreted to arise from the microscopic ionization of the substrate amino acid alpha-amino group, but also includes contributions from kinetic parameters.     

Experimental and computational investigations of oxygen reactivity in a heme and tyrosyl radical-containing fatty acid α-(di)oxygenase

Rice α-(di)oxygenase mediates the regio- and stereospecific oxidation of fatty acids using a persistent catalytic tyrosyl radical. Experiments conducted in the physiological O(2) concentration range, where initial hydrogen atom abstraction from the fatty acid occurs in a kinetically reversible manner, are described. Our findings indicate that O(2)-trapping of an α-carbon radical is likely to reversibly precede reduction of a 2-(R)-peroxyl radical intermediate in the first irreversible step. A mechanism of concerted proton-coupled electron transfer is proposed on the basis of natural abundance oxygen-18 kinetic isotope effects, deuterium kinetic isotope effects, and calculations at the density functional level of theory, which predict a polarized transition state in which electron transfer is advanced to a greater extent than proton transfer. The approach outlined should be useful for identifying mechanisms of concerted proton-coupled electron transfer in a variety of oxygen-utilizing enzymes.     

Nature of oxygen activation in glucose oxidase from Aspergillus niger: the importance of electrostatic stabilization in superoxide formation.

Glucose oxidase catalyzes the oxidation of glucose by molecular dioxygen, forming gluconolactone and hydrogen peroxide. A series of probes have been applied to investigate the activation of dioxygen in the oxidative half-reaction, including pH dependence, viscosity effects, 18O isotope effects, and solvent isotope effects on the kinetic parameter Vmax/Km(O2). The pH profile of Vmax/Km(O2) exhibits a pKa of 7.9 +/- 0.1, with the protonated enzyme form more reactive by 2 orders of magnitude. The effect of viscosogen on Vmax/Km(O2) reveals the surprising fact that the faster reaction at low pH (1.6 x 10(6) M-1 s-1) is actually less diffusion-controlled than the slow reaction at high pH (1.4 x 10(4) M-1 s-1); dioxygen reduction is almost fully diffusion-controlled at pH 9.8, while the extent of diffusion control decreases to 88% at pH 9.0 and 32% at pH 5.0, suggesting a transition of the first irreversible step from dioxygen binding at high pH to a later step at low pH. The puzzle is resolved by 18O isotope effects. 18(Vmax/Km) has been determined to be 1.028 +/- 0.002 at pH 5.0 and 1.027 +/- 0.001 at pH 9.0, indicating that a significant O-O bond order decrease accompanies the steps from dioxygen binding up to the first irreversible step at either pH. The results at high pH lead to an unequivocal mechanism; the rate-limiting step in Vmax/Km(O2) for the deprotonated enzyme is the first electron transfer from the reduced flavin to dioxygen, and this step accompanies binding of molecular dioxygen to the active site. In combination with the published structural data, a model is presented in which a protonated active site histidine at low pH accelerates the second-order rate constant for one electron transfer to dioxygen through electrostatic stabilization of the superoxide anion intermediate. Consistent with the proposed mechanisms for both high and low pH, solvent isotope effects indicate that proton transfer steps occur after the rate-limiting step(s). Kinetic simulations show that the model that is presented, although apparently in conflict with previous models for glucose oxidase, is in good agreement with previously published kinetic data for glucose oxidase. A role for electrostatic stabilization of the superoxide anion intermediate, as a general catalytic strategy in dioxygen-utilizing enzymes, is discussed.     

Transition-state structures for N-glycoside hydrolysis of AMP by acid and by AMP nucleosidase in the presence and absence of allosteric activator

The mechanism of acid and enzymatic hydrolysis of the N-glycosidic bond of AMP has been investigated by fitting experimentally observed kinetic isotope effects [Parkin, D. W., & Schramm, V. L. (1987) Biochemistry (preceding paper in this issue)] to calculated kinetic isotope effects for proposed transition-state structures. The sensitivity of the transition-state calculations was tested by "arying the transition-state structure and comparing changes in the calculated kinetic isotope effects with the experimental values of the isotope effect measurements. The kinetic isotope effects for the acid-catalyzed hydrolysis of AMP are best explained by a transition state with considerable oxycarbonium character in the ribose ring, significant bonding remaining to the departing adenine ring, participation of a water nucleophile, and protonation of the adenine ring. A transition-state structure without preassociation of the water nucleophile cannot be eliminated by the data. Enzymatic hydrolysis of the N-glycosidic bond of AMP by AMP nucleosidase from Azotobacter vinelandii was analyzed in the absence and presence of MgATP, the allosteric activator that increases Vmax approximately 200-fold. The transition states for enzyme-catalyzed hydrolysis that best explain the kinetic isotope effects involve early SN1 transition states with significant bond order in the glycosidic bond and protonation of the adenine base. The enzyme enforces participation of an enzyme-bound water molecule, which has weak bonding to C1' in the transition state. Activation of AMP nucleosidase by MgATP causes the bond order of the glycosidic bond in the transition state to increase significantly. Hyperconjugation in the ribosyl group is altered by enzymatic stabilization of the oxycarbonium ion. This change is consistent with the interaction of an amino acid on the enzyme. Together, these changes stabilize a carboxonium-like transition-state complex that occurs earlier in the reaction pathway than in the absence of allosteric activator. In addition to the allosteric changes that alter transition-state structure, the presence of other inductive effects that are unobserved by kinetic isotope measurements is also likely to increase the catalytic rate.     

Experimental Evidence for a Hydride Transfer Mechanism in Plant Glycolate Oxidase Catalysis

In plants, glycolate oxidase is involved in the photorespiratory cycle, one of the major fluxes at the global scale. To clarify both the nature of the mechanism and possible differences in glycolate oxidase enzyme chemistry from C₃ and C₄ plant species, we analyzed kinetic parameters of purified recombinant C₃ (Arabidopsis thaliana) and C₄ (Zea mays) plant enzymes and compared isotope effects using natural and deuterated glycolate in either natural or deuterated solvent. The ¹²C/¹³C isotope effect was also investigated for each plant glycolate oxidase protein by measuring the ¹³C natural abundance in glycolate using natural or deuterated glycolate as a substrate. Our results suggest that several elemental steps were associated with an hydrogen/deuterium isotope effect and that glycolate α-deprotonation itself was only partially rate-limiting. Calculations of commitment factors from observed kinetic isotope effect values support a hydride transfer mechanism. No significant differences were seen between C₃ and C₄ enzymes.     

Brassinin oxidase mediated transformation of the phytoalexin brassinin: structure of the elusive co-product, deuterium isotope effect and stereoselectivity

Brassinin oxidase, a fungal detoxifying enzyme that mediates the conversion of the phytoalexin brassinin into indole-3-carboxaldehyde, is the first enzyme described to date that catalyzes the transformation of a dithiocarbamate group into an aldehyde equivalent. Brassinin is an essential phytoalexin due to its antifungal activity and its role as biosynthetic precursor of other phytoalexins produced in plants of the family Brassicaceae (common name crucifer). In this report, the isolation, structure determination and synthesis of the elusive co-product of brassinin transformation by brassinin oxidase, S-methyl dithiocarbamate, the syntheses of dideuterated and (R) and (S) monodeuterated brassinins, kinetic analyses of isotope effects and chemical modifications of brassinin oxidase are described. The reaction of [1'-(2)H(2)]brassinin was found to be slowed by a kinetic isotope effect of 5.3 on the value of k(cat)/K(m). This result indicates that the hydride/hydrogen transfer step preceding brassinin transformation is rate determining in the overall reaction. In addition, the use of (R) and (S)-[1'-(2)H]brassinins as substrates indicated that the hydride/hydrogen transfer step is ca. 88% stereoselective for the pro-R hydrogen. A detailed chemical mechanism of the enzymatic transformation of brassinin is proposed.     

A general method for the analysis of random bisubstrate enzyme mechanisms

In the present communication, a general method for the kinetic analysis of random bisubstrate mechanisms is described. The method comprises a stepwise application of the following kinetic and ligand-binding experiments: determination of steady-state kinetic constants, product inhibition patterns, maximum rate relationships, application of alternate substrates, application of dead-end inhibitors, direct binding of substrates, kinetic isotope effects, and isotope exchange studies. This general method was applied to a practical example: a yeast alcohol dehydrogenase-catalyzed oxidation of 2-propanol by NAD⁺ at pH 7.0, 25°C. It was found that this fully reversible reaction proceeds by a steady-state random Bi-Bi mechanism, whereby both dead-end complexes are formed.     

Use of isotope effects to elucidate enzyme mechanisms

The chemical bond breaking steps are normally not rate limiting for enzymatic reactions. However, comparison of deuterium and tritium isotope effects on the same reaction, especially when coupled with 13C isotope effects for the same step measured with deuterated as well as unlabeled substrates, allows calculation of the intrinsic isotope effects on the bond breaking steps and thus a determination of the commitments to catalysis for the reactants. The variation in observed isotope effects as a function of reactant concentration can be used to determine kinetic mechanisms, while the pH variation of isotope effects can determine the stickiness of the reactants and which portions of the reactant mechanism are pH dependent. Finally the size of primary and secondary intrinsic isotope effects can be used to determine transition state structure.     

Stabilization of an intermediate in the oxidative half-reaction of human liver glycolate oxidase

Glycolate oxidase is a flavin-dependent enzyme that catalyzes the oxidation of α-hydroxy acids to the corresponding α-keto acids, with reduction of molecular oxygen to hydrogen peroxide. A number of probes have been used to investigate the oxidative half-reaction catalyzed by the enzyme, including steady state and rapid kinetics, pH studies, solvent kinetic isotope effects, and solvent viscosity effects. Here we present the first spectroscopic evidence of the formation of an intermediate with absorbance features resembling those of a flavosemiquinone in the oxidative half-reaction of glycolate oxidase.     

Deuterium isotope effects in norcamphor metabolism by cytochrome P-450cam: kinetic evidence for the two-electron reduction of a high-valent iron-oxo intermediate

The kinetics of NADH consumption, oxygen uptake, and hydrogen peroxide production have been studied for norcamphor metabolism by cytochrome P-450cam. The kinetic deuterium isotope effects on these processes, with specifically deuteriated norcamphor, are 0.77, 1.22, and 1.16, respectively. Steady-state UV-visible spectroscopy indicates that transfer of the second electron to the dioxy ferrous P-450 is the rate-limiting step, as it is when camphor is the substrate. The inverse deuterium isotope effect for NADH consumption is consistent with an isotope-dependent branching between monooxygenase and oxidase activity, where these reactivities differ in their NADH:oxygen stoichiometries. However, no isotope-dependent redistribution of steady-state intermediates was detected by isotopic difference UV-visible spectroscopy in the presence of norcamphor. The kinetic isotope effects and steady-state spectral results suggest that the high-valent iron-oxo hydroxylating intermediate [FeO]3+ is reduced by NADH and the physiological electron-transfer proteins to afford water.     

Isotopic probes yield microscopic constants: separation of binding energy from catalytic efficiency in the bovine plasma amine oxidase reaction

Bovine plasma amine oxidase catalyzes the oxidative deamination of primary amines. The reaction can be viewed as two half-reactions: enzyme reduction by substrate followed by enzyme reoxidation by dioxygen. Anaerobic stopped-flow kinetic measurements of the first half-reaction indicate very large deuterium isotope effects for benzylamine, m-tyramine, and dopamine, Dk = 13.5 +/- 1.3, which are ascribed to an intrinsic isotope effect. From the insensitivity of these isotope effects to amine concentration, stopped-flow data provide substrate dissociation constants, K1, and rate constants for the C-H bond cleavage step, k3, directly. Steady-state isotope effects have also been measured for benzylamine and six ring-substituted phenethylamines. Whereas a small range of values for kcat, 0.38-1.2 s-1, and Dkcat, 5.4-8.8, is observed, kcat/Km = 1.3 X 10(2) to 3.8 X 10(4) M-1 S-1 and D(kcat/Km) = 5.6-16.1 indicate a marked effect of ring substituent. As described earlier [Miller, S., & Klinman, J.P. (1982) Methods Enzymol. 87, 711], the availability of an intrinsic isotope effect for an enzymatic reaction permits calculation of microscopic constants from steady-state data. By employment of a minimal mechanism for bovine plasma amine oxidase involving a single precatalytic and multiple postcatalytic enzyme-substrate complexes, equations have been derived that allow calculation of k3 and K1 when DKeq congruent to 1 less than Dk. Unexpectedly, in the case of K1, we have shown that this parameter can be calculated from steady-state parameters without the requirement for an intrinsic isotope effect. This result should have general application to both ping-pong and sequential ternary-complex enzyme mechanisms. Of significance for future applications of steady-state isotope effects to the calculation of microscopic constants, values for K1 and k3 derived from steady-state parameters and single turnover measurements indicate excellent agreement. Compilation of parameters among six ring-substituted phenethylamines reveals alteration in delta G for enzyme-substrate complex formation by 2.8 kcal/mol, together with an essentially invariant rate constant for C-H bond activation. A detailed discussion of the relevance of these findings to the interrelationship of binding energy and catalytic efficiency in enzyme reactions is presented.     

Mycobacterium tuberculosis beta-ketoacyl-acyl carrier protein (ACP) reductase: kinetic and chemical mechanisms

Beta-ketoacyl-acyl carrier protein (ACP) reductase from Mycobacterium tuberculosis (MabA) is responsible for the second step of the type-II fatty acid elongation system of bacteria, plants, and apicomplexan organisms, catalyzing the NADPH-dependent reduction of beta-ketoacyl-ACP to generate beta-hydroxyacyl-ACP and NADP(+). In the present work, the mabA-encoded MabA has been cloned, expressed, and purified to homogeneity. Initial velocity studies, product inhibition, and primary deuterium kinetic isotope effects suggested a steady-state random bi-bi kinetic mechanism for the MabA-catalyzed reaction. The magnitudes of the primary deuterium kinetic isotope effect indicated that the C(4)-proS hydrogen is transferred from the pyridine nucleotide and that this transfer contributes modestly to the rate-limiting step of the reaction. The pH-rate profiles demonstrated groups with pK values of 6.9 and 8.0, important for binding of NADPH, and with pK values of 8.8 and 9.6, important for binding of AcAcCoA and for catalysis, respectively. Temperature studies were employed to determine the activation energy of the reaction. Solvent kinetic isotope effects and proton inventory analysis established that a single proton is transferred in a partially rate-limiting step and that the mechanism of carbonyl reduction is probably concerted. The observation of an inverse (D)2(O)V/K and an increase in (D)2(O)V when [4S-(2)H]NADPH was the varied substrate obscured the distinction between stepwise and concerted mechanisms; however, the latter was further supported by the pH dependence of the primary deuterium kinetic isotope effect. Kinetic and chemical mechanisms for the MabA-catalyzed reaction are proposed on the basis of the experimental data.     

Electronegativity of aromatic amines as a basis for the development of ground state inhibitors of lysyl oxidase

Benzylamine derivatives containing para substituents of differing electronegativity as well as isomers of aminomethylpyridine have been assessed for their substrate and inhibitor potentials toward lysyl oxidase. Substituted benzylamines with increasingly electronegative para substituents had the lowest KI values and thus were the most effective inhibitors of the oxidation of elastin by lysyl oxidase. The kcat values for these compounds as substrates of lysyl oxidase were also reduced with increasingly electronegative para substituents. Both the Dkcat and D(kcat/Km) kinetic isotope effects decreased with increasingly electronegative p-substituents in [alpha, alpha'-2H]benzylamines. In contrast, there was no Dkcat solvent isotope effect with [2H] H2O while the D(kcat/Km) solvent isotope effect tended to increase with increasingly electronegative p-substituents. These results are consistent with the stabilization of an enzyme-generated substrate carbanion and thus the retardation of substrate oxidation by electronegative substituents. Such ground state stabilization thus can result in compounds with increased potential for the inhibition of the oxidation of protein substrates of lysyl oxidase.     

Deuterium isotope effect of proton pumping in cytochrome c oxidase

In mitochondria and many aerobic bacteria cytochrome c oxidase is the terminal enzyme of the respiratory chain where it catalyses the reduction of oxygen to water. The free energy released in this process is used to translocate (pump) protons across the membrane such that each electron transfer to the catalytic site is accompanied by proton pumping. To investigate the mechanism of electron-proton coupling in cytochrome c oxidase we have studied the pH-dependence of the kinetic deuterium isotope effect of specific reaction steps associated with proton transfer in wild-type and structural variants of cytochrome c oxidases in which amino-acid residues in proton-transfer pathways have been modified. In addition, we have solved the structure of one of these mutant enzymes, where a key component of the proton-transfer machinery, Glu286, was modified to an Asp. The results indicate that the P3-->F3 transition rate is determined by a direct proton-transfer event to the catalytic site. In contrast, the rate of the F3-->O4 transition, which involves simultaneous electron transfer to the catalytic site and is characteristic of any transition during CytcO turnover, is determined by two events with similar rates and different kinetic isotope effects. These reaction steps involve transfer of protons, that are pumped, via a segment of the protein including Glu286 and Arg481.     

13C isotope effects as a probe of the kinetic mechanism and allosteric properties of Escherichia coli aspartate transcarbamylase.

13C kinetic isotope effects have been measured in carbamyl phosphate for the reaction catalyzed by aspartate transcarbamylase. For the holoenzyme, the value was 1.0217 at zero aspartate, but unity at infinite aspartate, with 4.8 mM aspartate eliminating half of the isotope effect. This pattern proves an ordered kinetic mechanism, with carbamyl phosphate adding before aspartate. The same parameters were observed in the presence of ATP or CTP, showing that there is only one form of active enzyme present, regardless of the presence or absence of allosteric modifiers. These data support the Monod model of allosteric behavior in which the equilibrium between fully active and inactive enzyme is perturbed by selective binding interactions of substrates and modifiers, and there are no enzyme forms having partial activity. Isolated catalytic subunits of the enzyme showed similar 13C isotope effects (1.0240 at zero aspartate, 1.0039 at infinite aspartate, 3.8 mM aspartate causing half of the change from one value to the other), but the finite isotope effect at infinite aspartate shows that the kinetic mechanism is now partly random. With the very slow and poorly bound aspartate analog cysteine sulfinate, the 13C isotope effects were 1.039 for both holoenzyme and catalytic subunits and were not decreased significantly by high levels of cysteine sulfinate. The value of 1.039 is probably close to the intrinsic isotope effect on the chemical reaction, while the kinetic mechanism with this substrate is now fully random because the chemistry is so much slower than release of either reactant from the enzyme.     

Evidence that both protium and deuterium undergo significant tunneling in the reaction catalyzed by bovine serum amine oxidase

The magnitudes of primary and secondary H/T and D/T kinetic isotope effects have been measured in the bovine serum amine oxidase catalyzed oxidation of benzylamine from 0 to 45 degrees C. Secondary H/T and D/T kinetic effects are small and in the range anticipated from equilibrium isotope effects; Arrhenius preexponential factors (AH/AT and AD/AT) determined from the temperature dependence of isotope effects also indicate semiclassical behavior. By contrast, primary H/T and D/T isotope effects, 35.2 +/- 0.8 and 3.07 +/- 0.07, respectively, at 25 degrees C, are larger than semiclassical values and give anomalously low preexponential factor ratios, AH/AT = 0.12 +/- 0.04 and AD/AT = 0.51 +/- 0.10. Stopped-flow studies indicate similar isotope effects on cofactor reduction as seen in the steady state, consistent with a single rate-limiting C-H bond cleavage step for Vmax/Km. The comparison of primary and secondary isotope effects allows us to rule out appreciable coupling between the primary and secondary hydrogens at C-1 of the substrate. From the properties of primary isotope effects, we conclude that both protium and deuterium undergo significant tunneling in the course of substrate oxidation. These findings represent the first example of quantum mechanical effects in an enzyme-catalyzed proton abstraction reaction.     

Computer simulation of primary kinetic isotope effects in the proposed rate-limiting step of the glyoxalase I catalyzed reaction

The proposed rate-limiting step of the glyoxalase I catalyzed reaction is the proton abstraction from the C1 carbon of the substrate by Glu(172). Here we examine primary kinetic isotope effects and the influence of quantum dynamics on this process by computer simulations. The calculations utilize the empirical valence bond method in combination with the molecular dynamics free energy perturbation technique and path integral simulations. For the enzyme-catalyzed reaction a H/D kinetic isotope effect of 5.0 +/- 1. 3 is predicted in reasonable agreement with the experimental result of about 3. Furthermore, the magnitude of quantum mechanical effects is found to be very similar for the enzyme reaction and the corresponding uncatalyzed process in solution, in agreement with other studies. The problems associated with attaining the required accuracy in order for the present approach to be useful as a diagnostic tool for the study of enzyme reactions are also discussed.