Change in Cofactor Specificity of Oxidoreductases by Adaptive Evolution of an Escherichia coli NADPH-Auxotrophic Strain

In the cell, NAD(H) and NADP(H) cofactors have different functions. The former mainly accepts electrons from catabolic reactions and carries them to respiration, while the latter provides reducing power for anabolism.

High-Yield ‘One-Pot’ Biosynthesis of Raspberry Ketone, a High-Value Fine Chemical

Cell-free extract and purified enzyme-based systems provide an attractive solution to study biosynthetic strategies towards a range of chemicals. 4-(4-hydroxyphenyl)-butan-2-one, also known as raspberry ketone, is the major fragrance component of raspberry fruit and is used as a natural additive in the food and sports industry. Current industrial processing of the natural form of raspberry ketone involves chemical extraction with a yield of ~1-4 mg kg-1 of fruit. Due to toxicity, microbial production provides only low yields of up to 5-100 mg L-1. Herein, we report an efficient cell-free strategy to probe a synthetic enzyme pathway that converts either L-tyrosine or the precursor, 4-(4-hydroxyphenyl)-buten-2-one (HBA), into raspberry ketone at up to 100% conversion. As part of this strategy, it is essential to recycle inexpensive cofactors. Specifically, the final enzyme step in the pathway is catalysed by raspberry ketone/zingerone synthase (RZS1), an NADPH-dependent double bond reductase. To relax cofactor specificity towards NADH, the preferred cofactor for cell-free biosynthesis, we identify a variant (G191D) with strong activity with NADH. We implement the RZS1 G191D variant within a ‘one-pot’ cell-free reaction to produce raspberry ketone at high-yield (61 mg L-1), which provides an alternative route to traditional microbial production. In conclusion, our cell-free strategy complements the growing interest in engineering synthetic enzyme cascades towards industrially relevant value-added chemicals.

Engineering Cofactor Specificity of a Thermostable Phosphite Dehydrogenase for a Highly Efficient and Robust NADPH Regeneration System

Nicotinamide adenine dinucleotide phosphate (NADP)-dependent dehydrogenases catalyze a range of chemical reactions useful for practical applications. However, their dependence on the costly cofactor, NAD(P)H remains a challenge which must be addressed. Here, we engineered a thermotolerant phosphite dehydrogenase from Ralstonia sp. 4506 (RsPtxD) by relaxing the cofactor specificity for a highly efficient and robust NADPH regeneration system. The five amino acid residues, Cys174–Pro178, located at the C-terminus of β7-strand region in the Rossmann-fold domain of RsPtxD, were changed by site-directed mutagenesis, resulting in four mutants with a significantly increased preference for NADP. The catalytic efficiency of mutant RsPtxDHARRA for NADP (Kcat/KM)NADP was 44.1 μM–1 min–1, which was the highest among the previously reported phosphite dehydrogenases. Moreover, the RsPtxDHARRA mutant exhibited high thermostability at 45°C for up to 6 h and high tolerance to organic solvents, when bound with NADP. We also demonstrated the applicability of RsPtxDHARRA as an NADPH regeneration system in the coupled reaction of chiral conversion of 3-dehydroshikimate to shikimic acid by the thermophilic shikimate dehydrogenase of Thermus thermophilus HB8 at 45°C, which could not be supported by the parent RsPtxD enzyme. Therefore, the RsPtxDHARRA mutant might be a promising alternative NADPH regeneration system for practical applications.

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

ABSTRACT Glucose-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.

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.

Emergence of NADP+-reducing enzymes in Escherichia coli central metabolism via adaptive evolution

The nicotinamide cofactor specificity of enzymes plays a key role in regulating metabolic processes and attaining cellular homeostasis. Multiple studies have used enzyme engineering tools or a directed evolution approach to switch the cofactor preference of specific oxidoreductases. However, whole-cell adaptation towards the emergence of novel cofactor regeneration routes was not previously explored. To address this challenge, we used an Escherichia coli NADPH-auxotroph strain. We continuously cultivated this strain under selective conditions. After 500-1100 generations of adaptive evolution using different carbon sources, we isolated several strains capable of growing without an external NADPH source. Most isolated strains were found to harbor a mutated NAD-dependent malic enzyme (MaeA). A single mutation in MaeA was found to switch cofactor specificity while lowering enzyme activity. Most mutated MaeA variants also harbored a second mutation that restored the catalytic efficiency of the enzyme. Remarkably, the best MaeA variants identified this way displayed overall superior kinetics relative to the wildtype variant with NAD+. In other evolved strains, the dihydrolipoamide dehydrogenase (Lpd) was mutated to accept NADP+ thus enabling the pyruvate dehydrogenase and 2-ketoglutarate dehydrogenase complexes to regenerate NADPH. Interestingly, no other central metabolism oxidoreductase seems to evolve towards reducing NADP+, which we attribute to several biochemical constraints such as unfavorable thermodynamics. This study demonstrates the potential and biochemical limits of evolving oxidoreductases within the cellular context towards changing cofactor specificity, further showing that long-term adaptive evolution can optimize enzyme activity beyond what is achievable via rational design or directed evolution using small libraries.ImportanceIn the cell, NAD(H) and NADP(H) cofactors have different functions. The former mainly accepts electrons from catabolic reactions and carries them to respiration, while the latter provides reducing power for anabolism. Correspondingly, the ratio of the reduced to the oxidized form differs for NAD (low) and NADP (high), reflecting their distinct roles. We challenged the flexibility of E. coli’s central metabolism in multiple adaptive evolution experiments using an NADPH-auxotroph strain. We found several mutations in two enzymes, changing the cofactor preference of malic enzyme and dihydrolipoamide dehydrogenase. Upon deletion of their corresponding genes we performed additional evolution experiments which did not lead to the emergence of any additional mutants. We attribute this restricted number of mutational targets to intrinsic thermodynamic barriers: The high ratio of NADPH to NADP+ limits metabolic redox reactions which can regenerate NADPH, mainly by mass action constraints.

Cofactor specificity of glucose-6-phosphate dehydrogenase isozymes in Pseudomonas putida reveals a general principle underlying glycolytic strategies in bacteria

Glucose-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 also provides reduced cofactors, thereby affecting redox balance. Although G6PDH is typically considered to display specificity towards nicotinamide adenine dinucleotide phosphate (NADP+), some variants accept nicotinamide 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 a genomic, genetic and biochemical perspective. P. putida represents an ideal model to tackle this endeavor, as its genome encodes numerous 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. Multiplication of G6PDH isoforms in these species goes hand-in-hand with low NADP+ affinity at least in one G6PDH isozyme. We propose that gene duplication and relaxation in cofactor specificity is an evolutionary strategy towards balancing the relative production of NADPH and NADH.

Exploration of the cofactor specificity of wild-type phosphite dehydrogenase and its mutant using molecular dynamics simulations

Phosphite dehydrogenase (Pdh) catalyzes the NAD-dependent oxidation of phosphite to phosphate with the formation of NADH.