Research Article  Open Access
Multiscale Modeling Indicates That Temperature Dependent [Ca^{2+}]_{i} Spiking in Astrocytes Is Quantitatively Consistent with Modulated SERCA Activity
Abstract
Changes in the cytosolic Ca^{2+} concentration () are the most predominant active signaling mechanism in astrocytes that can modulate neuronal activity and is assumed to influence neuronal plasticity. Although Ca^{2+} signaling in astrocytes has been intensively studied in the past, our understanding of the signaling mechanism and its impact on tissue level is still incomplete. Here we revisit our previously published data on the strong temperature dependence of Ca^{2+} signals in both cultured primary astrocytes and astrocytes in acute brain slices of mice. We apply multiscale modeling to test the hypothesis that the temperature dependent spiking is mainly caused by the increased activity of the sarcoendoplasmic reticulum ATPases (SERCAs) that remove Ca^{2+} from the cytosol into the endoplasmic reticulum. Quantitative comparison of experimental data with multiscale simulations supports the SERCA activity hypothesis. Further analysis of multiscale modeling and traditional rate equations indicates that the experimental observations are a spatial phenomenon where increasing pump strength leads to a decoupling of Ca^{2+} release sites and subsequently to vanishing spikes.
1. Introduction
The human brain is by far the most complex organ; it hosts our consciences, is the source of our inspiration, and controls many of our vegetative functions. The extraordinary complexity of cellular interactions in terms of structural as well as dynamic properties enables the huge flexibility of our behavior and the amazing plasticity required for learning [1]. While traditional neuroscience has focused mainly on neurons for decades, the importance of the more abundant glia cells has only become evident more recently [2, 3].
Glia cells are generally known for their supportive function for neurons. The minority of glia cells are oligodendrocytes that myelinate axons [4]. More recently, oligodendrocytes were also found to essentially support axon integrity by metabolic coupling through monocarboxylate transporters (MCTs) [5, 6]. Microglia represent the macrophages of the brain and play an important role in immune responses as well as in modifying the neuronal network structure by synapse pruning [7]. The majority of glia cells are astrocytes that are commonly considered as metabolic supporters of neurons. On the one hand astrocytes regulate the blood flow in response to neuronal activity [8]. Astrocytes take up approximately 80% of the metabolized glucose from the blood and perform glycolysis [9]. According to the astrocyteneuron lactate shuttle (ANLS) hypothesis [10], astrocytes generate lactate from pyruvate, the end product of glycolysis, by the lactate dehydrogenase (LDH) and export it into the extracellular space by MCTs. From the extracellular space, neurons take up the lactate and generate pyruvate by the reversible LDH reaction. The neuronal pyruvate is subsequently used as energy substrate for the TCA cycle in mitochondria to generate adenosine triphosphate (ATP) that is extensively needed to run the ion pumps in neurons to keep the large electrochemical potential across the plasma membrane. On the other hand, astrocytes play also an important function for neuronal proteostasis as they remove degraded material from the tissue either directly by intercellular transport towards the blood stream or by controlling the cerebral fluid flux [11].
Besides their supportive role in metabolism, astrocytes are now also known to modulate neuronal activity through active signaling mechanisms. Astrocytes are sealing up the dyadic cleft facilitating neurotransmitter diffusion and are essential for glutamate uptake and recycling [12]. Since astrocytes express also glutamate receptors, they can be stimulated by neuronal activity. Upon stimulation, astrocytes can secrete glutamate as well as ATP into the extracellular space including the dyadic cleft and thereby induce a local amplification mechanism. This neuronastrocyte crosstalk has led to the picture of the tripartite synapse [13, 14].
The most predominant and studied signaling mechanism in astrocytes is inositol 1,4,5trisphosphate (IP_{3}) mediated Ca^{2+} signaling. In general, Ca^{2+} signaling is a versatile and universal second messenger that translates extracellular signals into intracellular responses like molecular motor activation or modulation of gene expression [15–17]. Within the IP_{3} pathway, extracellular agonists are detected by plasma membrane receptors. Upon binding, these Gprotein coupled receptors activate the phospholipase C (PLC) that produces IP_{3} at the plasma membrane from which it diffuses into the cytosol. There IP_{3} can bind to IP_{3} receptors (IP_{3}Rs) localized in the membrane of the endoplasmic reticulum (ER). When IP_{3} and Ca^{2+} are bound to the receptor, it can be open and Ca^{2+} is released from the ER into the cytosol due to the large concentration differences between these two compartments of typically several orders of magnitude. From the cytosol the released Ca^{2+} is moved back into the extracellular space by plasma membrane calcium ATPase (PMCA) and into the ER by sacroendoplasmic reticulum calcium ATPases (SERCAs).
The collection of the Ca^{2+} tool kit [18] can generate a wide spectrum of cellular responses ranging from long lasting high cytosolic Ca^{2+} levels () to single transients. The most common Ca^{2+} signals in astrocytes are repetitive cytosolic transients often referred to as oscillations [18] that are typically elicited by neurotransmitters and hormones but astrocytes do also exhibit spontaneous spiking.
This oscillationlike behavior is rather common for many cell and stimulation types and has led to the general assumption that Ca^{2+} signaling is frequency coded [17]. This perspective is challenged by more recent studies including ours that have demonstrated that spiking is indeed stochastic and therefore not truly oscillatory [19–21]. By systematic analysis of single cell time courses for hundreds of cells of different cell types and stimulation antagonists, we found a preserved linear relation between the average period () and the standard deviation () of the interspike intervals. Applying methods from the theory of stochastic processes, we have shown that this pathway specific linear  relation characterizes the internal feedback strength [22] and represents an estimator for the information content [20, 23].
To understand the random signaling mechanism and underlying principles, we developed a computational multiscale model that takes into account the compartmentalization of the cell into cytosol and ER, the stochastic binding of signaling molecules to individual IP_{3}Rs, their spatial organization into channel clusters, and the spatial diffusion of Ca^{2+} and Ca^{2+} binding buffers [27]. The enabled mechanistic simulations have shown that the experimental observations are in accordance with the implemented hierarchical stochastic system where single channel fluctuations are carried onto the level of the cell by diffusion mediated Ca^{2+} induced Ca^{2+} release (CICR). The modeling based insights on the functional robustness have guided further experiments that have revealed a fold change behavior in Ca^{2+} signaling of stimulated HEK cells and hepatocytes and were used to develop and validate a universal encoding relation [28].
The establishment of the signaling mechanism was driven to a large extent by experiments with astrocytes and microglia [21]. In particular the spontaneous activity of cultured astrocytes was exploited to obtain a baseline activity of the Ca^{2+} signaling machinery without potential antagonist induced feedback mechanisms. The astrocytic base level activity is characterized by a slope of 1 in the linear  relation and corresponds to the theoretical minimal information gain. This is in accordance with the general assumption that stimulation should carry information that would correspond to a decreased slope of the  relation and is observed in stimulated cells [20].
The spontaneous Ca^{2+} transients in astrocytes occur independently of neuronal activity but are modulated by neuronal function like induced epileptic activity [29]. The astrocytic Ca^{2+} increase can in turn be sensed by neurons because it has been correlated with release of neuroactive substances [30, 31]. Thereby astrocytes can release glutamate [32, 33], DSerine [34, 35], ATP [36–39], chemokines [40], and nitric oxide (NO) [41]. The Ca^{2+} signal forms are rather heterogeneous and include random profiles and rhythmic oscillations as well as bursting activity.
Spontaneous activity has been reported for astrocytes in culture [42–44] and acutely isolated slices of several brain regions [31, 45–47]. Spontaneous astrocytic Ca^{2+} signaling is modulated by pathological events such as those in reactive astrocytes surrounding stab wounds in the neocortex that lack spontaneous Ca^{2+} activity or astrocytes in the neocortex under epileptic conditions that show more spontaneously active astrocytes [45].
Despite all these insights, the common mechanism of this spontaneous activity is still not fully understood. In particular, the higher abundance of spontaneous activity in cultured astrocytes compared to astrocytes within acute brain slices is only barely investigated [24]. We have shown in a comparison study that these differences can be explained to a large extent by the different temperatures at which the astrocytic activity is typically measured. Here we revisit the experimental data and integrate it with our mechanistic multiscale modeling framework to quantitatively test our hypothesis on the molecular foundation of the temperature dependency.
2. Results
Motivated by the typically observed and reported significant differences in Ca^{2+} singling between cultured astrocytes and astrocytes in acute brain slices, we reconsidered the experimental workflow [24]. The most dramatic difference between the two experimental protocols is the different temperatures these two cell systems are typically measured at [21, 48]. In thermodynamics the general impact of temperature on an individual reaction rate is well approximated by the Arrhenius law by , where denotes the approximate reaction rate, describes the activation energy of the reaction, and is the universal gas constant [49]. While the Arrhenius law describes the temperature dependency of an isolated reaction in a wellstirred reactor, the impact of temperature in more complex reaction systems may differ from the simple exponential relation. Alternatively it might be possible that the thermal energy in the measured temperature range is of the order of the activation energy, where the Arrhenius law shows a linear dependency. To explore if temperature could explain the Ca^{2+} signaling phenomenon, we summarize our careful analysis of temperature induced dynamics changes and use here computational modeling to investigate potential reasons for this phenomenon.
2.1. Experiments Show Strong Temperature Dependency in Cortical and Cultured Astrocytes
In the experimental study [24] we used cortical slices obtained from P 10–P 12 mouse brains and cultured astrocytes. The temperature was either changed abruptly or smoothly either by changing the heating temperature of the water bath for the perfusion buffer or by switching from one perfusion buffer to another, respectively. Thereby the temperature was measured within the recording chamber close to the sample cover slip.
2.1.1. Higher Temperatures Reversibly Suppress [Ca^{2+}]_{i} Spiking Instantaneously
Starting from the standard protocols, we first measured Ca^{2+} activity at room temperature and separately at more physiological temperatures with Fluo4AM loaded to the cells before (more detailed experimental description can be found in [24]). In cortical slices obtained from P 10–P 12 mouse brains measured at room temperature, astrocytes showed robust spontaneous signaling activity defined by long durations and large amplitudes of spikes. In contrast, when slices were measured at physiological temperatures, only very few cells showed spontaneous activity.
To compare this behavior with primary astrocytes cultured from mouse neocortex, we measured Ca^{2+} activity at room temperature (20°C) and then abruptly changed the perfusion buffer to one of 28°C or of 30°C in another setup. For these cultured astrocytes we observed a similar significant decrease of cells exhibiting Ca^{2+} signals already within the first minute after the temperature shift (Figure 1). Out of the active cells under 20°C (72 ± 24% of all cells in 4 independent experiments) more than a quarter of cells have not shown any spikes under 28°C and under 30°C only 20% exhibited at least a single spike. Remarkably, this dramatic decrease was fully reversible when perfusing the cells again with buffer at room temperature. Analogous results about the reversibility were also found in experiments with acute brain slices and the reversibility was independent of the experimental order of temperature changes.
(a)
(b)
(c)
(d)
(e)
(f)
The observed reduction in activity, starting at 28°C, was significantly smaller compared to recordings with the higher temperature of 30°C (Figure 1), indicating that between 28°C and 30°C is a thresholdlike temperature with regard to spike activity. The rather instantaneous response to temperature changes and the reversibility of the change in the activity pattern suggests that temperature is directly modifying the characteristic property of some molecular properties of the Ca^{2+} tool kit and not acting primarily on transcriptional regulation.
2.1.2. Temperature Changes Modulate Signal Form of [Ca^{2+}]_{i} Spike
Careful inspection of the time courses reveals also a strong temperature dependent on the duration and signal form of the Ca^{2+} transient in both cultured astrocytes and astrocytes in acute brain slices as shown in Figure 1. To quantify the effect on single transients, we measured the spike width (SW) of individual Ca^{2+} transients by the full width half maximum (FWHM) method. Despite the cell specific heterogeneity, spike durations of the Ca^{2+} transients of astrocytes in acute brain slices decreased significantly when the temperature was raised from 24°C to 37°C. The decrease from 45 ± 17 s to 13 ± 4 s was fully reversible, since lowering the temperature back to 24°C increased the spike duration again to 45 ± 18 s (10 slices, 42 cells at 70 spikes; 53 spikes; 71 spikes). We additionally determined the spike duration at 31°C which averaged at 22 ± 5 s (9 slices, 29 cells at 43 spikes) (Figure 1).
Within the analogous analysis of cultured astrocytes, we found that the spike duration changed from an average of 63 ± 18 s at 20°C to 21 ± 5 s at 30°C. This effect of temperature on spike duration was again reversible as spike duration increased again to 72 ± 20 s at 20°C after lowering the temperature from 30°C (4 independent experiments; at 15 cells, 26 spikes; at 30 cells, 56 spikes) (Figure 1). Within the other experimental setup, the average spike duration was 27 ± 6 s for 28°C as compared to 75 ± 16 s at 20°C (4 experiments; 30 cells at 103 spikes; 105 spikes).
The similarity of the signal form and particularly of the quantified spike widths of cells obtained from cultures and brain slices around 30°C points to temperature as the main reason for the typically reported differences in the signals. To analyze this behavior further we plotted the temperature dependence of the spike duration by pooling data from the two different experimental systems. The resulting data points can be fitted nicely by the exponential function with s and /°C indicating a functional relation between spike width and temperature.
Within the experimental comparison study we also investigated the role of extracellular Ca^{2+} influx and intracellular Ca^{2+} release from intracellular stores. While Ca^{2+} free buffer has reduced the total number of responsive cells only slightly to 78 ± 34% ( slices, averaged cell number/slice = 31 ± 15), blocking SERCA pumps by loading 1 μM thapsigargin for 15 min has abolished Ca^{2+} signals in more than 90% of cells at 31°C. Overall, these findings indicate a rapid change of the intracellular signaling mechanism. In particular, the astonishing reversibility of the Ca^{2+} activity and its signal forms when resetting the temperature to its previous value after tens of minutes suggests a molecular activity dependency on temperature.
2.2. Modeling the [Ca^{2+}]_{i} Temperature Dependence by Modulated SERCA Activity
As previously speculated, the experimentally reported temperature dependent activity of the SERCA pumps [26, 50] may explain the observed changes in signal forms. These studies report an increase of the pump strength with increasing temperature from 5°C to 38°C. Thus, the Ca^{2+} uptake at 5°C amounts only 5% of the uptake at 38°C. This can have a large effect on CICR that represents the spike generating mechanism because the IP_{3}Rs are not homogeneously distributed within the cell but form channel clusters that are typically separated by 1 μm or more. This separation may lead together with Ca^{2+} binding buffers and increasing SERCA activity to steeper gradients close to open channel clusters and subsequently to a weaker coupling between release sites and hence to a reduced CICR.
The data shown in Figure 1 support this perspective. The individual spike form is influenced by the stronger pumping as demonstrated experimentally by the decreased spike width for larger temperatures. Moreover, the spike widths of the two different cell types could be fitted nicely by one exponential function of the temperature. To test this SERCA based hypothesis we apply here our developed physiological multiscale model that is able to generate the whole spectrum of known signals and was used to investigate the experimentally quantified stochastic nature of spiking including the functional robustness of the  relation [27].
2.2.1. Multiscale Modeling Approach
Based on the observations that Ca^{2+} spiking is stochastic and obeys a nontrivial signaling signature, our multiscale model considers the molecular dynamics of individual IP_{3}Rs and Ca^{2+} diffusion as well as buffer dynamics [27]. Our modeling strategy for this coupled system is based on separation of the two involved length scales. On the microscopic scale we use a detailed model for the channels [25]. Those will be open and close in a stochastic manner due to stochastic binding of the signaling molecules IP_{3} and Ca^{2+}. The stochastic channel dynamics acts as a noisy source term that is driving the threedimensional reaction diffusion system (RDS) representing the cytosol. Thereby the model also takes reactions of Ca^{2+} with buffer molecules like BAPTA or calmodulin and the SERCA pumps into account. From this perspective the cellular dynamics is given by an in general nonlinear reaction diffusion system that considers the cytosolic and lumenal dynamics by where the first two equations describe the free Ca^{2+} concentration in the cytosol and the ER, respectively, and the other two equations correspond to a variety of cytosolic and lumenal buffers and . The first terms on the right hand side characterize diffusion with the individual diffusion coefficients . The following terms in (1a) and (1b) describe the leak flux (with permeability ) and release through channels (). The space dependence of reflects the location of channel clusters and the time dependent stochastic opening and closing of the channels. This time dependent source term is governed by stochastic simulation of the channel states. is the volume ratio . The third term results from Ca^{2+} uptake into the ER from the cytosol by SERCA pumps with a flux strength and the SERCA dissociation constant . The remaining terms in (1a)–(1d) describe the reactions of Ca^{2+} with cytosolic and lumenal buffers and , respectively, where denotes the capture rate and is the dissociation rate.
The cellular Ca^{2+} dynamics based on the CICR mechanism results from the local properties of the IP_{3}R channels and their spatial coupling by Ca^{2+} diffusion that is influenced by buffer reactions and SERCA pump activity. Thus, the cellular Ca^{2+} signaling is determined by the channel states and corresponding release fluxes that shape the signals in dependence on the properties of the cellular diffusion. In turn, the resulting Ca^{2+} concentration profiles influence the channel opening behavior due to its dependence on the local Ca^{2+} concentration.
Standard simulation techniques for reaction diffusion systems like finite difference or finite element methods have to calculate the concentration profiles for the complete numerical mesh for convergence criteria. The large concentration gradients typically observed close to open channels require a fine spatial discretization that leads in particular in three spatial dimensions to large computational costs [51]. To circumvent this restriction and enable simulations of biological relevant time courses, we developed an analytic solution of the reaction diffusion equations in terms of coupled Green’s function that allows for calculating at specific spatial positions only such as channel cluster locations. These local concentrations subsequently determine potential changes in the stochastic channel dynamics leading to modified source terms in the reaction diffusion system.
In order to solve the reaction diffusion system analytically we restrict the system to one mobile and one immobile buffer, introduce dimensionless variables, and linearize the corresponding equations [27]. To obtain an analytic solution of the resulting reaction diffusion system we have to determine the geometry and boundary conditions. As a natural choice we use a spherical geometry and, motivated by the observation that astrocytic Ca^{2+} signals do not change significantly in Ca^{2+} free medium [24], noflux boundary condition that reflects the properties of the plasma membrane for our in silico cell. With these specifications we could develop the analytic solution by coupled Green’s functions [52, 53] describing Ca^{2+} and buffer dynamics where open channels correspond to source terms. For mathematical details of the implemented Green’s cell algorithm (GCA) see [27] where we have shown how microscopic channel fluctuations are transmitted onto the level of the cell by diffusion thereby matching the experimental findings on cell variability and buffer influence.
For the IP_{3}R description within the GCA, we use the DeYoungKeizer model (DKM) that assumes three binding sites for each of the four subunits per channel [54]. Each of the three binding sites can be free or occupied leading to 2^{3} different states and 12 possible transitions in dependence on the Ca^{2+} concentration C and IP_{3} concentration I that can be visualized on a cube as shown in Figure 2(a). A subunit is active in the state only where the first index equals one if IP_{3} is bound, the second index corresponds to Ca^{2+} binding to an activating site, and the last index describes the occupation of the dominant inhibiting Ca^{2+} binding site. The activating binding sites for Ca^{2+} have a higher affinity compared to the dominant inhibiting site. This is a minimal choice to generate the nonlinearity in Ca^{2+}, which is the source of the bell shaped open probability (Figure 2(b)) and the basis of CICR [25].
(a)
(b)
2.2.2. Modeling the Impact of SERCA Pump Strength on Ca^{2+} Dynamics
Here we use our validated GCA implementation to test the hypothesis on the temperature dependent SERCA activity as a main reason for the observed differences in spiking. For the simulations, we use physiological cellular diffusion parameters listed in Table 1. Our in silico cell has 31 randomly distributed IP_{3}R clusters each of them having a random number of channels between 2 and 16, yielding 258 channels in total. The clusters are separated by at least 1.5 μm and only located within a sphere with 80% of the cell radius to guarantee convergence of Green’s function representation. To test the SERCA hypothesis, this prototype of a cell was simulated with different pump strengths .

