PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of ploscompComputational BiologyView this ArticleSubmit to PLoSGet E-mail AlertsContact UsPublic Library of Science (PLoS)
 
PLoS Comput Biol. 2010 August; 6(8): e1000870.
Published online 2010 August 5. doi:  10.1371/journal.pcbi.1000870
PMCID: PMC2917103

Calcium Signals Driven by Single Channel Noise

Jorg Stelling, Editor

Abstract

Usually, the occurrence of random cell behavior is appointed to small copy numbers of molecules involved in the stochastic process. Recently, we demonstrated for a variety of cell types that intracellular Ca2+ oscillations are sequences of random spikes despite the involvement of many molecules in spike generation. This randomness arises from the stochastic state transitions of individual Ca2+ release channels and does not average out due to the existence of steep concentration gradients. The system is hierarchical due to the structural levels channel - channel cluster - cell and a corresponding strength of coupling. Concentration gradients introduce microdomains which couple channels of a cluster strongly. But they couple clusters only weakly; too weak to establish deterministic behavior on cell level. Here, we present a multi-scale modelling concept for stochastic hierarchical systems. It simulates active molecules individually as Markov chains and their coupling by deterministic diffusion. Thus, we are able to follow the consequences of random single molecule state changes up to the signal on cell level. To demonstrate the potential of the method, we simulate a variety of experiments. Comparisons of simulated and experimental data of spontaneous oscillations in astrocytes emphasize the role of spatial concentration gradients in Ca2+ signalling. Analysis of extensive simulations indicates that frequency encoding described by the relation between average and standard deviation of interspike intervals is surprisingly robust. This robustness is a property of the random spiking mechanism and not a result of control.

Author Summary

The number of proteins organizing cellular processes is huge. The challenge for systems biology is to connect the properties of all these proteins to cellular behavior. Do individual state changes of molecules matter for cell behavior despite these large numbers? Recently, we have experimentally shown for four cell types that intracellular Ca2+ signalling is driven by single channel dynamics. Molecular fluctuations are used constructively for a stochastic spike generation mechanism. The hierarchical structure of Ca2+ signalling prevents averaging of fluctuations and, consequently, the sequence of global spikes still reflects this molecular noise. Here we present a stochastic 3-D multiscale modelling tool living up to these findings by following the consequences of individual channel state changes up to cell level. We simulate the variety of cell responses in different experiments. The stochastic spike generation mechanism is surprisingly robust, providing new insights into the relation of function and robustness. The modelling concept can be applied to a large class of reaction-diffusion processes including other pathways like cAMP.

Introduction

Cellular behavior is the dynamics emerging out of molecular properties and molecular interactions. Hence, cells are indispensably subject to intrinsic noise due to the randomness of diffusion and molecule state transitions in gene expression [1], [2], signaling pathways and control mechanisms. It drives noise induced cell differentiation [3], cell-to-cell variability of cloned cells [4] or second messenger dynamics [5]. While noise in gene expression can be attributed to small molecule numbers, we consider here noise in signalling pathways which occurs even in systems with large molecule numbers.

Molecular interactions create nonlinear feedback like substrate depletion and allosteric regulation in enzyme kinetics or mutual activation of ion channels in membrane potential dynamics. They also couple active molecules inside cells spatially by diffusion of product and substrate or electric currents. If this coupling is strong enough, cells respond spatially homogeneous. Otherwise, we observe dynamic spatial structures formed by concentrations of molecules in specific states. These structures are often called microdomains [6][9].

The existence of these dynamic structures determines in some systems whether the cell obeys deterministic or stochastic mechanisms. The dynamic compartmentalization of the cell by concentration gradients may prevent the establishment of deterministic dynamics by the law of large numbers even if the total number of molecules in the cell would suggest it otherwise. Microdomains are too small to behave deterministically. Not even the whole ensemble of microdomains will behave deterministically, if they are only weakly coupled or if there are only a few of them. Consequently, noise is not averaged out on cell level.

To determine whether we deal with a deterministic or stochastic system is important since these regimes may exhibit very different dependencies of behavior on system parameters [10]. For instance, repetitive spiking in intracellular An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e001.jpg signalling would be restricted to parameter values providing oscillatory dynamics with a deterministic mechanism [11], [12]. It may occur with a stochastic system also for parameters which would lead to bistable or excitable dynamics in the deterministic limit, i.e. for larger or different parameter ranges [13]. In the non-oscillatory parameter ranges, the mechanism creating almost regular spike sequences can be coherence resonance [14][16] rather than the existence of a limit cycle in phase space of the local dynamics. Noisy systems with gradients usually show also a dependency of system characteristics on parameters of spatial coupling which spatially homogeneous systems do not exhibit. An example is the dependency of the spiking frequency on diffusion properties (see below and [5]).

In summary, the interaction between noise and gradients determines parameter dependencies and mechanisms. Recent experimental and theoretical studies on intracellular An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e002.jpg dynamics taught us that cells may indeed work in this regime and may exhibit repetitive spiking with non-oscillatory local dynamics. Functionally relevant gradients are also observed with intracellular cAMP [8], [17][19], pH [20] and in phosphorylation/dephosphorylation dynamics [21], [22] suggesting that the lessons learned from An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e003.jpg dynamics may also apply to other systems.

One of these lessons is that the randomness of single molecule state changes is carried up from the molecular level to cell level [23], [24]. Cellular An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e004.jpg concentration spikes form random sequences of interspike intervals (ISIs) and that randomness arises from the randomness of single molecule state transitions [5], [25]. Consequently, the fluctuations of cellular signals contain information on single molecule behavior. It is a task for modelling now to establish the relation between these fluctuations and single molecule properties to decode this information.

Systems exhibiting the interaction between noise and gradients require modelling tools which can deal efficiently with the large concentration gradients and with the time scale range from molecular transitions to cell behavior. Here, we present such a modelling concept with the example of intracellular An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e005.jpg dynamics. It simulates all active molecules as stochastic Markov chains with all the individual state transitions and describes diffusion and some bulk reactions deterministically. Active molecules are those carrying the crucial feedbacks and nonlinearities. That allows for linearization of passive bulk reactions and the application of a multi-component Green's function to solve the partial differential equations in the cell analytically. We combine Green's functions with a local quasi-static approximation for the fast concentration changes and diffusion processes at the location of active molecules. That is possible due to the short diffusion time on the molecular length scale of a few nanometers. Since we use Green's functions for the long range concentration profiles we can restrict the calculation of concentration values to the location of active molecules. That renders this method extremely efficient even in 3 spatial dimensions.

We will apply this concept to intracellular An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e006.jpg dynamics and compare simulated time dependent concentrations with single cell time series obtained from cultured astrocytes all measured under the same condition without any stimulation. An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e007.jpg is a ubiquitous second messenger in eukaryotic cells that transmits a variety of extracellular signals to intracellular targets. An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e008.jpg controls fertilization, cell differentiation, gene expression, learning and memory [26]. It triggers secretion in glands, muscle contractions in the heart and transmits apoptosis signals [27], [28].

A main mechanism to increase the cytosolic An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e009.jpg concentration is release from intracellular stores, especially from the sarcoplasmic reticulum by ryanodine receptor channels (RyRs) or the endoplasmic reticulum (ER) by inositol 1,4,5-trisphosphate receptor channels (An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e010.jpg). These channels open in a An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e011.jpg dependent fashion - a self amplifying effect known as An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e012.jpg induced An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e013.jpg release (CICR) [27], [29]. If a single channel opens, An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e014.jpg is released into the cytosol, diffuses to adjacent channels and increases their open probability. Thus release may spread into the entire cell leading to a global cytosolic An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e015.jpg concentration spike.

