A membrane model for cytosolic calcium oscillations. A study using Xenopus oocytes.

Cytosolic calcium oscillations occur in a wide variety of cells and are involved in different cellular functions. We describe these calcium oscillations by a mathematical model based on the putative electrophysiological properties of the endoplasmic reticulum (ER) membrane. The salient features of our membrane model are calcium-dependent calcium channels and calcium pumps in the ER membrane, constant entry of calcium into the cytosol, calcium dependent removal from the cytosol, and buffering by cytoplasmic calcium binding proteins. Numerical integration of the model allows us to study the fluctuations in the cytosolic calcium concentration, the ER membrane potential, and the concentration of free calcium binding sites on a calcium binding protein. The model demonstrates the physiological features necessary for calcium oscillations and suggests that the level of calcium flux into the cytosol controls the frequency and amplitude of oscillations. The model also suggests that the level of buffering affects the frequency and amplitude of the oscillations. The model is supported by experiments indirectly measuring cytosolic calcium by calcium-induced chloride currents in Xenopus oocytes as well as cytosolic calcium oscillations observed in other preparations.     

Modeling of Glucose-Induced cAMP Oscillations in Pancreatic β Cells: cAMP Rocks when Metabolism Rolls

Recent advances in imaging technology have revealed oscillations of cyclic adenosine monophosphate (cAMP) in insulin-secreting cells. These oscillations may be in phase with cytosolic calcium oscillations or out of phase. cAMP oscillations have previously been modeled as driven by oscillations in calcium, based on the known dependence of the enzymes that generate cAMP (adenylyl cyclase) and degrade it (phosphodiesterase). However, cAMP oscillations have also been reported to occur in the absence of calcium oscillations. Motivated by similarities between the properties of cAMP and metabolic oscillations in pancreatic β cells, we propose here that in addition to direct control by calcium, cAMP is controlled by metabolism. Specifically, we hypothesize that AMP inhibits adenylyl cyclase. We incorporate this hypothesis into the dual oscillator model for β cells, in which metabolic (glycolytic) oscillations cooperate with modulation of ion channels and metabolism by calcium. We show that the combination of oscillations in AMP and calcium in the dual oscillator model can account for the diverse oscillatory patterns that have been observed, as well as for experimental perturbations of those patterns. Predictions to further test the model are provided.     

Spatial phase pattern in one-dimensional arrays of limit cycle oscillators with discrete coupling

We present a model of limit cycle oscillators for collective oscillations in intracellular calcium concentration in cell communities. A phase-dependent discrete coupling between nearest neighbors is introduced into the model on the basis of the experimental observation that intercellular transmission of calcium or calcium mobilizing messenger is effected by gap junction and gap junctional permeability is affected by intracellular calcium concentration. The spatial phase pattern of several clusters in which oscillations are in phase is found with the phase-dependent discrete coupling.     

Modelling mechanism of calcium oscillations in pancreatic acinar cells

We present a simple model for calcium oscillations in the pancreatic acinar cells. This model is based on the calcium release from two receptors, inositol trisphosphate receptors (IPR) and ryanodine receptors (RyR) through the process of calcium induced calcium release (CICR). In pancreatic acinar cells, when the Ca²⁺ concentration increases, the mitochondria uptake it very fast to restrict Ca²⁺ response in the cell. Afterwards, a much slower release of Ca²⁺ from the mitochondria serves as a calcium supply in the cytosol which causes calcium oscillations. In this paper we discuss a possible mechanism for calcium oscillations based on the interplay among the three calcium stores in the cell: the endoplasmic reticulum (ER), mitochondria and cytosol. Our model predicts that calcium shuttling between ER and mitochondria is a pacemaker role in the generation of Ca²⁺oscillations. We also consider the calcium dependent production and degradation of (1,4,5) inositol-trisphosphate (IP₃), which is a key source of intracellular calcium oscillations in pancreatic acinar cells. In this study we are able to predict the different patterns of calcium oscillations in the cell from sinusoidal to raised-baseline, high frequency and low-frequency baseline spiking.     

Modeling the interrelations between the calcium oscillations and ER membrane potential oscillations

