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.     

Calcium wave propagation by calcium-induced calcium release: An unusual excitable system

We discuss in detail the behaviour of a model, proposed by Goldbeteret al. (1990.Proc. natn. Acad. Sci. 87, 1461–1465), for intracellular calcium wave propagation by calcium-induced calcium release, focusing our attention on excitability and the propagation of waves in one spatial dimension. The model with no diffusion behaves like a generic excitable system, and threshold behaviour, excitability and oscillations can be understood within this general framework. However, when diffusion is included, the model no longer behaves like a generic excitable system; the fast and slow variables are not distinct and previous results on excitable systems do not necessarily apply. We consider a piecewise linear simplification of the model, and construct travelling pulse and periodic plane wave solutions to the simplified model. The analogous behaviour in the full model is studied numerically. Goldbeter's model for calciuminduced calcium release is an excitable system of a type not previously studied in detail.     

Modelling the spatio-temporal dynamics of multi-species host-parasitoid interactions: heterogeneous patterns and ecological implications

A mathematical model of the spatio-temporal dynamics of a two host, two parasitoid system is presented. There is a coupling of the four species through parasitism of both hosts by one of the parasitoids. The model comprises a system of four reaction-diffusion equations. The underlying system of ordinary differential equations, modelling the host-parasitoid population dynamics, has a unique positive steady state and is shown to be capable of undergoing Hopf bifurcations, leading to limit cycle kinetics which give rise to oscillatory temporal dynamics. The stability of the positive steady state has a fundamental impact on the spatio-temporal dynamics: stable travelling waves of parasitoid invasion exhibit increasingly irregular periodic travelling wave behaviour when key parameter values are increased beyond their Hopf bifurcation point. These irregular periodic travelling waves give rise to heterogeneous spatio-temporal patterns of host and parasitoid abundance. The generation of heterogeneous patterns has ecological implications and the concepts of temporary host refuge and niche formation are considered.     

On the roles of Ca2+ diffusion, Ca2+ buffers, and the endoplasmic reticulum in IP3-induced Ca2+ waves.

We have investigated the effects of Ca2+ diffusion, mobile and stationary Ca2+ buffers in the cytosol, and Ca2+ handling by the endoplasmic reticulum on inositol 1,4,5-trisphosphate-induced Ca2+ wave propagation. Rapid equilibration of free and bound Ca2+ is used to describe Ca2+ sequestration by buffers in both the cytosol and endoplasmic reticulum (ER) lumen. Cytosolic Ca2+ regulation is based on a kinetic model of the inositol 1,4,5-trisphosphate (IP3) receptor of De Young and Keizer that includes activation and inhibition of the IP3 receptor Ca2+ channel in the ER membrane and SERCA Ca2+ pumps in the ER. Diffusion of Ca2+ in the cytosol and the ER and the breakdown and diffusion of IP3 are also included in our calculations. Although Ca2+ diffusion is severely limited because of buffering, when conditions are chosen just below the threshold for Ca2+ oscillations, a pulse of IP3 or Ca2+ results in a solitary trigger wave that requires diffusion of Ca2+ for its propagation. In the oscillatory regime repetitive wave trains are observed, but for this type of wave neither the wave shape nor the speed is strongly dependent on the diffusion of Ca2+. Local phase differences lead to waves that are predominately kinematic in nature, so that the wave speed (c) is related to the wavelength (lambda) and the period of the oscillations (tau) approximately by the formula c = lambda/tau. The period is determined by features that control the oscillations, including [IP3] and pump activity, which are related to recent experiments. Both solitary waves and wave trains are accompanied by a Ca2+ depletion wave in the ER lumen, similar to that observed in cortical preparations from sea urchin eggs. We explore the effect of endogenous and exogenous Ca2+ buffers on wave speed and wave shape, which can be explained in terms of three distinct effects of buffering, and show that exogenous buffers or Ca2+ dyes can have considerable influence on the amplitude and width of the waves.     

Stimulus-dependent control of inositol 1,4,5-trisphosphate-induced Ca(2+) oscillation frequency by the endoplasmic reticulum Ca(2+)-ATPase

