Bioinformatics Advance Access originally published online on May 30, 2007
Bioinformatics 2007 23(15):1969-1977; doi:10.1093/bioinformatics/btm278
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Computational methods for diffusion-influenced biochemical reactions
ski 1,*1CWI (Center for Mathematics and Computer Science), Kruislaan 413 and 2Section Computational Science, Faculty of Science, University of Amsterdam, Kruislaan 403, 1098 SJ Amsterdam, The Netherlands
*To whom correspondence should be addressed.
| Abstract |
|---|
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/
Contact: m.dobrzynski{at}cwi.nl
Supplementary information: Supplementary data are available at Bioinformatics online.
Associate Editor: Jonathan Wren
Received on January 31, 2007; revised on April 20, 2007; accepted on May 17, 2007
This article has been cited by other articles:
![]() |
G. M. Dancik and K. S. Dorman mlegp: statistical analysis for computer models of biological systems using R Bioinformatics, September 1, 2008; 24(17): 1966 - 1967. [Abstract] [Full Text] [PDF] |
||||
