Stereoselective preparation of deuterated reduced nicotinamide adenine nucleotides and substrates by enzymatic synthesis

Methods for biochemical characterization of flavin-dependent N-monooxygenases involved in siderophore biosynthesis

Siderophores are essential molecules released by some bacteria and fungi in iron-limiting environments to sequester ferric iron, satisfying metabolic needs. Flavin-dependent N-hydroxylating monooxygenases (NMOs) catalyze the hydroxylation of nitrogen atoms to generate important siderophore functional groups such as hydroxamates. It has been demonstrated that the function of NMOs is essential for virulence, implicating these enzymes as potential drug targets. This chapter aims to serve as a resource for the characterization of NMO's enzymatic activities using several biochemical techniques. We describe assays that allow for the determination of steady-state kinetic parameters, detection of hydroxylated amine products, measurement of the rate-limiting step(s), and the application toward drug discovery efforts. While not exhaustive, this chapter will provide a foundation for the characterization of enzymes involved in siderophore biosynthesis, allowing for gaps in knowledge within the field to be addressed.

Photoinduced monooxygenation involving NAD(P)H-FAD sequential single-electron transfer

Light-dependent or light-stimulated catalysis provides a multitude of perspectives for implementation in technological or biomedical applications. Despite substantial progress made in the field of photobiocatalysis, the number of usable light-responsive enzymes is still very limited. Flavoproteins have exceptional potential for photocatalytic applications because the name-giving cofactor intrinsically features light-dependent reactivity, undergoing photoreduction with a variety of organic electron donors. However, in the vast majority of these enzymes, photoreactivity of the enzyme-bound flavin is limited or even suppressed. Here, we present a flavoprotein monooxygenase in which catalytic activity is controllable by blue light illumination. The reaction depends on the presence of nicotinamide nucleotide-type electron donors, which do not support the reaction in the absence of light. Employing various experimental approaches, we demonstrate that catalysis depends on a protein-mediated photoreduction of the flavin cofactor, which proceeds via a radical mechanism and a transient semiquinone intermediate.

The biosynthetic origin of psychoactive kavalactones in kava

Kava (Piper methysticum) is an ethnomedicinal shrub native to the Polynesian islands with well-established anxiolytic and analgesic properties. Its main psychoactive principles, kavalactones, form a unique class of polyketides that interact with the human central nervous system through mechanisms distinct from those of conventional psychiatric drugs. However, an unknown biosynthetic machinery and difficulty in chemical synthesis hinder the therapeutic use of kavalactones. In addition, kava also produces flavokavains, which are chalconoids with anticancer properties structurally related to kavalactones. Here, we report de novo elucidation of the key enzymes of the kavalactone and flavokavain biosynthetic network. We present the structural basis for the evolutionary development of a pair of paralogous styrylpyrone synthases that establish the kavalactone scaffold and the catalytic mechanism of a regio- and stereo-specific kavalactone reductase that produces a subset of chiral kavalactones. We further demonstrate the feasibility of engineering styrylpyrone production in heterologous hosts, thus opening a way to develop kavalactone-based non-addictive psychiatric therapeutics through synthetic biology.

Halophilic mechanism of the enzymatic function of a moderately halophilic dihydrofolate reductase from Haloarcula japonica strain TR-1

Dihydrofolate (DHF) reductase coded by a plasmid of the extremely halophilic archaeon Haloarcula japonica strain TR-1 (HjDHFR P1) shows moderate halophilicity on enzymatic activity at pH 6.0, although there is no significant effect of NaCl on its secondary structure. To elucidate the salt-activation and -inactivation mechanisms of this enzyme, we investigated the effects of pH and salt concentration, deuterium isotope effect, steady-state kinetics, and rapid-phase ligand-binding kinetics. Enzyme activity was increased eightfold by the addition of 500 mM NaCl at pH 6.0, fourfold by 250 mM at pH 8.0, and became independent of salt concentration at pH 10.0. Full isotope effects observed at pH 10.0 under 0-1000 mM NaCl indicated that the rate of hydride transfer, which was the rate-determining step at the basic pH region, was independent of salt concentration. Conversely, rapid-phase ligand-binding experiments showed that the amplitude of the DHF-binding reaction increased and the tetrahydrofolate (THF)-releasing rate decreased with increasing NaCl concentration. These results suggested that the salt-activation mechanism of HjDHFR P1 is via the population change of the anion-unbound and anion-bound conformers, which are binding-incompetent and -competent conformations for DHF, respectively, while that of salt inactivation is via deceleration of the THF-releasing rate, which is the rate-determining step at the neutral pH region.

Kinetics of NADP+/NADPH reduction-oxidation catalyzed by the ferredoxin-NAD(P)+ reductase from the green sulfur bacterium Chlorobaculum tepidum

Ferredoxin-NAD(P)⁺ oxidoreductase (FNR, [EC 1.18.1.2], [EC 1.18.1.3]) from the green sulfur bacterium Chlorobaculum tepidum (CtFNR) is a homodimeric flavoprotein with significant structural homology to bacterial NADPH-thioredoxin reductases. CtFNR homologs have been found in many bacteria, but only in green sulfur bacteria among photoautotrophs. In this work, we examined the reactions of CtFNR with NADP⁺, NADPH, and (4S-²H)-NADPD by stopped-flow spectrophotometry. Mixing CtFNR_ox with NADPH yielded a rapid decrease of the absorbance in flavin band I centered at 460 nm within 1 ms, and then the absorbance further decreased gradually. The magnitude of the decrease increased with increasing NADPH concentration, but even with ~50-fold molar excess NADPH, the absorbance change was only ~45 % of that expected for fully reduced protein. The absorbance in the charge transfer (CT) band centered around 600 nm increased rapidly within 1 ms, then slowly decreased to about 70 % of the maximum. When CtFNR_red was mixed with excess NADP⁺, the absorbance in the flavin band I increased to about 70 % of that of CtFNR_ox with an apparent rate of ~4 s^-1, whereas almost no absorption changes were observed in the CT band. Obtained data suggest that the reaction between CtFNR and NADP⁺/NADPH is reversible, in accordance with its physiological function.

Role of active site loop in coenzyme binding and flavin reduction in cytochrome P450 reductase

Cytochrome P450 reductase (CPR) contains a loop within the active site (comprising Asp(634), Ala(635), Arg(636) and Asn(637); human CPR numbering) that relocates upon NADPH binding. Repositioning of the loop triggers the reorientation of an FAD-shielding tryptophan (Trp(679)) to a partially stacked conformer, reducing the energy barrier for displacement of the residue by the NADPH nicotinamide ring: an essential step for hydride transfer. We used site-directed mutagenesis and kinetic analysis to investigate if the amino acid composition of the loop influences the catalytic properties of CPR. The D634A and D634N variants elicited a modest increase in coenzyme binding affinity coupled with a 36- and 10-fold reduction in cytochrome c(3+) turnover and a 17- and 3-fold decrease in the pre-steady state rate of flavin reduction. These results, in combination with a reduction in the kinetic isotope effect for hydride transfer, suggest that diminished activity is due to destabilization of the partially stacked conformer of Trp(677) and slower release of NADP(+). In contrast, R636A, R636S and an A635G/R636S double mutant led to a modest increase in cytochrome c(3+) reduction, which is linked to weaker coenzyme binding and faster interflavin electron transfer. A potential mechanism by which Arg(636) influences catalysis is discussed.

Pre-steady-state kinetic studies of redox reactions catalysed by Bacillus subtilis ferredoxin-NADP(+) oxidoreductase with NADP(+)/NADPH and ferredoxin

Ferredoxin-NADP(+) oxidoreductase ([EC1.18.1.2], FNR) from Bacillus subtilis (BsFNR) is a homodimeric flavoprotein sharing structural homology with bacterial NADPH-thioredoxin reductase. Pre-steady-state kinetics of the reactions of BsFNR with NADP(+), NADPH, NADPD (deuterated form) and B. subtilis ferredoxin (BsFd) using stopped-flow spectrophotometry were studied. Mixing BsFNR with NADP(+) and NADPH yielded two types of charge-transfer (CT) complexes, oxidized FNR (FNR(ox))-NADPH and reduced FNR (FNR(red))-NADP(+), both having CT absorption bands centered at approximately 600n m. After mixing BsFNR(ox) with about a 10-fold molar excess of NADPH (forward reaction), BsFNR was almost completely reduced at equilibrium. When BsFNR(red) was mixed with NADP(+), the amount of BsFNR(ox) increased with increasing NADP(+) concentration, but BsFNR(red) remained as the major species at equilibrium even with about 50-fold molar excess NADP(+). In both directions, the hydride-transfer was the rate-determining step, where the forward direction rate constant (~500 s(-1)) was much higher than the reverse one (<10 s(-1)). Mixing BsFd(red) with BsFNR(ox) induced rapid formation of a neutral semiquinone form. This process was almost completed within 1 ms. Subsequently the neutral semiquinone form was reduced to the hydroquinone form with an apparent rate constant of 50 to 70 s(-1) at 10°C, which increased as BsFd(red) increased from 40 to 120 μM. The reduction rate of BsFNR(ox) by BsFd(red) was markedly decreased by premixing BsFNR(ox) with BsFd(ox), indicating that the dissociation of BsFd(ox) from BsFNR(sq) is rate-limiting in the reaction. The characteristics of the BsFNR reactions with NADP(+)/NADPH were compared with those of other types of FNRs.

Probing the chemical mechanism of saccharopine reductase from Saccharomyces cerevisiae using site-directed mutagenesis

