Alteration of coenzyme specificity of malate dehydrogenase from Thermus flavus by site-directed mutagenesis

Enzymatic Characterization and Coenzyme Specificity Conversion of a Novel Dimeric Malate Dehydrogenase from Bacillus subtilis

Malate is an important material to various industrials and clinical applications. Bacillus subtilis is a widely used biocatalyst tool for chemical production. However, the specific enzymatic properties of malate dehydrogenase from Bacillus subtilis (BsMDH) remain largely unknown. In the present study, BsMDH was cloned, recombinantly expressed and purified to test its enzymatic properties. The molecular weight of single unit of BsMDH was 34,869.7 Da. Matrix-Assisted Laser-Desorption Ionization-Time-of-Flight Mass Spectrometry and gel filtration analysis indicated that the recombinant BsMDH could form dimers. The k_cat/K_m values of oxaloacetate and NADH were higher than those of malate and NAD⁺, respectively, indicating a better catalysis in the direction of malate synthesis than the reverse. Furthermore, six BsMDH mutants were constructed with the substitution of amino acids at the coenzyme binding site. Among them, BsMDH-T7 showed a greatly higher affinity and catalysis efficiency to NADPH than NADH with the degree of alteration of 2039, suggesting the shift of the coenzyme dependence from NADH to NADPH. In addition, BsMDH-T7 showed a relatively lower K_m value, but a higher k_cat and k_cat/K_m than NADPH-dependent MDHs from Thermus flavus and Corynebacterium glutamicum. Overall, these results indicated that BsMDH and BsMDH-T7 mutant might be promising enzymes for malate production.

Rossmann-toolbox: a deep learning-based protocol for the prediction and design of cofactor specificity in Rossmann fold proteins

The Rossmann fold enzymes are involved in essential biochemical pathways such as nucleotide and amino acid metabolism. Their functioning relies on interaction with cofactors, small nucleoside-based compounds specifically recognized by a conserved βαβ motif shared by all Rossmann fold proteins. While Rossmann methyltransferases recognize only a single cofactor type, the S-adenosylmethionine, the oxidoreductases, depending on the family, bind nicotinamide (nicotinamide adenine dinucleotide, nicotinamide adenine dinucleotide phosphate) or flavin-based (flavin adenine dinucleotide) cofactors. In this study, we showed that despite its short length, the βαβ motif unambiguously defines the specificity towards the cofactor. Following this observation, we trained two complementary deep learning models for the prediction of the cofactor specificity based on the sequence and structural features of the βαβ motif. A benchmark on two independent test sets, one containing βαβ motifs bearing no resemblance to those of the training set, and the other comprising 38 experimentally confirmed cases of rational design of the cofactor specificity, revealed the nearly perfect performance of the two methods. The Rossmann-toolbox protocols can be accessed via the webserver at https://lbs.cent.uw.edu.pl/rossmann-toolbox and are available as a Python package at https://github.com/labstructbioinf/rossmann-toolbox.

Phylogenetics and biochemistry elucidate the evolutionary link between l-malate and l-lactate dehydrogenases and disclose an intermediate group of sequences with mix functional properties

The NAD(P)-dependent malate dehydrogenases (MDH) (EC 1.1.1.37) and NAD-dependent lactate dehydrogenases (LDH) (EC. 1.1.1.27) form a large superfamily that has been characterized in organisms belonging to the three Domains of Life. MDH catalyzes the reversible conversion of the oxaloacetate into malate, while LDH operates at the late stage of glycolysis by converting pyruvate into lactate. Phylogenetic studies proposed that the LDH/MDH superfamily encompasses five main groups of enzymes. Here, starting from 16,052 reference proteomes, we reinvestigated the relationships between MDH and LDH. We showed that the LDH/MDH superfamily encompasses three main families: MDH1, MDH2, and a large family encompassing MDH3, LDH, and L-2-hydroxyisocaproate dehydrogenases (HicDH) sequences. An in-depth analysis of the phylogeny of the MDH3/LDH/HicDH family and of the nature of three important amino acids, located within the catalytic site and involved in binding and substrate discrimination, revealed a large group of sequences displaying unexpected combinations of amino acids at these three critical positions. This group branched in-between canonical MDH3 and LDH sequences. The functional characterization of several enzymes from this intermediate group disclosed a mix of functional properties, indicating that the MDH3/LDH/HicDH family is much more diverse than previously thought, and blurred the frontier between MDH3 and LDH enzymes. Present-days enzymes of the intermediate group are a valuable material to study the evolutionary steps that led to functional diversity and emergence of allosteric regulation within the LDH/MDH superfamily.

Enzyme-based amperometric biosensors for malic acid - A review

Malic acid is a key flavour component of many fruits and vegetables. There is significant interest in technologies for monitoring its concentration, particularly in winemaking. In this review we systematically and comprehensively chart progress in the development of enzyme-based amperometric biosensors for malic acid. We summarise the components and analytical parameters of malic acid sensors that have been reported over the past four decades, discussing their merits and pitfalls in terms of accuracy, sensitivity, linear range, response time and stability. We discuss how advances in electrode materials, electron mediators and the use of coupled enzymes have improved sensitivity and minimised interference, but also uncover a trade-off between sensitivity and linear range. A particular focus of our review is the three types of malate oxidoreductase enzyme that have been used in malic acid biosensors. We describe their different properties and conclude that identifying and/or engineering superior alternatives will be a key future direction for improving the commercial utility of malic acid biosensors.

Protein Engineering for Nicotinamide Coenzyme Specificity in Oxidoreductases: Attempts and Challenges

Oxidoreductases are ubiquitous enzymes that catalyze an extensive range of chemical reactions with great specificity, efficiency, and selectivity. Most oxidoreductases are nicotinamide cofactor-dependent enzymes with a strong preference for NADP or NAD. Because these coenzymes differ in stability, bioavailability and costs, the enzyme preference for a specific coenzyme is an important issue for practical applications. Different approaches for the manipulation of coenzyme specificity have been reported, with different degrees of success. Here we present various attempts for the switching of nicotinamide coenzyme preference in oxidoreductases by protein engineering. This review covers 103 enzyme engineering studies from 82 articles and evaluates the accomplishments in terms of coenzyme specificity and catalytic efficiency compared to wild type enzymes of different classes. We analyzed different protein engineering strategies and related them with the degree of success in inverting the cofactor specificity. In general, catalytic activity is compromised when coenzyme specificity is reversed, however when switching from NAD to NADP, better results are obtained. In most of the cases, rational strategies were used, predominantly with loop exchange generating the best results. In general, the tendency of removing acidic residues and incorporating basic residues is the strategy of choice when trying to change specificity from NAD to NADP, and vice versa. Computational strategies and algorithms are also covered as helpful tools to guide protein engineering strategies. This mini review aims to give a general introduction to the topic, giving an overview of tools and information to work in protein engineering for the reversal of coenzyme specificity.

In vitro Engineering of Novel Bioactivity in the Natural Enzymes

Enzymes catalyze various biochemical functions with high efficiency and specificity. In vitro design of the enzyme leads to novel bioactivity in this natural biomolecule that give answers of some vital questions like crucial residues in binding with substrate, molecular evolution, cofactor specificity etc. Enzyme engineering technology involves directed evolution, rational designing, semi-rational designing, and structure-based designing using chemical modifications. Similarly, combined computational and in vitro evolution approaches together help in artificial designing of novel bioactivity in the natural enzyme. DNA shuffling, error prone PCR and staggered extension process are used to artificially redesign active site of enzyme, which can alter its efficiency and specificity. Modifications of the enzyme can lead to the discovery of new path of molecular evolution, designing of efficient enzymes, locating active sites and crucial residues, shift in substrate, and cofactor specificity. The methods and thermodynamics of in vitro designing of the enzyme are also discussed. Similarly, engineered thermophilic and psychrophilic enzymes attain substrate specificity and activity of mesophilic enzymes that may also be beneficial for industry and therapeutics.

