Application of hybrid cybernetic model in simulating myeloma cell culture co-consuming glucose and glutamine with mixed consumption patterns

A dynamic model called hybrid cybernetic model (HCM) based on structured metabolic network is established for simulating mammalian cell metabolism featured with partially substitutable and partially complementary consumption patterns of two substrates, glucose and glutamine. Benefiting from the application of elementary mode analysis (EMA), the complicated metabolic network is decomposed into elementary modes (EMs) facilitating the employment of the hybrid cybernetic framework to investigate the external and internal flux distribution and the regulation mechanism among them. According to different substrate combination, two groups of EMs are obtained, i.e., EMs associated with glucose uptake and simultaneous uptake of glucose and glutamine. Uptake fluxes through various EMs are coupled together via cybernetic variables to maximize substrate uptake. External fluxes and internal fluxes could be calculated and estimated respectively, by the combination of the stoichiometrics of metabolic networks and fluxes through regulated EMs. The model performance is well validated via three sets of experimental data. Through parameter identification of limited number of experimental data, other external metabolites are precisely predicted. The obtained kinetic parameters of three experimental cultures have similar values, which indicates the robustness of the model. Furthermore, the prediction performance of the model is successfully validated based on identified parameters.     

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.     

Exacting predictions by cybernetic model confirmed experimentally: Steady state multiplicity in the chemostat

We demonstrate strong experimental support for the cybernetic model based on maximizing carbon uptake rate in describing the microorganism's regulatory behavior by verifying exacting predictions of steady state multiplicity in a chemostat. Experiments with a feed mixture of glucose and pyruvate show multiple steady state behavior as predicted by the cybernetic model. When multiplicity occurs at a dilution (growth) rate, it results in hysteretic behavior following switches in dilution rate from above and below. This phenomenon is caused by transient paths leading to different steady states through dynamic maximization of the carbon uptake rate. Thus steady state multiplicity is a manifestation of the nonlinearity arising from cybernetic mechanisms rather than of the nonlinear kinetics. The predicted metabolic multiplicity would extend to intracellular states such as enzyme levels and fluxes to be verified in future experiments. © 2012 American Institute of Chemical Engineers Biotechnol. Prog., 2012     

Dynamic modeling of the central carbon metabolism of Escherichia coli

Application of metabolic engineering principles to the rational design of microbial production processes crucially depends on the ability to describe quantitatively the systemic behavior of the central carbon metabolism to redirect carbon fluxes to the product-forming pathways. Despite the importance for several production processes, development of an essential dynamic model for central carbon metabolism of Escherichia coli has been severely hampered by the current lack of kinetic information on the dynamics of the metabolic reactions. Here we present the design and experimental validation of such a dynamic model, which, for the first time, links the sugar transport system (i.e., phosphotransferase system [PTS]) with the reactions of glycolysis and the pentose-phosphate pathway. Experimental observations of intracellular concentrations of metabolites and cometabolites at transient conditions are used to validate the structure of the model and to estimate the kinetic parameters. Further analysis of the detailed characteristics of the system offers the possibility of studying important questions regarding the stability and control of metabolic fluxes.     

Cybernetic modeling and regulation of metabolic pathways in multiple steady states of hybridoma cells

Hybridoma cells utilize a pair of complementary and partially substitutable substrates, glucose and glutamine, for growth. It has been shown that cellular metabolism shifts under different culture conditions. When those cultures at different metabolic states are switched to a continuous mode, they reach different steady states under the same operating conditions. A cybernetic model was constructed to describe the complementary and partial substitutable nature of substrate utilization. The model successfully predicted the metabolic shift and multiple steady-state behavior. The results are consistent with the experimental observation that the history of the culture affects the resulting steady state.     

Integrating cybernetic modeling with pathway analysis provides a dynamic, systems-level description of metabolic control