Saccharopine reductase catalyzes the reductive amination of l-α-aminoadipate-δ-semialdehyde with l-glutamate to give saccharopine. Two mechanisms have been proposed for the reductase, one that makes use of enzyme side chains as acid-base catalytic groups, and a second, in which the reaction is catalyzed by enzyme-bound reactants. Site-directed mutagenesis was used to change acid-base candidates in the active site of the reductase to eliminate their ionizable side chain. Thus, the D126A, C154S and Y99F and several double mutant enzymes were prepared. Kinetic parameters in the direction of glutamate formation exhibited modest decreases, inconsistent with the loss of an acid-base catalyst. The pH-rate profiles obtained with all mutant enzymes decrease at low and high pH, suggesting acid and base catalytic groups are still present in all enzymes. Solvent kinetic deuterium isotope effects are all larger than those observed for wild type enzyme, and approximately equal to one another, suggesting the slow step is the same as that of wild type enzyme, a conformational change to open the site and release products (in the direction of saccharopine formation). Overall, the acid-base chemistry is likely catalyzed by bound reactants, with the exception of deprotonation of the α-amine of glutamate, which likely requires an enzyme residue.

Practical aspects on the use of kinetic isotope effects as probes of flavoprotein enzyme mechanisms

The measurement of kinetic isotope effects (KIEs) has proved useful in many mechanistic studies of enzyme activity, most notably in enzyme-catalyzed hydrogen-transfer reactions. Primary KIEs (1° KIE) greater than unity indicate that transfer of the hydrogen species of interest is partially or fully rate limiting, and studies of the magnitude of the temperature and pressure dependence of these KIEs can inform on the mechanism of transfer. For example, KIE measurements have proved crucial in understanding the role of quantum mechanical tunneling in enzyme systems. The measurement of secondary KIEs (2° KIEs) is also informative and can be used to infer a significant tunneling contribution and details of transition state geometry. Here the deuterium label is introduced next to that of the transferred hydrogen. Measurements of 1° and 2° KIEs are being used increasingly in studies of H-transfer in flavoprotein enzymes and this requires the preparation of high purity and stereospecific labeled isotopologues. Strategies for the synthesis of labeled substrates are dependent on the enzyme system being studied. However, the nicotinamide coenzymes are often used in studies of flavoprotein enzyme mechanisms. Here, we provide practical details for the enzymatic synthesis of high purity deuterated isotopologues of the common biological coenzymes NADH and NADPH as well as the corresponding nonreactive mimics, tetrahydroNAD(P)H. Both forms of the coenzyme have proven useful in the study of mechanisms, particularly in relation to the involvement of quantum mechanical tunneling and dynamics in enzymatic H-transfer chemistry. The focus here is on practical considerations in the synthesis of these compounds. We also provide an abbreviated description of how measurements of KIEs can inform on flavoprotein mechanisms. The aim of this contribution is not to give a detailed description of the underlying theory (which has been reviewed extensively in the literature), but to provide a basic introduction and practical considerations for nonexpert readers who wish to incorporate such measurements in studies of enzyme mechanisms.

Artificial metalloenzymes for the diastereoselective reduction of NAD(+) to NAD(2)H

Stereoselectively labelled isotopomers of NAD(P)H are highly relevant for mechanistic studies of enzymes which utilize them as redox equivalents. Whereas several methods are firmly established for their generation in high diastereomeric purity by enzymatic methods, alternative methods have so far not been investigated. The article presents the stereoselective deuteration of NAD(+) at the 4-position (90% de) of the pyridinium-ring by means of an artificial metalloenzyme. The artificial metalloenzyme consists of a biotinylated iridium cofactor embedded in streptavidin isoforms and the resulting constructs have been previously shown to be compatible with natural enzymes. Alternative methods for stereoselective NAD(P)(+) reduction are expected to be of high interest for the mechanistic study of enzymes that accept NAD(P)H mimics and for the synthesis of structurally related fine chemicals.

Alteration of the flexible loop in 1-deoxy-D-xylulose-5-phosphate reductoisomerase boosts enthalpy-driven inhibition by fosmidomycin

1-Deoxy-d-xylulose-5-phosphate reductoisomerase (DXR), which catalyzes the first committed step in the 2-C-methyl-d-erythritol 4-phosphate pathway of isoprenoid biosynthesis used by Mycobacterium tuberculosis and other infectious microorganisms, is absent in humans and therefore an attractive drug target. Fosmidomycin is a nanomolar inhibitor of DXR, but despite great efforts, few analogues with comparable potency have been developed. DXR contains a strictly conserved residue, Trp203, within a flexible loop that closes over and interacts with the bound inhibitor. We report that while mutation to Ala or Gly abolishes activity, mutation to Phe and Tyr only modestly impacts k_cat and K_m. Moreover, pre-steady-state kinetics and primary deuterium kinetic isotope effects indicate that while turnover is largely limited by product release for the wild-type enzyme, chemistry is significantly more rate-limiting for W203F and W203Y. Surprisingly, these mutants are more sensitive to inhibition by fosmidomycin, resulting in K_m/K_i ratios up to 19-fold higher than that of wild-type DXR. In agreement, isothermal titration calorimetry revealed that fosmidomycin binds up to 11-fold more tightly to these mutants. Most strikingly, mutation strongly tips the entropy–enthalpy balance of total binding energy from 50% to 75% and 91% enthalpy in W203F and W203Y, respectively. X-ray crystal structures suggest that these enthalpy differences may be linked to differences in hydrogen bond interactions involving a water network connecting fosmidomycin’s phosphonate group to the protein. These results confirm the importance of the flexible loop, in particular Trp203, in ligand binding and suggest that improved inhibitor affinity may be obtained against the wild-type protein by introducing interactions with this loop and/or the surrounding structured water network.

Mechanistic studies on the flavin-dependent N⁶-lysine monooxygenase MbsG reveal an unusual control for catalysis

The mechanism of Mycobacterium smegmatis G (MbsG), a flavin-dependent l-lysine monooxygenase, was investigated under steady-state and rapid reaction conditions using primary and solvent kinetic isotope effects, substrate analogs, pH and solvent viscosity effects as mechanistic probes. The results suggest that l-lysine binds before NAD(P)H, which leads to a decrease in the rate constant for flavin reduction. l-lysine binding has no effect on the rate of flavin oxidation, which occurs in a one-step process without the observation of a C4a-hydroperoxyflavin intermediate. Similar effects were determined with several substrate analogs. Flavin oxidation is pH independent while the kcat/Km and kred/KD pH profiles for NAD(P)H exhibit single pKa values of ∼6.0, with increasing activity as the pH decreases. At lower pH, the enzyme becomes more uncoupled, producing more hydrogen peroxide and superoxide. Hydride transfer is partially rate-limiting at neutral pH and becomes more rate-limiting at low pH. An inverse solvent viscosity effect on kcat/Km for NAD(P)H was observed at neutral pH whereas a normal solvent viscosity effect was observed at lower pH. Together, the results indicate a unique mechanism where a rate-limiting and pH-sensitive conformational change occurs in the reductive half-reaction, which affects the efficiency of lysine hydroxylation.

Proximal FAD histidine residue influences interflavin electron transfer in cytochrome P450 reductase and methionine synthase reductase

Cytochrome P450 reductase (CPR) and methionine synthase reductase (MSR) transfer reducing equivalents from NADPH to FAD to FMN. In CPR, hydride transfer and interflavin electron transfer are kinetically coupled steps, but in MSR the two catalytic steps are represented by two distinct kinetic phases leading to transient formation of the FAD hydroquinone. In human CPR, His(322) forms a hydrogen-bond with the highly conserved Asp(677), a member of the catalytic triad. The catalytic triad is present in MSR, but Ala(312) replaces the histidine residue. To examine if this structural variation accounts for differences in their kinetic behavior, reciprocal substitutions were created. Substitution of His(322) for Ala in CPR does not affect the rate of NADPH hydride transfer or the FAD redox potentials, but does impede interflavin electron transfer. For MSR, swapping Ala(312) for a histidine residue resulted in the kinetic coupling of hydride and interflavin electron transfer, and eliminated the formation of the FAD hydroquinone intermediate. For both enzymes, placement of the His residue in the active site weakens coenzyme binding affinity. The data suggest that the proximal FAD histidine residue accelerates proton-coupled electron transfer from FADH2 to the higher potential FMN; a mechanism for this catalytic role is discussed.

Kinetic analysis of cytochrome P450 reductase from Artemisia annua reveals accelerated rates of NADH-dependent flavin reduction

Cytochrome P450 reductase from Artemisia annua (aaCPR) is a diflavin enzyme that has been employed for the microbial synthesis of artemisinic acid (a semi-synthetic precursor of the anti-malarial drug, artemisinin) based on its ability to transfer electrons to the cytochrome P450 monooxygenase, CYP71AV1. We have isolated recombinant aaCPR (with the N-terminal transmembrane motif removed) from Escherichia coli and compared its kinetic and thermodynamic properties with other CPR orthologues, most notably human CPR. The FAD and FMN redox potentials and the macroscopic kinetic constants associated with cytochrome c(3+) reduction for aaCPR are comparable to that of other CPR orthologues, with the exception that the apparent binding affinity for the oxidized coenzyme is ~ 30-fold weaker compared to human CPR. CPR from A. annua shows a 3.5-fold increase in uncoupled NADPH oxidation compared to human CPR and a strong preference (85 100-fold) for NADPH over NADH. Strikingly, reduction of the enzyme by the first and second equivalent of NADPH is much faster in aaCPR, with rates of > 500 and 17 s(-1) at 6 °C. We also optically detect a charge-transfer species that rapidly forms in < 3 ms and then persists during the reductive half reaction. Additional stopped-flow kinetic studies with NADH and (R)-[4-(2) H]NADPH suggest that the accelerated rate of flavin reduction is attributed to the relatively weak binding affinity of aaCPR for NADP(+) .

