A detailed metabolic flux analysis of an underdetermined network of CHO cells

CHOmpact: A reduced metabolic model of Chinese hamster ovary cells with enhanced interpretability

Metabolic modeling has emerged as a key tool for the characterization of biopharmaceutical cell culture processes. Metabolic models have also been instrumental in identifying genetic engineering targets and developing feeding strategies that optimize the growth and productivity of Chinese hamster ovary (CHO) cells. Despite their success, metabolic models of CHO cells still present considerable challenges. Genome-scale metabolic models (GeMs) of CHO cells are very large (>6000 reactions) and are difficult to constrain to yield physiologically consistent flux distributions. The large scale of GeMs also makes the interpretation of their outputs difficult. To address these challenges, we have developed CHOmpact, a reduced metabolic network that encompasses 101 metabolites linked through 144 reactions. Our compact reaction network allows us to deploy robust, nonlinear optimization and ensure that the computed flux distributions are physiologically consistent. Furthermore, our CHOmpact model delivers enhanced interpretability of simulation results and has allowed us to identify the mechanisms governing shifts in the anaplerotic consumption of asparagine and glutamate as well as an important mechanism of ammonia detoxification within mitochondria. CHOmpact, thus, addresses key challenges of large-scale metabolic models and will serve as a platform to develop dynamic metabolic models for the control and optimization of biopharmaceutical cell culture processes.

MYC overrides HIF-1α to regulate proliferating primary cell metabolism in hypoxia

Hypoxia requires metabolic adaptations to sustain energetically demanding cellular activities. While the metabolic consequences of hypoxia have been studied extensively in cancer cell models, comparatively little is known about how primary cell metabolism responds to hypoxia. Thus, we developed metabolic flux models for human lung fibroblast and pulmonary artery smooth muscle cells proliferating in hypoxia. Unexpectedly, we found that hypoxia decreased glycolysis despite activation of hypoxia-inducible factor 1α (HIF-1α) and increased glycolytic enzyme expression. While HIF-1α activation in normoxia by prolyl hydroxylase (PHD) inhibition did increase glycolysis, hypoxia blocked this effect. Multi-omic profiling revealed distinct molecular responses to hypoxia and PHD inhibition, and suggested a critical role for MYC in modulating HIF-1α responses to hypoxia. Consistent with this hypothesis, MYC knockdown in hypoxia increased glycolysis and MYC over-expression in normoxia decreased glycolysis stimulated by PHD inhibition. These data suggest that MYC signaling in hypoxia uncouples an increase in HIF-dependent glycolytic gene transcription from glycolytic flux.

A genome-scale dynamic constraint-based modelling (gDCBM) framework predicts growth dynamics, medium composition and intracellular flux distributions in CHO clonal variations

Optimizing mammalian cell growth and bioproduction is a tedious task. However, due to the inherent complexity of eukaryotic cells, heuristic experimental approaches such as, metabolic engineering and bioprocess design, are frequently integrated with mathematical models of cell culture to improve biological process efficiency and find paths for improvement. Constraint-based metabolic models have evolved over the last two decades to be used for dynamic modelling in addition to providing a linear description of steady-state metabolic systems. Formulation and implementation of the underlying optimization problems require special attention to the model's performance and feasibility, lack of defects in the definition of system components, and consideration of optimal alternate solutions, in addition to processing power limitations. Here, the time-resolved dynamics of a genome-scale metabolic network of Chinese hamster ovary (CHO) cell metabolism are shown using a genome-scale dynamic constraint-based modelling framework (gDCBM). The metabolic network was adapted from a reference model of CHO genome-scale metabolic model (GSMM), iCHO_DG44_v1, and dynamic restrictions were imposed to its exchange fluxes based on experimental results. We used this framework for predicting physiological changes in CHO clonal variants. Because of the methodical creation of the components for the flux balance analysis optimization problem and the integration of a switch time, this model can generate sequential predictions of intracellular fluxes during growth and non-growth phases (per hour of culture time) and transparently reveal the shortcomings in such practice. As a result of the differences exploited by various clones, we can understand the relevance of changes in intracellular flux distribution and exometabolomics. The integration of various omics data into the given gDCBM framework, as well as the reductionist analysis of the model, can further help bioprocess optimization.

Metabolic Flux Analysis of VERO Cells under Various Culture Conditions

