Dual sensitivity of sarcoplasmic/endoplasmic Ca2+-ATPase to cytosolic and endoplasmic reticulum Ca2+ as a mechanism of modulating cytosolic Ca2+ oscillations

Heat-hypersensitive mutants of ryanodine receptor type 1 revealed by microscopic heating

Malignant hyperthermia (MH) is a life-threatening disorder caused largely by mutations in ryanodine receptor type 1 (RyR1) Ca²⁺-release channels. Enhanced Ca²⁺ release through the mutant channels induces excessive heat development upon exposure to volatile anesthetics. However, the mechanism by which Ca²⁺ release is accelerated at an elevated temperature is yet to be identified. Fluorescence Ca²⁺ imaging with rapid heating by an infrared laser beam provides direct evidence that heat induces Ca²⁺ release through the RyR1 channel. And the mutant channels are more heat sensitive than the wild-type channels, thereby causing an increase in the cytosolic Ca²⁺ concentration in mutant cells. It is likely that the heat-induced Ca²⁺ release participates as an enhancer in the cellular mechanism of MH.

Thermoregulation is an important aspect of human homeostasis, and high temperatures pose serious stresses for the body. Malignant hyperthermia (MH) is a life-threatening disorder in which body temperature can rise to a lethal level. Here we employ an optically controlled local heat-pulse method to manipulate the temperature in cells with a precision of less than 1 °C and find that the mutants of ryanodine receptor type 1 (RyR1), a key Ca²⁺ release channel underlying MH, are heat hypersensitive compared with the wild type (WT). We show that the local heat pulses induce an intracellular Ca²⁺ burst in human embryonic kidney 293 cells overexpressing WT RyR1 and some RyR1 mutants related to MH. Fluorescence Ca²⁺ imaging using the endoplasmic reticulum–targeted fluorescent probes demonstrates that the Ca²⁺ burst originates from heat-induced Ca²⁺ release (HICR) through RyR1-mutant channels because of the channels’ heat hypersensitivity. Furthermore, the variation in the heat hypersensitivity of four RyR1 mutants highlights the complexity of MH. HICR likewise occurs in skeletal muscles of MH model mice. We propose that HICR contributes an additional positive feedback to accelerate thermogenesis in patients with MH.

IRAG2 Interacts with IP3-Receptor Types 1, 2, and 3 and Regulates Intracellular Ca2+ in Murine Pancreatic Acinar Cells

The inositol 1,4,5-triphosphate receptor-associated 2 (IRAG2) is also known as Jaw1 or lymphoid-restricted membrane protein (LRMP) and shares homology with the inositol 1,4,5-triphosphate receptor-associated cGMP kinase substrate 1 (IRAG1). IRAG1 interacts with inositol trisphosphate receptors (IP₃ receptors /IP₃R) via its coiled-coil domain and modulates Ca²⁺ release from intracellular stores. Due to the homology of IRAG1 and IRAG2, especially in its coiled-coil domain, it is possible that IRAG2 has similar interaction partners like IRAG1 and that IRAG2 also modulates intracellular Ca²⁺ signaling. In our study, we localized IRAG2 in pancreatic acinar cells of the exocrine pancreas, and we investigated the interaction of IRAG2 with IP₃ receptors and its impact on intracellular Ca²⁺ signaling and exocrine pancreatic function, like amylase secretion. We detected the interaction of IRAG2 with different subtypes of IP₃R and altered Ca²⁺ release in pancreatic acinar cells from mice lacking IRAG2. IRAG2 deficiency decreased basal levels of intracellular Ca²⁺, suggesting that IRAG2 leads to activation of IP₃R under unstimulated basal conditions. Moreover, we observed that loss of IRAG2 impacts the secretion of amylase. Our data, therefore, suggest that IRAG2 modulates intracellular Ca²⁺ signaling, which regulates exocrine pancreatic function.

Ca2+ dysregulation in the endoplasmic reticulum related to Alzheimer's disease: A review on experimental progress and computational modeling

