Calcium waves in agarose gel with cell organelles: implications of the velocity curvature relationship

Calcium oscillations and waves have been observed not only in several types of living cells but also in less complex systems of isolated cell organelles. Here we report the determination of apparent Ca2+ diffusion coefficients in a novel excitable medium of agarose gel with homogeneously distributed vesicles of skeletal sarcoplasmic reticulum. Spatiotemporal calcium patterns were visualized by confocal laser scanning fluorescence microscopy. To obtain characteristic parameters of the velocity curvature relationship, namely, apparent diffusion coefficient, velocity of plane calcium waves, and critical radius, positively and negatively curved wave fronts were analyzed. It is demonstrated that gel-immobilized cell organelles reveal features of an excitable medium. Apparent Ca2+ diffusion coefficients of the in vitro system, both in the absence or in the presence of mitochondria, were found to be higher than in cardiac myocytes and lower than in unbuffered agarose gel. Plane calcium waves propagated markedly slower in the in vitro system than in rat cardiac myocytes. Whereas mitochondria significantly reduced the apparent Ca2+ diffusion coefficient of the in vitro system, propagation velocity and critical size of calcium waves were found to be nearly unchanged. These results suggest that calcium wave propagation depends on the kinetics of calcium release rather than on diffusion.     

Two-dimensional model of calcium waves reproduces the patterns observed in Xenopus oocytes.

Biological excitability enables the rapid transmission of physiological signals over distance. Using confocal fluorescence microscopy, we previously reported circular, planar, and spiral waves of Ca2+ in Xenopus laevis oocytes that annihilated one another upon collision. We present experimental evidence that the excitable process underlying wave propagation depends on Ca2+ diffusion and does not require oscillations in inositol (1,4,5)trisphosphate (IP3) concentration. Extending an existing ordinary differential equation (ODE) model of Ca2+ oscillations to two spatial dimensions, we develop a partial differential equation (PDE) model of Ca2+ excitability. The model assumes that cytosolic Ca2+ couples neighboring Ca2+ release sites. This simple PDE model qualitatively reproduces our experimental observations.     

Migraine aura dynamics after reverse retinotopic mapping of weak excitation waves in the primary visual cortex

 A kinematical model for excitable wave propagation is analyzed to describe the dynamics of a typical neurological symptom of migraine. The kinematical model equation is solved analytically for a linear dependency between front curvature and velocity. The resulting wave starts from an initial excitation and moves in the medium that represents the primary visual cortex. Due to very weak excitability the wave propagates only across a confined area and eventually disappears. This cortical excitation pattern is projected onto a visual hemifield by reverse retinotopic mapping. Weak excitability explains the confined appearance of aura symptoms in time and sensory space. The affected area in the visual field matches in growth and form the one reported by migraine sufferers. The results can be extended from visual to tactile and to other sensory symptoms. If the spatiotemporal pattern from our model can be matched in future investigations with those from introspectives, it would allow one to draw conclusions on topographic mapping of sensory input in human cortex. Received: 25 April 2002 / Accepted: 20 February 2003 / Published online: 20 May 2003 RID="*" ID="*" Present address: M. A. Dahlem Leibniz-Institut für Neurobiologie, Brenneckestr. 6, 39118 Magdeburg, Germany Acknowledgements. We would like to thank V. Zykov for useful discussions on wave Propagation, and one of us (MAD) would like to thank Ed Chronicle for useful discussions on functional excitability. This project was supported by a scholarship Landesstipendium Sachsen-Anhalt to MAD. Correspondence to: M. A. Dahlem (e-mail: dahlem@ifn-magdeburg.de)     

The effect of heterogeneously-distributed RyR channels on calcium dynamics in cardiac myocytes

