Cofactor Preference Research Articles

New insights into the substrate inhibition of human 17β-hydroxysteroid dehydrogenase type 1

Human type 1 17β-hydroxysteroid dehydrogenase (17β-HSD1),a member of the short-chain dehydrogenase/reductase family, catalyzes the last step in the bioactivation of the most potent estrogen estradiol with high specificity and is thus involved in estrogen-dependent diseases. As an oxidoreductase, 17β-HSD1 can utilize both triphosphate and diphosphate cofactors in reaction at the molecular level, but more specific with triphosphate cofactor. The NADPH is much higher than NADP+ in living cells leading to preliminary reduction action. The enzyme also showed substrate-induced inhibition unprecedented in other members of 17β-HSDs. Our previous study elucidated the structural mechanism of substrate inhibition is due to the reversely bound estrone (E1) in the substrate-binding pocket of the enzyme resulting in a dead-end complex. However, the effect of the cofactor preference on the substrate inhibition of the enzyme is not yet clear. In the present study, we solved the ternary crystal structures of 17β-HSD1 in complex with E1 and cofactor analog NAD+ . Combined with molecular dynamics simulation using the enzyme with NADH/NADPH and different oriented E1 (normally oriented, E1N; reversely oriented, E1R), such ternary structure provides a complete picture of enzyme-substrate-cofactor interactions. The results reveal that different cofactors and substrate binding mode affect the allosteric effect between the two subunits of the enzyme. And the results from MD simulations confirmed that His221 plays a key role in the formation of dead-end complex in NADPH complex, and the absence of stable interaction between His221 and E1R in the NADH complex should be the main reason for its lack of substrate inhibition.

Viridiplantae-specific GLXI and GLXII isoforms co-evolved and detoxify glucosone in planta.

Reactive carbonyl species (RCS) such as methylglyoxal (MGO) and glyoxal (GO) are highly reactive, unwanted side-products of cellular metabolism maintained at harmless intracellular levels by specific scavenging mechanisms.MGO and GO are metabolized through the glyoxalase (GLX) system, which consists of two enzymes acting in sequence, GLXI and GLXII. While plant genomes encode a number of different GLX isoforms, their specific functions and how they arose during evolution are unclear. Here, we used Arabidopsis (Arabidopsis thaliana) as a model species to investigate the evolutionary history of GLXI and GLXII in plants and whether the GLX system can protect plant cells from the toxicity of RCS other than MGO and GO. We show that plants possess two GLX systems of different evolutionary origins and with distinct structural and functional properties. The first system is shared by all eukaryotes, scavenges MGO and GO, especially during seedling establishment, and features Zn2+-type GLXI proteins with a metal cofactor preference that were present in the last eukaryotic common ancestor. GLXI and GLXII of the second system, featuring Ni2+-type GLXI, were acquired by the last common ancestor of Viridiplantae through horizontal gene transfer from proteobacteria and can together metabolize keto-D-glucose (KDG, glucosone), a glucose-derived RCS, to D-gluconate. When plants displaying loss-of-function of a Viridiplantae-specific GLXI were grown in KDG, D-gluconate levels were reduced to 10%-15% of those in the wild type, while KDG levels showed an increase of 48%-67%. In contrast to bacterial GLXI homologs, which are active as dimers, plant Ni2+-type GLXI proteins contain a domain duplication, are active as monomers, and have a modified second active site. The acquisition and neofunctionalization of a structurally, biochemically, and functionally distinct GLX system indicates that Viridiplantae are under strong selection to detoxify diverse RCS.

Identification and investigation of a novel NADP+-dependent secoisolariciresinol dehydrogenase from Isatis indigotica.

Cofactors are crucial for the biosynthesis of natural compounds, and cofactor engineering is a useful strategy for enzyme optimization due to its potential to enhance enzyme efficiency. Secoisolariciresinol dehydrogenase (SIRD) was reported to convert secoisolariciresinol into matairesinol in an NAD+-dependent reaction. Here, a SIRD designated as IiSIRD2 identified from Isatis indigotica was found to utilize NADP+ as the cofactor. To explore the structural basis for this unique cofactor preference, model-based structural analysis was carried out, and it was postulated that a variation at the GXGGXG glycine-rich motif of IiSIRD2 alters its cofactor preference. This study paves way for future investigations on SIRD cofactor specificity and cofactor engineering to improve SIRD's catalytic efficiency.

