PubMed Health⌕ Search

Biomedical subjects

Wolfram Liebermeister

Publications and source records attributed to Wolfram Liebermeister.

11 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

A comprehensive library of fluorescent transcriptional reporters for Escherichia coli.

E. coli is widely used for systems biology research; there exists a need, however, for tools that can be used to accurately and comprehensively measure expression dynamics in individual living cells. To address this we present a library of transcriptional fusions of gfp to each of about 2,000 different promoters in E. coli K12, covering the great majority of the promoters in the organism. Each promoter fusion is expressed from a low-copy plasmid. We demonstrate that this library can be used to obtain highly accurate dynamic measurements of promoter activity on a genomic scale, in a glucose-lactose diauxic shift experiment. The library allowed detection of about 80 previously uncharacterized transcription units in E. coli, including putative internal promoters within previously known operons, such as the lac operon. This library can serve as a tool for accurate, high-resolution analysis of transcription networks in living E. coli cells.

Escherichia coli↗
PubMed2005

Response to temporal parameter fluctuations in biochemical networks.

Metabolic response coefficients describe how variables in metabolic systems, like steady state concentrations, respond to small changes of kinetic parameters. To extend this concept to temporal parameter fluctuations, we define spectral response coefficients that relate Fourier components of concentrations and fluxes to Fourier components of the underlying parameters. It is also straightforward to generalize other concepts from metabolic control theory, such as control coefficients with their summation and connectivity theorems. The first-order response coefficients describe forced oscillations caused by small harmonic oscillations of single parameters: they depend on the driving frequency and comprise the phases and amplitudes of the concentrations and fluxes. Close to a Hopf bifurcation, resonance can occur: as an example, we study the spectral densities of concentration fluctuations arising from the stochastic nature of chemical reactions. Second-order response coefficients describe how perturbations of different frequencies interact by mode coupling, yielding higher harmonics in the metabolic response. The temporal response to small parameter fluctuations can be computed by Fourier synthesis. For a model of glycolysis, this approximation remains fairly accurate even for large relative fluctuations of the parameters.

Animals↗
PubMed2005

Predicting physiological concentrations of metabolites from their molecular structure.

Physiological concentrations of metabolites can partly be explained by their molecular structure. We hypothesize that substances containing certain chemical groups show increased or decreased concentration in cells. We consider here, as chemical groups, local atomic configurations, describing an atom, its bonds, and its direct neighbor atoms. To test our hypothesis, we fitted a linear statistical model that relates experimentally determined logarithmic concentrations to feature vectors containing count numbers of the chemical groups. In order to determine chemical groups that have a clear effect on the concentration, we use a regularized (lasso) regression. In a dataset on 41 substances of central metabolism in different organisms, we found that the physical concentrations are increased by the occurrence of amino and hydroxyl groups, while aldehydes, ketones, and phosphates show decreased concentrations. The model explains about 22% of the variance of the logarithmic mean concentrations.

Animals↗
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

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

Linear modes of gene expression determined by independent component analysis.

MOTIVATION: The expression of genes is controlled by specific combinations of cellular variables. We applied Independent Component Analysis (ICA) to gene expression data, deriving a linear model based on hidden variables, which we term 'expression modes'. The expression of each gene is a linear function of the expression modes, where, according to the ICA model, the linear influences of different modes show a minimal statistical dependence, and their distributions deviate sharply from the normal distribution. RESULTS: Studying cell cycle-related gene expression in yeast, we found that the dominant expression modes could be related to distinct biological functions, such as phases of the cell cycle or the mating response. Analysis of human lymphocytes revealed modes that were related to characteristic differences between cell types. With both data sets, the linear influences of the dominant modes showed distributions with large tails, indicating the existence of specifically up- and downregulated target genes. The expression modes and their influences can be used to visualize the samples and genes in low-dimensional spaces. A projection to expression modes helps to highlight particular biological functions, to reduce noise, and to compress the data in a biologically sensible way.

Algorithms↗
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↗