In many cell types, receptor stimulation evokes cytosolic calcium oscillations with a frequency that depends on agonist dose. Previous studies demonstrated controversial effects of changing the activity of the endoplasmic reticulum Ca(2+)-ATPase upon the frequency of oscillations. By numerical simulations, we found that the model of De Young and Keizer (J. Keizer and G.W. De Young, 1994, J. Theor. Biol. 166: 431-442), unlike other models, can explain the observed discrepancies, assuming that the different experiments were performed at different stimulus levels. According to model predictions, partial inhibition of internal calcium pumps is expected to increase frequency at low stimulus strength and should have an opposite effect at strong stimuli. Similar results were obtained using an analytical estimation of oscillation period, based on calcium-dependent channel activation and inactivation. In experiments on HeLa cells, 4 nM thapsigargin increased the frequency of calcium oscillations induced by 1 and 2.5 microM histamine but had no effect on supramaximally stimulated cells. In HEp-2 cells, 2 nM thapsigargin slowed down the rapid, ATP-induced oscillations. Our results suggest that in the investigated cell types, the De Young-Keizer model based on inositol 1,4,5-trisphosphate-dependent calcium-induced calcium release can properly describe intracellular calcium oscillations.     

Receptors, sparks and waves in a fire-diffuse-fire framework for calcium release

Calcium ions are an important second messenger in living cells. Indeed calcium signals in the form of waves have been the subject of much recent experimental interest. It is now well established that these waves are composed of elementary stochastic release events (calcium puffs or sparks) from spatially localised calcium stores. The aim of this paper is to analyse how the stochastic nature of individual receptors within these stores combines to create stochastic behaviour on long time-scales that may ultimately lead to waves of activity in a spatially extended cell model. Techniques from asymptotic analysis and stochastic phase–plane analysis are used to show that a large cluster of receptor channels leads to a release probability with a sigmoidal dependence on calcium density. This release probability is incorporated into a computationally inexpensive model of calcium release based upon a stochastic generalisation of the fire-diffuse-fire (FDF) threshold model. Numerical simulations of the model in one and two dimensions (with stores arranged on both regular and disordered lattices) illustrate that stochastic calcium release leads to the spontaneous production of calcium sparks that may merge to form saltatory waves. Illustrations of spreading circular waves, spirals and more irregular waves are presented. Furthermore, receptor noise is shown to generate a form of array enhanced coherence resonance whereby all calcium stores release periodically and simultaneously.     

Global bifurcations in the FitzHugh-Nagumo equations for nerve wave propagation

The FitzHugh-Nagumo equations for action potential propagation along nerve axons and the corresponding ordinary differential equations for travelling waves are solved numerically. Above a critical value, a constant bias current can drive a wave-front solution. At the critical value, a global bifurcation occurs. As a result, the wave front switches into a pulse.Based on a thesis by one of the authors (H. F.).     

Travelling Waves in Hybrid Chemotaxis Models

Hybrid models of chemotaxis combine agent-based models of cells with partial differential equation models of extracellular chemical signals. In this paper, travelling wave properties of hybrid models of bacterial chemotaxis are investigated. Bacteria are modelled using an agent-based (individual-based) approach with internal dynamics describing signal transduction. In addition to the chemotactic behaviour of the bacteria, the individual-based model also includes cell proliferation and death. Cells consume the extracellular nutrient field (chemoattractant), which is modelled using a partial differential equation. Mesoscopic and macroscopic equations representing the behaviour of the hybrid model are derived and the existence of travelling wave solutions for these models is established. It is shown that cell proliferation is necessary for the existence of non-transient (stationary) travelling waves in hybrid models. Additionally, a numerical comparison between the wave speeds of the continuum models and the hybrid models shows good agreement in the case of weak chemotaxis and qualitative agreement for the strong chemotaxis case. In the case of slow cell adaptation, we detect oscillating behaviour of the wave, which cannot be explained by mean-field approximations.     

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.     

From periodic travelling waves to travelling fronts in the spike-diffuse-spike model of dendritic waves

In the vertebrate brain excitatory synaptic contacts typically occur on the tiny evaginations of neuron dendritic surface known as dendritic spines. There is clear evidence that spine heads are endowed with voltage-dependent excitable channels and that action potentials invade spines. Computational models are being increasingly used to gain insight into the functional significance of a spine with an excitable membrane. The spike-diffuse-spike (SDS) model is one such model that admits to a relatively straightforward mathematical analysis. In this paper we demonstrate that not only can the SDS model support solitary travelling pulses, already observed numerically in more detailed biophysical models, but that it has periodic travelling wave solutions. The exact mathematical treatment of periodic travelling waves in the SDS model is used, within a kinematic framework, to predict the existence of connections between two periodic spike trains of different interspike interval. The associated wave front in the sequence of interspike intervals travels with a constant velocity without degradation of shape, and might therefore be used for the robust encoding of information.     

