Broadening the cofactor specificity of a thermostable alcohol dehydrogenase using rational protein design introduces novel kinetic transient behavior

Cofactor specificity in the aldo‐keto reductase (AKR) superfamily has been well studied, and several groups have reported the rational alteration of cofactor specificity in these enzymes. Although most efforts have focused on mesostable AKRs, several putative AKRs have recently been identified from hyperthermophiles. The few that have been characterized exhibit a strong preference for NAD(H) as a cofactor, in contrast to the NADP(H) preference of the mesophilic AKRs. Using the design rules elucidated from mesostable AKRs, we introduced two site‐directed mutations in the cofactor binding pocket to investigate cofactor specificity in a thermostable AKR, AdhD, which is an alcohol dehydrogenase from Pyrococcus furiosus. The resulting double mutant exhibited significantly improved activity and broadened cofactor specificity as compared to the wild‐type. Results of previous pre‐steady‐state kinetic experiments suggest that the high affinity of the mesostable AKRs for NADP(H) stems from a conformational change upon cofactor binding which is mediated by interactions between a canonical arginine and the 2′‐phosphate of the cofactor. Pre‐steady‐state kinetics with AdhD and the new mutants show a rich conformational behavior that is independent of the canonical arginine or the 2′‐phosphate. Additionally, experiments with the highly active double mutant using NADPH as a cofactor demonstrate an unprecedented transient behavior where the binding mechanism appears to be dependent on cofactor concentration. These results suggest that the structural features involved in cofactor specificity in the AKRs are conserved within the superfamily, but the dynamic interactions of the enzyme with cofactors are unexpectedly complex. Biotechnol. Bioeng. 2010;107: 763–774. © 2010 Wiley Periodicals, Inc.     

Dissection of the physiological interconversion of 5alpha-DHT and 3alpha-diol by rat 3alpha-HSD via transient kinetics shows that the chemical step is rate-determining: effect of mutating cofactor and substrate-binding pocket residues on catalysis

3Alpha-hydroxysteroid dehydrogenases (3alpha-HSDs) catalyze the interconversion between 5alpha-dihydrotestosterone (5alpha-DHT), the most potent androgen, and 3alpha-androstanediol (3alpha-diol), a weak androgen metabolite. To identify the rate-determining step in this physiologically important reaction, rat liver 3alpha-HSD (AKR1C9) was used as the protein model for the human homologues in fluorescence stopped-flow transient kinetic and kinetic isotope effect studies. Using single and multiple turnover experiments to monitor the NADPH-dependent reduction of 5alpha-DHT, it was found that k(lim) and k(max) values were identical to k(cat), indicating that chemistry is rate-limiting overall. Kinetic isotope effect measurements, which gave (D)k(cat) = 2.4 and (D)2(O)k(cat) = 3.0 at pL 6.0, suggest that the slow chemical transformation is significantly rate-limiting. When the NADP(+)-dependent oxidation of 3alpha-diol was monitored, single and multiple turnover experiments showed a k(lim) and burst kinetics consistent with product release as being rate-limiting overall. When NAD(+) was substituted for NADP(+), burst phase kinetics was eliminated, and k(max) was identical to k(cat). Thus with the physiologically relevant substrates 5alpha-DHT plus NADPH and 3alpha-diol plus NAD(+), the slowest event is chemistry. R276 forms a salt-linkage with the phosphate of 2'-AMP, and when it is mutated, tight binding of NAD(P)H is no longer observed [Ratnam, K., et al. (1999) Biochemistry 38, 7856-7864]. The R276M mutant also eliminated the burst phase kinetics observed for the NADP(+)-dependent oxidation of 3alpha-diol. The data with the R276M mutant confirms that the release of the NADPH product is the slow event; and in its absence, chemistry becomes rate-limiting. W227 is a critical hydrophobic residue at the steroid binding site, and when it is mutated to alanine, k(cat)/K(m) for oxidation is significantly depressed. Burst phase kinetics for the NADP(+)-dependent turnover of 3alpha-diol by W227A was also abolished. In the W227A mutant, the slow release of NADPH is no longer observed since the chemical transformation is now even slower. Thus, residues in the cofactor and steroid-binding site can alter the rate-determining step in the NADP(+)-dependent oxidation of 3alpha-diol to make chemistry rate-limiting overall.     

Mutation of nicotinamide pocket residues in rat liver 3 alpha-hydroxysteroid dehydrogenase reveals different modes of cofactor binding

