The Community for Technology Leaders
RSS Icon
Subscribe
Issue No.01 - January-February (2011 vol.8)
pp: 27-35
Sean Mauch , California Institute of Technology, Pasadena
ABSTRACT
One can generate trajectories to simulate a system of chemical reactions using either Gillespie's direct method or Gibson and Bruck's next reaction method. Because one usually needs many trajectories to understand the dynamics of a system, performance is important. In this paper, we present new formulations of these methods that improve the computational complexity of the algorithms. We present optimized implementations, available from http://cain.sourceforge.net/>, that offer better performance than previous work. There is no single method that is best for all problems. Simple formulations often work best for systems with a small number of reactions, while some sophisticated methods offer the best performance for large problems and scale well asymptotically. We investigate the performance of each formulation on simple biological systems using a wide range of problem sizes. We also consider the numerical accuracy of the direct and the next reaction method. We have found that special precautions must be taken in order to ensure that randomness is not discarded during the course of a simulation.
INDEX TERMS
Biology and genetics, stochastic processes, algorithm design and analysis.
CITATION
Sean Mauch, "Efficient Formulations for Exact Stochastic Simulation of Chemical Systems", IEEE/ACM Transactions on Computational Biology and Bioinformatics, vol.8, no. 1, pp. 27-35, January-February 2011, doi:10.1109/TCBB.2009.47
REFERENCES
[1] L. Devroye, Non-Uniform Random Variate Generation. Springer, 1986.
[2] D.T. Gillespie, "A General Method for Numerically Simulating the Stochastic Time Evolution of Coupled Chemical Reactions," J. Computational Physics, vol. 22, no. 4, pp. 403-434, 1976.
[3] D.T. Gillespie, "Exact Stochastic Simulation of Coupled Chemical Reactions," J. Physical Chemistry, vol. 81, no. 25, pp. 2340-2361, 1977.
[4] M.A. Gibson and J. Bruck, "Efficient Exact Stochastic Simulation of Chemical Systems with Many Species and Many Channels," J. Physical Chemistry A, vol. 104, no. 9, pp. 1876-1889, 2000.
[5] D.T. Gillespie, "Stochastic Simulation of Chemical Kinetics," Ann. Rev. of Physical Chemistry, vol. 58, pp. 35-55, 2007.
[6] M. Matsumoto and T. Nishimura, "Mersenne Twister: A 623-Dimensionally Equidistributed Uniform Pseudorandom Number Generator," ACM Trans. Modelling and Computer Simulation, vol. 8, no. 1, pp. 3-30, http://www.math.sci.hiroshima-u.ac.jp/m-mat/ MTemt.html, Jan. 1998.
[7] G. Marsaglia and W.W. Tsang, "The Ziggurat Method for Generating Random Variables," J. Statistical Software, vol. 5, no. 8, pp. 1-7, http://www.jstatsoft.org/v05i08/, 2000.
[8] H. Rubin and B. Johnson, "Efficient Generation of Exponential and Normal Deviates," J. Statistical Computation and Simulation, vol. 76, no. 6, pp. 509-518, 2006.
[9] Y. Cao, H. Li, and L. Petzold, "Efficient Formulation of the Stochastic Simulation Algorithm for Chemically Reacting Systems," J. Chemical Physics, vol. 121, no. 9, pp. 4059-4067, 2004.
[10] J.M. McCollum, G.D. Peterson, C.D. Cox, M.L. Simpson, and N.F. Samatova, "The Sorting Direct Method for Stochastic Simulation of Biochemical Systems with Varying Reaction Execution Behavior," Computational Biology and Chemistry, vol. 30, no. 1, pp. 39-49, Feb. 2006.
[11] H. Li and L. Petzold, "Logarithmic Direct Method for Discrete Stochastic Simulation of Chemically Reacting Systems," technical report, Dept. of Computer Science, Univ. of California, Santa Barbara, http://www.engr.ucsb.edu~cse, 2006.
[12] A. Slepoy, A.P. Thompson, and S.J. Plimpton, "A Constant-Time Kinetic Monte Carlo Algorithm for Simulation of Large Biochemical Reaction Networks," J. Chemical Physics, vol. 128, no. 20,article no. 205101, 2008.
[13] T.H. Cormen, C.E. Leiserson, R.L. Rivest, and C. Stein, Introduction to Algorithms, second ed. MIT Press, 2001.
[14] D.J. Wilkinson, Stochastic Modelling for Systems Biology. CRC Press, 2006.
[15] Y. Cao and L. Petzold, "Accuracy Limitations and the Measurement of Errors in the Stochastic Simulation of Chemically Reacting Systems," J. Computational Physics, vol. 212, pp. 6-24, 2006.
[16] A.M. Kierzek, "STOCKS: STOChastic Kinetic Simulations of Biochemical Systems with the Gillespie Algorithm," Bioinformatics, vol. 18, no. 3, pp. 470-481, 2002.
[17] C.J. Proctor, C. Soti, R.J. Boys, C.S. Gillespie, D.P. Shanley, D.J. Wilkinson, and T.B.L. Kirkwood, "Modelling the Actions of Chaperones and Their Role in Ageing," Mechanisms of Ageing and Development, vol. 126, pp. 119-131, Jan. 2005.
[18] S. Ramsey, D. Orrell, and H. Bolouri, "Dizzy: Stochastic Simulation of Large-Scale Genetic Regulatory Networks," J. Bioinformatics and Computational Biology, vol. 3, pp. 415-436, Apr. 2005.
[19] S. Hoops, S. Sahle, R. Gauges, C. Lee, J. Pahle, N. Simus, M. Singhal, L. Xu, P. Mendes, and U. Kummer, "COPASI—A COmplex PAthway SImulator," Bioinformatics, vol. 22, pp. 3067-3074, 2006.
213 ms
(Ver 2.0)

Marketing Automation Platform Marketing Automation Tool