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.     

Ceramide increases mitochondrial free calcium levels via caspase 8 and Bid: role in initiation of cell death

We investigated how the mitochondrial phase of ceramide-mediated cell death is initiated in nerve growth factor (NGF)-differentiated PC12 cells. We distinguished three independent effects of ceramide: free radical production; a transient increase in cytosolic free calcium; and a long-lasting increase in mitochondrial free calcium. Only the latter led to cell death, which could be prevented by buffering of mitochondrial calcium with the calcium binding protein calbindin D-28K ectopically expressed in mitochondria. We showed that mitochondrial calcium did not increase as a result of the increase in cytosolic free calcium levels. Rather, it appears to derive from the endoplasmic reticulum (ER) since dantrolene, which inhibits release of calcium from ER into cytosol through ryanodine receptors, prevented the increase in cytosolic free calcium but potentiated the increase in mitochondrial free calcium. This suggests that a transfer of calcium occurs directly, or very locally, between the two organelles. This transfer implicated activation of caspase 8 and cleavage of its substrate Bid, a previously unknown function of these cell death intermediaries. The increase in mitochondrial free calcium was also responsible for the release of cytochrome c into the cytosol, underlining the critical role it plays in ceramide-mediated cell death.     

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.     

Spontaneous electrical and calcium oscillations in unstimulated pituitary gonadotrophs.

Single pituitary cells often fire spontaneous action potentials (APs), which are believed to underlie spiking fluctuations in cytosolic calcium concentration ([Ca2+]i). To address how these basal [Ca2+]i fluctuations depend on changes in plasma membrane voltage (V), simultaneous measurements of V and [Ca2+]i were performed in rat pituitary gonadotrophs. The data show that each [Ca2+]i spike is produced by the Ca2+ entry during a single AP. Using these and previously obtained patch-clamp data, we develop a quantitative mathematical model of this plasma membrane oscillator and the accompanying spatiotemporal [Ca2+]i oscillations. The model demonstrates that AP-induced [Ca2+]i spiking is prominent only in a thin shell layer neighboring the cell surface. This localized [Ca2+]i spike transiently activates the Ca2(+)- dependent K+ current resulting in a sharp afterhyperpolarization following each voltage spike. In accord with experimental observations, the model shows that the frequency and amplitude of the voltage spikes are highly sensitive to current injection and to the blocking of the Ca(2+)-sensitive current. Computations also predict that leaving the membrane channels intact, the firing rate can be modified by changing the Ca2+ handling parameters: the Ca2+ diffusion rate, the Ca2+ buffering capacity, and the plasma membrane Ca2+ pump rate. Finally, the model suggests reasons that spontaneous APs were seen in some gonadotrophs but not in others. This model provides a basis for further exploring how plasma membrane electrical activity is involved in the control of cytosolic calcium level in unstimulated as well as agonist-stimulated gonadotrophs.     

Uptake and release of Ca2+ by the endoplasmic reticulum contribute to the oscillations of the cytosolic Ca2+ concentration triggered by Ca2+ influx in the electrically excitable pancreatic B-cell.

The role of intracellular Ca2+ pools in oscillations of the cytosolic Ca2+ concentration ([Ca2+]c) triggered by Ca2+ influx was investigated in mouse pancreatic B-cells. [Ca2+]c oscillations occurring spontaneously during glucose stimulation or repetitively induced by pulses of high K+ (in the presence of diazoxide) were characterized by a descending phase in two components. A rapid decrease in [Ca2+]c coincided with closure of voltage-dependent Ca2+ channels and was followed by a slower phase independent of Ca2+ influx. Blocking the SERCA pump with thapsigargin or cyclopiazonic acid accelerated the rising phase of [Ca2+]c oscillations and increased their amplitude, which suggests that the endoplasmic reticulum (ER) rapidly takes up Ca2+. It also suppressed the slow [Ca2+]c recovery phase, which indicates that this phase corresponds to the slow release of Ca2+ that was taken up by the ER during the upstroke of the [Ca2+]c transient. Glucose promoted the buffering capacity of the ER and amplified the slow [Ca2+]c recovery phase. The slow phase induced by high K+ pulses was not affected by modulators of Ca2+- or inositol 1,4,5-trisphosphate-induced Ca2+ release, did not involve a depolarization-induced Ca2+ release, and was also observed at the end of a rapid rise in [Ca2+]c triggered from caged Ca2+. It is attributed to passive leakage of Ca2+ from the ER. We suggest that the ER displays oscillations of the Ca2+ concentration ([Ca2+]ER) concomitant and parallel to [Ca2+]c. The observation that thapsigargin depolarizes the membrane of B-cells supports the proposal that the degree of Ca2+ filling of the ER modulates the membrane potential. Therefore, [Ca2+]ER oscillations occurring during glucose stimulation are likely to influence the bursting behavior of B-cells and eventually [Ca2+]c oscillations.     

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.     

