A flexible state-space approach for the modeling of metabolic networks II: advanced interrogation of hybridoma metabolism

Having previously introduced the mathematical framework of topological metabolic analysis (TMA) - a novel optimization-based technique for modeling metabolic networks of arbitrary size and complexity - we demonstrate how TMA facilitates unique methods of metabolic interrogation. With the aid of several hybridoma metabolic investigations as case-studies (Bonarius et al., 1995, 1996, 2001), we first establish that the TMA framework identifies biologically important aspects of the metabolic network under investigation. We also show that the use of a structured weighting approach within our objective provides a substantial modeling benefit over an unstructured, uniform, weighting approach. We then illustrate the strength of TAM as an advanced interrogation technique, first by using TMA to prove the existence of (and to quantitatively describe) multiple topologically distinct configurations of a metabolic network that each optimally model a given set of experimental observations. We further show that such alternate topologies are indistinguishable using existing stoichiometric modeling techniques, and we explain the biological significance of the topological variables appearing within our model. By leveraging the manner in which TMA implements metabolite inputs and outputs, we also show that metabolites whose possible metabolic fates are inadequately described by a given network reconstruction can be quickly identified. Lastly, we show how the use of the TMA aggregate objective function (AOF) permits the identification of modeling solutions that can simultaneously consider experimental observations, underlying biological motivations, or even purely engineering- or design-based goals.     

Tendency modeling: a new approach to obtain simplified kinetic models of metabolism applied to Saccharomyces cerevisiae

A novel approach to construct kinetic models of metabolic pathways, to be used in metabolic engineering, is presented: the tendency modeling approach. This approach greatly facilitates the construction of these models and can easily be applied to complex metabolic networks. The resulting models contain a minimal number of parameters; identification of their values is straightforward. Use of in vitro obtained information in the identification of the kinetic equations is minimized. The tendency modeling approach has been used to derive a dynamic model of primary metabolism for aerobic growth of Saccharomyces cerevisiae on glucose, in which compartmentation is included. Simulation results obtained with the derived model are satisfying for most of the carbon metabolites that have been measured. Compared to a more detailed model, the simulations of our model are less accurate, but taking into account the much smaller number of kinetic parameters (35 instead of 84), the tendency the modeling approach is considered promising.     

Parameter estimation in models combining signal transduction and metabolic pathways: the dependent input approach

Biological complexity and limited quantitative measurements pose severe challenges to standard engineering methodologies for modelling and simulation of genes and gene products integrated in a functional network. In particular, parameter quantification is a bottleneck, and therefore parameter estimation, identifiability, and optimal experiment design are important research topics in systems biology. An approach is presented in which unmodelled dynamics are replaced by fictitious 'dependent inputs'. The dependent input approach is particularly useful in validation experiments, because it allows one to fit model parameters to experimental data generated by a reference cell type ('wild-type') and then test this model on data generated by a variation ('mutant'), so long as the mutations only affect the unmodelled dynamics that produce the dependent inputs. Another novel feature of the approach is in the inclusion of a priori information in a multi-objective identification criterion, making it possible to obtain estimates of parameter values and their variances from a relatively limited experimental data set. The pathways that control the nitrogen uptake fluxes in baker's yeast (Saccharomyces cerevisiae) have been studied. Well-defined perturbation experiments were performed on cells growing in steady-state. Time-series data of extracellular and intracellular metabolites were obtained, as well as mRNA levels. A nonlinear model was proposed and was shown to be structurally identifiable given data of its inputs and outputs. The identified model is a reliable representation of the metabolic system, as it could correctly describe the responses of mutant cells and different perturbations.     

A Green's function method for analysis of oxygen delivery to tissue by microvascular networks

A theoretical model is formulated for analyzing oxygen delivery from an arbitrary network configuration of cylindrical microvessels to a finite region of tissue. In contrast to models based on the classical Krogh cylinder approach, this model requires no a priori assumptions concerning the extent of the tissue region supplied with oxygen by each vessel segment. Steady-state conditions are assumed, and oxygen consumption in the tissue is assumed to be uniform. The nonlinear dissociation characteristics of oxyhemoglobin are taken into account. A computationally efficient Green's function approach is used, in which the tissue oxygen field is expressed in terms of the distribution of source strengths along each segment. The utility of the model is illustrated by analyses of oxygen delivery to a cuboidal tissue region by a single segment and by a six-segment network. It is found that the fractional contribution of the proximal segments to total oxygen delivery increases with decreasing flow rate and metabolic rate.     

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.     

Analysis and implementation of a neural extended Kalman filter for target tracking