Rat liver 3alpha-hydroxysteroid dehydrogenase (3alpha-HSD), an aldo-keto reductase, binds NADP(+) in an extended anti-conformation across an (alpha/beta)(8)-barrel. The orientation of the nicotinamide ring, which permits stereospecific transfer of the 4-pro-R hydride from NAD(P)H to substrate, is achieved by hydrogen bonds formed between the C3-carboxamide of the nicotinamide ring and Ser 166, Asn 167, and Gln 190 and by pi-stacking between this ring and Tyr 216. These residues were mutated to yield S166A, N167A, Q190A, and Y216S. In these mutants, K(d)(NADP(H)) increased by 2-11-fold but without a significant change in K(d)(NAD(H)). Steady-state kinetic parameters showed that K(m)(NADP)()+ increased 13-151-fold, and this was accompanied by comparable decreases in k(cat)/K(m)(NADP)()+. By contrast, K(m)(NAD)()+ increased 4-8-fold, but changes in k(cat)/K(m)(NAD)()+ were more dramatic and ranged from 23- to 930-fold. Corresponding changes in binding energies indicated that each residue contributed equally to the binding of NADP(H) in the ground and transition states. However, the same residues stabilized the binding of NAD(H) only in the transition state. These observations suggest that different modes of binding exist for NADP(H) and NAD(H). Importantly, these modes were revealed by mutating residues in the nicotinamide pocket indicating that direct interactions with the 2'-phosphate in the adenine mononucleotide is not the sole determinant of cofactor preference. The single mutations were unable to invert or racemize the stereochemistry of hydride transfer even though the nicotinamide pocket can accommodate both anti- and syn-conformers once the necessary hydrogen bonds are eliminated. When 4-pro-R-[(3)H]NADH was used to monitor incorporation into [(14)C]-5alpha-dihydrotestosterone, a decrease in the (3)H:(14)C ratio was observed in the mutants relative to wild-type enzyme reflecting a pronounced primary kinetic isotope effect. This observation coupled with the change in the binding energy for NAD(P)(H) in the transition state suggests that these mutants have altered the reaction trajectory for hydride transfer.     

Transient-phase kinetic studies on the nucleotide binding to 3alpha-hydroxysteroid dehydrogenase from Pseudomonas sp. B-0831 using fluorescence stopped-flow procedures.

The dual nucleotide cofactor-specific enzyme, 3alpha-hydroxysteroid dehydrogenase (3alpha-HSD) from Pseudomonas sp. B-0831, is a member of the short-chain dehydrogenase/reductase (SDR) superfamily. Transient-phase kinetic studies using the fluorescence stopped-flow method were conducted with 3alpha-HSD to characterize the nucleotide binding mechanism. The binding of oxidized nucleotides, NAD(+), NADP(+) and nicotinic acid adenine dinucleotide (NAAD(+)), agreed well with a one-step mechanism, while that of reduced nucleotide, NADH, showed a two-step mechanism. This difference draws attention to previous characteristic findings on rat liver 3alpha-HSD, which is a member of the aldo-keto reductase (AKR) superfamily. Although functionally similar, AKRs are structurally different from SDRs. The dissociation rate constants associated with the enzyme-nucleotide complex formation were larger than the k(cat) values for either oxidation or reduction of substrates, indicating that the release of cofactors is not rate-limiting overall. It should also be noted that k(cat) for a substrate, cholic acid, with NADP(+) was only 6% of that with NAD(+), and no catalytic activity was detectable with NAAD(+), despite the similar binding affinities of nucleotides. These results suggest that a certain type of nucleotide can modulate nucleotide-binding mode and further the catalytic function of the enzyme.     

The 2'-phosphate of NADP is responsible for proper orientation of the nicotinamide ring in the oxidative decarboxylation reaction catalyzed by sheep liver 6-phosphogluconate dehydrogenase

Sheep liver 6-phosphogluconate dehydrogenase shows a high specificity for NADP, with a much lower affinity for NAD. Discrimination between NADP and NAD suggests that the interactions between the 2'-phosphate and 6-phosphogluconate dehydrogenase contribute most of the binding energy for NADP. There are three active site residues, Asn-32, Arg-33, and Thr-34, that hydrogen-bond to the 2'-phosphate of NADP according to the crystal structure of the E.Nbr(8)ADP complex. In this study alanine mutagenesis was used to probe the contribution of each of the three residues to binding the cofactor and to catalysis. All mutant enzymes exhibit no significant change in V/E(t) or K(6PG) but an increase in K(NADP) that ranges from 6- to 80-fold. All mutant enzymes also exhibit at least a 7-fold increase in the primary kinetic (13)C-isotope effect-1, indicating that the decarboxylation step has become more rate-limiting. Data are consistent with significant roles for Asn-32, Arg-33, and Thr-34 in providing binding energy for NADP, and more importantly, the 2'-phosphate of NADP is required for proper orientation of the cofactor to allow rotation about the N-glycosidic bond as it is reduced in the hydride transfer step.     

