Abstract

BackgroundTranslating 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.ResultsWe 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.ConclusionWe 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.

Highlights

  • Translating a known metabolic network into a dynamic model requires reasonable guesses of all enzyme parameters

  • 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 structure of metabolic networks has been characterised for many organisms [1,2,3], and metabolic fluxes in large networks [4,5,6] are successfully described by pathway- or constraint-based methods [710]

Read more

Summary

Introduction

Translating a known metabolic network into a dynamic model requires reasonable guesses of all enzyme parameters. The structure of metabolic networks has been characterised for many organisms [1,2,3], and metabolic fluxes in large networks [4,5,6] are successfully described by pathway- or constraint-based methods [710]. Such methods do not explain how the fluxes are evoked by the activities of enzymes and how they respond to moderate perturbations. Parameter estimation is an integral part of kinetic modelling, and model fitting is currently receiving increasing attention [25,26,27,28,29]

Methods
Results
Discussion
Conclusion
Full Text
Published version (Free)

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call