Next Article in Journal
A Fast Factorisation of Semi-Primes Using Sum of Squares
Next Article in Special Issue
Exact Evaluation of Statistical Moments in Superradiant Emission
Previous Article in Journal
Combining the Magnetic Equivalent Circuit and Maxwell–Fourier Method for Eddy-Current Loss Calculation
Previous Article in Special Issue
Structures and Instabilities in Reaction Fronts Separating Fluids of Different Densities
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Functional Ca2+ Channels between Channel Clusters are Necessary for the Propagation of IP3R-Mediated Ca2+ Waves

by
Estefanía Piegari
and
Silvina Ponce Dawson
*
Departamento de Física FCEN-UBA and IFIBA (CONICET), Ciudad Universitaria, Pabellón I, Buenos Aires 1428, Argentina
*
Author to whom correspondence should be addressed.
Math. Comput. Appl. 2019, 24(2), 61; https://doi.org/10.3390/mca24020061
Submission received: 10 May 2019 / Revised: 8 June 2019 / Accepted: 10 June 2019 / Published: 11 June 2019
(This article belongs to the Special Issue Dynamics Days Latin America and the Caribbean 2018)

Abstract

:
The specificity and universality of intracellular Ca 2 + signals rely on the variety of spatio-temporal patterns that the Ca 2 + concentration can display. Ca 2 + release into the cytosol through inositol 1,4,5-trisphosphate receptors (IP 3 Rs) is key for this variety. The opening probability of IP 3 Rs depends on the cytosolic Ca 2 + concentration. All of the dynamics are then well described by an excitable system in which the signal propagation depends on the ability of the Ca 2 + released through one IP 3 R to induce the opening of other IP 3 Rs. In most cell types, IP 3 Rs are organized in clusters, i.e., the cytosol is a “patchy” excitable system in which the signals can remain localized (i.e., involving the release through one or more IP 3 Rs in a cluster), or become global depending on the efficiency of the Ca 2 + -mediated coupling between clusters. The spatial range over which the signals propagate determines the responses that the cell eventually produces. This points to the importance of understanding the mechanisms that make the propagation possible. Our previous qualitative comparison between experiments and numerical simulations seemed to indicate that Ca 2 + release not only occurs within the close vicinity of the clearly identifiable release sites (IP 3 R clusters) but that there are also functional IP 3 Rs in between them. In this paper, we present a quantitative comparison between experiments and models that corroborate this preliminary conclusion. This result has implications on how the Ca 2 + -mediated coupling between clusters works and how it can eventually be disrupted by the different Ca 2 + trapping mechanisms.

1. Introduction

Calcium (Ca 2 + ) signals are ubiquitous across cell types [1,2]. In many cases, they involve Ca 2 + release from the endoplasmic reticulum (ER) into the cytosol through Inositol 1,4,5-trisphosphate receptors (IP 3 Rs), which are Ca 2 + channels [3]. The opening probability of IP 3 Rs increases with the cytosolic Ca 2 + concentration [4], provided that this concentration is not too high. Thus, the very same Ca 2 + that is released through an open IP 3 R can induce the opening of neighboring IP 3 Rs. The combination of this Calcium Induced Calcium Release (CICR) [5] with the diffusion of Ca 2 + between its channels gives rise to propagating signals that can even embrace the whole cell [6,7,8]. IP 3 Rs become inhibited in the presence of high cytosolic Ca 2 + concentrations. From a physical/mathematical point, the dynamics that underlie these propagating Ca 2 + signals are well described by an excitable system. Now, in most cell types, IP 3 Rs are organized in clusters. Waves can then fail to propagate if the amount of Ca 2 + that reaches one cluster is not high enough to “cross” the excitability threshold [9]. We have recently studied the Ca 2 + -mediated coupling between neighboring IP 3 R-clusters by means of experiments in which we used two single-wavelength Ca 2 + -dyes [10]. Single-wavelength Ca 2 + -dyes are Ca 2 + indicators that increase their fluorescence enormously upon Ca 2 + binding, without changing their emission wavelength [11]. The presence of the dyes affects the elicited signals because they “trap” Ca 2 + and, in this way, can disrupt CICR. How the presence of different Ca 2 + trapping mechanisms (particularly, Ca 2 + -binding buffers) affects Ca 2 + signals had been studied with experiments in which varying quantities of exogenous Ca 2 + buffers had been introduced in the cells [12,13]. The experiments of [10] allowed us to make visible the invisible: by observing the signals simultaneously with a slow (Rhod-2) and a fast (Fluo-4) dye, we could infer directly the different ways in which slow and fast buffers affect the signals. Furthermore, performing experiments for various concentrations of the dyes and comparing qualitatively the spatio-temporal distributions of the Ca 2 + -bound to both dyes with those derived from numerical simulations, we inferred that there should be functional IP 3 Rs between the IP 3 R-clusters to explain the observations. This implied that Ca 2 + release not only occurred within the close vicinity of the centers of the clearly identifiable release sites (the clusters) but also in between.
In this paper, we present a more quantitative analysis of the experimental observations. To this end, we compute, from the experiments, the change with the slow dye concentration, Rhod-2, of the probability, P ( | n c ) that, given that there are n c simultaneously open channels in a primary cluster, the event does not induce the opening of IP 3 Rs in other (secondary) clusters (i.e., the event remains localized). We then use a simple model to compute numerically the probability that one IP 3 R located at a distance, d, from a primary IP 3 R-cluster becomes open (after a certain time) given that there are n c simultaneously open channels at the primary cluster. Changing the parameters of the simulation we conclude that the changes observed experimentally can only be explained if d 0.6 μ m (i.e., it is smaller than the typical inter-cluster distance, d 1.4 μ m ) and that the basal Ca 2 + concentration, [Ca] b , is reduced when the slow dye concentration, [ R ] T , is increased.
We present in what follows the experimental results that we use to estimate the changes in P ( | n c ) as [ R ] T is varied. We introduce in Section 2.2 the probabilistic model that we will then use to analyze the experimental data. In Section 2.3, we show the results of the numerical simulations with which we estimate some of the probabilities that enter the probabilistic model. In Section 2.4, we combine the experimental and numerical results and determine for what parameters of the simulations they are compatible within the framework probabilistic model. A discussion is included at the end.

2. Results

2.1. Experimental Results

