Stochastic reaction-diffusion modeling of calcium dynamics in 3D dendritic spines of Purkinje cells

Calcium (Ca²⁺) is a second messenger assumed to control changes in synaptic strength in the form of both long-term depression and long-term potentiation at Purkinje cell dendritic spine synapses via inositol trisphosphate (IP₃)-induced Ca²⁺ release. These Ca²⁺ transients happen in response to stimuli from parallel fibers (PFs) from granule cells and climbing fibers (CFs) from the inferior olivary nucleus. These events occur at low numbers of free Ca²⁺, requiring stochastic single-particle methods when modeling them. We use the stochastic particle simulation program MCell to simulate Ca²⁺ transients within a three-dimensional Purkinje cell dendritic spine. The model spine includes the endoplasmic reticulum, several Ca²⁺ transporters, and endogenous buffer molecules. Our simulations successfully reproduce properties of Ca²⁺ transients in different dynamical situations. We test two different models of the IP₃ receptor (IP₃R). The model with nonlinear concentration response of binding of activating Ca²⁺ reproduces experimental results better than the model with linear response because of the filtering of noise. Our results also suggest that Ca²⁺-dependent inhibition of the IP₃R needs to be slow to reproduce experimental results. Simulations suggest the experimentally observed optimal timing window of CF stimuli arises from the relative timing of CF influx of Ca²⁺ and IP₃ production sensitizing IP₃R for Ca²⁺-induced Ca²⁺ release. We also model ataxia, a loss of fine motor control assumed to be the result of malfunctioning information transmission at the granule to Purkinje cell synapse, resulting in a decrease or loss of Ca²⁺ transients. Finally, we propose possible ways of recovering Ca²⁺ transients under ataxia.