 Research
 Open Access
 Published:
What Is Required for Neuronal Calcium Waves? A Numerical Parameter Study
The Journal of Mathematical Neuroscience volume 8, Article number: 9 (2018)
Abstract
Neuronal calcium signals propagating by simple diffusion and reaction with mobile and stationary buffers are limited to cellular microdomains. The distance intracellular calcium signals can travel may be significantly increased by means of calciuminduced calcium release from internal calcium stores, notably the endoplasmic reticulum. The organelle, which can be thought of as a cellwithinacell, is able to sequester large amounts of cytosolic calcium ions via SERCA pumps and selectively release them into the cytosol through ryanodine receptor channels leading to the formation of calcium waves. In this study, we set out to investigate the basic properties of such dendritic calcium waves and how they depend on the three parameters dendrite radius, ER radius and ryanodine receptor density in the endoplasmic membrane. We demonstrate that there are stable and abortive regimes for calcium waves, depending on the above morphological and physiological parameters. In stable regimes, calcium waves can travel across long dendritic distances, similar to electrical action potentials. We further observe that abortive regimes exist, which could be relevant for spiketiming dependent plasticity, as travel distances and wave velocities vary with changing intracellular architecture. For some of these regimes, analytic functions could be derived that fit the simulation data. In parameter spaces, that are nontrivially influenced by the threedimensional calcium concentration profile, we were not able to derive such a functional description, demonstrating the mathematical requirement to model and simulate biochemical signaling in threedimensional space.
Introduction
Intracellular calcium signals define a transition point between electrical signals and biochemical responses in neurons. While basal calcium concentrations in the cytosol are very low, neurons can modulate local cytosolic calcium concentrations to induce microdomain calcium signals, which can integrate to produce longer ranging signal propagation towards the soma. There they reach the nucleus and trigger gene transcription responses relevant for learning [1–4] and neuroprotection [5–8]. Cellular calcium signals are shaped by calcium transport mechanisms embedded in the plasma membrane, across which calcium can be bidirectionally exchanged between cytosol and the extracellular space (see [9] for an overview). In addition, intracellular organelles like mitochondria and the endoplasmic reticulum (ER) function as large calcium stores [10–16]. Organelle membranes are equipped with calcium exchange mechanisms that allow transport of calcium from the cytosol into the organelles or vice versa. Some of these mechanisms—notably the ryanodine receptor channel (RyR) in the ER membrane, which is able to release large amounts of calcium from the ER into the cytosol—have a positive feedback property: Their opening is facilitated by the presence of cytosolic calcium and will trigger the release of even more calcium through surrounding channels. This calciuminduced calcium release mechanism (CICR) overcomes the limited reach of purely diffusive calcium signals in a buffered regime. CICR in neurons has been implicated in various neurodegenerative diseases [17–23] and has been studied with respect to CICR wavelike properties [24–38]. The importance of CICR wavelike dynamics therefore merits thorough investigation. One major question is how the cellular and intracellular architecture can shape and regulate calcium signals that are able to propagate over long distances, for instance in the context of synapsetonucleus communication. It was shown in [26] that variations in the activation of Ca^{2+} store release through inositol 1,4,5trisphosphate (IP_{3}) receptors were able to significantly change the calcium wave patterns. Additionally, the distances between distinct IP_{3} receptor clusters and the pump strength were shown to control wave stability and instability [37, 38]. Sequestration properties through SERCA pumps [36] and mitochondrially controlled CICR [27] were further implicated in CICR. Related work has shown that calcium signals can modulate the shape of the cell nucleus [39, 40]. Thus, there seem to exist mechanisms in which calcium signals shape the geometry of organelles and the intracellular architecture shapes calcium signals.
In this study, we therefore chose to investigate how the relationship between ER and dendrite geometry can influence the dynamics of CICR waves. The ER is typically thought of as a continuous extension of the outer nuclear envelope that forms the rough ER surrounding the cell nucleus. As smooth ER, it extends into dendrites, all the way into synaptic spines [41]. Reconstructions of portions of smooth ER in dendrites have shown the ER to define a moderately branched, continuous structure [41, 42]. From an idealized standpoint, the ER has previously been considered as a neuronwithinaneuron [43–45], approximated in form of a cable that is embedded in a larger cable (the dendrite). We decided to follow this path of morphology simplification and designed perfectly symmetrical ER/dendrite cableincable model domains (see Fig. 1).
The associated calcium model consists of a diffusionreaction partial differential equation for the propagation of calcium in the cytosol and the reaction of free calcium with the buffer calbindin. Endoplasmic calcium is modeled by a simple diffusion equation. Membrane exchange mechanisms are included via flux boundary conditions at the plasma membrane and the ER membrane: Ca^{2+}ATPase pumps, Na^{+}/Ca^{2+} exchangers, as well as a leakage term on the plasma membrane; sarco/endoplasmic reticulum Ca^{2+}ATPase pumps, ryanodine receptor channels and a leakage flux on the ER membrane.
Our study focused on: (1) how calcium dynamics depend on the parameters dendrite radius, ER radius, and ryanodine receptor density, (2) whether thresholds can be found for stable calcium waves and (3) on gathering information about wave velocities. Simulation results show how cellular architecture, even in this simplified scenario, can determine the shape and fate of calcium waves and their stability. These data were further used to derive empirical laws that explain the behavior of the observed calcium dynamics in distinct geometric regimes. As a side note and in agreement with [26], we also show that the higher spatial dimension (as compared to classical 1D cable models) is indeed required to adequately represent even the most basic properties of calcium wave propagation.
Model and Methods
Model Domain
We conducted all simulations on a perfectly cylindrical model dendrite with a fixed length of 50 μm and variable radius, containing a cylindrically shaped ER of variable radius positioned exactly at the center of the dendrite; see Fig. 1. We used the rotational symmetry to reduce the problem to two dimensions. Calcium signals were induced by a calcium influx density through the left boundary of the dendrite with an initial strength of \(2.5\times 10^{18} \ \mathrm{mol}\ \upmu \mathrm{m}^{2}\ \mathrm{s}^{1}\) that linearly decreased to zero within 1 ms.
Model Equations
Spatiotemporal calcium dynamics in the intracellular space are modeled by a system of diffusionreaction equations described in the following. Transmembrane currents through channels and pumps are incorporated into these equations as flux boundary conditions (cf. Sect. 2.3).
Our model includes cytosolic/endoplasmic calcium and, in the cytosol, calbindinD_{28k} (CalB) as a calcium buffer. In the equations, we will represent cytosolic and endoplasmic calcium concentrations by the symbols \(c_{c}\), \(c_{e}\), and the (unbound) buffer concentration by b, respectively. Mobility in the cytosol/ER is described by the diffusion equation
where \(u(x,t)\) represents both calcium and calbindin as functions in space and time.
The buffering reaction between cytosolic calcium and calbindin is described by
The fact that CalB has four distinct highaffinity calciumbinding sites [46] is taken into account by quadrupling its concentration for the simulation (thus implicitly assuming identical and independent binding sites). The full domain equations for cytosolic calcium and calbindin are thus given by
in the cytosolic domain, where the concentration of the CalBCa^{2+} compound is expressed by the difference of the total concentration of CalB present in the cytosol (\(b^{ \mathrm{tot}}\)) and free CalB, the former of which is assumed to be constant in space and time (this amounts to the assumption that free and calciumbinding CalB have the same diffusive properties). All diffusion and reaction parameters are listed in Table 1.
Membrane Transport Mechanisms
In order to study calcium waves, we included ryanodine receptors for calcium release from the ER. We also added sarco/endoplasmic reticulum Ca^{2+}ATPase pumps (SERCA) for reuptake as well as a leakage flux on the ER membrane. Moreover, we defined fluxes across the plasma membrane by addition of Ca^{2+}ATPase pumps (PMCA), Na^{+}/Ca^{2+} exchangers (NCX) and a plasma membrane leakage term. This amounts to the ER and plasma membrane flux density equations
Here, \(j_{R}\) is the RyR flux density, \(j_{S}\) the SERCA flux density and \(j_{l,e}\) the leakage flux density on the ER membrane. The flux densities of PMCA, NCX and the leakage flux density of the plasma membrane are denoted by \(j_{P}\), \(j_{N}\) and \(j_{l,p}\), respectively. All flux densities are oriented towards the cytosol so that positive signs mean influx into the cytosol. Besides RyR, ERmembranebound IP_{3} receptor (IP_{3}R) channels are another source for intracellular calcium. They are known to require calcium and, more importantly, IP_{3} in the cytosol to become permeable to luminal calcium. IP_{3} is produced at the plasma membrane by receptor activation of phospholipase C. Compared to RyR calcium efflux, IP_{3}R calcium dynamics take place on a slower time scale (governed by the diffusion of IP_{3}). Preliminary simulations for the present study have shown that IP_{3}Rs may contribute to the initiation of a traveling wave, when small calcium influxes are not enough to trigger a wave through RyRs alone. The wave dynamics are, however, dominated by RyR. We therefore decided to exclude IP_{3}Rmediated calcium dynamics.
To incorporate existing singlechannel or singlepump models in the partial differential equations formulation, we assume channels and pumps to be continuously distributed along the membranes, using their membrane densities to calculate the flux densities \(j_{\mathrm{ERM}}\) and \(j_{\mathrm{PM}}\).
RyR Channels
The calcium flux density through ryanodine receptor channels in the ER membrane is given by an expression of the form
where \(\rho_{R}\) is the density of RyR in the ER membrane, \(p^{o}_{R}\) is the open state probability of a single channel, and \(I_{R}\) the singlechannel Ca^{2+} current. We describe the single channel ionic current by
where the reference current \(I_{R}^{\mathrm{ref}}\) is approximated from data presented in [51].
The open probability for RyR channels is modeled according to [50] and can be calculated as the sum of the two open states \(o_{1}\) and \(o_{2}\) in the system of ordinary differential equations
emerging from a fourstate Markov model with the kinetic constants \(k_{a}^{\pm }\), \(k_{b}^{\pm }\) and \(k_{c}^{\pm }\).
SERCA Pumps
The current from sarco/endoplasmic reticulum Ca^{2+}ATPase pumps is described by a model from [24], which was adapted for the threedimensional case, and gives rise to the Ca^{2+} flux density
The model reflects the dependence of the Ca^{2+} current not only on the cytosolic concentration, but also on the endoplasmic saturation. The density \(\rho_{\mathrm{S}}\) of SERCA pumps in the ER membrane was adapted to the RyR channel density in each simulation to ensure a zero net flux through the membrane in equilibrium conditions.
An interesting question with regard to the modeling of transmembrane currents is their equilibration at stationary, i.e., resting concentrations. The usual approach is to calibrate a leakage current such that it exactly counters the combined net current of all other involved mechanisms at equilibrium. At the ER membrane, this is only possible as long as the resulting direction of the leakage current points outward. With varying RyR densities, however, this is not necessarily the case (as leakage and RyR current have the same direction). Since the SERCA current is the only inward current through the ER membrane that we consider in our model, we chose to calibrate their density to ensure a zero net flux at equilibrium. This results in variable SERCA densities (larger when the RyR density is larger), which appears plausible from an energy consumption point of view: An increase in leakage currents to counter a reduction in RyR density would only lead to an increase in “futile cycling” through leakage and SERCA pumps. Another approach to avoid unneccessary consumption of energy is to use a bidirectional SERCA model (cf. [55], also for a more detailed view on “futile cycling”). Our simulations showed that the exact equilibration of the membrane currents has only a minuscule impact on the results presented here—the threshold values presented in the Results section only change by about 1% with equilibration by leakage flux and constant SERCA density (data not shown). The reason is that the propagation of the wave front takes place on a much faster time scale than the reuptake of calcium through SERCA pumps. The same applies to the extrusion by PMCA and NCX pumps, which are therefore modeled in a basic way.
NCX and PMCA Pumps
For the Na^{+}/Ca^{2+} exchanger current, our model assumes a constant Na^{+} concentration at the plasma membrane. The current densities for both plasma membrane transport mechanisms are expressed as first and secondorder Hill equations, respectively (cf. [53]):
Leakage
Both ER and plasma membranes allow a leakage flux not accounted for by the above transport mechanisms. The plasma membrane leakage flux is calibrated to ensure zero membrane net flux in the equilibrium state. Leakage flux densities are modeled by
where \(c_{o}\) is the extracellular Ca^{2+} concentration, which is assumed to be constant throughout all simulations. Values for all model parameters are gathered in Table 1.
Numerical Methods
For numerical simulations, the three equations are discretized in space using a finite volumes method. This allows for a natural integration of current densities across ER and plasma membranes into the reactiondiffusion process.
The system of ordinary differential equations (in time) arising from this procedure is nonlinear (due to the nonlinear reaction term and, more importantly, the highly nonlinear transport terms across the membranes). As these equations, notably the RyR channel dynamics, require a precise calibration of the time step size used for numerical solution, we employed a linearlyimplicit extrapolation (LIMEX) scheme to solve the nonlinear system, since LIMEX offers automated errorestimatebased control of integration order and step size [56, 57].
For the results we present here, the linear problems emerging in the LIMEX method were solved using a BiCGSTAB [58] linear solver preconditioned by a geometric multigrid method using Gauss–Seidel smoothing and SuperLU [59] as base solver. Computations were facilitated by a domain decomposition parallelization approach and carried out using the UG 4 framework [60] on the JURECA computer system at the Jülich Supercomputing Centre [61].
Implementation
All model components were implemented in a NeuroBox project. NeuroBox [62] is a neuroscientific toolbox that combines 1D, 2D and 3D modeling and simulation of electrical and biochemical signaling in a visual workflow environment. Visual workflows are created with VRLStudio [63] and the generalpurpose numerical framework UG 4 [60] is used to solve the set of coupled nonlinear partial differential equations.
Results
Range of Calcium Waves
In a first series of simulations, we examined how calcium dynamics depend on the ER radius and the density of RyR channels in the ER membrane. To that end, we created five model dendrite domains with a fixed dendrite radius of 0.2 μm each and an ER radius of 40 nm, 50 nm, 60 nm, 70 nm, and 80 nm, respectively. We then simulated calcium dynamics with RyR densities ranging from 0.5 \(\upmu \mathrm{m} ^{2}\) to 4.0 \(\upmu \mathrm{m}^{2}\) on all five domains and measured how far the initial calcium signal was able to travel through the dendrite.
The propagation distance was defined as the distance between the left end of the dendrite (the activation point) and the rightmost position where the open probability of the RyR channels exceeded 0.1 at any point in time. In our study we observe both stable and abortive calcium waves. The qualitative differences are illustrated in Fig. 2. For a stable wave (depicted for dendrite radius 0.4 μm, ER radius 0.15 μm, RyR density 2.5 \(\upmu \mathrm{m}^{2}\) in Fig. 2A), the traveling wavefront is clearly visible and maintains a constant shape. Reduction of the ER radius to a value of 0.11 μm makes the wave regime abortive (Fig. 2B): The wave front does not propagate further than 20 μm and the concentration at the wave front diminishes over time. While the velocity of the abortive wave decays to zero, it remains constant in the stable case (Fig. 2C)—justifying the assumption that the stable wave can propagate arbitrarily far. Snapshots of a sample simulation of a stable wave are depicted in Fig. 2D).
We plotted the propagation distance as a function of RyR density in Fig. 3. The traces for ER radii of 50 nm, 60 nm, 70 nm and 80 nm rise slowly at first, but then very rapidly approach a singularity at a RyR density that is specific to each domain. This threshold RyR density decreases with increasing ER radius and also exists for an ER radius of 40 nm, but is outside the plotted RyR density range in Fig. 3.
We draw the following conclusions from these observations for our 0.2 μm model dendrite:

(i)
A minimal RyR density is required in all five domains to elicit a calcium wave. The sharp gradient in the curves depicted in Fig. 3 show that once a critical RyR density is surpassed, calcium waves become stable. Otherwise the initial calcium influx causes only a very localized transient that does not propagate along the dendrite. This result is comparable to minimal distance criteria for the IP_{3} receptor clusters studied in [37].

(ii)
Dendrites with smaller ER need a higher RyR density to elicit a wave, indicating that the ability to trigger a calcium wave scales with the rate of calcium release from the ER. This seems to agree with [37], where wave stability is controlled by IP_{3}R pump strength and, as mentioned above, by IP_{3}R cluster distances. The critical threshold for wave stability decreases with increasing ER radius. Since dendrite diameters typically decrease the further they are away from the soma, ER diameters are forced to decrease as well. This would lead to a higher probability for calcium waves to become abortive. Figure 3 demonstrates how such dendritic tapering can be compensated for by an increase in RyR density.

(iii)
Conversely, as calcium release from the ER increases both with ER radius and RyR density, there is also a threshold ER radius for any sufficiently large RyR density above which calcium waves are stable. Thus, larger ER compartments in dendrites will enable neurons to more readily induce longdistance calcium signals.

(iv)
There are regimes in which abortive calcium waves are initiated (closely below threshold density). Such threshold dependency has previously been introduced in [26]. Our results confirm the existence of such thresholded regimes and are later quantified by deriving empirical laws. While the general behavior will be either a microdomain calcium event or a calcium wave, there is a critical intermediate region in which abortive waves with variable travel distances can be elicited.
In summary, the interplay between ER radius and RyR density enables dendrites to modulate the stability and reach of calcium signals, which may be critical in synaptic crosscommunication and synapsetonucleus signaling.
Threshold for Stable Calcium Waves
Having found the existence of such a threshold, we examined the prerequisites for stable calcium waves more closely. The guiding questions for this investigation were:

(1)
Do such thresholds also exist for larger dendrites?

(2)
If so, what is the relationship between dendrite radius and ER radius (and RyR density) at the threshold?
To address these questions, we varied the dendrite radius between 0.1 μm and 1.0 μm and the RyR channel density in the ER membrane between 0.5 \(\upmu \mathrm{m}^{2}\) and 4.0 \(\upmu \mathrm{m}^{2}\). For various combinations of these two parameters, we determined the threshold ER radius, above which a stable calcium wave can be elicited, using a bisection technique: We first located an interval in the ER radius domain for which the lower bound did not allow for a stable calcium wave, while the upper bound did. Then we simulated the calcium dynamics with ER radius taken from the middle of the interval and repeated the procedure with the lower half of the interval, if a stable wave was detected, with the upper half otherwise, until an adequate precision was reached. For this purpose, we accepted as a stable wave any signal, which, upon activation by calcium influx through the lefthand side of the model dendrite, propagated all the way (i.e., 50 μm) along the dendrite and was detected, in a measurement zone located at the righthand side (see also Fig. 2). With regard to the steep, singularitylike gradients of the wave range/RyR density functions for ranges between 20 μm and 50 μm in Fig. 3, we are confident that this is an appropriate choice. Since even if a setting existed where a wave had a range of 50 μm but became abortive beyond that range, a minor increase in RyR density (or ER radius) would then substantially increase its range.
With a look at the results in Fig. 4, the first question can be answered affirmatively. In fact, for all combinations of dendrite radius and RyR density that we checked, we found an ER radius threshold above which stable calcium waves were produced. This threshold ER radius decreases with increasing RyR density (indicating that the ability to sustain a stable calcium wave scales with the rate of calcium release from the ER) and increases with growing dendrite radius (because there is more space that the released calcium can diffuse into, which reduces the effective nearERmembrane concentration). It is interesting to point out three more aspects:

(i)
The traces in Fig. 4A (threshold ER radius plotted against RyR density) appear to reach a lower limit for large RyR densities, meaning the ER needs to have a minimal size to produce stable calcium waves. This makes sense as it has to contain a minimal amount of calcium that can actually be released during a wave event. This result is in agreement with [36], where an increase in SERCA density led to an increase in luminal calcium. This was demonstrated to have a positive feedback on the wave dynamics.

(ii)
The traces in Fig. 4B (threshold ER radius plotted against dendrite radius) reach an upper limit for large dendrite radii. In other words, there is a limit to the thickness of the ER required to elicit stable calcium waves and the ER does not need to be bigger than that in, e.g., very thick proximal dendrites. Figure 5 depicts the limit ER threshold values as a function of RyR density. The data is fitted by a model function of the form
$$\begin{aligned} r_{\mathrm{lim}}(\rho ) = \frac{ac+b\rho }{\rho  c} \end{aligned}$$with \(a = 0.377\ \upmu \mathrm{m}\), \(b = 0.0115 \ \upmu \mathrm{m}\) and \(c = 0.637\ \upmu \mathrm{m}^{2}\).

(iii)
The trace corresponding to the lowest tested RyR density (0.5 \(\upmu \mathrm{m}^{2}\)) in Fig. 4B does not reach such an upper limit (the trace corresponding to a RyR density of 1.0 \(\upmu \mathrm{m}^{2}\) does reach a limit, but outside the depicted range). One can conclude that a minimal RyR density is required for stable calcium waves (see also [26]) if the ER is restricted to specific dendritic volumes. Further investigation of this observation showed that this is due to calcium buffering in the cytosol: If the RyRsupported calcium current through the membrane is too small, too much of the released calcium is buffered and there is not enough left to trigger further release through RyR channels. This results is in agreement with [37], where similar conclusions were made for IP_{3}Rinduced waves. Reducing the buffer concentration by a factor of ten allowed for a limit ER threshold radius for all tested RyR densities (cf. Fig. 5).
Empirical Threshold Laws
To complete the investigation of our second question, we tried to describe the relationship between threshold ER radius and dendrite radius/RyR density by deriving appropriate empirical laws for the traces in Fig. 4. The fact that we found a threshold ER radius for wave elicitation in all combinations of dendrite radius and RyR density suggests that there exists a function \(r(R,\rho )\) (r for ER radius, R for dendrite radius, ρ for RyR density) describing a threshold ER radius manifold in the threedimensional parameter space spanned by r, R and ρ.
The spatiotemporal patterns of calcium waves can be roughly separated into three groups by means of the distance between ER and plasma membrane. The first group is defined by morphological configurations in which this distance is small enough for the radial calcium concentration profiles to be approximately constant, the second group contains the cases in which the distance is large enough for it not to influence the wave propagation, and the third group comprises all intermediate cases. In this section, we will derive separate empirical threshold laws for the former two groups. To motivate this distinction, we introduce a little thought experiment (cf. Fig. 6A).
In a stable wave regime, the axial position x, where a cytosolic threshold calcium concentration \(c^{\star }\) is attained near the ER membrane, needs to shift by Δx to the right within a certain time Δt (governed by the wave velocity). The amount of calcium necessary to increase the original concentration at \(x + \Delta x\) to \(c^{\star }\) within this time span needs to be released from the ER. It is proportional to the number of releasing RyR channels (see also [26]), which, in turn, is proportional to: (i) the RyR density ρ and (ii) the radius r of the ER membrane, both assumed constant in time. The singlechannel current, however, is timedependent in that it depends on locally available endoplasmic calcium (see Eq. (8)), which decays rapidly during RyR release. While the amount of released calcium is proportional to r, the endoplasmic volume is proportional to \(r^{2}\). Assuming rapid radial equilibration within the ER, the decay of a singlechannel current is proportional to \(r^{1}\). Integrating the decreasing current from t to \(t+\Delta t\) shows that the amount of released calcium is also proportional to (iii) the factor \((1b/r)\) with some constant b, in a firstorder approximation. Putting (i) to (iii) together, an amount proportional to
is released into the cytosol, This is where we distinguish between regimes with small and large distances between plasma and ER membrane.
Small Distance Between Plasma and ER Membrane
When the distance between ER and plasma membrane is small, then there is only a negligible radial concentration gradient at all times, as diffusion will almost instantaneously produce constant radial concentration profiles (cf. Fig. 6B). With this in mind, the increase in cytosolic calcium concentration, caused by the amount of calcium released from the ER, is inversely proportional to \(R^{2}r^{2}\), assuming that a fixed percentage of the released calcium is bound by the mobile buffer, which leaves proportionality intact. As the required increase to bridge the gap between the concentrations at x and \(x+\Delta x\) is fix, a constant \(\frac{1}{2c}\) can be introduced to arrive at the relation
with the parameters b and c, usable in a least squares fitting, which very closely fits the 0.1 μm dendrite radius trace in Fig. 4A. The one for the 0.2 μm dendrite radius trace is also acceptable, but using R as an additional parameter in the fitting expression produces substantially better fitting results for the 0.3 μm and 0.4 μm dendrite radius traces, indicating that the model assumption of radially constant calcium concentrations is already largely violated in theses cases.
Large Distance Between Plasma and ER Membrane
When the distance between ER and plasma membrane is large, the wave front moves faster in the axial direction than its radial component reaches the plasma membrane, and the wave does not “sense” the presence of the plasma membrane (cf. Fig. 6C). The fact that the upper traces in Fig. 4A converge towards a limit trace and that the lower traces in Fig. 4B attain a limit is due to this circumstance. In this situation, radial calcium diffusion will be limited to a distance a from the ER membrane until the threshold concentration \(c^{\star }\) is reached at axial position \(x+\Delta x\). It is reasonable to assume that the calcium concentration near the ER membrane after release through RyR channels is inversely proportional to the volume it diffuses into, i.e., \((r+a)^{2}  r^{2}\). Like above, with a fixed percentage of the released calcium being bound by the buffer, we can introduce a constant \(\frac{c}{4a}\) (which depends on the buffer concentration) to arrive at the equation
The parameters a, b, and c can be determined to fit the limit threshold ER radii attained for large dendrite radii (Fig. 5, blue trace), also for buffer concentrations reduced by a factor of ten, in which case the parameter c is significantly smaller (Fig. 5, orange trace). Note that the buffer concentration puts a constraint on the RyR density: If the density falls below a threshold (the constant c), then there can be no calcium wave, regardless of how big the ER is.
Intermediate Distances
We were not able to find descriptive model functions to fit the threshold values for intermediate distances. In such regimes, the significant radial concentration gradient at the wave front, which varies strongly with the distance between the two membranes, makes it difficult to find expressions similar to the two presented ones.
Wave Velocity
During the described simulations, we observed that the velocity of calcium waves varies with ER radius and RyR density. This led us to finally examine the velocity of stable waves. Following the initial activation, all stable waves reached a constant velocity shortly after their formation (cf. Fig. 2C)). These constant velocities were measured using the rightmost axial coordinate where the RyR open probability exceeded 0.1 as the wave front position in every time step of the simulation. Results (shown for a dendrite radius of 0.2 μm in Fig. 7) indicate that calcium wave velocity increases both with increasing ER radius and RyR density. This can be expected given that the driving force for the movement of the wave is diffusion, which becomes faster when axial concentration gradients become bigger, which, in turn, is the case when the efflux density from the ER increases. The traces for wave velocity as a function of RyR density rise more rapidly in a small range above the threshold for stable waves and then exhibit a quasilinear behavior with a slope seemingly independent of the ER radius. Velocity traces for a larger dendrite of radius 1.0 μm show the same basic behavior (data not shown), though with a less pronounced dependency on the ER radius, which is in line with our findings for the ER radius threshold in the previous section.
Discussion
Our computational study finds that calcium waves can be triggered in an idealized morphological setting. If certain geometric and physiological conditions are met, wave propagation is a stable process, similar to the allornothing properties of electrical action potentials. A minimal ryanodine receptor density is required to guarantee sufficient calcium release from the ER and this density scales inversely with ER surface. We also observed abortive regimes in small parameter regions below thresholds for ER size and ryanodine receptor densities. When considering realistic morphologies, spatial variations in dendritic and ER morphology as well as changes in RyR density are likely to occur, thus, abortive regimes may be important when it comes to synaptic crosscommunication with respect to calcium signals.
We used spatially constant RyR densities and membrane radii throughout the study, since we were interested in obtaining functions expressing their relationships with respect to the potential for the elicitation of stable calcium waves. Of course, in real dendrites neither dendrite radius nor RyR density are spatially constant [64], nor can the ER be expected to be a perfectly shaped cylinder that is centrally positioned within the dendrite. Future research, which is outside the scope of this study, could integrate heterogeneous membrane composition and detailed threedimensional morphologies.
In the context of spiketimingdependent plasticity, calcium wave velocity may be a controlling parameter. Our study shows that calcium wave velocities vary depending on ER size and RyR density. This merits a closer investigation of the timing of calcium signals in connection to, e.g., backpropagating action potentials.
A more precise understanding of the intracellular architecture may be critical in wave timing and stability. For example, mitochondria are an important sink for excess cytosolic calcium. They can quickly absorb considerable amounts of calcium through mitochondrial calcium uniporter channels [65] and extend through dendrites as filamentous networks [41, 42]. It would therefore be interesting for future work to extend this study to a fully threedimensional, nonsymmetric case, to investigate the interplay between mitochondria and ER in the context of dendritic calcium waves.
In addition to the stability and velocity of dendritic calcium waves, changes in the cellular and intracellular architecture may promote direction selectivity of calcium signals. Since the radius of dendrites is variable, typically large in proximity of the soma and smaller in more distal regions, we wondered whether there could be some kind of direction selectivity in the propagation of calcium waves. Due to the increasing amount of releasable calcium in the growing ER, we suspected that stable waves might be better supported when propagating from thin dendrites to larger ones with larger ER than in the direction of decreasing radii. We tested this hypothesis with three simulations (data not shown), the first on a mediumsize dendrite with constant radius, the second on the same dendrite with increasing radius, and the last one on the same dendrite with decreasing radius. In all cases, the ER radius at each location was chosen to be just below the threshold for stable calcium waves at the chosen (constant) RyR density, so that one would expect the waves to terminate at some point. This indeed happened in all cases, however, the distance the waves were able to travel differed: shortest in the shrinking dendrite and longest in the expanding dendrite. The results confirmed our hypothesis. Yet, the effect was relatively small (approx. 11 μm difference in travel distance between smalltolarge and largetosmall scenarios) and required very precise calibration of the ER radius (too small and all waves would have terminated even earlier, too large and all waves would have been stable). We therefore believe that this effect is unlikely to be relevant under realistic conditions.
Our study further highlights the question of model dimensionality. While in many cases, in which electrical properties of neurons are studied, onedimensional multicompartment models are employed, this may not always be possible for intracellular biochemical modeling. The fact that we could not find a simple descriptive model for regimes with neither very small nor very large distance between ER and plasma membrane reflects the nontrivial dynamics of the calcium signal with respect to space. Even in this relatively simple model setup with rotational symmetry in the geometry, it is not generally possible to reduce dimensionality (by ignoring the radial concentration profile in the cytosol) to a model that is essentially 1D. A fortiori, such a model would not be able to generate the correct stability thresholds and wave velocities in mediumsize and large dendrites. For very thin dendrites, however, a simplification of 3D calcium models to 1D is conceivable provided the symmetry prerequisites of this study are met.
Abbreviations
 CalB:

calbindinD_{28k}
 CICR:

calciuminduced calcium release
 ER:

endoplasmic reticulum
 IP_{3} :

inositol 1,4,5trisphosphate
 IP_{3}R:

IP_{3} receptor
 LIMEX:

linearlyimplicit extrapolation
 NCX:

Na^{+}/Ca^{2+} exchanger
 PMCA:

plasma membrane Ca^{2+}ATPase
 RyR:

ryanodine receptor
 SERCA:

sarco/endoplasmic reticulum Ca^{2+}ATPase
References
 1.
Sheng M, McFadden G, Greenberg ME. Membrane depolarization and calcium induce cfos transcription via phosphorylation of transcription factor CREB. Neuron. 1990;4(4):571–82.
 2.
Sun P, Enslen H, Myung PS, Maurer RA. Differential activation of CREB by Ca^{2+}/calmodulindependent protein kinases type II and type IV involves phosphorylation of a site that negatively regulates activity. Genes Dev. 1994;8(21):2527–39.
 3.
Tao X, Finkbeiner S, Arnold DB, Shaywitz AJ, Greenberg ME. Ca^{2+} influx regulates BDNF transcription by a CREB family transcription factordependent mechanism. Neuron. 1998;20(4):709–26.
 4.
Shieh PB, Hu SC, Bobb K, Timmusk T, Ghosh A. Identification of a signaling pathway involved in calcium regulation of BDNF expression. Neuron. 1998;20(4):727–40.
 5.