Molecular determinants of steroid recognition and catalysis in aldo-keto reductases. Lessons from 3alpha-hydroxysteroid dehydrogenase.

Hydroxysteroid Dehydrogenases (HSDs) regulate the occupancy of steroid hormone receptors by converting active steroid hormones into their cognate inactive metabolites. HSDs belong to either the Short-chain Dehydrogenase/Reductases (SDRs) or the Aldo-Keto Reductases (AKRs). The AKRs include virtually all mammalian 3alpha-HSDs, Type 5 17beta-HSD, ovarian 20alpha-HSDs as well as the steroid 5beta-reductases. Selective inhibitors of 3alpha-HSD isoforms could control occupancy of the androgen and GABA(A) receptors, while broader based AKR inhibitors targeting 3alpha-HSD, 20alpha-HSD and prostaglandin F2alpha synthase could maintain pregnancy. We have determined three X-ray crystal structures of rat liver 3alpha-HSD, a representative AKR. These structures are of the apoenzyme (E), the binary-complex (E.NADP-), and the ternary complex (E.NADP+.testosterone). These structures are being used with site-directed mutagenesis to define the molecular determinants of steroid recognition and catalysis as a first step in rational inhibitor design. A conserved catalytic tetrad (Tyr55, Lys84, His117 and Asp50) participates in a 'proton-relay' in which Tyr55 acts as general acid/base catalyst. Its bifunctionality relies on contributions from His117 and Lys84 which alter the pKb and pKa, respectively of this residue. Point mutation of the tetrad results in different enzymatic activities. H117E mutants display 5beta-reductase activity while Y55F and Y55S mutants retain quinone reductase activity. Our results suggest that different transition states are involved in these reaction mechanisms. The ternary complex structure shows that the mature steroid binding pocket is comprised of ten residues recruited from five loops, and that there is significant movement of a C-terminal loop on binding ligand. Mutagenesis of pocket tryptophans shows that steroid substrates and classes of nonsteroidal inhibitors exhibit different binding modes which may reflect ligand-induced loop movement. Exploitation of these findings using steroidal and nonsteroidal mechanism based inactivators may lead to selective and broad based AKR inhibitors.     

Role of aspartic acid 38 in the cofactor specificity of Drosophila alcohol dehydrogenase.

Drosophila alcohol dehydrogenase (ADH), an NAD(+)-dependent dehydrogenase, shares little sequence similarity with horse liver ADH. However, these two enzymes do have substantial similarity in their secondary structure at the NAD(+)-binding domain [Benyajati, C., Place, A. P., Powers, D. A. & Sofer, W. (1981) Proc. Natl Acad. Sci. USA 78, 2717-2721]. Asp38, a conserved residue between Drosophila and horse liver ADH, appears to interact with the hydroxyl groups of the ribose moiety in the AMP portion of NAD+. A secondary-structure comparison between the nucleotide-binding domain of NAD(+)-dependent enzymes and that of NADP(+)-dependent enzymes also suggests that Asp38 could play an important role in cofactor specificity. Mutating Asp38 of Drosophila ADH into Asn38 decreases Km(app)NADP 62-fold and increases kcat/Km(app)NADP 590-fold at pH 9.8, when compared with wild-type ADH. These results suggest that Asp38 is in the NAD(+)-binding domain and its substituent, Asn38, allows Drosophila ADH to use both NAD+ and NADP+ as its cofactor. The observations from the experiments of thermal denaturation and kinetic measurement with pH also confirm that the repulsion between the negative charges of Asp38 and 2'-phosphate of NADP+ is the major energy barrier for NADP+ to serve as a cofactor for Drosophila ADH.     

Role of residues in the adenosine binding site of NAD of the Ascaris suum malic enzyme

