Design and thermodynamic analysis of a pathway enabling anaerobic production of poly-3-hydroxybutyrate in Escherichia coli

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.

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 (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.

An NADH preferring acetoacetyl-CoA reductase is engaged in poly-3-hydroxybutyrate accumulation in Escherichia coli

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.

Quantifying NAD(P)H production in the upper Entner-Doudoroff pathway from Pseudomonas putida KT2440

•

The first kinetic characterization of PputG6PDH-1 is presented.

•

The relative production of NADH and NADPH by PputG6PDH-1 is quantified.

•

The stoichiometric matrix of in silico metabolic models for Pseudomonas putida must be modified.

Despite the lack of biochemical information, all available in silico metabolic models of Pseudomonas putida KT2440 consider NADP as the only cofactor accepted by the glucose-6-phosphate dehydrogenases. Because the Entner–Doudoroff pathway is the main glycolytic route in this bacterium, determining how much NADH and NADPH are produced in the reaction catalyzed by these enzymes is very important for the correct interpretation of metabolic flux distributions. To determine the actual cofactor preference of the glucose-6-phosphate dehydrogenase encoded by the zwf-1 gene (PputG6PDH-1), the major isoform during growth on glucose, we purified this protein and studied its kinetic properties. Based on simple kinetic principles, we estimated the in vivo relative production of NADH and NADPH during the oxidation of glucose-6-phosphate (G6P). Contrary to the general assumption, our calculations showed that the reaction catalyzed by PputG6PDH-1 yields around 1/3 mol of NADPH and 2/3 mol of NADH per mol of oxidized G6P. Additionally, we obtained data suggesting that the reaction catalyzed by the 6-phosphogluconate dehydrogenase is active during growth on glucose, and it also produces NADH. These results indicate that the stoichiometric matrix of in silico models of P. putida KT2440 must be corrected and highlight the importance of considering the physiological concentrations of the involved metabolites to estimate the actual proportion of NADH and NADPH produced by a dehydrogenase.

Metabolic impact of an NADH-producing glucose-6-phosphate dehydrogenase in Escherichia coli

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 (Δpgi). 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 Δpgi 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 Δpgi 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.