Cybernetic modeling strives to uncover the inbuilt regulatory programs of biological systems and leverage them toward computational prediction of metabolic dynamics. Because of its focus on incorporating the global aims of metabolism, cybernetic modeling provides a systems-oriented approach for describing regulatory inputs and inferring the impact of regulation within biochemical networks. Combining cybernetic control laws with concepts from metabolic pathway analysis has culminated in a systematic strategy for constructing cybernetic models, which was previously lacking. The newly devised framework relies upon the simultaneous application of local controls that maximize the net flux through each elementary flux mode and global controls that modulate the activities of these modes to optimize the overall nutritional state of the cell. The modeling concepts are illustrated using a simple linear pathway and a larger network representing anaerobic E. coli central metabolism. The E. coli model successfully describes the metabolic shift that occurs upon deleting the pta-ackA operon that is responsible for fermentative acetate production. The model also furnishes predictions that are consistent with experimental results obtained from additional knockout strains as well as strains expressing heterologous genes. Because of the stabilizing influence of the included control variables, the resulting cybernetic models are more robust and reliable than their predecessors in simulating the network response to imposed genetic and environmental perturbations.     

Unveiling steady-state multiplicity in hybridoma cultures: the cybernetic approach

Mammalian cells grown in suspension produce waste metabolites such as lactate, alanine, and ammonia, which reduce the yield of cell mass and the desired product on the nutrients supplied. Previous studies (Cruz et al., 1999; Europa et al., 2000; Follstad et al., 1999) have shown that the cells can be made to alter their metabolism by starving them on their nutrients in continuous cultures at low dilution rates or starting the culture as a fed-batch. This leads to multiple steady states in continuous reactors, with some states being more favorable than others. Mathematical models that take into account the metabolic regulation that leads to these multiple steady states are invaluable tools for bioreactor control. In this article we present a cybernetic modeling strategy in which Metabolic Flux Analysis (MFA) is used to guide the cybernetic formulation. The hybridoma model presented as a result of this strategy considers the partially substitutable, partially complementary nature of glucose and glutamine. The choice of competitions within the network is guided by MFA and the model is successful in explaining the three multiple steady states observed. The cybernetic model though identified for the hybridoma experiments of Hu and others (Europa et al., 2000) seem generally applicable to mammalian systems as it captures the pathways that are common to mammalian cells grown in suspension. The model presented here could be used for start-up strategies for continuous reactors and model-based feedback control for maintaining high productivity of the reactor.     

The principle of flux minimization and its application to estimate stationary fluxes in metabolic networks.

Cellular functions are ultimately linked to metabolic fluxes brought about by thousands of chemical reactions and transport processes. The synthesis of the underlying enzymes and membrane transporters causes the cell a certain 'effort' of energy and external resources. Considering that those cells should have had a selection advantage during natural evolution that enabled them to fulfil vital functions (such as growth, defence against toxic compounds, repair of DNA alterations, etc.) with minimal effort, one may postulate the principle of flux minimization, as follows: given the available external substrates and given a set of functionally important 'target' fluxes required to accomplish a specific pattern of cellular functions, the stationary metabolic fluxes have to become a minimum. To convert this principle into a mathematical method enabling the prediction of stationary metabolic fluxes, the total flux in the network is measured by a weighted linear combination of all individual fluxes whereby the thermodynamic equilibrium constants are used as weighting factors, i.e. the more the thermodynamic equilibrium lies on the right-hand side of the reaction, the larger the weighting factor for the backward reaction. A linear programming technique is applied to minimize the total flux at fixed values of the target fluxes and under the constraint of flux balance (= steady-state conditions) with respect to all metabolites. The theoretical concept is applied to two metabolic schemes: the energy and redox metabolism of erythrocytes, and the central metabolism of Methylobacterium extorquens AM1. The flux rates predicted by the flux-minimization method exhibit significant correlations with flux rates obtained by either kinetic modelling or direct experimental determination. Larger deviations occur for segments of the network composed of redundant branches where the flux-minimization method always attributes the total flux to the thermodynamically most favourable branch. Nevertheless, compared with existing methods of structural modelling, the principle of flux minimization appears to be a promising theoretical approach to assess stationary flux rates in metabolic systems in cases where a detailed kinetic model is not yet available.     