Aromatic substitution of the FAD-shielding tryptophan reveals its differential role in regulating electron flux in methionine synthase reductase and cytochrome P450 reductase

Methionine synthase reductase (MSR) and cytochrome P450 reductase (CPR) transfer reducing equivalents from NADPH via an FAD and FMN cofactor to a redox partner protein. In both enzymes, hydride transfer from NADPH to FAD requires displacement of a conserved tryptophan that lies coplanar to the FAD isoalloxazine ring. Swapping the tryptophan for a smaller aromatic side chain revealed a distinct role for the residue in regulating MSR and CPR catalysis. MSR W697F and W697Y showed enhanced catalysis, noted by increases in kcat and k(cat)/K(m)(NADPH) for steady-state cytochrome c(3+) reduction and a 10-fold increase in the rate constant (k(obs1)) associated with hydride transfer. Elevated primary kinetic isotope effects on k(obs1) for W697F and W697Y suggest that preceding isotopically insensitive steps like displacement of W697 are less rate determining. MSR W697Y, but not MSR W697F, showed detectable formation of the disemiquinone intermediate, indicating that the polarity of the aromatic side chain influences the rate of interflavin electron transfer. By contrast, the CPR variants (W676F and W676Y) displayed modest decreases in cytochrome c(3+) reduction, a 30- and 3.5-fold decrease in the rate of FAD reduction, accumulation of a FADH2 -NADP(+) charge-transfer complex and dramatically suppressed rates of interflavin electron transfer. We conclude for MSR that hydride transfer is 'gated' by the free energy required to disrupt dispersion forces between the FAD isoalloxazine ring and W697. By contrast, the bulky indole ring of W676 accelerates catalysis in CPR by lowering the energy barrier for displacement of the oxidized nicotinamide ring coplanar with the FAD.

Stereochemical outcome at four stereogenic centers during conversion of prephenate to tetrahydrotyrosine by BacABGF in the bacilysin pathway

The first four enzymes of the bacilysin antibiotic pathway, BacABGF, convert prephenate to a tetrahydrotyrosine (H(4)Tyr) diastereomer on the way to the anticapsin warhead of the dipeptide antibiotic. BacB takes the BacA product endocyclic-Δ(4),Δ(8)-7R-dihydrohydroxyphenylpyruvate (en-H(2)HPP) and generates a mixture of 3E- and 3Z-olefins of the exocyclic-Δ(3),Δ(5)-dihydrohydroxyphenylpyruvate (ex-H(2)HPP). The NADH-utilizing BacG then catalyzes a conjugate reduction, adding a pro-S hydride equivalent to C(4) to yield tetrahydrohydroxyphenylpyruvate (H(4)HPP), a transamination away (via BacF) from 2S-H(4)Tyr. Incubations of the pathway enzymes in D(2)O yield deuterium incorporation at C(8) from BacA and then C(9) from BacB action. By (1)H NMR analysis of samples of H(4)Tyr, the stereochemistry at C(4), C(8), and C(9) can be assigned. BacG (followed by BacF) converts 3E-ex-H(2)HPP to 2S,4R,7R-H(4)Tyr. The 3Z isomer is instead reduced and transaminated to the opposite diastereomer at C(4), 2S,4S,7R-H(4)Tyr. Given that bacilysin has the 2S,4S stereochemistry in its anticapsin moiety, it is likely that the 2S,4S-H(4)Tyr is the diastereomer "on pathway". NMR determination of the stereochemistry of the CHD samples at C(8) and C(9) allows assignment of all stereogenic centers (except C(3)) in this unusual tetrahydro-aromatic amino acid building block, giving insights into and constraints on the BacA, BacB, and BacG mechanisms.

Supporting role of lysine 13 and glutamate 16 in the acid-base mechanism of saccharopine dehydrogenase from Saccharomyces cerevisiae

Saccharopine dehydrogenase (SDH) catalyzes the NAD+ dependent oxidative deamination of saccharopine to form lysine (Lys) and α-ketoglutarate (α-kg). The active site of SDH has a number of conserved residues that are believed important to the overall reaction. Lysine 13, positioned near the active site base (K77), forms a hydrogen bond to E78 neutralizing it, and contributing to setting the pKa of the catalytic residues to near neutral pH. Glutamate 16 is within hydrogen bond distance to the Nε atom of R18, which has strong H-bonding interactions with the α-carboxylate and α-oxo groups of α-kg. Mutation of K13 to M and E16 to Q decreased kcat by about 15-fold, and primary and solvent deuterium kinetic isotope effects measured with the mutant enzymes indicate hydride transfer is rate limiting for the overall reaction. The pH-rate profiles for K13M exhibited no pH dependence, consistent with an increase in negative charge in the active site resulting in the perturbation in the pKas of catalytic groups. Elimination of E16 affects optimal positioning of R18, which is involved in binding and holding α-kg in the correct conformation for optimum catalysis. In agreement, a ΔΔG°' of 2.60 kcal/mol is estimated from the change in Kα-kg for replacing E16 with Q.

Evidence in support of lysine 77 and histidine 96 as acid-base catalytic residues in saccharopine dehydrogenase from Saccharomyces cerevisiae

Saccharopine dehydrogenase (SDH) catalyzes the final reaction in the α-aminoadipate pathway, the conversion of l-saccharopine to l-lysine (Lys) and α-ketoglutarate (α-kg) using NAD⁺ as an oxidant. The enzyme utilizes a general acid-base mechanism to conduct its reaction with a base proposed to accept a proton from the secondary amine of saccharopine in the oxidation step and a group proposed to activate water to hydrolyze the resulting imine. Crystal structures of an open apo form and a closed form of the enzyme with saccharopine and NADH bound have been determined at 2.0 and 2.2 Å resolution, respectively. In the ternary complex, a significant movement of domain I relative to domain II that closes the active site cleft between the two domains and brings H96 and K77 into the proximity of the substrate binding site is observed. The hydride transfer distance is 3.6 Å, and the side chains of H96 and K77 are properly positioned to act as acid-base catalysts. Preparation of the K77M and H96Q single-mutant and K77M/H96Q double-mutant enzymes provides data consistent with their role as the general acid-base catalysts in the SDH reaction. The side chain of K77 initially accepts a proton from the ε-amine of the substrate Lys and eventually donates it to the imino nitrogen as it is reduced to a secondary amine in the hydride transfer step, and H96 protonates the carbonyl oxygen as the carbinolamine is formed. The K77M, H976Q, and K77M/H96Q mutant enzymes give 145-, 28-, and 700-fold decreases in V/E(t) and >10³-fold increases in V₂/K(Lys)E(t) and V₂/K(α-kg)E(t) (the double mutation gives >10⁵-fold decreases in the second-order rate constants). In addition, the K77M mutant enzyme exhibits a primary deuterium kinetic isotope effect of 2.0 and an inverse solvent deuterium isotope effect of 0.77 on V₂/K(Lys). A value of 2.0 was also observed for (D)(V₂/K(Lys))(D₂O) when the primary deuterium kinetic isotope effect was repeated in D₂O, consistent with a rate-limiting hydride transfer step. A viscosity effect of 0.8 was observed on V₂/K(Lys), indicating the solvent deuterium isotope effect resulted from stabilization of an enzyme form prior to hydride transfer. A small normal solvent isotope effect is observed on V, which decreases slightly when repeated with NADD, consistent with a contribution from product release to rate limitation. In addition, V₂/K(Lys)E(t) is pH-independent, which is consistent with the loss of an acid-base catalyst and perturbation of the pK(a) of the second catalytic group to a higher pH, likely a result of a change in the overall charge of the active site. The primary deuterium kinetic isotope effect for H96Q, measured in H₂O or D₂O, is within error equal to 1. A solvent deuterium isotope effect of 2.4 is observed with NADH or NADD as the dinucleotide substrate. Data suggest rate-limiting imine formation, consistent with the proposed role of H96 in protonating the leaving hydroxyl as the imine is formed. The pH-rate profile for V₂/K(Lys)E(t) exhibits the pK(a) for K77, perturbed to a value of ∼9, which must be unprotonated to accept a proton from the ε-amine of the substrate Lys so that it can act as a nucleophile. Overall, data are consistent with a role for K77 acting as the base that accepts a proton from the ε-amine of the substrate lysine prior to nucleophilic attack on the α-oxo group of α-ketoglutarate, and finally donating a proton to the imine nitrogen as it is reduced to give saccharopine. In addition, data indicate a role for H96 acting as a general acid-base catalyst in the formation of the imine between the ε-amine of lysine and the α-oxo group of α-ketoglutarate.

Contribution of K99 and D319 to substrate binding and catalysis in the saccharopine dehydrogenase reaction

Saccharopine dehydrogenase catalyzes the NAD-dependent oxidative deamination of saccharopine to l-lysine and α-ketoglutarate. Lysine 99 is within hydrogen-bond distance to the α-carboxylate of the lysine substrate and D319 is in the vicinity of the carboxamide side chain of NADH. Both are conserved and may be important to the overall reaction. Replacing K99 with M gives decreases of 110-, 80- and 20-fold in the V(2)/K(m) values for lysine, α-ketoglutarate and NADH, respectively. Deuterium isotope effects on V and V/K(Lys) increase, while the solvent deuterium isotope effects decrease compared to the C205S mutant enzyme. Data for K99M suggest a decreased affinity for all reactants and a change in the partition ratio of the imine intermediate to favor hydrolysis. A change in the bound conformation of the imine and/or the distance of the imine carbon to C4 of the nicotinamide ring of NADH is also suggested. Changing D319 to A decreases V(2)/K(NADH) by 33-fold. Primary deuterium and solvent deuterium isotope effects decrease compared to C205S suggesting a non-isotope sensitive step has become slower. NADH binds to enzyme first, and sets the site for binding of lysine and α-ketoglutarate. The slower step is likely the conformational change generated upon binding of NADH.