Calcium plays an essential role in excitation-contraction coupling in muscle, and derangements in calcium handling can produce a variety of potentially harmful conditions, especially in cardiac muscle. In cardiac tissue specialized invaginations of the sarcolemma, called T-tubules, penetrate deep into each sarcomere, and depolarization of the SL leads to an influx of calcium through voltage-sensitive channels in the T-tubules that in turn triggers further calcium release from the sarcoplasmic reticulum via ryanodine-sensitive calcium channels. Under certain conditions, such as elevated external Ca²⁺, cardiac cells can release calcium from the sarcoplasmic reticulum spontaneously, producing a calcium ’spark’ and propagating traveling waves of elevated Ca²⁺ concentration, without depolarization of the SL (Wier and Blatter, 1991a, Cell Calcium 12, 241–254; Williams, 1993, Cell Calcium 14, 724–735; Cheng et al., 1993a, Science 262, 740–744). However, under normal resting conditions these potentially harmful waves seldom occur. In this paper we investigate the role of the periodic distribution of ryanodine-sensitive channels in determining whether a spark can trigger a wave, using a modification of the kinetic model proposed by Tang and Othmer, 1994b, Biophys. J. 67, 2223–2235, for calcium-induced calcium release. We show that the spatial localization of these channels near the T-tubules has a significant effect on both wave propagation and the onset of oscillations in this system. Spatial localization provides a possible explanation for the differing effects of various experimental protocols on the system’s ability to propagate a traveling wave. Supported in part by NIH Grant GM29123.     

Do calcium buffers always slow down the propagation of calcium waves?

Calcium buffers are large proteins that act as binding sites for free cytosolic calcium. Since a large fraction of cytosolic calcium is bound to calcium buffers, calcium waves are widely observed under the condition that free cytosolic calcium is heavily buffered. In addition, all physiological buffered excitable systems contain multiple buffers with different affinities. It is thus important to understand the properties of waves in excitable systems with the inclusion of buffers. There is an ongoing controversy about whether or not the addition of calcium buffers into the system always slows down the propagation of calcium waves. To solve this controversy, we incorporate the buffering effect into the generic excitable system, the FitzHugh–Nagumo model, to get the buffered FitzHugh–Nagumo model, and then to study the effect of the added buffer with large diffusivity on traveling waves of such a model in one spatial dimension. We can find a critical dissociation constant ( $K=K(a)$ ) characterized by system excitability parameter $a$ such that calcium buffers can be classified into two types: weak buffers ( $K\in (K(a),\infty )$ ) and strong buffers ( $K\in (0,K(a))$ ). We analytically show that the addition of weak buffers or strong buffers but with its total concentration $b_0^1$ below some critical total concentration $b_{0,c}^1$ into the system can generate a traveling wave of the resulting system which propagates faster than that of the origin system, provided that the diffusivity $D_1$ of the added buffers is sufficiently large. Further, the magnitude of the wave speed of traveling waves of the resulting system is proportional to $\sqrt{D_1}$ as $D_1\rightarrow \infty $ . In contrast, the addition of strong buffers with the total concentration $b_0^1>b_{0,c}^1$ into the system may not be able to support the formation of a biologically acceptable wave provided that the diffusivity $D_1$ of the added buffers is sufficiently large.     

Molecular mechanisms of intracellular calcium excitability in X. laevis oocytes.

Following receptor activation in Xenopus oocytes, spiral waves of intracellular Ca2+ release were observed. We have identified key molecular elements in the pathway that give rise to Ca2+ excitability. The patterns of Ca2+ release produced by GTP-gamma-S and by inositol 1,4,5-trisphosphate (IP3) are indistinguishable from receptor-induced Ca2+ patterns. The regenerative Ca2+ activity is critically dependent on the presence of IP3 and on the concentration of intracellular Ca2+, but is independent of extracellular Ca2+. Broad regions of the intracellular milieu can be synchronously excited to initiate Ca2+ waves and produce pulsating foci of Ca2+ release. By testing the temperature dependence of wavefront propagation, we provide evidence for an underlying process limited by diffusion, consistent with the elementary theory of excitable media. We propose a model for intracellular Ca2+ signaling in which wave propagation is controlled by IP3-mediated Ca2+ release from internal stores, but is modulated by the cytoplasmic concentration and diffusion of Ca2+.     

Dynamics of excitable elements with time-delayed coupling

Motivated by recent experiments on intracellular calcium release we study the effects of different types of coupling on the dynamics of arrays of excitable elements. We intend to find a mechanism that produces a sustained activity of the elements following a spike. While instantaneous diffusive coupling does not exhibit this property, we show that, for a coupling term with temporal delay, signals from adjacent elements can serve as mutual excitations and thus prolong the duration of the signal. We propose that time delayed coupling is generated by diffusion between isolated clusters of calcium channels. Our model could thus provide an explanation for two different release modes observed in the Ca²⁺ system.     

Local positive feedback by calcium in the propagation of intracellular calcium waves.

