Stochastic simulation of the spatio-temporal dynamics of reaction-diffusion systems: the case for the bicoid gradient

More Information | Back to archive
Full Text of this article Full article [PDF] (2,31 MB)
doi doi:10.2390/biecoll-jib-2010-150
submission June 10, 2010
last revision October 26, 2010
published November 23, 2010
NCBI PubMed PubMed ID 21098882

Paola Lecca, Adaoha E. C. Ihekwaba, Lorenzo Dematté and Corrado Priami

Correspondence should be addressed to:
Paola Lecca
The Microsoft Research - University of Trento, Centre for Computational and Systems Biology, Povo (Trento), Italy
ue.ibsoc@nullaccel


Abstract

Reaction-diffusion systems are mathematical models that describe how the concentrations of substances distributed in space change under the influence of local chemical reactions, and diffusion which causes the substances to spread out in space. The classical representation of a reaction-diffusion system is given by semi-linear parabolic partial differential equations, whose solution predicts how diffusion causes the concentration field to change with time. This change is proportional to the diffusion coefficient. If the solute moves in a homogeneous system in thermal equilibrium, the diffusion coefficients are constants that do not depend on the local concentration of solvent and solute. However, in nonhomogeneous and structured media the assumption of constant intracellular diffusion coefficient is not necessarily valid, and, consequently, the diffusion coefficient is a function of the local concentration of solvent and solutes. In this paper we propose a stochastic model of reaction-diffusion systems, in which the diffusion coefficients are function of the local concentration, viscosity and frictional forces. We then describe the software tool Redi (REaction-DIffusion simulator) which we have developed in order to implement this model into a Gillespie-like stochastic simulation algorithm. Finally, we show the ability of our model implemented in the Redi tool to reproduce the observed gradient of the bicoid protein in the Drosophila Melanogaster embryo. With Redi, we were able to simulate with an accuracy of 1% the experimental spatio-temporal dynamics of the bicoid protein, as recorded in time-lapse experiments obtained by direct measurements of transgenic bicoidenhanced green fluorescent protein.

Reference

Paola Lecca, Adaoha E. C. Ihekwaba, Lorenzo Dematté and Corrado Priami. Stochastic simulation of the spatio-temporal dynamics of reaction-diffusion systems: the case for the bicoid gradient. Journal of Integrative Bioinformatics, 7(1):150, 2010. Online Journal: http://journal.imbio.de/index.php?paper_id=150
imprint | sitemap | credits | top