PubMed Health⌕ Search

Biomedical subjects

Edda Klipp

Publications and source records attributed to Edda Klipp.

16 recordsLinked to original sources
PubMed2006

Bringing metabolic networks to life: convenience rate law and thermodynamic constraints.

BACKGROUND: Translating a known metabolic network into a dynamic model requires rate laws for all chemical reactions. The mathematical expressions depend on the underlying enzymatic mechanism; they can become quite involved and may contain a large number of parameters. Rate laws and enzyme parameters are still unknown for most enzymes. RESULTS: We introduce a simple and general rate law called "convenience kinetics". It can be derived from a simple random-order enzyme mechanism. Thermodynamic laws can impose dependencies on the kinetic parameters. Hence, to facilitate model fitting and parameter optimisation for large networks, we introduce thermodynamically independent system parameters: their values can be varied independently, without violating thermodynamical constraints. We achieve this by expressing the equilibrium constants either by Gibbs free energies of formation or by a set of independent equilibrium constants. The remaining system parameters are mean turnover rates, generalised Michaelis-Menten constants, and constants for inhibition and activation. All parameters correspond to molecular energies, for instance, binding energies between reactants and enzyme. CONCLUSION: Convenience kinetics can be used to translate a biochemical network--manually or automatically--into a dynamical model with plausible biological properties. It implements enzyme saturation and regulation by activators and inhibitors, covers all possible reaction stoichiometries, and can be specified by a small number of parameters. Its mathematical form makes it especially suitable for parameter estimation and optimisation. Parameter estimates can be easily computed from a least-squares fit to Michaelis-Menten values, turnover rates, equilibrium constants, and other quantities that are routinely measured in enzyme assays and stored in kinetic databases.

Energy Metabolism↗
PubMed2006

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 heterogeneous, 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.

Bayes Theorem↗
PubMed2006

Mathematical modeling of intracellular signaling pathways.

Dynamic modeling and simulation of signal transduction pathways is an important topic in systems biology and is obtaining growing attention from researchers with experimental or theoretical background. Here we review attempts to analyze and model specific signaling systems. We review the structure of recurrent building blocks of signaling pathways and their integration into more comprehensive models, which enables the understanding of complex cellular processes. The variety of mechanisms found and modeling techniques used are illustrated with models of different signaling pathways. Focusing on the close interplay between experimental investigation of pathways and the mathematical representations of cellular dynamics, we discuss challenges and perspectives that emerge in studies of signaling systems.

Animals↗
PubMed2006

SBML-PET: a Systems Biology Markup Language-based parameter estimation tool.

UNLABELLED: The estimation of model parameters from experimental data remains a bottleneck for a major breakthrough in systems biology. We present a Systems Biology Markup Language (SBML) based Parameter Estimation Tool (SBML-PET). The tool is designed to enable parameter estimation for biological models including signaling pathways, gene regulation networks and metabolic pathways. SBML-PET supports import and export of the models in the SBML format. It can estimate the parameters by fitting a variety of experimental data from different experimental conditions. SBML-PET has a unique feature of supporting event definition in the SMBL model. SBML models can also be simulated in SBML-PET. Stochastic Ranking Evolution Strategy (SRES) is incorporated in SBML-PET for parameter estimation jobs. A classic ODE Solver called ODEPACK is used to solve the Ordinary Differential Equation (ODE) system. AVAILABILITY: http://sysbio.molgen.mpg.de/SBML-PET/. The website also contains detailed documentation for SBML-PET.

Algorithms↗
PubMed2006

A modelling approach to quantify dynamic crosstalk between the pheromone and the starvation pathway in baker's yeast.

Cells must be able to process multiple information in parallel and, moreover, they must also be able to combine this information in order to trigger the appropriate response. This is achieved by wiring signalling pathways such that they can interact with each other, a phenomenon often called crosstalk. In this study, we employ mathematical modelling techniques to analyse dynamic mechanisms and measures of crosstalk. We present a dynamic mathematical model that compiles current knowledge about the wiring of the pheromone pathway and the filamentous growth pathway in yeast. We consider the main dynamic features and the interconnections between the two pathways in order to study dynamic crosstalk between these two pathways in haploid cells. We introduce two new measures of dynamic crosstalk, the intrinsic specificity and the extrinsic specificity. These two measures incorporate the combined signal of several stimuli being present simultaneously and seem to be more stable than previous measures. When both pathways are responsive and stimulated, the model predicts that (a) the filamentous growth pathway amplifies the response of the pheromone pathway, and (b) the pheromone pathway inhibits the response of filamentous growth pathway in terms of mitogen activated protein kinase activity and transcriptional activity, respectively. Among several mechanisms we identified leakage of activated Ste11 as the most influential source of crosstalk. Moreover, we propose new experiments and predict their outcomes in order to test hypotheses about the mechanisms of crosstalk between the two pathways. Studying signals that are transmitted in parallel gives us new insights about how pathways and signals interact in a dynamical way, e.g., whether they amplify, inhibit, delay or accelerate each other.