The oxidation state of active site thiols determines activity of saccharopine dehydrogenase at low pH

Saccharopine dehydrogenase catalyzes the NAD-dependent conversion of saccharopine to generate L-lysine and α-ketoglutarate. A disulfide bond between cysteine 205 and cysteine 249, in the vicinity of the dinucleotide-binding site, is observed in structures of the apoenzyme, while a dithiol is observed in a structure with AMP bound, suggesting preferential binding of the dinucleotide to reduced enzyme. Mutation of C205 to S gave increased values of V/E(t) and V/KE(t) at pH 7 compared to wild type. Primary deuterium and solvent deuterium kinetic isotope effects suggest the catalytic pathway, which includes the hydride transfer and hydrolysis steps, contributes more to rate limitation in C205S, but the rates of the two steps relative to one another remain the same. There is a large increase in the rate constants V₁/E(t) and V₁/K(NAD)Et at pH values below 7 compared to WT. Data indicate the low pH increase in activity results from a decreased sensitivity of the C205S mutant enzyme to the protonation state of an enzyme group with a pK(a) of about 7, likely responsible for a pH-dependent conformational change. Reduction of WT and C205S mutant enzymes with TCEP gives equal activities at pH 6, consistent with the increased activity observed for the C205S mutant enzyme.

Three serendipitous pathways in E. coli can bypass a block in pyridoxal-5'-phosphate synthesis

Overexpression of seven different genes restores growth of a ΔpdxB strain of E. coli, which cannot make pyridoxal phosphate (PLP), on M9/glucose.

None of the enzymes encoded by these genes has a promiscuous 4-phosphoerythronate dehydrogenase activity that can replace the activity of PdxB.

Overexpression of these genes restores PLP synthesis by three different serendipitous pathways that feed into the normal PLP synthesis pathway downstream of the blocked step.

Reactions in one of these pathways are catalyzed by low-level activities of enzymes of unknown function and a promiscuous activity of an enzyme that normally has a role in another pathway; one reaction appears to be non-enzymatic.

Most metabolic enzymes are prodigious catalysts that have evolved to accelerate chemical reactions with high efficiency and specificity. However, many enzymes have inefficient promiscuous activities, as well, as a result of the assemblage of highly reactive catalytic residues and cofactors in active sites. Although promiscuous activities are generally orders of magnitude less efficient than well-evolved activities (O'Brien and Herschlag, 1998, 2001; Wang et al, 2003; Taylor Ringia et al, 2004), they often enhance reaction rates by orders of magnitude relative to those of uncatalyzed reactions (O'Brien and Herschlag, 1998, 2001). Thus, promiscuous activities provide a reservoir of novel catalytic activities that can be recruited to serve new functions.

The evolutionary potential of promiscuous enzymes extends beyond the recruitment of single enzymes to serve new functions. Microbes contain hundreds of enzymes—E. coli contains about 1700 (Freilich et al, 2005)—raising the possibility that promiscuous enzymes can be patched together to generate ‘serendipitous' pathways that are not part of normal metabolism. We distinguish serendipitous pathways from latent or cryptic pathways, which are bona fide pathways involving dedicated enzymes that are produced only under particular environmental circumstances. In contrast, serendipitous pathways are patched together from enzymes that normally serve other functions and are not regulated in a coordinated manner in response to the need to synthesize or degrade a metabolite.

In this study, we describe the discovery of three serendipitous pathways that allow synthesis of pyridoxal phosphate (PLP) in a strain of E. coli that lacks 4-phosphoerythronate dehydrogenase (PdxB) when one of the seven different genes is overexpressed. These genes were identified in a multicopy suppression experiment in which a library of E. coli genes (from the ASKA collection) was introduced into a ΔpdxB strain of E. coli that is unable to synthesize PLP. Surprisingly, none of the enzymes encoded by these genes has a promiscuous 4-phosphoerythronate (4PE) dehydrogenase activity that can substitute for the missing PdxB. Rather, overproduction of these enzymes appears to facilitate at least three serendipitous pathways that draw material from other metabolic pathways and feed into the normal PLP synthesis pathway downstream of the blocked step (Figure 1).

We have characterized one of these pathways in detail (Figure 3). The first step, dephosphorylation of 3-phosphohydroxypyruvate, is catalyzed by YeaB, a predicted NUDIX hydrolase of unknown function. Although catalysis is inefficient (k_cat=5.7×10⁻⁵ s⁻¹ and k_cat/K_M>0.028 M⁻¹ s⁻¹), the enzymatic rate is 4×10⁷-fold faster than the rate of the uncatalyzed reaction, and is sufficient to support PLP synthesis when YeaB is overproduced. The second step in the pathway is decarboxylation of 3-hydroxypyruvate (3HP). Although we found two enzymes (1-deoxyxylulose-5-phosphate synthase and the catalytic domain of α-ketoglutarate dehydrogenase) that catalyze this reaction with low but respectable activity in vitro, their involvement in pathway 1 was ruled out by genetic methods. Surprisingly, the non-enzymatic rate of decarboxylation of 3HP appears to be sufficient to support PLP synthesis. The third step in the pathway, condensation of glycolaldehyde and glycine to form 4-hydroxy-L-threonine, is catalyzed by LtaE, a low-specificity threonine aldolase whose physiological role is not known. The final step, phosphorylation of 4-hydroxy-L-threonine, is catalyzed by homoserine kinase (ThrB), which is required for synthesis of threonine. The promiscuous phosphorylation of 4-hydroxy-L-threonine is 80-fold slower than the physiological phosphorylation of homoserine. The involvement of LtaE and ThrB in pathway 1 was confirmed by genetic experiments showing that overexpression of yeaB no longer restores growth of ΔpdxB strains lacking either ltaE or thrB.

Although pathway 1 is inefficient, it provides the ΔpdxB strain with the ability to grow under conditions in which survival is otherwise impossible. In general, serendipitous assembly of an inefficient pathway from promiscuous activities of available enzymes will be important whenever the pathway provides increased fitness. This might occur when a critical metabolite is no longer available from the environment, and survival depends on assembly of a new biosynthetic pathway. A second circumstance in which an inefficient serendipitous pathway might improve fitness is the appearance of a novel compound in the environment that can be exploited as a source of carbon, nitrogen or phosphorous. Finally, chemotherapeutic agents that block metabolic pathways in bacteria or cancer cells could provide selective pressure for assembly of serendipitous pathways that allow synthesis of the end product of the blocked pathway and thus a previously unappreciated source of drug resistance. In all of these cases, even an inefficient pathway can provide a selective advantage over other cells in a particular environmental niche, allowing survival and subsequent mutations that elevate the efficiency of the pathway.

Our work is consistent with the hypothesis that the recognized metabolic network of E. coli is underlain by a denser network of reactions due to promiscuous enzymes that use and generate recognized metabolites, but also unusual metabolites that normally have no physiological role. The findings reported here highlight the abundance of cryptic capabilities in the E. coli proteome that can be drawn on to generate novel pathways. Such pathways could provide a starting place for assembly of more efficient pathways, both in nature and in the hands of metabolic engineers.

Bacterial genomes encode hundreds to thousands of enzymes, most of which are specialized for particular functions. However, most enzymes have inefficient promiscuous activities, as well, that generally serve no purpose. Promiscuous reactions can be patched together to form multistep metabolic pathways. Mutations that increase expression or activity of enzymes in such serendipitous pathways can elevate flux through the pathway to a physiologically significant level. In this study, we describe the discovery of three serendipitous pathways that allow synthesis of pyridoxal-5′-phosphate (PLP) in a strain of E. coli that lacks 4-phosphoerythronate (4PE) dehydrogenase (PdxB) when one of seven different genes is overexpressed. We have characterized one of these pathways in detail. This pathway diverts material from serine biosynthesis and generates an intermediate in the normal PLP synthesis pathway downstream of the block caused by lack of PdxB. Steps in the pathway are catalyzed by a protein of unknown function, a broad-specificity enzyme whose physiological role is unknown, and a promiscuous activity of an enzyme that normally serves another function. One step in the pathway may be non-enzymatic.

Glutamates 78 and 122 in the active site of saccharopine dehydrogenase contribute to reactant binding and modulate the basicity of the acid-base catalysts

Saccharopine dehydrogenase catalyzes the NAD-dependent oxidative deamination of saccharopine to give l-lysine and alpha-ketoglutarate. There are a number of conserved hydrophilic, ionizable residues in the active site, all of which must be important to the overall reaction. In an attempt to determine the contribution to binding and rate enhancement of each of the residues in the active site, mutations at each residue are being made, and double mutants are being made to estimate the interrelationship between residues. Here, we report the effects of mutations of active site glutamate residues, Glu(78) and Glu(122), on reactant binding and catalysis. Site-directed mutagenesis was used to generate E78Q, E122Q, E78Q/E122Q, E78A, E122A, and E78A/E122A mutant enzymes. Mutation of these residues increases the positive charge of the active site and is expected to affect the pK(a) values of the catalytic groups. Each mutant enzyme was completely characterized with respect to its kinetic and chemical mechanism. The kinetic mechanism remains the same as that of wild type enzymes for all of the mutant enzymes, with the exception of E78A, which exhibits binding of alpha-ketoglutarate to E and E.NADH. Large changes in V/K(Lys), but not V, suggest that Glu(78) and Glu(122) contribute binding energy for lysine. Shifts of more than a pH unit to higher and lower pH of the pK(a) values observed in the V/K(Lys) pH-rate profile of the mutant enzymes suggests that the presence of Glu(78) and Glu(122) modulates the basicity of the catalytic groups.