Alzheimer's disease (AD) is a devastating, incurable neurodegenerative disease affecting millions of people worldwide. Dysregulation of intracellular Ca(2+) signaling has been observed as an early event prior to the presence of clinical symptoms of AD and is believed to be a crucial factor contributing to its pathogenesis. The progressive and sustaining increase in the resting level of cytosolic Ca(2+) will affect downstream activities and neural functions. This review focuses on the issues relating to the increasing Ca(2+) release from the endoplasmic reticulum (ER) observed in AD neurons. Numerous research papers have suggested that the dysregulation of ER Ca(2+) homeostasis is associated with mutations in the presenilin genes and amyloid-β oligomers. These disturbances could happen at many different points in the signaling process, directly affecting ER Ca(2+) channels or interfering with related pathways, which makes it harder to reveal the underlying mechanisms. This review paper also shows that computational modeling is a powerful tool in Ca(2+) signaling studies and discusses the progress in modeling related to Ca(2+) dysregulation in AD research.

Altered Endoplasmic Reticulum Calcium Pump Expression during Breast Tumorigenesis

Endoplasmic reticulum calcium homeostasis is involved in several essential cell functions including cell proliferation, protein synthesis, stress responses or secretion. Calcium uptake into the endoplasmic reticulum is performed by Sarco/Endoplasmic Reticulum Calcium ATPases (SERCA enzymes). In order to study endoplasmic reticulum calcium homeostasis in situ in mammary tissue, in this work SERCA3 expression was investigated in normal breast and in its benign and malignant lesions in function of the cell type, degree of malignancy, and histological and molecular parameters of the tumors. Our data indicate, that although normal breast acinar epithelial cells express SERCA3 abundantly, its expression is strongly decreased already in very early non-malignant epithelial lesions such as adenosis, and remains low in lobular carcinomas. Whereas normal duct epithelium expresses significant amounts of SERCA3, its expression is decreased in several benign ductal lesions, as well as in ductal adenocarcinoma. The loss of SERCA3 expression is correlated with Elston-Ellis grade, negative hormone receptor expression or triple negative status in ductal carcinomas. The concordance between decreased SERCA3 expression and several histological, as well as molecular markers of ductal carcinogenesis indicates that endoplasmic reticulum calcium homeostasis is remodeled during tumorigenesis in the breast epithelium.

Highly cooperative dependence of sarco/endoplasmic reticulum calcium ATPase SERCA2a pump activity on cytosolic calcium in living cells

Sarco/endoplasmic reticulum (SR/ER) Ca(2+)-ATPase (SERCA) is an intracellular Ca(2+) pump localized on the SR/ER membrane. The role of SERCA in refilling intracellular Ca(2+) stores is pivotal for maintaining intracellular Ca(2+) homeostasis, and disturbed SERCA activity causes many disease phenotypes, including heart failure, diabetes, cancer, and Alzheimer disease. Although SERCA activity has been described using a simple enzyme activity equation, the dynamics of SERCA activity in living cells is still unknown. To monitor SERCA activity in living cells, we constructed an enhanced CFP (ECFP)- and FlAsH-tagged SERCA2a, designated F-L577, which retains the ATP-dependent Ca(2+) pump activity. The FRET efficiency between ECFP and FlAsH of F-L577 is dependent on the conformational state of the molecule. ER luminal Ca(2+) imaging confirmed that the FRET signal changes directly reflect the Ca(2+) pump activity. Dual imaging of cytosolic Ca(2+) and the FRET signals of F-L577 in intact COS7 cells revealed that SERCA2a activity is coincident with the oscillatory cytosolic Ca(2+) concentration changes evoked by ATP stimulation. The Ca(2+) pump activity of SERCA2a in intact cells can be expressed by the Hill equation with an apparent affinity for Ca(2+) of 0.41 ± 0.0095 μm and a Hill coefficient of 5.7 ± 0.73. These results indicate that in the cellular environment the Ca(2+) dependence of ATPase activation is highly cooperative and that SERCA2a acts as a rapid switch to refill Ca(2+) stores in living cells for shaping the intracellular Ca(2+) dynamics. F-L577 will be useful for future studies on Ca(2+) signaling involving SERCA2a activity.

