Bioinformatics Vol. 18 no. 3 2002
Pages 470-481
© 2002 Oxford University Press
STOCKS: STOChastic Kinetic Simulations of biochemical systems with Gillespie algorithm
Institute of Biochemistry and Biophysics, Polish Academy of Sciences, Pawinskiego 5a, 02-106 Warszawa, Poland
Received on June 5, 2001
; revised on August 22, 2001
; accepted on October 4, 2001
Motivation: The availability of a huge amount of molecular data concerning various biochemical reactions provoked numerous attempts to study the dynamics of cellular processes by means of kinetic models and computer simulations. Biochemical processes frequently involve small numbers of molecules (e.g. a few molecules of a transcriptional regulator binding to one ‘molecule’ of a DNA regulatory region). Such reactions are subject to significant stochastic fluctuations. Monte Carlo methods must be employed to study the functional consequences of the fluctuations and simulate processes that cannot be modelled by continuous fluxes of matter. This provides the motivation to develop software dedicated to Monte Carlo simulations of cellular processes with the rigorously proven Gillespie algorithm.
Results: STOCKS, software for the stochastic kinetic simulation of biochemical processes is presented. The program uses a rigorously derived Gillespie algorithm that has been shown to be applicable to the study of prokaryotic gene expression. Features dedicated to the study of cellular processes are implemented, such as the possibility to study a process in the range of several cell generations with the application of a simple cell division model. Taking expression of Escherichia coli beta-galactosidase as an example, it is shown that the program is able to simulate systems composed of reactions varying in several orders of magnitude by means of reaction rates and the numbers of molecules involved.
Availability: The software is available at ftp://ibbrain.ibb.waw.pl/stocksand http://www.ibb.waw.pl/stocks.
Supplementary information: Parameters of the model of prokaryotic gene expression are available in example files of software distribution.
Contact: andrzejk{at}ibb.waw.pl
CiteULike 
Connotea 
Del.icio.us What's this?
This article has been cited by other articles:
|
|
 |
|
  J. Pahle Biochemical simulations: stochastic, approximate stochastic and hybrid approaches Brief Bioinform, January 16, 2009; (2009) bbn050v1. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 T. Tian, S. Xu, J. Gao, and K. Burrage Simulated maximum likelihood method for estimating kinetic rates in gene expression Bioinformatics, January 1, 2007; 23(1): 84 - 91. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
K. Burrage, L. Hood, and M. A. Ragan Advanced computing for systems biology Brief Bioinform, December 1, 2006; 7(4): 390 - 398. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
A. Ng, B. Bursteinas, Q. Gao, E. Mollison, and M. Zvelebil Resources for integrative systems biology: from data through databases to networks and dynamic system models Brief Bioinform, December 1, 2006; 7(4): 318 - 330. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
C. Sanford, M. L.K. Yip, C. White, and J. Parkinson Cell++--simulating biochemical pathways Bioinformatics, December 1, 2006; 22(23): 2918 - 2925. [Abstract] [Full Text] [PDF]
|
|