Bioinformatics Advance Access originally published online on October 26, 2006
Bioinformatics 2007 23(1):84-91; doi:10.1093/bioinformatics/btl552
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Simulated maximum likelihood method for estimating kinetic rates in gene expression
1 Advanced Computational Modelling Centre, University of Queensland Brisbane, QLD 4072, Australia
2 Department of Mathematics, Hubei University of Technology Wuhan, Hubei 430068, P. R. China
3 School of Information Technology, Charles Sturt University Bathurst, NSW 2795, Australia
*To whom correspondence should be addressed.
Motivation: Kinetic rate in gene expression is a key measurement of the stability of gene products and gives important information for the reconstruction of genetic regulatory networks. Recent developments in experimental technologies have made it possible to measure the numbers of transcripts and protein molecules in single cells. Although estimation methods based on deterministic models have been proposed aimed at evaluating kinetic rates from experimental observations, these methods cannot tackle noise in gene expression that may arise from discrete processes of gene expression, small numbers of mRNA transcript, fluctuations in the activity of transcriptional factors and variability in the experimental environment.
Results: In this paper, we develop effective methods for estimating kinetic rates in genetic regulatory networks. The simulated maximum likelihood method is used to evaluate parameters in stochastic models described by either stochastic differential equations or discrete biochemical reactions. Different types of non-parametric density functions are used to measure the transitional probability of experimental observations. For stochastic models described by biochemical reactions, we propose to use the simulated frequency distribution to evaluate the transitional density based on the discrete nature of stochastic simulations. The genetic optimization algorithm is used as an efficient tool to search for optimal reaction rates. Numerical results indicate that the proposed methods can give robust estimations of kinetic rates with good accuracy.
Contact: tian{at}maths.uq.edu.au
Received on June 23, 2006; revised on October 5, 2006; accepted on October 23, 2006
CiteULike 
Connotea 
Del.icio.us What's this?
This article has been cited by other articles:
|
|
 |
|
  L. Zamparo and T. J. Perkins Statistical lower bounds on protein copy number from fluorescence expression images Bioinformatics, October 15, 2009; 25(20): 2670 - 2676. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
E. Cinquemani, A. Milias-Argeitis, S. Summers, and J. Lygeros Stochastic dynamics of genetic networks: modelling and parameter identification Bioinformatics, December 1, 2008; 24(23): 2748 - 2754. [Abstract] [Full Text] [PDF]
|
|