Function, kinetic properties, crystallization, and regulation of microbial malate dehydrogenase*

Malate dehydrogenase (MDH) is an enzyme widely distributed among living organisms and is a key protein in the central oxidative pathway. It catalyzes the interconversion between malate and oxaloacetate using NAD+ or NADP+ as a cofactor. Surprisingly, this enzyme has been extensively studied in eukaryotes but there are few reports about this enzyme in prokaryotes. It is necessary to review the relevant information to gain a better understanding of the function of this enzyme. Our review of the data generated from studies in bacteria shows much diversity in their molecular properties, including weight, oligomeric states, cofactor and substrate binding affinities, as well as differences in the direction of the enzymatic reaction. Furthermore, due to the importance of its function, the transcription and activity of this enzyme are rigorously regulated. Crystal structures of MDH from different bacterial sources led to the identification of the regions involved in substrate and cofactor binding and the residues important for the dimer-dimer interface. This structural information allows one to make direct modifications to improve the enzyme catalysis by increasing its activity, cofactor binding capacity, substrate specificity, and thermostability. A comparative analysis of the phylogenetic reconstruction of MDH reveals interesting facts about its evolutionary history, dividing this superfamily of proteins into two principle clades and establishing relationships between MDHs from different cellular compartments from archaea, bacteria, and eukaryotes.

Gradual neofunctionalization in the convergent evolution of trichomonad lactate and malate dehydrogenases

Lactate and malate dehydrogenases (LDH and MDH) are homologous, core metabolic enzymes common to nearly all living organisms. LDHs have evolved convergently from MDHs at least four times, achieving altered substrate specificity by a different mechanism each time. For instance, the LDH of anaerobic trichomonad parasites recently evolved independently from an ancestral trichomonad MDH by gene duplication. LDH plays a central role in trichomonad metabolism by catalyzing the reduction of pyruvate to lactate, thereby regenerating the NAD+ required for glycolysis. Using ancestral reconstruction methods, we identified the biochemical and evolutionary mechanisms responsible for this convergent event. The last common ancestor of these enzymes was a highly specific MDH, similar to modern trichomonad MDHs. In contrast, the LDH lineage evolved promiscuous activity by relaxing specificity in a gradual process of neofunctionalization involving one highly detrimental substitution at the “specificity residue” (R91L) and many additional mutations of small effect. L91 has different functional consequences in LDHs and in MDHs, indicating a prominent role for epistasis. Crystal structures of modern‐day and ancestral enzymes show that the evolution of substrate specificity paralleled structural changes in dimerization and α‐helix orientation. The relatively small “specificity residue” of the trichomonad LDHs can accommodate a range of substrate sizes and may permit solvent to access the active site, both of which promote substrate promiscuity. The trichomonad LDHs present a multi‐faceted counterpoint to the independent evolution of LDHs in other organisms and illustrate the diverse mechanisms by which protein function, structure, and stability coevolve.

PDB Code(s): 4UUL; 4UUM; 4UUN; 4UUO; 4UUP; 5A1T

Characterization of a thermostable glucose dehydrogenase with strict substrate specificity from a hyperthermophilic archaeon Thermoproteus sp. GDH-1

A hyperthermophilic archaeon was isolated from a terrestrial hot spring on Kodakara Island, Japan and designated as Thermoproteus sp. glucose dehydrogenase (GDH-1). Cell extracts from cells grown in medium supplemented with glucose exhibited NAD(P)-dependent glucose dehydrogenase activity. The enzyme (TgGDH) was purified and found to display a strict preference for d-glucose. The gene was cloned and expressed in Escherichia coli, resulting in the production of a soluble and active protein. Recombinant TgGDH displayed extremely high thermostability and an optimal temperature higher than 85 °C, in addition to its strict specificity for d-glucose. Despite its thermophilic nature, TgGDH still exhibited activity at 25 °C. We confirmed that the enzyme could be applied for glucose measurements at ambient temperatures, suggesting a potential of the enzyme for use in measurements in blood samples.

A computational strategy for altering an enzyme in its cofactor preference to NAD(H) and/or NADP(H)

Coenzyme engineering, especially for altered coenzyme specificity, has been a research hotspot for more than a decade. In the present study, a novel computational strategy that enhances the hydrogen-bond interaction between an enzyme and a coenzyme was developed and utilized to alter the coenzyme preference. This novel computational strategy only required the structure of the target enzyme. No other homologous enzymes were needed to achieve alteration in the coenzyme preference of a certain enzyme. Using our novel strategy, Gox2181 was reconstructed from exhibiting complete NADPH preference to exhibiting dual cofactor specificity for NADH and NADPH. Structure-guided Gox2181 mutants were designed in silico and molecular dynamics simulations were performed to evaluate the strength of hydrogen-bond interactions between the enzyme and the coenzyme NADPH. Three Gox2181 mutants displaying high structure stability and structural compatibility to NADH/NADPH were chosen for experimental confirmation. Among the three Gox2181 mutants, Gox2181-Q20R&D43S showed the highest enzymatic activity by utilizing NADPH as its coenzyme, which was even better than the wild-type enzyme. In addition, isothermal titration calorimetry analysis further verified that Gox2181-Q20R&D43S was able to interact with NADPH but the wild-type enzyme could not. This novel computational strategy represents an insightful approach for altering the cofactor preference of target enzymes.

Transhydrogenase promotes the robustness and evolvability of E. coli deficient in NADPH production

Metabolic networks revolve around few metabolites recognized by diverse enzymes and involved in myriad reactions. Though hub metabolites are considered as stepping stones to facilitate the evolutionary expansion of biochemical pathways, changes in their production or consumption often impair cellular physiology through their system-wide connections. How does metabolism endure perturbations brought immediately by pathway modification and restore hub homeostasis in the long run? To address this question we studied laboratory evolution of pathway-engineered Escherichia coli that underproduces the redox cofactor NADPH on glucose. Literature suggests multiple possibilities to restore NADPH homeostasis. Surprisingly, genetic dissection of isolates from our twelve evolved populations revealed merely two solutions: (1) modulating the expression of membrane-bound transhydrogenase (mTH) in every population; (2) simultaneously consuming glucose with acetate, an unfavored byproduct normally excreted during glucose catabolism, in two subpopulations. Notably, mTH displays broad phylogenetic distribution and has also played a predominant role in laboratory evolution of Methylobacterium extorquens deficient in NADPH production. Convergent evolution of two phylogenetically and metabolically distinct species suggests mTH as a conserved buffering mechanism that promotes the robustness and evolvability of metabolism. Moreover, adaptive diversification via evolving dual substrate consumption highlights the flexibility of physiological systems to exploit ecological opportunities.

Improving the NADH-cofactor specificity of the highly active AdhZ3 and AdhZ2 from Escherichia coli K-12

Biocatalysis is a promising tool for the sustainable production of chemicals. When cofactor depending enzymatic reactions are involved the applicability of the right cofactor is a central issue. One important example in this regard is the production of alcohols by nicotinamide cofactor (NAD(P)(+)) depending alcohol dehydrogenases. AdhZ3 from Escherichia coli, which is important for the production of alcohols from biomass, has a preference for NADPH as cofactor. We used a structure guided site-specific random approach, to change the cofactor preference towards NADH and to deduce more general rules for redesigning the cofactor specificity. Transfer of a triplet motif from NADH preferring horse liver ADH to AdhZ3 showed an insufficient switch in the preference towards NADH. A combinatorial site saturation mutagenesis altering three residues at once was applied. Library screening with two different cofactor concentrations (0.1 and 0.3mM) resulted in nine improved variants with AdhZ3-LND having the highest vmax and AdhZ3-CND having the lowest K(m). Asparagine was the most frequent amino acid found in eight of nine triplet motifs. To verify the triplet-motif, two variants of E. coli AdhZ2 DIN and LND were designed and confirmed for improved activity with NADH.

