Fire – diffuse – fire model of dynamics of intracellular calcium waves

When Ca21 is released from internal stores in living cells, the resulting wave of increased concentration can travel without deformation (continuous propagation) or with burst-like behavior (saltatory propagation). We analyze the ‘‘fire–diffuse–fire’’ model in order to illuminate the differences between these two modes of propagation. We show that the Ca21 release wave in immature Xenopus oocytes and cardiac myocytes is saltatory, whereas the fertilization wave in the mature oocyte is continuous. Traveling waves in living cells can vary greatly in their appearance. For example, the calcium wave in immature Xenopus frog eggs propagates as a sequence of bursts (1–4) (Fig. 1a), whereas the calcium wave that occurs during fertilization in mature Xenopus eggs appears to be continuous (5, 6) (Fig. 1b). It is commonly believed that information is encoded in the time course of the Ca21 signal (7–12). Thus, the distinction between these two modes of propagation likely has physiological significance. Remarkably, the waves in both cell types involve the release of calcium from internal stores via mechanisms that are physically and biologically similar. We report here a simple physical model that allows us to characterize these behaviors as limiting cases. A prominent feature of this model is the existence of a new time scale in addition to the usual chemical time scales, namely the intersite diffusion time, d2yD. (Here, D is the diffusion coefficient and d is the source separation.) If intersite diffusion is rate-limiting, the wave will exhibit burst-like (or saltatory) behavior (Figs. 1a and 2a). If chemical processes are rate-limiting (i.e., slow compared to intersite diffusion), then the resulting wave will be smooth (Figs. 1b and 2b). Although the intersite diffusion time is likely to be important for other wave phenomena, our analysis of recent experiments suggests that it is an essential feature for intracellular Ca21 waves. Ca21 is stored intracellularly in the endoplasmic or sarcoplasmic reticulum at 2–3 orders of magnitude greater than its concentration in the cytosol and is released by a process referred to as Ca21-induced Ca21 release (CICR). CICR involves Ca21 release through specialized Ca21 channels that are activated at slightly elevated levels of cytosolic Ca21 and then inactivated as the level of Ca21 rises further. After an open channel closes via inactivation, it cannot reopen for some time, during which it is in a ‘‘refractory’’ state. Calcium release provides a mechanism for wave propagation (10–17), whereas inactivation limits the amount released. High-resolution imaging of Ca21 in a variety of cell types reveals localized release events called ‘‘puffs’’ (1–4) (Fig. 1a), ‘‘sparks’’ (19), or ‘‘quantum emission domains’’ (8, 9). These events have been correlated with Ca21 release from either individual or small clusters of channels. In muscle cells, sparks are observed at the t-tubule structures in the sarcoplasmic reticulum (20), which are aligned in regular arrays with a characteristic separation of ;2 mm. In other cell types, inositol trisphosphate stimulates Ca21 release from clusters of inositol triphosphate receptors from the endoplasmic reticulum. Although these inositol triphosphate-dependent sites appear to be randomly distributed, they are characterized by a mean separation (1–4, 21). The release events are the elementary processes underlying Ca21 waves. Our description of these events is called the ‘‘fire–diffuse– fire’’ model. This model reproduces the full range of wave propagation, from saltatory (Fig. 2a) to continuous (Fig. 2b). In contrast, homogeneous reaction diffusion models predict only continuous propagation with a velocity (v) proportional to =Dyt, with t the time scale of a chemical process. When applied to Ca21 waves, the fire–diffuse–fire model predicts that the fertilization wave in Xenopus eggs is continuous. Moreover, it shows, in agreement with experiment, that waves in immature Xenopus oocytes (1–4) and in cardiac myocytes are saltatory. The fire–diffuse–fire model consists of a regular array of point-source release sites with spacing, d, embedded in a continuum in which calcium ions diffuse. A related model, in which the discrete entities are individual cells, has been used for modeling cAMP (22) waves. In cells in which release sites do not occur in regular arrays, d should be thought of as a mean spacing. The pumps that resequester the ions back into the stores operate on a somewhat slower time scale (23) and are neglected here for simplicity. Whenever the cytosolic Ca21 concentration in the vicinity of a release site reaches a threshold value, [Ca]T, the site begins releasing Ca21 ions at a rate syt. It remains open for a time t, and then closes after having released a quantity, s, of Ca21. Because the binding time for Ca21 with immobile buffers is faster than the other time scales in the problem (15), we treat the action of the Ca21 buffers as instantaneous. To account for the effect of buffering, we define s as the total amount of Ca21 released divided by a buffering factor of ;100. We use the data in ref. 24 for the buffered diffusion coefficient in Xenopus oocytes. Note that this means that both s and the diffusion coefficient, D, inherit the temperature dependence of the buffering factor. The processes of Ca21 release (or firing) and diffusion can be described by a single evolution equation described in detail in the Appendix and in ref. 25. The dynamics depends on two dimensionless parameters, G and b. The first is defined as G 5 ~syd! @Ca#T 2 @Ca#b , [1] where [Ca]b is the basal Ca21 concentration. The quantity syd3 is the release concentration, since it is the amount of Ca21 released per site after buffering (s) divided by the volume per site (1yd3). So G is the ratio of the release concentration to the difference between basal and threshold concentrations. The parameter b is the ratio of the time a site remains open, t, to the intersite diffusion time, d2yD, i.e. The publication costs of this article were defrayed in part by page charge payment. This article must therefore be hereby marked ‘‘advertisement’’ in accordance with 18 U.S.C. §1734 solely to indicate this fact. PNAS is available online at www.pnas.org. †To whom reprint requests should be addressed. e-mail: pearson@ lanl.gov.