The experiments analyzed here were presented and described in detail in [10]. A description of how they were performed is included in Materials and Methods. In brief, we elicited IP 3 -mediated Ca 2 + release events in Xenopus laevis oocytes that were previously injected with caged IP 3 , the Ca 2 + buffer, EGTA, and a fast, Fluo-4, and a slow, Rhod-2, Ca 2 + dye. The signals were elicited by uncaging the caged IP 3 with an UV flash. In the current paper, we analyze the changes observed in the distribution of localized Ca 2 + release events (puffs) elicited in this way as the concentration of the slow dye, [ R ] T = [ R h o d 2], is varied but all other experimental parameters are kept fixed ( [ E G T A ] = 90 μ M ; [ F l u o 4 ] = [ F ] T = 36 μ M , duration of the UV flash to uncage the IP 3 = ( 100 200 ) ms; see Table 1). In particular, we are interested in studying the changes in the “size” of the elicited release events that remain localized (Ca 2 + puffs), where by size we mean the number of IP 3 Rs that are simultaneously open at the release site (the cluster) during the release event. In order to compare the localized event size distributions obtained in experiments performed for different values of [ R ] T , we introduced in [10] a quantity, A l i b F , that is an increasing function of the Ca 2 + current that underlies the observed release event regardless of the value of [ R ] T . We describe in Section 4.3 how we compute A l i b F from the fluorescence emitted by the Ca 2 + -bound Fluo-4 molecules. We show in Figure 1 the cumulative distribution functions (CDFs) of A l i b F derived from the experiments performed for the conditions of Set III (dashed line), Set II (dotted line) and Set I (solid line). These CDFs were computed including only localized events, i.e., Ca 2 + puffs.
The Kolmogorov–Smirnov test rejects the null hypothesis that A l i b F from sets III and II and sets III and I come from the same continuous distribution with a 99% significance level ( p v a l u e = 2.1 × 10 4 and p v a l u e = 4.6 × 10 5 , respectively), but cannot reject that the data points from Set II and Set I come from the same distribution ( p v a l u e = 0.96 ). In any case, there is a tail in the CDF of A l i b F for Set I that is unobservable in that of Set II which is consistent with having more events with relatively larger underlying Ca 2 + currents in the former than in the latter ( A l i b F = 2.2 for Set III A l i b F = 2.7 for Set II and A l i b F = 2.9 for Set I). These comparisons indicate that puffs with relatively larger underlying Ca 2 + currents can be elicited as the concentration of the slow dye, Rhod-2, is increased. Although the transformation from puff amplitude to A l i b F involves certain uncertainties, the changes observed in the fluorescence rise time as [ R ] T is varied [10] support this conclusion.

2.2. Probabilistic Model to Analyze the Differences Observed in the Experimental Event Size Distributions for Different Values of [ R ] T

As analyzed in [14], being able to observe localized Ca 2 + release events (puffs) with larger Ca 2 + currents as a slow buffer concentration (in this case, Rhod-2) is increased can be due to a more efficient uncoupling between IP 3 R clusters due to the presence of the slow buffer. Namely, we have the hypothesis that the differences in the CDFs of localized release event sizes illustrated in Figure 1 occur because, as [ R ] T decreases, Ca 2 + release events with too many simultaneously open IP 3 Rs at the primary site can no longer remain localized, induce the opening of IP 3 Rs in neighboring (secondary) clusters and, thus, are not included to compute the CDF. We hereby introduce a way to analyze the experimental data to quantify what fraction of events that are localized for a given [ R ] T turn into waves as [ R ] T is decreased.
We define P ( n c ) as the probability that there are n c simultaneously open IP 3 Rs in a cluster for a given set of experimental conditions. Here, we will assume that all conditions remain the same except for the total slow dye concentration, [ R ] T . Thus, we will analyze the change of P ( n c ) with [ R ] T . Given that there are n c simultaneously open IP 3 Rs we want to distinguish whether this situation induces the opening of at least one IP 3 R in a neighboring (secondary) cluster (i.e., it initiates a wave) or it does not (i.e., the Ca 2 + release event due to the n c simultaneously open IP 3 Rs remains localized). We then write:
P ( n c ) = P ( n c & | [ R ] T ) + P ( n c & w | [ R ] T ) .
In Equation (1), P ( n c & | [ R ] T ) is the joint probability that n c channels are simultaneously open in a cluster and the event stays localized for a given value of [ R ] T . P ( n c & w | [ R ] T ) is the joint probability that n c channels are simultaneously open in a cluster and the event induces the opening of at least one IP 3 R in another cluster for a given value of [ R ] T . The symbol | means that these are two conditional probabilities for a given value of the slow dye concentration, [ R ] T . All the probabilities we work with here are defined over the set of events, i.e., for n c 1 . The aim of this calculation is to assess how the two joint probabilities of Equation (1) change with [ R ] T . Under the assumption that Rhod-2 is a slow buffer and, as such, does not affect CICR within the cluster [14], we consider that P ( n c ) does not depend on [ R ] T . What may change when varying Rhod-2 is whether the event with n c open channels in a cluster remains localized (stays as a puff) or elicits the opening of channels in a neighboring cluster (becomes a wave). We rewrite the two joint probabilities of interest as:
P ( n c & | [ R ] T ) = P ( n c | , [ R ] T ) P ( | [ R ] T ) = P ( n c | , [ R ] T ) 1 P ( w | [ R ] T ) ,
P ( n c & w | [ R ] T ) = P ( w | n c , [ R ] T ) P ( n c ) .
In these equations, P ( n c | , [ R ] T ) is the probability that a Ca 2 + release event that remains localized for a given [ R ] T corresponds to a situation with n c simultaneously open channels at the release site; P ( w | n c , [ R ] T ) is the probability that, for n c open channels in a cluster and a given [ R ] T , the resulting event induces the release of Ca 2 + from a neighboring cluster (i.e., generates a wave). P ( | n c , [ R ] T ) and P ( w | n c , [ R ] T ) , on the other hand, are the probabilities that a Ca 2 + release event obtained for given [ R ] T remains localized or initiates a wave, respectively. They satisfy: P ( | n c , [ R ] T ) + P ( w | n c , [ R ] T ) = 1 . We rewrite the latter as:
P ( w | [ R ] T ) = n c 1 P ( w | n c , [ R ] T ) P ( n c ) .
Combining Equations (1)–(3), we arrive at:
P ( n c ) = P ( n c | , [ R ] T ) 1 P ( w | [ R ] T ) + P ( w | n c , [ R ] T ) P ( n c ) .
Assuming that the Ca 2 + current through an open IP 3 R is approximately the same for all IP 3 Rs, we conclude that n c is proportional to the Ca 2 + current that underlies a Ca 2 + release event. The quantity A l i b F that we derive from the experimental data, on the other hand, is an increasing function of the underlying Ca 2 + current. We must point out that, if n c is large enough, puff amplitudes increase sublinearly with n c [15]. Assuming that n c and A l i b F are approximately linearly related, we can then use the experimental CDF of A l i b F , which we compute for the localized Ca 2 + release events, to estimate the CDF, F that can be computed from P ( n c | , [ R ] T ) :
F ( n , | , [ R ] T ) = n c = 1 n P ( n c | , [ R ] T ) .
The aim is to compare the distribution functions, F ( n , | , [ R ] T ) , for different values of [ R ] T using the corresponding experimental CDFs of A l i b F . In particular, we will compare the CDFs that are sufficiently different according to the K–S test: the ones with [ R ] T = 0 and with [ R ] T = 90 μ M . In what follows, we will drop the concentration units ( μ M ) from the expressions of the probabilities to simplify the notation. Defining P w ( n c ) P ( w | n c , 0 ) P ( w | n c , 90 ) and P w n c 1 P w ( n c ) P ( n c ) and using Equations (4) and (5), we obtain:
P w ( n c ) P ( n c ) = P ( n c | , 90 ) P ( n c | , 0 ) ( 1 P ( w | 90 ) ) + P ( n c | , 0 ) P w .
As illustrated in Figure 1, the experiments show that the difference between the two CDFs is more noticeable in the region of the largest size events, i.e., for the largest values of n c . We then compute:
n c n M P w ( n c ) P ( n c ) = F ( n M | , 0 ) F ( n M | , 90 ) ( 1 P ( w | 90 ) ) + P w ( 1 F ( n M | , 0 ) ) ,
where n M is the event size beyond which the CDFs start to differ more noticeably. As described later, the CDFs in the r.h.s. of Equation (8) can be estimated from the experimental CDFs. On the other hand, we estimate P ( w | n c , [ R ] T ) using the numerical simulations that we describe in the following section. Varying the parameters of the simulation, we determine the values for which we obtain estimates of the l.h.s. of this equation that are consistent with those of the r.h.s.