In many types of eukaryotic cells, the activation of surface receptors leads to the production of inositol 1,4,5-trisphosphate and calcium release from intracellular stores. Calcium release can occur in complex spatial patterns, including waves of release that traverse the cytoplasm. Fluorescence video microscopy was used to view calcium waves in single mouse neuroblastoma cells. The propagation of calcium waves was slowed by buffers that bind calcium quickly, such as BAPTA, but not by a buffer with slower on-rate, EGTA. This shows that a key feedback event in wave propagation is rapid diffusion of calcium occurring locally on a scale of < 1 micron. The length-speed product of wavefronts was used to determine that calcium acting in feedback diffuses at nearly the rate expected for free diffusion in aqueous solution. In cytoplasm, which contains immobile Ca2+ buffers, this rate of diffusion occurs only in the first 0.2 ms after release, within 0.4 micron of a Ca2+ release channel mouth. Calcium diffusion from an open channel to neighboring release sites is, therefore, a rate-determining regenerative step in calcium wave propagation. The theoretical limitations of the wave front analysis are discussed.     

Development of electrical excitability in embryonic neurons: Mechanisms and roles

Xenopus spinal neurons serve as a nearly ideal population of excitable cells for study of developmental regulation of electrical excitability. On the one hand, the firing properties of these neurons can be directly examined at early stages of differentiation and membrane excitability changes as neurons mature. Underlying changes in voltage-dependent ion channels have been characterized and the mechanisms that bring about these changes are being defined. On the other hand, these neurons have been shown to be spontaneously active at stages when action potentials provide significant calcium entry. Calcium entry provokes further elevation of intracellular calcium via release from intracellular stores. The resultant transient elevations of intracellular calcium encode differentiation in their frequency. Recent studies have shown that different neuronal subpopulations enlist distinct mechanisms for regulation of excitability and recruit specific programs of differentiation by particular patterns of activity. © 1998 John Wiley & Sons, Inc. J Neurobiol 37: 190–197, 1998     

The low conductance mitochondrial permeability transition pore confers excitability and CICR wave propagation in a computational model

Mitochondria have long been known to sequester cytosolic Ca²⁺ and even to shape intracellular patterns of endoplasmic reticulum-based Ca²⁺ signaling. Evidence suggests that the mitochondrial network is an excitable medium which can demonstrate independent Ca²⁺ induced Ca²⁺ release via the mitochondrial permeability transition. The role of this excitability remains unclear, but mitochondrial Ca²⁺ handling appears to be a crucial element in diverse diseases as diabetes, neurodegeneration and cardiac dysfunction that also have bioenergetic components. In this paper, we extend the modular Magnus-Keizer computational model for respiration-driven Ca²⁺ handling to include a permeability transition based on a channel-like pore mechanism. We demonstrate both excitability and Ca²⁺ wave propagation accompanied by depolarizations qualitatively similar to those reported in cell and isolated mitochondria preparations. These waves depend on the energy state of the mitochondria, as well as other elements of mitochondrial physiology. Our results support the concept that mitochondria can transmit state dependent signals about their function across the mitochondrial network. Our model provides the tools for predictions about the internal physiology that leads to this qualitatively different Ca²⁺ excitability seen in mitochondria.     

Saltatory propagation of Ca2+ waves by Ca2+ sparks.

Punctate releases of Ca2+, called Ca2+ sparks, originate at the regular array of t-tubules in cardiac myocytes and skeletal muscle. During Ca2+ overload sparks serve as sites for the initiation and propagation of Ca2+ waves in myocytes. Computer simulations of spark-mediated waves are performed with model release sites that reproduce the adaptive Ca2+ release observed for the ryanodine receptor. The speed of these waves is proportional to the diffusion constant of Ca2+, D, rather than D, as is true for reaction-diffusion equations in a continuous excitable medium. A simplified "fire-diffuse-fire" model that mimics the properties of Ca2+-induced Ca2+ release (CICR) from isolated sites is used to explain this saltatory mode of wave propagation. Saltatory and continuous wave propagation can be differentiated by the temperature and Ca2+ buffer dependence of wave speed.     

Analytical calculation of intracellular calcium wave characteristics