Although the culture of VERO cells in bioreactors is an important industrial bioprocess for the production of viruses and vaccines, surprisingly few reports on the analysis of the flux distribution in the cell metabolism have been published. In this study, an attempt is made to fill this gap by providing an analysis of relatively simple metabolic networks, which are constructed to describe the cell behavior in different culture conditions, e.g., the exponential growth phase (availability of glucose and glutamine), cell growth without glutamine, and cell growth without glucose and glutamine. The metabolic networks are kept as simple as possible in order to avoid underdeterminacy linked to the lack of extracellular measurements, and a unique flux distribution is computed in each case based on a mild assumption that the macromolecular composition of the cell is known. The result of this computation provides some insight into the metabolic changes triggered by the culture conditions, which could support the design of feedback control strategies in fed batch or perfusion bioreactors where the lactate concentration is measured online and regulated by controlling the delivery rates of glucose and, possibly, of some essential amino acids.

Metabolic trends of Chinese hamster ovary cells in biopharmaceutical production under batch and fed-batch conditions

Extensive knowledge of Chinese hamster ovary (CHO) cell metabolism is required to improve process productivity and culture performance in biopharmaceutical manufacturing. However, CHO cells show a dynamic metabolism during culturing in batch and fed-batch bioreactors. CHO cell metabolism is generally described as taking place in three stages: exponential growth phase, stationary phase, and death phase. This review aims to summarize the trends of central metabolism for CHO cells during each stage. Additional insights into how culture conditions are related to phase transitions and force metabolic rewiring are provided. Understanding of CHO cell metabolism lends itself to improving culture qualities by, for example, identifying sources of toxic byproducts and pathways for cellular engineering. In summary, this review describes the changes in CHO cell central metabolism over the course of the culture.

How to Tackle Underdeterminacy in Metabolic Flux Analysis? A Tutorial and Critical Review

Metabolic flux analysis is often (not to say almost always) faced with system underdeterminacy. Indeed, the linear algebraic system formed by the steady-state mass balance equations around the intracellular metabolites and the equality constraints related to the measurements of extracellular fluxes do not define a unique solution for the distribution of intracellular fluxes, but instead a set of solutions belonging to a convex polytope. Various methods have been proposed to tackle this underdeterminacy, including flux pathway analysis, flux balance analysis, flux variability analysis and sampling. These approaches are reviewed in this article and a toy example supports the discussion with illustrative numerical results.

Modelling Cell Metabolism: A Review on Constraint-Based Steady-State and Kinetic Approaches

Studying cell metabolism serves a plethora of objectives such as the enhancement of bioprocess performance, and advancement in the understanding of cell biology, of drug target discovery, and in metabolic therapy. Remarkable successes in these fields emerged from heuristics approaches, for instance, with the introduction of effective strategies for genetic modifications, drug developments and optimization of bioprocess management. However, heuristics approaches have showed significant shortcomings, such as to describe regulation of metabolic pathways and to extrapolate experimental conditions. In the specific case of bioprocess management, such shortcomings limit their capacity to increase product quality, while maintaining desirable productivity and reproducibility levels. For instance, since heuristics approaches are not capable of prediction of the cellular functions under varying experimental conditions, they may lead to sub-optimal processes. Also, such approaches used for bioprocess control often fail in regulating a process under unexpected variations of external conditions. Therefore, methodologies inspired by the systematic mathematical formulation of cell metabolism have been used to address such drawbacks and achieve robust reproducible results. Mathematical modelling approaches are effective for both the characterization of the cell physiology, and the estimation of metabolic pathways utilization, thus allowing to characterize a cell population metabolic behavior. In this article, we present a review on methodology used and promising mathematical modelling approaches, focusing primarily to investigate metabolic events and regulation. Proceeding from a topological representation of the metabolic networks, we first present the metabolic modelling approaches that investigate cell metabolism at steady state, complying to the constraints imposed by mass conservation law and thermodynamics of reactions reversibility. Constraint-based models (CBMs) are reviewed highlighting the set of assumed optimality functions for reaction pathways. We explore models simulating cell growth dynamics, by expanding flux balance models developed at steady state. Then, discussing a change of metabolic modelling paradigm, we describe dynamic kinetic models that are based on the mathematical representation of the mechanistic description of nonlinear enzyme activities. In such approaches metabolic pathway regulations are considered explicitly as a function of the activity of other components of metabolic networks and possibly far from the metabolic steady state. We have also assessed the significance of metabolic model parameterization in kinetic models, summarizing a standard parameter estimation procedure frequently employed in kinetic metabolic modelling literature. Finally, some optimization practices used for the parameter estimation are reviewed.

Dynamic Modeling of CHO Cell Metabolism Using the Hybrid Cybernetic Approach With a Novel Elementary Mode Analysis Strategy