Ascaris suum mitochondrial malic enzyme catalyzes the divalent metal ion dependent conversion of l-malate to pyruvate and CO(2), with concomitant reduction of NAD(P) to NAD(P)H. In this study, some of the residues that form the adenosine binding site of NAD were mutated to determine their role in binding of the cofactor and/or catalysis. D361, which is completely conserved among species, is located in the dinucleotide-binding Rossmann fold and makes a salt bridge with R370, which is also highly conserved. D361 was mutated to E, A and N. R370 was mutated to K and A. D361E and A mutant enzymes were inactive, likely a result of the increase in the volume in the case of the D361E mutant enzyme that caused clashes with the surrounding residues, and loss of the ionic interaction between D361 and R370, for D361A. Although the K(m) for the substrates and isotope effect values did not show significant changes for the D361N mutant enzyme, V/E(t) decreased by 1400-fold. Data suggested the nonproductive binding of the cofactor, giving a low fraction of active enzyme. The R370K mutant enzyme did not show any significant changes in the kinetic parameters, while the R370A mutant enzyme gave a slight change in V/E(t), contrary to expectations. Overall, results suggest that the salt bridge between D361 and R370 is important for maintaining the productive conformation of the NAD binding site. Mutation of residues involved leads to nonproductive binding of NAD. The interaction stabilizes one of the Rossmann fold loops that NAD binds. Mutation of H377 to lysine, which is conserved in NADP-specific malic enzymes and proposed to be a cofactor specificity determinant, did not cause a shift in cofactor specificity of the Ascaris malic enzyme from NAD to NADP. However, it is confirmed that this residue is an important second layer residue that affects the packing of the first layer residues that directly interact with the cofactor.     

Crystal structures of mouse 17alpha-hydroxysteroid dehydrogenase (apoenzyme and enzyme-NADP(H) binary complex): identification of molecular determinants responsible for the unique 17alpha-reductive activity of this enzyme

Very recently, the mouse 17alpha-hydroxysteroid dehydrogenase (m17alpha-HSD), a member of the aldo-keto reductase (AKR) superfamily, has been characterized and identified as the unique enzyme able to catalyze efficiently and in a stereospecific manner the conversion of androstenedione (Delta4) into epitestosterone (epi-T), the 17alpha-epimer of testosterone. Indeed, the other AKR enzymes that significantly reduce keto groups situated at position C17 of the steroid nucleus, the human type 3 3alpha-HSD (h3alpha-HSD3), the human and mouse type 5 17beta-HSD, and the rabbit 20alpha-HSD, produce only 17beta-hydroxy derivatives, although they possess more than 70% amino acid identity with m17alpha-HSD. Structural comparisons of these highly homologous enzymes thus offer an excellent opportunity of identifying the molecular determinants responsible for their 17alpha/17beta-stereospecificity. Here, we report the crystal structure of the m17alpha-HSD enzyme in its apo-form (1.9 A resolution) as well as those of two different forms of this enzyme in binary complex with NADP(H) (2.9 A and 1.35 A resolution). Interestingly, one of these binary complex structures could represent a conformational intermediate between the apoenzyme and the active binary complex. These structures provide a complete picture of the NADP(H)-enzyme interactions involving the flexible loop B, which can adopt two different conformations upon cofactor binding. Structural comparison with binary complexes of other AKR1C enzymes has also revealed particularities of the interaction between m17alpha-HSD and NADP(H), which explain why it has been possible to crystallize this enzyme in its apo form. Close inspection of the m17alpha-HSD steroid-binding cavity formed upon cofactor binding leads us to hypothesize that the residue at position 24 is of paramount importance for the stereospecificity of the reduction reaction. Mutagenic studies have showed that the m17alpha-HSD(A24Y) mutant exhibited a completely reversed stereospecificity, producing testosterone only from Delta4, whereas the h3alpha-HSD3(Y24A) mutant acquires the capacity to metabolize Delta4 into epi-T.     

Molecular determinants of the cofactor specificity of ribitol dehydrogenase, a short-chain dehydrogenase/reductase