A refined electrochemical model accounting for intracellular calcium oscillations and their interrelations with oscillations of the potential difference across the membrane of the endoplasmic reticulum (ER) or other intracellular calcium stores is established. The ATP dependent uptake of Ca2+ from the cytosol into the ER, the Ca2+ release from the ER through channels following a calcium-induced calcium release mechanism, and a potential-dependent Ca2+ leak flux out of the ER are included in the model and described by plausible rate laws. The binding of calcium to specific proteins such as calmodulin is taken into account. The quasi-electroneutrality condition allows us to express the transmembrane potential in terms of the concentrations of cytosolic calcium and free binding sites on proteins, which are the two independent variables of the model. We include monovalent ions in the model, because they make up a considerable portion in the balance of electroneutrality. As the permeability of the endoplasmic membrane for these ions is much higher than that for calcium ions, we assume the former to be in Nernst equilibrium. A stability analysis of the steady-state solutions (which are unique or multiple depending on parameter values) is carried out and the Hopf bifurcation leading from stable steady states to self-sustained oscillations is analysed with the help of appropriate mathematical techniques. The oscillations obtained by numerical integration exhibit the typical spike-like shape found in experiments and reasonable values of frequency and amplitude. The model describes the process of switching between stationary and pulsatile regimes as well as changes in oscillation frequency upon parameter changes. It turns out that calcium oscillations can arise without a permanent influx of calcium into the cell, when a calcium-buffering system such as calmodulin is included.     

Mechanical and biochemical modeling of cortical oscillations in spreading cells

Actomyosin-based cortical contractility is a common feature of eukaryotic cells and is involved in cell motility, cell division, and apoptosis. In nonmuscle cells, oscillations in contractility are induced by microtubule depolymerization during cell spreading. We developed an ordinary differential equation model to describe this behavior. The computational model includes 36 parameters. The values for all but two of the model parameters were taken from experimental measurements found in the literature. Using these values, we demonstrate that the model generates oscillatory behavior consistent with current experimental observations. The rhythmic behavior occurs because of the antagonistic effects of calcium-induced contractility and stretch-activated calcium channels. The model makes several experimentally testable predictions: 1), buffering intracellular calcium increases the period and decreases the amplitude of cortical oscillations; 2), increasing the number or activity of stretch activated channels leads to an increase in period and amplitude of cortical oscillations; 3), inhibiting Ca²⁺ pump activity increases the period and amplitude of oscillations; and 4), a threshold exists for the calcium concentration below which oscillations cease.     

Functional Modeling of the Shift in Cellular Calcium Dynamics at the Onset of Synchronization in Smooth Muscle Cells

In the present paper we address the nature of synchronization properties found in populations of mesenteric artery smooth muscle cells. We present a minimal model of the onset of synchronization in the individual smooth muscle cell that is manifested as a transition from calcium waves to whole-cell calcium oscillations. We discuss how different types of ion currents may influence both amplitude and frequency in the regime of whole-cell oscillations. The model may also explain the occurrence of mixed-mode oscillations and chaotic oscillations frequently observed in the experimental system.     

On the encoding and decoding of calcium signals in hepatocytes

Many different agonists use calcium as a second messenger. Despite intensive research in intracellular calcium signalling it is an unsolved riddle how the different types of information represented by the different agonists, is encoded using the universal carrier calcium. It is also still not clear how the information encoded is decoded again into the intracellular specific information at the site of enzymes and genes. After the discovery of calcium oscillations, one likely mechanism is that information is encoded in the frequency, amplitude and waveform of the oscillations. This hypothesis has received some experimental support. However, the mechanism of decoding of oscillatory signals is still not known. Here, we study a mechanistic model of calcium oscillations, which is able to reproduce both spiking and bursting calcium oscillations. We use the model to study the decoding of calcium signals on the basis of co-operativity of calcium binding to various proteins. We show that this co-operativity offers a simple way to decode different calcium dynamics into different enzyme activities.     

Hormone-induced calcium oscillations in liver cells can be explained by a simple one pool model