We present a theoretical analysis of intracellular calcium waves propagated by calcium feedback at the inositol 1,4,5-trisphosphate (IP3) receptor. The model includes essential features of calcium excitability, but is still analytically tractable. Formulas are derived for the wave speed, amplitude, and width. The calculations take into account cytoplasmic Ca buffering, the punctate nature of the Ca release channels, channel inactivation, and Ca pumping. For relatively fast buffers, the wave speed is well approximated by V(infinity) = (J(eff)D(eff)/C0)1/2, where J(eff) is an effective, buffered source strength; D(eff) is the effective, buffered diffusion constant of Ca; and C(0) is the Ca threshold for channel activation. It is found that the saturability and finite on-rate of buffers must be taken into account to accurately derive the wave speed and front width. The time scale governing Ca wave propagation is T(r), the time for Ca release to reach threshold to activate further release. Because IP3 receptor inactivation is slow on this time scale, channel inactivation does not affect the wave speed. However, inactivation competes with Ca removal to limit wave height and front length, and for biological parameter ranges, it is inactivation that determines these parameters. Channel discreteness introduces only small corrections to wave speed relative to a model in which Ca is released uniformly from the surface of the stores. These calculations successfully predict experimental results from basic channel and cell parameters and explain the slowing of waves by exogenous buffers.     

Traveling waves in the discrete fast buffered bistable system

We study the existence and uniqueness of traveling wave solutions of the discrete buffered bistable equation. Buffered excitable systems are used to model, among other things, the propagation of waves of increased calcium concentration, and discrete models are often used to describe the propagation of such waves across multiple cells. We derive necessary conditions for the existence of waves, and, under some restrictive technical assumptions, we derive sufficient conditions. When the wave exists it is unique and stable.      

On Propagation of Excitation Waves in Moving Media: The FitzHugh-Nagumo Model

Background

Existence of flows and convection is an essential and integral feature of many excitable media with wave propagation modes, such as blood coagulation or bioreactors.

Methods/Results

Here, propagation of two-dimensional waves is studied in parabolic channel flow of excitable medium of the FitzHugh-Nagumo type. Even if the stream velocity is hundreds of times higher that the wave velocity in motionless medium (), steady propagation of an excitation wave is eventually established. At high stream velocities, the wave does not span the channel from wall to wall, forming isolated excited regions, which we called “restrictons”. They are especially easy to observe when the model parameters are close to critical ones, at which waves disappear in still medium. In the subcritical region of parameters, a sufficiently fast stream can result in the survival of excitation moving, as a rule, in the form of “restrictons”. For downstream excitation waves, the axial portion of the channel is the most important one in determining their behavior. For upstream waves, the most important region of the channel is the near-wall boundary layers. The roles of transversal diffusion, and of approximate similarity with respect to stream velocity are discussed.

Conclusions

These findings clarify mechanisms of wave propagation and survival in flow.     

Modeling action potential generation and propagation in NRK fibroblasts

Normal rat kidney (NRK) fibroblasts change their excitability properties through the various stages of cell proliferation. The present mathematical model has been developed to explain excitability of quiescent (serum deprived) NRK cells. It includes as cell membrane components, on the basis of patch-clamp experiments, an inwardly rectifying potassium conductance (G_Kir), an L-type calcium conductance (G_CaL), a leak conductance (G_leak), an intracellular calcium-activated chloride conductance [G_Cl(Ca)], and a gap junctional conductance (G_gj), coupling neighboring cells in a hexagonal pattern. This membrane model has been extended with simple intracellular calcium dynamics resulting from calcium entry via G_CaL channels, intracellular buffering, and calcium extrusion. It reproduces excitability of single NRK cells and cell clusters and intercellular action potential (AP) propagation in NRK cell monolayers. Excitation can be evoked by electrical stimulation, external potassium-induced depolarization, or hormone-induced intracellular calcium release. Analysis shows the roles of the various ion channels in the ultralong (30 s) NRK cell AP and reveals the particular role of intracellular calcium dynamics in this AP. We support our earlier conclusion (De Roos A, Willems PH, van Zoelen EJ, and Theuvenet AP. Am J Physiol Cell Physiol 273: C1900–C1907, 1997) that AP generation and propagation may act as a rapid mechanism for the propagation of intracellular calcium waves, thus contributing to fast intercellular calcium signaling. The present model serves as a starting point to further analyze excitability changes during contact inhibition and cell transformation. Hodgkin-Huxley model; intracellular calcium dynamics; L-type calcium conductance; inward rectifier; calcium-activated chloride conductance; gap junctional coupling     

