Towards kinetic modeling of genome-scale metabolic networks without sacrificing stoichiometric,thermodynamic and physiological constraints

Mathematical modeling is an essential tool for the comprehensive understanding of cell metabolism and its interactions with the environmental and process conditions. Recent developments in the construction and analysis of stoichiometric models made it possible to define limits on steady-state metabolic behavior using flux balance analysis. However, detailed information on enzyme kinetics and enzyme regulation is needed to formulate kinetic models that can accurately capture the dynamic metabolic responses. The use of mechanistic enzyme kinetics is a difficult task due to uncertainty in the kinetic properties of enzymes. Therefore, the majority of recent works considered only mass action kinetics for reactions in metabolic networks. Herein, we applied the optimization and risk analysis of complex living entities (ORACLE) framework and constructed a large-scale mechanistic kinetic model of optimally grown Escherichia coli. We investigated the complex interplay between stoichiometry, thermodynamics, and kinetics in determining the flexibility and capabilities of metabolism. Our results indicate that enzyme saturation is a necessary consideration in modeling metabolic networks and it extends the feasible ranges of metabolic fluxes and metabolite concentrations. Our results further suggest that enzymes in metabolic networks have evolved to function at different saturation states to ensure greater flexibility and robustness of cellular metabolism.     

High-level heterologous production of propionate in engineered Escherichia coli

A propanologenic (i.e., 1-propanol-producing) bacterium Escherichia coli strain was previously derived by activating the genomic sleeping beauty mutase (Sbm) operon. The activated Sbm pathway branches out of the tricarboxylic acid (TCA) cycle at the succinyl-CoA node to form propionyl-CoA and its derived metabolites of 1-propanol and propionate. In this study, we targeted several TCA cycle genes encoding enzymes near the succinyl-CoA node for genetic manipulation to identify the individual contribution of the carbon flux into the Sbm pathway from the three TCA metabolic routes, that is, oxidative TCA cycle, reductive TCA branch, and glyoxylate shunt. For the control strain CPC-Sbm, in which propionate biosynthesis occurred under relatively anaerobic conditions, the carbon flux into the Sbm pathway was primarily derived from the reductive TCA branch, and both succinate availability and the SucCD-mediated interconversion of succinate/succinyl-CoA were critical for such carbon flux redirection. Although the oxidative TCA cycle normally had a minimal contribution to the carbon flux redirection, the glyoxylate shunt could be an alternative and effective carbon flux contributor under aerobic conditions. With mechanistic understanding of such carbon flux redirection, metabolic strategies based on blocking the oxidative TCA cycle (via ∆sdhA mutation) and deregulating the glyoxylate shunt (via ∆iclR mutation) were developed to enhance the carbon flux redirection and therefore propionate biosynthesis, achieving a high propionate titer of 30.9 g/L with an overall propionate yield of 49.7% upon fed-batch cultivation of the double mutant strain CPC-Sbm∆sdhA∆iclR under aerobic conditions. The results also suggest that the Sbm pathway could be metabolically active under both aerobic and anaerobic conditions.     

Metabolic adaptation of Escherichia coli during temperature-induced recombinant protein production: 2. Redirection of metabolic fluxes

The impact of temperature-induced synthesis of human basic fibroblast growth factor (hFGF-2) in high-cell-density cultures of recombinant Escherichia coli was studied by estimating metabolic flux variations. Metabolic flux distributions in E. coli were calculated by means of a stoichiometric network and linear programming. After the temperature upshift, a substantially elevated energy demand for synthesis of hFGF-2 and heat shock proteins resulted in a redirection of metabolic fluxes. Catabolic pathways like the Embden-Meyerhof-Parnas pathway and the tricarboxylic acid (TCA) cycle showed significantly enhanced activities, leading to reduced flux to growth-associated pathways like the pentose phosphate pathway and other anabolic pathways. Upon temperature upshift, an excess of NADPH was produced in the TCA cycle by isocitrate dehydrogenase. The metabolic model predicted the involvement of a transhydrogenase generating additional NADH from NADPH, thereby increasing ATP regeneration in the respiratory chain. The influence of the temperature upshift on the host's metabolism was investigated by means of a control strain harboring the "empty" parental expression vector. The metabolic fluxes after the temperature upshift were redirected similarly to the production strain; the effects, however, were observed to a lesser extent and with different time profiles.     