2.3. Numerical Simulations to Estimate the Probability That a Release Event from One (Primary) Cluster Induces the Release of Ca 2 + from Another (Secondary) Cluster

We compute the probability, P 0 ( t , d , n c , n s , [ R ] T ) , that n s IP 3 R located at a distance, d, from a (primary) cluster with n c IP 3 Rs that are simultaneously open at t = 0 , becomes open by a time, t. We want to compare how P 0 varies as [ R ] T is changed. We thus write explicitly its dependence on this variable. To compute P 0 ( t , d , n c , n s , [ R ] T ) , we proceed as explained in Materials and Methods (see also [10]) and the parameter values used are listed in Table 2. We show in Figure 2 the results obtained with n s = 1 . We show in Figs. Figure 2a–c the results obtained using the basal Ca 2 + concentration, [Ca] b = 0.1 μ M , for the conditions of Sets I, II and III. We show the results obtained at d = 0.6 μ m in Figure 2a and at d = 1.4 μ m (a typical inter-cluster distance) in Figure 2b,c. The number of simultaneously open channels is n c = 10 in Figure 2a,b and n c = 50 in Figure 2c. The change of P 0 with varying [ R ] T is unobservable for n c = 10 at d = 1.4 μ m (the difference is 0.004 for the times displayed in the figure) while it can be 0.085 at d = 0.6 μ m . Furthermore, it is P 0 ( t ) = P 0 ( t , d = 0.6 μ m , n c = 10 , [ R ] T = 90 μ M ) P 0 ( t , d = 0.6 μ m , n c = 10 , [ R ] T = 0 μ M ) 0.065 at t = d / V with V 10 μ ms 1 , a typical wave velocity. The maximum difference max t P 0 ( t ) increases with n c . This is shown in Figure 2c where n c = 50 , d = 1.4 μ m and max t P 0 ( t ) 0.026 . In Figure 2d, we show what happens when [Ca] b decreases. In this case, we compare P 0 at a distance d = 1.4 μ m from the source obtained for simulations performed with the concentrations of Set III (dashed line) and Set I (solid line) but with a different value of [Ca] b in each one (100 nM and 50 nM, respectively). A similar behaviour is obtained with n s = 5 (data not shown).
We interpret the [ R ] T -dependent changes of P 0 that are illustrated in Figure 2 as a sign of the change in the level of inter-cluster coupling (or, equivalently, disruption) that can be reached as the slow dye concentration is varied.
We further studied how sensitive is P 0 to changes in the distance to the source d. The results of Figure 3 are obtained using the basal Ca 2 + concentration, [Ca] b = 0.1 μ M . To illustrate the disruption when [ R ] T is increased, we show in Figure 3a P 0 = P 0 ( [ R ] T = 0 ) P 0 ( [ R ] T = 90 μ M ) computed with n c = 10 and n s = 1 as a function of the distance d for each time t = d / V (with V = 10 μ m / s ) and it can be observed that the probability of opening one IP 3 R in-between clusters decreases (from 0.4 μ m to 1.5 μ m ). When observing the probability of opening one IP 3 R as a function of d without adding the slow buffer ( R T = 0 , Set III) (solid line in Figure 3b), as d increases, this probability approximates to the basal probability (dotted line, P 0 computed as in Equation (14) but with no calcium dyes), almost no coupling can occur at the typical inter-cluster distance ( d = 1.4 μ m ). Thus, to explain the inter-cluster coupling, it is necessary to add a non-cluster IP 3 R in-between them. The optimal value of the parameter d should be on the order of 0.4–0.8 μ m (approximately the half distance between clusters). Not even adding n s = 5 sensing channels at the second cluster, the probability differs from the basal (dashed and dotted lines in Figure 3c, respectively). We choose d = 0.6 μ m to add an isolated IP 3 R in-between clusters (solid line in Figure 3c) and now the signal can propagate.
We now study whether the variations of Figure 2 can explain the changes in the distributions of Ca 2 + release during localized events observed in the experiments that are apparent in Figure 1.

2.4. Combining the Estimates Derived from the Experiments and from the Numerical Simulations to Interpret the Changes Observed Experimentally