Phe27Cys polymorphism of the human delta opioid receptor predisposes cells to compromised calcium signaling

A quarter of the human population with European background carries at least one allele of the OPRD1 gene that encodes the delta opioid receptor with cysteine at the amino acid position 27 (hδOR(Cys27)) instead of the evolutionary conserved phenylalanine (hδOR(Phe27)). The two variants have indistinguishable pharmacological properties but, importantly, hδOR(Cys27) differs from hδOR(Phe27) in having low maturation efficiency, lower stability at the cell surface and pronounced intracellular location. Both variants were previously shown to interact with the Sarco(endo)plasmic reticulum Ca²+ ATPase (SERCA) 2b in the early phase of their biosynthesis. We analyzed by pulse-chase assays, whether cellular signaling can affect hδOR(Cys27) maturation. Neither activation of the receptor by a δOR-specific agonist Leu-enkephalin, induction of intracellular calcium (Ca²+) release by ATP nor the direct stimulation of SERCA 2b by protein kinase C activation affected receptor maturation in HEK-293 cells. No signaling-mediated regulation of receptor maturation could therefore be demonstrated. Instead, we found by using single cell Ca²+ measurements that over-expression of hδOR(Cys27), but not hδOR(Phe27), compromised ATP-induced intracellular Ca²+-signaling. Furthermore, hδOR(Cys27) precursors showed slower dissociation from SERCA2b and hδOR(Cys27) expression caused down-regulation of the homocysteine-inducible endoplasmic reticulum-resident ubiquitin domain-like member 1 protein (HERP). We suggest that aging individuals with at least one hδOR(Cys27) encoding allele might have lowered threshold for Ca²+ dysregulation in neurons expressing hδOR.

The mechanism of agonist induced Ca2+ signalling in intact endothelial cells studied confocally in in situ arteries

In endothelial cells there remain uncertainties in the details of how Ca²⁺ signals are generated and maintained, especially in intact preparations. In particular the role of the sarco-endoplasmic reticulum Ca²⁺-ATPase (SERCA), in contributing to the components of agonist-induced signals is unclear.

The aim of this work was to increase understanding of the detailed mechanism of Ca²⁺ signalling in endothelial cells using real time confocal imaging of Fluo-4 loaded intact rat tail arteries in response to muscarinic stimulation. In particular we have focused on the role of SERCA, and its interplay with capacitative Ca²⁺ entry (CCE) and ER Ca²⁺ release and uptake. We have determined its contribution to the Ca²⁺ signal and how it varies with different physiological stimuli, including single and repeated carbachol applications and brief and prolonged exposures.

In agreement with previous work, carbachol stimulated a rise in intracellular Ca²⁺ in the endothelial cells, consisting of a rapid initial phase, then a plateau upon which oscillations of Ca²⁺ were superimposed, followed by a decline to basal Ca²⁺ levels upon carbachol removal. Our data support the following conclusions: (i) the size (amplitude and duration) of the Ca²⁺ spike and early oscillations are limited by SERCA activity, thus both are increased if SERCA is inhibited. (ii) SERCA activity is such that brief applications of carbachol do not trigger CCE, presumably because the fall in luminal Ca²⁺ is not sufficient to trigger it. However, longer applications sufficient to deplete the ER or even partial SERCA inhibition stimulate CCE. (iii) Ca²⁺ entry occurs via STIM-mediated CCE and SERCA contributes to the cessation of CCE. In conclusion our data show how SERCA function is crucial to shaping endothelial cell Ca signals and its dynamic interplay with both CCE and ER Ca releases.

Modulation of endoplasmic reticulum Ca2+ store filling by cyclic ADP-ribose promotes inositol trisphosphate (IP3)-evoked Ca2+ signals