Monte Carlo Method↗
PubMed2005

Alternative pathways as mechanism for the negative effects associated with overexpression of superoxide dismutase.

One of the most important antioxidant enzymes is superoxide dismutase (SOD), which catalyses the dismutation of superoxide radicals to hydrogen peroxide. The enzyme plays an important role in diseases like trisomy 21 and also in theories of the mechanisms of aging. But instead of being beneficial, intensified oxidative stress is associated with the increased expression of SOD and also studies on bacteria and transgenic animals show that high levels of SOD actually lead to increased lipid peroxidation and hypersensitivity to oxidative stress. Using mathematical models we investigate the question how overexpression of SOD can lead to increased oxidative stress, although it is an antioxidant enzyme. We consider the following possibilities that have been proposed in the literature: (i) Reaction of H(2)O(2) with CuZnSOD leading to hydroxyl radical formation. (ii) Superoxide radicals might reduce membrane damage by acting as radical chain breaker. (iii) While detoxifying superoxide radicals SOD cycles between a reduced and oxidized state. At low superoxide levels the intermediates might interact with other redox partners and increase the superoxide reductase (SOR) activity of SOD. This short-circuiting of the SOD cycle could lead to an increased hydrogen peroxide production. We find that only one of the proposed mechanisms is under certain circumstances able to explain the increased oxidative stress caused by SOD. But furthermore we identified an additional mechanism that is of more general nature and might be a common basis for the experimental findings. We call it the alternative pathway mechanism.

Humans↗
PubMed2005

Integrative model of the response of yeast to osmotic shock.

Integration of experimental studies with mathematical modeling allows insight into systems properties, prediction of perturbation effects and generation of hypotheses for further research. We present a comprehensive mathematical description of the cellular response of yeast to hyperosmotic shock. The model integrates a biochemical reaction network comprising receptor stimulation, mitogen-activated protein kinase cascade dynamics, activation of gene expression and adaptation of cellular metabolism with a thermodynamic description of volume regulation and osmotic pressure. Simulations agree well with experimental results obtained under different stress conditions or with specific mutants. The model is predictive since it suggests previously unrecognized features of the system with respect to osmolyte accumulation and feedback control, as confirmed with experiments. The mathematical description presented is a valuable tool for future studies on osmoregulation in yeast and-with appropriate modifications-other organisms. It also serves as a starting point for a comprehensive description of cellular signaling.

Cell Membrane↗
PubMed2005

Minimum information requested in the annotation of biochemical models (MIRIAM).

Most of the published quantitative models in biology are lost for the community because they are either not made available or they are insufficiently characterized to allow them to be reused. The lack of a standard description format, lack of stringent reviewing and authors' carelessness are the main causes for incomplete model descriptions. With today's increased interest in detailed biochemical models, it is necessary to define a minimum quality standard for the encoding of those models. We propose a set of rules for curating quantitative models of biological systems. These rules define procedures for encoding and annotating models represented in machine-readable form. We believe their application will enable users to (i) have confidence that curated models are an accurate reflection of their associated reference descriptions, (ii) search collections of curated models with precision, (iii) quickly identify the biological phenomena that a given curated model or model constituent represents and (iv) facilitate model reuse and composition into large subcellular models.

Biochemistry↗
PubMed2005

Biochemical network models simplified by balanced truncation.

Modelling of biochemical systems usually focuses on certain pathways, while the concentrations of so-called external metabolites are considered fixed. This approximation ignores feedback loops mediated by the environment, that is, via external metabolites and reactions. To achieve a more realistic, dynamic description that is still numerically efficient, we propose a new methodology: the basic idea is to describe the environment by a linear effective model of adjustable dimensionality. In particular, we (a) split the entire model into a subsystem and its environment, (b) linearize the environment model around a steady state, and (c) reduce its dimensionality by balanced truncation, an established method for large-scale model reduction. The reduced variables describe the dynamic modes in the environment that dominate its interaction with the subsystem. We compute metabolic response coefficients that account for complexity-reduced dynamics of the environment. Our simulations show that a dynamic environment model can improve the simulation results considerably, even if the environment model has been drastically reduced and if its kinetic parameters are only approximately known. The speed-up in computation gained by model reduction may become vital for parameter estimation in large cell models.

Kinetics↗
PubMed2004

Modelling the dynamics of the yeast pheromone pathway.

We present a mathematical model of the dynamics of the pheromone pathways in haploid yeast cells of mating type MATa after stimulation with pheromone alpha-factor. The model consists of a set of differential equations and describes the dynamics of signal transduction from the receptor via several steps, including a G protein and a scaffold MAP kinase cascade, up to changes in the gene expression after pheromone stimulation in terms of biochemical changes (complex formations, phosphorylations, etc.). The parameters entering the models have been taken from the literature or adapted to observed time courses or behaviour. Using this model we can follow the time course of the various complex formation processes and of the phosphorylation states of the proteins involved. Furthermore, we can explain the phenotype of more than a dozen well-characterized mutants and also the graded response of yeast cells to varying concentrations of the stimulating pheromone.