We first estimate the r.h.s. of Equation (8) assuming that F ( n , | , [ R ] T ) is given by the experimental CDF of A l i b F for the same [ R ] T and some unknown factor between n and A l i b F . We recall here that the experimental CDF corresponds only to localized events (i.e., events at a primary cluster that do not induce the opening of IP 3 Rs at another secondary IP 3 R cluster). In order to estimate this unknown factor, we associate, n M , (the value after which the differences in the CDFs become more noticeable) to a value, A l i b F , for which F ( n M | , 0 ) is sufficiently close to 1. The basic assumption here is that, for [ R ] T = 0 , almost all primary events with n n M simultaneously open IP 3 Rs initiate waves in which case they do not remain localized and are, therefore, not included in the computation of the CDF. We choose A l i b F = 4.8 for which, according to the experimental data, it is F ( n M | , 0 ) 0.98 (see Figure 1). For this value, it is F n M | , 90 ) 0.90 . Thus, we estimate F ( n M | , 0 ) F ( n M | , 90 ) 0.08 . We do not have a direct estimation of P ( w | [ R ] T ) . Assuming that n c = 10 is the most probable value for the number of simultaneously open IP 3 Rs in a cluster, we approximate P ( w | [ R ] T ) = n c 1 P ( w | n c , [ R ] T ) P ( n c ) P ( w | n c = 10 , [ R ] T ) P 0 ( t = d / V , d , n c = 10 , [ R ] T ) with P 0 the open probability computed numerically that we introduced in the previous section and V a typical Ca 2 + wave velocity. For V, we try two values, V = 20 μ m / s and V = 10 μ m / s . For d, we try the typical inter-cluster distance, d = 1.4 μ m and the closer distance, d = 0.6 μ m that was probed in the previous section. Using d = 1.4 μ m and V = 20 μ m / s , the simulations give P ( w | 90 ) 0.14 and P w 0 . The estimate of the r.h.s. of Equation (8) then results equal to 0.07 . This value changes to 0.06 if we use V = 10 μ m / s . Using d = 0.6 μ m and V = 20 μ m / s the simulations give P ( w | 90 ) 0.12 and P w = 0.028 . The estimate of the r.h.s. of Equation (8) then results as 0.07 . This value changes to 0.06 if we use V = 10 μ m / s .
We now use the simulations of the previous section to put an upper bound, P w , m a x , on P w ( n c ) in the l.h.s of Equation (8). With such an upper bound, we can write n c n M P w ( n c ) P ( n c ) P w , m a x n c n M P ( n c ) = P w , m a x ( 1 F ( n M ) ) where F ( n M ) = n c n M P ( n c ) . Similarly to the way we have followed estimating P ( w | [ R ] T ) , we compute P ( w | n c , [ R ] T ) P 0 ( t = d / V , d , n c , [ R ] T ) , with P 0 the open probability of the previous section illustrated in Figure 2. In particular, using the results of these simulations, we conclude that P w ( n c ) P ( w | n c , 0 ) P ( w | n c , 90 μ M ) is larger the larger the value of n c . Thus, we obtain the upper bound, P w , m a x , using similar simulations to those of Figure 2 but for n c = 50 (a very large number of simultaneously open IP 3 Rs). Namely, we estimate P w , m a x P 0 ( t = d / V , d , n c = 50 , 0 ) P 0 ( t = d / V , d , n c = 50 , 90 ) . In order to put an upper bound on the l.h.s. of Equation (8), we need a bound for F ( n M ) , the CDF of all the (primary) event sizes at n = n M . As already explained, we assume that F ( n ) = n c n P ( n c ) does not depend on [ R ] T . Given our interpretation of the results, we assume that the difference between the CDF of all the (primary) event sizes, F ( n ) , and the CDF of the primary event sizes that remain localized for a given value of [ R ] T , F n | , [ R ] T , is due to the existence of primary events (of large enough size) that initiate waves for that value of [ R ] T . Assuming that the fraction of primary events that initiate waves for [ R ] T = 90 μ M is negligible, we can approximate F ( n M ) F n M | , 90 0.90 . If we do not want to use this approximation, then we can use the bound F ( n M ) < 0.90 . In what follows, we mostly use F ( n M ) = 0.90 , but we repeat some computations changing it to 0.80 to see how much the estimates could change. Proceeding as just explained, for d = 1.4 μ m , we obtain P w ( n c ) P w , m a x 0.003 for V = 20 μ m / s and P w , m a x 0.026 for V = 10 μ m / s . Using F ( n M ) = 0.9 , we then obtain 0.0003 and 0.0026 as upper bounds of the l.h.s. of Equation (8) for V = 20 μ m / s and V = 10 μ m / s , respectively. These two upper bounds are at least one order of magnitude smaller than the values obtained for the r.h.s. of Equation (8). If we use F ( n M ) = 0.8 to compute the l.h.s. of this equation, the latter estimate doubles with respect to the previous one. Thus, the order of magnitude difference between the left and right estimates for d = 1.4 μ m remains the same. Repeating the computations for d = 0.6 μ m , we obtain P w ( n c ) P w , m a x 0.19 for V = 20 μ m / s and P w , m a x 0.15 for V = 10 μ m / s . Using F ( n M ) = 0.9 , we then get 0.019 and 0.015 as upper bounds of the l.h.s. of Equation (8) for V = 20 μ m / s and V = 10 μ m / s , respectively. In this case, the values of the left- and right-hand sides are of the same order of magnitude. These estimates come closer together if we use F ( n M ) = 0.8 in the l.h.s. of the equation. In such a case, we obtain 0.04 and 0.03 for the l.h.s. estimate using V = 20 μ m / s and V = 10 μ m / s , respectively, two values that are pretty similar to the r.h.s. estimates, 0.07 and 0.06.

2.5. Changes in Basal Calcium Concentration, [Ca] b

As illustrated in Figure 2d, decreasing basal [Ca] with increasing [ R ] T changes the open probability at the distance, d = 1.4 μ m , in the direction that is needed to explain the observed changes in the event size distributions. We analyze here whether there is any evidence of a decreasing basal Ca 2 + with increasing [ R ] T in the experimental data. We show in Figure 4 the cumulative density functions of the mean basal fluorescence emitted by the Ca 2 + -bound Fluo-4 molecules, f 0 , F , in (a) and of the mean basal fluorescence emitted by the Ca 2 + -bound Rhod-2 molecules, f 0 , R in (b) for the experiments with [ R ] T = 36 μ M (dotted line) and with [ R ] T = 90 μ M (solid line). In Figure 4b, we rescaled f 0 R by 90 / 36 = 2.5 in the case of Set I to make the distributions of experiments of Set I (which has [ R ] T = 90 μ M ) and II (which has [ R ] T = 36 μ M ) readily comparable. The values of f 0 , F and f 0 , R were derived from the fluorescence observations as explained in Materials and Methods. We observe that the CDFs move to smaller values of their arguments with increasing [ R ] T . As the mean basal fluorescence is an increasing function of [Ca] b (see Equation (12)), this observation supports the idea that, on average, [Ca] b decreases with increasing Rhod-2.
In order to estimate the variation in [Ca] b with increasing [ R ] T , we compare f 0 , D ¯ for sets I and II. In particular, we obtain f 0 , F ¯ = 6.1 a.u., f 0 , R ¯ = 13.2 a.u. and f 0 , F ¯ = 7.2 a.u., f 0 , R ¯ = 6.4 a.u. for sets I and II, respectively. Inserting these values into Equation (13), using that N R = 32 for set II and N R = 80 for set I, and assuming that [Ca] b = 100 nM for set II, we obtain [Ca] b = 60–80 nM for set I, depending on whether we use the mean Rhod-2 or mean Fluo-4 basal fluorescence values.
We now repeat the calculations of the previous section but using the values of P 0 prescribed by the simulations with [Ca] b = 100 nM for set II and [Ca] b = 50 nM for set I (Figure 2d). In this case, the r.h.s. estimates do not change much from the previous calculations. The l.h.s. estimates, on the other hand, change slightly coming closer together with the r.h.s. estimates. For example, using V = 20 μ m / s , we obtain l.h.s. 0.023 for d = 0.6 μ m and l.h.s. 0.007 for d = 1.4 μ m . If we use V = 10 μ m / s , we obtain l.h.s. 0.019 for d = 0.6 μ m and l.h.s. 0.013 for d = 1.4 μ m .