Ribitol dehydrogenase from Zymomonas mobilis (ZmRDH) catalyzes the conversion of ribitol to d-ribulose and concomitantly reduces NAD(P)(+) to NAD(P)H. A systematic approach involving an initial sequence alignment-based residue screening, followed by a homology model-based screening and site-directed mutagenesis of the screened residues, was used to study the molecular determinants of the cofactor specificity of ZmRDH. A homologous conserved amino acid, Ser156, in the substrate-binding pocket of the wild-type ZmRDH was identified as an important residue affecting the cofactor specificity of ZmRDH. Further insights into the function of the Ser156 residue were obtained by substituting it with other hydrophobic nonpolar or polar amino acids. Substituting Ser156 with the negatively charged amino acids (Asp and Glu) altered the cofactor specificity of ZmRDH toward NAD(+) (S156D, [k(cat)/K(m)(,NAD)]/[k(cat)/K(m)(,NADP)] = 10.9, where K(m)(,NAD) is the K(m) for NAD(+) and K(m)(,NADP) is the K(m) for NADP(+)). In contrast, the mutants containing positively charged amino acids (His, Lys, or Arg) at position 156 showed a higher efficiency with NADP(+) as the cofactor (S156H, [k(cat)/K(m)(,NAD)]/[k(cat)/K(m)(,NADP)] = 0.11). These data, in addition to those of molecular dynamics and isothermal titration calorimetry studies, suggest that the cofactor specificity of ZmRDH can be modulated by manipulating the amino acid residue at position 156.     

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.     

Human 20alpha-hydroxysteroid dehydrogenase: crystallographic and site-directed mutagenesis studies lead to the identification of an alternative binding site for C21-steroids

Human 20alpha-hydroxysteroid dehydrogenase (h20alpha-HSD; AKR1C1) catalyzes the transformation of progesterone (Prog) into 20alpha-hydroxy-progesterone (20alpha-OHProg). Although h20alpha-HSD shares 98% sequence identity with human type 3 3alpha-HSD (h3alpha-HSD3, AKR1C2), these two enzymes differ greatly in their activities. In order to explain these differences, we have solved the crystal structure of h20alpha-HSD in a ternary complex with NADP(+) and 20alpha-OHProg at 1.59A resolution. The steroid is stabilized by numerous hydrophobic interactions and a hydrogen bond between its O20 and the N(epsilon ) atom of His222. This new interaction prevents the formation of a hydrogen bond with the cofactor, as seen in h3alpha-HSD3 ternary complexes. By combining structural, direct mutagenesis and kinetic studies, we found that the H(222)I substitution decreases the K(m) value for the cofactor 95-fold. With these results, we hypothesize that the rotation of the lateral chain of His222 could be a mediating step between the transformation of Prog and the release of the cofactor. Moreover, crystal structure analysis and direct mutagenesis experiments lead us to identify a new residue involved in the binding of Prog. Indeed, the R(304)L substitution leads to a 65-fold decrease in the K(m) value for Prog reduction. We thus propose that Prog is maintained in a new steroid-binding site composed mainly of residues found in the carboxy-terminal region of the protein.     

Mechanism of coenzyme binding to human methionine synthase reductase revealed through the crystal structure of the FNR-like module and isothermal titration calorimetry

Human methionine synthase reductase (MSR) is a 78 kDa flavoprotein that regenerates the active form of cobalamin-dependent methionine synthase (MS). MSR contains one FAD and one FMN cofactor per polypeptide and functions in the sequential transfer of reducing equivalents from NADPH to MS via its flavin centers. We report the 1.9 A crystal structure of the NADP+-bound FNR-like module of MSR that spans the NADP(H)-binding domain, the FAD-binding domain, the connecting domain, and part of the extended hinge region, a feature unique to MSR. The overall fold of the protein is similar to that of the corresponding domains of the related diflavin reductase enzymes cytochrome P450 reductase and neuronal nitric oxide synthase (NOS). However, the extended hinge region of MSR, which is positioned between the NADP(H)/FAD- and FMN-binding domains, is in an unexpected orientation with potential implications for the mechanism of electron transfer. Compared with related flavoproteins, there is structural variation in the NADP(H)-binding site, in particular regarding those residues that interact with the 2'-phosphate and the pyrophosphate moiety of the coenzyme. The lack of a conserved binding determinant for the 2'-phosphate does not weaken the coenzyme specificity for NADP(H) over NAD(H), which is within the range expected for the diflavin oxidoreductase family of enzymes. Isothermal titration calorimetry reveals a binding constant of 37 and 2 microM for binding of NADP+ and 2',5'-ADP, respectively, for the ligand-protein complex formed with full-length MSR or the isolated FNR module. These values are consistent with Ki values (36 microM for NADP+ and 1.4 microM for 2',5'-ADP) obtained from steady-state inhibition studies. The relatively weaker binding of NADP+ to MSR compared with other members of the diflavin oxidoreductase family might arise from unique electrostatic repulsive forces near the 5'-pyrophosphate moiety and/or increased hydrophobic stacking between Trp697 and the re face of the FAD isoalloxazine ring. Small structural permutations within the NADP(H)-binding cleft have profound affects on coenzyme binding, which likely retards catalytic turnover of the enzyme in the cell. The biological implications of an attenuated mechanism of MS reactivation by MSR on methionine and folate metabolism are discussed.     