Computational design of short-chain dehydrogenase Gox2181 for altered coenzyme specificity

Short-chain dehydrogenase Gox2181 from Gluconobacter oxydans catalyzes the reduction of 2,3-pentanedione by using NADH as the physiological electron donor. To realize its synthetic biological application for coenzyme recycling use, computational design and site-directed mutagenesis have been used to engineer Gox2181 to utilize not only NADH but also NADPH as the electron donor. Single and double mutations at residues Q20 and D43 were made in a recombinant expression system that corresponded to Gox2181-D43Q and Gox2181-Q20R&D43Q, respectively. The design of mutant Q20R not only resolved the hydrogen bond interaction and electrostatic interaction between R and 2'-phosphate of NADPH, but also could enhance the binding with 2'-phophated of NADPH by combining with D43Q. Molecular dynamics simulation has been carried out to testify the hydrogen bond interactions between mutation sites and 2'-phosphate of NADPH. Steady-state turnover measurement results indicated that Gox2181-D43Q could use both NADH and NADPH as its coenzyme, and so could Gox2181-Q20R&D43Q. Meanwhile, compared to the wild-type enzyme, Gox2181-D43Q exhibited dramatically reduced enzymatic activity while Gox2181-Q20R&D43Q successfully retained the majority of enzymatic activity.

Reversal of the extreme coenzyme selectivity of Clostridium symbiosum glutamate dehydrogenase

Active-site mutants of glutamate dehydrogenase from Clostridium symbiosum have been designed and constructed and the effects on coenzyme preference evaluated by detailed kinetic measurements. The triple mutant F238S/P262S/D263K shows complete reversal in coenzyme selectivity from NAD(H) to NADP(H) with retention of high levels of catalytic activity for the new coenzyme. For oxidized coenzymes, k(cat) /K(m) ratios of the wild-type and triple mutant enzyme indicate a shift in preference of approximately 1.6 × 10(7) -fold, from ∼ 80,000-fold in favour of NAD(+) to ∼ 200-fold in favour of NADP(+). For reduced coenzymes the corresponding figure is 1.7 × 10(4) -fold, from ∼ 1000-fold in favour of NADH to ∼ 17-fold in favour of NADPH. A fourth mutation (N290G), previously identified as having a potential bearing on coenzyme specificity, did not engender any further shift in preference when incorporated into the triple mutant, despite having a significant effect when expressed as a single mutant.

An Examination by Site-Directed Mutagenesis of Putative Key Residues in the Determination of Coenzyme Specificity in Clostridial NAD-Dependent Glutamate Dehydrogenase

Sequence and structure comparisons of various glutamate dehydrogenases (GDH) and other nicotinamide nucleotide-dependent dehydrogenases have potentially implicated certain residues in coenzyme binding and discrimination. We have mutated key residues in Clostridium symbiosum NAD⁺-specific GDH to investigate their contribution to specificity and to enhance acceptance of NADPH. Comparisons with E. coli NADPH-dependent GDH prompted design of mutants F238S, P262S, and F238S/P262S, which were purified and assessed at pH 6.0, 7.0, and 8.0. They showed markedly increased catalytic efficiency with NADPH, especially at pH 8.0 (∼170-fold for P262S and F238S/P262S with relatively small changes for NADH). A positive charge introduced through the D263K mutation also greatly increased catalytic efficiency with NADPH (over 100-fold at pH 8) and slightly decreased activity with NADH. At position 242, “P6” of the “core fingerprint,” where NAD⁺- and NADP⁺-dependent enzymes normally have Gly or Ala, respectively, clostridial GDH already has Ala. Replacement with Gly produced negligible shift in coenzyme specificity.

Re-engineering the discrimination between the oxidized coenzymes NAD+ and NADP+ in clostridial glutamate dehydrogenase and a thorough reappraisal of the coenzyme specificity of the wild-type enzyme

Clostridial glutamate dehydrogenase mutants, designed to accommodate the 2'-phosphate of disfavoured NADPH, showed the expected large specificity shifts with NAD(P)H. Puzzlingly, similar assays with oxidized cofactors initially revealed little improvement with NADP(+) , although rates with NAD(+) were markedly diminished. This article reveals that the enzyme's discrimination in favour of NAD(+) and against NADP(+) had been greatly underestimated and has indeed been abated by a factor of > 16,000 by the mutagenesis. Initially, stopped-flow studies of the wild-type enzyme showed a burst increase of A(340) with NADP(+) but not NAD(+), with amplitude depending on the concentration of the coenzyme, rather than enzyme. Amplitude also varied with the commercial source of the NADP(+). FPLC, HPLC and mass spectrometry identified NAD(+) contamination ranging from 0.04 to 0.37% in different commercial samples. It is now clear that apparent rates of NADP(+) utilization mainly reflected the reduction of contaminating NAD(+), creating an entirely false view of the initial coenzyme specificity and also of the effects of mutagenesis. Purification of the NADP(+) eliminated the burst. With freshly purified NADP(+), the NAD(+) : NADP(+) activity ratio under standard conditions, previously estimated as 300 : 1, is 11,000. The catalytic efficiency ratio is even higher at 80,000. Retested with pure cofactor, mutants showed marked specificity shifts in the expected direction, for example, 16 200 fold change in catalytic efficiency ratio for the mutant F238S/P262S, confirming that the key structural determinants of specificity have been successfully identified. Of wider significance, these results underline that, without purification, even the best commercial coenzyme preparations are inadequate for such studies.

Identification and biochemical characterization of a thermostable malate dehydrogenase from the mesophile Streptomyces coelicolor A3(2)

We identified and characterized a malate dehydrogenase from Streptomyces coelicolor A3(2) (ScMDH). The molecular mass of ScMDH was 73,353.5 Da with two 36,675.0 Da subunits as analyzed by matrix-assisted laser-desorption ionization-time-of-flight mass spectrometry (MALDI-TOF-MS). The detailed kinetic parameters of recombinant ScMDH are reported here. Heat inactivation studies showed that ScMDH was more thermostable than most MDHs from other organisms, except for a few extremely thermophile bacteria. Recombinant ScMDH was highly NAD(+)-specific and displayed about 400-fold (k(cat)) and 1,050-fold (k(cat)/K(m)) preferences for oxaloacetate reduction over malate oxidation. Substrate inhibition studies showed that ScMDH activity was inhibited by excess oxaloacetate (K(i)=5.8 mM) and excess L-malate (K(i)=12.8 mM). Moreover, ScMDH activity was not affected by most metal ions, but was strongly inhibited by Fe(2+) and Zn(2+). Taken together, our findings indicate that ScMDH is significantly thermostable and presents a remarkably high catalytic efficiency for malate synthesis.

Relaxing the coenzyme specificity of 1,3-propanediol oxidoreductase from Klebsiella pneumoniae by rational design