LimbäckStokin K, Korzus E, NagaokaYasuda R, Mayford M. Nuclear calcium/calmodulin regulates memory consolidation. J Neurosci. 2004;24(48):10858–67.
 6.
Papadia S, Stevenson P, Hardingham NR, Bading H, Hardingham GE. Nuclear Ca^{2+} and the cAMP response elementbinding protein family mediate a late phase of activitydependent neuroprotection. J Neurosci. 2005;25(17):4279–87.
 7.
Zhang SJ, Zou LL, Lau D, Ditzel DA, DelucingeVivier C, Aso Y, Descombes P, Bading H. Nuclear calcium signaling controls expression of a large gene pool: identification of a gene program for acquired neuroprotection induced by synaptic activity. PLoS Genet. 2009;5(8):e1000604. https://doi.org/10.1371/journal.pgen.1000604.
 8.
Zhang SJ, Buchthal B, Lau D, Hayer S, Dick O, Schwaninger M, Veltkamp R, Zou M, Weiss U, Bading H. A signaling cascade of nuclear calciumCREBATF3 activated by synaptic NMDA receptors defines a gene repression module that protects against extrasynaptic NMDA receptorinduced neuronal cell death and ischemic brain damage. J Neurosci. 2011;31(13):4978–90.
 9.
Brini M, Carafoli E. The plasma membrane Ca^{2+} ATPase and the plasma membrane sodium calcium exchanger cooperate in the regulation of cell calcium. Cold Spring Harb Perspect Biol. 2011;3(2):a004168. https://doi.org/10.1101/cshperspect.a004168.
 10.
Görlach A, Klappa P, Kietzmann T. The endoplasmic reticulum: folding, calcium homeostasis, signaling, and redox control. Antioxid Redox Signal. 2006;8(9–10):1391–418.
 11.
Berridge MJ. Calcium microdomains: organization and function. Cell Calcium. 2006;40(5–6):405–12.
 12.
Foskett JK, White C, Cheung KH, Mak DO. Inositol trisphosphate receptor Ca^{2+} release channels. Physiol Rev. 2007;87(2):593–658.
 13.
Sutko JL, Airey JA. Ryanodine receptor Ca^{2+} release channels: does diversity in form equal diversity in function? Physiol Rev. 1996;76(4):1027–71.
 14.
Rizzuto R, Brini M, Murgia M, Pozzan T. Microdomains with high Ca^{2+} close to IP_{3}sensitive channels that are sensed by neighboring mitochondria. Science. 1993;262(5134):744–7.
 15.
Walsh C, Barrow S, Voronina S, Chvanov M, Petersen OH, Tepikin A. Modulation of calcium signalling by mitochondria. Biochim Biophys Acta. 2009;1787(11):1374–82.
 16.
Goetz JG, Genty H, StPierre P, Dang T, Joshi B, Sauvé R, Vogl W, Nabi IR. Reversible interactions between smooth domains of the endoplasmic reticulum and mitochondria are regulated by physiological cytosolic Ca^{2+} levels. J Cell Sci. 2007;120(20):3553–64.
 17.
Bezprozvanny I, Hayden MR. Deranged neuronal calcium signaling and Huntington’s disease. Biochem Biophys Res Commun. 2004;322:1310–7.
 18.
Nimmrich V, Grimm C, Draguhn A, Barghorn S, Lehmann A, Schoemaker H, Hillen H, Gross G, Ebert U, Bruehl C. Amyloid β oligomers (A \(\beta_{1{}42}\) globulomer) suppress spontaneous synaptic activity by inhibition of P/Qtype calcium currents. J Neurosci. 2008;28(4):788–97.
 19.
Ito E, Oka K, Etcheberrigaray R, Nelson TJ, McPhie DL, TofelGrehl B, Gibson GE, Alkon DL. Internal Ca^{2+} mobilization is altered in fibroblasts from patients with Alzheimer disease. Proc Natl Acad Sci USA. 1994;91(2):534–8.
 20.
Yoo AS, Cheng I, Chung S, Grenfell TZ, Lee H, PackChung E, Handler M, Shen J, Xia W, Tesco G, Saunders AJ, Ding K, Frosch MP, Tanzi RE, Kim TW. Presenilinmediated modulation of capacitative calcium entry. Neuron. 2000;27(3):561–72.
 21.
Stutzmann GE, Smith I, Caccamo A, Oddo S, Laferla FM, Parker I. Enhanced ryanodine receptor recruitment contributes to Ca^{2+} disruptions in young, adult, and aged Alzheimer’s disease mice. J Neurosci. 2006;26(19):5180–9.
 22.
Cheung KH, Shineman D, Müller M, Cárdenas C, Mei L, Yang J, Tomita T, Iwatsubo T, Lee VMY, Foskett JK. Mechanism of Ca^{2+} disruption in Alzheimer’s disease by presenilin regulation of InsP_{3} receptor channel gating. Neuron. 2008;58(6):871–83. https://doi.org/10.1016/j.neuron.2008.04.015.
 23.
Green KN, Demuro A, Akbari Y, Hitt BD, Smith IF, Parker I, LaFerla FM. SERCA pump activity is physiologically regulated by presenilin and regulates amyloid β production. J Cell Biol. 2008;181(7):1107–16.
 24.
Sneyd J, TsanevaAtanasova K, Bruce JIE, Straub SV, Giovannucci DR, Yule DI. A model of calcium waves in pancreatic and parotid acinar cells. Biophys J. 2003;85(3):1392–405. https://doi.org/10.1016/S00063495(03)74572X.
 25.
Ross WN. Understanding calcium waves and sparks in central neurons. Nat Rev Neurosci. 2012;13:157–68. https://doi.org/10.1038/nrn3168.
 26.