Figure 3 shows the cell wide dynamics, where panel (a) exhibits the number of open channels and panel (b) shows the resulting cytosolic Ca^{2+} concentration for pump strengths varying between 22 s^{−1} and 200 s^{−1} what is around the estimated physiological reference value of 86 s^{−1} [55]. A first conspicuous property shown in panel (a) is that the amplitudes of open channels during a spike are rather similar and vary between 25 and 30, rather independently of the pump strengths . Nevertheless, has a huge influence on the time course.
For low pump strengths, Ca^{2+} is removed too slowly from the cytosol for a coordinated signal. During the initial activation from base level the amplitude of open channels becomes maximal and is limited by inhibition according to the channel properties explained in Figure 2. When a channel recovers from inhibition due to the stochastic unbinding of Ca^{2+} from a dominantly inhibiting site of a subunit, the cytosolic Ca^{2+} concentration is still in the Ca^{2+} concentration range where the higher affinity for the Ca^{2+} activating binding site leads to a large open probability shown in Figure 2(b). Subsequently, the channel will open again, increase the local Ca^{2+} concentration, and reinduce channel inhibition. Repeating this scenario, the channels will mainly switch between activation and inhibition back and forth where a rather constant proportion of up to 50% of subunits are inhibited until randomly no reactivation occurs and the cell is reaching its resting state again. This permanent switching between activation and inhibition leads to the observed plateaus for the two lowest pump strengths.
This behavior is also reflected by the resulting cytosolic Ca^{2+} concentrations. For s^{−1}, the cell exhibits a plateau response with some superimposed spiking, which are often observed in experiments with astrocytes from acute brain slices [24, 56–58]. For s^{−1}, the shorter plateaus of open channels are translated into a more oscillationlike behavior of the cytosolic Ca^{2+} concentration exhibiting rather diverse spike amplitudes. The concentration is in a similar range as the cell with the lower pump strength, indicating a regime where local channel inhibition and spatial coupling determined by pumps have an equal influence on the cellular dynamics.
For increasing , the cells exhibit an even more spiking like behavior since long lasting channel activity is absent, and therefore shows more pronounced peaks. This can be again explained by the interplay of inhibition and Ca^{2+} removal. After a spike has occurred and most channels are inhibited, the SERCAs decrease the cytosolic Ca^{2+} concentration sufficiently fast leading to base level when channels recover from inhibition. Subsequently, the cell relaxes to its resting level, from which it can be activated again by CICR, that is, by triggering opening of a channel cluster. For very large pump strengths, the spikes of the number of open channels have still a similar height to those for smaller , but their widths decrease drastically. This induces low and slim peaks of the cytosolic Ca^{2+} concentration.
So far, the simulations seem to support the hypothesis of the decreased SERCA activity leading to an increased activity at lower temperatures. In terms of the experimental findings, a temperature increase in experiments and observed changes in the Ca^{2+} signals correspond to an increase of the pump rate from the region from 42 s^{−1} to 62 s^{−1} to 152 s^{−1} or higher.
For a further detailed analysis, we determine the spike widths of the simulated cytosolic Ca^{2+} peaks as done for the experimental data in Figure 1. Figure 4(a) exhibits the dependence of the spike width on the pump strength . The analysis reveals an exponential dependence on the pump rate as shown by the line defined by . At second glance we observe that the spike width becomes most regular for s^{−1} as depicted by the small error bar that exhibits the smallest relative variation even compared to the smaller error bars for large . This might indicate an optimal dynamical regime.
(a)
(b)
An analogue behavior can be found for the average amplitudes in Figure 4(b). The dependence on can be fitted by the exponential relation , and again the mean amplitude for s^{−1} exhibits the smallest relative variations. The amplitude behavior obtained from simulations has a similar trend like those seen in the experiments but contradicts former theoretical results where the amplitude stayed constant [59].
To analyze the potential optimal regime, we determine the oscillation characteristics and . Indeed, exhibits a minimum of 42 s for a pump rate of 62 s^{−1} as shown in Figure 5(a). Panel (b) displays the dependence of on showing a slope close to one as in earlier experiments and simulations and indicates a deterministic time of 20 s. The spread of the data points is exclusively caused by the different pump strengths. Similar simulations with smaller buffer concentrations have shown an analogue behavior with a shift of the minimum of to larger . This is in line with our consideration, since a lower buffer concentration leads to an increased spatial coupling that can be compensated by higher pump rates.
(a)
(b)
(c)
For an additional test of our conception based on the spatial coupling of release sites, we analyze the local Ca^{2+} concentrations at channel clusters in Figure 5(c). For all pump strengths the concentration peaks are in the range of 100 μM independently of . Hence, the different cell wide cytosolic Ca^{2+} signals shown in Figure 3(b) are not caused by different local properties but are caused by the spatial coupling of release sites. Thus, the results of the temperature experiments are consistent with the hypothesized spatial control mechanism.
2.2.3. Rate Equation Based Modeling
To further investigate the role of the spatial aspect in Ca^{2+} signaling, we model the effect of SERCA activity on the Ca^{2+} profile by a traditional rate equation model implemented by ordinary differential equations (ODEs). These types of models do not consider any spatial aspects but treat the cell as a wellstirred homogeneous entity and were extensively used in the establishment of the Ca^{2+} signaling mechanism [60, 61]. A variety of models have been proposed and applied to diverse aspects of Ca^{2+} signaling where the common mechanism is Ca^{2+} cycling between the ER and the cytosol. Here we use a model based on previous studies which reflects the main Ca^{2+} fluxes according to Figure 2 as well as Ca^{2+} buffering by mitochondria. The governing rate equations for the cytosolic and ER Ca^{2+} concentration readwhere and are the fractions of free Ca^{2+} within the cytosol and the ER [62], is again the volume ratio between cytosol and ER, denotes the mitochondria to cytosol volume ratio, and the corresponding fluxes reflect typical physiological characteristics. The 2 most important fluxes are the Ca^{2+} release from the ER by IP_{3}R () and the SERCA mediated Ca^{2+} removal from the cytosol into the ER () given by respectively. The dominant flux for the spiking is the nonlinear release flux inducing the CICR by the bell shaped dependence on the Ca^{2+} concentration [63] here modeled in analogy to previous approaches [64]. The form of the complementary SERCA mediated Ca^{2+} flux into the ER () reflects the fact that two Ca^{2+} ions are pumped into the ER by the use of one ATP molecule where the ATP concentration is assumed to be constant. The remaining fluxes are given in Table 2 where and describe Ca^{2+} flux from extracellular space and mitochondria, respectively. corresponds to the Ca^{2+} leakage from ER and is the flux into the extracellular space through plasma membrane Ca^{2+} ATPases. The corresponding parameters are also listed in Table 2.