1,3-Propanediol has wide applications for large volume markets, particularly in the polymer business. Microbial production of 1,3-propanediol has been considered as a competitor to the traditional petrochemical routes. However, the formation of 1,3-propanediol is limited by the amount of NADH supplied by the oxidative pathway of glycerol dismutation. Previous metabolic flux analysis revealed that relaxation of the coenzyme specificity of 1,3-propanediol oxidoreductase for both NADH and NADPH would increase the production of 1,3-propanediol as well as maintaining the NADH-NAD(+) circle. This work tried to accomplish such a relaxation by rational protein design. Overall binding free energy indicated that the electrostatic energy was the major force discriminating NADH from NADPH. Computational alanine-scanning mutagenesis of the active site residues illustrated that Asp41 was the key residue responsible for the coenzyme specificity. Compared with Asp41Ala, Asp41Gly could further weaken the repulsion between Asp41 and the phosphate group esterified to the 2'-hydroxyl group of the ribose at the adenine end of NADPH. Site-directed mutagenesis was conducted and the relaxation was successfully realized.

Crystal structure of archaeal highly thermostable L-aspartate dehydrogenase/NAD/citrate ternary complex

The crystal structure of the highly thermostable L-aspartate dehydrogenase (L-aspDH; EC 1.4.1.21) from the hyperthermophilic archaeon Archaeoglobus fulgidus was determined in the presence of NAD and a substrate analog, citrate. The dimeric structure of A. fulgidus L-aspDH was refined at a resolution of 1.9 A with a crystallographic R-factor of 21.7% (R(free) = 22.6%). The structure indicates that each subunit consists of two domains separated by a deep cleft containing an active site. Structural comparison of the A. fulgidus L-aspDH/NAD/citrate ternary complex and the Thermotoga maritima L-aspDH/NAD binary complex showed that A. fulgidus L-aspDH assumes a closed conformation and that a large movement of the two loops takes place during substrate binding. Like T. maritima L-aspDH, the A. fulgidus enzyme is highly thermostable. But whereas a large number of inter- and intrasubunit ion pairs are responsible for the stability of A. fulgidus L-aspDH, a large number of inter- and intrasubunit aromatic pairs stabilize the T. maritima enzyme. Thus stabilization of these two L-aspDHs appears to be achieved in different ways. This is the first detailed description of substrate and coenzyme binding to L-aspDH and of the molecular basis of the high thermostability of a hyperthermophilic L-aspDH.

Structural basis for the alteration of coenzyme specificity in a malate dehydrogenase mutant

To elucidate the structural basis for the alteration of coenzyme specificity from NADH toward NADPH in a malate dehydrogenase mutant EX7 from Thermus flavus, we determined the crystal structures at 2.0 A resolution of EX7 complexed with NADPH and NADH, respectively. In the EX7-NADPH complex, Ser42 and Ser45 form hydrogen bonds with the 2'-phosphate group of the adenine ribose of NADPH, although the adenine moiety is not seen in the electron density map. In contrast, although Ser42 and Ser45 occupy a similar position in the EX7-NADH complex structure, both the adenine and adenine ribose moieties of NADH are missing in the map. These results and kinetic analysis of site-directed mutant enzymes indicate (1) that the preference of EX7 for NADPH over NADH is ascribed to the recognition of the 2'-phosphate group by two Ser and Arg44, and (2) that the adenine moiety of NADPH is not recognized in this mutant.

Complete reversal of coenzyme specificity of isocitrate dehydrogenase from Haloferax volcanii

Haloferax volcanii Ds-threo-isocitrate dehydrogenase (ICDH) was highly expressed in bacteria as inclusion bodies. The recombinant enzyme was refolded, purified and characterized, and was found to be NADP-dependent like the wild-type protein. Sequence alignment of several isocitrate dehydrogenases from evolutionarily divergent organisms including H. volcanii revealed that the amino acid residues involved in coenzyme specificity are highly conserved. Our objective was to switch the coenzyme specificity of halophilic ICDH by altering these conserved amino acids. We were able to switch coenzyme specificity from NADP+ to NAD+ by changing five amino acids by site-directed mutagenesis (Arg291, Lys343, Tyr344, Val350 and Tyr390). The five mutants of ICDH were overexpressed in Escherichia coli as inclusion bodies and each recombinant ICDH protein was refolded and purified, and its kinetic parameters were determined. Coenzyme specificity did not switch until all five amino acids were substituted.

Crystal structure of NAD-dependent malate dehydrogenase complexed with NADP(H)

For better understanding of the coenzyme specificity in NAD-dependent MDH (tMDH) from Thermus flavus AT-62, we determined the crystal structures of tMDH-NADP(H) complex at maximally 1.65 A resolution. The overall structure is almost the same as that of the tMDH-NADH complex. However, NADP(H) binds to tMDH in the reverse orientation, where adenine occupies the position near the catalytic center and nicotinamide is positioned at the adenine binding site of the tMDH-NADH complex. Consistent with this, kinetic analysis of the malate-oxidizing reaction revealed that NADP(+) inhibited tMDH at high concentrations. This has provided the first evidence for the alternative binding mode of the nicotinamide coenzyme, that has pseudo-symmetry in its structure, in a single enzyme.

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.

The first crystal structure of hyperthermostable NAD-dependent glutamate dehydrogenase from Pyrobaculum islandicum

The extremely thermostable NAD-dependent glutamate dehydrogenase (NAD-GluDH) from Pyrobaculum islandicum, a member of the Crenarchaeota, was crystallized, and its 3D structure has been determined by X-ray diffraction methods. The homohexameric structure of Pb. islandicum glutamate dehydrogenase (Pis-GluDH) was solved and refined at a resolution of 2.9A with a crystallographic R-factor of 19.9% (Rfree 26.0%). The structure indicates that each subunit consists of two domains separated by a deep cleft containing an active site. The secondary structural elements and catalytically important residues of the enzyme were highly conserved among the NAD(P)-dependent GluDHs from other sources. A structural comparison of Pis-GluDH with other NAD(P)-dependent GluDHs suggests that a significant difference in the alpha8-loop-alpha9 region of this enzyme is associated with its coenzyme specificity. From the analysis of the 3D structure, hydrophobic interactions between intersubunits were found to be important features for the enzyme oligomerization. It has been reported that Pis-GluDH is highly thermostable, like the GluDH of the hyperthermophilic archaeum Pyrococcus furiosus, and the increase in the intersubunit ion pair networks is responsible for the extreme thermostability of the Pc. furiosus enzyme. However, the number of intersubunit ion pairs in the Pis-GluDH molecules is much smaller than those of the Pc. furiosus GluDH. The number of hydrophobic interactions at the intersubunit interfaces were increased and responsible for the extremely high thermostability. This indicates that the major molecular strategy for high thermostability of the GluDHs may be different for each hyperthermophile.

Complete reversal of coenzyme specificity of xylitol dehydrogenase and increase of thermostability by the introduction of structural zinc

Pichia stipitis NAD(+)-dependent xylitol dehydrogenase (XDH), a medium-chain dehydrogenase/reductase, is one of the key enzymes in ethanol fermentation from xylose. For the construction of an efficient biomass-ethanol conversion system, we focused on the two areas of XDH, 1) change of coenzyme specificity from NAD(+) to NADP(+) and 2) thermostabilization by introducing an additional zinc atom. Site-directed mutagenesis was used to examine the roles of Asp(207), Ile(208), Phe(209), and Asn(211) in the discrimination between NAD(+) and NADP(+). Single mutants (D207A, I208R, F209S, and N211R) improved 5 approximately 48-fold in catalytic efficiency (k(cat)/K(m)) with NADP(+) compared with the wild type but retained substantial activity with NAD(+). The double mutants (D207A/I208R and D207A/F209S) improved by 3 orders of magnitude in k(cat)/K(m) with NADP(+), but they still preferred NAD(+) to NADP(+). The triple mutant (D207A/I208R/F209S) and quadruple mutant (D207A/I208R/F209S/N211R) showed more than 4500-fold higher values in k(cat)/K(m) with NADP(+) than the wild-type enzyme, reaching values comparable with k(cat)/K(m) with NAD(+) of the wild-type enzyme. Because most NADP(+)-dependent XDH mutants constructed in this study decreased the thermostability compared with the wild-type enzyme, we attempted to improve the thermostability of XDH mutants by the introduction of an additional zinc atom. The introduction of three cysteine residues in wild-type XDH gave an additional zinc-binding site and improved the thermostability. The introduction of this mutation in D207A/I208R/F209S and D207A/I208R/F209S/N211R mutants increased the thermostability and further increased the catalytic activity with NADP(+).