Sensing and refilling calcium stores in an excitable cell.

Inositol 1,4,5-trisphosphate (IP3)-induced Ca2+ mobilization leads to depletion of the endoplasmic reticulum (ER) and an increase in Ca2+ entry. We show here for the gonadotroph, an excitable endocrine cell, that sensing of ER Ca2+ content can occur without the Ca2+ release-activated Ca2+ current (Icrac), but rather through the coupling of IP3-induced Ca2+ oscillations to plasma membrane voltage spikes that gate Ca2+ entry. Thus we demonstrate that capacitative Ca2+ entry is accomplished through Ca(2+)-controlled Ca2+ entry. We develop a comprehensive model, with parameter values constrained by available experimental data, to simulate the spatiotemporal behavior of agonist-induced Ca2+ signals in both the cytosol and ER lumen of gonadotrophs. The model combines two previously developed models, one for ER-mediated Ca2+ oscillations and another for plasma membrane potential-driven Ca2+ oscillations. Simulations show agreement with existing experimental records of store content, cytosolic Ca2+ concentration ([Ca2+]i), and electrical activity, and make a variety of new, experimentally testable predictions. In particular, computations with the model suggest that [Ca2+]i in the vicinity of the plasma membrane acts as a messenger for ER content via Ca(2+)-activated K+ channels and Ca2+ pumps in the plasma membrane. We conclude that, in excitable cells that do not express Icrac, [Ca2+]i profiles provide a sensitive mechanism for regulating net calcium flux through the plasma membrane during both store depletion and refilling.     

Computational model of excitable cell indicates ATP free energy dynamics in response to calcium oscillations are undampened by cytosolic ATP buffers

Mitochondria in excitable cells are recurrently exposed to pulsatile calcium gradients that activate cell function. Rapid calcium uptake by the mitochondria has previously been shown to cause uncoupling of oxidative phosphorylation. To test (i) if periodic nerve firing may cause oscillation of the cytosolic thermodynamic potential of ATP hydrolysis and (ii) if cytosolic adenylate (AK) and creatine kinase (CK) ATP buffering reactions dampen such oscillations, a lumped kinetic model of an excitable cell capturing major aspects of the physiology has been developed. Activation of ATP metabolism by low-frequency calcium pulses caused large oscillation of the cytosolic, but not mitochondrial ATP/ADP, ratio. This outcome was independent of net ATP synthesis or hydrolysis during mitochondrial calcium uptake. The AK/CK ATP buffering reactions dampened the amplitude and rate of cytosolic ATP/ADP changes on a timescale of seconds, but not milliseconds. These model predictions suggest that alternative sources of capacitance in neurons and striated muscles should be considered to protect ATP-free energy-driven cell functions.     

Involvement of endoplasmic reticulum in paclitaxel-induced apoptosis