Chinese hamster ovary (CHO) cell culture has a major importance on the production of biopharmaceuticals, including recombinant therapeutic proteins such as monoclonal antibodies (MAb). Mathematical modeling of biological systems can successfully assess metabolism complexity while providing logical and systematic methods for relevant genetic target and culture parameter identification toward cell growth and productivity improvements. Most modeling approaches on CHO cells have been performed under stationary constraints, and only a few dynamic models have been presented on simplified reaction sets, due to substantial overparameterization problems. The hybrid cybernetic modeling (HCM) approach has been recently used to describe the dynamic behavior by incorporating regulation between different metabolic states by elementary mode participation control, with sets of equations evaluated by objective functions. However, as metabolic networks evaluated are constructed toward a genomic scale, and cell compartmentalization is considered, identification of the active set becomes more difficult as EM number exponentially grows. Thus, the development of robust approaches for EM active set selection and analysis with smaller computational requirements is required to impulse the use of cybernetic modeling on larger up to genome-scale networks. In this report, a novel elementary mode selection strategy, based on a polar representation of the convex solution space is presented and coupled to a cybernetic approach to model the dynamic physiologic and metabolic behavior of CHO-S cell cultures. The proposed Polar Space Yield Analysis (PSYA) was compared to other reported elementary mode selection approaches derived from Common Metabolic Objective Analysis (CMOA) used in Flux Balance Analysis (FBA), Yield Space Analysis (YSA), and Lumped Yield Space Analysis (LYSA). For this purpose, exponential growth phase dynamic metabolic models were calculated using kinetic rate equations based on previously modeled growth parameters. Finally, complete culture dynamic metabolic flux models were constructed using the HCM approach with selected elementary mode sets. The yield space elementary mode- and the polar space elementary mode- hybrid cybernetic models presented the best fits and performances. Also, a flux reaction perturbation prediction approach based on the polar yield solution space resulted useful for metabolic network flux distribution capability analysis and identification of potential genetic modifications targets.

Novel column generation-based optimization approach for poly-pathway kinetic model applied to CHO cell culture

Mathematical modelling can provide precious tools for bioprocess simulation, prediction, control and optimization of mammalian cell-based cultures. In this paper we present a novel method to generate kinetic models of such cultures, rendering complex metabolic networks in a poly-pathway kinetic model. The model is based on subsets of elementary flux modes (EFMs) to generate macro-reactions. Thanks to our column generation-based optimization algorithm, the experimental data are used to identify the EFMs, which are relevant to the data. Here the systematic enumeration of all the EFMs is eliminated and a network including a large number of reactions can be considered. In particular, the poly-pathway model can simulate multiple metabolic behaviors in response to changes in the culture conditions.

We apply the method to a network of 126 metabolic reactions describing cultures of antibody-producing Chinese hamster ovary cells, and generate a poly-pathway model that simulates multiple experimental conditions obtained in response to variations in amino acid availability. A good fit between simulated and experimental data is obtained, rendering the variations in the growth, product, and metabolite uptake/secretion rates. The intracellular reaction fluxes simulated by the model are explored, linking variations in metabolic behavior to adaptations of the intracellular metabolism.

•

Novel method to model multiple states by a poly-pathway kinetic model.

•

EFMs relevant to data identified by column generation (CG)-based optimization.

•

CG optimization enables use of networks much larger than systematic enumeration.

•

A kinetic model simulates changes in metabolic rates linked to available amino acids.

•

The flux distribution of each metabolic state is visualized in the original network.

Theoretical model and characteristics of mitochondrial thermogenesis

Based on the first law of thermodynamics and the thermal diffusion equation, the deduced theoretical model of mitochondrial thermogenesis satisfies the Laplace equation and is a special case of the thermal diffusion equation. The model settles the long-standing question of the ability to increase cellular temperature by endogenous thermogenesis and explains the thermogenic characteristics of brown adipocytes. The model and calculations also suggest that the number of free available protons is the major limiting factor for endogenous thermogenesis and its speed.

The best models of metabolism

Biochemical systems are among of the oldest application areas of mathematical modeling. Spanning a time period of over one hundred years, the repertoire of options for structuring a model and for formulating reactions has been constantly growing, and yet, it is still unclear whether or to what degree some models are better than others and how the modeler is to choose among them. In fact, the variety of options has become overwhelming and difficult to maneuver for novices and experts alike. This review outlines the metabolic model design process and discusses the numerous choices for modeling frameworks and mathematical representations. It tries to be inclusive, even though it cannot be complete, and introduces the various modeling options in a manner that is as unbiased as that is feasible. However, the review does end with personal recommendations for the choices of default models. WIREs Syst Biol Med 2017, 9:e1391. doi: 10.1002/wsbm.1391 For further resources related to this article, please visit the WIREs website.

Selection of chemically defined media for CHO cell fed-batch culture processes