Hormone-induced oscillations of the free intracellular calcium concentration are thought to be relevant for frequency encoding of hormone signals. In liver cells, such Ca2+ oscillations occur in response to stimulation by hormones acting via phosphoinositide breakdown. This observation may be explained by cooperative, positive feedback of Ca2+ on its own release from one inositol 1,4,5-trisphosphate-sensitive pool, obviating oscillations of inositol 1,4,5-trisphosphate. The kinetic rate laws of the associated model have a mathematical structure reminiscent of the Brusselator, a hypothetical chemical model involving a rather improbable trimolecular reaction step, thus giving a realistic biological interpretation to this hallmark of dissipative structures. We propose that calmodulin is involved in mediating this cooperativity and positive feedback, as suggested by the presented experiments. For one, hormone-induced calcium oscillations can be inhibited by the (nonphenothiazine) calmodulin antagonists calmidazolium or CGS 9343 B. Alternatively, in cells overstimulated by hormone, as characterized by a non-oscillatory elevated Ca2+ concentration, these antagonists could again restore sustained calcium oscillations. The experimental observations, including modulation of the oscillations by extracellular calcium, were in qualitative agreement with the predictions of our mathematical model.     

Ca2+ dynamics in a population of smooth muscle cells: modeling the recruitment and synchronization

Many experimental studies have shown that arterial smooth muscle cells respond with cytosolic calcium rises to vasoconstrictor stimulation. A low vasoconstrictor concentration gives rise to asynchronous spikes in the calcium concentration in a few cells (asynchronous flashing). With a greater vasoconstrictor concentration, the number of smooth muscle cells responding in this way increases (recruitment) and calcium oscillations may appear. These oscillations may eventually synchronize and generate arterial contraction and vasomotion. We show that these phenomena of recruitment and synchronization naturally emerge from a model of a population of smooth muscle cells coupled through their gap junctions. The effects of electrical, calcium, and inositol 1,4,5-trisphosphate coupling are studied. A weak calcium coupling is crucial to obtain a synchronization of calcium oscillations and the minimal required calcium permeability is deduced. Moreover, we note that an electrical coupling can generate oscillations, but also has a desynchronizing effect. Inositol 1,4,5-trisphosphate diffusion does not play an important role to achieve synchronization. Our model is validated by published in vitro experiments obtained on rat mesenteric arterial segments.     

Switching from simple to complex oscillations in calcium signaling

We present a new model for calcium oscillations based on experiments in hepatocytes. The model considers feedback inhibition on the initial agonist receptor complex by calcium and activated phospholipase C, as well as receptor type-dependent self-enhanced behavior of the activated G(alpha) subunit. It is able to show simple periodic oscillations and periodic bursting, and it is the first model to display chaotic bursting in response to agonist stimulations. Moreover, our model offers a possible explanation for the differences in dynamic behavior observed in response to different agonists in hepatocytes.     

Ryanodine block of calcium oscillations in heart muscle and the sodium-tension relationship

The control of tension is examined in cardiac Purkinje fibers. We show, in accordance with earlier results (14, 15), that tension produced by this preparation is a steep power function of intracellular sodium. With the aid of ryanodine, a pharmacological agent that blocks spontaneous and spatially asynchronous calcium release from the sarcoplasmic reticulum (SR), we investigate the influence that such calcium fluctuations have on tension. We find that even when we control for alterations of intracellular pH, the presence of such fluctuations reduces the dependence of tension on intracellular sodium. We present a simple model that can explain how the presence of oscillations of intracellular calcium leads to a reduction of the slope of the tension-calcium relationship. We show, furthermore, that when the oscillations are spatially asynchronous, this reduction of slope is even greater. The modeling takes account of the known relationships between tension and calcium and tension and sarcomere length. We conclude that the effect of ryanodine to steepen the tension-sodium relationship can be explained by ryanodine's blocking calcium release from the SR, thereby abolishing oscillations of intracellular calcium.     

Modelling the transition from simple to complex Ca2+oscillations in pancreatic acinar cells

