Utilizing anaerobic metabolisms for the production of biotechnologically relevant products presents potential advantages, such as increased yields and reduced energy dissipation. However, lower energy dissipation may indicate that certain reactions are operating closer to their thermodynamic equilibrium. While stoichiometric analyses and genetic modifications are frequently employed in metabolic engineering, the use of thermodynamic tools to evaluate the feasibility of planned interventions is less documented. In this study, we propose a novel metabolic engineering strategy to achieve an efficient anaerobic production of poly-(R)-3-hydroxybutyrate (PHB) in the model organism Escherichia coli. Our approach involves re-routing of two-thirds of the glycolytic flux through non-oxidative glycolysis and coupling PHB synthesis with NADH re-oxidation. We complemented our stoichiometric analysis with various thermodynamic approaches to assess the feasibility and the bottlenecks in the proposed engineered pathway. According to our calculations, the main thermodynamic bottleneck are the reactions catalyzed by the acetoacetyl-CoA β-ketothiolase (EC 2.3.1.9) and the acetoacetyl-CoA reductase (EC 1.1.1.36). Furthermore, we calculated thermodynamically consistent sets of kinetic parameters to determine the enzyme amounts required for sustaining the conversion fluxes. In the case of the engineered conversion route, the protein pool necessary to sustain the desired fluxes could account for 20% of the whole cell dry weight.

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 (AAR^Cn1). Aiming to have a deeper understanding of the structural determinants of the cofactor preference in AAR^Cn1, 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 AAR^Cn1. One of these engineered enzymes, Chimera 5, experimentally showed a selectivity ratio ((k_cat/K_M)^NADH/(k_cat/K_M)^NADPH) ≈ 18, which is 160 times higher than the selectivity ratio experimentally observed in the parental AAR^Cn1. 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 AAR^Cn1, 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.

Oxygen supply implies higher production cost and reduction of maximum theoretical yields. Thus, generation of fermentation products is more cost-effective. Aiming to find a key piece for the production of (poly)-3-hydroxybutyrate (PHB) as a fermentation product, here we characterize an acetoacetyl-CoA reductase, isolated from a Candidatus Accumulibacter phosphatis-enriched mixed culture, showing a (k_cat^NADH/K_M^NADH)/(k_cat^NADPH/K_M^NADPH)>500. Further kinetic analyses indicate that, at physiological concentrations, this enzyme clearly prefers NADH, presenting the strongest NADH preference so far observed among the acetoacetyl-CoA reductases. Structural and kinetic analyses indicate that residues between E37 and P41 have an important role for the observed NADH preference. Moreover, an operon was assembled combining the phaCA genes from Cupriavidus necator and the gene encoding for this NADH-preferring acetoacetyl-CoA reductase. Escherichia coli cells expressing that assembled operon showed continuous accumulation of PHB under oxygen limiting conditions and PHB titer increased when decreasing the specific oxygen consumption rate. Taken together, these results show that it is possible to generate PHB as a fermentation product in E. coli, opening opportunities for further protein/metabolic engineering strategies envisioning a more efficient anaerobic production of PHB.