Biocatalytic Cascade of Sebacic Acid Production with In Situ Co-Factor Regeneration Enabled by Engineering of an Alcohol Dehydrogenase

Sebacic acid (1,10-decanedioic acid) is an important chemical intermediate. Traditional chemical oxidation methods for sebacic acid production do not conform with “green” manufacturing. With the rapid development of enzymatic technologies, a biocatalytic cascade method based on the Baeyer–Villiger monooxygenase was developed. The most attractive point of the method is the oleic acid that can be utilized as raw material, which is abundant in nature. However, this bio-catalysis process needs co-factor electron carriers, and the high cost of the co-factor limits its progress. In this piece of work, a co-factor in situ regeneration system between ADH from Micrococcus luteus WIUJH20 (MlADH) and BVMO is proposed. Since the co-factors of both enzymes are different, switching the co-factor preference of native MlADH from NAD+ to NADP+ is necessary. Switching research was carried out based on in silico simulation, and the sites of Tyr36, Asp 37, Ala38, and Val39 were selected for mutation investigation. The experimental results demonstrated that mutants of MlADH_D37G and MlADH_D37G/A38T/V39K would utilize NADP+ efficiently, and the mutant of MlADH_D37G/A38T/V39K demonstrated the highest sebacic acid yield with the combination of BVMO. The results indicated that the in situ co-factor generation system is successfully developed, which would improve the efficiency of the biocatalytic cascade for sebacic acid production and is helpful for simplifying product isolation, thus, reducing the cost of the enzymatic transformations process.

Rationally Engineering the Cofactor Specificity of LfSDR1 for Biocatalytic Synthesis of the Key Intermediate of Telotristat Ethyl

AbstractSwitching cofactor preference of oxidoreductases from NADPH to NADH by rational engineering, replacing the expensive cofactor NADP+ with the cheap cofactor NAD+, is a focus of attention in the industrial application of oxidoreductases. This study focuses on the reversal of cofactor preference for short‐chain dehydrogenases/reductases (SDRs). Combined with bioinformatics analyses and in silico analyses, a small and smart mutant library (Mu1‐Mu3) of LfSDR1 was rationally designed and constructed. Thus, the excellent NADH‐dependent recombinant LfSDR1‐V186A/G92V/E141L/G38D/T15A variant (Mu2) was obtained. Meanwhile, novel enzymatic processes for synthesis of the key intermediates [(R)‐2 and (S)‐4] of telotristat ethyl and crizotinib were successfully created, which mainly relied on Mu2 coupled with an FDH‐catalyzed cofactor regeneration system. A co‐expressed E. coli whole‐cell biocatalyst containing the genes of Mu2 and PpFDH was developed to reduce ketones 1 and 3. Finally, ketone 1 was almost completely converted into the product (R)‐2 with a space‐time yield of 115.7 g⋅L−1⋅d−1 and a 98.8 % ee value.

A phosphite-based screening platform for identification of enzymes favoring nonnatural cofactors

Enzymes with dedicated cofactor preference are essential for advanced biocatalysis and biomanufacturing, especially when employing nonnatural nicotinamide cofactors in redox reactions. However, directed evolution of an enzyme to switch its cofactor preference is often hindered by the lack of efficient and affordable method for screening as the cofactor per se or the substrate can be prohibitively expensive. Here, we developed a growth-based selection platform to identify nonnatural cofactor-dependent oxidoreductase mutants. The growth of bacteria depended on the nicotinamide cytosine dinucleotide (NCD) mediated conversion of non-metabolizable phosphite into phosphate. The strain BW14329 lacking the ability to oxidize phosphite was suitable as host, and NCD-dependent phosphite dehydrogenase (Pdh*) is essential to the selection platform. Previously confirmed NCD synthetase with NCD synthesis capacity and NCD-dependent malic enzyme were successfully identified by using the platform. The feasibility of this strategy was successfully demonstrated using derived NCD-active malic enzyme as well as for the directed evolution of NCD synthetase in Escherichia coli. A phosphite-based screening platform was built for identification of enzymes favoring nonnatural cofactor NCD. In the future, once Pdh variants favoring other biomimetic or nonnatural cofactors are available this selection platform may be readily redesigned to attain new enzyme variants with anticipated cofactor preference, providing opportunities to further expand the chemical space of redox cofactors in chemical biology and synthetic biology.

The catalytic mechanism of the mitochondrial methylenetetrahydrofolate dehydrogenase/cyclohydrolase (MTHFD2)

