Modeling stochastic calcium dynamics in the dendritic spines: a hybrid algorithm
© Kalantzis et al; licensee BioMed Central Ltd. 2008
Published: 11 July 2008
Stochastic fluctuation plays an important role in biology. In a small subcellular compartment like dendritic spines of pyramidal hippocampal cells, the basal level of Ca2+ concentration (~100 nM) entails six Ca2+ ions. In such a small volume, the stochastic nature of chemical reactions considerably influences intracellular signal transductions . In fact, Ca+2 is an important intracellular second messenger regulating many postsynaptic enzymes which in turn trigger modifications of synaptic strength or structure. Most of these interactions are highly nonlinear. For these reasons, we need an efficient and accurate stochastic simulation method to study synaptic functions.
One solution to this problem is a stochastic simulation algorithm (SSA) developed by Gillespie for simulating spatially homogeneous chemical systems . For spatially inhomogeneous system such as intracellular environments, a major modification of the Gillespie algorithm is required. For example, Bernstein  subdivides the simulated volume into small cells and models the diffusion of molecules between neighboring elements as a Markov transitions process treated as a special kind of "reaction". This method was shown to accurately reproduce the macroscopic limit of reaction diffusion system with cubic autocatalytic reactions. However, a much higher simulation speed is desired for practical applications. Hybrid algorithms were proposed to solve this problem [4, 5]. In these algorithms, the reactions are partitioned into fast and slow based on the propensities and the population of the different species: the former (fast component) is simulated deterministically and the latter stochastically.
However, such a hybrid scheme becomes problematic and produces inaccurate results especially in a stiff system. For example, dendritic spines contain NMDA receptors whose kinetic transitions are relatively slow but inject fast diffusing Ca2+ ions into the spines which control the subsequent biochemical signaling cascades. Because of this stiffness of our system, the propensities of the "slow" subset of reactions do not change smoothly in time and consequently the fast-slow partitioning failed to describe the calcium dynamics. In this paper, we describe a novel hybrid stochastic algorithm to overcome this problem. In short, we developed a different partitioning strategy to avoid the conflicts between fast and slow systems in the simulation. We compare our results with the fully stochastic algorithm and the fully deterministic methods. Using this hybrid algorithm, we show how the stochastic fluctuations of Ca2+ dynamics influence the induction of synaptic plasticity.
- Shouval HZ, Kalantzis G: Stochastic properties of synaptic transmission affect the shape of spike time dependent plasticity curves. J Neurophysiol. 2005, 93: 1069-1073. 10.1152/jn.00504.2004.View ArticlePubMedGoogle Scholar
- Gillespie DT: A general method for numerically simulating the stochastic evolution of coupled chemical reactions. J Comput Phys. 1976, 22: 403-10.1016/0021-9991(76)90041-3.View ArticleGoogle Scholar
- Bernstein D: Simulating mesoscopic reaction-diffusion systems using the Gillespie algorithm. Phys Rev E. 2005, 71: 041103-10.1103/PhysRevE.71.041103.View ArticleGoogle Scholar
- Chiam KH, Tan C, Bhargava V, Rajagopal G: Hybrid simulations of stochastic reaction-diffusion processes for modeling intracellular signaling pathways. Phys Rev E. 2006, 74: 51910-10.1103/PhysRevE.74.051910.View ArticleGoogle Scholar
- Alfonsi A, et al: Adaptive simulation of hybrid stochastic and deterministic models for biochemical systems. ESAIM Proceedings. 2005, 14: 1-13.Google Scholar
This article is published under license to BioMed Central Ltd.