Reconstitution of oxaloacetate metabolism in the tricarboxylic acid cycle in Synechocystis sp. PCC 6803: discovery of important factors that directly affect the conversion of oxaloacetate

The tricarboxylic acid (TCA) cycle is one of the most important metabolic pathways in nature. Oxygenic photoautotrophic bacteria, cyanobacteria, have an unusual TCA cycle. The TCA cycle in cyanobacteria contains two unique enzymes that are not part of the TCA cycle in other organisms. In recent years, sustainable metabolite production from carbon dioxide using cyanobacteria has been looked at as a means to reduce the environmental burden of this gas. Among cyanobacteria, the unicellular cyanobacterium Synechocystis sp. PCC 6803 (Synechocystis 6803) is an optimal host for sustainable metabolite production. Recently, metabolite production using the TCA cycle in Synechocystis 6803 has been carried out. Previous studies revealed that the branch point of the oxidative and reductive TCA cycles, oxaloacetate metabolism, plays a key role in metabolite production. However, the biochemical mechanisms regulating oxaloacetate metabolism in Synechocystis 6803 are poorly understood. Concentrations of oxaloacetate in Synechocystis 6803 are extremely low, such that in vivo analysis of oxaloacetate metabolism does not seem realistic. Therefore, using purified enzymes, we reconstituted oxaloacetate metabolism in Synechocystis 6803 in vitro to reveal the regulatory mechanisms involved. Reconstitution of oxaloacetate metabolism revealed that pH, Mg²⁺ and phosphoenolpyruvate are important factors affecting the conversion of oxaloacetate in the TCA cycle. Biochemical analyses of the enzymes involved in oxaloacetate metabolism in this and previous studies revealed the biochemical mechanisms underlying the effects of these factors on oxaloacetate conversion. In addition, we clarified the function of two l- malate dehydrogenase isozymes in oxaloacetate metabolism. These findings serve as a basis for various applications of the cyanobacterial TCA cycle.     

Metabolic analysis of antibody producing CHO cells in fed‐batch production

Chinese hamster ovary (CHO) cells are commonly used for industrial production of recombinant proteins in fed batch or alternative production systems. Cells progress through multiple metabolic stages during fed‐batch antibody (mAb) production, including an exponential growth phase accompanied by lactate production, a low growth, or stationary phase when specific mAb production increases, and a decline when cell viability declines. Although media composition and cell lineage have been shown to impact growth and productivity, little is known about the metabolic changes at a molecular level. Better understanding of cellular metabolism will aid in identifying targets for genetic and metabolic engineering to optimize bioprocess and cell engineering. We studied a high expressing recombinant CHO cell line, designated high performer (HP), in fed‐batch productions using stable isotope tracers and biochemical methods to determine changes in central metabolism that accompany growth and mAb production. We also compared and contrasted results from HP to a high lactate producing cell line that exhibits poor growth and productivity, designated low performer (LP), to determine intrinsic metabolic profiles linked to their respective phenotypes. Our results reveal alternative metabolic and regulatory pathways for lactate and TCA metabolite production to those reported in the literature. The distribution of key media components into glycolysis, TCA cycle, lactate production, and biosynthetic pathways was shown to shift dramatically between exponential growth and stationary (production) phases. We determined that glutamine is both utilized more efficiently than glucose for anaplerotic replenishment and contributes more significantly to lactate production during the exponential phase. Cells shifted to glucose utilization in the TCA cycle as growth rate decreased. The magnitude of this metabolic switch is important for attaining high viable cell mass and antibody titers. We also found that phosphoenolpyruvate carboxykinase (PEPCK1) and pyruvate kinase (PK) are subject to differential regulation during exponential and stationary phases. The concomitant shifts in enzyme expression and metabolite utilization profiles shed light on the regulatory links between cell metabolism, media metabolites, and cell growth. Biotechnol. Bioeng. 2013; 110: 1735–1747. © 2013 Wiley Periodicals, Inc.     

Reverse TCA cycle flux through isocitrate dehydrogenases 1 and 2 is required for lipogenesis in hypoxic melanoma cells

