Structural kinetic modeling of metabolic networks

The stability and robustness of metabolic states: identifying stabilizing sites in metabolic networks

The dynamic behavior of metabolic networks is governed by numerous regulatory mechanisms, such as reversible phosphorylation, binding of allosteric effectors or temporal gene expression, by which the activity of the participating enzymes can be adjusted to the functional requirements of the cell. For most of the cellular enzymes, such regulatory mechanisms are at best qualitatively known, whereas detailed enzyme-kinetic models are lacking. To explore the possible dynamic behavior of metabolic networks in cases of lacking or incomplete enzyme-kinetic information, we present a computational approach based on structural kinetic modeling. We derive statistical measures for the relative impact of enzyme-kinetic parameters on dynamic properties (such as local stability) and apply our approach to the metabolism of human erythrocytes. Our findings show that allosteric enzyme regulation significantly enhances the stability of the network and extends its potential dynamic behavior. Moreover, our approach allows to differentiate quantitatively between metabolic states related to senescence and metabolic collapse of the human erythrocyte. We think that the proposed method represents an important intermediate step on the long way from topological network analysis to detailed kinetic modeling of complex metabolic networks.

Bistability and oscillations in the Huang-Ferrell model of MAPK signaling

Physicochemical models of signaling pathways are characterized by high levels of structural and parametric uncertainty, reflecting both incomplete knowledge about signal transduction and the intrinsic variability of cellular processes. As a result, these models try to predict the dynamics of systems with tens or even hundreds of free parameters. At this level of uncertainty, model analysis should emphasize statistics of systems-level properties, rather than the detailed structure of solutions or boundaries separating different dynamic regimes. Based on the combination of random parameter search and continuation algorithms, we developed a methodology for the statistical analysis of mechanistic signaling models. In applying it to the well-studied MAPK cascade model, we discovered a large region of oscillations and explained their emergence from single-stage bistability. The surprising abundance of strongly nonlinear (oscillatory and bistable) input/output maps revealed by our analysis may be one of the reasons why the MAPK cascade in vivo is embedded in more complex regulatory structures. We argue that this type of analysis should accompany nonlinear multiparameter studies of stationary as well as transient features in network dynamics.

Molecular studies of cell communication systems lead to models with multiple free parameters. Analysis of dynamical behavior of these models presents considerable challenge. We have developed a computational approach for the efficient exploration of dynamic behavior in such models and applied this method to the model of the Mitogen Activated Protein Kinase cascade, a signaling network conserved in all eukaryotes. Previous analysis of this model suggested that it works as a reversible switch. We have shown that it can also function as an irreversible switch and as a clock.

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.

Experimental and mathematical approaches to modeling plant metabolic networks

To support their sessile and autotrophic lifestyle higher plants have evolved elaborate networks of metabolic pathways. Dynamic changes in these metabolic networks are among the developmental forces underlying the functional differentiation of organs, tissues and specialized cell types. They are also important in the various interactions of a plant with its environment. Further complexity is added by the extensive compartmentation of the various interconnected metabolic pathways in plants. Thus, although being used widely for assessing the control of metabolic flux in microbes, mathematical modeling approaches that require steady-state approximations are of limited utility for understanding complex plant metabolic networks. However, considerable progress has been made when manageable metabolic subsystems were studied. In this article, we will explain in general terms and using simple examples the concepts underlying stoichiometric modeling (metabolic flux analysis and metabolic pathway analysis) and kinetic approaches to modeling (including metabolic control analysis as a special case). Selected studies demonstrating the prospects of these approaches, or combinations of them, for understanding the control of flux through particular plant pathways are discussed. We argue that iterative cycles of (dry) mathematical modeling and (wet) laboratory testing will become increasingly important for simulating the distribution of flux in plant metabolic networks and deriving rational experimental designs for metabolic engineering efforts.

Analysis of structures causing instabilities

We present a simple new method to systematically identify all topological structures (e.g., positive feedback loops) potentially leading to locally unstable steady states: ICSA-The instability causing structure analysis. Systems without any instability causing structure (i.e., not fulfilling the necessary topological condition for instabilities) cannot have unstable steady states. It follows that common bistability or multistability and Hopf bifurcations are excluded and sustained oscillations and deterministic chaos are most unlikely. The ICSA leads to new insights into the topological organization of chemical and biochemical systems, such as metabolic, gene regulatory, and signal transduction networks.

From structure to dynamics of metabolic pathways: application to the plant mitochondrial TCA cycle

MOTIVATION

Mitochondrial metabolism, dominated by the reactions of the tricarboxylic acid (TCA) cycle, is of vital importance for a wide range of metabolic processes. In particular for autotrophic tissue, such as plant leaves, the TCA cycle marks the point of divergence of anabolic pathways and plays an essential role in biosynthesis. However, despite extensive knowledge about its stoichiometric properties, the function and the dynamical capabilities of the TCA cycle remain largely unknown.

METHODS AND RESULTS

Based on a recently proposed formalism, we investigate the dynamic and functional properties of the mitochondrial TCA cycle of plants. Starting with the structural properties, as described by the elementary flux modes of the system, we aim for the transition from structure to the dynamics of the TCA cycle. Using a parametric description of the system, encompassing all possible differential equations and parameter values, we detect and quantify regimes of different dynamic behavior. Optimizing the system with respect to dynamic stability, we demonstrate that maximal stability is associated with specific (relative) metabolite concentrations and flux values that are subsequently compared to the experimental literature. Our analysis also serves as a general example how to elucidate the transition from the structure to the dynamics of metabolic pathways.

Bringing metabolic networks to life: integration of kinetic, metabolic, and proteomic data

Background

Translating a known metabolic network into a dynamic model requires reasonable guesses of all enzyme parameters. In Bayesian parameter estimation, model parameters are described by a posterior probability distribution, which scores the potential parameter sets, showing how well each of them agrees with the data and with the prior assumptions made.

Results

We compute posterior distributions of kinetic parameters within a Bayesian framework, based on integration of kinetic, thermodynamic, metabolic, and proteomic data. The structure of the metabolic system (i.e., stoichiometries and enzyme regulation) needs to be known, and the reactions are modelled by convenience kinetics with thermodynamically independent parameters. The parameter posterior is computed in two separate steps: a first posterior summarises the available data on enzyme kinetic parameters; an improved second posterior is obtained by integrating metabolic fluxes, concentrations, and enzyme concentrations for one or more steady states. The data can be heterogenous, incomplete, and uncertain, and the posterior is approximated by a multivariate log-normal distribution. We apply the method to a model of the threonine synthesis pathway: the integration of metabolic data has little effect on the marginal posterior distributions of individual model parameters. Nevertheless, it leads to strong correlations between the parameters in the joint posterior distribution, which greatly improve the model predictions by the following Monte-Carlo simulations.

Conclusion

We present a standardised method to translate metabolic networks into dynamic models. To determine the model parameters, evidence from various experimental data is combined and weighted using Bayesian parameter estimation. The resulting posterior parameter distribution describes a statistical ensemble of parameter sets; the parameter variances and correlations can account for missing knowledge, measurement uncertainties, or biological variability. The posterior distribution can be used to sample model instances and to obtain probabilistic statements about the model's dynamic behaviour.