Determinants of the dual cofactor specificity and substrate cooperativity of the human mitochondrial NAD(P)+-dependent malic enzyme: functional roles of glutamine 362

The human mitochondrial NAD(P)+-dependent malic enzyme (m-NAD-ME) is a malic enzyme isoform with dual cofactor specificity and substrate binding cooperativity. Previous kinetic studies have suggested that Lys362 in the pigeon cytosolic NADP+-dependent malic enzyme has remarkable effects on the binding of NADP+ to the enzyme and on the catalytic power of the enzyme (Kuo, C. C., Tsai, L. C., Chin, T. Y., Chang, G.-G., and Chou, W. Y. (2000) Biochem. Biophys. Res. Commun. 270, 821-825). In this study, we investigate the important role of Gln362 in the transformation of cofactor specificity from NAD+ to NADP+ in human m-NAD-ME. Our kinetic data clearly indicate that the Q362K mutant shifted its cofactor preference from NAD+ to NADP+. The Km(NADP) and kcat(NADP) values for this mutant were reduced by 4-6-fold and increased by 5-10-fold, respectively, compared with those for the wild-type enzyme. Furthermore, up to a 2-fold reduction in Km(NADP)/Km(NAD) and elevation of kcat(NADP)/kcat(NAD) were observed for the Q362K enzyme. Mutation of Gln362 to Ala or Asn did not shift its cofactor preference. The Km(NADP)/Km(NAD) and kcat(NADP)/kcat(NAD) values for Q362A and Q362N were comparable with those for the wild-type enzyme. The DeltaG values for Q362A and Q362N with either NAD+ or NADP+ were positive, indicating that substitution of Gln with Ala or Asn at position 362 brings about unfavorable cofactor binding at the active site and thus significantly reduces the catalytic efficiency. Our data also indicate that the cooperative binding of malate became insignificant in human m-NAD-ME upon mutation of Gln362 to Lys because the sigmoidal phenomenon appearing in the wild-type enzyme was much less obvious that that in Q362K. Therefore, mutation of Gln362 to Lys in human m-NAD-ME alters its kinetic properties of cofactor preference, malate binding cooperativity, and allosteric regulation by fumarate. However, the other Gln362 mutants, Q362A and Q362N, have conserved malate binding cooperativity and NAD+ specificity. In this study, we provide clear evidence that the single mutation of Gln362 to Lys in human m-NAD-ME changes it to an NADP+-dependent enzyme, which is characteristic because it is non-allosteric, non-cooperative, and NADP+-specific.     

Mechanistic investigation of a highly active phosphite dehydrogenase mutant and its application for NADPH regeneration

NAD(P)H regeneration is important for biocatalytic reactions that require these costly cofactors. A mutant phosphite dehydrogenase (PTDH-E175A/A176R) that utilizes both NAD and NADP efficiently is a very promising system for NAD(P)H regeneration. In this work, both the kinetic mechanism and practical application of PTDH-E175A/A176R were investigated for better understanding of the enzyme and to provide a basis for future optimization. Kinetic isotope effect studies with PTDH-E175A/A176R showed that the hydride transfer step is (partially) rate determining with both NAD and NADP giving (D)V values of 2.2 and 1.7, respectively, and (D)V/K(m,phosphite) values of 1.9 and 1.7, respectively. To better comprehend the relaxed cofactor specificity, the cofactor dissociation constants were determined utilizing tryptophan intrinsic fluorescence quenching. The dissociation constants of NAD and NADP with PTDH-E175A/A176R were 53 and 1.9 microm, respectively, while those of the products NADH and NADPH were 17.4 and 1.22 microm, respectively. Using sulfite as a substrate mimic, the binding order was established, with the cofactor binding first and sulfite binding second. The low dissociation constant for the cofactor product NADPH combined with the reduced values for (D)V and k(cat) implies that product release may become partially rate determining. However, product inhibition does not prevent efficient in situ NADPH regeneration by PTDH-E175A/A176R in a model system in which xylose was converted into xylitol by NADP-dependent xylose reductase. The in situ regeneration proceeded at a rate approximately fourfold faster with PTDH-E175A/A176R than with either WT PTDH or a NADP-specific Pseudomonas sp.101 formate dehydrogenase mutant with a total turnover number for NADPH of 2500.     