Site-directed mutagenesis as a probe of the acid-base catalytic mechanism of homoisocitrate dehydrogenase from Saccharomyces cerevisiae

Homoisocitrate dehydrogenase (HIcDH) catalyzes the Mg2+- and K+-dependent oxidative decarboxylation of homoisocitrate to alpha-ketoadipate using NAD as the oxidant. A recent consideration of the structures of enzymes in the same family as HIcDH, including isopropylmalate and isocitrate dehydrogenases, suggests all of the family members utilize a Lys-Tyr pair to catalyze the acid-base chemistry of the reaction [Aktas, D. F., and Cook, P. F. (2009) Biochemistry 48, 3565-3577]. Multiple-sequence alignment indicates the active site Lys-Tyr pair consists of lysine 206 and tyrosine 150. Therefore, the K206M and Y150F mutants of HIcDH were prepared and characterized to test the potential roles of these residues as acid-base catalysts. The V/Et values of the K206M and Y150F mutant enzymes at pH 7.5 are decreased by approximately 2400- and approximately 680-fold, respectively, compared to that of wild-type HIcDH; the K(m) for HIc does not change significantly. V/Et and V/K(MgHIc)Et for the K206M mutant enzyme are pH-independent below pH 6 and decrease to a constant value above pH 7, while V/K(NAD)Et is independent over the pH range from 6.2 to 9.5. In the case of the Y150F mutant enzyme, V/Et and V/K(NAD)Et are pH-independent above pH 9.5 and decrease to a constant value below pH 8. This behavior can be compared to that of the wild-type enzyme, where V/Et decreases at high and low pH, giving pKa values of approximately 6.5 and approximately 9.5, respectively. Data were interpreted in terms of a group with a pKa of 6.5 that acts as a general base in the hydride transfer step and a group with a pKa of 9.5 that acts as a general acid to protonate C3 in the tautomerization reaction [Lin, Y., Volkman, J., Nicholas, K. M., Yamamoto, T., Eguchi, T., Nimmo, S. L., West, A. H., and Cook, P. F. (2008) Biochemistry 47, 4169-4180]. Solvent deuterium isotope effects on V and V/K(MgHIc) were near unity for the K206M mutant enzyme but approximately 2.2 for the Y150F mutant enzyme. The dramatic decreases in activity, the measured solvent deuterium isotope effects, and changes in the pH dependence of kinetic parameters compared to that of the wild type are consistent with K206 acting as a general base in the hydride transfer step of the wild-type enzyme but as a general acid in the Y150F mutant enzyme, replacing Y150 in the tautomerization reaction. In addition, Y150 acts as a general acid in the tautomerization reaction of the wild-type enzyme and replaces K206 as the general base in the hydride transfer step of the K206M mutant enzyme.

Correction of pre-steady-state KIEs for isotopic impurities and the consequences of kinetic isotope fractionation

We show, both experimentally and by kinetic modeling, that enzymatic single-turnover (pre-steady-state) H-transfer reactions can be significantly complicated by kinetic isotope fractionation. This fractionation results in the formation of more protiated than deuterated product and is a unique problem for pre-steady-state reactions. When observed rate constants are measured using rapid-mixing (e.g., stopped flow) methodologies, kinetic isotope fractionation can lead to a large underestimation of both the magnitude and temperature dependence of kinetic isotope effects (KIEs). This fractionation is related to the isotopic purity of the substrates used and highlights a major problem with experimental studies which measure KIEs with substrates that are not isotopically pure. As it is not always possible to prepare isotopically pure substrates, we describe two general methods for the correction, for known isotope impurities, of KIEs calculated from pre-steady-state measurements.

A stopped flow transient kinetic analysis of substrate binding and catalysis in Escherichia coli D-3-phosphoglycerate dehydrogenase

Pre-steady state, stopped flow analysis of Escherichia coli D-3-phosphoglycerate dehydrogenase was performed by following the fluorescence of protein tryptophan and the fluorescence resonance energy transfer from protein tryptophan to bound NADH. The results indicate that binding of substrates is ordered, with coenzyme, NADH, binding first. Furthermore, the analysis indicated that there are two sets of sites on the tetrameric enzyme that can be differentiated by their kinetic behavior. NADH binding was consistent with an initial binding event followed by a slow conformational change for each site. The slow conformational change is responsible for the apparent tight binding of NADH to the apoenzyme but is too slow to participate in the catalytic cycle when the enzyme is rapidly turning over. Subsequent binding of the substrate, alpha-ketoglutarate, was characterized by a rapid equilibrium binding event followed by a conformational change for each site. Catalysis in the direction of NAD(+) reduction showed a distinct burst of activity followed by a slow rate of turnover, indicating that the rate-limiting step is after hydride transfer. Catalysis in the direction of NADH oxidation did not display burst kinetics, indicating that the rate-limiting step is at or before the hydride transfer step. The burst data indicated that the rate of NAD(+) reduction (3.8 s(-1)) is similar to the k(cat) of the enzyme (2-3 s(-1)) in that direction. However, analysis of the reaction with deuterated NADH failed to show an effect on the velocity of the reaction with a V(H)/V(D)=1.07+/-0.06. None of the other rates determined by stopped flow analysis could account for the k(cat) of the enzyme in either direction (forward k(cat)=0.01 s(-1), reverse k(cat)=2-3 s(-1)), suggesting that the rate-limiting step in both directions is a conformational change in the enzyme that is not detected optically.

Chemical mechanism of homoisocitrate dehydrogenase from Saccharomyces cerevisiae

Homoisocitrate dehydrogenase (HIcDH, 3-carboxy-2-hydroxyadipate dehydrogenase) catalyzes the fourth reaction of the alpha-aminoadipate pathway for lysine biosynthesis, the conversion of homoisocitrate to alpha-ketoadipate using NAD as an oxidizing agent. A chemical mechanism for HIcDH is proposed on the basis of the pH dependence of kinetic parameters, dissociation constants for competitive inhibitors, and isotope effects. According to the pH-rate profiles, two enzyme groups act as acid-base catalysts in the reaction. A group with a p K a of approximately 6.5-7 acts as a general base accepting a proton as the beta-hydroxy acid is oxidized to the beta-keto acid, and this residue participates in all three of the chemical steps, acting to shuttle a proton between the C2 hydroxyl and itself. The second group acts as a general acid with a p K a of 9.5 and likely catalyzes the tautomerization step by donating a proton to the enol to give the final product. The general acid is observed in only the V pH-rate profile with homoisocitrate as a substrate, but not with isocitrate as a substrate, because the oxidative decarboxylation portion of the isocitrate reaction is limiting overall. With isocitrate as the substrate, the observed primary deuterium and (13)C isotope effects indicate that hydride transfer and decarboxylation steps contribute to rate limitation, and that the decarboxylation step is the more rate-limiting of the two. The multiple-substrate deuterium/ (13)C isotope effects suggest a stepwise mechanism with hydride transfer preceding decarboxylation. With homoisocitrate as the substrate, no primary deuterium isotope effect was observed, and a small (13)C kinetic isotope effect (1.0057) indicates that the decarboxylation step contributes only slightly to rate limitation. Thus, the chemical steps do not contribute significantly to rate limitation with the native substrate. On the basis of data from solvent deuterium kinetic isotope effects, viscosity effects, and multiple-solvent deuterium/ (13)C kinetic isotope effects, the proton transfer step(s) is slow and likely reflects a conformational change prior to catalysis.

Mechanism of proton transfer in the 3alpha-hydroxysteroid dehydrogenase/carbonyl reductase from Comamonas testosteroni

3alpha-hydroxysteroid dehydrogenase/carbonyl reductase from Comamonas testosteroni catalyzes the oxidation of androsterone with NAD(+) to form androstanedione and NADH with a concomitant releasing of protons to bulk solvent. To probe the proton transfer during the enzyme reaction, we used mutagenesis, chemical rescue, and kinetic isotope effects to investigate the release of protons. The kinetic isotope effects of (D)V and (D(2)O)V for wild-type enzyme are 1 and 2.1 at pL 10.4 (where L represents H, (2)H), respectively, and suggest a rate-limiting step in the intramolecular proton transfer. Substitution of alanine for Lys(159) changes the rate-limiting step to the hydride transfer, evidenced by an equal deuterium isotope effect of 1.8 on V(max) and V/K(androsterone) and no solvent kinetic isotope effect at saturating 3-(cyclohexylamino)propanesulfonic acid (CAPS). However, a value of 4.4 on V(max) is observed at 10 mm CAPS at pL 10.4, indicating a rate-limiting proton transfer. The rate of the proton transfer is blocked in the K159A and K159M mutants but can be rescued using exogenous proton acceptors, such as buffers, small primary amines, and azide. The Brønsted relationship between the log(V/K(d)(-base)Et) of the external amine (corrected for molecular size effects) and pK(a) is linear for the K159A mutant-catalyzed reaction at pH 10.4 (beta = 0.85 +/- 0.09) at 5 mm CAPS. These results show that proton transfer to the external base with a late transition state occurred in a rate-limiting step. Furthermore, a proton inventory on V/Et is bowl-shaped for both the wild-type and K159A mutant enzymes and indicates a two-proton transfer in the transition state from Tyr(155) to Lys(159) via 2'-OH of ribose.