Thul R, Rietdorf K, Bootman MD, Coombes S. Unifying principles of calcium wave propagation—insights from a threedimensional model for atrial myocytes. Biochim Biophys Acta, Mol Cell Res. 2015;1853(9):2131–43. https://doi.org/10.1016/j.bbamcr.2015.02.019. 13th European Symposium on Calcium.
 27.
Falcke M. Deterministic and stochastic models of intracellular Ca^{2+} waves. New J Phys. 2003;5:96.1–96.28.
 28.
Kaźmierczak B, Peradzyński Z. Calcium waves with fast buffers and mechanical effects. J Math Biol. 2011;62(1):1–38. https://doi.org/10.1007/s0028500903232.
 29.
Dani JW, Chernjavsky A, Smith SJ. Neuronal activity triggers calcium waves in hippocampal astrocyte networks. Neuron. 1992;8:429–40.
 30.
Smith N. Generation of calcium waves in living cells by pulsedlaserinduced photodisruption. Appl Phys Lett. 2001;79:1208–10. https://doi.org/10.1063/1.1397255.
 31.
Yang X. Computational modelling of nonlinear calcium waves. Appl Math Model. 2006;30:200–8.
 32.
Thul R, Smith GD, Coombes S. A bidomain threshold model of propagating calcium waves. J Math Biol. 2008;56:435–63.
 33.
Timofeeva Y, Coombes S. Wave bifurcation and propagation failure in a model of Ca^{2+} release. J Math Biol. 2003;47:249–69.
 34.
Petrovic P, Valent I, Cocherova E, Pavelkova J, Zahradnikova A. Ryanodine receptor gating controls generation of diastolic calcium waves in cardiac myocytes. J Gen Physiol. 2015;145:489–511.
 35.
Thurley K, Skupin A, Thul R, Falcke M. Fundamental properties of Ca^{2+} signals. Biochem Biophys Acta. 2012;1820:1185–94. https://doi.org/10.1016/j.bbagen.2011.10.007.
 36.
Falcke M, Li Y, Lechleiter JD, Camacho P. Modeling the dependence of the period of intracellular Ca^{2+} waves on SERCA expression. Biophys J. 2003;85:1474–81. https://doi.org/10.1016/S00063495(03)745809.
 37.
Guisoni N, Ferrero P, Layana C, Diambra L. Abortive and propagating intracellular calcium waves: analysis from a hybrid model. PLoS ONE. 2015;10:1–15. https://doi.org/10.1371/journal.pone.0115187.
 38.
Thul R, Falcke M. Stability of membrane bound reactions. Phys Rev Lett. 2004;93:188103. https://doi.org/10.1103/PhysRevLett.93.188103.
 39.
Wittmann M, Queisser G, Eder A, Wiegert JS, Bengtson CP, Hellwig A, Wittum G, Bading H. Synaptic activity induces dramatic changes in the geometry of the cell nucleus: interplay between nuclear structure, histone H3 phosphorylation, and nuclear calcium signaling. J Neurosci. 2009;29(47):14687–700.
 40.
Queisser G, Wiegert JS, Bading H. Structural dynamics of the cell nucleus: basis for morphology modulation of nuclear calcium signaling and gene transcription. Nucleus. 2011;2(2):98–104.
 41.
Wu Y, Whiteus C, Xu CS, Hayworth KJ, Weinberg RJ, Hess HF, De Camilli P. Contacts between the endoplasmic reticulum and other membranes in neurons. Proc Natl Acad Sci. 2017;114(24):4859–67. https://doi.org/10.1073/pnas.1701078114. http://www.pnas.org/content/114/24/E4859.full.pdf.
 42.
Popov V, Medvedev NI, Davies HA, Stewart MG. Mitochondria form a filamentous reticular network in hippocampal dendrites but are present as discrete bodies in axons: a threedimensional ultrastructural study. J Comp Neurol. 2005;492(1):50–65. https://doi.org/10.1002/cne.20682.
 43.
Terasaki M, Slater NT, Fein A, Schmidek A, Reese TS. Continuous network of endoplasmic reticulum in cerebellar Purkinje neurons. Proc Natl Acad Sci. 1994;91(16):7510–4. http://www.pnas.org/content/91/16/7510.full.pdf.
 44.
Berridge MJ. Neuronal calcium signaling. Neuron. 1998;21(1):13–26.
 45.
Shemer I, Brinne B, Tegnér J, Grillner S. Electrotonic signals along intracellular membranes may interconnect dendritic spines and nucleus. PLoS Comput Biol. 2008;4(3):1–19. https://doi.org/10.1371/journal.pcbi.1000036.
 46.
Veenstra TD, Johnson KL, Tomlinson AJ, Naylor S, Kumar R. Determination of calciumbinding sites in rat brain calbindin D_{28k} by electrospray ionization mass spectrometry. Biochemistry. 1997;36(12):3535–42. https://doi.org/10.1021/bi9628329. PMID: 9132004.
 47.
Müller A, Kukley M, Stausberg P, Beck H, Müller W, Dietrich D. Endogenous Ca^{2+} buffer concentration and Ca^{2+} microdomains in hippocampal neurons. J Neurosci. 2005;25(3):558–65. https://doi.org/10.1523/JNEUROSCI.379904.2005. http://www.jneurosci.org/content/25/3/558.full.pdf+html.
 48.
Allbritton NL, Meyer T, Stryer L. Range of messenger action of calcium ion and inositol 1,4,5trisphosphate. Science. 1992;258(5089):1812–5. https://doi.org/10.1126/science.1465619. http://www.sciencemag.org/content/258/5089/1812.full.pdf.
 49.
Schmidt H, Stiefel KM, Racay P, Schwaller B, Eilers J. Mutational analysis of dendritic Ca^{2+} kinetics in rodent Purkinje cells: role of parvalbumin and calbindin D_{28k}. J Physiol. 2003;551(1):13–32. https://doi.org/10.1113/jphysiol.2002.035824. http://jp.physoc.org/content/551/1/13.full.pdf+html.
 50.
Keizer J, Levine L. Ryanodine receptor adaptation and Ca^{2+}induced Ca^{2+} releasedependent Ca^{2+} oscillations. Biophys J. 1996;71(6):3477–87. https://doi.org/10.1016/S00063495(96)795437. http://download.cell.com/biophysj/pdf/PIIS0006349596795437.pdf.
 51.
Tinker A, Lindsay ARG, Williams AJ. Cation conduction in the calcium release channel of the cardiac sarcoplasmic reticulum under physiological and pathophysiological conditions. Cardiovasc Res. 1993;27(10):1820–5. https://doi.org/10.1093/cvr/27.10.1820. http://cardiovascres.oxfordjournals.org/content/27/10/1820.full.pdf.
 52.