A mathematical model is proposed which systematically investigates complex calcium oscillations in pancreatic acinar cells. This model is based on calcium-induced calcium release via inositol trisphosphate receptors (IPR) and ryanodine receptors (RyR) and includes calcium modulation of inositol (1,4,5) trisphosphate (IP₃) levels through feedback regulation of degradation and production. In our model, the apical and the basal regions are separated by a region containing mitochondria, which is capable of restricting Ca²⁺ responses to the apical region. We were able to reproduce the observed oscillatory patterns, from baseline spikes to sinusoidal oscillations. The model predicts that calcium-dependent production and degradation of IP₃ is a key mechanism for complex calcium oscillations in pancreatic acinar cells. A partial bifurcation analysis is performed which explores the dynamic behaviour of the model in both apical and basal regions.     

Agonist-dependent phosphorylation of the inositol 1,4,5-trisphosphate receptor: A possible mechanism for agonist-specific calcium oscillations in pancreatic acinar cells.

The properties of inositol 1,4,5-trisphosphate (IP3)-dependent intracellular calcium oscillations in pancreatic acinar cells depend crucially on the agonist used to stimulate them. Acetylcholine or carbachol (CCh) cause high-frequency (10-12-s period) calcium oscillations that are superimposed on a raised baseline, while cholecystokinin (CCK) causes long-period (>100-s period) baseline spiking. We show that physiological concentrations of CCK induce rapid phosphorylation of the IP3 receptor, which is not true of physiological concentrations of CCh. Based on this and other experimental data, we construct a mathematical model of agonist-specific intracellular calcium oscillations in pancreatic acinar cells. Model simulations agree with previous experimental work on the rates of activation and inactivation of the IP3 receptor by calcium (DuFour, J.-F., I.M. Arias, and T.J. Turner. 1997. J. Biol. Chem. 272:2675-2681), and reproduce both short-period, raised baseline oscillations, and long-period baseline spiking. The steady state open probability curve of the model IP3 receptor is an increasing function of calcium concentration, as found for type-III IP3 receptors by Hagar et al. (Hagar, R.E., A.D. Burgstahler, M.H. Nathanson, and B.E. Ehrlich. 1998. Nature. 396:81-84). We use the model to predict the effect of the removal of external calcium, and this prediction is confirmed experimentally. We also predict that, for type-III IP3 receptors, the steady state open probability curve will shift to lower calcium concentrations as the background IP3 concentration increases. We conclude that the differences between CCh- and CCK-induced calcium oscillations in pancreatic acinar cells can be explained by two principal mechanisms: (a) CCK causes more phosphorylation of the IP3 receptor than does CCh, and the phosphorylated receptor cannot pass calcium current; and (b) the rate of calcium ATPase pumping and the rate of calcium influx from the outside the cell are greater in the presence of CCh than in the presence of CCK.     

Calcium oscillations in a triplet of pancreatic acinar cells

We use a mathematical model of calcium dynamics in pancreatic acinar cells to investigate calcium oscillations in a ring of three coupled cells. A connected group of cells is modeled in two different ways: 1), as coupled point oscillators, each oscillator being described by a spatially homogeneous model; and 2), as spatially distributed cells coupled along their common boundaries by gap-junctional diffusion of inositol trisphosphate and/or calcium. We show that, although the point-oscillator model gives a reasonably accurate general picture, the behavior of the spatially distributed cells cannot always be predicted from the simpler analysis; spatially distributed diffusion and cell geometry both play important roles in determining behavior. In particular, oscillations in which two cells are in synchrony, with the third phase-locked but not synchronous, appears to be more dominant in the spatially distributed model than in the point-oscillator model. In both types of model, intercellular coupling leads to a variety of synchronous, phase-locked, or asynchronous behaviors. For some parameter values there are multiple, simultaneous stable types of oscillation. We predict 1), that intercellular calcium diffusion is necessary and sufficient to coordinate the responses in neighboring cells; 2), that the function of intercellular inositol trisphosphate diffusion is to smooth out any concentration differences between the cells, thus making it easier for the diffusion of calcium to synchronize the oscillations; 3), that groups of coupled cells will tend to respond in a clumped manner, with groups of synchronized cells, rather than with regular phase-locked periodic intercellular waves; and 4), that enzyme secretion is maximized by the presence of a pacemaker cell in each cluster which drives the other cells at a frequency greater than their intrinsic frequency.     