The participation of the mitochondrial pathway in paclitaxel-induced apoptosis has been well documented. After addition of paclitaxel to U937 cells, however, we observed an early expression of five endoplasmic reticulum (ER) stress response genes that preceded the release of cytochrome c from the mitochondria and the cleavage of the caspases. Involvement of the ER was supported by the following evidence. Paclitaxel treatment not only activated calpain and caspase-4, but also induced a gradual increase in the cytosolic Ca(2+) concentration at 3-6 h. Paclitaxel-induced apoptosis can be inhibited by the calpain inhibitor calpeptin and IP(3) receptor inhibitors. Either buffering of the cytosolic Ca(2+) or inhibition of mitochondrial calcium uptake reduced BiP expression. These inhibitors also reduced mitochondrial apoptotic signals, such as mitochondrion membrane potential disruption, cytochrome c release and eventually reduced the death of U937 cells. Paclitaxel-induced Bax/Bak translocation to the ER and Bax dimerization on the ER membrane occurred within 3 h, which led to a Ca(2+) efflux into cytosol. Moreover, we found that cytochrome c translocated to the ER after releasing from mitochondria and then interacted with the IP(3) receptor at 12-15 h. This phenomenon has been known to amplify apoptotic signaling. Taken together, ER would seem to contribute to paclitaxel-induced apoptosis via both the early release of Ca(2+) and the late amplification of mitochondria-mediated apoptotic signals.     

Contribution of the endoplasmic reticulum to the glucose-induced [Ca(2+)](c) response in mouse pancreatic islets

Thapsigargin (TG), a blocker of Ca(2+) uptake by the endoplasmic reticulum (ER), was used to evaluate the contribution of the organelle to the oscillations of cytosolic Ca(2+) concentration ([Ca(2+)](c)) induced by repetitive Ca(2+) influx in mouse pancreatic beta-cells. Because TG depolarized the plasma membrane in the presence of glucose alone, extracellular K(+) was alternated between 10 and 30 mM in the presence of diazoxide to impose membrane potential (MP) oscillations. In control islets, pulses of K(+), mimicking regular MP oscillations elicited by 10 mM glucose, induced [Ca(2+)](c) oscillations whose nadir remained higher than basal [Ca(2+)](c). Increasing the depolarization phase of the pulses while keeping their frequency constant (to mimic the effects of a further rise of the glucose concentration on MP) caused an upward shift of the nadir of [Ca(2+)](c) oscillations that was reproduced by raising extracellular Ca(2+) (to increase Ca(2+) influx) without changing the pulse protocol. In TG-pretreated islets, the imposed [Ca(2+)](c) oscillations were of much larger amplitude than in control islets and occurred on basal levels. During intermittent trains of depolarizations, control islets displayed mixed [Ca(2+)](c) oscillations characterized by a summation of fast oscillations on top of slow ones, whereas no progressive summation of the fast oscillations was observed in TG-pretreated islets. In conclusion, the buffering capacity of the ER in pancreatic beta-cells limits the amplitude of [Ca(2+)](c) oscillations and may explain how the nadir between oscillations remains above baseline during regular oscillations or gradually increases during mixed [Ca(2+)](c) oscillations, two types of response observed during glucose stimulation.     

Effects of inositol 1,4,5-trisphosphate receptor-mediated intracellular stochastic calcium oscillations on activation of glycogen phosphorylase

In various cell types cytosolic calcium (Ca(2+)) is an important regulator. The possible role of Ca(2+) release from the inositol 1,4,5-trisphosphate (IP(3)) receptor channel in the regulation of the phosphorylation-dephosphorylation cycle process involved in glycogen degradation by glycogen phosphorylase have theoretically investigated by using the Li-Rinzel model for cytosolic Ca(2+) oscillations. For the case of deterministic cytosolic Ca(2+) oscillations, there exists an optimal frequency of cytosolic Ca(2+) oscillations at which the average fraction of active glycogen phosphorylase reaches a maximum value, and a mutation for the average fraction of active glycogen phosphorylase occurs at the higher bifurcation point of Ca(2+) oscillations. For the case of stochastic cytosolic Ca(2+) oscillations, the fraction of active phosphorylase is strongly affected by the number of IP(3) receptor channels and the level of IP(3) concentration. Small number of IP(3) receptor channels can potentiate the sensitivity of the activity of glycogen phosphorylase. The average frequency and amplitude of active phosphorylase stochastic oscillations are increased with the level of increasing IP(3) stimuli. The various distributions for the amplitude of active glycogen phosphorylase oscillations in parameters plane are discussed.     

Interaction of calmodulin with Sec61α limits Ca2+ leakage from the endoplasmic reticulum