Having a better motion model in the state estimator is one way to improve target tracking performance. Since the motion model of the target is not known a priori, either robust modeling techniques or adaptive modeling techniques are required. The neural extended Kalman filter is a technique that learns unmodeled dynamics while performing state estimation in the feedback loop of a control system. This coupled system performs the standard estimation of the states of the plant while estimating a function to approximate the difference between the given state-coupling function model and the behavior of the true plant dynamics. At each sample step, this new model is added to the existing model to improve the state estimate. The neural extended Kalman filter has also been investigated as a target tracking estimation routine. Implementation issues for this adaptive modeling technique, including neural network training parameters, were investigated and an analysis was made of the quality of performance that the technique can have for tracking maneuvering targets.     

Phenotype-centric modeling for rational metabolic engineering

Phenotype-centric modeling enables a paradigm shift in the analysis of mechanistic models. It brings the focus to a network's biochemical phenotypes and their relationship with measurable traits (e.g., product yields, system dynamics, signal amplification factors, etc.) and away from computationally intensive simulation-centric modeling. Here, we explore applications of this new modeling strategy in the field of rational metabolic engineering using the amorphadiene biosynthetic network as a case study. This network has previously been studied using a mechanistic model and the simulation-centric strategy, and thus provides an excellent means to compare and contrast results obtained from these two very different strategies. We show that the phenotype-centric strategy, without values for the parameters, not only identifies beneficial intervention strategies obtained with the simulation-centric strategy, but it also provides an understanding of the mechanistic context for the validity of these predictions. Additionally, we propose a set of hypothetical strains with the potential to outperform reported production strains and to enhance the mechanistic understanding of the amorphadiene biosynthetic network. Further, we identify the landscape of possible intervention strategies for the given model. We believe that phenotype-centric modeling can advance the field of rational metabolic engineering by enabling the development of next generation kinetics-based algorithms and methods that do not rely on a priori knowledge of kinetic parameters but allow a structured, global analysis of system design in the parameter space.     

A phenomenological model for atherosclerotic plaque growth and rupture

The objective of this communication is to develop a computer-based framework for the overall coupled phenomena leading to growth and rupture of atherosclerotic plaques. The modeling is purposely simplified to expose the dominant phenomenological controlling mechanisms, and their coupled interaction. The main ingredients of the present simplified modeling approach, describing the events that occur due to the presence and oxidation of excess low-density lipoprotein (LDL) in the intima, are: (i) adhesion of monocytes to the endothelial surface, which is controlled by the intensity of the blood flow and the adhesion molecules stimulated by the excess LDL, (ii) penetration of the monocytes into the intima and subsequent inflammation of the tissue, and (iii) rupture of the plaque accompanied with some degree of thrombus formation or even subsequent occlusive thrombosis. The set of resulting coupled equations, each modeling entirely different physical events, is solved using an iterative staggering scheme, which allows the equations to be solved in a computationally convenient decoupled fashion. Theoretical convergence properties of the scheme are given as a function of physical parameters involved. A numerical example is given to illustrate the modeling approach and an a priori prediction for time to rupture as a function of arterial geometry, diameter of the monocyte, adhesion stress, bulk modulus of the ruptured wall material, blood viscosity, flow rate and mass density of the monocytes.     

基于约束的建模方法在代谢网络中的应用研究

代谢网络在各种细胞功能和生命过程中发挥着至关重要的作用。随着细胞网络重建工程的迅速发展,可用的基因组水平代谢网络越来越多,因而计算方法在这些网络的结构功能分析中越来越重要。基于约束的建模方法不像图论方法那样仅考虑代谢模型的纯拓扑结构,也不像各种动力学建模方法那样需求详尽的热力学参数,因而极具优势。采用基于约束的建模方法对一个含619个基因,655个代谢物和743个代谢反应的金黄色葡萄球菌(Staphylococcusaureus)代谢网络进行了分析,主要研究了该模型的网络结构特征,以及其最优生长率、动态生长情况和基因删除学习等。本研究提供了一个对金黄色葡萄球菌代谢网络进行约束建模分析的初步框架。     

On‐line adaptive metabolic flux analysis: Application to PHB production by mixed microbial cultures