Relaxing the nicotinamide cofactor specificity of phosphite dehydrogenase by rational design

Homology modeling was used to identify two particular residues, Glu175 and Ala176, in Pseudomonas stutzeri phosphite dehydrogenase (PTDH) as the principal determinants of nicotinamide cofactor (NAD(+) and NADP(+)) specificity. Replacement of these two residues by site-directed mutagenesis with Ala175 and Arg176 both separately and in combination resulted in PTDH mutants with relaxed cofactor specificity. All three mutants exhibited significantly better catalytic efficiency for both cofactors, with the best kinetic parameters displayed by the double mutant, which had a 3.6-fold higher catalytic efficiency for NAD(+) and a 1000-fold higher efficiency for NADP(+). The cofactor specificity was changed from 100-fold in favor of NAD(+) for the wild-type enzyme to 3-fold in favor of NADP(+) for the double mutant. Isoelectric focusing of the proteins in a nondenaturing gel showed that the replacement with more basic residues indeed changed the effective pI of the protein. HPLC analysis of the enzymatic products of the double mutant verified that the reaction proceeded to completion using either substrate and produced only the corresponding reduced cofactor and phosphate. Thermal inactivation studies showed that the double mutant was protected from thermal inactivation by both cofactors, while the wild-type enzyme was protected by only NAD(+). The combined results provide clear evidence that Glu175 and Ala176 are both critical for nicotinamide cofactor specificity. The rationally designed double mutant might be useful for the development of an efficient in vitro NAD(P)H regeneration system for reductive biocatalysis.

Engineering and characterization of a NADPH-utilizing cytochrome b5 reductase