The tricarboxylic acid (TCA) cycle is the central hub of oxidative metabolism, running in the classic forward direction to provide carbon for biosynthesis and reducing agents for generation of ATP. Our metabolic tracer studies in melanoma cells showed that in hypoxic conditions the TCA cycle is largely disconnected from glycolysis. By studying the TCA branch point metabolites, acetyl CoA and citrate, as well as the metabolic endpoint glutamine and fatty acids, we developed a comprehensive picture of the rewiring of the TCA cycle that occurs in hypoxia. Hypoxic tumor cells maintain proliferation by running the TCA cycle in reverse. The source of carbon for acetyl CoA, citrate, and fatty acids switches from glucose in normoxia to glutamine in hypoxia. This hypoxic flux from glutamine into fatty acids is mediated by reductive carboxylation. This reductive carboxylation is catalyzed by two isocitrate dehydrogenases, IDH1 and IDH2. Their combined action is necessary and sufficient to effect the reverse TCA flux and maintain cellular viability.     

Metabolic engineering of Mannheimia succiniciproducens for malic acid production using dimethylsulfoxide as an electron acceptor

Microbial production of various TCA intermediates and related chemicals through the reductive TCA cycle has been of great interest. However, rumen bacteria that naturally possess strong reductive TCA cycle have been rarely studied to produce these chemicals, except for succinic acid, due to their dependence on fumarate reduction to transport electrons for ATP synthesis. In this study, malic acid (MA), a dicarboxylic acid of industrial importance, was selected as a target chemical for mass production using Mannheimia succiniciproducens, a rumen bacterium possessing a strong reductive branch of the TCA cycle. The metabolic pathway was reconstructed by eliminating fumarase to prevent MA conversion to fumarate. The respiration system of M. succiniciproducens was reconstructed by introducing the Actinobacillus succinogenes dimethylsulfoxide (DMSO) reductase to improve cell growth using DMSO as an electron acceptor. Also, the cell membrane was engineered by employing Pseudomonas aeruginosa cis-trans isomerase to enhance MA tolerance. High inoculum fed-batch fermentation of the final engineered strain produced 61 g/L of MA with an overall productivity of 2.27 g/L/h, which is the highest MA productivity reported to date. The systems metabolic engineering strategies reported in this study will be useful for developing anaerobic bioprocesses for the production of various industrially important chemicals.     

Kinetic-model-based pathway optimization with application to reverse glycolysis in mammalian cells

Over the last two decades, model-based metabolic pathway optimization tools have been developed for the design of microorganisms to produce desired metabolites. However, few have considered more complex cellular systems such as mammalian cells, which requires the use of nonlinear kinetic models to capture the effects of concentration changes and cross-regulatory interactions. In this study, we develop a new two-stage pathway optimization framework based on kinetic models that incorporate detailed kinetics and regulation information. In Stage 1, a set of optimization problems are solved to identify and rank the enzymes that contribute the most to achieving the metabolic objective. Stage 2 then determines the optimal enzyme interventions for specified desired numbers of enzyme adjustments. It also incorporates multi-scenario optimization, which allows the simultaneous consideration of multiple physiological conditions. We apply the proposed framework to find enzyme adjustments that enable a reverse glucose flow in cultured mammalian cells, thereby eliminating the need for glucose feed in the late culture stage and enhancing process robustness. The computational results demonstrate the efficacy of the proposed approach; it not only captures the important regulations and key enzymes for reverse glycolysis but also identifies differences and commonalities in the metabolic requirements for different carbon sources.     

Genome-scale analysis of Mannheimia succiniciproducens metabolism

Mannheimia succiniciproducens MBEL55E isolated from bovine rumen is a capnophilic gram-negative bacterium that efficiently produces succinic acid, an industrially important four carbon dicarboxylic acid. In order to design a metabolically engineered strain which is capable of producing succinic acid with high yield and productivity, it is essential to optimize the whole metabolism at the systems level. Consequently, in silico modeling and simulation of the genome-scale metabolic network was employed for genome-scale analysis and efficient design of metabolic engineering experiments. The genome-scale metabolic network of M. succiniciproducens consisting of 686 reactions and 519 metabolites was constructed based on reannotation and validation experiments. With the reconstructed model, the network structure and key metabolic characteristics allowing highly efficient production of succinic acid were deciphered; these include strong PEP carboxylation, branched TCA cycle, relative weak pyruvate formation, the lack of glyoxylate shunt, and non-PTS for glucose uptake. Constraints-based flux analyses were then carried out under various environmental and genetic conditions to validate the genome-scale metabolic model and to decipher the altered metabolic characteristics. Predictions based on constraints-based flux analysis were mostly in excellent agreement with the experimental data. In silico knockout studies allowed prediction of new metabolic engineering strategies for the enhanced production of succinic acid. This genome-scale in silico model can serve as a platform for the systematic prediction of physiological responses of M. succiniciproducens to various environmental and genetic perturbations and consequently for designing rational strategies for strain improvement.     