In this work, an algorithm for on‐line adaptive metabolic flux analysis (MFA) is proposed and applied to polyhydroxybutyrate (PHB) production by mixed microbial cultures (MMC). In this process, population dynamics constitutes an important source of perturbation to MFA calculations because some stoichiometric and energetic parameters of the underlying metabolic network are continuously changing over time. The proposed algorithm is based on the application of the observer‐based estimator (OBE) to the central MFA equation, whereby the role of the OBE is to force the accumulation of intracellular metabolites to converge to zero by adjusting the values of unknown network parameters. The algorithm was implemented in a reactor equipped with on‐line analyses of dissolved oxygen and carbon dioxide through respirometric and titrimetric measurements. The oxygen and carbon dioxide fluxes were measured directly, whereas acetate, PHB, and sludge production fluxes were estimated indirectly using a projection of latent structures model calibrated a priori with off‐line measurements. The algorithm was implemented in a way that the network parameters associated with biosynthesis were adjusted on‐line. The algorithm proofed to converge exponentially with the steady state error always below 1 mmol/L. The estimated fluxes passed the consistency index test for experimental error variances as low as 1%. The comparison of measured and estimated respiratory coefficient and of the theoretical and estimated yield of sludge on acetate further confirmed the metabolic consistency of the parameters that were estimated on‐line. © 2009 American Institute of Chemical Engineers Biotechnol. Prog., 2009     

Dynamics of muscle glycogenolysis modeled with pH time course computation and pH-dependent reaction equilibria and enzyme kinetics

Cellular metabolites are moieties defined by their specific binding constants to H+, Mg2+, and K+ or anions without ligands. As a consequence, every biochemical reaction in the cytoplasm has an associated proton stoichiometry that is generally noninteger- and pH-dependent. Therefore, with metabolic flux, pH is altered in a medium with finite buffer capacity. Apparent equilibrium constants and maximum enzyme velocities, which are functions of pH, are also altered. We augmented an earlier mathematical model of skeletal muscle glycogenolysis with pH-dependent enzyme kinetics and reaction equilibria to compute the time course of pH changes. Analysis shows that kinetics and final equilibrium states of the closed system are highly constrained by the pH-dependent parameters. This kinetic model of glycogenolysis, coupled to creatine kinase and adenylate kinase, simulated published experiments made with a cell-free enzyme mixture to reconstitute the network and to synthesize PCr and lactate in vitro. Using the enzyme kinetic and thermodynamic data in the literature, the simulations required minimal adjustments of parameters to describe the data. These results show that incorporation of appropriate physical chemistry of the reactions with accurate kinetic modeling gives a reasonable simulation of experimental data and is necessary for a physically correct representation of the metabolic network. The approach is general for modeling metabolic networks beyond the specific pathway and conditions presented here.     

An integrated modeling-experimental strategy for the analysis of metabolic pathways

An inherent problem in studying the behavior of a metabolic pathway is the impossibility of developing a complete, detailed model that includes all the cellular processes that have an impact on the set of fluxes in such a pathway. Lacking this, one requires some means of modeling the interactions between a metabolic pathway and other cellular processes for the purpose of analyzing pathway characteristics within the cell (e.g., determining sensitivity coefficients for various steps in the pathway) with a minimal amount of time and effort. A general framework is developed for studying these issues in a rigorous manner. Using this framework, detailed knowledge about a metabolic pathway (i.e., a set of rate expressions for steps in the pathway) can be combined with the results from a relatively simple set of experiments in order to obtain estimates for the sensitivity of the pathway to enzyme activities, inhibition constants, and other parameters that determine the pathway's behavior, while accounting for the pathway's interaction with the rest of the cellular metabolism. A model system representing amino acid production is used to illustrate the problem and to provide results based on computational experiments. The modeling strategy described here should be useful in genetic design to improve pathway fluxes and metabolic network selectivity.     

Indistinguishability and identifiability analysis of linear compartmental models.

Two compartmental model structures are said to be indistinguishable if they have the same input-output properties. In cases in which available a priori information is not sufficient to specify a unique compartmental model structure, indistinguishable model structures may have to be generated and their attributes examined for relevance. An algorithm is developed that, for a given compartmental model, investigates the complete set of models with the same number of compartments and the same input-output structure as the original model, applies geometrical rules necessary for indistinguishable models, and test models meeting the geometrical criteria for equality of transfer functions. Identifiability is also checked in the algorithm. The software consists of three programs. Program 1 determines the number of locally identifiable parameters. Program 2 applies several geometrical rules that eliminate many (generally most) of the candidate models. Program 3 checks the equality between system transfer functions of the original model and models being tested. Ranks of Jacobian matrices and submatrices and other criteria are used to check patterns of moment invariants and local identifiability. Structural controllability and structural observability are checked throughout the programs. The approach was successfully used to corroborate results from examples investigated by others.     

Hybrid dynamic modeling of Escherichia coli central metabolic network combining Michaelis–Menten and approximate kinetic equations

