A bidomain threshold model of propagating calcium waves

We present a bidomain fire-diffuse-fire model that facilitates mathematical analysis of propagating waves of elevated intracellular calcium (Ca) in living cells. Modelling Ca release as a threshold process allows the explicit construction of travelling wave solutions to probe the dependence of Ca wave speed on physiologically important parameters such as the threshold for Ca release from the endoplasmic reticulum (ER) to the cytosol, the rate of Ca resequestration from the cytosol to the ER, and the total [Ca] (cytosolic plus ER). Interestingly, linear stability analysis of the bidomain fire-diffuse-fire model predicts the onset of dynamic wave instabilities leading to the emergence of Ca waves that propagate in a back-and-forth manner. Numerical simulations are used to confirm the presence of these so-called "tango waves" and the dependence of Ca wave speed on the total [Ca]. The original publication is available at www.springerlink.com (Journal of Mathematical Biology)

Receptors, sparks and waves in a fire-diffuse-fire framework for calcium release

Calcium ions are an important second messenger in living cells. Indeed calcium signals in the form of waves have been the subject of much recent experimental interest. It is now well established that these waves are composed of elementary stochastic release events (calcium puffs or sparks) from spatially localised calcium stores. The aim of this paper is to analyse how the stochastic nature of individual receptors within these stores combines to create stochastic behaviour on long timescales that may ultimately lead to waves of activity in a spatially extended cell model. Techniques from asymptotic analysis and stochastic phase-plane analysis are used to show that a large cluster of receptor channels leads to a release probability with a sigmoidal dependence on calcium density. This release probability is incorporated into a computationally inexpensive model of calcium release based upon a stochastic generalization of the Fire-Diffuse-Fire (FDF) threshold model. Numerical simulations of the model in one and two dimensions (with stores arranged on both regular and disordered lattices) illustrate that stochastic calcium release leads to the spontaneous production of calcium sparks that may merge to form saltatory waves. Illustrations of spreading circular waves, spirals and more irregular waves are presented. Furthermore, receptor noise is shown to generate a form of array enhanced coherence resonance whereby all calcium stores release periodically and simultaneously.

Sensitisation waves in a bidomain fire-diffuse-fire model of intracellular Ca²⁺ dynamics

We present a bidomain threshold model of intracellular calcium (Ca²⁺) dynamics in which, as suggested by recent experiments, the cytosolic threshold for Ca²⁺ liberation is modulated by the Ca²⁺ concentration in the releasing compartment. We explicitly construct stationary fronts and determine their stability using an Evans function approach. Our results show that a biologically motivated choice of a dynamic threshold, as opposed to a constant threshold, can pin stationary fronts that would otherwise be unstable. This illustrates a novel mechanism to stabilise pinned interfaces in continuous excitable systems. Our framework also allows us to compute travelling pulse solutions in closed form and systematically probe the wave speed as a function of physiologically important parameters. We find that the existence of travelling wave solutions depends on the time scale of the threshold dynamics, and that facilitating release by lowering the cytosolic threshold increases the wave speed. The construction of the Evans function for a travelling pulse shows that of the co-existing fast and slow solutions the slow one is always unstable.