3. Discussion

Intracellular Ca 2 + signals are ubiquitous across cell types. The spatial range over which they spread is key to determining the eventual end responses. This points to the importance of understanding how intracellular Ca 2 + signals propagate inside the cells. To this end, Ca 2 + release from the endoplasmic reticulum into the cytosol through IP 3 Rs plays a major role. IP 3 Rs are Ca 2 + channels whose open probability depends on the cytosolic Ca 2 + concentration. Therefore, the Ca 2 + released through an open IP 3 R can induce the opening of nearby IP 3 Rs. High Ca 2 + concentrations, on the other hand, lead to the inhibition of the channels. This dynamic is clearly excitable. In this regard, however, the excitability of the cytosol is “patchy”: IP 3 Rs tend to be organized in clusters separated by 1.4 μ m . This may lead to propagation failure when the Ca 2 + released from one IP 3 R reaches the vicinity of another one at a concentration that is not enough to induce its opening. The inter-cluster Ca 2 + -mediated coupling can be interfered by means of Ca 2 + buffers. This is used as an experimental tool, but the cells can do it as well.
In [10], we presented experimental results in which we studied how the presence of competing Ca 2 + trapping mechanisms of different kinetics altered the resulting intracellular signals. Differently from previous studies [12], in [10], we made visible the invisible by using two dyes of different kinetics as the Ca 2 + trapping mechanisms. The work of [10] not only allowed us to draw conclusions on how the signals were reshaped by the presence of the different buffers but also gave some indications on the spatial distribution of the IP 3 Rs involved in the signals. In particular, based on a qualitative comparison between experiments and numerical simulations, we concluded in [10] that Ca 2 + release seemed to occur not only from the clearly identifiable release sites (IP 3 R clusters), but also from some functional, probably isolated, IP 3 Rs in between them. In this paper, we have presented a quantitative analysis of the experiments of [10] that corroborated this conclusion.
For the quantitative comparison between experiments and models, in this paper, we have focused on the size distribution of the localized Ca 2 + release events (puffs) that were obtained with the experiments of [10] for [ E G T A ] = 90 μ M , [ F ] T = 36 μ M and two concentrations of the slow dye Rhod-2, [ R ] T = 36 μ M and 90 μ M . Given that the fluorescence amplitude observed for a given release event could change with varying [ R ] T even if the underlying Ca 2 + current remained the same, we characterized the observed puffs by the quantity, A l i b F (Equation (11)) that we introduced in [10] to overcome this problem. The A l i b F distributions obtained for the analyzed experiments showed a shift towards larger values of A l i b F as [ R ] T was increased (see Figure 1). This shift agrees with previous observations and analyses according to which the increase of a slow Ca 2 + buffer concentration (in this case, Rhod-2) disrupts the Ca 2 + -mediated coupling between clusters (Figure 3a) [10,12,14,16]. Namely, we interpret this shift as reflecting the fact that events that are characterized by a certain number of simultaneously open IP 3 Rs at a primary cluster and remain localized for a given value of [ R ] T can induce the opening of IP 3 Rs at other (secondary) clusters for smaller values of [ R ] T . While the former events are puffs and would then be considered for the computation of the localized event size distribution, the latter would not because they correspond to waves.
We observed that clusters can become coupled when adding an IP 3 R in-between them (Figure 3c). We introduced a probabilistic model in order to analyze quantitatively whether the differences observed for the experiments performed for [ R ] T = 0 and [ R ] T = 90 μ M could be explained if Ca 2 + release only occurred through IP 3 R-clusters separated by 1.4 μ m or not. Within the framework of the probabilistic model, we then combined the analysis of the experimental data with some probability estimates derived from numerical simulations similar to those presented in [10]. We determined in this way that the numerically estimated values were not compatible with the differences observed experimentally if the only Ca 2 + release sites involved were 1.4 μ m apart from one another. The experimental and numerical results were more compatible if we assumed that there was Ca 2 + release from at least one IP 3 R at a distance 0.6 μ m from the primary Ca 2 + release cluster. The presence of the slow dye, on the other hand, could reduce the basal Ca 2 + concentration. We analyzed that possibility in the experimental data (Figure 4) and estimated that [Ca] b could have been reduced by half when [ R ] T was changed from 0 to 90 μ M . The numerical simulations, on the other hand, showed that a decreasing value of [Ca] b with increasing [ R ] T gave better results for d = 1.4 μ m (Figure 2d) in terms of their compatibility with their experiments. We then re-analyzed the experimental data but using numerical simulations that included this change in [Ca] b with varying [ R ] T . The best situation to explain Figure 1 was obtained with simulations that combined a change in [Ca] b with [ R ] T and the presence of a functional IP 3 R at a shorter distance ( d 0.6 μ m ) than the typical inter-cluster one.
Our quantitative analysis of the experiments of [10] presented in this paper confirms that the spatial landscape over which intracellular Ca 2 + signals propagate do not consist solely of patches of excitability that are 1.4–2 μ m apart from one another but that there are also “relay stations” (isolated functional IP 3 Rs) in between. Probably, the existence of these in-between IP 3 Rs is necessary for the propagation of Ca 2 + waves.

4. Materials and Methods

4.1. Oocyte Preparation

Experiments were performed on Xenopus laevis immature oocytes previously treated with collagenase. Oocytes were loaded by intracellular microinjection with different compounds. Two calcium dyes Fluo-4 dextran high affinity ( K d = 0.8 μ M ) and Rhod-2 ( K d = 2 μ M ) were used to probe cytosolic [Ca]. Caged InsP 3 (D-Myo-Inositol 1,4,5-Triphosphate,P4(5)-(1-(2-Nitrophenyl)ethyl) Ester) was used to induce IP 3 R opening. The exogenous Ca 2 + buffer EGTA was also used. Final intracellular concentrations of the different compounds were calculated assuming a 1 μ L cytosolic volume. Final intracellular concentration of InsP 3 was 9 μ M in all of the experiments. The different concentrations used in each experiment are detailed in Table 1 where we classify the experiments in three sets. Fluo-4, Rhod-2 and InsP 3 were from Molecular Probes Inc.; EGTA was from Sigma Aldrich. Recordings were made at room temperature.

