Computational methods for diffusion-influenced biochemical reactions

MOTIVATION We compare stochastic computational methods accounting for space and discrete nature of reactants in biochemical systems. Implementations based on Brownian dynamics (BD) and the reaction-diffusion master equation are applied to a simplified gene expression model and to a signal transduction pathway in Escherichia coli. RESULTS In the regime where the number of molecules is small and reactions are diffusion-limited predicted fluctuations in the product number vary between the methods, while the average is the same. Computational approaches at the level of the reaction-diffusion master equation compute the same fluctuations as the reference result obtained from the particle-based method if the size of the sub-volumes is comparable to the diameter of reactants. Using numerical simulations of reversible binding of a pair of molecules we argue that the disagreement in predicted fluctuations is due to different modeling of inter-arrival times between reaction events. Simulations for a more complex biological study show that the different approaches lead to different results due to modeling issues. Finally, we present the physical assumptions behind the mesoscopic models for the reaction-diffusion systems. AVAILABILITY Input files for the simulations and the source code of GMP can be found under the following address: http://www.cwi.nl/projects/sic/bioinformatics2007/

[1]  E. Cox,et al.  Real-Time Kinetics of Gene Activity in Individual Bacteria , 2005, Cell.

[2]  Bastien Chopard,et al.  Multiparticle Lattice Gas Automata For Reaction Diffusion Systems , 1994 .

[3]  A. van Oudenaarden,et al.  Noise Propagation in Gene Networks , 2005, Science.

[4]  D. Bray,et al.  Simulated Diffusion of Phosphorylated CheY through the Cytoplasm of Escherichia coli , 2005, Journal of bacteriology.

[5]  T. Elston,et al.  Stochasticity in gene expression: from theories to phenotypes , 2005, Nature Reviews Genetics.

[6]  D. Bray,et al.  Stochastic simulation of chemical reactions with spatial resolution and single molecule detail , 2004, Physical biology.

[7]  K. Kaneko,et al.  Transitions induced by the discreteness of molecules in a small autocatalytic system. , 2000, Physical review letters.

[8]  P. R. ten Wolde,et al.  Green's-function reaction dynamics: a particle-based approach for simulating biochemical networks in time and space. , 2005, The Journal of chemical physics.

[9]  D. Gillespie.  A rigorous derivation of the chemical master equation , 1992 .

[10]  D. Gillespie.  A General Method for Numerically Simulating the Stochastic Time Evolution of Coupled Chemical Reactions , 1976 .

[11]  N. Agmon,et al.  Theory of reversible diffusion‐influenced reactions , 1990 .

[12]  Raymond Kapral,et al.  Diffusion and reaction in crowded environments , 2007 .

[13]  P. T. Wolde,et al.  Simulating biochemical networks at the particle level and in time and space: Green's function reaction dynamics. , 2005 .

[14]  C. Brooks.  Computer simulation of liquids , 1989 .

[15]  A. Kierzek,et al.  The Effect of Transcription and Translation Initiation Frequencies on the Stochastic Fluctuations in Prokaryotic Gene Expression* , 2001, The Journal of Biological Chemistry.

[16]  D. Gillespie.  Exact Stochastic Simulation of Coupled Chemical Reactions , 1977 .

[17]  David Fange,et al.  Noise-Induced Min Phenotypes in E. coli , 2006, PLoS Comput. Biol..

[18]  A. Arkin,et al.  It's a noisy business! Genetic regulation at the nanomolar scale. , 1999, Trends in genetics : TIG.

[19]  Mansour,et al.  Reaction-diffusion master equation: A comparison with microscopic simulations. , 1996, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[20]  Martin Howard,et al.  Stochastic model for Soj relocation dynamics in Bacillus subtilis. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[21]  Sidney Redner,et al.  A guide to first-passage processes , 2001 .

[22]  C. Lumsden,et al.  Stochastic Simulation of Coupled Reaction-Diffusion Processes , 1996 .