Quantitative prediction of cellular metabolism with constraint-based models: the COBRA Toolbox

The manner in which microorganisms utilize their metabolic processes can be predicted using constraint-based analysis of genome-scale metabolic networks. Herein, we present the constraint-based reconstruction and analysis toolbox, a software package running in the Matlab environment, which allows for quantitative prediction of cellular behavior using a constraint-based approach. Specifically, this software allows predictive computations of both steady-state and dynamic optimal growth behavior, the effects of gene deletions, comprehensive robustness analyses, sampling the range of possible cellular metabolic states and the determination of network modules. Functions enabling these calculations are included in the toolbox, allowing a user to input a genome-scale metabolic model distributed in Systems Biology Markup Language format and perform these calculations with just a few lines of code. The results are predictions of cellular behavior that have been verified as accurate in a growing body of research. After software installation, calculation time is minimal, allowing the user to focus on the interpretation of the computational results.     

Metabolic design based on a coupled gene expression-metabolic network model of tryptophan production in Escherichia coli

The presumably high potential of a holistic design approach for complex biochemical reaction networks is exemplified here for the network of tryptophan biosynthesis from glucose, a system whose components have been investigated thoroughly before. A dynamic model that combines the behavior of the trp operon gene expression with the metabolic network of central carbon metabolism and tryptophan biosynthesis is investigated. This model is analyzed in terms of metabolic fluxes, metabolic control, and nonlinear optimization. We compare two models for a wild-type strain and another model for a tryptophan producer. An integrated optimization of the whole network leads to a significant increase in tryptophan production rate for all systems under study. This enhancement is well above the increase that can be achieved by an optimization of subsystems. A constant ratio of control coefficients on tryptophan synthesis rate has been identified for the models regarding or disregarding trp operon expression. Although we found some examples where flux control coefficients even contradict the trends of enzyme activity changes in an optimized profile, flux control can be used as an indication for enzymes that have to be taken into account in optimization.     

Computational approaches to the topology, stability and dynamics of metabolic networks

Cellular metabolism is characterized by an intricate network of interactions between biochemical fluxes, metabolic compounds and regulatory interactions. To investigate and eventually understand the emergent global behavior arising from such networks of interaction is not possible by intuitive reasoning alone. This contribution seeks to describe recent computational approaches that aim to asses the topological and functional properties of metabolic networks. In particular, based on a recently proposed method, it is shown that it is possible to acquire a quantitative picture of the possible dynamics of metabolic systems, without assuming detailed knowledge of the underlying enzyme-kinetic rate equations and parameters. Rather, the method builds upon a statistical exploration of the comprehensive parameter space to evaluate the dynamic capabilities of a metabolic system, thus providing a first step towards the transition from topology to function of metabolic pathways. Utilizing this approach, the role of feedback mechanisms in the maintenance of stability is discussed using minimal models of cellular pathways.     

Systematic Construction of Kinetic Models from Genome-Scale Metabolic Networks

The quantitative effects of environmental and genetic perturbations on metabolism can be studied in silico using kinetic models. We present a strategy for large-scale model construction based on a logical layering of data such as reaction fluxes, metabolite concentrations, and kinetic constants. The resulting models contain realistic standard rate laws and plausible parameters, adhere to the laws of thermodynamics, and reproduce a predefined steady state. These features have not been simultaneously achieved by previous workflows. We demonstrate the advantages and limitations of the workflow by translating the yeast consensus metabolic network into a kinetic model. Despite crudely selected data, the model shows realistic control behaviour, a stable dynamic, and realistic response to perturbations in extracellular glucose concentrations. The paper concludes by outlining how new data can continuously be fed into the workflow and how iterative model building can assist in directing experiments.     

