Search

In the last years, Salmonella has been extensively studied not only due to its importance as a pathogen, but also as a host to produce pharmaceutical compounds. However, the full exploitation of Salmonella as a platform for bioproduct delivery has been hampered by the lack of information about its metabolism. Genome-scale metabolic models can be valuable tools to delineate metabolic engineering strategies as long as they closely represent the actual metabolism of the target organism. In the present study, a 13C-MFA approach was applied to map the fluxes at the central carbon pathways of S. typhimurium LT2 growing at glucose-limited chemostat cultures. The experiments were carried out in a 2L bioreactor, using defined medium enriched with 20% 13C-labeled glucose. Metabolic flux distributions in central carbon pathways of S. typhimurium LT2 were estimated using OpenFLUX2 based on the labeling pattern of biomass protein hydrolysates together with biomass composition. The results suggested that pentose phosphate is used to catabolize glucose, with minor fluxes through glycolysis. In silico simulations, using Optflux and pFBA as simulation method, allowed to study the performance of the genome-scale metabolic model. In general, the accuracy of in silico simulations was improved by the superimposition of estimated intracellular fluxes to the existing genome-scale metabolic model, showing a better fitting to the experimental extracellular fluxes, whereas the intracellular fluxes of pentose phosphate and anaplerotic reactions were poorly described.

Volatile fatty acids (VFAs) are an inexpensive and renewable carbon source that can be generated from gas fermentation and anaerobic digestion of fermentable wastes. The oleaginous yeast Yarrowia lipolytica is a promising biocatalyst that can utilize VFAs and convert them into triacylglycerides (TAGs). However, currently there is limited knowledge on the metabolism of Y. lipolytica when cultured on VFAs. To develop a better understanding, we used acetate as the sole carbon source to culture two strains, a control strain and a previously engineered strain for lipid overaccumulation. For both strains, metabolism during the growth phase and lipid production phase were investigated by metabolic flux analysis using two parallel sodium acetate tracers. The resolved flux distributions demonstrate that the glyoxylate shunt pathway is constantly active and the flux through gluconeogenesis varies depending on strain and phase. In particular, by regulating the activities of malate transport and pyruvate kinase, the cells divert only a portion of the glyoxylate shunt flux required to satisfy the needs for anaplerotic reactions and NADPH production through gluconeogenesis and the oxidative pentose phosphate pathway (PPP). Excess flux flows back to the tricarboxylic acid (TCA) cycle for energy production. As with the case of glucose as the substrate, the primary source for lipogenic NADPH is derived from the oxidative PPP.

In contrast to stoichiometric-based models, the development of large-scale kinetic models of metabolism has been hindered by the challenge of identifying kinetic parameter values and kinetic rate laws applicable to a wide range of environmental and/or genetic perturbations. The recently introduced ensemble modeling (EM) procedure provides a promising remedy to address these challenges by decomposing metabolic reactions into elementary reaction steps and incorporating all phenotypic observations, upon perturbation, in its model parameterization scheme. Here, we present a kinetic model of Escherichia coli core metabolism that satisfies the fluxomic data for wild-type and seven mutant strains by making use of the EM concepts. This model encompasses 138 reactions, 93 metabolites and 60 substrate-level regulatory interactions accounting for glycolysis/gluconeogenesis, pentose phosphate pathway, TCA cycle, major pyruvate metabolism, anaplerotic reactions and a number of reactions in other parts of the metabolism. Parameterization is performed using a formal optimization approach that minimizes the discrepancies between model predictions and flux measurements. The predicted fluxes by the model are within the uncertainty range of experimental flux data for 78% of the reactions (with measured fluxes) for both the wild-type and seven mutant strains. The remaining flux predictions are mostly within three standard deviations of reported ranges. Converting the EM-based parameters into a Michaelis-Menten equivalent formalism revealed that 35% of Km and 77% of kcat parameters are within uncertainty range of the literature-reported values. The predicted metabolite concentrations by the model are also within uncertainty ranges of metabolomic data for 68% of the metabolites. A leave-one-out cross-validation test to evaluate the flux prediction performance of the model showed that metabolic fluxes for the mutants located in the proximity of mutations used for training the model can be predicted more accurately. The constructed model and the parameterization procedure presented in this study pave the way for the construction of larger-scale kinetic models with more narrowly distributed parameter values as new metabolomic/fluxomic data sets are becoming available for E. coli and other organisms.

Corynebacterium glutamicum possesses high in vivo activity of the gluconeogenic phosphoenolpyruvate carboxykinase (PEPCk) during growth on glucose, resulting together with anaplerotic carboxylation reactions in a PEP/pyruvate/oxaloacetate substrate cycle. The present study investigated the changes in intracellular fluxes and metabolite concentrations that are caused by altered PEPCk activity in l -lysine-producing C. glutamicum MH20-22B, applying a recently developed 13C labeling-based strategy for anaplerotic flux resolution and quantification. Abolition of PEPCk activity by deletion of the respective pck gene resulted in increased intracellular concentrations of oxaloacetate, l -aspartate, α-ketoglutarate, pyruvate, and l -lysine and in a 60% enhanced flux toward l -lysine biosynthesis, whereas increasing the PEPCk activity by pck overexpression had opposite effects. The results of the combined measurements of enzyme activities, in vivo fluxes, and metabolite concentrations were exploited to elucidate the in vivo regulation of anaplerotic reactions in C. glutamicum, and implications for the metabolic engineering of amino-acid-producing strains are discussed.

The extensive use of13C enrichments in precursor metabolites for flux quantification does not rely on NADPH stoichiometries and can therefore be used to quantify reducing power fluxes. As an application of this concept, the NADPH fluxes were quantified in an l -lysine producer ofCorynebacterium glutamicumgrown into metabolic and isotopic steady state with [1-13C]glucose. In this case, where the organism's NADPH-dependent glutamate dehydrogenase consumes reducing power, the NADPH flux generated is 210% (molar flux relative to glucose uptake rate) with its major part (72% of the total) generated via the pentose phosphate pathway activity. An isogenic strain in which the glutamate dehydrogenase ofC. glutamicumwas replaced by the NADH-dependent glutamate dehydrogenase ofPeptostreptococcus asaccharolyticuswas made and the metabolite fluxes were again estimated. The major response to this local perturbation is a drastically reduced NADPH generation of only 139%. Most of the NADPH (62% of the total) is now generated via the tricarboxylic acid cycle activity. This shows the extraordinary flexibility of the central metabolism and provides a picture of the global regulatory properties of the central metabolism. Furthermore, a detailed analysis of the fluxes and exchange fluxes within the anaplerotic reactions is given. It is hypothesized that these reactions might also serve to balance the total reducing power budget as well as the energy budget within the cell.