4.2. Confocal Microscopy

Confocal imaging was performed using a spectral confocal scanning microscope Olympus FluoView1000 that has a spectral scan unit connected to an inverted microscope IX81. The caged compound was photolyzed with the UV part of the spectrum of a mercury lamp that comes with the microscope using the modification introduced in [17]. Fluo-4 was excited with the 488 nm line of a multiline Argon laser, Rhod-2 was excited using the 543 nm line of a He–Ne laser. Both lasers were focused on the oocyte with a 60 × oil immersion objective (NA 1.35). The Fluo-4 and Rhod-2 emitted fluorescences were simultaneously detected in the 500–600 nm and the 600–630 ranges, respectively, with PMT detectors. All the experiments were performed in the linescan imaging mode to improve the temporal resolution. Linescan images were obtained by scanning along a fixed line (250 px) within the oocyte. The acquisition rate was fixed at 10 μ s per pixel resulting in a scan rate of 3.26 ms per line. The caged compound was photo-released approximately 3 s after the linescan acquisition started.

4.3. Image Analysis

All images were analyzed using routines written in MATLAB. In the experiments where we simultaneously acquired the fluorescence coming from two channels (around 510 nm for Fluo-4 and 570 nm for Rhod-2), we used a linear unmmixing method to minimize the effect of the spectral bleed-through ( R = 0.1626 was the linear unmixing coefficient used). The images were also smoothed by averaging over the eight nearest pixels.
The events were identified and the images were processed as explained in [10]. From the fluorescence distribution, f D ( x i , t j ) , collected in each of the channels ( D = R for Rhod-2 and D = F for Fluo-4) at each pixel, ( x i , t j ) , and the relative increase in fluorescence at the peak of the signal, f r , D = max x i , t j ( ( f r , D ( x i , t j ) f 0 , D ( x i ) ) / f 0 , D ( x i ) ) , with f 0 , D ( x i ) the mean basal fluorescence at x i observed with D before the UV flash, we computed the corresponding Ca 2 + -bound dye and (maximum) relative Ca 2 + bound dye concentrations ( [ C a D ] and [ C a D ] r max ( [ C a D ] [ C a D ] b ) / [ C a D ] b , respectively, with [ C a D ] b the basal Ca 2 + -bound dye concentration). To estimate the Ca 2 + -bound dye concentration, we followed [18] neglecting fluctuations in the number of dye molecules that contribute to the fluorescence at each pixel, N D , (for more details, see [10]):
[ C a D ] = [ D ] T q 1 , D q 2 , D f D γ N D q 2 , D , D = R , F ,
[ C a D ] r = f r , D 1 + q 2 , D / q 1 , D ( 1 q 2 , D q 1 , D ) [ C a ] b [ C a ] b + K d , D , D = R , F .
To compute these quantities we followed [18] and used N F = 32 , N R = 32 for [ R ] T = 36 μ M , N R = 80 for [ R ] T = 90 μ M , [Ca] b = 100 nM , q 1 , F = 0.45 , q 2 , F = 0.01 , K d , F = 0.8 μ M , q 1 , R = 0.36 , q 2 , R = 0.02 and K d , R = 2 μ M .
In this paper, we only analyzed the events observed in the Fluo-4 channel, i.e., for D = F . The event size of each analyzed puff was then characterized by the maximum value of the relative increase in the Ca 2 + -bound Fluo-4 concentration, [ C a F ] r that we derived from the observed fluorescence. As done in [10], we then used the total Rhod-2 concentration, [ R ] T , of the experiment to obtain estimates of the maximum values, [ C a F ] r ( [ R ] = 0 ) , that would have been attained for the same release event if only the dye, Fluo-4, had been present. As discussed in [10], this estimate that we call A l i b F is an increasing function of the Ca 2 + current that underlies the release event regardless of the value of [ R ] T used if the Ca 2 + current arises from a very localized spatial region (the cluster). As done in [10], we computed it as:
A l i b F [ C a F ] r + α F , R [ R ] T ,
with α F , R = 4.58 × 10 3 .

4.4. Basal Calcium Estimation

In order to study the behavior of the mean basal Ca 2 + concentration, [ C a ] b , for each experiment type probed in the paper, we follow some of the steps of the method introduced in [18]. We work with linescan images obtained before any UV flash has been applied, i.e., we analyze basal fluorescence. On these images, we get rid of the horizontal lines that are persistently dark, which correspond to the cortical granules. We then compute the mean basal fluorescence for each linescan image as:
f 0 , D = 1 N i b f j = 1 j U V f D ( x i , t j ) , D = R , F ,
where the sum over i runs over the horizontal lines that are not persistently dark and the subscript, D, denotes whether the fluorescence comes from the Fluo-4 ( D = F ) or Rhod-2 ( D = R ) molecules. Using the values, f 0 , D , obtained for each experiment type, we compute the corresponding cumulative distribution functions of the mean basal fluorescence. To transform from basal fluorescence to basal Ca 2 + , we use the following expression derived from [18]:
f 0 , D ¯ = γ D ( q 1 , D q 2 , D ) [ C a ] b [ C a ] b + K d , D + q 2 , D N D ¯ ,
which takes into account the contributions to the fluorescence from the free and the Ca 2 + -bound dye molecules with brightness q 2 , D and q 1 , D , respectively. In Equation (13) K d , D is the dissociation constant of the Ca 2 + -dye reaction, N D ¯ is the mean number of dye molecules that contribute to the fluorescence collected at a pixel and γ D is a multiplying factor introduced by the detector ( γ R = 6 and γ F = 5 [18]).

4.5. Numerical Simulations