Two CHO cell clones derived from the same parental CHOBC® cell line and producing the same monoclonal antibody (BC-G, a low producing clone; BC-P, a high producing clone) were tested in four basal media in all possible combinations with three feeds (=12 conditions) in fed-batch cultures. Higher amino acid feeding did not always lead to higher mAb production. The two clones showed differences in cell physiology, metabolism and optimal medium-feed combinations. During the phase transitions of all cultures, cell metabolism showed a shift represented by lower specific consumption and production rates, except for the specific glucose consumption rate in cultures fed by Actifeed A/B. The BC-P clone fed by Actifeed A/B showed a threefold cell volume increase and an increase of the specific consumption rate of glucose in the stationary phase. Since feeding was based on glucose this resulted in accumulation of amino acids for this feed, while this did not occur for the poorer feed (EFA/B). The same feed also led to an increase of cell size for the BC-G clone, but to a lesser extent.

What can mathematical modelling say about CHO metabolism and protein glycosylation?

Chinese hamster ovary cells have been in the spotlight for process optimization in recent years, due to being the major, long established cell factory for the production of recombinant proteins. A deep, quantitative understanding of CHO metabolism and mechanisms involved in protein glycosylation has proven to be attainable through the development of high throughput technologies. Here we review the most notable accomplishments in the field of modelling CHO metabolism and protein glycosylation.

Selection of chemically defined media for CHO cell fed-batch culture processes

Two CHO cell clones derived from the same parental CHO^BC® cell line and producing the same monoclonal antibody (BC-G, a low producing clone; BC-P, a high producing clone) were tested in four basal media in all possible combinations with three feeds (=12 conditions) in fed-batch cultures. Higher amino acid feeding did not always lead to higher mAb production. The two clones showed differences in cell physiology, metabolism and optimal medium-feed combinations. During the phase transitions of all cultures, cell metabolism showed a shift represented by lower specific consumption and production rates, except for the specific glucose consumption rate in cultures fed by Actifeed A/B. The BC-P clone fed by Actifeed A/B showed a threefold cell volume increase and an increase of the specific consumption rate of glucose in the stationary phase. Since feeding was based on glucose this resulted in accumulation of amino acids for this feed, while this did not occur for the poorer feed (EFA/B). The same feed also led to an increase of cell size for the BC-G clone, but to a lesser extent.

A Consensus Genome-scale Reconstruction of Chinese Hamster Ovary Cell Metabolism

Chinese hamster ovary (CHO) cells dominate biotherapeutic protein production and are widely used in mammalian cell line engineering research. To elucidate metabolic bottlenecks in protein production and to guide cell engineering and bioprocess optimization, we reconstructed the metabolic pathways in CHO and associated them with >1,700 genes in the Cricetulus griseus genome. The genome-scale metabolic model based on this reconstruction, iCHO1766, and cell-line-specific models for CHO-K1, CHO-S, and CHO-DG44 cells provide the biochemical basis of growth and recombinant protein production. The models accurately predict growth phenotypes and known auxotrophies in CHO cells. With the models, we quantify the protein synthesis capacity of CHO cells and demonstrate that common bioprocess treatments, such as histone deacetylase inhibitors, inefficiently increase product yield. However, our simulations show that the metabolic resources in CHO are more than three times more efficiently utilized for growth or recombinant protein synthesis following targeted efforts to engineer the CHO secretory pathway. This model will further accelerate CHO cell engineering and help optimize bioprocesses.

PFA toolbox: a MATLAB tool for Metabolic Flux Analysis

BACKGROUND

Metabolic Flux Analysis (MFA) is a methodology that has been successfully applied to estimate metabolic fluxes in living cells. However, traditional frameworks based on this approach have some limitations, particularly when measurements are scarce and imprecise. This is very common in industrial environments. The PFA Toolbox can be used to face those scenarios.

RESULTS

Here we present the PFA (Possibilistic Flux Analysis) Toolbox for MATLAB, which simplifies the use of Interval and Possibilistic Metabolic Flux Analysis. The main features of the PFA Toolbox are the following: (a) It provides reliable MFA estimations in scenarios where only a few fluxes can be measured or those available are imprecise. (b) It provides tools to easily plot the results as interval estimates or flux distributions. (c) It is composed of simple functions that MATLAB users can apply in flexible ways. (d) It includes a Graphical User Interface (GUI), which provides a visual representation of the measurements and their uncertainty. (e) It can use stoichiometric models in COBRA format. In addition, the PFA Toolbox includes a User's Guide with a thorough description of its functions and several examples.

CONCLUSIONS

The PFA Toolbox for MATLAB is a freely available Toolbox that is able to perform Interval and Possibilistic MFA estimations.