The inositol 1,4,5-trisphosphate (An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e016.jpg) pathway initiates An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e017.jpg release from the ER in many cell types (including astrocytes [30]), since binding of An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e018.jpg to the An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e019.jpg primes them for activation by An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e020.jpg (Figure 1 in Text S1). The spatial arrangement of An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e021.jpg in channel clusters leads to a hierarchical system with the structural levels channel, channel cluster and cluster array, which is the cell level. An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e022.jpg pumps and buffers generate large gradients close to open channel clusters. Thus, channels within a cluster are strongly coupled and the coupling between clusters is only weak - the geometrical hierarchy entails a hierarchy of coupling strengths.

Stochastic binding of An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e023.jpg and An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e024.jpg to the binding sites of An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e025.jpg leads to random opening of a single channel in a cluster [31], [32]. This causes other channels of the same cluster to open also leading to a puff. An individual cluster is stochastic due to the small number of An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e026.jpg per cluster [33][35]. The opening of a single cluster can only be detected by adjacent clusters due to the strong An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e027.jpg gradients [23], [24], [27], [36], [37]. Since they are again only a few, it remains random whether they are opened by the initial puff. If a supercritical number of puffs arises, release spreads into the whole cell causing a global spike. Thus, due to the hierarchy of coupling strength, randomness is carried up from the channel level to the cell level.

In order to model the hierarchical system, we have to consider the stochastic behavior of individual An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e028.jpg and the spatial heterogeneity of cells induced by An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e029.jpg clustering. That leads to a reaction diffusion system (RDS) with local stochastic source terms. For sufficient fast simulations, we decompose the system into local stochastic dynamics comprising channel state transitions and fast local concentration changes and a deterministic global dynamics for which we derive an analytical solution in form of a three component Green's function (Text S1). The solution is driven by stochastic channel behavior described by a hybrid deterministic-stochastic algorithm. We apply the model to a variety of experiments to demonstrate its potential.

Results

Multi-scale modelling exploiting the hierarchical organization of Ca2+ signals

Our modelling concept simulates active molecules individually by Markov chains, the concentration dynamics in the range of the molecule locally quasi-statically and the diffusional long range coupling by Green's functions. Simulations are orders of magnitude faster than numerical schemes based on spatial grids. Their efficiency derives from the methods which we apply. The use of hybrid deterministic-stochastic algorithms for the Markov chains allows for time steps much larger than traditional Gillespie algorithms. In between stochastic molecule state transitions, we integrate the concentration dynamics. The local quasi-static approximation reduces clusters to spatial An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e030.jpg-function sources which turns integrals into sums. It also substantially reduces the number of modes to be used in the Green's function. And finally Green's function enables us to restrict the calculation of concentration values to the locations of active molecules.

Channel and cluster level

An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e031.jpg dynamics and spatial channel clustering lead to the hierarchical system depicted in Figure 1. An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e032.jpg channels are tetrameres [38]. A single channel opens and closes in dependence on binding and dissociation of An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e033.jpg and An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e034.jpg to the binding sites of its subunits (see below). An open channel conducts a An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e035.jpg current from the ER into the cytosol which is due to the huge concentration difference of up to 4 orders of magnitude across the ER membrane.

Figure 1
IP3R properties and clustering generate a hierarchical system.

An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e045.jpg form clusters on the membrane of the ER consisting of 1 to 10 channels [33], [35]. They physically interact within a cluster and are consequently separated by a few nanometers only [35]. The An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e046.jpg in a cluster are strongly coupled by the large local An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e047.jpg concentration close to open channels.

Typical inter-cluster distances found experimentally are in the range of 1–7 An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e048.jpg [39]. Figure 1A shows a representative example of cluster arrangement used in simulations. Due to cytosolic buffers and SERCAs, the local An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e049.jpg concentrations close to an open channel cluster exhibit large gradients such that coupling between clusters is weak compared to intra-cluster coupling. This leads to the hierarchical organization of An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e050.jpg signals. Stochastic opening of a single channel (blip) is locally amplified by CICR leading to a puff (Figure 1B and D). The concentration gradients keep the probability for activation of adjacent clusters small and only a fraction of puffs activates several neighboring clusters. Once a supercritical number of open clusters is reached, more of them open forming a global signal. In that way, the triggering random opening of a single An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e051.jpg is carried up to the macroscopic scale. The mechanism transforms the fast noise of channel state changes on a millisecond time scale into fluctuations of interspike intervals of tens of seconds as shown in Figure 1D.

An early and widely used channel state model is the DeYoung-Keizer model [40], [41]. It assumes independent subunit dynamics and allocates three binding sites to each subunit as shown in Figure 1C. One site for An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e052.jpg and one for An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e053.jpg that cooperatively activate the subunit. Another binding site with lower affinity for An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e054.jpg inhibits the subunit dominantly. These two different affinities lead to a biphasic dependence of the stationary open probability on the An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e055.jpg concentration (see Figure 1 in Text S1). Only the state An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e056.jpg out of the 8 possible subunit states An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e057.jpg corresponds to an active subunit (Figure 1C), where the first index refers to the An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e058.jpg binding and is 1, if An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e059.jpg is bound and 0 otherwise. Analogously, the second and third index describe An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e060.jpg binding to the activating and inhibiting site, respectively. A channel opens, if at least 3 subunits are in the active state.

The 12 possible transitions between the 8 subunit states correspond to transitions in a state scheme forming a cube (Figure 1C). Some of the transition probabilities depend on the local An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e061.jpg and An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e062.jpg concentrations (Figure 1 in Text S1). In simulations, the transitions are realized by a hybrid deterministic-stochastic algorithm [42], which uses the local An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e063.jpg concentrations and the dissociation rates and binding rate constants given in Table 1 in Text S1.

Since An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e064.jpg within one cluster are close to each other, a cluster can be approximated by one spatial An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e065.jpg-source for the purpose of simulating the cluster current in the long range cellular dynamics. The current depends on the number of open channels An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e066.jpg, the time course of which comes out of the stochastic simulation of channel states. It is proportional to the concentration difference An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e067.jpg across the ER membrane at the location of the channel molecule. Hence, we actually need to solve the complete reaction-diffusion problem to determine it. But the concentration difference at the cluster is not well defined with a An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e068.jpg-source term. Therefore, we calculate the cluster current using a spatially extended cluster with radius An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e069.jpg as described in detail in Ref. [43]. The solution of that problem converges within fractions of a millisecond to its stationary state in the range of the channel molecule [43]. That part of the solution is all we need to calculate the current of the An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e070.jpgth cluster. Using the stationary concentration profiles we obtain:

equation image
(1)

with An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e072.jpg where An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e073.jpg denotes the channel flux constant. An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e074.jpg and An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e075.jpg are the diffusion coefficients of An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e076.jpg in the ER and the cytosol. The cluster radius An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e077.jpg depends on the number of open channels An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e078.jpg and the single channel radius An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e079.jpg. The advantage of the approximation is that it takes local ER depletion into account but only depends on the the spatially averaged concentrations An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e080.jpg and An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e081.jpg, which form the boundary conditions for the local quasi-static approximation (see [43] for details). If channel distances within a cluster are of the order of magnitude of the diffusion length of free An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e082.jpg, the internal cluster geometry becomes relevant. In that case, several An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e083.jpg-functions can be used for one cluster.

The approximation allows as well for determination of the local An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e084.jpg concentration at an open channel cluster resulting from its own current (1) as

equation image
(2)