The construction of dynamic metabolic models at reaction network level requires the use of mechanistic enzymatic rate equations that comprise a large number of parameters. The lack of knowledge on these equations and the difficulty in the experimental identification of their associated parameters, represent nowadays the limiting factor in the construction of such models. In this study, we compare four alternative modeling approaches based on Michaelis–Menten kinetics for the bi-molecular reactions and different types of simplified rate equations for the remaining reactions (generalized mass action, convenience kinetics, lin-log and power-law). Using the mechanistic model for Escherichia coli central carbon metabolism as a benchmark, we investigate the alternative modeling approaches through comparative simulations analyses. The good dynamic behavior and the powerful predictive capabilities obtained using the hybrid model composed of Michaelis–Menten and the approximate lin-log kinetics indicate that this is a possible suitable approach to model complex large-scale networks where the exact rate laws are unknown.     

Validation-based model selection for 13C metabolic flux analysis with uncertain measurement errors

Accurate measurements of metabolic fluxes in living cells are central to metabolism research and metabolic engineering. The gold standard method is model-based metabolic flux analysis (MFA), where fluxes are estimated indirectly from mass isotopomer data with the use of a mathematical model of the metabolic network. A critical step in MFA is model selection: choosing what compartments, metabolites, and reactions to include in the metabolic network model. Model selection is often done informally during the modelling process, based on the same data that is used for model fitting (estimation data). This can lead to either overly complex models (overfitting) or too simple ones (underfitting), in both cases resulting in poor flux estimates. Here, we propose a method for model selection based on independent validation data. We demonstrate in simulation studies that this method consistently chooses the correct model in a way that is independent on errors in measurement uncertainty. This independence is beneficial, since estimating the true magnitude of these errors can be difficult. In contrast, commonly used model selection methods based on the χ²-test choose different model structures depending on the believed measurement uncertainty; this can lead to errors in flux estimates, especially when the magnitude of the error is substantially off. We present a new approach for quantification of prediction uncertainty of mass isotopomer distributions in other labelling experiments, to check for problems with too much or too little novelty in the validation data. Finally, in an isotope tracing study on human mammary epithelial cells, the validation-based model selection method identified pyruvate carboxylase as a key model component. Our results argue that validation-based model selection should be an integral part of MFA model development.     

Network modeling of dynamic brain interactions predicts emergence of neural information that supports human cognitive behavior

How cognitive task behavior is generated by brain network interactions is a central question in neuroscience. Answering this question calls for the development of novel analysis tools that can firstly capture neural signatures of task information with high spatial and temporal precision (the “where and when”) and then allow for empirical testing of alternative network models of brain function that link information to behavior (the “how”). We outline a novel network modeling approach suited to this purpose that is applied to noninvasive functional neuroimaging data in humans. We first dynamically decoded the spatiotemporal signatures of task information in the human brain by combining MRI-individualized source electroencephalography (EEG) with multivariate pattern analysis (MVPA). A newly developed network modeling approach—dynamic activity flow modeling—then simulated the flow of task-evoked activity over more causally interpretable (relative to standard functional connectivity [FC] approaches) resting-state functional connections (dynamic, lagged, direct, and directional). We demonstrate the utility of this modeling approach by applying it to elucidate network processes underlying sensory–motor information flow in the brain, revealing accurate predictions of empirical response information dynamics underlying behavior. Extending the model toward simulating network lesions suggested a role for the cognitive control networks (CCNs) as primary drivers of response information flow, transitioning from early dorsal attention network-dominated sensory-to-response transformation to later collaborative CCN engagement during response selection. These results demonstrate the utility of the dynamic activity flow modeling approach in identifying the generative network processes underlying neurocognitive phenomena.

How is cognitive task behavior generated by brain network interactions? This study describes a novel network modeling approach and applies it to source electroencephalography data. The model accurately predicts future information dynamics underlying behavior and (via simulated lesioning) suggests a role for cognitive control networks as key drivers of response information flow.     

Cybernetic modeling based on pathway analysis for Penicillium chrysogenum fed-batch fermentation

A macrokinetic model employing cybernetic methodology is proposed to describe mycelium growth and penicillin production. Based on the primordial and complete metabolic network of Penicillium chrysogenum found in the literature, the modeling procedure is guided by metabolic flux analysis and cybernetic modeling framework. The abstracted cybernetic model describes the transients of the consumption rates of the substrates, the assimilation rates of intermediates, the biomass growth rate, as well as the penicillin formation rate. Combined with the bioreactor model, these reaction rates are linked with the most important state variables, i.e., mycelium, substrate and product concentrations. Simplex method is used to estimate the sensitive parameters of the model. Finally, validation of the model is carried out with 20 batches of industrial-scale penicillin cultivation.