Email updates

Keep up to date with the latest news and content from BMC Neuroscience and BioMed Central.

This article is part of the supplement: Seventeenth Annual Computational Neuroscience Meeting: CNS*2008

Open Access Poster presentation

Modeling stochastic calcium dynamics in the dendritic spines: a hybrid algorithm

Georgios Kalantzis1, Yoshihisa Kubota1 and Harel Z Shouval123*

Author Affiliations

1 Department of Neurobiology and Anatomy, The University of Texas Medical School, Houston, Texas, USA

2 Department of Biomedical Engineering, The University of Texas, Austin, Texas, USA

3 Institute for Brain and Neural Systems, Brown University, Providence, Rhode Island, USA

For all author emails, please log on.

BMC Neuroscience 2008, 9(Suppl 1):P86  doi:10.1186/1471-2202-9-S1-P86

The electronic version of this article is the complete one and can be found online at:

Published:11 July 2008

© 2008 Kalantzis et al; licensee BioMed Central Ltd.

Poster presentation

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 [1]. 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 [2]. For spatially inhomogeneous system such as intracellular environments, a major modification of the Gillespie algorithm is required. For example, Bernstein [3] 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.


  1. Shouval HZ, Kalantzis G: Stochastic properties of synaptic transmission affect the shape of spike time dependent plasticity curves.

    J Neurophysiol 2005, 93:1069-1073. PubMed Abstract | Publisher Full Text OpenURL

  2. Gillespie DT: A general method for numerically simulating the stochastic evolution of coupled chemical reactions.

    J Comput Phys 1976, 22:403. Publisher Full Text OpenURL

  3. Bernstein D: Simulating mesoscopic reaction-diffusion systems using the Gillespie algorithm.

    Phys Rev E 2005, 71:041103. Publisher Full Text OpenURL

  4. 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. Publisher Full Text OpenURL

  5. Alfonsi A, et al.: Adaptive simulation of hybrid stochastic and deterministic models for biochemical systems.

    ESAIM Proceedings 2005, 14:1-13. OpenURL