Early development and quorum sensing in bacterial biofilms

 We develop mathematical models to examine the formation, growth and quorum sensing activity of bacterial biofilms. The growth aspects of the model are based on the assumption of a continuum of bacterial cells whose growth generates movement, within the developing biofilm, described by a velocity field. A model proposed in Ward et al. (2001) to describe quorum sensing, a process by which bacteria monitor their own population density by the use of quorum sensing molecules (QSMs), is coupled with the growth model. The resulting system of nonlinear partial differential equations is solved numerically, revealing results which are qualitatively consistent with experimental ones. Analytical solutions derived by assuming uniform initial conditions demonstrate that, for large time, a biofilm grows algebraically with time; criteria for linear growth of the biofilm biomass, consistent with experimental data, are established. The analysis reveals, for a biologically realistic limit, the existence of a bifurcation between non-active and active quorum sensing in the biofilm. The model also predicts that travelling waves of quorum sensing behaviour can occur within a certain time frame; while the travelling wave analysis reveals a range of possible travelling wave speeds, numerical solutions suggest that the minimum wave speed, determined by linearisation, is realised for a wide class of initial conditions. Received: 10 February 2002 / Revised version: 29 October 2002 / Published online: 19 March 2003 Key words or phrases: Bacterial biofilm – Quorum sensing – Mathematical modelling – Numerical solution – Asymptotic analysis – Travelling wave analysis     

Properties of intracellular Ca2+ waves generated by a model based on Ca(2+)-induced Ca2+ release.

Cytosolic Ca2+ waves occur in a number of cell types either spontaneously or after stimulation by hormones, neurotransmitters, or treatments promoting Ca2+ influx into the cells. These waves can be broadly classified into two types. Waves of type 1, observed in cardiac myocytes or Xenopus oocytes, correspond to the propagation of sharp bands of Ca2+ throughout the cell at a rate that is high enough to permit the simultaneous propagation of several fronts in a given cells. Waves of type 2, observed in hepatocytes, endothelial cells, or various kinds of eggs, correspond to the progressive elevation of cytosolic Ca2+ throughout the cell, followed by its quasi-homogeneous return down to basal levels. Here we analyze the propagation of these different types of intracellular Ca2+ waves in a model based on Ca(2+)-induced Ca2+ release (CICR). The model accounts for transient or sustained waves of type 1 or 2, depending on the size of the cell and on the values of the kinetic parameters that measure Ca2+ exchange between the cytosol, the extracellular medium, and intracellular stores. Two versions of the model based on CICR are considered. The first version involves two distinct Ca2+ pools sensitive to inositol 1,4,5-trisphosphate (IP3) and Ca2+, respectively, whereas the second version involves a single pool sensitive both to Ca2+ and IP3 behaving as co-agonists for Ca2+ release. Intracellular Ca2+ waves occur in the two versions of the model based on CICR, but fail to propagate in the one-pool model at subthreshold levels of IP3. For waves of type 1, we investigate the effect of the spatial distribution of Ca(2+)-sensitive Ca2+ stores within the cytosol, and show that the wave fails to propagate when the distance between the stores exceeds a critical value on the order of a few microns. We also determine how the period and velocity of the waves are affected by changes in parameters measuring stimulation, Ca2+ influx into the cell, or Ca2+ pumping into the stores. For waves of type 2, the numerical analysis indicates that the best qualitative agreement with experimental observations is obtained for phase waves. Finally, conditions are obtained for the occurrence of "echo" waves that are sometimes observed in the experiments.     

Simulation of waves in calcium models with 3D spherical geometry

Waves of calcium ions are present in fertilized eggs of many species. Models for pulse and tidal wave propagation have usually been studied in one or two spatial coordinates only. We examine in three spatial coordinates some established models, based on Ca(2+)-induced Ca(2+)-release from both (assumed) continuously or heterogeneously distributed stores of endoplasmic reticulum (ER) through channels activated by inositol triphosphate (IP(3)). With continuous IP(3) distribution decreasing radially towards the interior, we obtain concave pulse shapes for waves penetrating the interior. Concave waves are also recorded in systems with ER confined to distributions of small spheres (microdomains) inside the cell, which we simulate for front waves (tides) in bistable systems.     

Coexistence of Tonic Spiking Oscillations in a Leech Neuron Model