Metabolic control at the cytosol–mitochondria interface in different growth phases of CHO cells

Metabolism at the cytosol–mitochondria interface and its regulation is of major importance particularly for efficient production of biopharmaceuticals in Chinese hamster ovary (CHO) cells but also in many diseases. We used a novel systems-oriented approach combining dynamic metabolic flux analysis and determination of compartmental enzyme activities to obtain systems level information with functional, spatial and temporal resolution. Integrating these multiple levels of information, we were able to investigate the interaction of glycolysis and TCA cycle and its metabolic control. We characterized metabolic phases in CHO batch cultivation and assessed metabolic efficiency extending the concept of metabolic ratios. Comparing in situ enzyme activities including their compartmental localization with in vivo metabolic fluxes, we were able to identify limiting steps in glycolysis and TCA cycle. Our data point to a significant contribution of substrate channeling to glycolytic regulation. We show how glycolytic channeling heavily affects the availability of pyruvate for the mitochondria. Finally, we show that the activities of transaminases and anaplerotic enzymes are tailored to permit a balanced supply of pyruvate and oxaloacetate to the TCA cycle in the respective metabolic states. We demonstrate that knowledge about metabolic control can be gained by correlating in vivo metabolic flux dynamics with time and space resolved in situ enzyme activities.     

Ensemble modeling for strain development of l-lysine-producing Escherichia coli

One of the main strategies to improve the production of relevant metabolites has been the manipulation of single or multiple key genes in the metabolic pathways. This kind of strategy requires several rounds of experiments to identify enzymes that impact either yield or productivity. The use of mathematical tools to facilitate this process is desirable. In this work, we apply the Ensemble Modeling (EM) framework, which uses phenotypic data (effects of enzyme overexpression or genetic knockouts on the steady-state production rate) to screen for potential models capable of describe existing data and thus gaining insight to improve strains for l-lysine production. Described herein is a strategy to generate a set of kinetic models that describe a set of enzyme overexpression phenotypes previously determined in an Escherichia coli strain that produces increased levels of l-lysine in an industrial laboratory. This final ensemble of models captures the kinetic characteristics of the cell through screening of phenotypes after sequential overexpression of enzymes. Furthermore, these models demonstrate some predictive capability, as starting from the reference producing strain (overexpressing desensitized dihydrodipicolinate synthetase (dapA*)) this set of models is able to predict that the desensitization of aspartate kinase (lysC*) is the next rate-controlling step in the l-lysine pathway. Moreover, this set of models allows for the generation of further targets for testing, for example, phosphoenolpyruvate (Ppc), aspartate aminotransferase (AspC), and glutamate dehydrogenase (GdhA). This work demonstrates the usefulness, applicability, and scope that the Ensemble Modeling framework offers to build production strains.     

Enhancement of recombinant protein production in Escherichia coli by coproduction of aspartase

As commonly recognized, the excretion of acetate by the aerobic growth of Escherichia coli on glucose is a manifestation of imbalanced flux between glycolysis and the tricarboxylic acid (TCA) cycle. Accordingly, this may restrict the production of recombinant proteins in E. coli, due to the limited amounts of precursor metabolites produced in TCA cycle. To approach this issue, an extra supply of intermediate metabolites in TCA cycle was made by conversion of aspartate to fumarate, a reaction mediated by the activity of L-aspartate ammonia-lyase (aspartase). As a result, in the glucose minimal medium containing aspartate, the production of two recombinant proteins, beta-galactosidase and green fluorescent protein, in the aspartase-producing strain was substantially increased by 5-fold in association with 30-40% more biomass production. This preliminary study illustrates the great promise of this approach used to enhance the production of these two recombinant proteins.