Methylenetetrahydrofolate dehydrogenase/cyclohydrolase (MTHFD2) is a new drug target that is expressed in cancer cells but not in normal adult cells, which provides an Achilles heel to selectively kill cancer cells. Despite the availability of crystal structures of MTHFD2 in the inhibitor- and cofactor-bound forms, key information is missing due to technical limitations, including (a) the location of absolutely required Mg2+ ion, and (b) the substrate-bound form of MTHFD2. Using computational modeling and simulations, we propose that two magnesium ions are present at the active site whereby (i) Arg233, Asp225, and two water molecules coordinate , while together with Arg233 stabilize the inorganic phosphate (Pi); (ii) Asp168 and three water molecules coordinate , and further stabilizes Pi by forming a hydrogen bond with two oxygens of Pi; (iii) Arg201 directly coordinates the Pi; and (iv) through three water-mediated interactions, Asp168 contributes to the positioning and stabilization of , and Pi. Our computational study at the empirical valence bond level allowed us also to elucidate the detailed reaction mechanisms. We found that the dehydrogenase activity features a proton-coupled electron transfer with charge redistribution connected to the reorganization of the surrounding water molecules which further facilitates the subsequent cyclohydrolase activity. The cyclohydrolase activity then drives the hydration of the imidazoline ring and the ring opening in a concerted way. Furthermore, we have uncovered that two key residues, Ser197/Arg233, are important factors in determining the cofactor (NADP+/NAD+) preference of the dehydrogenase activity. Our work sheds new light on the structural and kinetic framework of MTHFD2, which will be helpful to design small molecule inhibitors that can be used for cancer treatment.

Altering NAD(P)H cofactor specificity by structure‐guided modification of class II HMG‐CoA reductase

HMG‐CoA reductase (HMGR) is a central enzyme in metabolism and natural product biosynthesis, as it catalyzes the committed and rate‐determining step of the mevalonate pathway, which produces steroids and isoprenoids, the largest and most diverse class of natural products. Here, HMGR converts HMG‐CoA to mevalonate using two equivalents of a redox cofactor, either NADH or NADPH. Class I HMGRs are highly specific for NADPH, including human HMGR, the target of cholesterol‐lowering statin drugs. However, cofactor specificity of microbial class II HMGR varies widely, as some homologs exhibit a strong preference for either NADH or NADPH, while others homologs are less specific. We have recently identified a key structural feature termed the “cofactor helix” that may govern cofactor specificity in class II HMGR, but until now this helix’s role in conferring NAD(P)H preference has not been tested. In this work, we performed structure‐guided enzyme modification with the aim of altering the cofactor preferences of various class II HMGR homologs. Through kinetic and X‐ray crystallographic studies on these enzyme variants, we have found that the cofactor helix plays a critical but not exclusive role in controlling NAD(P)H preference, expanding our understanding of the structural determinants of cofactor specificity in class II HMGR.

Engineering an acetoacetyl-CoA reductase from Cupriavidus necator toward NADH preference under physiological conditions

The coupling of PHB generation with NADH reoxidation is required to generate PHB as a fermentation product. A fundamental trait to accomplish this feature is to express a functional NADH-preferring acetoacetyl-CoA reductase, engaged in PHB accumulation. One way to obtain such a reductase is by engineering the cofactor preference of the acetoacetyl-CoA reductase encoded by the phaB1 gene from Cupriavidus necator (AARCn1). Aiming to have a deeper understanding of the structural determinants of the cofactor preference in AARCn1, and to obtain an NADH-preferring acetoacetyl-CoA reductase derived from this protein, some engineered enzymes were expressed, purified and kinetically characterized, together with the parental AARCn1. One of these engineered enzymes, Chimera 5, experimentally showed a selectivity ratio ((kcat/KM)NADH/(kcat/KM)NADPH) ≈ 18, which is 160 times higher than the selectivity ratio experimentally observed in the parental AARCn1. A thermodynamic-kinetic approach was employed to estimate the cofactor preference and flux capacity of Chimera 5 under physiological conditions. According to this approach, Chimera 5 could prefer NADH over NADPH between 25 and 150 times. Being a derivative of AARCn1, Chimera 5 should be readily functional in Escherichia coli and C. necator. Moreover, with the expected expression level, its activity should be enough to sustain PHB accumulation fluxes similar to the fluxes previously observed in these biotechnologically relevant cell factories.