The leech neuron model studied here has a remarkable dynamical plasticity. It exhibits a wide range of activities including various types of tonic spiking and bursting. In this study we apply methods of the qualitative theory of dynamical systems and the bifurcation theory to analyze the dynamics of the leech neuron model with emphasis on tonic spiking regimes. We show that the model can demonstrate bi-stability, such that two modes of tonic spiking coexist. Under a certain parameter regime, both tonic spiking modes are represented by the periodic attractors. As a bifurcation parameter is varied, one of the attractors becomes chaotic through a cascade of period-doubling bifurcations, while the other remains periodic. Thus, the system can demonstrate co-existence of a periodic tonic spiking with either periodic or chaotic tonic spiking. Pontryagins averaging technique is used to locate the periodic orbits in the phase space.     

ATP- and Gap Junction–dependent Intercellular Calcium Signaling in Osteoblastic Cells

Many cells coordinate their activities by transmitting rises in intracellular calcium from cell to cell. In nonexcitable cells, there are currently two models for intercellular calcium wave propagation, both of which involve release of inositol trisphosphate (IP₃)- sensitive intracellular calcium stores. In one model, IP₃ traverses gap junctions and initiates the release of intracellular calcium stores in neighboring cells. Alternatively, calcium waves may be mediated not by gap junctional communication, but rather by autocrine activity of secreted ATP on P₂ purinergic receptors. We studied mechanically induced calcium waves in two rat osteosarcoma cell lines that differ in the gap junction proteins they express, in their ability to pass microinjected dye from cell to cell, and in their expression of P2Y₂ (P_2U) purinergic receptors. ROS 17/2.8 cells, which express the gap junction protein connexin43 (Cx43), are well dye coupled, and lack P_2U receptors, transmitted slow gap junction-dependent calcium waves that did not require release of intracellular calcium stores. UMR 106-01 cells predominantly express the gap junction protein connexin 45 (Cx45), are poorly dye coupled, and express P_2U receptors; they propagated fast calcium waves that required release of intracellular calcium stores and activation of P_2U purinergic receptors, but not gap junctional communication. ROS/P_2U transfectants and UMR/Cx43 transfectants expressed both types of calcium waves. Gap junction–independent, ATP-dependent intercellular calcium waves were also seen in hamster tracheal epithelia cells. These studies demonstrate that activation of P_2U purinergic receptors can propagate intercellular calcium, and describe a novel Cx43-dependent mechanism for calcium wave propagation that does not require release of intracellular calcium stores by IP₃. These studies suggest that gap junction communication mediated by either Cx43 or Cx45 does not allow passage of IP₃ well enough to elicit release of intracellular calcium stores in neighboring cells.     

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.     

A PKC wave follows the calcium wave after activation of Xenopus eggs

Calcium waves are well-known hallmarks of egg activation that trigger resumption of the cell cycle and development of the embryo. These waves rapidly and efficiently assure that activation signals are transmitted to all regions of the egg. Although the mechanism by which the calcium wave propagates across an egg as large as that of Xenopus is not known, two models prevail. One model is a wave of calcium-induced calcium release (CICR) and the other is propagation by inositol-induced calcium release (IICR). IICR requires a wave of phosphatidylinositol 4,5-bisphosphate (PIP2) hydrolysis, generating two second messengers, IP3, which then releases calcium and DAG, which activates protein kinase C (PKC). We show here that a wave of PKC-green fluorescent protein travels across the egg immediately following, and at the same velocity as, the calcium wave. This is the first example of a PKC wave in a vertebrate egg and supports the IICR model of wave propagation.     

Slow and fast invasion waves in a model of acid-mediated tumour growth

This work is concerned with a reaction-diffusion system that has been proposed as a model to describe acid-mediated cancer invasion. More precisely, we consider the properties of travelling waves that can be supported by such a system, and show that a rich variety of wave propagation dynamics, both fast and slow, is compatible with the model. In particular, asymptotic formulae for admissible wave profiles and bounds on their wave speeds are provided.     

Periodic travelling waves in cyclic predator–prey systems

Predation is an established cause of cycling in prey species. Here, the ability of predation to explain periodic travelling waves in prey populations, which have recently been found in a number of spatiotemporal field studies, is examined. The nature of periodic waves in these systems, and the way in which they can be generated by the invasion of predators into a prey population is discussed. A theoretical calculation that predicts, as a function of two parameter ratios, whether such an invasion will lead to a stable periodic travelling wave that would be observed in practice is presented ‐ the alternative outcome is spatiotemporal chaos. The calculation also predicts quantitative details of the periodic waves, such as speed and amplitude. The results give new insights into the types of predator‐prey systems in which one would expect to see periodic travelling waves following an invasion by predators.