Dynamic estimation of specific fluxes in metabolic networks using non-linear dynamic optimization

BACKGROUND

Metabolic network models describing the biochemical reaction network and material fluxes inside microorganisms open interesting routes for the model-based optimization of bioprocesses. Dynamic metabolic flux analysis (dMFA) has lately been studied as an extension of regular metabolic flux analysis (MFA), rendering a dynamic view of the fluxes, also in non-stationary conditions. Recent dMFA implementations suffer from some drawbacks, though. More specifically, the fluxes are not estimated as specific fluxes, which are more biologically relevant. Also, the flux profiles are not smooth, and additional constraints like, e.g., irreversibility constraints on the fluxes, cannot be taken into account. Finally, in all previous methods, a basis for the null space of the stoichiometric matrix, i.e., which set of free fluxes is used, needs to be chosen. This choice is not trivial, and has a large influence on the resulting estimates.

RESULTS

In this work, a new methodology based on a B-spline parameterization of the fluxes is presented. Because of the high degree of non-linearity due to this parameterization, an incremental knot insertion strategy has been devised, resulting in a sequence of non-linear dynamic optimization problems. These are solved using state-of-the-art dynamic optimization methods and tools, i.e., orthogonal collocation, an interior-point optimizer and automatic differentiation. Also, a procedure to choose an optimal basis for the null space of the stoichiometric matrix is described, discarding the need to make a choice beforehand. The proposed methodology is validated on two simulated case studies: (i) a small-scale network with 7 fluxes, to illustrate the operation of the algorithm, and (ii) a medium-scale network with 68 fluxes, to show the algorithm's capabilities for a realistic network. The results show an accurate correspondence to the reference fluxes used to simulate the measurements, both in a theoretically ideal setting with no experimental noise, and in a realistic noise setting.

CONCLUSIONS

Because, apart from a metabolic reaction network and the measurements, no extra input needs to be given, the resulting algorithm is a systematic, integrated and accurate methodology for dynamic metabolic flux analysis that can be run online in real-time if necessary.

On dynamically generating relevant elementary flux modes in a metabolic network using optimization

Elementary flux modes (EFMs) are pathways through a metabolic reaction network that connect external substrates to products. Using EFMs, a metabolic network can be transformed into its macroscopic counterpart, in which the internal metabolites have been eliminated and only external metabolites remain. In EFMs-based metabolic flux analysis (MFA) experimentally determined external fluxes are used to estimate the flux of each EFM. It is in general prohibitive to enumerate all EFMs for complex networks, since the number of EFMs increases rapidly with network complexity. In this work we present an optimization-based method that dynamically generates a subset of EFMs and solves the EFMs-based MFA problem simultaneously. The obtained subset contains EFMs that contribute to the optimal solution of the EFMs-based MFA problem. The usefulness of our method was examined in a case-study using data from a Chinese hamster ovary cell culture and two networks of varied complexity. It was demonstrated that the EFMs-based MFA problem could be solved at a low computational cost, even for the more complex network. Additionally, only a fraction of the total number of EFMs was needed to compute the optimal solution.

DRUM: a new framework for metabolic modeling under non-balanced growth. Application to the carbon metabolism of unicellular microalgae

Metabolic modeling is a powerful tool to understand, predict and optimize bioprocesses, particularly when they imply intracellular molecules of interest. Unfortunately, the use of metabolic models for time varying metabolic fluxes is hampered by the lack of experimental data required to define and calibrate the kinetic reaction rates of the metabolic pathways. For this reason, metabolic models are often used under the balanced growth hypothesis. However, for some processes such as the photoautotrophic metabolism of microalgae, the balanced-growth assumption appears to be unreasonable because of the synchronization of their circadian cycle on the daily light. Yet, understanding microalgae metabolism is necessary to optimize the production yield of bioprocesses based on this microorganism, as for example production of third-generation biofuels. In this paper, we propose DRUM, a new dynamic metabolic modeling framework that handles the non-balanced growth condition and hence accumulation of intracellular metabolites. The first stage of the approach consists in splitting the metabolic network into sub-networks describing reactions which are spatially close, and which are assumed to satisfy balanced growth condition. The left metabolites interconnecting the sub-networks behave dynamically. Then, thanks to Elementary Flux Mode analysis, each sub-network is reduced to macroscopic reactions, for which simple kinetics are assumed. Finally, an Ordinary Differential Equation system is obtained to describe substrate consumption, biomass production, products excretion and accumulation of some internal metabolites. DRUM was applied to the accumulation of lipids and carbohydrates of the microalgae Tisochrysis lutea under day/night cycles. The resulting model describes accurately experimental data obtained in day/night conditions. It efficiently predicts the accumulation and consumption of lipids and carbohydrates.