In eukaryotes, protein transport into the endoplasmic reticulum (ER) is facilitated by a protein-conducting channel, the Sec61 complex. The presence of large, water-filled pores with uncontrolled ion permeability, as formed by Sec61 complexes in the ER membrane, would seriously interfere with the regulated release of calcium from the ER lumen into the cytosol, an essential mechanism for intracellular signalling. We identified a calmodulin (CaM)-binding motif in the cytosolic N-terminus of mammalian Sec61α that bound CaM but not Ca2+-free apocalmodulin with nanomolar affinity and sequence specificity. In single-channel measurements, CaM potently mediated Sec61-channel closure in Ca2+-dependent manner. At the cellular level, two different CaM antagonists stimulated calcium release from the ER through Sec61 channels. However, protein transport into microsomes was not modulated by Ca2+-CaM. Molecular modelling of the ribosome/Sec61/CaM complexes supports the view that simultaneous ribosome and CaM binding to the Sec61 complex may be possible. Overall, CaM is involved in limiting Ca2+ leakage from the ER.     

Characterization of estrogen and antiestrogen binding to the cytosol and microsomes of breast tumors

The binding of [3H]estradiol and [3H]hydroxytamoxifen to the cytosol and microsomal fractions of several human breast tumors was investigated. By washing microsomal membranes with a KCl-free or a KCl-containing medium we could distinguish between intrinsic, extrinsic and contaminant estradiol binding sites in these membranes. We observed that treatment of the microsomes with low salt medium removes about 80% of the total estradiol binding sites, whereas 20% are not extractable. The concentration of unextractable [3H]estradiol binding sites in the microsomes varies in proportion to the level of cytosolic estrogen receptors (ER). About 10% of the total extranuclear specific estrogen binding sites was consistently found tightly associated to the microsomal fraction, which displays an affinity for estradiol (Kd = 0.1-0.6 nM) similar to that of the cytosolic ER. The displacement of [3H]estradiol with unlabeled hormone or with the antiestrogens, nafoxidine, enclomiphene and tamoxifen (TAM) exhibits identical IC50 values either in the cytosol or in the microsomal membranes. On the other hand, the microsomal fraction of breast tumors also binds [3H]hydroxyTAM, but with higher capacity and lower affinity than those of the cytosolic fraction. Furthermore, we did not observe correlation between the concentrations of ER and of antiestrogen binding sites (AEBS) in the tumors. These results indicate that microsomal membranes of human breast tumors contain estrogen binding sites which may be related to the cytosol ER recycling and that specific AEBS are predominantly localized in this membrane system. Furthermore, it is shown that the magnitude of estradiol binding to microsomes depends on the ER positive degree of the tumors, whereas the magnitude of the antiestrogen binding to the microsomes is independent of the ER status of the tumors.     

Modulation of neuronal excitability by intracellular calcium buffering: from spiking to bursting

We have investigated the detailed regulation of neuronal firing pattern by the cytosolic calcium buffering capacity using a combination of mathematical modeling and patch-clamp recording in acute slice. Theoretical results show that a high calcium buffer concentration alters the characteristic regular firing of cerebellar granule cells and that a transition to various modes of oscillations occurs, including bursting. Using bifurcation analysis, we show that this transition from spiking to bursting is a consequence of the major slowdown of calcium dynamics. Patch-clamp recordings on cerebellar granule cells loaded with a high concentration of the fast calcium buffer BAPTA (15 mM) reveal dramatic alterations in their excitability as compared to cells loaded with 0.15 mM BAPTA. In high calcium buffering conditions, granule cells exhibit all bursting behaviors predicted by the model whereas bursting is never observed in low buffering. These results suggest that cytosolic calcium buffering capacity can tightly modulate neuronal firing patterns leading to generation of complex patterns and therefore that calcium-binding proteins may play a critical role in the non-synaptic plasticity and information processing in the central nervous system.     

A mathematical model for LH release in response to continuous and pulsatile exposure of gonadotrophs to GnRH