Elementary mode analysis: a useful metabolic pathway analysis tool for characterizing cellular metabolism

Elementary mode analysis is a useful metabolic pathway analysis tool to identify the structure of a metabolic network that links the cellular phenotype to the corresponding genotype. The analysis can decompose the intricate metabolic network comprised of highly interconnected reactions into uniquely organized pathways. These pathways consisting of a minimal set of enzymes that can support steady state operation of cellular metabolism represent independent cellular physiological states. Such pathway definition provides a rigorous basis to systematically characterize cellular phenotypes, metabolic network regulation, robustness, and fragility that facilitate understanding of cell physiology and implementation of metabolic engineering strategies. This mini-review aims to overview the development and application of elementary mode analysis as a metabolic pathway analysis tool in studying cell physiology and as a basis of metabolic engineering.     

定量分析丝氨酸/甘氨酸转换影响胃癌细胞增殖的动力学机制

目的 胃癌（GC）严重影响人类的健康生活,研究表明其与丝氨酸/甘氨酸代谢密切相关。丝氨酸/甘氨酸代谢对于肿瘤细胞的增殖能力具有重要影响。本文的研究目的是探究丝氨酸/甘氨酸代谢能够影响胃癌细胞增殖能力的分子机制。方法 本文通过一种基于随机和非梯度系统的势能景观所建立的大型代谢网络动力学建模方法,构建了一个稳定的胃癌细胞代谢动力学模型。基于对模型的调控,定量分析丝氨酸/甘氨酸代谢影响胃癌细胞增殖的动力学机制。对一般代谢网络动力学方程添加随机噪声,通过随机动力学分解得到代谢网络参数空间的李雅普诺夫（Lyapunov）函数。进一步减少与随机波动相关的Lyapunov函数变化,从而得到稳定的代谢网络。结果 在动力学参数不足的情况下,成功构建了胃癌细胞代谢网络的动力学模型。当胞外丝氨酸可用时,模型优先消耗丝氨酸;当甘氨酸生成丝氨酸的速率增加时,模型显著上调生成S-腺苷甲硫氨酸（SAM）和S-腺苷同型半胱氨酸（SAH）的稳态通量。结论 本文证明了胃癌细胞对于丝氨酸的优先摄取以及丝氨酸/甘氨酸转化速率对SAM生成的重要作用,其可能通过调节细胞甲基化进程影响胃癌细胞的增殖能力,为靶向丝氨酸/甘氨酸代谢的癌症治疗提供了新的思路和方向。     

Something from nothing: bridging the gap between constraint-based and kinetic modelling

Two divergent modelling methodologies have been adopted to increase our understanding of metabolism and its regulation. Constraint-based modelling highlights the optimal path through a stoichiometric network within certain physicochemical constraints. Such an approach requires minimal biological data to make quantitative inferences about network behaviour; however, constraint-based modelling is unable to give an insight into cellular substrate concentrations. In contrast, kinetic modelling aims to characterize fully the mechanics of each enzymatic reaction. This approach suffers because parameterizing mechanistic models is both costly and time-consuming. In this paper, we outline a method for developing a kinetic model for a metabolic network, based solely on the knowledge of reaction stoichiometries. Fluxes through the system, estimated by flux balance analysis, are allowed to vary dynamically according to linlog kinetics. Elasticities are estimated from stoichiometric considerations. When compared to a popular branched model of yeast glycolysis, we observe an excellent agreement between the real and approximate models, despite the absence of (and indeed the requirement for) experimental data for kinetic constants. Moreover, using this particular methodology affords us analytical forms for steady state determination, stability analyses and studies of dynamical behaviour.     

Dynamic metabolic modeling for a MAB bioprocess