In addition to its well established function in activating Ca(2+) release from the endoplasmic reticulum (ER) through ryanodine receptors (RyR), the second messenger cyclic ADP-ribose (cADPR) also accelerates the activity of SERCA pumps, which sequester Ca(2+) into the ER. Here, we demonstrate a potential physiological role for cADPR in modulating cellular Ca(2+) signals via changes in ER Ca(2+) store content, by imaging Ca(2+) liberation through inositol trisphosphate receptors (IP(3)R) in Xenopus oocytes, which lack RyR. Oocytes were injected with the non-metabolizable analog 3-deaza-cADPR, and cytosolic [Ca(2+)] was transiently elevated by applying voltage-clamp pulses to induce Ca(2+) influx through expressed plasmalemmal nicotinic channels. We observed a subsequent potentiation of global Ca(2+) signals evoked by strong photorelease of IP(3), and increased numbers of local Ca(2+) puffs evoked by weaker photorelease. These effects were not evident with cADPR alone or following cytosolic Ca(2+) elevation alone, indicating that they did not arise through direct actions of cADPR or Ca(2+) on the IP(3)R, but likely resulted from enhanced ER store filling. Moreover, the appearance of a new population of puffs with longer latencies, prolonged durations, and attenuated amplitudes suggests that luminal ER Ca(2+) may modulate IP(3)R function, in addition to simply determining the size of the available store and the electrochemical driving force for release.

Human delta opioid receptor biogenesis is regulated via interactions with SERCA2b and calnexin

Sarco(endo)plasmic reticulum calcium ATPase (SERCA)2b maintains the cellular Ca(2+) homeostasis by transferring Ca(2+) from the cytosol to the lumen of the endoplasmic reticulum (ER). Recently, SERCA2b has also been shown to be involved in the biosynthesis of secreted and membrane proteins via direct protein-protein interactions, involving components of the ER folding and quality-control machinery, as well as newly synthesized G protein-coupled receptors. Here we demonstrate that the human delta opioid receptor (hdeltaOR) exists in a ternary complex with SERCA2b and the ER molecular chaperone calnexin. The interaction between SERCA2b and hdeltaOR in vivo did not require calnexin as it was independent of the C-terminal calnexin-interacting domain of SERCA2b. However, the receptor was able to mediate co-immunoprecipitation of calnexin with the C-terminally truncated SERCA2b. The association of SERCA2b with hdeltaOR was regulated in vitro by Ca(2+) and ATP in a manner that was opposite to the calnexin-hdeltaOR interaction. Importantly, co-expression of the catalytically inactive SERCA2b(D351A) or calnexin binding-compromised SERCA2bDeltaC mutants with the receptor decreased the expression of mature receptors in a manner that did not directly relate to changes in the ER Ca(2+) concentration. We conclude that dynamic interactions among SERCA2b, calnexin and the hdeltaOR precursor orchestrate receptor biogenesis and are regulated by Ca(2+) and ATP. We further hypothesize that the primary role of SERCA2b in this process is to act as a Ca(2+) sensor in the vicinity of active translocons, integrating protein folding with local fluctuations of ER Ca(2+) levels.

Calibration of dynamic models of biological systems with KInfer

Methods for parameter estimation that are robust to experimental uncertainties and to stochastic and biological noise and that require a minimum of a priori input knowledge are of key importance in computational systems biology. The new method presented in this paper aims to ensure an inference model that deduces the rate constants of a system of biochemical reactions from experimentally measured time courses of reactants. This new method was applied to some challenging parameter estimation problems of nonlinear dynamic biological systems and was tested both on synthetic and real data. The synthetic case studies are the 12-state model of the SERCA pump and a model of a genetic network containing feedback loops of interaction between regulator and effector genes. The real case studies consist of a model of the reaction between the inhibitor kappaB kinase enzyme and its substrate in the signal transduction pathway of NF-kappaB, and a stiff model of the fermentation pathway of Lactococcus lactis.

Modelling sarcoplasmic reticulum calcium ATPase and its regulation in cardiac myocytes