Cryo-EM structure of transcription termination factor Rho from Mycobacterium tuberculosis reveals bicyclomycin resistance mechanism

The bacterial Rho factor is a ring-shaped motor triggering genome-wide transcription termination and R-loop dissociation. Rho is essential in many species, including in Mycobacterium tuberculosis where rho gene inactivation leads to rapid death. Yet, the M. tuberculosis Rho [MtbRho] factor displays poor NTPase and helicase activities, and resistance to the natural Rho inhibitor bicyclomycin [BCM] that remain unexplained. To address these issues, we solved the cryo-EM structure of MtbRho at 3.3 Å resolution. The MtbRho hexamer is poised into a pre-catalytic, open-ring state wherein specific contacts stabilize ATP in intersubunit ATPase pockets, thereby explaining the cofactor preference of MtbRho. We reveal a leucine-to-methionine substitution that creates a steric bulk in BCM binding cavities near the positions of ATP γ-phosphates, and confers resistance to BCM at the expense of motor efficiency. Our work contributes to explain the unusual features of MtbRho and provides a framework for future antibiotic development.

Debottlenecking 4-hydroxybenzoate hydroxylation in Pseudomonas putida KT2440 improves muconate productivity from p-coumarate

The transformation of 4-hydroxybenzoate (4-HBA) to protocatechuate (PCA) is catalyzed by flavoprotein oxygenases known as para-hydroxybenzoate-3-hydroxylases (PHBHs). In Pseudomonas putida KT2440 (P. putida) strains engineered to convert lignin-related aromatic compounds to muconic acid (MA), PHBH activity is rate-limiting, as indicated by the accumulation of 4-HBA, which ultimately limits MA productivity. Here, we hypothesized that replacement of PobA, the native P. putida PHBH, with PraI, a PHBH from Paenibacillus sp. JJ-1b with a broader nicotinamide cofactor preference, could alleviate this bottleneck. Biochemical assays confirmed the strict preference of NADPH for PobA, while PraI can utilize either NADH or NADPH. Kinetic assays demonstrated that both PobA and PraI can utilize NADPH with comparable catalytic efficiency and that PraI also efficiently utilizes NADH at roughly half the catalytic efficiency. The X-ray crystal structure of PraI was solved and revealed absolute conservation of the active site architecture to other PHBH structures despite their differing cofactor preferences. To understand the effect in vivo, we compared three P. putida strains engineered to produce MA from p-coumarate (pCA), showing that expression of praI leads to lower 4-HBA accumulation and decreased NADP+/NADPH ratios relative to strains harboring pobA, indicative of a relieved 4-HBA bottleneck due to increased NADPH availability. In bioreactor cultivations, a strain exclusively expressing praI achieved a titer of 40 g/L MA at 100% molar yield and a productivity of 0.5 g/L/h. Overall, this study demonstrates the benefit of sampling readily available natural enzyme diversity for debottlenecking metabolic flux in an engineered strain for microbial conversion of lignin-derived compounds to value-added products.

Rational engineering of cofactor specificity of glutamate dehydrogenase for poly-γ-glutamic acid synthesis in Bacillus licheniformis

Poly-γ-glutamic acid (γ-PGA) is a multifunctional biopolymer mainly produced by Bacillus. The cofactor specificity of enzymes plays a critical role in regulating metabolic process and metabolite production. Here, we present a novel approach for switching cofactor specificity of glutamate dehydrogenase RocG from nicotinamide adenine dinucleotide phosphate (NADPH) to nicotinamide adenine dinucleotide (NADH) to improve γ-PGA production. Firstly, 3D structural modeling and molecular docking were performed to predict the binding modes of NADH and NADPH. Several site-specific mutants based on the conventional and Random Accelerated Molecular Dynamics simulations were obtained to alter cofactor specificity. Then, the effects of RocG variants overexpressions on γ-PGA production were evaluated. Compared to the wild-type, the mutant RocGD276E showed highest increase in γ-PGA yield, increased by 40.50%. Meanwhile, yields of main by-products acetoin and 2,3-butandieol were decreased by 21.70% and 16.53%, respectively. Finally, the results of enzymatic properties confirmed that glutamate dehydrogenase mutant RocGD276E exhibited the higher affinity for NADH, caused a shift in coenzyme preference from NADPH to NADH, with a catalytic efficiency comparable with NADPH-dependent RocG. Taken together, this research demonstrated that switching the cofactor preference of glutamate dehydrogenase via rational design was an effective strategy for high-level production of γ-PGA in Bacillus licheniformis.