characteristic separation of ϳ2 m.In other cell types, inositol trisphosphate stimulates Ca 2ϩ release from clusters of inositol triphosphate receptors from the endoplasmic reticulum.Although these inositol triphosphate-dependent sites appear to be randomly distributed, they are characterized by a mean separation (1)(2)(3)(4)21).The release events are the elementary processes underlying Ca 2ϩ waves.
Our description of these events is called the ''fire-diffusefire'' model.This model reproduces the full range of wave propagation, from saltatory (Fig. 2a) to continuous (Fig. 2b).In contrast, homogeneous reaction diffusion models predict only continuous propagation with a velocity (v) proportional to ͌ D͞, with the time scale of a chemical process.When applied to Ca 2ϩ waves, the fire-diffuse-fire model predicts that the fertilization wave in Xenopus eggs is continuous.Moreover, it shows, in agreement with experiment, that waves in immature Xenopus oocytes (1)(2)(3)(4) and in cardiac myocytes are saltatory.
The fire-diffuse-fire model consists of a regular array of point-source release sites with spacing, d, embedded in a continuum in which calcium ions diffuse.A related model, in which the discrete entities are individual cells, has been used for modeling cAMP (22) waves.In cells in which release sites do not occur in regular arrays, d should be thought of as a mean spacing.The pumps that resequester the ions back into the stores operate on a somewhat slower time scale (23) and are neglected here for simplicity.Whenever the cytosolic Ca 2ϩ concentration in the vicinity of a release site reaches a threshold value, [Ca 2ϩ ] T , the site begins releasing Ca 2ϩ ions at a rate ͞.It remains open for a time , and then closes after having released a quantity, , of Ca 2ϩ .Because the binding time for Ca 2ϩ with immobile buffers is faster than the other time scales in the problem (15), we treat the action of the Ca 2ϩ buffers as instantaneous.To account for the effect of buffering, we define as the total amount of Ca 2ϩ released divided by a buffering factor of ϳ100.We use the data in ref. 24 for the buffered diffusion coefficient in Xenopus oocytes.Note that this means that both and the diffusion coefficient, D, inherit the temperature dependence of the buffering factor.
The processes of Ca 2ϩ release (or firing) and diffusion can be described by a single evolution equation described in detail in the Appendix and in ref. 25.The dynamics depends on two dimensionless parameters, ⌫ and ␤.The first is defined as where [Ca 2ϩ ] b is the basal Ca 2ϩ concentration.The quantity ͞d 3 is the release concentration, since it is the amount of Ca 2ϩ released per site after buffering () divided by the volume per site (1͞d 3 ).So ⌫ is the ratio of the release concentration to the difference between basal and threshold concentrations.The parameter ␤ is the ratio of the time a site remains open, , to the intersite diffusion time, d 2 ͞D, i.e.
If release is the rate-limiting process, then ␤ Ͼ Ͼ 1, and the waves are continuous.If intersite diffusion is the rate-limiting process, then ␤ Ͻ Ͻ 1, and the waves are saltatory (25).
A traveling saltatory wave solution (␤ ϭ 0.08 and ⌫ ϭ 3.05) to the fire-diffuse-fire model is illustrated in Fig. 2a.The vertical axis is the Ca 2ϩ concentration, and the horizontal axis is position in site spacing units.The parameters were estimated from the literature on the immature oocyte (1)(2)(3)(4)26).The wave is traveling to the right with a velocity of 20 m/s.Well behind the front (i.e., for x Ͻ ϳ2), the concentration is [Ca 2ϩ ](x, t) ϭ ͞d 3 ϩ [Ca 2ϩ ] b ϭ 0.66 M. The time at which site n begins to fire (i.e., to release Ca 2ϩ ) is denoted by t n .The first profile in the sequence is the solid line with the short spike above site 5.This profile occurs just after the Ca 2ϩ concentration at site 5 crosses threshold at time t ϭ t 5 .The second profile in the sequence is the long-dashed line with the tall peak above site 5.This is the profile at time t ϭ t 5 ϩ .The peak is wider because of diffusion and higher because the total amount () of Ca 2ϩ ions have been released at this time.The third profile is the short-dashed line and corresponds to a time halfway between t 5 and t 6 (after site 5 has completed firing, but before site 6 has reached threshold).The fourth profile is the solid line with the short spike above site 6.This is the profile just after the Ca 2ϩ concentration reaches threshold at time t ϭ t 6 .We see that saltatory propagation consists of a distinct emission of Ca 2ϩ when a particular site fires.Obviously, both Ca 2ϩ release and Ca 2ϩ diffusion occur simultaneously, but the firing time is so short compared to the intersite diffusion time, d 2 ͞D, that during the time that a single site is firing, very few ions diffuse to the next site.The shorthand description of these processes, fire-diffuse-fire, is thus an accurate description of the sequence of events shown in Fig. 2a.
Continuous propagation (␤ ϭ 50 and ⌫ ϭ 5) is illustrated in Fig. 2b.The parameters in this case are appropriate for the fertilization wave in mature Xenopus eggs (5).For D, we used a higher value (D ϭ 50 m 2 ͞s because the Ca 2ϩ concentration is higher, thus implying less buffering in the egg (24).The wave is traveling to the right with a velocity of 5.2 m͞s.The experimental observations provide the concentration in the wave back, from which it is possible only to estimate ͞d 3 but not and d separately.Here, D ϭ 450 m 2 .Thus, ␤ ϭ D͞d 2 will be in the range 28-450 for d in the range 1-4 m.In Fig. 2b, we used d ϭ 3 m.These values of ␤ are large enough so that the propagation is continuous.In this regime, the velocity does not depend strongly on d provided ⌫ ϭ 5 is fixed.Decreasing d to 1 m (while decreasing by a factor of 27, so that ͞d 3 remains constant) results in a solution that is indistinguishable from the one in Fig. 2b, which shows the wave profile at two different times.It travels nearly without deformation.
For steady propagation, we have derived a relationship (25) of the form where and are dimensionless velocities and ϵ vd͞D. [5] Combining Eqs. 2, 4, and 5, it follows that ϭ ͱ␤.
[6] It can be shown that F Ն 1 so that there is no solution to Eq. 3 if ⌫ Ͻ 1.This corresponds to the intuition that waves will not propagate if the release concentration is not sufficient to bring the neighboring sites up to threshold.In myocytes, the diffusion constant of Ca 2ϩ is of the order of D ϭ 30 m 2 s Ϫ1 , d Ϸ 2 m, Ϸ 15 ms, and [Ca 2ϩ ] T Ϸ 0.2 M (24, 28).Under conditions in which isolated sparks are found experimentally, the unitary Ca 2ϩ current at a release site is estimated to be 2pA (13,16), so that including the buffering, we find Ϸ 1.5 ϫ 10 Ϫ15 mol.Using Eqs. 1 and 2, this gives ⌫ Ϸ 1 and ␤ Ϸ 0.1, in agreement with the observation that waves do not propagate under these conditions.We have shown for the extreme saltatory limit ( 3 0) that the fire-diffuse-fire model undergoes a series of bifurcations and crises that lead to chaos as propagation failure is approached (i.e., as ⌫ is decreased) (25,28).The general feature of complex dynamics should be observed whenever such a system is driven towards propagation failure.On the other hand, waves do propagate in myocytes under conditions of Ca 2ϩ overload (20), which increases the unitary current, thus increasing ⌫.Simulations and more detailed calculations (27,28) suggest that ⌫ Ϸ 5.This generates waves with the experimentally observed velocity of about 65 m/s, and since ␤ Ϸ 0.1, these waves are saltatory, in agreement with experiment.
Although the explicit form of F in terms of and is complicated (25), the resulting function is simple.Fig. 3 is a contour plot of F(, ) upon which are superimposed straight lines ϭ ͌ ␤ for several values of ␤.Any two of the quantities ⌫, ␤, , and , when different from zero, determine the other two.
If 2 ͞ Ͻ Ͻ 1 (which holds for small enough ␤) the propagation is saltatory.In this case, F can be approximated by a function of only so that Eq. 3 reduces to F(, 0) ϵ g() ϭ ⌫ from which we find that ϭ g Ϫ1 (⌫) so that Here, g Ϫ1 denotes the inverse function of g.If ⌫ is large enough, we get an approximate explicit formula for g Ϫ1 which gives v Ϸ 4D͞dlog ⌫.If 2 ͞ Ͼ Ͼ 1 and ␤ Ͼ Ͼ 1, propagation is continuous.In this case, F can be approximated by a function of only so that Eq. 3 reduces to F(0, ) ϵ f() ϭ ⌫, from which we find that ϭ f Ϫ1 (⌫) so that Here, if ⌫ is large enough, we get an approximate explicit formula for f Ϫ1 that gives v Ϸ ͌ ⌫ ͌ D͞.This is Luther's equation (15,16) (v ϭ ␣ ͌ D͞) with ␣ ϭ ͌ ⌫.For the continuous-fertilization Ca 2ϩ wave in Fig. 1b, ⌫ ϭ 5, and this approximation is accurate to within a few percent.It is, however, completely incorrect for saltatory waves.Typically, propagation is saltatory if ␤ Ͻ Ͻ 1 and continuous if ␤ Ͼ Ͼ 1.Given a value of ⌫, it is clear from Fig. 3 that for small enough ␤ the intersection between the line ϭ ͌ ␤ and the contour will occur at values that satisfy 2 ͞ Ͻ Ͻ 1.This is the saltatory case.On the other hand, for large enough ␤, the intersection will occur at values that satisfy 2 ͞ Ͼ Ͼ 1.This is the continuous case.One of the most distinctive physical differences between continuous and saltatory waves is the number of sites that are simultaneously firing.This is given by 1 ϩ v͞d ϭ 1 ϩ 2 ͞.The curve ϭ 2 is shown as the dashed line in Fig. 3. Data points lying below this curve have only a single site firing at a time and those above have two or more.As illustrated by the vertical lines in Fig. 2b, many sites fire simultaneously during continuous propagation.On the other hand, as shown in Fig. 2a, only a single site fires at a time during saltatory propagation.More detailed models will share this and other general features of the fire-diffuse-fire model.Thus, whether a wave is saltatory or continuous will be determined by ␤.The trend of increasing velocity with increasing release concentration and the approach to propagation failure as the release density is decreased will both hold.The parameter values used here for oocytes, myocytes, and eggs are within experimentally accepted ranges.However, the specific values were chosen to reproduce the experimental velocities.Other values within the experimentally observed range give velocities that are off by as much as a factor of 5.The match (or mismatch) between experimentally observed and theoretically predicted velocities cannot be used either to validate or to invalidate a particular model.For example, using the Luther equation as in refs.15 and 16 (v ϭ ͌ D/) gives a velocity of 22 m/s for the waves in the immature oocyte as in Fig. 2a.These values (D ϭ 25 m 2 /s, ϭ 0.05s, and v ϭ 22͞s) are all within the experimental uncertainty.However, the Luther equation is manifestly wrong since the wave is saltatory (1,2).Decreasing by a factor of 50 in Fig. 2a results in a change in v of only a few percent, whereas the Luther equation yields a velocity of 160 m/s!The fact that ͌ D͞ is approximately equal to the velocity of the model saltatory wave is purely coincidental.
Whether the mode of propagation is continuous or saltatory is at least as important a fact as the numerical value of the velocity.Our analysis provides a simple, robust prescription for predicting the mode of propagation: estimate ␤.If estimates of ␤ are all Ͻ0.1, it is safe to assume that propagation is saltatory; if they are Ͼ10, propagation is continuous.For intermediate values of ␤ (0.1 Յ ␤ Յ 10), a more accurate model (namely, ref. 25) would be necessary to determine the mode of propagation.As it is widely held that information is encoded in the Ca 2ϩ signal (7)(8)(9)(10)(11)(12), the distinction between saltatory and continuous waves likely has considerable physiological significance.
The fire-diffuse-fire model of calcium propagation is described by the single-evolution equation where [Ca 2ϩ ](x, t) is the average concentration of calcium in the directions perpendicular to that of propagation (x), ␦() is the ␦ function, H() is the Heaviside step function [H() ϭ 0 for Ͻ 0, H() ϭ 1 for Ն 0], D is the diffusion coefficient of calcium, t i is the first time at which the ith site takes on the threshold value, [Ca 2ϩ ] T , is the total number of ions released by a site, and is the ''rise time'' for the receptor (i.e., the length of time the receptor is open during a release event).The ␦ function density for the receptors means we are treating the receptor͞channels as point sources.As stated in the text, the binding to buffers is assumed instantaneous.Thus, the effect of buffering is reflected in the values used for D and .Details can be found at http:͞͞www-xdiv.lanl.gov͞XCM͞pearson͞fdf.

FIG. 1 .
FIG. 1.(a) Space-time (or ''waterfall'') plot of the calcium wave in an immature Xenopus oocyte.The image was obtained with high-resolution confocal microscopy (1).''Warmer'' colors indicate elevated calcium levels.The individual release sites are clearly visible as streaks.(b) Relatively low-resolution snapshot of the calcium fertilization wave in mature Xenopus oocyte(5).The lighter region corresponds to higher calcium concentration, which is propagating into the region of lower calcium concentration, i.e., the wave is propagating to the left and up.The diameter is 1.2 mm.The front is ϳ50 m deep.If the mean site spacing is 2 m, then the zone in which calcium is actively being released (i.e., the front) is ϳ25 sites deep.The resolution is too low to see individual sites.One of the predictions of the fire-diffuse-fire model is that even with the spatial resolution of a, it would not be possible to see the individual sites because the release rate is so slow that the calcium diffuses away before significant local buildup can occur.The tiny local buildups that do occur are statistical in nature.

FIG. 2 .
FIG. 2. (a) Four snapshots of a rightward saltatory wave solution of the fire-diffuse-fire model.The dashed line denotes the threshold concentration, [Ca 2ϩ ]T.The values used here were taken from the literature on the immature Xenopus oocyte (1-4): d ϭ 3.9 m, ϭ 0.05 s, D ϭ 25 m 2 ͞s [Ca 2ϩ ]T ϭ 0.25 M, [Ca 2ϩ ]b ϭ 0.05 M, and ϭ 3.5 ϫ 10 Ϫ14 mol.The wave is traveling to the right with a velocity of 20 m͞s.Note that only one site fires at a time and that the shape of the wave changes with time.(b) Two snapshots of a rightward continuous wave.The vertical lines indicate the sites that are simultaneously firing at the time of the first snapshot (solid line).The dashed line corresponds to a later time.The parameters used here were taken from the literature on the fertilization wave in Xenopus eggs (5).The values we used correspond to the interior region of the egg: d ϭ 3 m, D ϭ 50 m 2 ͞s, ϭ 9 s, ϭ 1.35 ϫ 10 Ϫ14 mol, [Ca 2ϩ ]T ϭ 0.4 M, and [Ca 2ϩ ]b ϭ 0.3 M.The wave is traveling to the right with a velocity of 5.2 m͞s.Note that many sites are firing simultaneously and that the wave travels without observable deformation.Despite the larger diffusion coefficient, the continuum wave is slower than the saltatory wave.This is because of the difference in release rates ().