Structural studies of the pigeon cytosolic NADP(+)-dependent malic enzyme.

Malic enzymes are widely distributed in nature, and have important biological functions. They catalyze the oxidative decarboxylation of malate to produce pyruvate and CO(2) in the presence of divalent cations (Mg(2+), Mn(2+)). Most malic enzymes have a clear selectivity for the dinucleotide cofactor, being able to use either NAD(+) or NADP(+), but not both. Structural studies of the human mitochondrial NAD(+)-dependent malic enzyme established that malic enzymes belong to a new class of oxidative decarboxylases. Here we report the crystal structure of the pigeon cytosolic NADP(+)-dependent malic enzyme, in a closed form, in a quaternary complex with NADP(+), Mn(2+), and oxalate. This represents the first structural information on an NADP(+)-dependent malic enzyme. Despite the sequence conservation, there are large differences in several regions of the pigeon enzyme structure compared to the human enzyme. One region of such differences is at the binding site for the 2'-phosphate group of the NADP(+) cofactor, which helps define the cofactor selectivity of the enzymes. Specifically, the structural information suggests Lys362 may have an important role in the NADP(+) selectivity of the pigeon enzyme, confirming our earlier kinetic observations on the K362A mutant. Our structural studies also revealed differences in the organization of the tetramer between the pigeon and the human enzymes, although the pigeon enzyme still obeys 222 symmetry.     

Crystal structure of the MJ0490 gene product of the hyperthermophilic archaebacterium Methanococcus jannaschii, a novel member of the lactate/malate family of dehydrogenases

The MJ0490 gene, one of the only two genes of Methanococcus jannaschii showing sequence similarity to the lactate/malate family of dehydrogenases, was classified initially as coding for a putative l-lactate dehydrogenase (LDH). It has been re-classified as a malate dehydrogenase (MDH) gene, because it shows significant sequence similarity to MT0188, MDH II from Methanobacterium thermoautotrophicum strain DeltaH. The three-dimensional structure of its gene product has been determined in two crystal forms: a "dimeric" structure in the orthorhombic crystal at 1.9 A resolution and a "tetrameric" structure in the tetragonal crystal at 2.8 A. These structures share a similar subunit fold with other LDHs and MDHs. The tetrameric structure resembles typical tetrameric LDHs. The dimeric structure is equivalent to the P-dimer of tetrameric LDHs, unlike dimeric MDHs, which correspond to the Q-dimer. The structure reveals that the cofactor NADP(H) is bound at the active site, despite the fact that it was not intentionally added during protein purification and crystallization. The preference of NADP(H) over NAD(H) has been supported by activity assays. The cofactor preference is explained by the presence of a glycine residue in the cofactor binding pocket (Gly33), which replaces a conserved aspartate (or glutamate) residue in other NAD-dependent LDHs or MDHs. Preference for NADP(H) is contributed by hydrogen bonds between the oxygen atoms of the monophosphate group and the ribose sugar of adenosine in NADP(H) and the side-chains of Ser9, Arg34, His36, and Ser37. The MDH activity of MJ0490 is made possible by Arg86, which is conserved in MDHs but not in LDHs. The enzymatic assay showed that the MJ0490 protein possesses the fructose-1,6-bisphosphate-activated LDH activity (reduction). Thus the MJ0490 gene product appears to be a novel member of the lactate/malate dehydrogenase family, displaying an LDH scaffold and exhibiting a relaxed substrate and cofactor specificities in NADP(H) and NAD(H)-dependent malate and lactate dehydrogenase reactions.     

Kinetic studies of dogfish liver glutamate dehydrogenase.