When developing large-scale mathematical models of physiology, some reduction in complexity is necessarily required to maintain computational efficiency. A prime example of such an intricate cell is the cardiac myocyte. For the predictive power of the cardiomyocyte models, it is vital to accurately describe the calcium transport mechanisms, since they essentially link the electrical activation to contractility. The removal of calcium from the cytoplasm takes place mainly by the Na(+)/Ca(2+) exchanger, and the sarcoplasmic reticulum Ca(2+) ATPase (SERCA). In the present study, we review the properties of SERCA, its frequency-dependent and beta-adrenergic regulation, and the approaches of mathematical modelling that have been used to investigate its function. Furthermore, we present novel theoretical considerations that might prove useful for the elucidation of the role of SERCA in cardiac function, achieving a reduction in model complexity, but at the same time retaining the central aspects of its function. Our results indicate that to faithfully predict the physiological properties of SERCA, we should take into account the calcium-buffering effect and reversible function of the pump. This 'uncomplicated' modelling approach could be useful to other similar transport mechanisms as well.

SERCA pump activity is physiologically regulated by presenilin and regulates amyloid beta production

In addition to disrupting the regulated intramembraneous proteolysis of key substrates, mutations in the presenilins also alter calcium homeostasis, but the mechanism linking presenilins and calcium regulation is unresolved. At rest, cytosolic Ca²⁺ is maintained at low levels by pumping Ca²⁺ into stores in the endoplasmic reticulum (ER) via the sarco ER Ca²⁺-ATPase (SERCA) pumps. We show that SERCA activity is diminished in fibroblasts lacking both PS1 and PS2 genes, despite elevated SERCA2b steady-state levels, and we show that presenilins and SERCA physically interact. Enhancing presenilin levels in Xenopus laevis oocytes accelerates clearance of cytosolic Ca²⁺, whereas higher levels of SERCA2b phenocopy PS1 overexpression, accelerating Ca²⁺ clearance and exaggerating inositol 1,4,5-trisphosphate–mediated Ca²⁺ liberation. The critical role that SERCA2b plays in the pathogenesis of Alzheimer's disease is underscored by our findings that modulating SERCA activity alters amyloid β production. Our results point to a physiological role for the presenilins in Ca²⁺ signaling via regulation of the SERCA pump.

Intraluminal calcium as a primary regulator of endoplasmic reticulum function

The concentration of Ca2+ inside the lumen of endoplasmic reticulum (ER) regulates a vast array of spatiotemporally distinct cellular processes, from intracellular Ca2+ signals to intra-ER protein processing and cell death. This review summarises recent data on the mechanisms of luminal Ca2+-dependent regulation of Ca2+ release and uptake as well as ER regulation of cellular adaptive processes. In addition we discuss general biophysical properties of the ER membrane, as trans-endomembrane Ca2+ fluxes are subject to basic electrical forces, determined by factors such as the membrane potential of the ER and the ease with which Ca2+ fluxes are able to change this potential (i.e. the resistance of the ER membrane). Although these electrical forces undoubtedly play a fundamental role in shaping [Ca2+](ER) dynamics, at present there is very little direct experimental information about the biophysical properties of the ER membrane. Further studies of how intraluminal [Ca2+] is regulated, best carried out with direct measurements, are vital for understanding how Ca2+ orchestrates cell function. Direct monitoring of [Ca2+](ER) under conditions where the cytosolic [Ca2+] is known may also help to capture elusive biophysical information about the ER, such as the potential difference across the ER membrane.

Ca2+ signalling and Ca2+-activated ion channels in exocrine acinar cells

The development of the calcium signalling field, from its early beginnings some 40 years ago to the present, is described. Calcium signalling in exocrine gland acinar cells and the effects of neurotransmitter- or hormone-elicited rises in the cytosolic calcium ion concentration on ion channel gating are reviewed. The highly polarized arrangement of the organelle systems in living acinar cells is described as well as its importance for the physiologically relevant local and polarized calcium signalling events.

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(2+) pump activity increases the period and amplitude of oscillations; and 4), a threshold exists for the calcium concentration below which oscillations cease.

Calcium oscillations