In a previous study, a model was developed to investigate the release of luteinizing hormone (LH) from pituitary cells in response to a short pulse of gonadotropin-releasing hormone (GnRH). The model included: binding of GnRH to its receptor (R), dimerization and internalization of the hormone receptor complex, interaction with a G protein, production of inositol 1,4,5-trisphosphate (IP₃), release of calcium from the endoplasmic reticulum (ER), entrance of calcium into the cytosol via voltage gated membrane channels, pumping of calcium out of the cytosol via membrane and ER pumps, and release of LH. The extended model, presented in this paper, also includes the following physiologically important phenomena: desensitization of calcium channels; internalization of the dimerized receptors and recycling of some of the internalized receptors; an increase in G _q concentration near the plasma membrane in response to receptor dimerization; and basal rates of synthesis and degradation of the receptors. With suitable choices of the parameters, good agreement with a variety of experimental data of the LH release pattern in response to pulses of various durations, repetition rates, and concentrations of GnRH were obtained. The mathematical model allows us to assess the effects of internalization and desensitization on the shapes and time courses of LH response curves.     

Calmodulin regulation of the calcium-leak channel Sec61 is unique to vertebrates

In eukaryotes, protein transport into the endoplasmic reticulum (ER) is facilitated by a protein-conducting channel, the Sec61 complex. The presence of large, water-filled pores with uncontrolled ion permeability, such as those formed by Sec61 complexes in the ER membrane, would interfere with the regulated release of calcium from the ER lumen into the cytosol, an essential mechanism of intracellular signaling. We identified a calmodulin (CaM) binding motif in the cytosolic N-terminus of Sec61α from Canis familiaris that binds CaM, but not Ca(2+)-free apo-CaM, with nanomolar affinity and sequence specificity. In single channel lipid bilayer measurements, CaM potently mediated Sec61-channel closure in a Ca(2+)-dependent manner. No functional CaM binding motif was identified in the corresponding region of Sec61p from Saccharomyces cerevisiae, and no channel closure occurred in the presence of CaM and Ca(2+). Therefore, CaM binding to the cytosolic N-terminus of Sec61α is involved in limiting Ca(2+)-leakage from the ER in C. familiaris but not S. cerevisiae.     

A Simulation Study of Calcium Dynamics Features Caused by Exchange between the Cytosol and Organellar Stores of Neurons

The objects of the study were single-compartment mathematical models corresponding to a fragment of the dendrite of a cerebellar Purkinje neuron. The fragments contained the mitochondria (model 1) or a cistern of the endoplasmic reticulum, ER (model 2), functioning as calcium stores. With simulating single excitatory synaptic actions, we examined the dependence of the dynamics of intracellular Ca²⁺ levels on the maximum rate of Ca²⁺ exchange between the cytosol and these stores, as well as on the intensity of the diffusion flow into adjacent organelle-free regions. The plasma membrane of the compartment had ion channels (including those of the synaptic current) and the calcium pump characteristic of the mentioned neurons. The model equations took into account Ca²⁺ exchange between the cytosol, extracellular environment, and organellar stores, as well as the diffusion process. In model 1, the mitochondria exchanged Ca²⁺ with the cytosol through the uniporter and sodium-calcium exchanger; mitochondrial processes, such as the tricarboxylic acid cycle and aerobic cellular respiration, were also included. In model 2, the ER membrane had the calcium pump, leak channels, and channels of calcium-induced and inositol-3-phosphate-dependent Ca²⁺ release. The stores (mitochondria or ER) occupied 36% of the total volume of the compartment. An increase in the maximum rate of calcium exchange with the stores led to a proportional decrease in the peak Ca²⁺ concentrations in the cytosol ([Ca²⁺]_i), more pronounced in the case of the ER; the Ca²⁺ concentration in both types of stores increased significantly. Due to the higher storage rate, the ER was able to absorb several times greater amounts of Ca²⁺ than the mitochondria did. With smaller diffusion flux (e.g., similarly to the case of diffusion from a larger-sized head into the neck of the dendritic spine), the intensity of cytosolic transients increased at fixed kinetics of flux exchange with the stores. Therefore, the organellar stores can significantly modulate not only the intensity but also the time course of changes in the intracellular Ca²⁺ levels.     

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.