A Redox-Neutral, Two-Enzyme Cascade for the Production of Malate and Gluconate from Pyruvate and Glucose

A triple mutant of NADP(H)-dependent malate dehydrogenase from thermotolerant Thermococcus kodakarensis has an altered cofactor preference for NAD+, as well as improved malate production compared to wildtype malate dehydrogenase. By combining mutant malate dehydrogenase with glucose dehydrogenase from Sulfolobus solfataricus and NAD+/NADH in a closed reaction environment, gluconate and malate could be produced from pyruvate and glucose. After 3 h, the yield of malate was 15.96 mM. These data demonstrate the feasibility of a closed system capable of cofactor regeneration in the production of platform chemicals.

Engineering Promiscuous Alcohol Dehydrogenase Activity of a Reductive Aminase AspRedAm for Selective Reduction of Biobased Furans.

Catalytic promiscuity is a promising starting point for improving the existing enzymes and even creating novel enzymes. In this work, site-directed mutagenesis was performed to improve promiscuous alcohol dehydrogenase activity of reductive aminase from Aspergillus oryzae (AspRedAm). AspRedAm showed the cofactor preference toward NADPH in reductive aminations, while it favored NADH in the reduction reactions. Some key amino acid residues such as N93, I118, M119, and D169 were identified for mutagenesis by molecular docking. Variant N93A showed the optimal pH and temperature of 8 and 30°C, respectively, in the reduction of 5-hydroxymethylfurfural (HMF). The thermostability was enhanced upon mutation of N93 to alanine. The catalytic efficiency of variant N93A (k cat/K m, 23.6 mM−1 s−1) was approximately 2-fold higher compared to that of the wild-type (WT) enzyme (13.1 mM−1 s−1). The improved catalytic efficiency of this variant may be attributed to the reduced steric hindrance that stems from the smaller side chain of alanine in the substrate-binding pocket. Both the WT enzyme and variant N93A had broad substrate specificity. Escherichia coli (E. coli) cells harboring plain vector enabled selective reduction of biobased furans to target alcohols, with the conversions of 35–95% and the selectivities of >93%. The introduction of variant N93A to E. coli resulted in improved substrate conversions (>98%) and selectivities (>99%).

Cofactor Specificity of Glucose-6-Phosphate Dehydrogenase Isozymes in Pseudomonas putida Reveals a General Principle Underlying Glycolytic Strategies in Bacteria.

ABSTRACTGlucose-6-phosphate dehydrogenase (G6PDH) is widely distributed in nature and catalyzes the first committing step in the oxidative branch of the pentose phosphate (PP) pathway, feeding either the reductive PP or the Entner-Doudoroff pathway. Besides its role in central carbon metabolism, this dehydrogenase provides reduced cofactors, thereby affecting redox balance. Although G6PDH is typically considered to display specificity toward NADP+, some variants accept NAD+ similarly or even preferentially. Furthermore, the number of G6PDH isozymes encoded in bacterial genomes varies from none to more than four orthologues. On this background, we systematically analyzed the interplay of the three G6PDH isoforms of the soil bacterium Pseudomonas putida KT2440 from genomic, genetic, and biochemical perspectives. P. putida represents an ideal model to tackle this endeavor, as its genome harbors gene orthologues for most dehydrogenases in central carbon metabolism. We show that the three G6PDHs of strain KT2440 have different cofactor specificities and that the isoforms encoded by zwfA and zwfB carry most of the activity, acting as metabolic “gatekeepers” for carbon sources that enter at different nodes of the biochemical network. Moreover, we demonstrate how multiplication of G6PDH isoforms is a widespread strategy in bacteria, correlating with the presence of an incomplete Embden-Meyerhof-Parnas pathway. The abundance of G6PDH isoforms in these species goes hand in hand with low NADP+ affinity, at least in one isozyme. We propose that gene duplication and relaxation in cofactor specificity is an evolutionary strategy toward balancing the relative production of NADPH and NADH.IMPORTANCE Protein families have likely arisen during evolution by gene duplication and divergence followed by neofunctionalization. While this phenomenon is well documented for catabolic activities (typical of environmental bacteria that colonize highly polluted niches), the coexistence of multiple isozymes in central carbon catabolism remains relatively unexplored. We have adopted the metabolically versatile soil bacterium Pseudomonas putida KT2440 as a model to interrogate the physiological and evolutionary significance of coexisting glucose-6-phosphate dehydrogenase (G6PDH) isozymes. Our results show that each of the three G6PDHs in this bacterium display distinct biochemical properties, especially at the level of cofactor preference, impacting bacterial physiology in a carbon source-dependent fashion. Furthermore, the presence of multiple G6PDHs differing in NAD+ or NADP+ specificity in bacterial species strongly correlates with their predominant metabolic lifestyle. Our findings support the notion that multiplication of genes encoding cofactor-dependent dehydrogenases is a general evolutionary strategy toward achieving redox balance according to the growth conditions.