Changes in cellular Ca2+ concentration control a wide range of physiological processes, from the subsecond release of synaptic neurotransmitters, to the regulation of gene expression over months or years. Ca2+ can also trigger cell death through both apoptosis and necrosis, and so the regulation of cellular Ca2+ concentration must be tightly controlled through the concerted action of pumps, channels and buffers that transport Ca2+ into and out of the cell cytoplasm. A hallmark of cellular Ca2+ signalling is its spatiotemporal complexity: stimulation of cells by a hormone or neurotransmitter leads to oscillations in cytoplasmic Ca2+ concentration that can vary markedly in time course, amplitude, frequency, and spatial range. In this chapter we review some of the biological roles of Ca2+, the experimental characterisation of complex dynamic changes in Ca2+ concentration, and attempts to explain this complexity using computational models. We consider the 'toolkit' of cellular proteins which influence Ca2+ concentrarion, describe mechanistic models of key elements of the toolkit, and fit these into the framework of whole cell models of Ca2+ oscillations and waves. Finally, we will touch on recent efforts to use stochastic modelling to elucidate elementary Ca2+ signal events, and how these may evolve into global signals.

A three-dimensional model of intercellular calcium signaling in epithelial cells

We have developed a fully three-dimensional (3D) model of calcium signaling in epithelial cells based on a set of reaction diffusion equations that are solved on a large-scale finite-element code in three dimensions. We have explicitly included the cellular compartments including the cell nucleus, cytoplasm, and gap junctions. The model allows for buffering of free Ca2+, calcium-induced calcium release, and the explicit inclusion of mobile buffers. To make quantitative comparisons to experimental results, we used fluorescence microscopy images of cells to generate an accurate mesh describing cell morphology. We found that Ca2+ wave propagation through the tissue is a function of both initial conditions used to start the wave and various geometrical parameters that affect propagation such as gap junction density and distribution, and the presence of nuclei. The exogenous dyes used in experimental imaging also affect wave propagation.

The endoplasmic reticulum: folding, calcium homeostasis, signaling, and redox control

The endoplasmic reticulum (ER) plays a major role in regulating synthesis, folding, and orderly transport of proteins. It is also essentially involved in various cellular signaling processes, primarily by its function as a dynamic Ca(2+) store. Compared to the cytosol, oxidizing conditions are found in the ER that allow oxidation of cysteine residues in nascent polypeptide chains to form intramolecular disulfide bonds. However, compounds and enzymes such as PDI that catalyze disulfide bonds become reduced and have to be reoxidized for further catalytic cycles. A number of enzymes, among them products of the ERO1 gene, appear to provide oxidizing equivalents, and oxygen appears to be the final oxidant in aerobic living organisms. Thus, protein oxidation in the ER is connected with generation of reactive oxygen species (ROS). Changes in the redox state and the presence of ROS also affect the Ca(2+) homeostasis by modulating the functionality of ER-based channels and buffering chaperones. In addition, a close relationship exists between oxidative stress and ER stress, which both may activate signaling events leading to a rebalance of folding capacity and folding demand or to cell death. Thus, redox homeostasis appears to be a prerequisite for proper functioning of the ER.

Ca2+ refilling of the endoplasmic reticulum is largely preserved albeit reduced Ca2+ entry in endothelial cells

In this study the relationship between the efficiency of endoplasmic reticulum (ER) Ca2+ refilling and the extent of Ca2+ entry was investigated in endothelial cells. ER and mitochondrial Ca2+ concentration were measured using genetically encoded Ca2+ sensors, while the amount of entering Ca2+ was controlled by varying either the extracellular Ca2+ or the electrical driving force for Ca2+ by changing the plasma membrane potential. In the absence of an agonist, ER Ca2+ replenishment was fully accomplished even if the Ca2+ concentration applied was reduced from 2 to 0.5mM. A similar strong efficiency of ER Ca2+ refilling was obtained under condition of plasma membrane depolarization. However, in the presence of histamine, ER Ca2+ refilling depended on mitochondrial Ca2+ transport and was more susceptible to membrane depolarization. Store-operated Ca2+ entry (SOCE), was strongly reduced under low Ca2+ and depolarizing conditions but increased if ER Ca2+ uptake was blocked or if ER Ca2+ was released continuously by IP(3). A correlation of the kinetics of ER Ca2+refilling with cytosolic Ca2+ signals revealed that termination of SOCE is a rapid event that is not delayed compared to ER refilling. Our data indicate that ER refilling occurs in priority to, and independently from the cytosolic Ca2+ elevation upon Ca2+ entry and that this important process is widely achieved even under conditions of diminished Ca2+entry.

