Bioinformatics Advance Access originally published online on March 20, 2008
Bioinformatics 2008 24(9):1191-1197; doi:10.1093/bioinformatics/btn103
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Bayesian inference of the sites of perturbations in metabolic pathways via Markov chain Monte Carlo
1Manchester Interdisciplinary Biocentre, The University of Manchester, Manchester M1 7DN, 2School of Chemistry, The University of Manchester, Manchester M13 9PL and 3School of Computer Science, The University of Manchester, Manchester M13 9PL, UK
*To whom correspondence should be addressed.
| Abstract |
|---|
Motivation: Genetic modifications or pharmaceutical interventions can influence multiple sites in metabolic pathways, and often these are distant from the primary effect. In this regard, the ability to identify target and off-target effects of a specific compound or gene therapy is both a major challenge and critical in drug discovery.
Results: We applied Markov Chain Monte Carlo (MCMC) for parameter estimation and perturbation identification in the kinetic modeling of metabolic pathways. Variability in the steady-state measurements in cells taken from a population can be caused by differences in initial conditions within the population, by variation of parameters among individuals and by possible measurement noise. MCMC-based parameter estimation is proposed as a method to help in inferring parameter distributions, taking into account uncertainties in the initial conditions and in the measurement data. The inferred parameter distributions are then used to predict changes in the network via a simple classification method. The proposed technique is applied to analyze changes in the pathways of pyruvate metabolism of mutants of Lactococcus lactis, based on previously published experimental data.
Availability: MATLAB code used in the simulations is available from ftp://anonymous@dbkweb.mib.man.ac.uk/pub/Bioinformatics_BJ.zip
Contact: bayujw{at}ieee.org
Supplementary information: Supplementary data are available at Bioinformatics online.
Associate Editor: Thomas Lengauer
1They are set such that the true value lies within the interval (exp(–m
), exp(m
)) where
is the log-variance and m is the log-mean. Note that exp(
0.4805)
2 and exp(
7.6872)
16.
Received on January 29, 2008; revised on February 22, 2008; accepted on March 15, 2008