Initial-rate studies were made of the oxidation of L-glutamate by NAD+ and NADP+ catalysed by highly purified preparations of dogfish liver glutamate dehydrogenase. With NAD+ as coenzyme the kinetics show the same features of coenzyme activation as seen with the bovine liver enzyme [Engel & Dalziel (1969) Biochem. J. 115, 621--631]. With NADP+ as coenzyme, initial rates are much slower than with NAD+, and Lineweaver--Burk plots are linear over extended ranges of substrate and coenzyme concentration. Stopped-flow studies with NADP+ as coenzyme give no evidence for the accumulation of significant concentrations of NADPH-containing complexes with the enzyme in the steady state. Protection studies against inactivation by pyridoxal 5'-phosphate indicate that NAD+ and NADP+ give the same degree of protection in the presence of sodium glutarate. The results are used to deduce information about the mechanism of glutamate oxidation by the enzyme. Initial-rate studies of the reductive amination of 2-oxoglutarate by NADH and NADPH catalysed by dogfish liver glutamate dehydrogenase showed that the kinetic features of the reaction are very similar with both coenzymes, but reactions with NADH are much faster. The data show that a number of possible mechanisms for the reaction may be discarded, including the compulsory mechanism (previously proposed for the enzyme) in which the sequence of binding is NAD(P)H, NH4+ and 2-oxoglutarate. The kinetic data suggest either a rapid-equilibrium random mechanism or the compulsory mechanism with the binding sequence NH4+, NAD(P)H, 2-oxoglutarate. However, binding studies and protection studies indicate that coenzyme and 2-oxoglutarate do bind to the free enzyme.     

The cofactor preference of glucose-6-phosphate dehydrogenase from Escherichia coli--modeling the physiological production of reduced cofactors

In Escherichia coli, the pentose phosphate pathway is one of the main sources of NADPH. The first enzyme of the pathway, glucose-6-phosphate dehydrogenase (G6PDH), is generally considered an exclusive NADPH producer, but a rigorous assessment of cofactor preference has yet to be reported. In this work, the specificity constants for NADP and NAD for G6PDH were determined using a pure enzyme preparation. Absence of the phosphate group on the cofactor leads to a 410-fold reduction in the performance of the enzyme. Furthermore, the contribution of the phosphate group to binding of the transition state to the active site was calculated to be 3.6 kcal·mol(-1). In order to estimate the main kinetic parameters for NAD(P) and NAD(P)H, we used the classical initial-rates approach, together with an analysis of reaction time courses. To achieve this, we developed a new analytical solution to the integrated Michaelis-Menten equation by including the effect of competitive product inhibition using the ω-function. With reference to relevant kinetic parameters and intracellular metabolite concentrations reported by others, we modeled the sensitivity of reduced cofactor production by G6PDH as a function of the redox ratios of NAD/NADH (rR(NAD)) and NADP/NADPH (rR(NADP)). Our analysis shows that NADPH production sharply increases within the range of thermodynamically feasible values of rR(NADP), but NADH production remains low within the range feasible for rR(NAD). Nevertheless, we show that certain combinations of rR(NADP) and rR(NAD) sustain greater levels of NADH production over NADPH.     

Mouse 17alpha-hydroxysteroid dehydrogenase (AKR1C21) binds steroids differently from other aldo-keto reductases: identification and characterization of amino acid residues critical for substrate binding

The mouse 17alpha-hydroxysteroid dehydrogenase (m17alpha-HSD) is the unique known member of the aldo-keto reductase (AKR) superfamily able to catalyze efficiently and in a stereospecific manner the conversion of androstenedione (Delta4) into epi-testosterone (epi-T), the 17alpha-epimer of testosterone. Structural and mutagenic studies had already identified one of the residues delineating the steroid-binding cavity, A24, as the major molecular determinant for the stereospecificity of m17alpha-HSD. We report here a ternary complex crystal structure (m17alpha-HSD:NADP(+):epi-T) determined at 1.85 A resolution that confirms this and reveals a unique steroid-binding mode for an AKR enzyme. Indeed, in addition to the interactions found in all other AKRs (van der Waals contacts stabilizing the core of the steroid and the hydrogen bonds established at the catalytic site by the Y55 and H117 residues with the oxygen atom of the ketone group to be reduced), m17alpha-HSD establishes with the other extremity of the steroid nucleus an additional interaction involving K31. By combining direct mutagenesis and kinetic studies, we found that the elimination of this hydrogen bond did not affect the affinity of the enzyme for its steroid substrate but led to a slight but significant increase of its catalytic efficiency (k(cat)/K(m)), suggesting a role for K31 in the release of the steroidal product at the end of the reaction. This previously unobserved steroid-binding mode for an AKR is similar to that adopted by other steroid-binding proteins, the hydroxysteroid dehydrogenases of the short-chain dehydrogenases/reductases (SDR) family and the steroid hormone nuclear receptors. Mutagenesis and structural studies made on the human type 3 3alpha-HSD, a closely related enzyme that shares 73% amino acids identity with the m17alpha-HSD, also revealed that the residue at position 24 of these two enzymes directly affects the binding and/or the release of NADPH, in addition to its role in their 17alpha/17beta stereospecificity.