the validity of which had been shown for the buffer concentrations used here [43]. Note that the total concentration at a cluster is the sum of the concentration (2) and the concentrations induced by currents of other open channel clusters. After closing, the An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e086.jpg concentration is determined by the cellular concentration dynamics (see below) 10 nm apart from the release site.

Cellular concentration dynamics

The modelling strategy for the cellular An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e087.jpg dynamics is based on the separation of two length scales. On the microscopic scale of channel clusters, we use a detailed and stochastic channel model to determine local An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e088.jpg currents. On the macroscopic scale of the cell, we use a linearized spatial bi-domain model, and Green's function to integrate it. The microscopic scale determines the currents representing the An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e089.jpg sources of the macroscopic scale. We implement ideas proposed in [43] and use the currents An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e090.jpg of Eq. (1) as the amplitudes of the spatial An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e091.jpg-functions representing the cluster source terms in Eqs. (3). A similar approach was taken by Solovey et al. [44]. We circumvent the concentration divergence at An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e092.jpg-function sources by using Eq. (2) for the value of the local concentration at open clusters. Vice versa, the macroscopic scale affects the concentration values entering the transition rates of the microscopic state schemes.

The ER is a tubular network spreading throughout the cell [45]. Diffusion in such a geometry can be described by diffusion in unrestricted space with a decreased diffusion coefficient [46]. Subsequently, we can superimpose the ER and the cytosol leading to a bi-domain model. Due to the quasi-static approximation (Eq. 1), we do not need to determine the spatially resolved concentration in the ER. Lumenal and cytosolic domains are coupled by a homogeneous An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e093.jpg leak flux An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e094.jpg through the ER membrane, An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e095.jpg re-uptake An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e096.jpg of the ER by SERCA pumps and by the stochastic channel currents An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e097.jpg. Within the cytosol we take free An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e098.jpg, one mobile buffer An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e099.jpg and one immobile buffer An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e100.jpg with the total concentrations An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e101.jpg and An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e102.jpg into account leading to the reaction diffusion equations

equation image
(3a)

equation image
(3b)

equation image
(3c)

where we used buffer conservation and linear pump and leak fluxes with the flux constants An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e106.jpg and An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e107.jpg. An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e108.jpg is the stochastic channel cluster current of the An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e109.jpgth cluster with strength defined by Equation (1).

Scaling concentrations, space and time with typical values reveals the number of independent parameters. It entails the definitions of Table 2. We linearize Eqs. (3), since we would like to use Green's function to solve them. Our parameter values are in the range of the applicability of the linearization to the buffer dynamics as described by Smith et al. [47] for the stationary profiles. We additionally have linearized the pump dynamics. The linearization does not exhibit saturation, which is relevant for calcium concentrations above An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e110.jpg, with An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e111.jpg being the dissociation constant of the pump (Figure 2 in Text S1). These concentrations occur close to open clusters. In that area, the dynamics are dominated by the diffusion term and the channel term, which reduces the relative error due to the linearization of pump and buffer rates substantially. However, if precise knowledge of concentration values close to open channels or clusters is required, the complete non-linear reaction diffusion equations must be solved like e.g. in [42]. The scaled linear reaction diffusion system (Text S1) describes the spatially resolved concentration dynamics by:

equation image
(4a)

equation image
(4b)

equation image
(4c)

where the leak flux depends on the average lumenal concentration, only. All the reaction rate constants depend on the resting state concentration An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e115.jpg, An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e116.jpg and An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e117.jpg due to the linearization: An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e118.jpg, An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e119.jpg and An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e120.jpg. For simplicity we subsumed also An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e121.jpg and An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e122.jpg under An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e123.jpg.

Table 2
Definition of scaling factors and non-dimensional parameters.

The cytosolic concentrations An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e124.jpg are determined by the 3-component Green's function with An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e125.jpg clusters localized at An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e126.jpg (see also Figure 3 in Text S1)

equation image
(5)

with the Bessel function of the first kind An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e128.jpg and the Legendre polynomial An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e129.jpg, where An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e130.jpg is the angle between the source location An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e131.jpg and the point An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e132.jpg given by

equation image
(6)

The An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e134.jpg are determined by the boundary conditions at the plasma membrane (see Text S1).

The three-component response functions An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e135.jpg and An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e136.jpg include the time integration over the source history, i.e. the time dependent channel flux strength An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e137.jpg, and take the buffer reactions as well as the coupling with the ER into account:

equation image
(7a)

equation image
(7b)

with the dimensionless cell radius An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e140.jpg and the normalization factors An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e141.jpg and An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e142.jpg given in the Text S1. The coupling between the cytosol and the ER by An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e143.jpg and An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e144.jpg as well as the reaction rates of An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e145.jpg with the two buffers determine the time constants An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e146.jpg of the response functions (0), which are implicitly given by the roots of the determinant of the coupling matrix

equation image
(8)

The method allows for spatially resolved concentration dynamics as shown in Figure 2 and in the Video S1 by an iso-concentration surface of 2 An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e148.jpg. An initially opening cluster increases the open probability of adjacent An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e149.jpg clusters and release is spreading through the cell until inhibition stops release.

Figure 2
Spatially resolved Ca2+ dynamics.

For the global An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e155.jpg dynamics, the average concentrations are obtained by spatial integration of the analytical solution (9) as

equation image
(9)

where An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e157.jpg denotes the cell radius. The first component of An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e158.jpg describes the cytosolic average concentration An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e159.jpg. With this, the lumenal average An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e160.jpg concentration An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e161.jpg in dimensionless units is determined by

equation image
(10)

which takes into account the leak, pump and channel fluxes, and An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e163.jpg is the volume ratio An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e164.jpg of the cytosol and the ER. An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e165.jpg denotes the equilibrium average lumenal concentration at An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e166.jpg. The difference between the average cytosolic and lumenal concentration An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e167.jpgAn external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e168.jpg determines the cluster current according to Eq. (1) (see Text S1).

The two main approximations of our method are the local quasi-static approximation and the linearization of the passive bulk processes. These assumptions do not allow for a precise study of the intra-cluster concentration dynamics. That can be done with finite element methods like in ref. [42]. The structure of the Green's function solution enables an elegant parallel algorithm that we call the Green's cell. It is orders of magnitude faster than finite element methods and able to simulate long lasting whole cell dynamics in feasible computing time. In the Green's cell algorithm the actual concentration of each cluster is calculated with the Green's function and local quasi-static approximation in dependence on the source history of all clusters by a single process. The concentrations are sent to the master process which determines the corresponding state transition and reaction time by the hybrid algorithm and also calculates the average concentrations. The transition times are re-distributed to the cluster processes where they are used to update the concentrations. For further details see Figure 4 in Text S1.

Stochasticity in measured and simulated Ca2+ signals

Our recent experimental investigation started from the assumption of a random spike generation by wave nucleation followed by a deterministic refractory time. This prediction yields in a linear dependence of the standard deviation on the average period which was also experimentally confirmed [5]. Previous studies report a possible feedback of An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e169.jpg on PKC activity in glutamate stimulated rat astrocytes [48][50]. This may lead to a positive feedback on the An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e170.jpg level by activation of PLCAn external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e171.jpg. The measured relation between standard deviation and average of interspike intervals for spontaneous spiking has a slope equal to 1 [5], demonstrating that spike generation is poissonian and the spike generation probability is constant on the time scale of ISI. Clearly, there is no feedback on that time scale.

