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

Dendritic spines act as computational units and must adapt their responses according to their activation history. Calcium influx acts as the first signaling step during postsynaptic activation and is a determinant of synaptic weight change. Dendritic spines also come in a variety of sizes and shapes. To probe the relationship between calcium dynamics and spine morphology, we used a stochastic reaction-diffusion model of calcium dynamics in idealized and realistic geometries. We show that despite the stochastic nature of the various calcium channels, receptors, and pumps, spine size and shape can separately modulate calcium dynamics and subsequently synaptic weight updates in a deterministic manner. The relationships between calcium dynamics and spine morphology identified in idealized geometries also hold in realistic geometries suggesting that there are geometrically determined deterministic relationships that may modulate synaptic weight change.

Calcium (Ca2+) imaging reveals a variety of correlated firing in cultures of dissociated hippocampal neurons, pinpointing the non‐synaptic paracrine release of glutamate as a possible mediator for such firing patterns, although the biophysical underpinnings remain unknown. An intriguing possibility is that extracellular glutamate could bind metabotropic receptors linked with inositol trisphosphate (IP3) mediated release of Ca2+ from the endoplasmic reticulum of individual neurons, thereby modulating neural activity in combination with sarco/endoplasmic reticulum Ca2+ transport ATPase (SERCA) and voltage‐gated Ca2+ channels (VGCC). However, the possibility that such release may occur in different neuronal compartments and can be inherently stochastic poses challenges in the characterization of such interplay between various Ca2+ channels. Here we deploy biophysical modeling in association with Monte Carlo parameter sampling to characterize such interplay and successfully predict experimentally observed Ca2+ patterns. The results show that the neurotransmitter level at the plasma membrane is the extrinsic source of heterogeneity in somatic Ca2+ transients. Our analysis, in particular, identifies the origin of such heterogeneity to an intrinsic differentiation of hippocampal neurons in terms of multiple cellular properties pertaining to intracellular Ca2+ signaling, such as VGCC, IP3 receptor, and SERCA expression. In the future, the biophysical model and parameter estimation approach used in this study can be upgraded to predict the response of a system of interconnected neurons.

The modeling and simulation of stochastic reaction–diffusion processes is a topic of steady interest that is approached with a wide range of methods. At the level of particle-resolved descriptions, where chemical reactions are coupled to the spatial diffusion of individual particles, there exist comprehensive numerical simulation schemes, while the corresponding mathematical formalization is relatively underdeveloped. The aim of this paper is to provide a framework to systematically formulate the probabilistic evolution equation, termed chemical diffusion master equation (CDME), that governs particle-based stochastic reaction–diffusion processes. To account for the non-conserved and unbounded particle number of this type of open systems, we employ a classical analogue of the quantum mechanical Fock space that contains the symmetrized probability densities of the many-particle configurations in space. Following field-theoretical ideas of second quantization, we introduce creation and annihilation operators that act on single-particle densities and provide natural representations of symmetrized probability densities as well as of reaction and diffusion operators. These operators allow us to consistently and systematically formulate the CDME for arbitrary reaction schemes. The resulting form of the CDME further serves as the foundation to derive more coarse-grained descriptions of reaction–diffusion dynamics. In this regard, we show that a discretization of the evolution equation by projection onto a Fock subspace generated by a finite set of single-particle densities leads to a generalized form of the well-known reaction–diffusion master equation, which supports non-local reactions between grid cells and which converges properly in the continuum limit.