The avian cell line AGE1.CR.pIX characterized by metabolic flux analysis

Background

In human vaccine manufacturing some pathogens such as Modified Vaccinia Virus Ankara, measles, mumps virus as well as influenza viruses are still produced on primary material derived from embryonated chicken eggs. Processes depending on primary cell culture, however, are difficult to adapt to modern vaccine production. Therefore, we derived previously a continuous suspension cell line, AGE1.CR.pIX, from muscovy duck and established chemically-defined media for virus propagation.

Results

To better understand vaccine production processes, we developed a stoichiometric model of the central metabolism of AGE1.CR.pIX cells and applied flux variability and metabolic flux analysis. Results were compared to literature dealing with mammalian and insect cell culture metabolism focusing on the question whether cultured avian cells differ in metabolism. Qualitatively, the observed flux distribution of this avian cell line was similar to distributions found for mammalian cell lines (e.g. CHO, MDCK cells). In particular, glucose was catabolized inefficiently and glycolysis and TCA cycle seem to be only weakly connected.

Conclusions

A distinguishing feature of the avian cell line is that glutaminolysis plays only a minor role in energy generation and production of precursors, resulting in low extracellular ammonia concentrations. This metabolic flux study is the first for a continuous avian cell line. It provides a basis for further metabolic analyses to exploit the biotechnological potential of avian and vertebrate cell lines and to develop specific optimized cell culture processes, e.g. vaccine production processes.

An in-silico study of the regulation of CHO cells glycolysis

In this work, a kinetic-metabolic model previously developed for CHO cells is used to study glycolysis regulation. The model is assessed for its biological relevance by analyzing its ability to simulate metabolic events induced following a hypoxic perturbation. Feedback and feedforward regulatory mechanisms known to occur to either inhibit or activate fluxes of glycolysis, are implemented in various combined scenarios and their effects on the metabolic response were analyzed. This study aims at characterizing the role of intermediates of glycolysis and of the cell energetic state, described as the AMP-to-ATP ratio, as inhibitors and activators of glycolysis pathway. In addition to the glycolysis pathway, we here describe the transient metabolic response of pathways that are connected to glycolysis, such as the pentose phosphate pathway, TCA cycle, cell bioenergetics system, glutamine and amino acids metabolisms. Taken individually, each regulatory mechanism leads to an oscillatory behavior in response to a hypoxic perturbation, while their combination clearly damps oscillations. However, only the addition of the cell energetic state to the regulatory mechanisms results in a non-oscillating response leading to metabolic flux rate rearrangement corresponding to the anaerobic metabolism expected to prevail under hypoxic conditions. We thus demonstrate in this work, from model simulations, that the robustness of a cell energetic metabolism can be described from a combination of feedback and feedforward inhibition and activation regulatory mechanisms of glycolysis fluxes, involving intermediates of glycolysis and the cell energetic state itself.

Non-stationary 13C metabolic flux analysis of Chinese hamster ovary cells in batch culture using extracellular labeling highlights metabolic reversibility and compartmentation

BACKGROUND

Mapping the intracellular fluxes for established mammalian cell lines becomes increasingly important for scientific and economic reasons. However, this is being hampered by the high complexity of metabolic networks, particularly concerning compartmentation.

RESULTS

Intracellular fluxes of the CHO-K1 cell line central carbon metabolism were successfully determined for a complex network using non-stationary 13C metabolic flux analysis. Mass isotopomers of extracellular metabolites were determined using [U-13C6] glucose as labeled substrate. Metabolic compartmentation and extracellular transport reversibility proved essential to successfully reproduce the dynamics of the labeling patterns. Alanine and pyruvate reversibility changed dynamically even if their net production fluxes remained constant. Cataplerotic fluxes of cytosolic phosphoenolpyruvate carboxykinase and mitochondrial malic enzyme and pyruvate carboxylase were successfully determined. Glycolytic pyruvate channeling to lactate was modeled by including a separate pyruvate pool. In the exponential growth phase, alanine, glycine and glutamate were excreted, and glutamine, aspartate, asparagine and serine were taken up; however, all these amino acids except asparagine were exchanged reversibly with the media. High fluxes were determined in the pentose phosphate pathway and the TCA cycle. The latter was fueled mainly by glucose but also by amino acid catabolism.

CONCLUSIONS

The CHO-K1 central metabolism in controlled batch culture proves to be robust. It has the main purpose to ensure fast growth on a mixture of substrates and also to mitigate oxidative stress. It achieves this by using compartmentation to control NADPH and NADH availability and by simultaneous synthesis and catabolism of amino acids.