To show that the experimental findings are indeed consistent with our ideas of spike generation, we use our modelling tool to study how molecular noise of single channels can be translated into global signals and whether it is sufficient to cause the observed randomness of spike sequences. Figure 3A shows an example of single cell measurements, where the upper panel exhibits the fluorescent signal An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e172.jpg related to the cytosolic An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e173.jpg concentration and the lower panel the individual ISIs. It demonstrates the stochasticity of spiking, since variations in ISIs are in the range of their average. Simulations of a cell with 47 clusters each containing a random number of An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e174.jpg between 4 and 16 exhibit a behavior very similar to experiments showing that single channel noise can lead to time varying ISIs, since there are not any other sources of noise in the simulations (Figure 3B and C). The simulated An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e175.jpg oscillations exhibit in accordance with experimental observations different flavors ranging from rare and irregular spiking to fast and more periodic spiking. The standard deviation An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e176.jpg depends linearly on the average period An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e177.jpg [5]. Recently we have shown that this linear dependence is not a self-evident relation [51]. In particular, it was found that self-sustained oscillatory systems exhibit a different relation than the one observed in An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e178.jpg spiking experiments. The dependence of An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e179.jpg on An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e180.jpg obtained here in simulations is shown in Figure 3D and exhibits a linear dependence with a slope of 1 which was found in experiments for spontaneous oscillations [5], [52]. The offset of the regression line on the An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e181.jpg -axis of about 20 s is the deterministic recovery time.

Figure 3
Stochasticity of Ca2+ oscillations.

Dependence on IP3 and Ca2+ concentrations

The different An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e195.jpgAn external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e196.jpg data points in Figure 3D result from different combinations of the An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e197.jpg and An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e198.jpg base level concentrations, which are both parameters in the model. In vivo the An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e199.jpg concentration is related to the stimulation level by activation of Phospholipase C and An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e200.jpg production. The An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e201.jpg base level is determined by the leak and the pump flux through the ER membrane. In simulations, we adjust the leak flux according to An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e202.jpg and the pump strength. If both concentrations are rather high in the range of An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e203.jpg no spiking occurs since channels are activated as soon as they are in the excitable state (Figure 5 in Text S1). We observe fast and regular spiking (Figure 3C,E and F and Figure 5 in Text S1) for intermediate concentrations. The ISIs have a An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e204.jpg close to the deterministic refractory time, since a new spike is initiated as soon as the recovery time has elapsed. Regular spiking corresponds to cells with small An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e205.jpg in Figure 3D. A further decrease in one of the concentrations increases An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e206.jpg and An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e207.jpg, in a way depending on the other concentration (Figure 3B,E and F). If both concentrations are small, global spiking vanishes and the signal consists of uncorrelated blips.

Different Ca2+ signals in dependence on physiologic parameters

In the previous analysis of the dependence of oscillations on the concentrations, we have already seen that the modelling tool can generate a large spectrum of An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e208.jpg signals ranging from stochastic spiking to almost periodic oscillations. Here, we show that the model can produce all known An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e209.jpg -induced forms of An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e210.jpg signals in dependence on physiologic parameters. Figure 4 exhibits different experimental signal forms and the corresponding simulation results for a cell with 32 clusters. The variety of signals is achieved by varying cell parameters leading to distinct cell responses as shown by the behavior of open channels (black) and number of inhibited subunits (magenta) as well as by the resulting average An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e211.jpg concentration in the cytosol (red) and in the ER (blue). Fast and rather regular oscillations occur by the interplay of activation and inhibition leading to array enhanced coherence resonance as was hypothesized before [5]. This can be seen in the behavior of the channel state dynamics. The number of inhibited subunits (magenta) increases dramatically during a spike and finally inhibition terminates it (Figure 4A). In the following the number of inhibited subunits relaxes slowly towards its resting level. Only very few channels open directly after a spike and these openings do not initiate a new spike, since the number of inhibited subunits is still to high (higher than approximately 220). That causes the deterministic time An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e212.jpg also observed experimentally [5], [52]. But a spike occurs very soon after the number of inhibited subunits has fallen below a critical range since the open probability is rather high with these parameter values. That keeps the stochastic part of the ISI small and spike sequences regular. Moreover, the amplitude of the spike of open channels seems to be smaller, if the spike is initiated at times where the number of inhibited subunits is still high.

Figure 4
Spontaneous Ca2+ signals in individual astrocytes measured under identical conditions (upper row) and simulations of a cell with 32 clusters with different parameters (red line, middle row) exhibit good agreement in the cytosolic Ca2+ ...

We find longer and more irregular ISIs for decreased An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e223.jpg and An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e224.jpg base level concentrations, since the probability of a channel opening is decreased. As a consequence, the cell relaxes to a resting state between spikes with only a few inhibited subunits (Figure 4B). The spike amplitudes of both the number of open channels and of the average An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e225.jpg concentration are slightly increased compared to the regular spiking.

SERCA pumps also shape An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e226.jpg signals. Recent studies have shown that different phenotypes of cloned cells with regard to An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e227.jpg signalling occur due to small variations in SERCA expression levels and activity of RyR [4]. Here, we find that a decreased SERCA activity leads to a burst like behavior (Figure 4C), since An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e228.jpg is removed slower from the cytosol and thus can activate channels which have recovered early from inhibition or channels which have not been activated before.

For even smaller SERCA activity, cells exhibit long lasting plateau An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e229.jpg signals often with superimposed oscillations (Figure 4D). In these cases, released An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e230.jpg stays within the cytosol and reactivates An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e231.jpg again and again. Cooperativeness induced by inhibition leads to superimposed oscillations on the high An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e232.jpg level. The panels of Fig. 4 provide also an idea of cell variability within one cell type and even within one experiment.

Increased randomness by Ca2+ buffers

A direct consequence of the diffusion mediated signal mechanism is the dependence on the strength of spatial coupling by An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e233.jpg diffusion. That coupling strength can be modulated by exogenous An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e234.jpg buffers, since they reduce the diffusion length of free An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e235.jpg. We took advantage of this property of buffers to demonstrate the spatial character of An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e236.jpg oscillations [5]. Note that we used concentrations of An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e237.jpg buffers much smaller than usually applied in order to suppress any kind of An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e238.jpg signal. We measured spiking for several minutes to obtain reference values for ISIs, loaded additional buffer and continued measuring (see Figure 5A). The individual ISIs (blue crosses) are increased and exhibit a larger variability after buffer loading.

Figure 5
Buffers render spiking more irregular by decreasing spatial coupling.

To understand the experimental observation in more detail, we use simulations to analyze the response to additional buffer. Analogously to the experiment, we simulate a fixed cellular arrangement with different mobile buffer concentrations. Figure 5B shows a representative example, where the red and the blue parts correspond to 25 An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e249.jpg and 250 An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e250.jpg EGTA, respectively. Like in the experiment, larger buffer concentration leads to less and more irregular spiking. In the part with the higher buffer concentration, we observe isolated events which do not lead to global waves since coupling of clusters is too weak. These local events are rare in the reference measurements, since a triggering event initiates a global wave very likely.

From population simulations, where individual cells differ in their spatial arrangement of clusters, initial buffer and An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e251.jpg base level concentrations, we obtain the An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e252.jpgAn external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e253.jpg relation shown in Figure 5C, where cells are shifted by an increase of 10 An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e254.jpg in the EGTA concentration. Similar to experiment [52], cells exhibit individual increases of An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e255.jpg and An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e256.jpg with a slope of the shift close to 1 comparable with the population slopes for the two measuring periods.

Influence of buffer kinetics