[23]  Johan Hattne,et al.  Stochastic reaction-diffusion simulation with MesoRD , 2005, Bioinform..

[24]  Jaap A. Kaandorp,et al.  Spatial stochastic modelling of the phosphoenolpyruvate-dependent phosphotransferase (PTS) pathway in Escherichia coli , 2006, Bioinform..

[25]  R Metzler,et al.  The future is noisy: the role of spatial fluctuations in genetic switching. , 2001, Physical review letters.

[26]  U. Bhalla.  Signaling in small subcellular volumes. I. Stochastic and diffusion effects on individual pathways. , 2004, Biophysical journal.

[27]  J. Stoyanov.  A Guide to First‐passage Processes , 2003 .

[28]  Kunihiko Kaneko,et al.  Discreteness-induced stochastic steady state in reaction diffusion systems : self-consistent analysis and stochastic simulations , 2004, physics/0409027.

[29]  P. Swain,et al.  Gene Regulation at the Single-Cell Level , 2005, Science.

[30]  Jeroen S. van Zon,et al.  Diffusion of transcription factors can drastically enhance the noise in gene expression. , 2006, Biophysical journal.

[31]  Dan ie l T. Gil lespie.  A rigorous derivation of the chemical master equation , 1992 .

[32]  Karen Lipkow,et al.  Changing Cellular Location of CheZ Predicted by Molecular Simulations , 2006, PLoS Comput. Biol..

[33]  P. R. ten Wolde,et al.  Simulating biochemical networks at the particle level and in time and space: Green's function reaction dynamics. , 2005, Physical review letters.

[34]  J. Elf,et al.  Spontaneous separation of bi-stable biochemical systems into spatial domains of opposite phases. , 2004, Systems biology.

[35]  N. Kampen,et al.  Stochastic processes in physics and chemistry , 1981 .

[36]  R. Metzler.  The Future is Noisy , 2001 .

[37]  Eric Renshaw,et al.  Spatial heterogeneity and the stability of reaction states in autocatalysis. , 2002, Physical review. E, Statistical, nonlinear, and soft matter physics.

[38]  Sandeep Krishna,et al.  Stochastic simulations of the origins and implications of long-tailed distributions in gene expression. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[39]  Kunihiko Kaneko,et al.  Molecular discreteness in reaction-diffusion systems yields steady states not seen in the continuum limit. , 2003, Physical review. E, Statistical, nonlinear, and soft matter physics.

[40]  C. Gardiner.  Handbook of Stochastic Methods , 1983 .

[41]  J. Blom,et al.  Why the phosphotransferase system of Escherichia coli escapes diffusion limitation. , 2003, Biophysical journal.

[42]  A. Arkin,et al.  Stochastic kinetic analysis of developmental pathway bifurcation in phage lambda-infected Escherichia coli cells. , 1998, Genetics.

[43]  Masaru Tomita,et al.  Space in systems biology of signaling pathways – towards intracellular molecular crowding in silico , 2005, FEBS letters.

[44]  Michael A. Gibson,et al.  Efficient Exact Stochastic Simulation of Chemical Systems with Many Species and Many Channels , 2000 .

[45]  A. Oudenaarden,et al.  Enhancement of cellular memory by reducing stochastic transitions , 2005, Nature.

[46]  S. Solomon,et al.  The importance of being discrete: life always wins on the surface. , 1999, Proceedings of the National Academy of Sciences of the United States of America.

[47]  W. Ebeling.  Stochastic Processes in Physics and Chemistry , 1995 .

[48]  J. Marko,et al.  How do site-specific DNA-binding proteins find their targets? , 2004, Nucleic acids research.

[49]  P. R. ten Wolde,et al.  Green's-function reaction dynamics: a particle-based approach for simulating biochemical networks in time and space. , 2005, The Journal of chemical physics.

[50]  Benjamin B. Kaufmann,et al.  Contributions of low molecule number and chromosomal positioning to stochastic gene expression , 2005, Nature Genetics.