Microsomal cytochrome b(5) reductase (EC 1.6.2.2) catalyzes the reduction of ferricytochrome b(5) using NADH as the physiological electron donor. Site-directed mutagenesis has been used to engineer the soluble rat cytochrome b(5) reductase diaphorase domain to utilize NADPH as the preferred electron donor. Single and double mutations at residues D239 and F251 were made in a recombinant expression system that corresponded to D239E, S and T, F251R, and Y, D239S/F251R, D239S/F251Y, and D239T/F251R, respectively. Steady-state turnover measurements indicated that D239S/F251Y was bispecific while D239T, D239S/F251R, and D239T/F251R were each NADPH-specific. Wild-type (WT) cytochrome b(5) reductase showed a 3700-fold preference for NADH whereas the mutant with the highest NADPH efficiency, D239T, showed an 11-fold preference for NADPH, a 39200-fold increase. Wild-type cytochrome b(5) reductase only formed a stable charge-transfer complex with NADH while D239T formed complexes with both NADH and NADPH. The rates of hydride ion transfer, determined by stopped-flow kinetics, were k(NADH-WT) = 130 s(-1), k(NADPH-WT) = 5 s(-1), k(NADH-D239T) = 180 s(-1), and k(NADPH-D239T) = 73 s(-1). K(s) determinations by differential spectroscopy demonstrated that D239T could bind nonreducing pyridine nucleotides with a phosphate or a hydroxyl substituent at the 2' position, whereas wild-type cytochrome b(5) reductase would only bind 2' hydroxylated molecules. Oxidation-reduction potentials (E degrees ', n = 2) for the flavin cofactor were WT = -268 mV, D239T = -272 mV, WT+NAD(+) = -190 mV, D239T+NAD(+) = -206 mV, WT+NADP(+) = -253 mV, and D239T+NADP(+) = -215 mV, which demonstrated the thermodynamic contribution of NADP(+) binding to D239T. The crystal structures of D239T and D239T in complex with NAD(+) indicated that the loss of the negative electrostatic surface that precluded 2' phosphate binding in the wild-type enzyme was primarily responsible for the observed improvement in the use of NADPH by the D239T mutant.

The B' helix determines cytochrome P450nor specificity for the electron donors NADH and NADPH

Nitric oxide reductase (Nor) cytochrome P450nor (P450nor) is unique because it is catalytically self-sufficient, receiving electrons directly from NADH or NADPH. However, little is known about the direct binding of NADH to cytochrome. Here, we report that oxidized pyridine nucleotides (NAD(+) and NADP(+)) and an analogue induce a spectral perturbation in bound heme when mixed with P450nor. The P450nor isoforms are classified according to electron donor specificity for NADH or NADPH. One type (Fnor, a P450nor of Fusarium oxysporum) utilizes only NADH. We found that NAD(+) induced a type I spectral change in Fnor, whereas NADP(+) induced a reverse type I spectral change, although the K(d) values for both were comparable. In contrast, NADP(+) as well as NAD(+) caused a type I spectral change in Tnor, a P450nor isozyme from Trichosporon cutaneum that utilizes both NADH and NADPH as electron donors. The B' helix region of Tnor ((73)SAGGKAAA(80)) contains some Ala and Gly residues, whereas the sequence is replaced at a few sites with more bulky amino acid residues in Fnor ((73)SASGKQAA(80)). A single mutation (S75G) significantly improved the NADPH- dependent Nor activity of Fnor, and the overall activity was accelerated via the NADPH-enhanced reduction step. These results showed that pyridine nucleotide cofactors can bind P450nor and that only a few residues in the B' helix region determine cofactor specificity. We further showed that a poor electron donor (NADPH) could also bind Fnor, but an appropriate configuration for electron transfer is blocked by steric hindrance mainly by Ser(75) against the 2'-phosphate moiety. The present results along with previous observations together revealed a novel motif for cofactor binding.

The commercial production of chemicals using pathway engineering

Integration of metabolic pathway engineering and fermentation production technologies is necessary for the successful commercial production of chemicals. The 'toolbox' to do pathway engineering is ever expanding to enable mining of biodiversity, to maximize productivity, enhance carbon efficiency, improve product purity, expand product lines, and broaden markets. Functional genomics, proteomics, fluxomics, and physiomics are complementary to pathway engineering, and their successful applications are bound to multiply product turnover per cell, channel carbon efficiently, shrink the size of factories (i.e., reduce steel in the ground), and minimize product development cycle times to bring products to market.

A highly specific monomeric isocitrate dehydrogenase from Corynebacterium glutamicum

The monomeric isocitrate dehydrogenase (IDH) of Corynebacterium glutamicum is compared to the topologically distinct dimeric IDH of Escherichia coli. Both IDHs have evolved to efficiently catalyze identical reactions with similar pH optimum as well as striking specificity toward NADP and isocitrate. However, the monomeric IDH is 10-fold more active (calculated as kcat/Km.isocitrate/Km.NADP) and 7-fold more NADP-specific than the dimeric enzyme, favoring NADP over NAD by a factor of 50,000. Such an extraordinary coenzyme specificity is not rivaled by any other characterized dehydrogenases. In addition, the monomeric enzyme is 10-fold more specific for isocitrate. The spectacular substrate specificity may be predominantly attributed to the isocitrate-assisted stabilization of catalytic complex during hydride transfer. No significant overall sequence identity is found between the monomeric and dimeric enzymes. However, structure-based alignment leads to the identification of three regions in the monomeric enzyme that match closely the three motifs located in the central region of dimeric IDHs and the homologous isopropylmalate dehydrogenases. The role of Lys253 as catalytic residue has been demonstrated by site-directed mutagenesis. Our results suggest that monomeric and dimeric forms of IDHs are functionally and structurally homologous.

Cloning and characterization of ftsZ and pyrF from the archaeon Thermoplasma acidophilum

To characterize cytoskeletal components of archaea, the ftsZ gene from Thermoplasma acidophilum was cloned and sequenced. In T. acidophilum ftsZ, which is involved in cell division, was found to be in an operon with the pyrF gene, which encodes orotidine-5'-monophosphate decarboxylase (ODC), an essential enzyme in pyrimidine biosynthesis. Both ftsZ and pyrF from T. acidophilum were expressed in Escherichia coli and formed functional proteins. FtsZ expression in wild-type E. coli resulted in the filamentous phenotype characteristic of ftsZ mutants. T. acidophilum pyrF expression in an E. coli mutant lacking pyrF complemented the mutation and rescued the strain. Sequence alignments of ODCs from archaea, bacteria, and eukarya reveal five conserved regions, two of which have homology to 3-hexulose-6-phosphate synthase (HPS), suggesting a common substrate recognition and binding motif.

Inhibition of the thioredoxin-dependent activation of the NADP-malate dehydrogenase and cofactor specificity

The chloroplastic NADP-malate dehydrogenase is activated by reduction of its N- and C-terminal disulfides by reduced thioredoxin. The activation is inhibited by NADP(+), the oxidized form of the cofactor. Previous studies suggested that the C-terminal disulfide was involved in this process. Recent structural data pointed toward a possible direct interaction between the C terminus of the oxidized enzyme and the cofactor. In the present study, the relationship between the cofactor specificity for catalysis and for inhibition of activation has been investigated by changing the cofactor specificity of the enzyme by substitution of selected residues of the cofactor-binding site. An NAD-specific thiol-regulated MDH was engineered. Its activation was inhibited by NAD(+) but no longer by NADP(+). These results demonstrate that the oxidized cofactor is bound at the same site as the reduced cofactor and support the idea of a direct interaction between the negatively charged C-terminal end of the enzyme and the positively charged nicotinamide ring of the cofactor, in agreement with the structural data. The structural requirements for cofactor specificity are modeled and discussed.

Competition between C-terminal tyrosine and nicotinamide modulates pyridine nucleotide affinity and specificity in plant ferredoxin-NADP(+) reductase

Chloroplast ferredoxin-NADP(+) reductase has a 32,000-fold preference for NADPH over NADH, consistent with its main physiological role of NADP(+) photoreduction for de novo carbohydrate biosynthesis. Although it is distant from the 2'-phosphoryl group of NADP(+), replacement of the C-terminal tyrosine (Tyr(308) in the pea enzyme) by Trp, Phe, Gly, and Ser produced enzyme forms in which the preference for NADPH over NADH was decreased about 2-, 10-, 300-, and 400-fold, respectively. Remarkably, in the case of the Y308S mutant, the k(cat) value for the NADH-dependent activity approached that of the NADPH-dependent activity of the wild-type enzyme. Furthermore, difference spectra of the NAD(+) complexes revealed that the nicotinamide ring of NAD(+) binds at nearly full occupancy in the active site of both the Y308G and Y308S mutants. These results correlate well with the k(cat) values obtained with these mutants in the NADH-ferricyanide reaction. The data presented support the hypothesis that specific recognition of the 2'-phosphate group of NADP(H) is required but not sufficient to ensure a high degree of discrimination against NAD(H) in ferredoxin-NADP(+) reductase. Thus, the C-terminal tyrosine enhances the specificity of the reductase for NADP(H) by destabilizing the interaction of a moiety common to both coenzymes, i.e. the nicotinamide.

Cloning, overexpression, and mutagenesis of the Sporobolomyces salmonicolor AKU4429 gene encoding a new aldehyde reductase, which catalyzes the stereoselective reduction of ethyl 4-chloro-3-oxobutanoate to ethyl (S)-4-chloro-3-hydroxybutanoate

We cloned and sequenced the gene encoding an NADPH-dependent aldehyde reductase (ARII) in Sporobolomyces salmonicolor AKU4429, which reduces ethyl 4-chloro-3-oxobutanoate (4-COBE) to ethyl (S)-4-chloro-3-hydroxybutanoate. The ARII gene is 1,032 bp long, is interrupted by four introns, and encodes a 37,315-Da polypeptide. The deduced amino acid sequence exhibited significant levels of similarity to the amino acid sequences of members of the mammalian 3beta-hydroxysteroid dehydrogenase-plant dihydroflavonol 4-reductase superfamily but not to the amino acid sequences of members of the aldo-keto reductase superfamily or to the amino acid sequence of an aldehyde reductase previously isolated from the same organism (K. Kita, K. Matsuzaki, T. Hashimoto, H. Yanase, N. Kato, M. C.-M. Chung, M. Kataoka, and S. Shimizu, Appl. Environ. Microbiol. 62:2303-2310, 1996). The ARII protein was overproduced in Escherichia coli about 2, 000-fold compared to the production in the original yeast cells. The enzyme expressed in E. coli was purified to homogeneity and had the same catalytic properties as ARII purified from S. salmonicolor. To examine the contribution of the dinucleotide-binding motif G(19)-X-X-G(22)-X-X-A(25), which is located in the N-terminal region, during ARII catalysis, we replaced three amino acid residues in the motif and purified the resulting mutant enzymes. Substrate inhibition of the G(19)-->A and G(22)-->A mutant enzymes by 4-COBE did not occur. The A(25)-->G mutant enzyme could reduce 4-COBE when NADPH was replaced by an equimolar concentration of NADH.

Probing the affinity and specificity of yeast alcohol dehydrogenase I for coenzymes

Yeast (Saccharomyces cerevisiae) alcohol dehydrogenase I (SceADH) binds NAD+ and NADH less tightly and turns over substrates more rapidly than does horse (Equus caballus) liver alcohol dehydrogenase E isoenzyme (EcaADH), and neither enzyme uses NADP efficiently. Amino acid residues in the proposed adenylate binding pocket of SceADH were substituted in attempts to improve affinity for coenzymes or reactivity with NADP. Substitutions in SceADH (Gly202Ile or Ser246Ile) with the corresponding residues in the adenine binding site of the homologous EcaADH have modest effects on coenzyme binding and other kinetic constants, but the Ser246Ile substitution decreases turnover numbers by 350-fold. The Ser176Phe substitution (also near adenine site) significantly decreases affinity for coenzymes and turnover numbers. In the consensus nucleotide-binding betaalphabeta fold sequence, SceADH has two alanine residues (177-GAAGGLG-183) instead of the Leu200 in EcaADH (199-GLGGVG-204); the Ala178-Ala179 to Leu substitution significantly decreases affinity for coenzymes and turnover numbers. Some NADP-dependent enzymes have an Ala corresponding to Gly183 in SceADH; the Gly183Ala substitution significantly decreases affinity for coenzymes and turnover numbers. NADP-dependent enzymes usually have a neutral residue instead of the Asp (Asp201 in SceADH) that interacts with the hydroxyl groups of the adenosine ribose, along with a basic residue (at position 202 or 203) to stabilize the 2'-phosphate of NADP. The Gly203Arg change in SceADH does not significantly affect the kinetics. The Gly183Ala or Gly203Arg substitutions do not enable SceADH to use NADP+ as coenzyme. SceADH with the single Asp201Gly or double Asp201Gly:Gly203Arg substitutions have similar, low activity with NADP+. The results suggest that several of the amino acid residues participate in coenzyme binding and that conversion of specificity for coenzyme requires multiple substitutions.

A crystallographic comparison between mutated glyceraldehyde-3-phosphate dehydrogenases from Bacillus stearothermophilus complexed with either NAD+ or NADP+

Mutations have been introduced in the cytosolic glyceraldehyde-3-phosphate dehydrogenase (GAPDH) from Bacillus stearothermophilus in order to convert its cofactor selectivity from a specificity towards NAD into a preference for NADP. In the B-S mutant, five mutations (L33T, T34G, D35G, L187A, P188S) were selected on the basis of a sequence alignment with NADP-dependent chloroplastic GAPDHs. In the D32G-S mutant, two of the five mutations mentioned above (L187A, P188S) have been used in combination with another one designed from electrostatic considerations (D32G). Both mutants exhibit a dual-cofactor selectivity at the advantage of either NAD (B-S) or NADP (D32G-S). In order to analyse the cofactor-binding site plasticity at the molecular level, crystal structures of these mutants have been solved, when complexed with either NAD+ (D32G-Sn, resolution 2.5 A, R = 13.9%; B-Sn, 2.45 A, 19.3%) or NADP+ (D32G-Sp, 2.2 A, 19.2%; B-Sp, 2.5 A, 14.4%). The four refined models are very similar to that of the wild-type GAPDH and as expected resemble more closely the holo form than the apo form. In the B-S mutant, the wild-type low affinity for NADP+ seems to be essentially retained because of repulsive electrostatic contacts between the extra 2'-phosphate and the unchanged carboxylate group of residue D32. Such an antideterminant effect is not well compensated by putative attractive interactions which had been expected to arise from the newly-introduced side-chains. In this mutant, recognition of NAD+ is slightly affected with respect to that known on the wild-type, because mutations only weakly destabilize hydrogen bonds and van der Waals contacts originally present in the natural enzyme. Thus, the B-S mutant does not mimic efficiently the chloroplastic GAPDHs, and long-range and/or second-layer effects, not easily predictable from visual inspection of three-dimensional structures, need to be taken into account for designing a true "chloroplastic-like" mutant of cytosolic GAPDH. In the case of the D32G-S mutant, the dissociation constants for NAD+ and NADP+ are practically reversed with respect to those of the wild-type. The strong alteration of the affinity for NAD+ obviously proceeds from the suppression of the two wild-type hydrogen bonds between the adenosine 2'- and 3'-hydroxyl positions and the D32 carboxylate group. As expected, the efficient recognition of NADP+ is partly promoted by the removal of intra-subunit electrostatic repulsion (D32G) and inter-subunit steric hindrance (L187A, P188S). Another interesting feature of the reshaped NADP+-binding site is provided by the local stabilization of the extra 2'-phosphate which forms a hydrogen bond with the side-chain hydroxyl group of the newly-introduced S188. When compared to the presently known natural NADP-binding clefts, this result clearly demonstrates that an absolute need for a salt-bridge involving the 2'-phosphate is not required to switch the cofactor selectivity from NAD to NADP. In fact, as it is the case in this mutant, only a moderately polar hydrogen bond can be sufficient to make the extra 2'-phosphate of NADP+ well recognized by a protein environment.

Malate dehydrogenase from the mesophile Chlorobium vibrioforme and from the mild thermophile Chlorobium tepidum: molecular cloning, construction of a hybrid, and expression in Escherichia coli

The genes (mdh) encoding malate dehydrogenase (MDH) from the mesophile Chlorobium vibrioforme and the moderate thermophile C. tepidum were cloned and sequenced, and the complete amino acid sequences were deduced. When the region upstream of mdh was analyzed, a sequence with high homology to an operon encoding ribosomal proteins from Escherichia coli was found. Each mdh gene consists of a 930-bp open reading frame and encodes 310 amino acid residues, corresponding to a subunit weight of 33,200 Da for the dimeric enzyme. The amino acid sequence identity of the two MDHs is 86%. Homology searches using the primary structures of the two MDHs revealed significant sequence similarity to lactate dehydrogenases. A hybrid mdh was constructed from the 3' part of mdh from C. tepidum and the 5' part of mdh from C. vibrioforme. The thermostabilities of the hybrid enzyme and of MDH from C. vibrioforme and C. tepidum were compared.

Structural determinants of nucleotide coenzyme specificity in the distinctive dinucleotide binding fold of HMG-CoA reductase from Pseudomonas mevalonii

The 102-residue small domain of the 428-residue NAD(H)-dependent HMG-CoA reductase of Pseudomonas mevalonii (EC 1.1.1.88) binds NAD(H) at a distinctive, non-Rossmann dinucleotide binding fold. The three-dimensional structure reveals that Asp146 lies close to the 2'-OH of NAD-. To investigate the role of this residue in determination of coenzyme specificity, Asp146 was mutated to Ala, Gly, Ser, and Asn. The mutant enzymes were analyzed for their ability to catalyze the oxidative acylation of mevalonate to HMG-CoA using either the natural coenzyme NAD+ or the alternate coenzyme NADP+. Mutation of Asp146 to Ala or Gly increased the specificity for NADP+, expressed as the ratio of kcat/K(m) for NADP+ to kcat/K(m) for NAD+, 1200-fold (enzyme D146G) and 6700-fold (enzyme D146A). Mutation of Asp146 was accompanied by 565-fold (D146G) and 330-fold (D146A) increases in kcat/K(m) for NADP+ and 2-fold (D146G) and 20-fold (D146A) decreases in kcat/K(m) for NAD+. To further improve NADP+ specificity, Gln147, Leu148, Leu149, or Thr192 of enzyme D146G or D146A was replaced by lysine or arginine, which could stabilize the 2'-phosphate of NADP+. Enzymes D146G/T192K, D146G/T192R, D146G/L148K, D146A/L148K, and D146A/L148R exhibited 3200-, 4500-, 56000-, 72000-, and 83000-fold increases in the specificity for NADP+ relative to the wild-type enzyme.

Phosphorus-31 nuclear magnetic resonance studies on coenzyme binding and specificity in glyceraldehyde-3-phosphate dehydrogenase

Binding of NAD(P)+ to wild type and a series of mutants of the glycolytic NAD-dependent glyceraldehyde-3-phosphate dehydrogenase (GAPDH) from Bacillus stearothermophilus designed to alter the cofactor specificity [Clermont, S., Corbier, C., Mely, Y., Gerard, D., Wonacott, A., & Branlant, G. (1993) Biochemistry 21, 10178-10184] has been studied by 31P NMR. In the mutants with the L187A and P188S substitutions, the pyrophosphate signals are split, and the upfield resonance has been assigned to the P(a) phosphate. Titration of the NADP+ 2'-phosphate pKa deduced from its chemical shift shows that the electrostatic environment in the binding site is largely affected by the single point mutations. pKas ranging from 7.7 for the L187A-P188S mutant to < 5.7 for the D32G-L187A-P188S and D32A-L187A-P188S mutants have been observed, thus indicating that the binding of NADP+ is modulated by the ionization state of its 2'-phosphate. In the quintuple mutant L33T-T34G-D35G-L187A-P188S, designed in comparison with the photosynthetic NAD(P)-dependent GAPDH of the chloroplast, the 2'-phosphate has a pKa of 6.8. As further stabilizing interactions like hydrogen bonds or positively charged side chains would lower this pKa, it is suggested that the 2'-phosphate ionization state of bound NADP+ in chloroplastic GAPDH is dianionic. The NADP+ dissociation rate constants (k(off)) of the three mutants D32G, L187A-P188S, and D32G-L187A-P188S, are higher at pH 6.1 than at pH 8.1 and are similar at the same pH, indicating that the difference in binding affinity between these three mutants results from the molecular recognition step or conformational change upon binding.

Determinants of cofactor specificity in isocitrate dehydrogenase: structure of an engineered NADP+ --> NAD+ specificity-reversal mutant

The 7-fold mutation Cys201Met/Cys332Tyr/Lys344Asp/Tyr345Ile/Val35 1Ala/Tyr391Lys/Arg395Ser converts the cofactor specificity of Escherichia coli isocitrate dehydrogenase from a 7000-fold preference for NADP+ to a 200-fold preference for NAD+, with overall activity comparable to that of wild-type NAD+-dependent isocitrate dehydrogenases. The structure of the NAD+-dependent mutant has been determined and refined to a working R-factor of 0.186 at 1.9 A resolution. The structure shows that NADP+ affinity is destroyed by removing favorable interactions between the 2'-phosphate and Tyr345, Tyr391, and Arg395 and by adding a repulsive interaction with Asp344. NAD+ affinity is enhanced by adding hydrogen bonds between Asp344 and the free 2'-hydroxyl. The favorable Asp344-2'-OH interaction requires a change in the pucker of the ribose to C2' endo and a shift in the adenine ring. The ring shift is made possible by a series of changes in steric packing interactions. The linchpin for repacking in the adenosine binding site is residue 351. The side chain of this "second layer" residue dictates packing of the surrounding "first layer" residues which interact with the 2' moiety and, in turn, directly determine specificity.

Site-directed mutagenesis of a serine residue in cinnamyl alcohol dehydrogenase, a plant NADPH-dependent dehydrogenase, affects the specificity for the coenzyme

Using recombinant cinnamyl alcohol dehydrogenase isoform 2 (CAD2, EC 1.1.1.195), an NADPH-dependent aromatic alcohol dehydrogenase involved in lignification in vascular plants, we have investigated the detailed steady-state kinetic mechanism of CAD2 and the role of a serine residue in determining the cofactor specificity of CAD2. Site-directed mutagenesis (S212D) and overexpression of the WT and mutant S212D forms of CAD2 in Escherichia coli, followed by kinetic studies on the purified WT and mutant proteins, confirmed the involvement of S212D in recognizing the phosphate group of NADPH and provided information on the structural requirements for NADPH specificity. From substrate kinetic patterns and product inhibition studies both WT and S212D mutant forms of CAD2 have been shown to follow rapid equilibrium random bireactant kinetics with the value of the interaction factor (alpha) for WT (0.25) being significantly less than that for S212D CAD2 (0.45). The changes in binding energy arising from the mutation on the binding of the 2'-phosphate site of the coenzyme were assessed. A marked degree of physical interaction was detected between the enzymatic binding sites of the coniferyl alcohol substrate and the 2'-phosphate binding region, which are quite distant in the three-dimensional structure. The inhibition by 2',5'-ADP and 5'-AMP was found to be weak for both WT and S212D CAD2. Strong substrate inhibition was detected for CAD2, and its implications for plant physiological studies were assessed.(ABSTRACT TRUNCATED AT 250 WORDS)

Biological production of semisynthetic opiates using genetically engineered bacteria

Semisynthetic derivatives of morphine and related alkaloids are in widespread clinical use. Due to the complexity of these molecules, however, chemical transformations are difficult to achieve in high yields. We recently identified the powerful analgesic hydromorphone as an intermediate in the metabolism of morphine by Pseudomonas putida M10. Here we describe the construction of recombinant strains of Escherichia coli that express morphine dehydrogenase and morphinone reductase. These strains are capable of efficiently transforming the naturally occurring alkaloids morphine and codeine to hydromorphone and the antitussive hydrocodone, respectively. Our results demonstrate the potential for recombinant DNA technology to provide biological routes for the synthesis of known and novel semisynthetic opiate drugs.

Three-dimensional structure prediction of the NAD binding site of proton-pumping transhydrogenase from Escherichia coli

A three-dimensional structure of the NAD site of Escherichia coli transhydrogenase has been predicted. The model is based on analysis of conserved residues among the transhydrogenases from five different sources, homologies with enzymes using NAD as cofactors or substrates, hydrophilicity profiles, and secondary structure predictions. The present model supports the hypothesis that there is one binding site, located relatively close to the N-terminus of the alpha-subunit. The proposed structure spans residues alpha 145 to alpha 287, and it includes five beta-strands and five alpha-helices oriented in a typical open twisted alpha/beta conformation. The amino acid sequence following the GXGXXG dinucleotide binding consensus sequence (residues alpha 172 to alpha 177) correlates exactly to a typical fingerprint region for ADP binding beta alpha beta folds in dinucleotide binding enzymes. In the model, aspartic acid alpha 195 forms hydrogen bonds to one or both hydroxyl groups on the adenosine ribose sugar moiety. Threonine alpha 196 and alanine alpha 256, located at the end of beta B and beta D, respectively, create a hydrophobic sandwich with the adenine part of NAD buried inside. The nicotinamide part is located in a hydrophobic cleft between alpha A and beta E. Mutagenesis work has been carried out in order to test the predicted model and to determine whether residues within this domain are important for proton pumping directly. All data support the predicted structure, and no residue crucial for proton pumping was detected. Since no three-dimensional structure of transhydrogenase has been solved, a well based tertiary structure prediction is of great value for further experimental design in trying to elucidate the mechanism of the energy-linked proton pump.

Malate dehydrogenase: a model for structure, evolution, and catalysis

Malate dehydrogenases are widely distributed and alignment of the amino acid sequences show that the enzyme has diverged into 2 main phylogenetic groups. Multiple amino acid sequence alignments of malate dehydrogenases also show that there is a low degree of primary structural similarity, apart from in several positions crucial for nucleotide binding, catalysis, and the subunit interface. The 3-dimensional structures of several malate dehydrogenases are similar, despite their low amino acid sequence identity. The coenzyme specificity of malate dehydrogenase may be modulated by substitution of a single residue, as can the substrate specificity. The mechanism of catalysis of malate dehydrogenase is similar to that of lactate dehydrogenase, an enzyme with which it shares a similar 3-dimensional structure. Substitution of a single amino acid residue of a lactate dehydrogenase changes the enzyme specificity to that of a malate dehydrogenase, but a similar substitution in a malate dehydrogenase resulted in relaxation of the high degree of specificity for oxaloacetate. Knowledge of the 3-dimensional structures of malate and lactate dehydrogenases allows the redesign of enzymes by rational rather than random mutation and may have important commercial implications.

Chemical modification of arginine and lysine residues in coenzyme-binding domain of carbonyl reductase from rabbit kidney: indomethacin affords a significant protection against inactivation of the enzyme by phenylglyoxal

Carbonyl reductase from rabbit kidney was inactivated by phenylglyoxal (PGO) and 2,4,6-trinitrobenzenesulfonate sodium (TNBS). NADP+ protected the enzyme from the inactivations by PGO and TNBS, suggesting that essential arginine and lysine residues are located in coenzyme-binding domain of the enzyme. Judging from the effects of PGO-treated enzymes in the presence and in the absence of NADP+ on the fluorescence intensity of NADPH, one essential arginine residue in coenzyme-binding domain was found to have a role in the binding of NADPH to the enzyme. Indomethacin afforded a significant protection against inactivation of the enzyme by PGO, whereas it could not protect the enzyme from the inactivation by TNBS. It is reasonable to postulate that indomethacin interacts at least in part with or near one essential arginine residue in coenzyme-binding domain of carbonyl reductase from rabbit kidney.