BAPTA and EGTA are common An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e257.jpg buffers to suppress An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e258.jpg signals and we have used both in experiments [5]. Cells responded much more sensitive to BAPTA than to EGTA. BAPTA has much larger binding and dissociation rate constants than EGTA (Table 1). A disadvantage of the experiment is that the buffer is loaded into a cell by its esterificated form and the total amount that has entered is unknown and difficult to control. Here, we use modelling to illuminate the influence of the different buffer kinetics and concentrations of EGTA and BAPTA.

Table 1
Physiologic standard parameters used in simulation if not stated otherwise.

Figure 5D shows the dependence of An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e330.jpg for fixed cell parameters on the buffer concentration in magenta for EGTA and in black for BAPTA, where squares denote simulations with a single channel current of 0.12 pA and the dots correspond to 1.2 pA. The larger current was achieved by an increased lumenal An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e331.jpg concentration. Cells only differ in the buffer type. We see that increasing BAPTA has a stronger effect than EGTA, which is mainly caused by the larger capture rate. Moreover, we observe a nonlinear dependence of An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e332.jpg on the buffer concentration. The nonlinearity explains the individual shifts of cells in the An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e333.jpgAn external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e334.jpg plane shown Figure 5D. The comparison of the two different current strengths for BAPTA (black) indicates the role of spatial coupling. Higher currents lead to stronger coupling, and subsequently increasing buffer concentrations have a smaller effect on An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e335.jpg.

Cell characteristics in dependence on single channel currents

From the buffer simulations, we can determine the An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e336.jpgAn external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e337.jpg relation shown in Figure 5E. For the smaller currents, there is no qualitative difference between EGTA and BAPTA. Both exhibit a slope close to 1 as shown by the regression lines and an estimated deterministic time of 20 s. The simulations with higher cluster currents indicate a similar deterministic refractory period but the slope of the An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e338.jpgAn external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e339.jpg relation decreases to approximately 0.6. This might explain the experimentally found differences between cell types. Larger currents lead to stronger coupling on the macroscopic length scale and hence to smaller variations.

To confirm these findings and to test the dependency of the slope on other parameters, we analyze spiking of cells for the two different single channel currents. In each simulation set the cells have identical properties and differ only with respect to the buffer content leading to the distinct An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e340.jpg and An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e341.jpg values in Figure 5E (see also Section 6 in Text S1). From these values we determine the population slopes An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e342.jpg. Figure 5F shows An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e343.jpg averaged over different spatial arrangements, An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e344.jpg concentrations (stimulation levels) and pump strengths (see Figure 6 in Text S1). Analogously, we investigated An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e345.jpg, An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e346.jpg and An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e347.jpg (data not shown). The results are very similar to those with An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e348.jpg. For smaller single channel current we obtain always a slope close to 1 when varying all 4 cell properties and for the larger current a slope to 0.6. Varying the buffer concentration, spatial arrangement of clusters, An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e349.jpg concentration or pump strength (within certain limits) does not change the An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e350.jpgAn external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e351.jpg relation but only the position of the system on it.

Discussion

We have presented here an efficient modelling concept for An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e352.jpg dynamics in 3 spatial dimensions. It simulates cell behavior starting from individual channels in full detail. Using Green's function and multiscale techniques allow for taking concentration gradients into account and thus for capturing the hierarchy of coupling strengths. The method can simulate up to 4000 seconds real time within 24 h on 8 CPUs for a cell with 32 clusters and 10 channels per cluster. In comparison to grid-based numerical methods, its main advantage is a gain of computational speed of several orders of magnitude, which enables us to simulate whole spike sequences. We demonstrate the potential of this modelling concept by simulating a variety of experiments. We compare the in silico data with time series obtained from spontaneous oscillations in cultured astrocytes, but several of the results will also apply to other cell types like those analyzed in [5].

These recent experiments showed for 4 different cell types that the sequences of interspike intervals in An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e353.jpg signalling are random [5]. In line with the ideas on the An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e354.jpg signalling mechanisms, we assumed single molecule state transitions to be a sufficient source of noise. We confirm this assumptions with our simulations here in which these state transitions are the only source of randomness. The fluctuations are carried up through the structural levels due to the existence of concentration gradients and hierarchies of coupling strength.

With our bottom-up modelling approach, we were able to generate all experimentally known An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e355.jpg signal types in dependence on physiologic parameters. Spiking exhibits the random ISI sequences observed experimentally with fast regular sequences and slow irregular ones. In particular, the dependency on parameters of spatial coupling observed in experiments is reproduced. We find a sigmoidal response of the An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e356.jpg concentration upon very strong stimulation or strong spatial coupling, which is well known as over stimulation. We observe also bursting. We do not compare our bursting simulations with specific experiments here, but we would like to mention a general aspect. This signal type is usually ascribed to the existence of a dynamic feedback like store depletion or inhibition of An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e357.jpg production which terminates bursts. Such a feedback is not required with a stochastic model. The random length of bursts in our stochastic model offers also a simple explanation for the irregular burst length observed in experiments.

With our method we are able to follow the An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e358.jpg dynamics from the molecular to the cellular level. The single molecule fluctuations determine the global behavior, since they initiate cellular signals. Simultaneously, the local rough channel signal is smoothed on the cell level by the hierarchical system due to diffusion. The universality and variety of signalling cross talks between An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e359.jpg signalling and other pathways render An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e360.jpg a potential source of noise in cellular systems. The fluctuations can be used for cell variability [4] with regards to gene regulation [53], [54] and cell differentiation [3] and provides a flexibility to changing external conditions which is needed during evolution [55].

The An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e361.jpgAn external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e362.jpg relation and functional robustness

Both the experiments and simulations show a simple linear relation between the standard deviation of ISI An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e363.jpg and the average ISI An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e364.jpg. The existence of this linear relation turned out to be surprisingly robust. It survives even an increase of the single channel current by an order of magnitude. This relation describes for each individual cell the response to stimulation changes. Cells shift the spike pattern from slow and irregular to fast and regular along the An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e365.jpgAn external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e366.jpg relation when we increase stimulation. That supplements the current ideas on frequency encoding [54], [56].

At the same time, the An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e367.jpgAn external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e368.jpg relation describes the outcome of spiking experiments with a group of cells. In the experiments, we subjected a sample of cells to the same protocol, and we obtained as many different responses as there are cells in the sample [5]. That set of responses is not arbitrarily scattered across the An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e369.jpgAn external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e370.jpg -plane but aligns along the An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e371.jpgAn external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e372.jpg relation. All the variability among individual cells with respect to expression levels of pathway components, cell volume, ER volume, shape, ion concentration, etc. does not lead to severe deviations from this An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e373.jpgAn external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e374.jpg relation. An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e375.jpg spiking is robust against variability of many pathway components in the sense that the An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e376.jpgAn external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e377.jpg relation is robust. We learn from the simulations here, that it is rather the stochastic spike generation mechanism than control and regulation which provides for this robustness.

If we call the An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e378.jpgAn external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e379.jpg relation from a single cell obtained by parameter changes individual relation and that obtained from a sample of cells population relation, we can describe our findings as identity of individual and population relation.

We could reproduce the variability within a population of cells in simulations by varying cluster array geometry, pump strength, stimulation or buffering conditions. Changing these parameter values simply shifted the system on the An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e380.jpgAn external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e381.jpg relation and did not modify the relation. But changing the single channel current by one order of magnitude did change the slope of the An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e382.jpgAn external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e383.jpg relation.