Alpha-secondary isotope effects as probes of "tunneling-ready" configurations in enzymatic H-tunneling: insight from environmentally coupled tunneling models

Using alpha-secondary kinetic isotope effects (2 degrees KIEs) in conjunction with primary (1 degrees ) KIEs, we have investigated the mechanism of environmentally coupled hydrogen tunneling in the reductive half-reactions of two homologous flavoenzymes, morphinone reductase (MR) and pentaerythritol tetranitrate reductase (PETNR). We find exalted 2 degrees KIEs (1.17-1.18) for both enzymes, consistent with hydrogen tunneling. These 2 degrees KIEs, unlike 1 degrees KIEs, are independent of promoting motions-a nonequilibrium pre-organization of cofactor and active site residues that is required to bring the reactants into a "tunneling-ready" configuration. That these 2 degrees KIEs are identical suggests the geometries of the "tunneling-ready" configurations in both enzymes are indistinguishable, despite the fact that MR, but not PETNR, has a clearly temperature-dependent 1 degrees KIE. The work emphasizes the benefit of combining studies of 1 degrees and 2 degrees KIEs to report on pre-organization and local geometries within the context of contemporary environmentally coupled frameworks for H-tunneling.

Elucidation of a complete kinetic mechanism for a mammalian hydroxysteroid dehydrogenase (HSD) and identification of all enzyme forms on the reaction coordinate: the example of rat liver 3alpha-HSD (AKR1C9)

Hydroxysteroid dehydrogenases (HSDs) are essential for the biosynthesis and mechanism of action of all steroid hormones. We report the complete kinetic mechanism of a mammalian HSD using rat 3alpha-HSD of the aldo-keto reductase superfamily (AKR1C9) with the substrate pairs androstane-3,17-dione and NADPH (reduction) and androsterone and NADP(+) (oxidation). Steady-state, transient state kinetics, and kinetic isotope effects reconciled the ordered bi-bi mechanism, which contained 9 enzyme forms and permitted the estimation of 16 kinetic constants. In both reactions, loose association of the NADP(H) was followed by two conformational changes, which increased cofactor affinity by >86-fold. For androstane-3,17-dione reduction, the release of NADP(+) controlled k(cat), whereas the chemical event also contributed to this term. k(cat) was insensitive to [(2)H]NADPH, whereas (D)k(cat)/K(m) and the (D)k(lim) (ratio of the maximum rates of single turnover) were 1.06 and 2.06, respectively. Under multiple turnover conditions partial burst kinetics were observed. For androsterone oxidation, the rate of NADPH release dominated k(cat), whereas the rates of the chemical event and the release of androstane-3,17-dione were 50-fold greater. Under multiple turnover conditions full burst kinetics were observed. Although the internal equilibrium constant favored oxidation, the overall K(eq) favored reduction. The kinetic Haldane and free energy diagram confirmed that K(eq) was governed by ligand binding terms that favored the reduction reactants. Thus, HSDs in the aldo-keto reductase superfamily thermodynamically favor ketosteroid reduction.

Evidence from Raman spectroscopy that InhA, the mycobacterial enoyl reductase, modulates the conformation of the NADH cofactor to promote catalysis

InhA, the enoyl reductase from Mycobacterium tuberculosis, catalyzes the NADH-dependent reduction of trans-2-enoyl-ACPs. In the present work, Raman spectroscopy has been used to identify catalytically relevant changes in the conformation of the nicotinamide ring that occur when NADH binds to InhA. For 4(S)-NADD, there is an 11 cm-1 decrease in the wavenumber of the C4-D stretching band (nuC-D) and a 50% decrease in the width of this band upon binding to InhA. While a similar reduction in line width is observed for the corresponding band arising from 4(R)-NADD, nuC-D for this isomer increases 34 cm-1 upon binding to InhA. These changes in nuC-D indicate that the nicotinamide ring adopts a bound conformation in which the 4(S)C-D bond is in a pseudoaxial orientation. Mutagenesis of F149, a conserved active site residue close to the cofactor, demonstrates that this enzyme-induced modulation in cofactor structure is directly linked to catalysis. In contrast to the wild-type enzyme, Raman spectra of NADD bound to F149A InhA resemble those of NADD in solution. Consequently, F149A is no longer able to optimally position the cofactor for hydride transfer, which correlates with the 30-fold decrease in kcat and 2-fold increase in D(V/KNADH) caused by this mutation. These studies thus substantiate the proposal that hydride transfer is promoted by pseudoaxial positioning of the NADH pro-4S bond, and indicate that catalysis of substrate reduction by InhA results, in part, from correct orientation of the cofactor in the ground state.

A proposed proton shuttle mechanism for saccharopine dehydrogenase from Saccharomyces cerevisiae

Saccharopine dehydrogenase [N6-(glutaryl-2)-L-lysine:NAD oxidoreductase (L-lysine forming)] catalyzes the final step in the alpha-aminoadipate pathway for lysine biosynthesis. It catalyzes the reversible pyridine nucleotide-dependent oxidative deamination of saccharopine to generate alpha-Kg and lysine using NAD+ as an oxidizing agent. The proton shuttle chemical mechanism is proposed on the basis of the pH dependence of kinetic parameters, dissociation constants for competitive inhibitors, and isotope effects. In the direction of lysine formation, once NAD+ and saccharopine bind, a group with a pKa of 6.2 accepts a proton from the secondary amine of saccharopine as it is oxidized. This protonated general base then does not participate in the reaction again until lysine is formed at the completion of the reaction. A general base with a pKa of 7.2 accepts a proton from H2O as it attacks the Schiff base carbon of saccharopine to form the carbinolamine intermediate. The same residue then serves as a general acid and donates a proton to the carbinolamine nitrogen to give the protonated carbinolamine. Collapse of the carbinolamine is then facilitated by the same group accepting a proton from the carbinolamine hydroxyl to generate alpha-Kg and lysine. The amine nitrogen is then protonated by the group that originally accepted a proton from the secondary amine of saccharopine, and products are released. In the reverse reaction direction, finite primary deuterium kinetic isotope effects were observed for all parameters with the exception of V2/K(NADH), consistent with a steady-state random mechanism and indicative of a contribution from hydride transfer to rate limitation. The pH dependence, as determined from the primary isotope effect on DV2 and D(V2/K(Lys)), suggests that a step other than hydride transfer becomes rate-limiting as the pH is increased. This step is likely protonation/deprotonation of the carbinolamine nitrogen formed as an intermediate in imine hydrolysis. The observed solvent isotope effect indicates that proton transfer also contributes to rate limitation. A concerted proton and hydride transfer is suggested by multiple substrate/solvent isotope effects, as well as a proton transfer in another step, likely hydrolysis of the carbinolamine. In agreement, dome-shaped proton inventories are observed for V2 and V2/K(Lys), suggesting that proton transfer exists in at least two sequential transition states.

Promoting motions in enzyme catalysis probed by pressure studies of kinetic isotope effects

Use of the pressure dependence of kinetic isotope effects, coupled with a study of their temperature dependence, as a probe for promoting motions in enzymatic hydrogen-tunneling reactions is reported. Employing morphinone reductase as our model system and by using stopped-flow methods, we measured the hydride transfer rate (a tunneling reaction) as a function of hydrostatic pressure and temperature. Increasing the pressure from 1 bar (1 bar = 100 kPa) to 2 kbar accelerates the hydride transfer reaction when both protium (from 50 to 161 s(-1) at 25 degrees C) and deuterium (12 to 31 s(-1) at 25 degrees C) are transferred. We found that the observed primary kinetic isotope effect increases with pressure (from 4.0 to 5.2 at 25 degrees C), an observation incompatible with the Bell correction model for hydrogen tunneling but consistent with a full tunneling model. By numerical modeling, we show that both the pressure and temperature dependencies of the reaction rates are consistent with the framework of the environmentally coupled tunneling model of Kuznetsov and Ulstrup [Kuznetsov AM, Ulstrup J (1999) Can J Chem 77:1085-1096], providing additional support for the role of a promoting motion in the hydride tunneling reaction in morphinone reductase. Our study demonstrates the utility of "barrier engineering" by using hydrostatic pressure as a probe for tunneling regimes in enzyme systems and provides added and independent support for the requirement of promoting motions in such tunneling reactions.

Secoisolariciresinol dehydrogenase: mode of catalysis and stereospecificity of hydride transfer in Podophyllum peltatum

Secoisolariciresinol dehydrogenase (SDH) catalyzes the NAD+ dependent enantiospecific conversion of secoisolariciresinol into matairesinol. In Podophyllum species, (-)-matairesinol is metabolized into the antiviral compound, podophyllotoxin, which can be semi-synthetically converted into the anticancer agents, etoposide, teniposide and Etopophos. Matairesinol is also a precursor of the cancer-preventative "mammalian" lignan, enterolactone, formed in the gut following ingestion of, for example, various high fiber dietary foods, as well as being an intermediate to numerous defense compounds in vascular plants. This study investigated the mode of enantiospecific Podophyllum SDH catalysis, the order of binding, and the stereospecificity of hydride abstraction/transfer from secoisolariciresinol to NAD+. SDH contains a highly conserved catalytic triad (Ser153, Tyr167 and Lys171), whose activity was abolished with site-directed mutagenesis of Tyr167Ala and Lys171Ala, whereas mutagenesis of Ser153Ala only resulted in a much reduced catalytic activity. Isothermal titration calorimetry measurements indicated that NAD+ binds first followed by the substrate, (-)-secoisolariciresinol. Additionally, for hydride transfer, the incoming hydride abstracted from the substrate takes up the pro-S position in the NADH formed. Taken together, a catalytic mechanism for the overall enantiospecific conversion of (-)-secoisolariciresinol into (-)-matairesinol is proposed.