Calcium oscillations and waves generated by multiple release mechanisms in pancreatic acinar cells

We explore the dynamic behavior of a model of calcium oscillations and wave propagation in the basal region of pancreatic acinar cells [Sneyd, J., et al., Biophys. J. 85: 1392-1405, 2003]. Since it is known that two principal calcium release pathways are involved, inositol trisphosphate receptors (IPR) and ryanodine receptors (RyR), we study how the model behavior depends on the density of each receptor type. Calcium oscillations can be mediated either by IPR or RyR. Continuous increases in either RyR or IPR density can lead to the appearance and disappearance of oscillations multiple times, and the two receptor types interact via their common effect on cytoplasmic calcium concentration and the subsequent effect on the total amount of calcium inside the cell. Increases in agonist concentration can stimulate oscillations via the RyR by increasing calcium influx. Using a two time-scale approach, we explain these complex behaviors by treating the total amount of cellular calcium as a slow parameter. Oscillations are controlled by the shape of the slow manifold and where it intersects the nullcline of the slow variable. When calcium diffusion is included, the existence of traveling waves in the model equation is strongly dependent on the interplay between the total amount of calcium in the cell and membrane transport, a feature that can be experimentally tested. Our results help us understand the behavior of a model that includes both receptors in comparison to the properties of each receptor type in isolation.

A buffering SERCA pump in models of calcium dynamics

Many cells use oscillations in calcium concentration to transmit messages. The oscillations largely result from an influx of calcium into the cytosol from the endoplasmic reticulum (ER), followed by an efflux of calcium from the cytosol back into the ER. The sarcoplasmic/endoplasmic reticulum calcium ATPase (SERCA) pump pumps calcium into the ER. It binds calcium on the cytosolic side and releases it on the ER side and in the delay between binding and release, calcium is buffered by the pump. We developed a model of a buffering SERCA pump and investigated whether including this in a model of calcium oscillations has any significant effects. We found that the oscillations produced when using the SERCA pump, which does not buffer calcium, have a larger amplitude and a slightly smaller period than when using the buffering SERCA pump. We show that the buffering SERCA pump shows adaptation to a stimulus, and we demonstrate that, by using a bidirectional SERCA pump, we are able to eliminate futile cycling of calcium between the cytosol and ER when the cell is at rest.

Biophysical re-equilibration of Ca2+ fluxes as a simple biologically plausible explanation for complex intracellular Ca2+ release patterns

Physiological regulation of Ca(2+) release from the endoplasmic reticulum (ER) is critical for cell function. Recent direct measurements of free [Ca(2+)] inside the ER ([Ca(2+)](ER)) revealed that [Ca(2+)](ER) itself is a key regulator of ER Ca(2+) handling. However, the role of this new regulatory process in generating various patterns of Ca(2+) release remains to be elucidated in detail. Here, we incorporate the recently quantified experimental correlations between [Ca(2+)](ER) and Ca(2+) movements across the ER membrane into a mathematical model ER Ca(2+) handling. The model reproduces basic experimental dynamics of [Ca(2+)](ER). Although this was not goal in model design, the model also exhibits mechanistically unclear experimental phenomena such as "quantal" Ca(2+) release, and "store charging" by increasing resting cytosolic [Ca(2+)]. While more complex explanations cannot be ruled out, on the basis of our data we propose that "quantal release" and "store charging" could be simple re-equilibration phenomena, predicted by the recently quantified biophysical dynamics of Ca(2+) movements across the ER membrane.