Non-stationary 13C metabolic flux analysis of Chinese hamster ovary cells in batch culture using extracellular labeling highlights metabolic reversibility and compartmentation

BACKGROUND

Mapping the intracellular fluxes for established mammalian cell lines becomes increasingly important for scientific and economic reasons. However, this is being hampered by the high complexity of metabolic networks, particularly concerning compartmentation.

RESULTS

Intracellular fluxes of the CHO-K1 cell line central carbon metabolism were successfully determined for a complex network using non-stationary 13C metabolic flux analysis. Mass isotopomers of extracellular metabolites were determined using [U-13C6] glucose as labeled substrate. Metabolic compartmentation and extracellular transport reversibility proved essential to successfully reproduce the dynamics of the labeling patterns. Alanine and pyruvate reversibility changed dynamically even if their net production fluxes remained constant. Cataplerotic fluxes of cytosolic phosphoenolpyruvate carboxykinase and mitochondrial malic enzyme and pyruvate carboxylase were successfully determined. Glycolytic pyruvate channeling to lactate was modeled by including a separate pyruvate pool. In the exponential growth phase, alanine, glycine and glutamate were excreted, and glutamine, aspartate, asparagine and serine were taken up; however, all these amino acids except asparagine were exchanged reversibly with the media. High fluxes were determined in the pentose phosphate pathway and the TCA cycle. The latter was fueled mainly by glucose but also by amino acid catabolism.

CONCLUSIONS

The CHO-K1 central metabolism in controlled batch culture proves to be robust. It has the main purpose to ensure fast growth on a mixture of substrates and also to mitigate oxidative stress. It achieves this by using compartmentation to control NADPH and NADH availability and by simultaneous synthesis and catabolism of amino acids.

Dynamics of growth and metabolism controlled by glutamine availability in Chinese hamster ovary cells

The physiology of animal cells is characterized by constantly changing environmental conditions and adapting cellular responses. Applied dynamic metabolic flux analysis captures metabolic dynamics and can be applied to industrially relevant cultivation conditions. We investigated the impact of glutamine availability or limitation on the physiology of CHO K1 cells in eight different batch and fed-batch cultivations. Varying glutamine availability resulted in global metabolic changes. We observed dose-dependent effects of glutamine in batch cultivation. Identifying metabolic links from the glutamine metabolism to specific metabolic pathways, we show that glutamine feeding results in its coupling to tricarboxylic acid cycle fluxes and in its decoupling from metabolic waste production. We provide a mechanistic explanation of the cellular responses upon mild or severe glutamine limitation and ammonia stress. The growth rate of CHO K1 decreased with increasing ammonia levels in the supernatant. On the other hand, growth, especially culture longevity, was stimulated at mild glutamine-limiting conditions. Flux rearrangements in the pyruvate and amino acid metabolism compensate glutamine limitation by consumption of alternative carbon sources and facilitating glutamine synthesis and mitigate ammonia stress as result of glutamine abundance.

Factors causing compositional changes in soy protein hydrolysates and effects on cell culture functionality

Soy protein hydrolysates significantly enhance cell growth and recombinant protein production in cell cultures. The extent of this enhancement in cell growth and IgG production is known to vary from batch to batch. This can be due to differences in the abundance of different classes of compounds (e.g., peptide content), the quality of these compounds (e.g., glycated peptides), or the presence of specific compounds (e.g., furosine). These quantitative and qualitative differences between batches of hydrolysates result from variation in the seed composition and seed/meal processing. Although a considerable amount of literature is available that describes these factors, this knowledge has not been combined in an overview yet. The aim of this review is to identify the most dominant factors that affect hydrolysate composition and functionality. Although there is a limited influence of variation in the seed composition, the overview shows that the qualitative changes in hydrolysate composition result in the formation of minor compounds (e.g., Maillard reaction products). In pure systems, these compounds have a profound effect on the cell culture functionality. This suggests that the presence of these compounds in soy protein hydrolysates may affect hydrolysate functionality as well. This influence on the functionality can be of direct or indirect nature. For instance, some minor compounds (e.g., Maillard reaction products) are cytotoxic, whereas other compounds (e.g., phytates) suppress protein hydrolysis during hydrolysate production, resulting in altered peptide composition, and, thus, affect the functionality.

Isotopically nonstationary 13C flux analysis of Myc-induced metabolic reprogramming in B-cells