Synthesis and utility of 14C-labeled nicotinamide cofactors

A new method for the synthesis of the reduced form of beta-nicotinamide [U-14C]adenine dinucleotide 2(')-phosphate([Ad-14C]NADPH) is presented. The present synthesis results in a radioactive material with a specific activity that is greater than 220 mCi/mmol. This method could easily be adapted for syntheses of 14C-labeled NADH, NADP(+), or any nicotinamide cofactors with radiolabels in other positions. Since these cofactors are so ubiquitous, the use and applications of such labeled material has broad implications. The utility of the labeled cofactor for determination of substrates for nicotinamide-dependent enzymes in the nano- to femtomole scale, in alternative enzymatic assays, and in kinetic isotope effect studies is discussed.

H-tunneling in the multiple H-transfers of the catalytic cycle of morphinone reductase and in the reductive half-reaction of the homologous pentaerythritol tetranitrate reductase

The mechanism of flavin reduction in morphinone reductase (MR) and pentaerythritol tetranitrate (PETN) reductase, and flavin oxidation in MR, has been studied by stopped-flow and steady-state kinetic methods. The temperature dependence of the primary kinetic isotope effect for flavin reduction in MR and PETN reductase by nicotinamide coenzyme indicates that quantum mechanical tunneling plays a major role in hydride transfer. In PETN reductase, the kinetic isotope effect (KIE) is essentially independent of temperature in the experimentally accessible range, contrasting with strongly temperature-dependent reaction rates, consistent with a tunneling mechanism from the vibrational ground state of the reactive C-H/D bond. In MR, both the reaction rates and the KIE are dependent on temperature, and analysis using the Eyring equation suggests that hydride transfer has a major tunneling component, which, unlike PETN reductase, is gated by thermally induced vibrations in the protein. The oxidative half-reaction of MR is fully rate-limiting in steady-state turnover with the substrate 2-cyclohexenone and NADH at saturating concentrations. The KIE for hydride transfer from reduced flavin to the alpha/beta unsaturated bond of 2-cyclohexenone is independent of temperature, contrasting with strongly temperature-dependent reaction rates, again consistent with ground-state tunneling. A large solvent isotope effect (SIE) accompanies the oxidative half-reaction, which is also independent of temperature in the experimentally accessible range. Double isotope effects indicate that hydride transfer from the flavin N5 atom to 2-cyclohexenone, and the protonation of 2-cyclohexenone, are concerted and both the temperature-independent KIE and SIE suggest that this reaction also proceeds by ground-state quantum tunneling. Our results demonstrate the importance of quantum tunneling in the reduction of flavins by nicotinamide coenzymes. This is the first observation of (i) three H-nuclei in an enzymic reaction being transferred by tunneling and (ii) the utilization of both passive and active dynamics within the same native enzyme.

The crystal structure and stereospecificity of levodione reductase from Corynebacterium aquaticum M-13

The (6R)-2,2,6-trimethyl-1,4-cyclohexanedione (levodione) reductase (LVR) of the soil isolate bacterium Corynebacterium aquaticum M-13 is a NAD(H)-linked enzyme that catalyzes reversible oxidoreduction between (4R)-hydroxy-(6R)-2,2,6-trimethylcyclohexanone (actinol) and levodione. Here the crystal structure of a ternary complex of LVR with NADH and its inhibitor 2-methyl-2,4-pentanediol has been determined by molecular replacement and refined at 1.6-A resolution with a crystallographic R factor of 0.199. The overall structure is similar to those of other short-chain alcohol dehydrogenase/reductase enzymes. The positions of NADH and 2-methyl-2,4-pentanediol indicate the binding site of the substrate and identify residues that are likely to be important in the catalytic reaction. Modeling of the substrate binding in the active site suggests that the specificity of LVR is determined by electrostatic interactions between the negatively charged surface of Glu-103 of LVR and the positively charged surface on the re side of levodione. Mutant LVR enzymes in which Glu-103 is substituted with alanine (E103A), glutamine (E103Q), asparagines (E103N), or aspartic acid (E103D) show a 2-6-fold increase in Km values as compared with wild-type LVR and a much lower enantiomeric excess of the reaction products (60%) than the wild-type enzyme (95%). Together, these data indicate that Glu-103 has an important role in determining the stereospecificity of LVR.

Biosynthesis of 4-methylproline in cyanobacteria: cloning of nosE and nosF genes and biochemical characterization of the encoded dehydrogenase and reductase activities

The biosynthesis of the unusual amino acid 4-methylproline in the Nostoc genus of cyanobacteria was investigated on the genetic and enzymatic level. Two genes involved in the biosynthesis were cloned and the corresponding enzymes, a zinc-dependent long-chain dehydrogenase and a Delta(1)-pyrroline-5-carboxylic acid (P5C) reductase homologue, were overexpressed in Escherichia coli and biochemically characterized. Putative substrates were synthesized to test enzyme substrate specificities, and deuterium labeling studies were carried out to reveal the stereospecificities of the enzymatic reactions with respect to the substrates as well as to the coenzymes.

Stopped-flow kinetic studies of electron transfer in the reductase domain of neuronal nitric oxide synthase: re-evaluation of the kinetic mechanism reveals new enzyme intermediates and variation with cytochrome P450 reductase

The reduction by NADPH of the FAD and FMN redox centres in the isolated flavin reductase domain of calmodulin-bound rat neuronal nitric oxide synthase (nNOS) has been studied by anaerobic stopped-flow spectroscopy using absorption and fluorescence detection. We show by global analysis of time-dependent photodiode array spectra, single wavelength absorption and NADPH fluorescence studies, that at least four resolvable steps are observed in stopped-flow studies with NADPH and that flavin reduction is reversible. The first reductive step represents the rapid formation of an equilibrium between an NADPH-enzyme charge-transfer species and two-electron-reduced enzyme bound to NADP(+). The second and third steps represent further reduction of the enzyme flavins and NADP(+) release. The fourth step is attributed to the slow accumulation of an enzyme species that is inferred not to be relevant catalytically in steady-state reactions. Stopped-flow flavin fluorescence studies indicate the presence of slow kinetic phases, the timescales of which correspond to the slow phase observed in absorption and NADPH fluorescence transients. By analogy with stopped-flow studies of cytochrome P450 reductase, we attribute these slow fluorescence and absorption changes to enzyme disproportionation and/or conformational change. Unlike for the functionally related cytochrome P450 reductase, transfer of the first hydride equivalent from NADPH to nNOS reductase does not generate the flavin di-semiquinoid state. This indicates that internal electron transfer is relatively slow and is probably gated by NADP(+) release. Release of calmodulin from the nNOS reductase does not affect the kinetics of inter-flavin electron transfer under stopped-flow conditions, although the observed rate of formation of the equilibrium between the NADPH-oxidized enzyme charge-transfer species and two-electron-reduced enzyme bound to NADP(+) is modestly slower in calmodulin-depleted enzyme. Our studies indicate the need for significant re-interpretation of published kinetic data for electron transfer in the reductase domain of neuronal nitric oxide synthase.

Fast hydride transfer in proton-translocating transhydrogenase revealed in a rapid mixing continuous flow device

Transhydrogenase couples the redox reaction between NAD(H) and NADP(H) to proton translocation across a membrane. Coupling is achieved through changes in protein conformation. Upon mixing, the isolated nucleotide-binding components of transhydrogenase (dI, which binds NAD(H), and dIII, which binds NADP(H)) form a catalytic dI(2).dIII(1) complex, the structure of which was recently solved by x-ray crystallography. The fluorescence from an engineered Trp in dIII changes when bound NADP(+) is reduced. Using a continuous flow device, we have measured the Trp fluorescence change when dI(2).dIII(1) complexes catalyze reduction of NADP(+) by NADH on a sub-millisecond scale. At elevated NADH concentrations, the first-order rate constant of the reaction approaches 21,200 s(-1), which is larger than that measured for redox reactions of nicotinamide nucleotides in other, soluble enzymes. Rather high concentrations of NADH are required to saturate the reaction. The deuterium isotope effect is small. Comparison with the rate of the reverse reaction (oxidation of NADPH by NAD(+)) reveals that the equilibrium constant for the redox reaction on the complex is >36. This high value might be important in ensuring high turnover rates in the intact enzyme.

Characterization of hydride transfer to flavin adenine dinucleotide in neuronal nitric oxide synthase reductase domain: geometric relationship between the nicotinamide and isoalloxazine rings

Based on the similarity in both structure and function of the reductase domain of neuronal nitric oxide synthase (nNOSred) to that of NADPH-cytochrome P450 reductase (CPR), we determined whether the characteristics of hydride transfer from NADPH to flavin adenine dinucleotide (FAD) were similar for both proteins. Secondly, we questioned whether hydride transfer from NADPH to either nNOSred or holo-nNOS was rate limiting for reactions catalyzed by these two proteins. Utilizing 500 MHz proton NMR and deuterated substrate, we determined that the stereospecificity of hydride transfer from NADPH and the conformation of the nicotinamide ring around the glycosidic bond were similar between CPR and nNOSred. Specifically, nNOSred abstracts the A-side hydrogen from NADPH, and the nicotinamide ring is in the anti conformation. We determined that the rate of hydride transfer to FAD appears to become partially rate limiting only for exceptionally good electron acceptors such as cytochrome c. Hydride transfer is not rate limiting for NO. production under any conditions used in this study. Interestingly, the deuterium isotope effect was decreased in the cytochrome c reductase assay with both nNOS and nNOSred when the assays were conducted in high ionic strength buffer, suggesting an increase in the rate of hydride transfer to FAD. These results are in stark contrast to results obtained with CPR (D. S. Sem and C. B. Kasper, 1995, Biochemistry 34, 3391-3398) whereby hydride transfer is partially rate limiting at high, but not at low, ionic strength. The seemingly opposite results in deuterium isotope effect observed with CPR and nNOSred, under conditions of high and low ionic strength, suggest differences in structure and/or regulation of these important flavoproteins.

Structure-function aspects and inhibitor design of type 5 17beta-hydroxysteroid dehydrogenase (AKR1C3)

17beta-Hydroxysteroid dehydrogenase (17beta-HSD) type 5 has been cloned from human prostate and is identical to type 2 3alpha-HSD and is a member of the aldo-keto reductase (AKR) superfamily; it is formally AKR1C3. In vitro the homogeneous recombinant enzyme expressed in Escherichia coli functions as a 3-keto-, 17-keto- and 20-ketosteroid reductase and as a 3alpha-, 17beta- and 20alpha-hydroxysteroid oxidase. The enzyme will reduce 5alpha-DHT, Delta(4)-androstene-3,17-dione, estrone and progesterone to produce 3alpha-androstanediol, testosterone, 17beta-estradiol and 20alpha-hydroxprogesterone, respectively. It will also oxidize 3alpha-androstanediol, testosterone, 17beta-estradiol and 20alpha-hydroxyprogesterone to produce 5alpha-androstane-3,17-dione, Delta(4)-androstene-3,17-dione, and progesterone, respectively. Many of these properties are shared by the related AKR1C1, AKR1C2 and AKR1C4 isoforms. RT-PCR shows that AKR1C3 is dominantly expressed in the human prostate and mammary gland. Examination of k(cat)/K(m) for these reactions indicates that as a reductase it prefers 5alpha-dihydrotestosterone and 5alpha-androstane-3,17-dione as substrates to Delta(4)-androstene-3,17-dione, suggesting that in the prostate it favors the formation of inactive androgens. Its concerted reductase activity may, however, lead to a pro-estrogenic state in the breast since it will convert estrone to 17beta-estradiol; convert Delta(4)-androstene-3,17-dione to testosterone (which can be aromatized to 17beta-estradiol); and it will reduce progesterone to its inactive metabolite 20alpha-hydroxyprogesterone. Drawing on detailed structure-function analysis of the related rat 3alpha-HSD (AKR1C9), which shares 69% sequence identity with AKR1C3, it is predicted that AKR1C3 catalyzes an ordered bi bi mechanism, that the rate determining step is k(chem), and that an oxyanion prevails in the transition state. Based on these relationships steroidal-based inhibitors that compete with the steroid product would be desirable since they would act as uncompetitive inhibitors. With regards to transition state analogs steroid carboxylates and pyrazoles may be preferred while 3alpha, 17beta or 20alpha-spiro-oxiranes may act as mechanism-based inactivators.

Mechanism of human aldehyde reductase: characterization of the active site pocket

Human aldehyde reductase is a NADPH-dependent aldo-keto reductase that is closely related (65% identity) to aldose reductase, an enzyme involved in the pathogenesis of some diabetic and galactosemic complications. In aldose reductase, the active site residue Tyr48 is the proton donor in a hydrogen-bonding network involving residues Asp43/Lys77, while His110 directs the orientation of substrates in the active site pocket. Mutation of the homologous Tyr49 to phenylalamine or histidine (Y49F or Y49H) and of Lys79 to methionine (K79M) in aldehyde reductase yields inactive enzymes, indicating similar roles for these residues in the catalytic mechanism of aldehyde reductase. A H112Q mutant aldehyde reductase exhibited a substantial decrease in catalytic efficiency (kcat/Km) for hydrophilic (average 150-fold) and aromatic substrates (average 4200-fold) and 50-fold higher IC50 values for a variety of inhibitors than that of the wild-type enzyme. The data suggest that His112 plays a major role in determining the substrate specificity of aldehyde reductase, similar to that shown earlier for the homologous His110 in aldose reductase [Bohren, K. M., et. al. (1994) Biochemistry 33, 2021-2032]. Mutation of Ile298 or Val299 affected the kinetic parameters to a much lesser degree. Unlike native aldose reductase, which contains a thiol-sensitive Cys298, neither the I298C or V299C mutant exhibited any thiol sensitivity, suggesting a geometry of the active site pocket different from that in aldose reductase. Also different from aldose reductase, the detection of a significant primary deuterium isotope effect on kcat (1.48 +/- 0.02) shows that nucleotide exchange is only partially rate-limiting. Primary substrate and solvent deuterium isotope effects on the H112Q mutant suggest that hydride and proton transfers occur in two discrete steps with hydride transfer taking place first. Dissociation constants and spectroscopic and fluorimetric properties of nucleotide complexes with various mutants suggest that, in addition to Tyr49 and His112, Lys79 plays a hitherto unappreciated role in nucleotide binding. The mode of inhibition of aldehyde reductase by aldose reductase inhibitors (ARIs) is generally similar to that of aldose reductase and involves binding to the E:NADP+ complex, as shown by kinetic and direct inhibitor-binding experiments. The order of ARI potency was AL1576 (Ki = 60 nM) > tolrestat > ponalrestat > sorbinil > FK366 > zopolrestat > alrestatin (Ki = 148 microM). Our data on aldehyde reductase suggest that the active site pocket significantly differs from that of aldose reductase, possibly due to the participation of the C-terminal loop in its formation.

The involvement of NADP(H) binding and release in energy transduction by proton-translocating nicotinamide nucleotide transhydrogenase from Escherichia coli

Proton-translocating transhydrogenase was solubilised and purified from membranes of Escherichia coli. Consistent with recent evidence [Hutton, M., Day, J., Bizouarn, T. and Jackson, J.B. (1994) Eur. J. Biochem. 219, 1041-1051], at low pH and salt concentration, the enzyme catalysed rapid reduction of the NAD+ analogue AcPdAD+ by a combination of NADH and NADPH. At saturating concentrations of NADPH, the dependence of the steady-state rate on the concentrations of NADH and AcPdAD+ indicated that, with respect to these two nucleotides, the reaction proceeds by a ping-pong mechanism. High concentrations of either NADH or AcPdAD+ led to substrate inhibition. These observations support the view that, in this reaction, NADP(H) remains bound to the enzyme: AcPdAD+ is reduced by enzyme-bound NADPH, and NADH is oxidised by enzyme-bound NADP+, in a cyclic process. When this reaction was carried out with [4A-2H]NADH replacing [4A-1H]NADH, the rate was decreased by 46%, suggesting that the H- transfer steps are rate-limiting. In simple 'reverse' transhydrogenation, the reduction of AcPdAD+ was slower with [4B-2H]NADPH than with [4B-1H]NADPH when the reaction was performed at pH 8.0, but there was no deuterium isotope effect at pH 6.0. This indicates that H- transfer is rate-limiting at pH 8.0 and supports our earlier suggestion that NADP+ release from the enzyme is rate-limiting at low pH. The lack of a deuterium isotope effect in the reduction of thio-NADP+ by NADH at low pH is also consistent with the view that NADPH release from the enzyme is slow under these conditions. A steady-state rate equation is derived for the reduction of AcPdAD+ by NADPH plus NADH, assuming operation of the cyclic pathway. It adequately accounts for the pH dependence of the enzyme, for the features described above and for kinetic characteristics of E. coli transhydrogenase described in the literature.

Kinetic mechanism for the model reaction of NADPH-cytochrome P450 oxidoreductase with cytochrome c

The kinetic mechanism of NADPH-cytochrome P450 oxidoreductase (P450R) has been determined for the model reaction with cytochrome c3+. Although initial velocity studies show parallel patterns, consistent with a classical (one-site) ping-pong mechanism that precludes the formation of a ternary NADPH-P450R-cytochrome c3+ complex, product and dead-end inhibition results suggest a nonclassical (two-site) ping-pong mechanism [Northrop, D. B. (1969) J. Biol. Chem. 244, 5808-5819]. This mechanism is a hybrid of the random sequential (ternary complex) and ping-pong mechanisms, since ternary complexes can form as well as intermediate, modified forms of the enzyme that can be present in the absence of any bound substrate. The complete rate equation is derived for this mechanism, and values for Vmax, (V/K)NADPH, (V/K)cytc, and the corresponding Michaelis constants are presented in terms of microscopic rate constants along with the expected product inhibition patterns (Appendix). Inhibition by NADP+ is competitive versus NADPH and uncompetitive versus cytochrome c3+, while inhibition by cytochrome c2+ is competitive versus cytochrome c3+ and noncompetitive versus NADPH. These inhibition patterns are consistent with the proposed two-site mechanism. This mechanism would give the same initial velocity patterns as the classical one-site ping-pong mechanism, but it allows for the formation of a ternary complex, with NADPH and cytochrome c3+ reacting independently at two separate sites on P450R. The D(V/K)NADPH isotope effect is not affected by cytochrome c3+ concentration, consistent with our assumption (in deriving the rate equation) that binding at the two sites is independent. At the high ionic strength used in this study (850 mM), the mechanism is two-site ping-pong, with the electron acceptor site itself reacting with cytochrome c3+ in a tetra uni ping-pong manner.