Alteration of cofactor specificity of the acrylyl-CoA reductase from Escherichia coli.

Alteration of the cofactor specificity of acrylyl-CoA reductase (AcuI) catalyzing the NAD(P)H-dependent reduction of acrylyl-CoA to propionyl-CoA is often desirable for designing of artificial metabolic pathways of various appointments. Several variants of AcuIs from Escherichia coli K-12 with multiple amino acid substitutions to alter the cofactor preference were obtained by site directed mutagenesis and the modified enzymes as His6-tagged proteins were characterized. The simultaneous substitutions of arginine-180, arginine-198 and serine-178 residues by alanine in the enzyme pocket sequence as well as other amino acid changes decreased both NADPH- and NADH-dependent activities in comparison to the wild-type enzyme. The replacement of serine-156 by glutamic acid decreased NADPH-dependent activity at least 7000-fold but NADH-dependent activity only by threefold. The replacement of serine-156 by aspartic acid decreased NADPH-dependent activity 70-fold with fair preservation of activity and specificity to NADH. These results demonstrated a relevance of Asp156 in the interaction of AcuI from E. coli K-12 with NADH as a coenzyme. These findings may provide reference information for shifting coenzyme specificity of acrylyl-CoA reductases.

Structural Insights into Malic Enzyme Variants Favoring an Unnatural Redox Cofactor.

The use of nicotinamide cytosine dinucleotide (NCD), a biocompatible nicotinamide adenosine dinucleotide (NAD) analogue, is of great scientific and biotechnological interest. Several redox enzymes have been devised to favor NCD, and have been successfully applied in creating NCD-dependent redox systems. However, molecular interactions between cofactor and protein have still to be disclosed in order to guide further engineering efforts. Here we report the structural analysis of an NCD-favoring malic enzyme (ME) variant derived from Escherichia coli. The X-ray crystal structure data revealed that the residues located at position 346 and 401 in ME acted as the "gatekeepers" of the adenine moiety binding cavity. When Arg346 was substituted with either acidic or aromatic residues, the corresponding mutants showed substantially reduced NCD preference. Inspired by these observations, we generated Lactobacillus helveticus derived d-lactate dehydrogenase variants at Ile177, the counterpart to Arg346 in ME, and found a similar trend in terms of cofactor preference changes. As many NAD-dependent oxidoreductases share key structural features, our results provide guidance for protein engineering to obtain more NCD-favoring variants.

Crystal structure of the 6-phosphogluconate dehydrogenase from Gluconobacter oxydans reveals tetrameric 6PGDHs as the crucial intermediate in the evolution of structure and cofactor preference in the 6PGDH family

Background: The enzyme 6-phosphogluconate dehydrogenase (6PGDH) is the central enzyme of the oxidative pentose phosphate pathway. Members of the 6PGDH family belong to different classes: either homodimeric enzymes assembled from long-chain subunits or homotetrameric ones assembled from short-chain subunits. Dimeric 6PGDHs bear an internal duplication absent in tetrameric 6PGDHs and distant homologues of the β-hydroxyacid dehydrogenase (βHADH) superfamily. Methods: We use X-ray crystallography to determine the structure of the apo form of the 6PGDH from Gluconobacter oxydans (Go6PGDH). We carried out a structural and phylogenetic analysis of short and long-chain 6PGDHs. We put forward an evolutionary hypothesis explaining the differences seen in oligomeric state vs. dinucleotide preference of the 6PGDH family. We determined the cofactor preference of Go6PGDH at different 6-phosphogluconate concentrations, characterizing the wild-type enzyme and three-point mutants of residues in the cofactor binding site of Go6PGDH. Results: The structural comparison suggests that the 6PG binding site initially evolved by exchanging C-terminal α-helices between subunits. An internal duplication event changed the quaternary structure of the enzyme from a tetrameric to a dimeric arrangement. The phylogenetic analysis suggests that 6PGDHs have spread from Bacteria to Archaea and Eukarya on multiple occasions by lateral gene transfer. Sequence motifs consistent with NAD+- and NADP+-specificity are found in the β2-α2 loop of dimeric and tetrameric 6PGDHs. Site-directed mutagenesis of Go6PGDH inspired by this analysis fully reverses dinucleotide preference. One of the mutants we engineered has the highest efficiency and specificity for NAD+ so far described for a 6PGDH. Conclusions: The family 6PGDH comprises dimeric and tetrameric members whose active sites are conformed by a C-terminal α-helix contributed from adjacent subunits. Dimeric 6PGDHs have evolved from the duplication-fusion of the tetrameric C-terminal domain before independent transitions of cofactor specificity. Changes in the conserved β2-α2 loop are crucial to modulate the cofactor specificity in Go6PGDH.