Chiu VCK, Haynes DH. Rapid kinetic studies of active Ca^{2+} transport in sarcoplasmic reticulum. J Membr Biol. 1980;56:219–39. https://doi.org/10.1007/BF01869478.
 53.
Graupner M. A theory of plasma membrane calcium pump function and its consequences for presynaptic calcium dynamics. Diploma thesis. Technische Universität Dresden; 2003.
 54.
Elwess NL, Filoteo AG, Enyedi A, Penniston JT. Plasma membrane Ca^{2+} pump isoforms 2a and 2b are unusually responsive to calmodulin and Ca^{2+}. J Biol Chem. 1997;272(29):17981–6. https://doi.org/10.1074/jbc.272.29.17981. http://www.jbc.org/content/272/29/17981.full.pdf+html.
 55.
Higgins ER, Cannell MB, Sneyd J. A buffering SERCA pump in models of calcium dynamics. Biophys J. 2006;91(1):151–63. https://doi.org/10.1529/biophysj.105.075747.
 56.
Deuflhard P, Hairer E, Zugck J. Onestep and extrapolation methods for differentialalgebraic systems. Numer Math. 1987;51(5):501–16. https://doi.org/10.1007/BF01400352.
 57.
Deuflhard P. Order and stepsize control in extrapolation methods. Numer Math. 1983;41(3):399–422. https://doi.org/10.1007/BF01418332.
 58.
van der Vorst HA. BiCGSTAB: a fast and smoothly converging variant of BiCG for the solution of nonsymmetric linear systems. SIAM J Sci Stat Comput. 1992;13(2):631–44. https://doi.org/10.1137/0913035. http://epubs.siam.org/doi/abs/10.1137/0913035.
 59.
Demmel JW, Eisenstat SC, Gilbert JR, Li XS, Liu JWH. A supernodal approach to sparse partial pivoting. SIAM J Matrix Anal Appl. 1999;20(3):720–55.
 60.
Vogel A, Reiter S, Rupp M, Nägel A, Wittum G. UG 4: a novel flexible software system for simulating PDE based models on high performance computers. Comput Vis Sci. 2013;16(4):165–79. https://doi.org/10.1007/s0079101402329.
 61.
Jülich Supercomputing Centre. JURECA: generalpurpose supercomputer at Jülich Supercomputing Centre. J LargeScale Res Facil. 2016;2:A62. https://doi.org/10.17815/jlsrf2121.
 62.
Breit M, Stepniewski M, Grein S, Gottmann P, Reinhardt L, Queisser G. Anatomically detailed and largescale simulations studying synapse loss and synchrony using NeuroBox. Front Neuroanat. 2016;10:8. https://doi.org/10.3389/fnana.2016.00008.
 63.
Hoffer M, Poliwoda C, Wittum G. Visual reflection library: a framework for declarative GUI programming on the Java platform. Comput Vis Sci. 2013;16(4):181–92. https://doi.org/10.1007/s007910140230y.
 64.
Hertle DN, Yeckel MF. Distribution of inositol1,4,5trisphosphate receptor isotypes and ryanodine receptor isotypes during maturation of the rat hippocampus. Neuroscience. 2007;150(3):625–38. https://doi.org/10.1016/j.neuroscience.2007.09.058.
 65.
Kirichok Y, Krapivinsky G, Clapham DE. The mitochondrial calcium uniporter is a highly selective ion channel. Nature. 2004;427:360–4. https://doi.org/10.1038/nature02246.
Acknowledgements
We wish to thank Arne Nägel for the implementation of the LIMEX scheme in UG 4 and for his advice on how to use it. We gratefully acknowledge the computing time granted by the John von Neumann Institute for Computing (NIC) and provided on the supercomputer JURECA at Jülich Supercomputing Centre (JSC).
Availability of data and materials
Please contact the authors for data requests.
Funding
This research was funded by BMBF (Collaborative Research in Computational Neuroscience 01GQ1410B to GQ).
Author information
Affiliations
Contributions
MB and GQ formulated the scientific questions and designed the simulations. The simulations were executed and analyzed by MB. The research project was supervised by GQ. MB and GQ wrote the manuscript. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Not applicabale.
Competing interests
The authors declare that they have no competing interests.
Consent for publication
Not applicabale.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Electronic Supplementary Material
Below are the links to the electronic supplementary material.
13408_2018_64_MOESM1_ESM.mkv
Video of wave propagation in thin dendrite: Typical wave in a dendrite of radius 0.2 μm, with an ER radius of 50 nm and a ryanodine receptor density in the ER membrane of 3.5 \(\upmu \mathrm{m}^{2}\). One second in the video corresponds to 10 ms of simulated time. Calcium concentrations are colorcoded blue (\(c _{c}\leq 50\ \mathrm{n}\mathrm{M}\), \(c_{e}\leq 25 \ \upmu \mathrm{M} \)) to red (\(c_{c}\geq 1 \ \upmu \mathrm{M} \), \(c_{e}\geq 250\ \upmu \mathrm{M} \)). Radial coordinate scaled by a factor of 40 to enhance visibility (MKV 74 kB)
13408_2018_64_MOESM2_ESM.mkv
Video of wave propagation in thick dendrite: Typical wave in a dendrite of radius 1.0 μm, with an ER radius of 0.2 μm and a ryanodine receptor density in the ER membrane of 2.5 \(\upmu \mathrm{m}^{2}\). One second in the video corresponds to 10 ms of simulated time. Calcium concentrations are colorcoded blue (\(c_{c}\leq 50 \ \mathrm{n}\mathrm{M}\), \(c_{e}\leq 25 \ \upmu \mathrm{M} \)) to red (\(c_{c}\geq 1\ \upmu \mathrm{M} \), \(c_{e}\geq 250\ \upmu \mathrm{M}\)). Radial coordinate scaled by a factor of 8 to enhance visibility (MKV 88 kB)
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.
About this article
Cite this article
Breit, M., Queisser, G. What Is Required for Neuronal Calcium Waves? A Numerical Parameter Study. J. Math. Neurosc. 8, 9 (2018). https://doi.org/10.1186/s134080180064x
Received:
Accepted:
Published:
Keywords
 Calcium waves
 Endoplasmic reticulum
 Ryanodine receptors
 3D modeling
 Structurefunction interplay
 Numerical simulation