Image processing techniques applied to excitation waves in the chicken retina

Image processing techniques are described in detail that are used to gain information about the dynamics of wave propagation in excitable media. We focus on a phenomenon called spreading depression (SD) observed in the chicken retina, but the techniques described here concern a large variety of excitable systems. Despite the impressive progress both in SD research of the past 50 years and, during nearly the same period, in the theory of self-organization of wave patterns, there is still little mutual overlap. However, the increasing demands for understanding complex systems, like neuronal tissue, require such theoretical concepts. Arguments are given why the chicken retina is a nearly perfect experimental system for assessing and further developing these concepts.     

A mode control model of a neuron's axon and dendrites

The ability of a neuron network to process information depends upon the ability of the individual neurons to transport impulses and to control the signal transport process in other neurons. The transport process for the action potential seen at the axon depends upon the excitable characteristic of the neural membrane. Propagation of signals in the dendrites, where synaptic imputs are most likely processed, is not clearly understood. Extracellular recordings of dendritic systems indicate that the dendrites are partially excitable and can conduct spikes. Further, electrical stimulation of the reticular formation or specific thalamic nuclei suggest that the conduction process can be modified in the dendrites of cortical cells.A Mode Control model is described which demonstrates many of the observed transport and control properties of dendrite and axon membrane. The model is based upon a simple extension of Fitzhugh's BVP model. Lateral transport over the membrane has been introduced by applying Kirchhoff's laws. Reinterpreting the variables, the influence of membrane potential, pH, and calcium ions can be identified. Modification of the voltage-current characteristic of the membrane model can change the axon model to a dendrite model. The dendrite model possesses a diffusion equation mode, a wave equation mode and a pulse mode. Signals are transferred in the wave and pulse mode and blocked in the diffusion mode. The dendrite's mode is controlled by the resting depolarization level. Experimental evidence tends to confirm these phenomena.The work described in this paper was performed while attending the University of California, Berkeley, under a National Institutes of Health Traineeship.     

The effect of ion pumps on the speed of travelling waves in the fire-diffuse-fire model of Ca2+ release

The fire-diffuse-fire model provides an idealized model of Ca²⁺ release within living cells. The effect of calcium pumps, which drive Ca²⁺ back into internal stores, is often neglected for mathematical simplicity. Here we show how to explicitly analyse such effects by extending the work of Keizer et al. [Keizer, J. E., G. D. Smith, S. Ponce Dawson and J. Pearson (1998). Saltatory propagation of Ca²⁺ waves by Ca²⁺ sparks. Biophys. J. 75, 595–600.]. For travelling waves, in which release events occur sequentially, we construct the speed of waves in terms of the time-scale at which pumps operate. An immediate consequence of this analysis is that the inclusion of calcium pumps leads to multiple solutions. A linear stability analysis determines those solution branches in parameter space which are stable. Numerical continuation is used to provide explicit examples of the bifurcation diagrams of the speed of waves as a function of physiologically significant system parameters.     

Intercellular Ca2+ wave propagation through gap-junctional Ca2+ diffusion: a theoretical study

Intercellular regenerative calcium waves in systems such as the liver and the blowfly salivary gland have been hypothesized to spread through calcium-induced calcium release (CICR) and gap-junctional calcium diffusion. A simple mathematical model of this mechanism is developed. It includes CICR and calcium removal from the cytoplasm, cytoplasmic and gap-junctional calcium diffusion, and calcium buffering. For a piecewise linear approximation of the calcium kinetics, expressions in terms of the cellular parameters are derived for 1) the condition for the propagation of intercellular waves, and 2) the characteristic time of the delay of a wave encountered at the gap junctions. Intercellular propagation relies on the local excitation of CICR in the perijunctional space by gap-junctional calcium influx. This mechanism is compatible with low effective calcium diffusivity, and necessitates that CICR can be excited in every cell along the path of a wave. The gap-junctional calcium permeability required for intercellular waves in the model falls in the range of reported gap-junctional permeability values. The concentration of diffusive cytoplasmic calcium buffers and the maximal rate of CICR, in the case of inositol 1,4,5-trisphosphate (IP3) receptor calcium release channels set by the IP(3) concentration, are shown to be further determinants of wave behavior.