To assess the rate of CICR-mediated coupling between neighboring clusters, we compute the probability that an IP 3 R that n s IP 3 R located at a distance, d, from a Ca 2 + point source becomes open during a time interval, t , since the start of the release by means of:
P 0 ( t , d , n c , [ R ] T ) = 1 exp 0 t k o n n s [ C a 2 + ] ( d , t ) d t
with k o n = 20 μ M 1 s 1 the rate of Ca 2 + binding to the activating site of an IP 3 R of the DeYoung–Keizer model [19]. We compute [ C a 2 + ] ( d , t ) simulating a set of reaction-diffusion equations in a spherical volume (assuming spherical symmetry with r the radial coordinate) for: Ca 2 + , an immobile endogenous buffer (S), two cytosolic indicators (F and R) and an exogenous mobile buffer ( E G T A ). A point source located at the origin and pumps (P) that remove Ca 2 + uniformly in space are also included. For the source, we assume that it consists of n c channels that open simultaneously at t = 0 , each of which becomes close after a time that is drawn from an exponential distribution of mean t o p e n = 20 ms [19]. For the Ca 2 + -buffer or dye reactions we consider that a single Ca 2 + ion binds to a single buffer or dye molecule (X) according to:
C a 2 + + X k o n X k o f f X [ C a X ] ,
where X represents F, R, E G T A , or S and k o n X and k o f f X are the forward and backward binding rate constants of the corresponding reaction, respectively. We assume that the total concentrations of dyes and buffers ( [ F ] T , [ R ] T , [ E G T A ] T , and [ S ] T ) are spatially uniform at t = 0 so that they remain uniform and constant for all times. We also assume that [Ca 2 + ] is initially uniform, equal to its basal value and in equilibrium with the buffers and dyes. The parameter values used are listed in Table 3.

Author Contributions

E.P. conducted the experiments. E.P. and S.P.D. performed the simulations. E.P. and S.P.D. analyzed the experiments and simulations and wrote the paper. S.P.D. conceived the work.

Funding

This research has been supported by UBA (UBACyT 20020170100482BA) and ANPCyT (PICT 2015-3824).

Conflicts of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

  1. Berridge, M.J.; Bootman, M.D.; Lipp, P. Calcium—A life and death signal. Nature 1998, 395, 645–648. [Google Scholar] [CrossRef] [PubMed]
  2. Bootman, M.D.; Collins, T.J.; Peppiatt, C.M.; Prothero, L.S.; MacKenzie, L.; Smet, P.D.; Travers, M.; Tovey, S.C.; Seo, J.T.; Berridge, M.J.; et al. Calcium signalling—An overview. Semin. Cell Dev. Biol. 2001, 12, 3–10. [Google Scholar] [CrossRef] [PubMed]
  3. Choe, C.U.; Ehrlich, B.E. The Inositol 1,4,5-Trisphosphate Receptor (IP3R) and Its Regulators: Sometimes Good and Sometimes Bad Teamwork. Sci. Signal. 2006, 2006, re15. [Google Scholar] [CrossRef] [PubMed]
  4. Foskett, J.K.; White, C.; Cheung, K.H.; Mak, D.O.D. Inositol Trisphosphate Receptor Ca2+ Release Channels. Physiol. Rev. 2007, 87, 593–658. [Google Scholar] [CrossRef] [PubMed]
  5. Fabiato, A. Calcium-induced release of calcium from the cardiac sarcoplasmic reticulum. Am. J. Physiol. 1983, 245, 1–15. [Google Scholar] [CrossRef] [PubMed]
  6. Sun, X.P.; Callamara, N.; Marchant, J.S.; Parker, I. A continuum of InsP3-mediated elementary Ca2+ signalling events in Xenopus oocyte. J. Physiol. 1998, 509, 67–80. [Google Scholar] [CrossRef] [PubMed]
  7. Smith, I.F.; Parker, I. Imaging the quantal substructure of single IP3R channel activity during Ca2+ puffs in intact mammalian cells. Proc. Natl. Acad. Sci. USA 2009, 106, 6404–6409. [Google Scholar] [CrossRef] [PubMed]
  8. Solovey, G.; Dawson, S.P. Intra-cluster percolation of calcium signals. PLoS ONE 2010, 5, 1–8. [Google Scholar] [CrossRef] [PubMed]
  9. Keizer, J.; Smith, G.D.; Ponce-Dawson, S.; Pearson, J.E. Saltatory Propagation of Ca2+ Waves by Ca2+ Sparks. Biophys. J. 1998, 75, 595–600. [Google Scholar] [CrossRef]
  10. Piegari, E.; Lopez, L.F.; Dawson, S.P. Using two dyes to observe the competition of Ca2+ trapping mechanisms and their effect on intracellular Ca2+ signals. Phys. Biol. 2018, 15, 066006. [Google Scholar] [CrossRef] [PubMed]
  11. Paredes, R.M.; Etzler, J.C.; Watts, L.T.; Zheng, W.; Lechleiter, J.D. Chemical calcium indicators. Methods 2008, 46, 143–151. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Dargan, S.L.; Parker, I. Buffer kinetics shape the spatiotemporal patterns of IP3-evoked Ca2+ signals. J. Physiol. 2003, 553, 775–788. [Google Scholar] [CrossRef] [PubMed]
  13. Dargan, S.L.; Schwaller, B.; Parker, I. Spatiotemporal patterning of IP3-mediated Ca2+ signals in Xenopus oocytes by Ca2+-binding proteins. J. Physiol. 2004, 556, 447–461. [Google Scholar] [CrossRef] [PubMed]
  14. Piegari, E.; Sigaut, L.; Ponce Dawson, S. Ca2+ images obtained in different experimental conditions shed light on the spatial distribution of IP3 receptors that underlie Ca2+ puffs. Cell Calcium 2015, 57, 109–119. [Google Scholar] [CrossRef] [PubMed]
  15. Solovey, G.; Fraiman, D.; Dawson, S.P. Mean field strategies induce unrealistic nonlinearities in calcium puffs. Front. Physiol. 2011, 2, 1–11. [Google Scholar] [CrossRef]
  16. Callamaras, N.; Parker, I. Phasic characteristic of elementary Ca2+ release sites underlies quantal responses to IP3. EMBO J. 2000, 19, 3608–3617. [Google Scholar] [CrossRef] [PubMed]
  17. Sigaut, L.; Barella, M.; Espada, R.; Ponce, M.L.; Dawson, S.P. Custom-made modification of a commercial confocal microscope to photolyze caged compounds using the conventional illumination module and its application to the observation of Inositol 1,4,5-trisphosphate-mediated calcium signals. J. Biomed. Opt. 2011, 16, 066013. [Google Scholar] [CrossRef] [PubMed]
  18. Piegari, E.; Lopez, L.; Perez Ipiña, E.; Ponce Dawson, S. Fluorescence fluctuations and equivalence classes of Ca2+ imaging experiments. PLoS ONE 2014, 9, e95860. [Google Scholar] [CrossRef] [PubMed]
  19. Youngt, G.W.D.E.; Keizer, J. A single pool. Nature 1992, 89, 9895–9899. [Google Scholar]