That suggests a mathematical definition of robustness which accounts for the fact that cells should be able to execute certain functions (e.g. to spike with a range of ISI), but not necessarily at the same strength of stimulation or normalized values of other parameters. We denote with An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e384.jpg and An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e385.jpg two variables describing the function (e.g. An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e386.jpg and An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e387.jpg), and with An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e388.jpg,…,An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e389.jpg and An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e390.jpg,…,An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e391.jpg two sets of parameters (e.g. stimulation strength, temperature, cell volume). The relation between An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e392.jpg and An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e393.jpg is robust with respect to value changes of parameters An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e394.jpg, if it has the structure An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e395.jpg. The parameters An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e396.jpg change only the value of An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e397.jpg while the An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e398.jpg control also the properties of An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e399.jpg, i.e. the properties of the pathway. We call this robustness of the function An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e400.jpg functional robustness (in difference to the robustness of the value of An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e401.jpg). If we identify the stimulation strength with An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e402.jpg, all cells distinguished by the values of An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e403.jpg only can realize frequency encoding with the same An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e404.jpgAn external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e405.jpg relation by varying An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e406.jpg. They can realize this function also by varying another An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e407.jpg-parameter or several of them: function and functional robustness are closely related.

The statement on robustness can also be interpreted with respect to identity of pathways converging onto An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e408.jpg spiking. An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e409.jpg signals can be caused by many different stimuli. The pathways upstream from An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e410.jpg responding to the stimuli must differ with respect to their value of the An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e411.jpg, in order to be distinguishable by pathways downstream from An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e412.jpg.

In summary, cells realize frequency encoding - the function of An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e413.jpg spiking - by mainly moving up and down the relation between standard deviation and average of ISI and to some degree by modulating the deterministic part of the ISI [52]. The An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e414.jpgAn external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e415.jpg relation exists for a stochastic process only, since An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e416.jpg holds for deterministic systems. The An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e417.jpgAn external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e418.jpg relation turned out to be functionally robust with respect to changes of values of one set of parameters. That set may describe cell variability within one cell type or pathway. Changing substantially another set of parameters modified the An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e419.jpgAn external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e420.jpg relation. That set appears rather to specify the identity of pathways converging on An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e421.jpg spiking.

The role of IP3R clusters for astrocyte Ca2+ signalling

Our model predicts that close proximity of An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e422.jpg clusters is a prerequisite for a spontaneous An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e423.jpg response to spread throughout a cell. Indeed, there are types of astrocytes in which An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e424.jpg responses spread within the cell and those, such as Bergmann glia where this is not observed. Interestingly local, subcellular spontaneous An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e425.jpg responses have been recorded which represent functional microdomains [57]. Complementary to the functional units, morphological units are described which are separated from each other by fine processes [58]. It is well conceivable that these thin processes separate endoplasmic reticulum between microdomains by more than 2 An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e426.jpg and according to our model this separation would prevent the spread of a local An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e427.jpg signal to other parts of the cell. In contrast, in cultured astrocytes, the endoplasmic reticulum is preferentially arranged around the cell center without apparent discontinuity [59] and these cells frequently exhibit spontaneous An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e428.jpg responses. In situ, spontaneous An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e429.jpg responses are reported for hippocampal astrocytes and these astrocytes are less polarized as compared to Bergmann glial cells and we would predict that they are less compartimentalized. Indeed, morphological studies indicate that hippocampal astrocytes have five to ten main processes from which smaller extensions branch off [60]. The synchronized activity obviously can spread within the volume of the main processes and soma of hippocampal astrocytes. Moreover, in contrast to culture, the endoplasmic reticulum in astrocytes in hippocampus tissue is preferentially located close to the plasma membrane [59]. These different morphological arrangements result in distinct patterns of An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e430.jpg responses and as a consequence in different gene expression patterns [53].

Do we need such a modelling tool beyond intracellular Ca2+ dynamics

The rise of cell imaging during the last decades illustrated the spatial structure of cells and protein localization. Obviously, cells are not homogeneous and active molecules coupled by diffusional transport are very common. Concentration gradients are functionally relevant, if they create microdomains inside which a pathway is in a state different from its state at other locations in the cell. They have been shown to exist for ‘the other’ fast diffusing intracellular messenger cAMP and in phosphorylation/dephosphorylation dynamics.

Hence, the need for spatially resolved cell models exists and we can apply the modelling concept, if all essential non-linearities are in the discrete active molecules or the boundary conditions and we can linearize remaining bulk reactions. The excellent validity of the linearization for the buffer reactions of An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e431.jpg dynamics has been shown by Smith et al. [47]. We expect a degradation reaction like the cAMP degradation by PDEs also to be linearizable in good approximation. If local concentrations at active molecules should be outside the range of validity of the linearization, that can be fixed by the choice of the local quasi-static approximation of the diffusion process there in many cases. The non-linearities of cAMP production by membrane-bound adenylyl cyclase can be formulated as boundary condition and Green's function must then be used iteratively with an update of the boundary condition in each time step. These remarks illustrate that there is flexibility in the choice of reactions to be linearized which crucially expands the applicability of the concept.

Methods

Cell preparation

Astrocyte cell cultures were prepared from cortex of newborn NMRI mice [61]. Briefly, brain tissue was freed from blood vessels and meninges, trypsinised and gently triturated with a fire-polished pipette in the presence of 0.05% DNAase (Worthington Biochem. Corp., Freehold, NY, USA). Cells were washed twice and plated directly on poly-L-lysine (PLL; 100 An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e432.jpg ; Sigma, Deisenhofen, Germany) coated glass coverslips (An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e433.jpg) at densities of 3 to An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e434.jpg cells/coverslip, kept in An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e435.jpg-10-cm-dishes using Dulbecco's modified Eagle's medium (DMEM) supplemented with 10% fetal calf serum (FCS), 2 mM L-glutamine, 100 units/ml penicillin, and 100 An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e436.jpg streptomycin. One day later, cultures were washed twice with Hank's balanced salt solution (HBSS).

Cells were maintained for at least 4 days and after reaching a subconfluent state, microglial cells and oligodendrocytes as well as their early precursors were dislodged by manual shaking and removed by washing with HBSS. The purity of the astrocytes was routinely determined by immunofluorescence using an antibody against glial fibrillary acidic protein (GFAP, Sigma), a specific astrocytic marker. The cultures typically exhibited more than 90% cells positive for GFAP.

Cell imaging

Cultured cells plated on glass coverslips were measured between p4 and p6. Cells were incubated with the An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e437.jpg indicator dye Fluo-4-acetoxymethyl-ester (Fluo-4 AM, 5 An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e438.jpg, Molecular Probes, Eugene, USA) for 30 min at room temperature in HEPES buffer (148.9 mM NaCl, 5.4 mM KCl, 1 mM An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e439.jpg, 10 mM An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e440.jpg, 10 mM HEPES, 5 mM D-glucose, pH 7.4) containing 0.01% Pluronic-127 (Molecular Probes). Subsequently cells were washed and kept in HEPES buffer for 15–20 min prior to the measurements with the conventional imaging system at a frequency of 0.33 Hz. Cultures were fixed within the microscope chamber of an upright microscope (Axioskop FS, Zeiss, Oberkochen, Germany) equipped with a 20× water immersion objective (UMPlanFl, numeric aperture: 0.5, Olympus, Hamburg, Germany) by a U-shaped platinum wire and superfused with HEPES buffer at 20An external file that holds a picture, illustration, etc.
Object name is pcbi.1000870.e441.jpg. Substances were applied by changing the perfusate. Cells were illuminated (495 nm) from a monochromator (T.I.L.L. Photonics) and fluorescent images (515–545 nm) collected every 3 s with a 12 bit camera (SensiCam) on an upright microscope. At this state, no intercellular waves were observed. Single cell time series were extracted from these images with ImagingCellsEasily software.