To investigate the influence of SERCA pump strength on Ca^{2+} dynamics, we simulate the ODEs for a wide range of pump strength and three different values of the Ca^{2+} release flux () from IP_{3}R channels and analyze the spike characteristics (Figure 6). This parameter scan shows that the Ca^{2+} period within the rate equation model exhibits a similar behavior as the spatial GCA simulations in terms of a minimal period in dependence on the pump strength (Figure 6(a)). While this finding is not in clear favor for the perspective of spatial coupling of release sites as the spike generating mechanism, the results obtained with the rate equation model given by (2a) and (2b) about the spike width and amplitude exhibit a different behavior than the experimental data.
(a)
(b)
(c)
In the ODE model a higher SERCA pump rate causes more Ca^{2+} within the ER and also a higher Ca^{2+} release into the cytosol because IP_{3}R flux is controlled to a large proportion by the concentration difference and thus by the ER concentration. A high SERCA pump rate leads to a larger spike width as observed in Figure 6(b). These higher fluxes in and out of the ER also yield higher spike amplitudes (Figure 6(c)). Both observations are in contrast to the experimental observations and the GCA simulations, potentially indicating the fact that a spatial distribution of channels is a necessity for the experimentally observed SERCA pump effects.
3. Discussion
Previous experimental findings have shown a clear temperature dependency of Ca^{2+} signals in astrocytes and indicated the different temperatures at which cultured astrocytes and astrocytes in acute brain slices are typically measured as a potential explanation for the reported differences in Ca^{2+} signaling [24]. Based on these observation and experimental studies reporting a temperature dependent SERCA activity, we used here multiscale modeling to investigate if this activity change may explain the experimental finding.
For this purpose we simulated in silico cells with our GCA implementation that allows for controlling all physiological parameters and obtaining mechanistic insights. From the activity patterns of identical in silico cells that only differ in the SERCA activity shown in Figure 3, we can directly conclude that the interplay of the IP_{3}R properties with respect to CICR and inhibition together with the SERCA pumps can generate a wide spectrum of cellular Ca^{2+} signals. In particular, spiking occurs for an intermediate pump rate where Ca^{2+} is removed fast enough from the cytosol to prevent too frequent repetitive stimulation but where the spatial coupling between the release sites is preserved to allow for cell wide coordinated signals.
As a more quantitative analysis we analyzed the spike width in the simulation data in analogy to the experimental findings where we found that the experimental spike width followed an exponential dependency on temperature independently if the width was determined from cultured astrocytes or from astrocytes in acute brain slices (Figure 1). Remarkably, we obtained again an exponential dependency of the spike width from the simulations as shown in Figure 4. The modeling based exponential relation does not depend explicitly on the temperature but depends on the pump strength that was varied in simulations. From the obtained two exponential dependencies of the spike width on the temperature in experiments and on the SERCA activity in simulations we can infer the relation between these two quantities by setting as
This relation predicts a linear relation between the SERCA activity and the temperature as shown in Figure 7 where the model inference is shown by the line. Note that the proportional assumption can lead to an additional scaling factor in the second term of (4) but not in the first one that describes the slope of the relation. A linear relation was also measured experimentally in tuna by others [26]. Within this study no absolute values of the pump strength were measured but relative changes for different temperatures compared to a reference at 30°C were determined as shown in Figure 7 by the data points. The qualitative agreement between the directly measured relative changes (data points in Figure 7) and the spike width based model inference (line in Figure 7) is another strong evidence that modified SERCA activity may induce the change in astrocytic Ca^{2+} signals at different temperatures.
It is worth emphasizing that we varied in our multiscale simulations the pump strength motivated by the experimentally reported temperature dependent SERCA activity but did not include any specific relation. From the unbiased parameter scan on we first observed the general tendency of less coordinated Ca^{2+} signals for higher temperature in accordance with the experimental data. More strikingly, we found an exponential dependency on of the spike width in the simulation data in analogy to the exponential dependency of the spike width on the temperature in experiments. The exponential relation was not included explicitly within the model but followed naturally from the physiological simulations indicating the importance of the spatial signal organization. With the two obtained exponential relations for the spike width from experiments and from simulations, we then inferred the linear relation between SERCA activity and temperature shown by the line in Figure 7. This inferred linear relation is in great qualitative accordance with the previous directly measured temperature dependency of the SERCA activity shown by the data dots in Figure 7. While the previous study has motivated the modeling strategy of modulating the pump strength , no specific relation has been assumed but is only caused by the spatial signal mechanism.
Strikingly, the two independently obtained relations for the SERCA pump activity both exhibit a linear dependency on the temperature questioning the applicability of the Arrhenius law in this case.
The finding either indicates an activation energy in the order of the thermal energy or a more complex reaction system than what the Arrhenius law is considering. While the activation energy of SERCA is larger than the thermal energy as indicated by the use of ATP for the pump reaction, the underlying effective reaction scheme is more complex. The pump reaction consists of several (potentially interdependent) binding steps of Ca^{2+} and ATP to the SERCA protein. Since each of the binding and unbinding processes will be affected by temperature changes, the resulting reaction may exhibit already different characteristics from the collision of free particles assumed in Arrhenius law. More importantly, we determine here the cell wide averaged SERCA activity. The modulated SERCA activity is changing the spatial coupling between the release sites which leads to a modified hierarchical signal organisation by changing diffusion properties within the diffusion mediated CICR mechanism. This renders the reaction scheme even more complex and may be the main reason for the deviation from Arrhenius law.
When comparing the different scales for the pump strengthtemperature dependence in Figure 7, it becomes obvious that the used SERCA activity parameters in the simulations were smaller than those we would expect from the inference. Furthermore, the slope of the relation from inference is twice as large as the linear fit to the experimental data. This quantitative discrepancy is potentially caused by slightly different spatial characteristic like Ca^{2+} buffer content and physiologic settings of the cells. Nevertheless, the integrative approach shows how molecular properties can be characterized from global cell behavior by mechanistic modeling.
To further evaluate the spatial characteristic of the phenomenon, we were also analyzing local concentrations at channel clusters. Their amplitudes were rather unaffected by the pump strength but the collective IP_{3}R activity was significant modulated. This is a further evidence that it is the cell level that shapes the signals and not local properties. In line with these findings, we found that the rate equation based model (2a) and (2b) is unable to recapitulate the experimental finding when modulating the pump activity. This model approach neglects the spatial scale and thus the spatial coupling between release sites. The failure to model the experimental observations can be also seen as an indirect support for the spatial characteristic of Ca^{2+} signals. In general, the observed experimental behaviors might be also obtained by specific realizations of rate equation models that might consider slightly different flux definitions or integrate further pathways like ATP consumption by SERCA pumps and facilitated ATP production by Ca^{2+} triggered mitochondrial activity. These models are typically rather focused on specific scenarios and often not able to generate the wide spectrum of observed Ca^{2+} signals. In contrast, we have shown that the multiscale approach of the GCA is able to generate all these experimental observations in dependence on reasonable physiological parameters [27].
Interestingly, we found in the multiscale GCA simulations a regime of optimal signaling in terms of regularity of Ca^{2+} spikes. The nonlinear dependence of on the pump strength might indicate a possible control mechanism. A recent study [65] has reported how frequency modulation may control gene regulation. If we therefore assume that the ability to spike and to use frequency coding is the purpose of the Ca^{2+} signaling pathway, cells can control this behavior by the expression level of SERCA. From our findings here, we would expect a negative feedback of SERCA expression on fast oscillations and a positive feedback on slow oscillations with high concentration peaks and plateaus.
The impact of temperature on Ca^{2+} signaling is probably not an astrocytespecific phenomenon. In cardiac muscle the frequency of Ca^{2+} spikes decreases dramatically at 37°C compared with 22°C. In other cell types like in rabbit renal tubules, hepatocytes, and parenchymal and endothelial cells acute hypothermia affects the intracellular Ca^{2+} homeostasis and goes along with a rise in cytosolic Ca^{2+} levels [66–68]. Due to the wide distribution of SERCA and the universality of Ca^{2+} signaling, this mechanism might be used in many temperature sensitive cellular processes, which could make SERCA a key temperature sensor.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This work was supported by the International Research Program of the University Luxembourg by the I1REXBPUL13MANI grant and by funding from the National Biomedical Computation Resource (NBCR) through the NIH P41 GM103426 grant from the National Institutes of Health.
References
 C. Cherniak, “The bounded brain: toward quantitative neuroanatomy,” Journal of Cognitive Neuroscience, vol. 2, no. 1, pp. 58–68, 1990. View at: Publisher Site  Google Scholar
 N. J. Allen and B. A. Barres, “Neuroscience: Glia—more than just brain glue,” Nature, vol. 457, no. 7230, pp. 675–677, 2009. View at: Publisher Site  Google Scholar
 F. A. C. Azevedo, L. R. B. Carvalho, L. T. Grinberg et al., “Equal numbers of neuronal and nonneuronal cells make the human brain an isometrically scaledup primate brain,” Journal of Comparative Neurology, vol. 513, no. 5, pp. 532–541, 2009. View at: Publisher Site  Google Scholar
 N. Baumann and D. PhamDinh, “Biology of oligodendrocyte and myelin in the mammalian central nervous system,” Physiological Reviews, vol. 81, no. 2, pp. 871–927, 2001. View at: Google Scholar
 U. Fünfschilling, L. M. Supplie, D. Mahad et al., “Glycolytic oligodendrocytes maintain myelin and longterm axonal integrity,” Nature, vol. 485, no. 7399, pp. 517–521, 2012. View at: Publisher Site  Google Scholar
 Y. Lee, B. M. Morrison, Y. Li et al., “Oligodendroglia metabolically support axons and contribute to neurodegeneration,” Nature, vol. 487, no. 7408, pp. 443–448, 2012. View at: Publisher Site  Google Scholar
 G. W. Kreutzberg, “Microglia: a sensor for pathological events in the CNS,” Trends in Neurosciences, vol. 19, no. 8, pp. 312–318, 1996. View at: Publisher Site  Google Scholar
 T. Takano, G.F. Tian, W. Peng et al., “Astrocytemediated control of cerebral blood flow,” Nature Neuroscience, vol. 9, no. 2, pp. 260–267, 2006. View at: Publisher Site  Google Scholar
 L. Pellerin, G. Pellegri, P. G. Bittar et al., “Evidence supporting the existence of an activitydependent astrocyteneuron lactate shuttle,” Developmental Neuroscience, vol. 20, no. 45, pp. 291–299, 1998. View at: Publisher Site  Google Scholar
 L. Hertz, “The astrocyteneuron lactate shuttle: a challenge of a challenge,” Journal of Cerebral Blood Flow and Metabolism, vol. 24, no. 11, pp. 1241–1248, 2004. View at: Google Scholar
 M. Nedergaard, “Garbage truck of the brain,” Science, vol. 340, no. 6140, pp. 1529–1530, 2013. View at: Publisher Site  Google Scholar
 N. J. Allen, “Synaptic plasticity: astrocytes wrap it up,” Current Biology, vol. 24, no. 15, pp. R697–R699, 2014. View at: Publisher Site  Google Scholar
 H. R. Parri, T. M. Gould, and V. Crunelli, “Spontaneous astrocytic Ca^{2+} oscillations in situ drive NMDARmediated neuronal excitation,” Nature Neuroscience, vol. 4, no. 8, pp. 803–812, 2001. View at: Publisher Site  Google Scholar
 M. M. Halassa, T. Fellin, and P. G. Haydon, “The tripartite synapse: roles for gliotransmission in health and disease,” Trends in Molecular Medicine, vol. 13, no. 2, pp. 54–63, 2007. View at: Publisher Site  Google Scholar
 M. J. Berridge, “Inositol trisphosphate and diacylglycerol: two interacting second messengers,” Annual Review of Biochemistry, vol. 56, pp. 159–193, 1987. View at: Publisher Site  Google Scholar
 M. J. Berridge, “Inositol trisphosphate and calcium signalling,” Nature, vol. 361, no. 6410, pp. 315–325, 1993. View at: Publisher Site  Google Scholar
 M. J. Berridge, P. Lipp, and M. D. Bootman, “The versatility and universality of calcium signalling,” Nature Reviews Molecular Cell Biology, vol. 1, no. 1, pp. 11–21, 2000. View at: Google Scholar
 M. D. Bootman and M. J. Berridge, “The elemental principles of calcium signaling,” Cell, vol. 83, no. 5, pp. 675–678, 1995. View at: Publisher Site  Google Scholar
 G. Dupont, A. AbouLovergne, and L. Combettes, “Stochastic aspects of oscillatory Ca^{2+} dynamics in hepatocytes,” Biophysical Journal, vol. 95, no. 5, pp. 2193–2202, 2008. View at: Publisher Site  Google Scholar
 A. Skupin and M. Falcke, “Statistical properties and information content of Ca2^{+} oscillations,” Genome Informatics, vol. 19, pp. 69–79, 2008. View at: Google Scholar
 A. Skupin, H. Kettenmann, U. Winkler et al., “How does intracellular Ca^{2+} oscillate: by chance or by the clock?” Biophysical Journal, vol. 94, no. 6, pp. 2404–2411, 2008. View at: Publisher Site  Google Scholar
 A. Skupin and M. Falcke, “Statistical analysis of calcium oscillations,” European Physical Journal: Special Topics, vol. 187, no. 1, pp. 231–240, 2010. View at: Publisher Site  Google Scholar
 A. Skupin and M. Falcke, “From puffs to global Ca^{2+} signals: how molecular properties shape global signals,” Chaos, vol. 19, no. 3, Article ID 037111, 2009. View at: Publisher Site  Google Scholar
 C. G. Schipke, A. Heidemann, A. Skupin, O. Peters, M. Falcke, and H. Kettenmann, “Temperature and nitric oxide control spontaneous calcium transients in astrocytes,” Cell Calcium, vol. 43, no. 3, pp. 285–295, 2008. View at: Publisher Site  Google Scholar
 N. Christian, A. Skupin, S. Morante, K. Jansen, G. Rossi, and O. Ebenhöh, “Mesoscopic behavior from microscopic Markov dynamics and its application to calcium release channels,” Journal of Theoretical Biology, vol. 343, pp. 102–112, 2014. View at: Publisher Site  Google Scholar
 A. M. LandeiraFernandez, J. M. Morrissette, J. M. Blank, and B. A. Block, “Temperature dependence of the Ca2+atpase (serca2) in the ventricles of tuna and mackerel,” The American Journal of Physiology—Regulatory Integrative and Comparative Physiology, vol. 286, no. 2, pp. R398–R404, 2004. View at: Google Scholar
 A. Skupin, H. Kettenmann, and M. Falcke, “Calcium signals driven by single channel noise,” PLoS Computational Biology, vol. 6, no. 8, Article ID e1000870, 13 pages, 2010. View at: Publisher Site  Google Scholar  MathSciNet
 K. Thurley, S. C. Tovey, G. Moenke et al., “Reliable encoding of stimulus intensities within random sequences of intracellular Ca^{2+} spikes,” Science Signaling, vol. 7, no. 331, pp. ra59–ra59, 2014. View at: Publisher Site  Google Scholar
 F. Aguado, J. F. EspinosaParrilla, M. A. Carmona, and E. Soriano, “Neuronal activity regulates correlated network properties of spontaneous calcium transients in astrocytes in situ,” The Journal of Neuroscience, vol. 22, no. 21, pp. 9430–9444, 2002. View at: Google Scholar
 L. Pasti, A. Volterra, T. Pozzan, and G. Carmignoto, “Intracellular calcium oscillations in astrocytes: a highly plastic, bidirectional form of communication between neurons and astrocytes in situ,” The Journal of Neuroscience, vol. 17, no. 20, pp. 7817–7830, 1997. View at: Google Scholar
 M. Zonta, M. C. Angulo, S. Gobbo et al., “Neurontoastrocyte signaling is central to the dynamic control of brain microcirculation,” Nature Neuroscience, vol. 6, no. 1, pp. 43–50, 2003. View at: Publisher Site  Google Scholar
 P. Bezzi, G. Carmignoto, L. Pasti et al., “Prostaglandins stimulate calciumdependent glutamate release in astrocytes,” Nature, vol. 391, no. 6664, pp. 281–285, 1998. View at: Publisher Site  Google Scholar
 A. Araque, N. Li, R. T. Doyle, and P. G. Haydon, “SNARE proteindependent glutamate release from astrocytes,” The Journal of Neuroscience, vol. 20, no. 2, pp. 666–673, 2000. View at: Google Scholar
 L. Pasti, M. Zonta, S. Vicini, and G. Carmignoto, “Cytosolic calcium oscillations in astrocytes may regulate exocytotic release of glutamate,” Journal of Neuroscience, vol. 21, no. 2, pp. 477–484, 2001. View at: Google Scholar
 A. Panatier, D. T. Theodosis, J.P. Mothet et al., “Gliaderived dserine controls nmda receptor activity and synaptic memory,” Cell, vol. 125, no. 4, pp. 775–784, 2006. View at: Publisher Site  Google Scholar
 G. Queiroz, P. J. GebickeHaerter, A. Schobert, K. Starke, and I. Von Kügelgen, “Release of ATP from cultured rat astrocytes elicited by glutamate receptor activation,” Neuroscience, vol. 78, no. 4, pp. 1203–1208, 1997. View at: Publisher Site  Google Scholar
 M. L. Cotrina, J. H.C. Lin, A. AlvesRodrigues et al., “Connexins regulate calcium signaling by controlling ATP release,” Proceedings of the National Academy of Sciences of the United States of America, vol. 95, no. 26, pp. 15735–15740, 1998. View at: Publisher Site  Google Scholar
 P. B. Guthrie, J. Knappenberger, M. Segal, M. V. L. Bennett, A. C. Charles, and S. B. Kater, “ATP released from astrocytes mediates glial calcium waves,” The Journal of Neuroscience, vol. 19, no. 2, pp. 520–528, 1999. View at: Google Scholar
 C. M. Anderson, J. P. Bergher, and R. A. Swanson, “Atpinduced atp release from astrocytes,” Journal of Neurochemistry, vol. 88, no. 1, pp. 246–256, 2004. View at: Google Scholar
 M. Eddelston and L. Mucke, “Molecular profile of reactive astrocytes—implications for their role in neurologic disease,” Neuroscience, vol. 54, no. 1, pp. 15–36, 1993. View at: Publisher Site  Google Scholar
 N. Li, J.Y. Sul, and P. G. Haydon, “A calciuminduced calcium inux factor, nitric oxide, modulates the refilling of calcium stores in astrocytes,” Journal of Neuroscience, vol. 23, no. 32, pp. 10302–10310, 2003. View at: Google Scholar
 A. N. van den Pol, S. M. Finkbeiner, and A. H. CornellBell, “Calcium excitability and oscillations in suprachiasmatic nucleus neurons and glia in vitro,” The Journal of Neuroscience, vol. 12, no. 7, pp. 2648–2664, 1992. View at: Google Scholar
 A. C. Charles, “Glianeuron intercellular calcium signaling,” Developmental Neuroscience, vol. 16, no. 34, pp. 196–206, 1994. View at: Publisher Site  Google Scholar
 T. J. Manning and H. Sontheimer, “Spontaneous intracellular calcium oscillations in cortical astrocytes from a patient with intractable childhood epilepsy (rasmussen's encephalitis),” Glia, vol. 21, no. 3, pp. 332–337, 1997. View at: Google Scholar
 A. Tashiro, J. Goldberg, and R. Yuste, “Calcium oscillations in neocortical astrocytes under epileptiform conditions,” Journal of Neurobiology, vol. 50, no. 1, pp. 45–55, 2002. View at: Publisher Site  Google Scholar
 H. R. Parri and V. Crunelli, “The role of Ca^{2+} in the generation of spontaneous astrocytic Ca^{2+} oscillations,” Neuroscience, vol. 120, no. 4, pp. 979–992, 2003. View at: Publisher Site  Google Scholar
 R. Zur Nieden and J. W. Deitmer, “The role of metabotropic glutamate receptors for the generation of calcium oscillations in rat hippocampal astrocytes in situ,” Cerebral Cortex, vol. 16, no. 5, pp. 676–687, 2006. View at: Publisher Site  Google Scholar
 G. Gimpl, F. Kirchhoff, R. E. Lang, and H. Kettenmann, “Identification of neuropeptide Y receptors in cultured astrocytes from neonatal rat brain,” Journal of Neuroscience Research, vol. 34, no. 2, pp. 198–205, 1993. View at: Publisher Site  Google Scholar
 R. Heinrich and S. Schuster, The Regulation of Cellular Systems, Springer, Berlin, Germany, 1996. View at: Publisher Site
 L. Dode, K. Van Baelen, F. Wuytack, and W. L. Dean, “Low temperature molecular adaptation of the skeletal muscle sarco(endo)plasmic reticulum Ca^{2+}atpase 1 (serca 1) in the wood frog (Rana sylvatica),” Journal of Biological Chemistry, vol. 276, no. 6, pp. 3911–3919, 2001. View at: Publisher Site  Google Scholar
 S. Rüdiger, J. W. Shuai, W. Huisinga et al., “Hybrid stochastic and deterministic simulations of calcium blips,” Biophysical Journal, vol. 93, no. 6, pp. 1847–1857, 2007. View at: Publisher Site  Google Scholar
 M. Ozisik, Heat Conduction, John Wiley & Sons, New York, NY, USA, 1993.
 J. Beck, K. Cole, A. HajiSheikh, and B. Litkouhi, Heat Conduction Using Green's Function, Series in Computational and Physical Processes in Mechanics and Thermal Sciences, Taylor & Francis, New York, NY, USA, 1992.
 G. W. De Young and J. Keizer, “A singlepool inositol 1,4,5trisphosphatereceptorbased model for agoniststimulated oscillations in Ca^{2+} concentration,” Proceedings of the National Academy of Sciences of the United States of America, vol. 89, no. 20, pp. 9895–9899, 1992. View at: Publisher Site  Google Scholar
 K. Bentele and M. Falcke, “Quasisteady approximation for ion channel currents,” Biophysical Journal, vol. 93, no. 8, pp. 2597–2608, 2007. View at: Publisher Site  Google Scholar
 A. Araque, G. Carmignoto, and P. G. Haydon, “Dynamic signaling between astrocytes and neurons,” Annual Review of Physiology, vol. 63, pp. 795–813, 2001. View at: Publisher Site  Google Scholar
 L. Pasti, T. Pozzan, and G. Carmignoto, “Longlasting changes of calcium oscillations in astrocytes: a new form of glutamatemediated plasticity,” The Journal of Biological Chemistry, vol. 270, no. 25, pp. 15203–15210, 1995. View at: Publisher Site  Google Scholar
 A. D. Laskey, B. J. Roth, P. B. Simpson, and J. T. Russell, “Images of Ca^{2+} flux in astrocytes: evidence for spatially distinct sites of Ca^{2+} release and uptake,” Cell Calcium, vol. 23, no. 6, pp. 423–432, 1998. View at: Publisher Site  Google Scholar
 M. Falcke, Y. Li, J. D. Lechleiter, and P. Camacho, “Modeling the dependence of the period of intracellular Ca^{2+} waves on SERCA expression,” Biophysical Journal, vol. 85, no. 3, pp. 1474–1481, 2003. View at: Publisher Site  Google Scholar
 T. Meyer, D. Holowka, and L. Stryer, “Highly cooperative opening of calcium channels by inositol 1,4,5trisphosphate,” Science, vol. 240, no. 4852, pp. 653–656, 1988. View at: Publisher Site  Google Scholar
 G. Dupont and A. Goldbeter, “Theoretical insights into the origin of signalinduced calcium oscillations,” in Cell to Cell Signalling, A. Goldbeter, Ed., pp. 461–474, Academic Press, 1989. View at: Google Scholar
 R. Bertram, L. S. Satin, M. G. Pedersen, D. S. Luciani, and A. Sherman, “Interaction of glycolysis and mitochondrial respiration in metabolic oscillations of pancreatic islets,” Biophysical Journal, vol. 92, no. 5, pp. 1544–1555, 2007. View at: Publisher Site  Google Scholar
 M. Falcke, “Reading the patterns in living cells—the physics of Ca^{2+} signaling,” Advances in Physics, vol. 53, no. 3, pp. 255–440, 2004. View at: Publisher Site  Google Scholar
 D.O. D. Mak, S. M. J. McBride, and J. K. Foskett, “Spontaneous channel activity of the inositol 1,4,5trisphosphate (InsP3) receptor (InsP3R). Application of allosteric modeling to calcium and InsP3 regulation of the InsP3R singlechannel gating,” Journal of General Physiology, vol. 122, no. 5, pp. 583–603, 2003. View at: Publisher Site  Google Scholar
 L. Cai, C. K. Dalal, and M. B. Elowitz, “Frequencymodulated nuclear localization bursts coordinate gene regulation,” Nature, vol. 455, no. 7212, pp. 485–490, 2008. View at: Publisher Site  Google Scholar
 J. F. McAnulty, M. S. Ametani, J. H. Southard, and F. O. Belzer, “Effect of hypothermia on intracellular Ca^{2+} in rabbit renal tubules suspended in UWgluconate preservation solution,” Cryobiology, vol. 33, no. 1, pp. 196–204, 1996. View at: Publisher Site  Google Scholar
 J.S. Kim and J. H. Southard, “Alteration in cellular calcium and mitochondrial functions in rat liver during cold preservation,” Transplantation, vol. 65, no. 3, pp. 369–375, 1998. View at: Publisher Site  Google Scholar
 P. Haddad, J.C. Cabrillac, D. Piche, L. Musallam, and P.M. Huet, “Changes in intracellular calcium induced by acute hypothermia in parenchymal, endothelial, and Kupffer cells of the rat liver,” Cryobiology, vol. 39, no. 1, pp. 69–79, 1999. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2015 Niko Komin et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.