Biosynthesis of poly-3-hydroxybutyrate (PHB) as a fermentation product enables the coupling of growth and product generation. Moreover, the reduction of oxygen supply should reduce operative cost and increase product yield. Generation of PHB as a fermentation product depends on the in vivo activity of an NADH-preferring acetoacetyl-CoA reductase. Proof of this concept requires (i) quantification of the cofactor preference, in physiologically relevant conditions, of a putative NADH-preferring acetoacetyl-CoA reductase and (ii) verification of PHB accumulation using an NADH-preferring acetoacetyl-CoA reductase in a species naturally incapable of doing so, for example, Escherichia coli. This dataset contains kinetic data obtained by spectrophotometry and data from a continuous culture of an engineered E. coli strain accumulating PHB under oxygen-limiting conditions. In this dataset it is possible to find (1) enzyme stability assays; (2) initial rates and progress curves from reactions catalyzed by two acetoacetyl-CoA reductases; (3) estimations of the relative use of NADH and NADPH by two acetoacetyl-CoA reductases; (4) estimations of the flux capacity of the reaction catalyzed by an acetoacetyl-CoA reductase; (5) biomass composition of an engineered E. coli strain transformed with a plasmid; (6) calculation of reconciled specific rates of this engineered strain growing on sucrose as the sole carbon source under oxygen limitation and (7) metabolic fluxes distributions during the continuous growth of this engineered strain. Because a relatively small number of acetoacetyl-CoA reductases have been kinetically characterized, data and scripts here provided could be useful for further kinetic characterizations. Moreover, the procedure described to estimate biomass composition could be interesting to estimate plasmid and protein burden in other strains. Application of data reconciliation to fermentations should help to obtain specific rates consistent with the principle of mass and electron conservation. All the required data and scripts to perform these analyses are deposited in a Mendeley Data repository. This article was co-submitted with the manuscript entitled “An NADH preferring acetoacetyl-CoA reductase is engaged in poly-3-hydroxybutyrate accumulation in Escherichiasia. coli”.

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 nonavailability of isolates has led to controversial conclusions on the metabolic pathways used. In this study, we experimentally determined the redox cofactor preferences of different oxidoreductases in the central carbon metabolism of a highly enriched "Ca Accumulibacter phosphatis" culture. Remarkably, we observed that the acetoacetyl coenzyme A reductase engaged in polyhydroxyalkanoate (PHA) synthesis is NADH preferring instead of showing the generally assumed NADPH dependency. This allows rethinking of the ecological role of PHA accumulation as a fermentation product under anaerobic conditions and not just a stress response. Based on previously published metaomics 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 tricarboxylic acid (TCA) cycle. Thus, optimal metabolic flux strategies will depend on the environment (acetate uptake) and on intracellular storage compound availability (polyphosphate/glycogen). This NADH-related metabolic flexibility is enabled by the NADH-driven PHA synthesis. It allows for maintaining metabolic activity under various environmental substrate conditions, with high carbon conservation and lower energetic costs than for NADPH-dependent PHA synthesis. Such (flexible) metabolic redox coupling can explain the competitiveness of PAOs 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 (GAOs). These observations stem from slightly different experimental conditions, and controversy arises only when one assumes that metabolism can operate only in a 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.

Biotechnological industry strives to develop anaerobic bioprocesses fueled by abundant and cheap carbon sources, like sucrose. However, oxygen-limiting conditions often lead to by-product formation and reduced ATP yields. While by-product formation is typically decreased by gene deletion, the breakdown of oligosaccharides with inorganic phosphate instead of water could increment the ATP yield. To observe the effect of oxygen limitation during sucrose consumption, a non-fermentative Escherichia coli K-12 strain was transformed with genes enabling sucrose assimilation. It was observed that the combined deletion of the genes adhE, adhP, mhpF, ldhA, and pta abolished the anaerobic growth using sucrose. Therefore, the biomass-specific conversion rates were obtained using oxygen-limited continuous cultures. Strains performing the breakdown of the sucrose by hydrolysis (SUC-HYD) or phosphorolysis (SUC-PHOSP) were studied in such conditions. An experimentally validated in silico model, modified to account for plasmid and protein burdens, was employed to calculate carbon and electron consistent conversion rates. In both strains, the biomass yields were lower than expected and, strikingly, SUC-PHOSP showed a yield lower than SUC-HYD. Flux balance analyses indicated a significant increase in the non-growth-associated ATP expenses by comparison with the growth on glucose. The observed fructose-1,6-biphosphatase and phosphoglucomutase activities, as well as the concentrations of glycogen, suggest the operation of ATP futile cycles triggered by a combination of the oxygen limitation and the metabolites released during the sucrose breakdown.

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.