Complex calcium oscillations and the role of mitochondria and cytosolic proteins

Intracellular calcium oscillations, which are oscillatory changes of cytosolic calcium concentration in response to agonist stimulation, are experimentally well observed in various living cells. Simple calcium oscillations represent the most common pattern and many mathematical models have been published to describe this type of oscillation. On the other hand, relatively few theoretical studies have been proposed to give an explanation of complex intracellular calcium oscillations, such as bursting and chaos. In this paper, we develop a new possible mechanism for complex calcium oscillations based on the interplay between three calcium stores in the cell: the endoplasmic reticulum (ER), mitochondria and cytosolic proteins. The majority ( approximately 80%) of calcium released from the ER is first very quickly sequestered by mitochondria. Afterwards, a much slower release of calcium from the mitochondria serves as the calcium supply for the intermediate calcium exchanges between the ER and the cytosolic proteins causing bursting calcium oscillations. Depending on the permeability of the ER channels and on the kinetic properties of calcium binding to the cytosolic proteins, different patterns of complex calcium oscillations appear. With our model, we are able to explain simple calcium oscillations, bursting and chaos. Chaos is also observed for calcium oscillations in the bursting mode.     

Birhythmicity, trirhythmicity and chaos in bursting calcium oscillations

We have analyzed various types of complex calcium oscillations. The oscillations are explained with a model based on calcium-induced calcium release (CICR). In addition to the endoplasmic reticulum as the main intracellular Ca2+ store, mitochondrial and cytosolic Ca2+ binding proteins are also taken into account. This model was previously proposed for the study of the physiological role of mitochondria and the cytosolic proteins in gene rating complex Ca2+ oscillations [1]. Here, we investigated the occurrence of different types of Ca2+ oscillations obtained by the model, i.e. simple oscillations, bursting, and chaos. In a bifurcation diagram, we have shown that all these various modes of oscillatory behavior are obtained by a change of only one model parameter, which corresponds to the physiological variability of an agonist. Bursting oscillations were studied in more detail because they express birhythmicity, trirhythmicity and chaotic behavior. Two different routes to chaos are observed in the model: in addition to the usual period doubling cascade, we also show intermittency. For the characterization of the chaotic behavior, we made use of return maps and Lyapunov exponents. The potential biological role of chaos in intracellular signaling is discussed.     

Model of the Ca2+ oscillator for shuttle streaming in Physarum polycephalum.

We propose a mechanism for the cytoplasmic Ca++ oscillator which is thought to power shuttle streaming in strands of the slime-mold Physarum polycephalum. The mechanism uses a phosphorylation-dephosphorylation cycle of myosin light chain kinase. This kinase is bistable if the kinase phosphorylation chain, through adenylate cyclase and cAMP, is activated by calcium. Relaxation oscillations can then occur if calcium is exchanged between the cytoplasm and internal vacuoles known to exist in physarum. As contractile activity in physarum myosin is inhibited by calcium, this model can give calcium oscillations 180 degrees out of phase with actin filament tension as observed. Oscillations of ATP concentration are correctly predicted to be in phase with the tension, provided the actomyosin cycling rate is comparable with ATPase rates for phosphorylation of the myosin light chain and its kinase.     

Intercellular communication via intracellular calcium oscillations

In this letter, we present the results of a simple model for intercellular communication via calcium oscillations, motivated in part by a recent experimental study. The model describes two cells (a "donor" and "sensor") whose intracellular dynamics involve a calcium-induced, calcium release process. The cells are coupled by assuming that the input of the sensor cell is proportional to the output of the donor cell. As one varies the frequency of calcium oscillations of the donor cell, the sensor cell passes through a sequence of N : M phase-locked regimes and exhibits a "Devil's staircase" behavior. Such a phase-locked response has been seen experimentally in pulsatile stimulation of single cells. We also study a stochastic version of the coupled two-cell model. We find that phase locking holds for realistic choices for the cell volume.