Cell Cycle Proteins↗
PubMed2004

Finding kinetic parameters using text mining.

The mathematical modeling and description of complex biological processes has become more and more important over the last years. Systems biology aims at the computational simulation of complex systems, up to whole cell simulations. An essential part focuses on solving a large number of parameterized differential equations. However, measuring those parameters is an expensive task, and finding them in the literature is very laborious. We developed a text mining system that supports researchers in their search for experimentally obtained parameters for kinetic models. Our system classifies full text documents regarding the question whether or not they contain appropriate data using a support vector machine. We evaluated our approach on a manually tagged corpus of 800 documents and found that it outperforms keyword searches in abstracts by a factor of five in terms of precision.

Computational Biology↗
PubMed2004

Alternative pathways might mediate toxicity of high concentrations of superoxide dismutase.

One of the most important antioxidant enzymes is superoxide dismutase (SOD), which catalyzes the dismutation of superoxide radicals to peroxide. The gene for CuZnSOD lies in humans on chromosome 21, and its activity is increased in patients with Down syndrome. However, instead of being beneficial, increased lipid peroxidation is associated with this increased expression, and also studies on bacteria and transgenic animals show that high levels of SOD actually lead to increased lipid peroxidation and hypersensitivity to oxidative stress. Using mathematical models, we investigated the question of how overexpression of SOD can lead to increased oxidative stress, although it is an antioxidant enzyme. We considered several possibilities that have been proposed in the literature, such as CuZnSOD-catalyzed hydroxyl radical formation, superoxide-mediated inhibition of membrane peroxidation, and short-circuiting of the Cu(I)ZnSOD/Cu(II)ZnSOD redox cycle. We found that one of the proposed mechanisms under certain circumstances is able to explain the increased oxidative stress caused by SOD. Furthermore, we identified an additional mechanism that agrees well with experimental observations. We call it the "alternative pathway" mechanism, because it depends on superoxide radicals having alternative pathways besides their reaction with SOD. The alternative pathway mechanism is a very general explanation for SOD-associated oxidative stress, because it does not depend on the specific type of SOD, nor on the redox status of the cell. We therefore think that it might be the common mechanism for the detrimental effects seen in cells and organisms with increased levels of the different forms of superoxide dismutase.

Antioxidants↗
PubMed2004

Inferring dynamic properties of biochemical reaction networks from structural knowledge.

Functional properties of biochemical networks depend on both the network structure and the kinetic parameters. Extensive data on metabolic network topologies have been collected in databases, but much less information is available about the kinetic constants or metabolite concentrations. Depending on the values of these parameters, metabolic fluxes and control coefficients may vary within a wide range. Nevertheless, some of the parameters may have little influence on the observables of interest. We address the question whether, despite uncertainty about kinetic parameters, probabilistic statements can be made about dynamic network features. To this end, we perform a variability analysis of the parameters: assuming that the parameters follow statistical distributions, we compute the resulting distributions of the network properties like metabolic fluxes, concentrations, or control coefficients by Monte Carlo simulation. In this manner, we study systematically the possible distributions arising from typical topologies of biochemical networks such as linear chains, branched networks, and signaling and gene expression cascades. This analysis reveals how much information about dynamic behavior can be drawn from structural knowledge.

Biochemistry↗
PubMed2002

Prediction of temporal gene expression. Metabolic opimization by re-distribution of enzyme activities.

A computational approach is used to analyse temporal gene expression in the context of metabolic regulation. It is based on the assumption that cells developed optimal adaptation strategies to changing environmental conditions. Time-dependent enzyme profiles are calculated which optimize the function of a metabolic pathway under the constraint of limited total enzyme amount. For linear model pathways it is shown that wave-like enzyme profiles are optimal for a rapid substrate turnover. For the central metabolism of yeast cells enzyme profiles are calculated which ensure long-term homeostasis of key metabolites under conditions of a diauxic shift. These enzyme profiles are in close correlation with observed gene expression data. Our results demonstrate that optimality principles help to rationalize observed gene expression profiles.

Databases as Topic↗
PubMedDate not supplied

A theory of optimal differential gene expression.

We investigate a model of optimal regulation, intended to describe large-scale differential gene expression. Relations between the optimal expression patterns and the function of genes are deduced from an optimality principle: the regulators have to maximise a fitness function which they influence directly via a cost term, and indirectly via their control on important cell variables, such as metabolic fluxes. According to the model, the optimal linear response to small perturbations reflects the regulators' functions, namely their linear influences on the cell variables. The optimal behaviour can be realised by a linear feedback mechanism. Known or assumed properties of response coefficients lead to predictions about regulation patterns. A symmetry relation predicted for deletion experiments is verified with gene expression data. Where the optimality assumption is valid, our results justify the use of expression data for functional annotation and for pathway reconstruction and suggest the use of linear factor models for the analysis of gene expression data.

Adaptation, Physiological↗