Supporting Information

Text S1

Detailed mathematical model and supporting results.

(0.75 MB PDF)

Video S1

The movie shows the free cytosolic calcium concentration during a spike lasting 15 s by an iso-concentration surface of 2 µM. The initial puff activates adjacent channel clusters by increasing their open probability. The clusters open and close randomly until inhibition terminates the release. Parameter values are in Table 1, the spatial arrangement of clusters is shown in Figure 1.

(9.94 MB AVI)

Acknowledgments

We thank Irene Haupt for excellent cell preparation and Carola Schipke for experimental advice.

Footnotes

The authors have declared that no competing interests exist.

A.S. was supported by the International Research Training Group (IRTG) “Genomics and Systems Biology of Molecular Networks” of the German Research Association (DFG) and by the GoFORSYS project no. 0313924 of the Federal Ministry of Education and Research (BMBF). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

References

1. Elowitz M, Surrette M, Hsing W, Swain P. Stochastic gene expression in a single cell. Science. 2002;297:1183–1186. [PubMed]
2. Suel GM, Kulkarni RP, Dworkin J, Garcia-Ojalvo J, Elowitz MB. Tunability and noise dependence in differentiation dynamics. Science. 2007;315:1716–1719. [PubMed]
3. Chang H, Hemberg M, Barahona B, Ingber D, Huang S. Transcriptome-wide noise controls lineage choice in mammalian progenitor cells. Nature. 2008;453:554–557. [PubMed]
4. Nakumura N, Yamazawa T, Okubo Y, Tanabe M, Takeshima H, et al. Temporal switching and cell-to-cell variability in Ca2+ release activity in mammalian cells. Mol Syst Biol. 2009;5:247. [PMC free article] [PubMed]
5. Skupin A, Kettenmann H, Winkler U, Wartenberg M, Sauer H, et al. How does intracellular Ca2+ oscillate: By chance or by the clock? Biophys J. 2008;94:2404–2411. [PubMed]
6. Berridge M. Calcium microdomains: Organization and function. Cell Calcium. 2006;40:405–412. [PubMed]
7. Bootman M, Lipp P, Berridge M. The organisation and functions of local Ca2+ signals. J Cell Sci. 2002;114:2213–2222. [PubMed]
8. Zaccolo M, Pozzan T. Discrete microdomains with high concentrations of cAMP in stimulated rat neonatale cardiac myocytes. Science. 2002;295:1711–1715. [PubMed]
9. Tovey S, Deodes S, Taylor E, Church J, Tayor C. Selective coupling of type 6 adenylyl cyclase with type 2 IP3 receptors mediates direct sensitization of IP3 receptors by cAMP. J Cell Biol. 2008;183:297–311. [PMC free article] [PubMed]
10. Kummer U, Krajnc B, Pahle J, Green AK, Dixon CJ, et al. Transition from stochastic to deterministic behavior in calcium oscillations. Biophys J. 2005;89:1603–1611. [PubMed]
11. Thul R, Falcke M. Stability of membrane bound reactions. Phys Rev Lett. 2004;93:188103-1–188103-4. [PubMed]
12. Thul R, Thurley K, Falcke M. Toward a predictive model of Ca2+ puffs. Chaos. 2009;19:037108. [PubMed]
13. Falcke M, Bär M, Lechleiter J, Hudson J. Spiral breakup and defect dynamics in a model for intracellular Ca2+ dynamics. Physica D. 1999;129:236–252.
14. Pikovsky AS, Kurths J. Coherence resonance in a noise-driven excitable system. Phys Rev Lett. 1997;78:775–778.
15. Falcke M, Li Y, Lechleiter J, Camacho P. Modeling the dependence of the period of intracellular Ca2+ waves on SERCA expression. Biophys J. 2003;85:1474–1481. [PubMed]
16. Lindner B, Neiman A, Schimansky-Geier L. Effects of noise in exitable systems. Phys Rep. 2004;392:321–424.
17. Karpen J, Rich T. Resolution of cAMP signals in three-dimensional microdomains using novel, real-time sensors. Proc West Pharmacol Soc. 2003;47:1–5. [PubMed]
18. Willoughby D, Cooper D. Organization and Ca2+ regulation of adenylyl cyclase in cAMP microdomains. Physiol Rev. 2007;87:965–1010. [PubMed]
19. Leroy J, Abi-Gerges A, Nikolaev V, Richter W, Lechene P, et al. Spatiotemporal dynamics of beta-Adrenergic cAMP signals and l-type Ca2+ channels regulation in adult rat ventricular myocytes: Role of phophodiesterase. Circ Res. 2008;102:1091–1100. [PubMed]
20. Cardone R, Casavola V, Reshkin S. The role of disturbed pH dynamics and the Na+/H+ exchanger in metastasis. Nature Rev Cancer. 2005;5:786–795. [PubMed]
21. Kholodenko B. Four-dimensional organization of protein singnaling-cascades: the role of diffusion, endocytosis and molecular motors. J Exp Biol. 2003;206:2073–2082. [PubMed]
22. Kholodenko B. Cell-signalling dynamics in time and space. Nat Rev Mol Cell Biol. 2006;7:165–176. [PMC free article] [PubMed]
23. Marchant J, Parker I. Role of elementary Ca2+ puffs in generating repetitive Ca2+ oscillations. EMBO Journal. 2001;20:65–76. [PubMed]
24. Marchant J, Callamaras N, Parker I. Initiation of IP3-mediated Ca2+ waves in Xenopus oocytes. EMBO J. 1999;18:5285–5299. [PubMed]
25. Dupont G, Abou-Lovergne A, Combettes L. Stochastic aspects of oscillatory calcium dynamics in hepatocytes. Biophys J. 2008;95:2193–2202. [PubMed]
26. Berridge M, Bootman M, Lipp P. Calcium - a life and death signal. Nature. 1998;395:645–648. [PubMed]
27. Berridge M. Elementary and global aspects of calcium signalling. J Physiol. 1997;499:291–306. [PubMed]
28. Orrenius S, Zhivotovsky B, Nicotera P. Regulation of cell death: the calcium-apoptosis link. Nat Rev Mol Cell Biol. 2003;4:552–565. [PubMed]
29. Bootman M, Berridge M, Lipp P. Cooking with calcium: The recipes for composing global signals from elementary events. Cell. 1997;91:367–373. [PubMed]
30. Fiacco T, McCarthy K. Astrocyte calcium elevations: Properties, propagation, and effects on brain signaling. Glia. 2006;54:676–690. [PubMed]
31. Falcke M. On the role of stochastic channel behavior in intracellular Ca2+ dynamics. Biophys J. 2003;84:42–56. [PubMed]
32. Shuai J, Pearson JE, Parker I. Modeling Ca2+ feedback on a single inositol 1,4,5-trisphosphate receptor and its modulation by Ca2+ buffers. Biophys J. 2008;95:3738–3752. [PubMed]
33. Smith I, Witgen S, Parker I. Localization of puff sites adjacent to the plasma membrane: Functional and spatial characterization of Ca2+ signaling in SH-SY5Y cells utiliting membrane-permeant IP3. Cell Calcium. 2009;45:65–77. [PMC free article] [PubMed]
34. Smith I, Parker I. Imaging and quantal substructure of single IP3R channel activity during Ca2+ puffs in intact mammalian cells. Proc Natl Acad Sci USA. 2009;106:6406–6409. [PubMed]
35. Taufiq-Ur-Rahman, Skupin A, Falcke M, Taylor C. Clustering of IP3 receptors by IP3 retunes their regulation by IP3 and Ca2+. Nature. 2009;458:655–659. [PMC free article] [PubMed]
36. Yao Y, Parker I. Quantal puffs of intracellular Ca2+ evoked by inositol trisphosphate in Xenopus oocytes. J Physiol (Cambridge) 1995;482:533–553. [PubMed]
37. Falcke M. Reading the patterns in living cells - the Physics of Ca2+ signaling. Advances in Physics. 2004;53:255–440.
38. Nucifora F, Sharp A, Milgram S, Ross C. Inositol 1,4,5-trisphosphate receptors in endocrine cells: localization and association in hetero- and homotetramers. Mol Biol Cell. 1996;7:949–960. [PMC free article] [PubMed]
39. Tateishi Y, Hattori M, Nakayama T, Iwai M, Bannai H, et al. Cluster formation of inositol 1,4,5-trisphosphate receptor requires its transition to open state. J Biol Chem. 2004;280:6816–6822. [PubMed]
40. DeYoung G, Keizer J. A single-pool inositol 1,4,5-trisphosphate-receptor-based model for agonist-stimulated oscillations in Ca2+ concentration. Proc Natl Acad Sci USA. 1992;89:9895–9899. [PubMed]
41. Sneyd J, Falcke M, Dufour JF, Fox C. Mathematical models of the inositol trisphosphate receptor. Prog Biophys Mol Bio. 2004;85:121–140. [PubMed]
42. Rüdiger S, Shuai JW, Huisinga W, Nagaiah C, Warnecke G, et al. Hybrid Stochastic and Deterministic Simulations of Calcium Blips. Biophys J. 2007;93:1847–1857. [PubMed]
43. Bentele K, Falcke M. Quasi-steady approximation for ion channel currents. Biophys J. 2007;93:2597–2608. [PubMed]
44. Solovey G, Fraiman D, Pando B, Dawson SP. Simplified model of cytosolic Ca2+ dynamics in the presence of one or several clusters of Ca2+-release channels. Phys Rev E Stat Nonlin Soft Matter Phys. 2008;78:041915. [PubMed]
45. Roderick H, Berridge M, Bootman M. Understanding Calcium Dynamics - Experiments and Theory. Berlin Heidelberg New York: Springer; 2003. pp. 17–35. volume 623 of Lecture Notes in Physics, chapter 2.
46. Ölveczky B, Verkman A. Monte carlo analysis of obstructed diffusion in three dimensions: Application to molecular diffusion in organelles. Biophys J. 1998;74:2722–2730. [PubMed]
47. Smith G, Dai L, Miura R, Sherman A. Asymptotic analysis of buffered calcium diffusion near a point source. SIAM J Appl Math. 2001;61:1816–1838.
48. Codazzi F, Teruel MN, Meyer T. Control of astrocyte Ca2+ oscillations and waves by oscillating translocation and activation of protein kinase C. Curr Biol. 2001;11:1089–1097. [PubMed]
49. Venance L, Stella N, Glowinsji J, Giaume J. Mechanism involved in initiation and propagation of receptor-induced intercellular calcium signaling in cultured rat astrocytes. J Neurosci. 1997;17:1981–1992. [PubMed]
50. Pasti L, Pozzan T, Carmignoto G. Long-lasting changes of calcium oscillations in astrocytes. a new form of glutamate-mediated plasticity. J Biol Chem. 1995;270:15203–15210. [PubMed]
51. Skupin A, Falcke M. From puffs to global Ca2+ signals: how molecular properties shape global signals. Chaos. 2009;19:037111. [PubMed]
52. Skupin A, Falcke M. Statistical properties and information content of Ca2+ oscillations. Gen Inform. 2008;19:69–79.
53. Capite J, Ng S, Parekh A. Decoding of cytoplasmic Ca2+ oscillations through the spatial signature drives gene expression. Curr Biology. 2009;19:853–858. [PubMed]
54. Cai L, Dalal C, Elowitz M. Frequency-modulated nuclear localization bursts coordinate gene regulation. Nature. 2008;455:485–490. [PMC free article] [PubMed]
55. Reanney D. Genetic noise in evolution? Nature. 1984;307:318–319. [PubMed]
56. Dolmetsch R, Xu K, Lewis R. Calcium oscillations increase the efficiency and specificity of gene expression. Nature. 1998;392:933–936. [PubMed]
57. Grosche J, Matyash V, Möller T, Verkhratsky A, Reichenbach A, et al. Microdomains for neuron-glia interaction: parallel fiber signaling to Bergman glia cells. Nature Neurosci. 1999;2:139–143. [PubMed]
58. Grosche J, Kettenmann H, Reichenbach A. Bergman glia cells form distinct morphological structures to interact with cerebellar neurons. J Neurosci. 2002;68:138–149. [PubMed]
59. Pivneva T, Haas B, Reyes-Haro D, Laube G, Veh R, et al. Store-operated Ca2+ entry in astrocytes: Different spatial arrangement of endoplasmic reticulum explains functional diversity in vitro and in situ. Cell Calcium. 2008;43:591–601. [PubMed]
60. Bushhong E, Martone M, Jones Y, Elissman M. Protoplasmic astrocytes in CA1 stratum radiatum occupy seperate anatomical domains. J Neurosci. 2002;22:183–192. [PubMed]
61. Lyon S, Kettenmann H. Oligodendrocytes and microglia are selectively vulnerable to combined hypoxia and hypoglycemia injury in vitro. J Cereb Blood Flow Metab. 1998;18:521–530. [PubMed]
62. Cornell-Bell A, Finkbeiner S, Cooper M, Smith S. Glutamate induces calcium waves in cultured astrocytes: Long range glial signaling. Science. 1990;247:470–473. [PubMed]
63. Suhara W, Kobayashi M, Sagara H, Hamadad K, Goto T, et al. Visualization of inositol 1,4,5-trisphosphate receptor by atomic force microscopy. Neurosci Lett. 2006;391:102–107. [PubMed]
64. Allbritton N, Meyer T, Sryer L. Range of messenger action of calcium ion and inositol 1,4,5 trisphosphate. Science. 1992;258:1812–1815. [PubMed]
65. Thul R, Falcke M. Release currents of IP3 receptor channel clusters and concentration profiles. Biophys J. 2004;86:2660–2673. [PubMed]
66. Jafri M, Keizer J. On the roles of Ca2+ diffusion, Ca2+ buffers and the endoplasmatic reticulum in IP3-induced Ca2+ waves. Biophys J. 1995;69:2139–2153. [PubMed]
67. Irvine R, Letcher A, Heslop J, Berridge M. The inositol tris/tetrakisphosphate pathway –demons of Ins(1,4,5)P3 3 kinase activity in mammel tissue. Nature. 1986;320:631–634. [PubMed]
68. Meldolesi J, Pozzan T. The endoplasmic reticulum Ca2+ store: a view from the lumen. Trends Biochem Sci. 1998;23:10–14. [PubMed]
69. Pape P, Jong D, Chandler W. Calcium release and its voltage dependence in frog cut muscle fibers equilibrated with 20 mM EGTA. J Gen Physiol. 1995;106:259–336. [PMC free article] [PubMed]
70. Richardson A, Taylor C. Effects of Ca2+ chelators on purified inositol 1,4,5-trisphoasphate (Ins(1,4,5)P3) receptors and Ins(1,4,5)P3-stimulated Ca2+-mobilization. J Biol Chem. 1993;268:11528–11533. [PubMed]

Articles from PLoS Computational Biology are provided here courtesy of Public Library of Science