Dehydroquinate dehydratase/shikimate dehydrogenases involved in gallate biosynthesis of the aluminum-tolerant tree species Eucalyptus camaldulensis

Main conclusionEucalyptus camaldulensis EcDQD/SDH2 and 3 combine gallate formation, dehydroquinate dehydratase, and shikimate dehydrogenase activities. They are candidates for providing the essential gallate for the biosynthesis of the aluminum-detoxifying metabolite oenothein B.The tree species Eucalyptus camaldulensis shows exceptionally high tolerance against aluminum, a widespread toxic metal in acidic soils. In the roots of E. camaldulensis, aluminum is detoxified via the complexation with oenothein B, a hydrolyzable tannin. In our approach to elucidate the biosynthesis of oenothein B, we here report on the identification of E. camaldulensis enzymes that catalyze the formation of gallate, which is the phenolic constituent of hydrolyzable tannins. By systematical screening of E. camaldulensis dehydroquinate dehydratase/shikimate dehydrogenases (EcDQD/SDHs), we found two enzymes, EcDQD/SDH2 and 3, catalyzing the NADP+-dependent oxidation of 3-dehydroshikimate to produce gallate. Based on extensive in vitro assays using recombinant EcDQD/SDH2 and 3 enzymes, we present for the first time a detailed characterization of the enzymatic gallate formation activity, including the cofactor preferences, pH optima, and kinetic constants. Sequence analyses and structure modeling suggest the gallate formation activity of EcDQD/SDHs is based on the reorientation of 3-dehydroshikimate in the catalytic center, which facilitates the proton abstraction from the C5 position. Additionally, EcDQD/SDH2 and 3 maintain DQD and SDH activities, resulting in a 3-dehydroshikimate supply for gallate formation. In E. camaldulensis, EcDQD/SDH2 and 3 are co-expressed with UGT84A25a/b and UGT84A26a/b involved in hydrolyzable tannin biosynthesis. We further identified EcDQD/SDH1 as a “classical” bifunctional plant shikimate pathway enzyme and EcDQD/SDH4a/b as functional quinate dehydrogenases of the NAD+/NADH-dependent clade. Our data indicate that in E. camaldulensis the enzymes EcDQD/SDH2 and 3 provide the essential gallate for the biosynthesis of the aluminum-detoxifying metabolite oenothein B.

Structure-Guided Design of Formate Dehydrogenase for Regeneration of a Non-Natural Redox Cofactor.

Formate dehydrogenase (FDH) has been widely used for the regeneration of the reduced nicotinamide adenine dinucleotide (NADH). To utilize nicotinamide cytosine dinucleotide (NCD) as a non-natural redox cofactor, it remains challenging as NCDH, the reduced form of NCD, has to be efficiently regenerated. Here we demonstrate successful engineering of FDH for NCDH regeneration. Guided by the structural information of FDH from Pseudomonas sp. 101 (pseFDH) and the NAD-pseFDH complex, semi-rational strategies were applied to design mutant libraries and screen for NCD-linked activity. The most active mutant reached a cofactor preference switch from NAD to NCD by 3700-fold. Homology modeling analysis showed that these mutants had reduced cofactor binding pockets and dedicated hydrophobic interactions for NCD. Efficient regeneration of NCDH was implemented by powering an NCD-dependent D-lactate dehydrogenase for stoichiometric and stereospecific reduction of pyruvate to D-lactate at the expense of formate.