Halomonas bluephagenesis has been developed as a platform strain for the next generation industrial biotechnology (NGIB) with advantages of resistances to microbial contamination and high cell density growth (HCD), especially for production of polyhydroxyalkanoates (PHA) including poly(3-hydroxybutyrate) (PHB), poly(3-hydroxybutyrate-co-4-hydroxybutyrate) (P34HB) and poly(3-hydroxybutyrate-co-3-hydroxyvalerate) (PHBV). However, little is known about the mechanism behind PHA accumulation under oxygen limitation. This study for the first time found that H. bluephagenesis utilizes NADH instead of NADPH as a cofactor for PHB production, thus revealing the rare situation of enhanced PHA accumulation under oxygen limitation. To increase NADH/NAD ratio for enhanced PHA accumulation under oxygen limitation, an electron transport pathway containing electron transfer flavoprotein subunits α and β encoded by etf operon was blocked to increase NADH supply, leading to 90% PHB accumulation in the cell dry weight (CDW) of H. bluephagenesis compared with 84% by the wild type. Acetic acid, a cost-effective carbon source, was used together with glucose to balance the redox state and reduce inhibition on pyruvate metabolism, resulting in 22% more CDW and 94% PHB accumulation. The cellular redox state changes induced by the addition of acetic acid increased 3HV ratio in its copolymer PHBV from 4% to 8%, 4HB in its copolymer P34HB from 8% to 12%, respectively, by engineered H. bluephagenesis. The strategy of systematically modulation on the redox potential of H. bluephagenesis led to enhanced PHA accumulation and controllable monomer ratios in PHA copolymers under oxygen limitation, reducing energy consumption and scale-up complexity.

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

In Escherichia coli, the oxidative branch of the pentose phosphate pathway (oxPPP) is one of the major sources of NADPH when glucose is the sole carbon nutrient. However, unbalanced NADPH production causes growth impairment as observed in a strain lacking phosphoglucoisomerase (Dpgi). In this work, we studied the metabolic response of this bacterium to the replacement of its glucose-6-phosphate dehydrogenase (G6PDH) by an NADH-producing variant. The homologous enzyme from Leuconostoc mesenteroides was studied by molecular dynamics and site-directed mutagenesis to obtain the NAD-preferring LmG6PDH R46E,Q47E . Through homologous recombination, the zwf loci (encoding G6PDH) in the chromosomes of WT and Dpgi E. coli strains were replaced by DNA encoding LmG6PDH R46E,Q47E . Contrary to some predictions performed with flux balance analysis, the replacements caused a substantial effect on the growth rates, increasing 59 % in the Dpgi strain, while falling 44 % in the WT. Quantitative PCR (qPCR) analysis of the zwf locus showed that the expression level of the mutant enzyme was similar to the native enzyme and the expression of genes encoding key enzymes of the central pathways also showed moderate changes among the studied strains. The phenotypic and qPCR data were integrated into in silico modelling, showing an operative G6PDH flux contributing to the NADH pool. Our results indicated that, in vivo, the generation of NADH by G6PDH is beneficial or disadvantageous for growth depending on the operation of the upper Embden-Meyerhof pathway. Interestingly, a genomic database search suggested that in bacteria lacking phosphofructokinase, the G6PDHs tend to have similar preferences for NAD and NADP. The importance of the generation of NADPH in a pathway such as the oxPPP is discussed.

Environmental fluctuations in the availability of nutrients lead to intricate metabolic strategies. Candidatus Accumulibacter phosphatis, a polyphosphate accumulating organism (PAO) responsible for enhanced biological phosphorus removal (EBPR) from wastewater treatment systems, is prevalent in aerobic/anaerobic environments. While the overall metabolic traits of these bacteria are well described, the non-availability of isolates has led to controversial conclusions on the metabolic pathways used. Here, we experimentally determined the redox cofactor preference of different oxidoreductases in the central carbon metabolism of a highly enriched Ca. A. phosphatis culture. Remarkably, we observed that the acetoacetyl-CoA reductase engaged in polyhydroxyalkanoates (PHA) synthesis is NADH-preferring instead of the generally assumed NADPH dependency. This allows re-thinking the ecological role of PHA accumulation as a fermentation product under anaerobic conditions and not just a stress response. Based on previously published meta-omics data and the results of enzymatic assays, a reduced central carbon metabolic network was constructed and used for simulating different metabolic operating modes. In particular, scenarios with different acetate-to-glycogen consumption ratios were simulated, which demonstrated optima using different combinations of glycolysis, glyoxylate shunt or branches of the TCA cycle. Thus, optimal metabolic flux strategies will depend on the environment (acetate uptake) and on intracellular storage compounds availability (polyphosphate/glycogen). This NADH-related metabolic flexibility is enabled by the NADH-driven PHA synthesis. It allows for maintaining metabolic activity under varying environmental substrate conditions, with high carbon conservation and lower energetic costs compared to NADPH dependent PHA synthesis. Such (flexible) metabolic redox coupling can explain PAOs' competitiveness under oxygen-fluctuating environments. IMPORTANCE Here we demonstrate how microbial storage metabolism can adjust to a wide range of environmental conditions. Such flexibility generates a selective advantage under fluctuating environmental conditions. It can also explain the different observations reported in PAO literature, including the capacity of Ca. Accumulibacter phosphatis to act like glycogen accumulating organisms (GAO). These observations stem from slightly different experimental conditions and controversy only arises when one assumes metabolism can only operate in one single mode. Furthermore, we also show how the study of metabolic strategies is possible when combining -omics data with functional cofactor assays and modeling. Genomic information can only provide the potential of a microorganism. The environmental context and other complementary approaches are still needed to study and predict the functional expression of such metabolic potential.

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.