Production of monoclonal antibodies (MAb) for diagnostic or therapeutic applications has become an important task in the pharmaceutical industry. The efficiency of high-density reactor systems can be potentially increased by model-based design and control strategies. Therefore, a reliable kinetic model for cell metabolism is required. A systematic procedure based on metabolic modeling is used to model nutrient uptake and key product formation in a MAb bioprocess during both the growth and post-growth phases. The approach combines the key advantages of stoichiometric and kinetic models into a complete metabolic network while integrating the regulation and control of cellular activity. This modeling procedure can be easily applied to any cell line during both the cell growth and post-growth phases. Quadratic programming (QP) has been identified as a suitable method to solve the underdetermined constrained problem related to model parameter identification. The approach is illustrated for the case of murine hybridoma cells cultivated in stirred spinners.     

Prediction of dynamic behavior of mutant strains from limited wild-type data

Metabolic engineering is the field of introducing genetic changes in organisms so as to modify their function towards synthesizing new products of high impact to society. However, engineered cells frequently have impaired growth rates thus seriously limiting the rate at which such products are made. The problem is attributable to inadequate understanding of how a metabolic network functions in a dynamic sense. Predictions of mutant strain behavior in the past have been based on steady state theories such as flux balance analysis (FBA), minimization of metabolic adjustment (MOMA), and regulatory on/off minimization (ROOM). Such predictions are restricted to product yields and cannot address productivity, which is of focal interest to applications. We demonstrate that our framework ( [Song and Ramkrishna, 2010] and [Song and Ramkrishna, 2011]), based on a “cybernetic” view of metabolic systems, makes predictions of the dynamic behavior of mutant strains of Escherichia coli from a limited amount of data obtained from the wild-type. Dynamic frameworks must necessarily address the issue of metabolic regulation, which the cybernetic approach does by postulating that metabolism is an optimal dynamic response of the organism to the environment in driving reactions towards ensuring survival. The predictions made in this paper are without parallel in the literature and lay the foundation for rational metabolic engineering.     

Model of Tryptophan Metabolism,Readily Scalable Using Tissue-specific Gene Expression Data

Tryptophan is utilized in various metabolic routes including protein synthesis, serotonin, and melatonin synthesis and the kynurenine pathway. Perturbations in these pathways have been associated with neurodegenerative diseases and cancer. Here we present a comprehensive kinetic model of the complex network of human tryptophan metabolism based upon existing kinetic data for all enzymatic conversions and transporters. By integrating tissue-specific expression data, modeling tryptophan metabolism in liver and brain returned intermediate metabolite concentrations in the physiological range. Sensitivity and metabolic control analyses identified expected key enzymes to govern fluxes in the branches of the network. Combining tissue-specific models revealed a considerable impact of the kynurenine pathway in liver on the concentrations of neuroactive derivatives in the brain. Moreover, using expression data from a cancer study predicted metabolite changes that resembled the experimental observations. We conclude that the combination of the kinetic model with expression data represents a powerful diagnostic tool to predict alterations in tryptophan metabolism. The model is readily scalable to include more tissues, thereby enabling assessment of organismal tryptophan metabolism in health and disease.     

Development of a mathematical model for evaluating the dynamics of normal and apoptotic Chinese hamster ovary cells

A metabolic flux based methodology was developed for modeling the metabolism of a Chinese hamster ovary cell line. The elimination of insignificant fluxes resulted in a simplified metabolic network which was the basis for modeling the significant metabolites. Employing kinetic rate expressions for growing and non-growing subpopulations, a logistic model was developed for cell growth and dynamic models were formulated to describe culture composition and monoclonal antibody (MAb) secretion. The model was validated for a range of nutrient concentrations. Good agreement was obtained between model predictions and experimental data. The ultimate goal of this study is to establish a comprehensive dynamic model which may be used for model-based optimization of the cell culture for MAb production in both batch and fed-batch systems.