We assessed several methods of (13)C metabolic flux analysis (MFA) and found that isotopically nonstationary MFA achieved maximum flux resolution in cultured P493-6 B-cells, which have been engineered to provide tunable expression of the Myc oncoprotein. Comparison of metabolic flux maps obtained under oncogenic (High) and endogenous (Low) Myc expression levels revealed network-wide reprogramming in response to ectopic Myc expression. High Myc cells relied more heavily on mitochondrial oxidative metabolism than Low Myc cells and globally upregulated their consumption of amino acids relative to glucose. TCA cycle and amphibolic mitochondrial pathways exhibited 2- to 4-fold flux increases in High Myc cells, in contrast to modest increases in glucose uptake and lactate excretion. Because our MFA approach relied exclusively upon isotopic measurements of protein-bound amino acids and RNA-bound ribose, it is readily applicable to more complex tumor models that are not amenable to direct extraction and isotopic analysis of free intracellular metabolites.

Dynamic metabolic models of CHO cell cultures through minimal sets of elementary flux modes

The concept of Elementary Flux Modes (EFMs) has been of central importance in a number of studies involving the analysis of metabolism. In Provost and Bastin (2007) this concept is used to translate the metabolic networks of the different phases of CHO cell cultures into macroscopic bioreactions linking extracellular substrates to products. However, a critical issue concerns the calculation of these elementary flux vectors, as their number combinatorially increases with the size of the metabolic network. In this study, a detailed metabolic network of CHO cells is considered, where the above-mentioned combinatorial explosion makes the computation of the elementary flux modes impossible. To alleviate this problem, a methodology proposed in Jungers et al. (2011) is used to compute a decomposition of admissible flux vectors in a minimal number of elementary flux modes without explicitly enumerating all of them. As a result, a set of macroscopic bioreactions linking the extracellular measured species is obtained at a very low computational expense. The procedure is repeated for the several cell culture phases and a global model is built using a multi-model approach, which is able to successfully predict the evolution of experimental data.

Modeling metabolic networks for mammalian cell systems: general considerations, modeling strategies, and available tools

Over the past decades, the availability of large amounts of information regarding cellular processes and reaction rates, along with increasing knowledge about the complex mechanisms involved in these processes, has changed the way we approach the understanding of cellular processes. We can no longer rely only on our intuition for interpreting experimental data and evaluating new hypotheses, as the information to analyze is becoming increasingly complex. The paradigm for the analysis of cellular systems has shifted from a focus on individual processes to comprehensive global mathematical descriptions that consider the interactions of metabolic, genomic, and signaling networks. Analysis and simulations are used to test our knowledge by refuting or validating new hypotheses regarding a complex system, which can result in predictive capabilities that lead to better experimental design. Different types of models can be used for this purpose, depending on the type and amount of information available for the specific system. Stoichiometric models are based on the metabolic structure of the system and allow explorations of steady state distributions in the network. Detailed kinetic models provide a description of the dynamics of the system, they involve a large number of reactions with varied kinetic characteristics and require a large number of parameters. Models based on statistical information provide a description of the system without information regarding structure and interactions of the networks involved. The development of detailed models for mammalian cell metabolism has only recently started to grow more strongly, due to the intrinsic complexities of mammalian systems, and the limited availability of experimental information and adequate modeling tools. In this work we review the strategies, tools, current advances, and recent models of mammalian cells, focusing mainly on metabolism, but discussing the methodology applied to other types of networks as well.

Optical sensors for monitoring dynamic changes of intracellular metabolite levels in mammalian cells

Knowledge of the in vivo levels, distribution and flux of ions and metabolites is crucial to our understanding of physiology in both healthy and diseased states. The quantitative analysis of the dynamics of ions and metabolites with subcellular resolution in vivo poses a major challenge for the analysis of metabolic processes. Genetically encoded Förster resonance energy transfer (FRET) sensors can be used for real-time in vivo detection of metabolites. FRET sensor proteins, for example, for glucose, can be targeted genetically to any cellular compartment, or even to subdomains (e.g., a membrane surface), by adding signal sequences or fusing the sensors to specific proteins. The sensors can be used for analyses in individual mammalian cells in culture, in tissue slices and in intact organisms. Applications include gene discovery, high-throughput drug screens or systematic analysis of regulatory networks affecting uptake, efflux and metabolism. Quantitative analyses obtained with the help of FRET sensors for glucose or other ions and metabolites provide valuable data for modeling of flux. Here we provide a detailed protocol for monitoring glucose levels in the cytosol of mammalian cell cultures through the use of FRET glucose sensors; moreover, the protocol can be used for other ions and metabolites and for analyses in other organisms, as has been successfully demonstrated in bacteria, yeast and even intact plants. The whole procedure typically takes ∼4 d including seeding and transfection of mammalian cells; the FRET-based analysis of transfected cells takes ∼5 h.