Figure 1. CDF of event sizes (as measured by A l i b F ) for experiments performed with [ F ] T = 36 μ M and [ E G T A ] = 90 μ M and different values of [ R ] T (the solid line corresponds to Set I, the dotted line to Set II and the dashed line to Set III (see Table 1).
Figure 1. CDF of event sizes (as measured by A l i b F ) for experiments performed with [ F ] T = 36 μ M and [ E G T A ] = 90 μ M and different values of [ R ] T (the solid line corresponds to Set I, the dotted line to Set II and the dashed line to Set III (see Table 1).
Mca 24 00061 g001
Figure 2. P 0 dependence with distance d from the source and [Ca] b . (a) performed at d = 0.6 μ m with n c = 10 ; (b) at d = 1.4 μ m with n c = 10 ; (c) at d = 1.4 μ m and n c = 50 and (d) d = 1.4 μ m with n c = 10 while varying [Ca] b from 100 nM (Set III) to 50 nM (Set I). The dashed line corresponds to Set III, the dotted line corresponds to Set II, and the solid line corresponds to Set I.
Figure 2. P 0 dependence with distance d from the source and [Ca] b . (a) performed at d = 0.6 μ m with n c = 10 ; (b) at d = 1.4 μ m with n c = 10 ; (c) at d = 1.4 μ m and n c = 50 and (d) d = 1.4 μ m with n c = 10 while varying [Ca] b from 100 nM (Set III) to 50 nM (Set I). The dashed line corresponds to Set III, the dotted line corresponds to Set II, and the solid line corresponds to Set I.
Mca 24 00061 g002
Figure 3. Existence of IP 3 Rs in-between clusters is necessary to explain the observations. (a) P 0 = P 0 ( [ R ] T = 0 ) P 0 ( [ R ] T = 90 μ M ) with n c = 10 and n s = 1 as a function of the distance d for each time t = d / V (with V = 10 μ m / s ); (b) P 0 computed as in (a) for the Set III (solid line) and basal (dotted line) conditions; (c) P 0 computed in the condition of Set III with n c = 10 and n s = 5 at a distance d = 1.4 μ m from the Ca 2 + point source as a function of time when an IP 3 R is added in-between cluster (solid line), with no additional IP 3 R (dashed line) and at the basal condition (dotted lined).
Figure 3. Existence of IP 3 Rs in-between clusters is necessary to explain the observations. (a) P 0 = P 0 ( [ R ] T = 0 ) P 0 ( [ R ] T = 90 μ M ) with n c = 10 and n s = 1 as a function of the distance d for each time t = d / V (with V = 10 μ m / s ); (b) P 0 computed as in (a) for the Set III (solid line) and basal (dotted line) conditions; (c) P 0 computed in the condition of Set III with n c = 10 and n s = 5 at a distance d = 1.4 μ m from the Ca 2 + point source as a function of time when an IP 3 R is added in-between cluster (solid line), with no additional IP 3 R (dashed line) and at the basal condition (dotted lined).
Mca 24 00061 g003
Figure 4. CDF of basal fluorescence. f 0 , D for experiments with [ E G T A ] = 90 μ M are shown in the Fluo-4 ( D = F , (a)) and Rhod-2 ( D = R , (b)) channels. In both cases, Set II is depicted with dotted lines and Set I with solid lines. In (b), the values, f 0 , R , of Set I are divided by 90 / 36 = 2.5 to make both distributions comparable.
Figure 4. CDF of basal fluorescence. f 0 , D for experiments with [ E G T A ] = 90 μ M are shown in the Fluo-4 ( D = F , (a)) and Rhod-2 ( D = R , (b)) channels. In both cases, Set II is depicted with dotted lines and Set I with solid lines. In (b), the values, f 0 , R , of Set I are divided by 90 / 36 = 2.5 to make both distributions comparable.
Mca 24 00061 g004
Table 1. Combinations of the dyes and EGTA concentrations used in the different experiments.
Table 1. Combinations of the dyes and EGTA concentrations used in the different experiments.
Experiment[Fluo-4] ( μ M )[Rhod-2] ( μ M )[EGTA] ( μ M )
Set I369090
Set II363690
Set III36090
Table 2. Value of the parameters varied to compute P 0 .
Table 2. Value of the parameters varied to compute P 0 .
ParameterAbbreviatureValues
Number of IP 3 Rs in the source n c 1, 10, 50
Distance to the Ca 2 + sourced(0.4–1.5) μ m
Number of sensing IP 3 Rs n s 1, 5
Rhod-2 concentration [ R ] T 0, 90 μ M
Velocity of propagationV10, 20 μ m/s
Table 3. Parameter values used to solve the simulations.
Table 3. Parameter values used to solve the simulations.
ParameterValueUnits
Free Calcium
D C a 220 μ m 2 s 1
[ C a ] b 0.05–0.1 μ M
Calcium dye Fluo-4-dextran
D F 15 μ m 2 s 1
k o n F 240 μ M 1 s 1
k o f f F 180 s 1
[ F ] T 36 μ M
Calcium dye Rhod-2-dextran
D R 15 μ m 2 s 1
k o n R 70 μ M 1 s 1
k o f f R 130 s 1
[ R ] T 0, 36, 90 μ M
Exogenous buffer EGTA
D E G T A 80 μ m 2 s 1
k o n E G T A 5 μ M 1 s 1
k o f f E G T A 0,75 s 1
[ D ] T 90 μ M
Endogenous immobile buffer
D S 0 μ m 2 s 1
k o n S 400 μ M 1 s 1
k o f f S 800 s 1
[ S ] T 300 μ M
Pump
k p 0.1 s 1
v p 0.9 μ Ms 1
Source
n c 1, 10, 50-
t o p e n 20ms
I C a 0.1pA

Share and Cite

MDPI and ACS Style

Piegari, E.; Ponce Dawson, S. Functional Ca2+ Channels between Channel Clusters are Necessary for the Propagation of IP3R-Mediated Ca2+ Waves. Math. Comput. Appl. 2019, 24, 61. https://doi.org/10.3390/mca24020061

AMA Style

Piegari E, Ponce Dawson S. Functional Ca2+ Channels between Channel Clusters are Necessary for the Propagation of IP3R-Mediated Ca2+ Waves. Mathematical and Computational Applications. 2019; 24(2):61. https://doi.org/10.3390/mca24020061

Chicago/Turabian Style

Piegari, Estefanía, and Silvina Ponce Dawson. 2019. "Functional Ca2+ Channels between Channel Clusters are Necessary for the Propagation of IP3R-Mediated Ca2+ Waves" Mathematical and Computational Applications 24, no. 2: 61. https://doi.org/10.3390/mca24020061

APA Style

Piegari, E., & Ponce Dawson, S. (2019). Functional Ca2+ Channels between Channel Clusters are Necessary for the Propagation of IP3R-Mediated Ca2+ Waves. Mathematical and Computational Applications, 24(2), 61. https://doi.org/10.3390/mca24020061

Article Metrics

Back to TopTop