Metabolite interactions in the bacterial Calvin cycle and implications for flux regulation.

Metabolite-level regulation of enzyme activity is important for microbes to cope with environmental shifts. Knowledge of such regulations can also guide strain engineering for biotechnology. Here we apply limited proteolysis-small molecule mapping (LiP-SMap) to identify and compare metabolite-protein interactions in the proteomes of two cyanobacteria and two lithoautotrophic bacteria that fix CO2 using the Calvin cycle. Clustering analysis of the hundreds of detected interactions shows that some metabolites interact in a species-specific manner. We estimate that approximately 35% of interacting metabolites affect enzyme activity in vitro, and the effect is often minor. Using LiP-SMap data as a guide, we find that the Calvin cycle intermediate glyceraldehyde-3-phosphate enhances activity of fructose-1,6/sedoheptulose-1,7-bisphosphatase (F/SBPase) from Synechocystis sp. PCC 6803 and Cupriavidus necator in reducing conditions, suggesting a convergent feed-forward activation of the cycle. In oxidizing conditions, glyceraldehyde-3-phosphate inhibits Synechocystis F/SBPase by promoting enzyme aggregation. In contrast, the glycolytic intermediate glucose-6-phosphate activates F/SBPase from Cupriavidus necator but not F/SBPase from Synechocystis. Thus, metabolite-level regulation of the Calvin cycle is more prevalent than previously appreciated.

Thermodynamic limitations of PHB production from formate and fructose in Cupriavidus necator.

The chemolithotroph Cupriavidus necator H16 is known as a natural producer of the bioplastic-polymer PHB, as well as for its metabolic versatility to utilize different substrates, including formate as the sole carbon and energy source. Depending on the entry point of the substrate, this versatility requires adjustment of the thermodynamic landscape to maintain sufficiently high driving forces for biological processes. Here we employed a model of the core metabolism of C. necator H16 to analyze the thermodynamic driving forces and PHB yields from formate for different metabolic engineering strategies. For this, we enumerated elementary flux modes (EFMs) of the network and evaluated their PHB yields as well as thermodynamics via Max-min driving force (MDF) analysis and random sampling of driving forces. A heterologous ATP:citrate lyase reaction was predicted to increase driving force for producing acetyl-CoA. A heterologous phosphoketolase reaction was predicted to increase maximal PHB yields as well as driving forces. These enzymes were then verified experimentally to enhance PHB titers between 60 and 300% in select conditions. The EFM analysis also revealed that PHB production from formate may be limited by low driving forces through citrate lyase and aconitase, as well as cofactor balancing, and identified additional reactions associated with low and high PHB yield. Proteomics analysis of the engineered strains confirmed an increased abundance of aconitase and cofactor balancing. The findings of this study aid in understanding metabolic adaptation. Furthermore, the outlined approach will be useful in designing metabolic engineering strategies in other non-model bacteria.

Protein allocation and utilization in the versatile chemolithoautotroph Cupriavidus necator.

Bacteria must balance the different needs for substrate assimilation, growth functions, and resilience in order to thrive in their environment. Of all cellular macromolecules, the bacterial proteome is by far the most important resource and its size is limited. Here, we investigated how the highly versatile 'knallgas' bacterium Cupriavidus necator reallocates protein resources when grown on different limiting substrates and with different growth rates. We determined protein quantity by mass spectrometry and estimated enzyme utilization by resource balance analysis modeling. We found that C. necator invests a large fraction of its proteome in functions that are hardly utilized. Of the enzymes that are utilized, many are present in excess abundance. One prominent example is the strong expression of CBB cycle genes such as Rubisco during growth on fructose. Modeling and mutant competition experiments suggest that CO2-reassimilation through Rubisco does not provide a fitness benefit for heterotrophic growth, but is rather an investment in readiness for autotrophy.

Kinetic modeling of the Calvin cycle identifies flux control and stable metabolomes in Synechocystis carbon fixation.

Biological fixation of atmospheric CO2 via the Calvin-Benson-Bassham cycle has massive ecological impact and offers potential for industrial exploitation, either by improving carbon fixation in plants and autotrophic bacteria, or by installation into new hosts. A kinetic model of the Calvin-Benson-Bassham cycle embedded in the central carbon metabolism of the cyanobacterium Synechocystis sp. PCC 6803 was developed to investigate its stability and underlying control mechanisms. To reduce the uncertainty associated with a single parameter set, random sampling of the steady-state metabolite concentrations and the enzyme kinetic parameters was employed, resulting in millions of parameterized models which were analyzed for flux control and stability against perturbation. Our results show that the Calvin cycle had an overall high intrinsic stability, but a high concentration of ribulose 1,5-bisphosphate was associated with unstable states. Low substrate saturation and high product saturation of enzymes involved in highly interconnected reactions correlated with increased network stability. Flux control, that is the effect that a change in one reaction rate has on the other reactions in the network, was distributed and mostly exerted by energy supply (ATP), but also by cofactor supply (NADPH). Sedoheptulose 1,7-bisphosphatase/fructose 1,6-bisphosphatase, fructose-bisphosphate aldolase, and transketolase had a weak but positive effect on overall network flux, in agreement with published observations. The identified flux control and relationships between metabolite concentrations and system stability can guide metabolic engineering. The kinetic model structure and parameterizing framework can be expanded for analysis of metabolic systems beyond the Calvin cycle.

Thermodynamic analysis of computed pathways integrated into the metabolic networks of E. coli and Synechocystis reveals contrasting expansion potential.

Introducing biosynthetic pathways into an organism is both reliant on and challenged by endogenous biochemistry. Here we compared the expansion potential of the metabolic network in the photoautotroph Synechocystis with that of the heterotroph E. coli using the novel workflow POPPY (Prospecting Optimal Pathways with PYthon). First, E. coli and Synechocystis metabolomic and fluxomic data were combined with metabolic models to identify thermodynamic constraints on metabolite concentrations (NET analysis). Then, thousands of automatically constructed pathways were placed within each network and subjected to a network-embedded variant of the max-min driving force analysis (NEM). We found that the networks had different capabilities for imparting thermodynamic driving forces toward certain compounds. Key metabolites were constrained differently in Synechocystis due to opposing flux directions in glycolysis and carbon fixation, the forked tri-carboxylic acid cycle, and photorespiration. Furthermore, the lysine biosynthesis pathway in Synechocystis was identified as thermodynamically constrained, impacting both endogenous and heterologous reactions through low 2-oxoglutarate levels. Our study also identified important yet poorly covered areas in existing metabolomics data and provides a reference for future thermodynamics-based engineering in Synechocystis and beyond. The POPPY methodology represents a step in making optimal pathway-host matches, which is likely to become important as the practical range of host organisms is diversified.