This Article Abstract Full Text (PDF) Alert me when this article is cited Alert me if a correction is posted Services Similar articles in this journal Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via HighWire Citing Articles via Google Scholar Google Scholar Articles by Baran, I. Search for Related Content PubMed PubMed Citation Articles by Baran, I.

Integrated Luminal and Cytosolic Aspects of the Calcium Release Control

Biophysics Department, Faculty of Medicine, "Carol Davila" University of Medicine and Pharmacology, 76241 Bucharest, Romania

Correspondence: Address reprint requests to Irina Baran, Bd. Eroii Sanitari 8, 76241 Bucharest, Romania. Tel.: 040-21-3125955; Fax: 040-21-3125955; E-mail: baran{at}theorl.theory.nipne.ro.

	ABSTRACT

TOP ABSTRACT INTRODUCTION INTRALUMINAL CALCIUM DYNAMICS MODEL FOR FAST GATING... SIMULATION OF CALCIUM SIGNALS... DISCUSSION APPENDIX 1: STEADY-STATE OPEN... APPENDIX 2. DWELL-TIME... APPENDIX 3: SIMULATION PROCEDURE REFERENCES

 
We propose here a unitary approach to the luminal and cytosolic control of calcium release. A minimal number of model elements that realistically describe different data sets are combined and adapted to correctly respond to various physiological constraints. We couple the kinetic properties of the inositol 1,4,5 trisphosphate receptor/calcium channel with the dynamics of Ca²⁺ and K⁺ in both the lumen and cytosol, and by using a detailed simulation approach, we propose that local (on a radial distance 2 µm) calcium oscillations in permeabilized cells are driven by the slow inactivation of channels organized in discrete clusters composed of between six and 15 channels. Moreover, the character of these oscillations is found to be extremely sensitive to K⁺, so that the cytosolic and luminal calcium variations are in or out of phase if the store at equilibrium has tens or hundreds µM Ca²⁺, respectively, depending on the K⁺ gradient across the reticulum membrane. Different patterns of calcium signals can be reproduced through variation of only a few parameters.

	INTRODUCTION

TOP ABSTRACT INTRODUCTION INTRALUMINAL CALCIUM DYNAMICS MODEL FOR FAST GATING... SIMULATION OF CALCIUM SIGNALS... DISCUSSION APPENDIX 1: STEADY-STATE OPEN... APPENDIX 2. DWELL-TIME... APPENDIX 3: SIMULATION PROCEDURE REFERENCES

 
Release of calcium ions from specialized compartments is a key process in numerous cellular events (Berridge, 1993; Clapham, 1995). The complex nature of intracellular calcium dynamics relies on the interplay between specifically induced calcium fluxes, such as release from the endo- or sarcoplasmic reticulum through the inositol 1,4,5 trisphosphate (IP₃) receptor or the ryanodine receptor (RyR), uptake by the endoplasmic reticulum (ER) through the smooth ER Ca²⁺ pump, or calcium influx from the outside through Ca²⁺ channels on the plasma membrane, among which store-operated calcium channels are thought to directly sense the ER calcium content. Exchange with mitochondria plays also a significant role in preventing prolonged increases in cytosolic Ca²⁺ (Hajnóczky et al., 1999).

Nevertheless, the coupled control on both channel gating and store filling plays the key role in shaping cytosolic calcium signals. Studies on the IP₃ receptor/Ca²⁺ channel in its native environment or reconstituted into lipid bilayers have revealed a biphasic effect of Ca²⁺ on the channel, given by ion binding to antagonist—activating/inhibitory—classes of sites, whereas IP₃'s binding properties and effects have received various interpretation. The global picture is even more complexed by parallel variations of the intraluminal source of calcium during release. The ER is known to have a high buffering capacity for Ca²⁺, but little is known about the exact dynamics and role, either direct or indirect, of store filling in calcium signaling. Recently an important feature of calcium release regulation has been characterized (Nguyen et al., 1998) through direct involvement of potassium ions toward mobilization of Ca²⁺ from the luminal matrix. Moreover, besides the dynamic Ca²⁺/K⁺ exchange during release, there exists the additional possibility of a direct effect of luminal calcium at the receptor level. Although this mechanism gained convincing evidence (Györke and Györke, 1998) in the case of the RyR, the corresponding effect on the IP₃ receptor (IP₃R) still awaits clarification (Taylor and Traynor, 1995).

The many facets of intracellular calcium signaling (Berridge et al., 2000) have proved to be rather difficult to incorporate into a consistent, unitary description. Different proposed models (see Dupont, 1999, for review) involving calcium release properties are able to explain quantal release or incremental detection when applied to permeabilized cell-like conditions, but cannot reproduce IP₃-induced oscillations unless applied to intact cell-like conditions, where a conservation rule is generally used for intracellular calcium. The few current quantitative models (De Young and Keizer, 1992; Bezprozvanny and Erlich, 1994; Kaftan et al., 1997; Swillens et al., 1998; Moraru et al., 1999; Swillens et al., 1999) on regulation of IP₃R activity have been mainly used to explain the Ca²⁺ and IP₃ dependence of the steady-state open probability, or to simulate Ca²⁺ release at discrete sites, at constant luminal or cytosolic Ca²⁺. Moreover, these models are based on data obtained with channels incorporated into reconstituted bilayers, which are quantitatively largely different from data on the outer nuclear membrane.

The present study explores a new dimension of the calcium release dynamics and intends to answer two main questions: 1), to what extent can calcium signals depend on K⁺, and 2), how can in- or out-of-phase local calcium oscillations be explained in permeabilized cells? To this end we construct a basis for quantifying the dynamical relation between store and cytosolic calcium by coupling it to the IP₃R/Ca²⁺ channel kinetics and to the control of the luminal calcium content by potassium. The model we propose here is based on the derivation of basic parameters that define both Ca²⁺/K⁺ exchange inside the store and gating properties of the ER calcium channel in nuclear membranes. With these parameters in hand, we simulate Ca²⁺ release through a cluster of channels in cylindrical geometry, taking into consideration ion diffusion and binding to endo- and exogenous buffers within the lumen and cytosol, in conditions that reproduce the main features of the physiologic subcellular environment. The model accurately describes a series of disparate experimental data that have not been theoretically considered to the moment and gives a more complete and detailed simulation approach to the calcium release process than previously done. Various approximations currently used (such as buffering approximations, spherical symmetry, constant luminal Ca²⁺, constant channel current, and neglect of intraluminal ion diffusion) are replaced here by direct approach to the underlying phenomena.

We discuss two possible scenarios for the K⁺ role in calcium release and show that the previously proposed (Nguyen et al., 1998) K⁺-dependent release cycle leads to parallel, not antiparallel, local calcium oscillations. Interestingly, we find that local calcium oscillations observed on a seconds scale in permeabilized cells might be driven by the slow inactivation of channels organized in discrete clusters composed of between 6 and 15 channels, and out-of-phase calcium variations appear at rather constant luminal K⁺. Contrary to the expectation, active implication of K⁺ actually appears as a key factor in inducing parallel oscillations.

	INTRALUMINAL CALCIUM DYNAMICS

TOP ABSTRACT INTRODUCTION INTRALUMINAL CALCIUM DYNAMICS MODEL FOR FAST GATING... SIMULATION OF CALCIUM SIGNALS... DISCUSSION APPENDIX 1: STEADY-STATE OPEN... APPENDIX 2. DWELL-TIME... APPENDIX 3: SIMULATION PROCEDURE REFERENCES

 
In the following we consider the observation (Nguyen et al., 1998) that the glycoprotein matrix inside calcium stores functions as a continuous Ca²⁺/K⁺ exchanger. We have used a simple model for describing this commerce, whereby Ca²⁺ and K⁺ compete in binding to a common matrix domain M, as shown in Fig. 1. Ion binding kinetics is assumed to be first-order for calcium and h_K-order for potassium. All over the paper h denotes Hill coefficients; k, rate constants, K, apparent dissociation constants, and [X], concentration of the X specie. Then the equations governing the dynamical changes in the luminal concentrations read:

(1)

(2)

(3)

(4)

where [M] + [MCa] + [MK] = M_max is the total common Ca²⁺/K⁺ binding domains concentration.

View larger version (6K): [in this window] [in a new window]   FIGURE 1  Reactions assumed in Ca2+ and K+ competing binding to the same luminal matrix domain M.

(5)

where is the total luminal calcium (free + bound) content of resting (unstimulated) permeabilized cells in 140 mM K⁺ incubating medium. Here the apparent dissociation constants are defined as and

The best fit (shown in Fig. 2) of the data (Nguyen et al., 1998) to Eq. 5 was obtained with the following parameter values:

1). For the endoplasmic reticulum: and and

2). For isolated secretory vesicles: , and

View larger version (17K): [in this window] [in a new window]   FIGURE 2  Derivation of intraluminal Ca2+-K+ exchange parameters. Fit of data (Nguyen et al., 1998) to Eq. 5 is obtained by least-square method. Either plot ( filled circles, or open circles) is drawn by using the corresponding parameters specified in the text.

Regulation of the steady free calcium level (100 µM–300 mM) inside the ER is thought to be achieved through binding of Ca²⁺ ions to luminal glycoproteins as calsequestrin or calreticulin (Nguyen et al., 1998). Purified calsequestrin (with 109 acidic amino acids) binds 40 mol Ca²⁺/mol of calsequestrin with K_D = 0.4–1 mM at 150 mM KCl (Mitchell et al., 1988), and 1 mol of calreticulin binds 20 mol Ca²⁺. A quick inspection of our results shows that a mean buffer protein concentration of 360 mM/40 10 mM would rely inside the ER of cilliated cells. Biochemical analysis should reveal the value of M_max in the absence of K⁺. In favor of this, a decrease in the K⁺ concentration from 150 mM to 20 mM was shown to induce a 6- to 10-fold increase in Ca²⁺ binding, meaning that a single protein molecule can bind 400 Ca²⁺ ions at 20 mM K⁺, thus supplying a large number of Ca²⁺ binding sites of the order 400 x 10 mM = 4 M. The 6- to 10-fold increase in calcium binding also agrees with our result

In the mucin granules of goblet cells, the resting Ca²⁺ level is 20 µM (Nguyen et al., 1998). The same calculation as above yields a quantity of 500 nM luminal binding protein. The different K_Ca and h_K obtained with the two organelles suggest differences in intraluminal conditions (pH, chemical composition, or homogeneity degree) between different calcium stores. Nevertheless further diversity may be due to the cell type (cardiac/cilliated/goblet) or environmental (in vitro/in vivo) conditions (pH, ionic strength) of analysis. This, and the small amount of data, could explain why such a large capacity of intraluminal proteins to bind Ca²⁺ has not been detected yet. Quantitative measure of the K⁺ binding capacity of ER buffers would be extremely valuable in determining the exact mechanisms involved in the luminal Ca²⁺/K⁺ ion exchange.

Two other points are to be stressed here. One, it is possible that at very low calcium concentration, the ER membrane gets very leaky, similarly to the plasma membrane (Hille, 1992), so we preferred removing the low calcium data, which otherwise alter the fit. Second, the above parameters constrain a minimal free resting which seems to disagree with the 25–35 µM values observed (Nguyen et al., 1998). We distinguish between two possible explanations for this apparent discrepancy. First, the data on steady-state Ca²⁺/K⁺ ion exchange in both the ER and secretory vesicles are obtained with pH 5.5, whereas pH_ER = 7.1 has been determined in vivo in HeLa cells (Kim et al., 1998). The most reasonable assumption is that a higher, physiological pH might change the Ca²⁺ and K⁺ affinities for the luminal matrix. Indeed, a steady can be readily obtained with calculation based on the assumption that and are reduced by the increased pH through a factor f_dis > 1 (discussed later on). As a second possibility, we assume f_dis = 1 (the pH has no effect on and ) and consider calcium outflow at rest. Inasmuch as a steady state is reached in the Ca²⁺/K⁺ exchange experiments, we infer that the corresponding flux is flowing through the IP₃Rs; otherwise the lumen eventually gets empty. This scenario works if IP₃Rs have nonzero open probability in resting cells, which implies a basal, endogenous [IP₃] > 0. In the absence of any other data that could clarify these aspects, we shall consider both alternatives in our simulations on calcium release.

	MODEL FOR FAST GATING OF THE IP3R/Ca2+ CHANNEL

TOP ABSTRACT INTRODUCTION INTRALUMINAL CALCIUM DYNAMICS MODEL FOR FAST GATING... SIMULATION OF CALCIUM SIGNALS... DISCUSSION APPENDIX 1: STEADY-STATE OPEN... APPENDIX 2. DWELL-TIME... APPENDIX 3: SIMULATION PROCEDURE REFERENCES

 
IP₃ binds to a specific protein complex—the IP₃R/Ca²⁺ channel, located on intracellular Ca²⁺ stores where it promotes and regulates release of Ca²⁺ into the cytosol. Each IP₃R is a tetramer having a well-defined IP₃ binding site per subunit near its NH₂-terminus (Marshall and Taylor, 1993; Taylor and Traynor, 1995). Two additional sites, with low and high Ca²⁺ affinity have been described (Kaftan et al., 1997; Moraru et al., 1999), but their location is uncertain. Among all IP₃ binding sites, the 1-nM site is the least (1%) frequent.

Biochemical studies indicated that Ca²⁺ binds to the same receptor at up to eight sites per monomer (Mignery et al., 1992; Sienaert et al., 1996, ¹⁹⁹⁷), with seven sites residing on the cytosolic region of the receptor (Sienaert et al., 1997). The Ca²⁺ affinity has been determined for only one of these sites (K_d = 0.8 µM, Hill coefficient 1). The luminal Ca²⁺ binding site (Sienaert et al., 1996) could modulate the channel activity as in the case of the homologous Ca²⁺ release channel of skeletal and cardiac muscle, the RyR, where the existence of such site has clearly been evidenced (Györke and Györke, 1998). The IP₃R is currently considered to be regulated by two functional classes of cytosolic calcium binding sites, of which one is activating (we denote it S_a) and one inhibitory (S_i). With the channel in its ER membrane environment, the open probability (P_o) data in the presence of 10 µM IP₃ are well described by a biphasic Hill formula (see Eq. 1 in Mak et al., 1998), whereby Ca²⁺ binds to S_a with activation constant K_a = 210 ± 20 nM and Hill coefficient h_a = 1.9 ± 0.3, and to S_i with K_i = 54 ± 3 µM and Hill coefficient h_a = 3.9 ± 0.7. As for bilayer experiments, with 2 µM IP₃ applied, the channel's activity is maximal in the presence of 0.25 µM cytosolic Ca²⁺ (Bezprozvanny et al., 1991; Kaftan et al., 1997) and completely inhibited at 5 µM Ca²⁺, showing both positive and negative feedback of Ca²⁺ in a bell-shaped P_o([Ca²⁺]) curve, whereas with 180 µM IP₃ there is maximal activity at 1 µM Ca²⁺, with a shallow decline at Ca²⁺ up to 40 µM. The predictions here (Bezprozvanny et al., 1991; Kaftan et al., 1997) for the K_d's are in the range 0.03–0.2 µM for the activating site and 0.2–1 µM for the inhibitory site, with two or three Ca²⁺ ions required to bind to the channel for each class of sites.

The regulatory domain of the receptor also includes two–three putative adenine nucleotide-binding sites (Marshall and Taylor, 1993; Taylor and Traynor, 1995). ATP bound to these sites stimulates the receptor, but the way it does so has to be clarified. Direct and systematic investigations (Mak et al., 1999) of the mechanism involved in this regulation showed that ATP allosterically regulates the Ca²⁺ sensitivity of the Ca²⁺ activating sites with a functional dissociation constant K_ATP = 0.27 mM and Hill coefficient of 1, meaning that only one ATP molecule binds to the receptor to stimulate it. The same study suggest that ATP is also involved in the regulation of the Ca²⁺ affinity of the inhibitory sites, but the data are not sufficient to draw a firm conclusion.

The type I IP₃ receptor has been shown to have a high open probability, (Mak et al., 1998, 1999), in its native membrane environment, at the physiologic IP₃ (1 µM) and cytosolic Ca²⁺ concentration ([Ca²⁺]_cyt 0.1 µM), although displaying a greatly reduced activity when incorporated into lipid bilayers, (Kaftan et al.,1997), respectively. In either case the P_o dependence on calcium (Bezprozvanny et al., 1991) is bell-shaped, as has been shown also for IP₃-induced calcium-45 release in brain synaptosoms (Finch et al., 1991). The data used and analyzed so far in various theoretical models are those obtained with the channel incorporated in artificial membranes. To describe its activity, 8- (De Young and Keizer, 1992), 16- (Kaftan et al., 1997), 18- (Swillens et al., 1999) and recently more complex, 125-state (Moraru et al., 1999) models have been proposed, out of which only the last one describes features of both P_o and open-time distributions.

Here we use and analyze data exclusively obtained with channels in their native environment and construct a simple operational model to accurately describe not only P_o but also _o (open time) and _c (close time), which is quite a difficult task. Our first steps were to consider first a two-site (four-state) and then three-site (eight-state) model. However, the quantitative results were unsatisfactory, for the models in no way could simultaneously fit the three data sets (Mak et al., 1998). We then added a fourth regulatory site and obtained good agreement to the data.

So, our minimal model treats the tetramer receptor as a complex unit, equivalent to a 16-state system at thermal equilibrium. Kinetics of calcium binding at either site are independent of each other, but dependent on binding to the receptor of a specific factor, f_a or f_i, respectively, as shown in Fig. 3 A. Hence, we are introducing two other regulatory sites, S_fa and S_fi. Because the open dwell time histograms of this channel generally show two distinct components, we assume there are two open states, characterized by bound S_a and both S_i and S_fi unbound.

View larger version (20K): [in this window] [in a new window]   FIGURE 3  Four-state channel representation and transitions within basic regulatory modules of IP3R fast and slow gating, as assumed in the model. In A, x = a or i and X = A or I characterize the activating and inhibitory module, respectively.

(6)

where (x = a, i) is the equilibrium constant of the transition q_(f)x = 0 to q_(f)x = 1. Here q_(f)x is the site S_(f)x occupancy (0 or 1) and

Based on estimation (see Appendix 2) of errors introduced by events undetected during recording of single-channel activity, the data obtained with the Xenopus oocyte type I IP₃R/calcium channel in patches of the outer nuclear membrane are well (least-square method) fitted with the parameter values given in Table 1, with arbitrary and Mean values of the parameters are obtained by best fit with _d = 0.3 ms, the approximate experimental value (Mak et al., 1998; Mak, personal communication). Parameter deviations from mean values are determined such that theoretical curves obtained with 0.2 ms _d 0.4 ms remain within experimental error bars. The model simultaneously comes in good agreement (shown in Figs. 4 and 5) with the measured open/close time and P_o for the specific conditions of the electrophysiological experiments (Mak et al., 1998), namely: 1.5-µM trans Ca²⁺, 10-µM IP₃, and 500-µM ATP on the cytosolic IP₃R side. We have to mention that we obtained a good agreement with the triplet data, i.e., and only when a high affinity calcium binding site (K_a = 12 nM) was afforded for. The 500-µM ATP data could be fitted as well with K_a0 = 0.3 µM and h_a0 = 2, with essentially the same results.

View larger version (33K): [in this window] [in a new window]   FIGURE 4  [Ca2+]cyt dependence of the mean open- and closed-channel duration at 10 µM IP3 and 500 µM free ATP. Data are from Mak et al. (1998) and Mak (personal communication). Continuous lines are drawn according to Eqs. 17–19 with parameter values given in the text. The curves in each panel correspond respectively, from the lowest to the upper, to different values of d: 0, 0.2, 0.3, and 0.4 ms in the lower panel and 0, 0.2, and 0.4 ms in the upper panel.

View larger version (25K): [in this window] [in a new window]   FIGURE 5  Open probability of IP3R at 10 µM IP3 and different cytosolic ATP concentrations: 0–12 µM (small squares), 500 µM (filled circles), 4.8 mM (large squares), and 9.5 mM (open circles). Data are from Mak et al. (1998, 1999). Curves represent calculated Actual Po differs by <5% (not shown). Dashed line corresponds to predicted Po at 1 mM [Ca2+]lum in the absence of ATP, assuming

Our calculation (not shown) based on the procedure described (Appendices 1 and 2) proves that simpler gating models (with two or three regulatory sites) cannot adequately fit into the experimental data. We want to stress that the two gating modules discussed here have been developed from the two standard calcium sites. The hypothesis is made that kinetics of Ca²⁺ binding at each of these two sites are regulated by a distinct factor. At first we had devoted a lot of work to the identification of the ligand f_a as IP₃ but this resulted (not shown) in a marked dependence of the activation domain of the curve on the concentration of IP₃, contrary to the data (Mak et al., 1998). Another good candidate, namely ATP, for the same factor f_a has been suggested before (Mak et al., 1999) and we have proved here that interpreting ATP as the factor implicated in the activation module regulation very well can explain the related features of channel activity.

As for the inhibiting module, evidence is less conclusive. A good choice initially seemed to be IP₃, as in the model of De Young and Keizer (1992), with the IP₃ binding transition adapted to the present situation, but the results (not shown) came in disagreement with the IP₃ effects on P_o (Mak et al., 1998), mainly because P_o,max greatly depended on IP₃ concentration, and the activation domain of the curve was shifted to the right by decreasing [IP₃], which is not the case. Neither could ATP provide good explanation for this behavior (from our investigations, not shown). Consequently, we hypothesize here that the factor f_i might be luminal Ca²⁺. In this way the model can explain key features of channel kinetical behavior, such as inhibitory effects of high (0.6–1 mM) [Ca²⁺]_lum in the absence of ATP (Thrower et al., 2000); apparent invariability of P_o at low [Ca²⁺]_lum (Mak et al., 1998); decreased open channel duration observed in the presence of calcium in the trans chamber (Bezprozvanny and Erlich, 1994); or the shallow decrease (Thrower et al., 2000) of the longest open time (here _o1) with increasing [Ca²⁺]_lum. The implicated mechanism considered here is closing of the IP₃R by luminal Ca²⁺, whereas for the RyR the opposite effect, namely opening of the luminal gate by calcium, can account for the data obtained (Györke and Györke, 1998). More precisely, for the IP₃R, the transition from the (0,0) to the (0,1) state of the inhibition module is assumed to be driven by calcium binding to the luminal site, thus closing the channel (if open). It is quite difficult to make an exact correspondence between data obtained in different experimental conditions. Channel inactivation by 0.6–1 mM [Ca²⁺]_lum in absence of ATP agrees well with (our model predicts—not shown—that this inactivation is only apparent, with all open channel events being too short for detection with 0.5-ms time resolution (Thrower et al., 2000); in addition, removal of 1 mM [Ca²⁺]_lum inhibition by ATP would agree with an increase of with ATP, which pushes _o above the detection limit). Instead, for specific experimental conditions (Mak and Foskett, 1997), the value gives and values, close to those obtained by experiment ( and <5 ms, respectively), but implies a lower threshold for the inhibitory [Ca²⁺]_lum (50 µM). The calculated subunitary Hill coefficient could result from negative cooperativity in calcium binding at the luminal site. It yields a dissociation constant for and for

The analysis here does not uniquely declare as the factor f_i, but provides a series of arguments in favor of this equivalence. In further simulations we shall consider three possible alternatives on this matter (see below). As a final point here, it is worth mentioning that the IP₃ dependence of the curve can be well reproduced (not shown) by considering the proposed allosteric regulation of the calcium affinity for the S_i site by IP₃ binding to its 50-nM site with a Hill coefficient of 4 (see Eq. 2, Mak et al., 1998), where K_max = 54 µM is replaced by the 45 µM obtained here. This dependence will be used in all our simulations (see below).

	SIMULATION OF CALCIUM SIGNALS IN PERMEABILIZED CELLS

TOP ABSTRACT INTRODUCTION INTRALUMINAL CALCIUM DYNAMICS MODEL FOR FAST GATING... SIMULATION OF CALCIUM SIGNALS... DISCUSSION APPENDIX 1: STEADY-STATE OPEN... APPENDIX 2. DWELL-TIME... APPENDIX 3: SIMULATION PROCEDURE REFERENCES

 
Because the calcium release apparatus seems to behave differently in various cell types, generating different patterns of Ca²⁺ signals, a better understanding of the molecular and cellular mechanisms involved is of great interest. Given the large differences between the channel's properties considered here as compared to previous theoretical studies, we have been looking for defining their effects and determining how the 16-state model could integrate into an operational model of calcium release. We expected largely different features of release kinetics because the channel is inhibited now at higher [Ca²⁺]_cyt (100 µM vs. 1 µM). Furthermore, calcium release regulation by K⁺ appears to be crucial in the oscillatory Ca²⁺ response to agonist stimulation, inasmuch as the calcium signal disappears (Nguyen et al., 1998) in the presence of apamin, a blocker of the small conductance calcium-activated potassium (ASK) channels (see below).

To quantitatively study the effects of potassium ions in calcium release regulation, we coupled the IP₃ gating model analyzed above to the dynamics of calcium and potassium in both the lumen and cytosol. Numerical and biophysical aspects of our simulation procedure are detailed in the Appendix 3. The model was tested first on a documented scenario (Nguyen et al., 1998) and improved to come in agreement with two basic processes: 1); calcium oscillations are triggered locally in permeabilized cells; and 2), K⁺ is implicated in the (in- or out-of-phase) character of these oscillations. It is worth mentioning that a rise of such oscillations in permeabilized cells has not been theoretically explained yet. We have examined the system's behavior under various physiological constraints and inspected the sensitivity of the calcium signals to several parameters. To remain in agreement with the experimental evidence, further improvement of the model has emerged and critical ranges of the free parameters in the model have been determined.

Slow kinetics of calcium release
Simulations based on the 16-state model
At first we strictly used the scenario proposed (Nguyen et al., 1998) and considered (Scenario I) that the equilibrium resting state has low (25–35 µM) [Ca²⁺]_ER. Accordingly, we had to consider a potassium extrusion term, which could be due to the K⁺-H⁺ exchanger across the reticulum membrane (Nguyen et al., 1998), to restore the chemical equilibrium (see Eq. 5) inside the store. As mentioned before, we were led to lower the dissociation constants and by a factor f_dis. After IP₃ stimulation, the IP₃R/calcium channels open and the resulting increase in [Ca²⁺]_cyt activates the ASK channels on the ER membrane (Nguyen et al., 1998). Consequently, K⁺ inflow liberates Ca²⁺ from the luminal matrix. In all our simulations, the system first displayed a pronounced but transient increase in local [Ca²⁺]_cyt (3 µM) and ran thereafter into a quasi steady state with high [Ca²⁺]_ER (between 2 and 5 mM) and slightly elevated [Ca²⁺]_cyt (0.6 µM). Calcium release is silenced by the inhibitory effect of the high microdomain (Naraghi and Neher, 1997) calcium concentration built up by a current as high as 0.4–1 pA flowing through the channel from the submembrane [Ca²⁺]_ER of 2–5 mM.

We then undertook the possibility that gating kinetics is modified in vivo (Mak and Foskett, 1997), and assumed 100 times slower channel gating. An example, obtained with N_Ca = 2, is shown in Fig. 6 A. Increasing N_Ca produces a similar first spike, followed by more frequent, smaller spikes. However, in all the cases we have investigated (for large variation of the free parameters (not shown)), the resulting variations in intra- and extraluminal free calcium concentrations were not out of phase as expected, but always parallel. In addition, these variations are irregular, not oscillatory-type.

View larger version (50K): [in this window] [in a new window]   FIGURE 6  Simulated responses after stimulation of permeabilized cells with 3 µM IP3 at t = 0 in A and B and at moments indicated by arrows in C, D, and E. In C and D, a basal [IP3]b = 30 nM and 100 µM EGTA are considered; otherwise EGTA is present with 40 µM.

A slow gating component is required
From the results obtained so far, we conclude that the simple model for IP₃R gating is not sufficient to explain the observed patterns of calcium signals. So, we suspected that a slow kinetic component of the Xenopus oocyte IP₃R has not manifested in vitro (Mak et al., 1998), or was too slow to be detected (see the case of the RyR that adapts 100 times slower in vitro than in vivo (Sachs et al., 1995). Recordings of the type I IP₃R activity in its native environment show indeed that an additional, slow kinetic component is active in vivo. However, analysis of inactivation at the receptor level is rendered extremely difficult by its slow kinetics, which reduces the number of events below a reasonable limit appropriate for a standard statistical analysis. When analyzed in the "on-nucleus" configuration (Mak and Foskett, 1997), the channel preserved its fast kinetics features (discussed above), but displayed an additional, slow kinetic component.

All these, along with observations that Ca²⁺ and IP₃ inactivate IP₃R and IP₃-mediated calcium release on a seconds scale (Parker and Ivorra, 1990; Finch et al., 1991; Hajnóczky and Thomas, 1994; Oancea and Meyer, 1996), have determined us to introduce a Ca²⁺- and IP₃-dependent slow regulatory module. The reaction scheme, shown in Fig. 3 B, was first proposed to explain adaptive behavior and can account for channel inactivation after ligand binding (Sachs et al., 1995). The current active substates in Fig. 3 B are (1,0) and (0,1), so the channel has now four open states, which may correspond to the four different subconductance states observed experimentally (Bezprozvanny et al., 1991, Mak and Foskett, 1997; Thrower et al., 2000). We have chosen the functional dependence consistent with the interpretation (Hajnóczky and Thomas, 1994) that IP₃ may aid Ca²⁺ to access its inhibitory site (here by fastening Ca²⁺-binding kinetics), and with three other observations: 1), the initial release rate increases (Hajnóczky and Thomas, 1994) with IP₃ in the range 1.5–7.5 µM, so we assume that the IP₃ site involved is the low affinity one, i.e., 2), the channel recovers from inactivation after IP₃ removal or chelation of calcium in the presence of saturating IP₃ (Hajnóczky and Thomas, 1994); and 3), half-maximal channel inactivation appears at 0.2 or 0.6 µM [Ca²⁺]_cyt (Finch et al., 1991; Hajnóczky and Thomas, 1994); we use then K_I = 0.2 µM and h_I = 2, as proposed for channels in reconstituted bilayers (Bezprozvanny et al., 1991). The thermodynamical condition was used as before to derive Other parameters were chosen to yield two inactivation components that differ by at least one order (Hajnóczky and Thomas, 1994) and to match as possible the results of Mak and Foskett (1997). For the conditions therein (10 µM IP₃, 1 mM ATP, 200 nM [Ca²⁺]_cyt, and 200 nM [Ca²⁺]_lum), event distributions generated by the full 64-state model satisfactorily agree (not shown) with both P_o and the mean open and closed time histograms for two sets of kinetic parameters, as specified in Table 3. Analytical and simulated activity of the 64-state, but not 16-state, channel shows inactivation (not shown) of biexponential time course with rate constants of the order of those determined by experiment (Hajnóczky and Thomas, 1994).

Because the present data on the dependence of IP₃R gating on luminal calcium are not conclusive, we relaxed the assumption (see above) that by considering either a weak dependence or independence (i.e., ) of the inhibitory module on Both choices were not critical to the results. On the contrary, in several simulations (not shown) with the same conditions but the channels' P_o decreases considerably at so that the system remained in a quasi steady state with high and low We have not followed this line of investigation, inasmuch as it would be only speculative at the moment. In numerous simulations we have varied the EGTA concentration as specified in Table 2, but its effect on the slow timescale was not critical, as it mainly operates on the signal's amplitude, which remains in the noninhibitory [Ca²⁺]_cyt domain (not shown).

Efficiency of interchannel communication was shown (Swillens et al., 1999) to decrease with the distance between receptors (d), and a compact (d = 12 nm) cluster of 20–30 channels could explain the currently observed distributions of puffs and blips. Due to the high complexity of the system, our spatial grid was always d 35 nm. As an alternative (Swillens et al., 1999) we also performed several simulations where a highly compact cluster of 25 channels was comprised in a unique volume element, with 50-nm gridding. The open probabilities were independently simulated for all channels, which are facing now the same microdomain [Ca²⁺]_cyt. The results obtained with this procedure are not qualitatively different from those shown in Fig. 6 D.

Extremely high density of IP₃Rs corroborated by rapid Ca²⁺ mobilization inside the ER can produce sustained, high global [Ca²⁺]_cyt after stimulation, without requiring external calcium influx as currently invoked for explaining a high Ca²⁺ level in some cells after stimulation (Dupont, 1999). For example, a formation of six clusters, 50-nm interspaced, each consisting of 20, highly packed IP₃Rs, yields (not shown) an elevated global [Ca²⁺]_cyt fluctuating between 0.83 and 0.94 µM, whereas the radial average varies between 29–38 µM and [Ca²⁺]_ER remains at 100–113 µM, for

Use of the full model according to Scenario I generates parallel variations in local intra- and extraluminal calcium (an example is shown in Fig. 6 B, obtained with N_Ca = 9 and the Set 2 of inactivation parameters). Differences between simulations using either of the two sets of inactivation parameters were not crucial for the dynamical behavior of the system. We also have checked that with K_I = 0.66 or 0.8 µM (Finch et al., 1991; Moraru et al., 1999) and h_I = 1 or 2, oscillatory evolution can as well be reproduced within both scenarios, in conditions mentioned above.

From a careful inspection by simulation (not shown) we conclude that variation of the standard parameters within ranges given in Table 2 do not result in qualitative changes in the shape of calcium signals on the seconds timescale. Instead, the completely free parameters left (namely N_Ca, N_K, f_dis, and ) appear to be crucial in their modulation. Generally we used N_Ca = 1–25 (Swillens et al., 1999) and N_K = 50–250. We have examined the system's response to large variations in these parameters (not shown). Depending on the scenario, Table 4 indicates parameter values constrained by the requirement to obtain realistic results. For example, within Scenario I, the values therein give fast Ca²⁺ mobilization and results relevant to the situation. Lower values greatly enlarge the dynamic profiles (not shown). In summary, the system's behavior was found to be extremely sensitive to three factors:

1. The rapidity of Ca²⁺ and K⁺ ions in dissociating from the luminal glycoprotein matrix.
   
2. The number of Ca²⁺ channels in the cluster.
   
3. The difference between the luminal and cytosolic K⁺ concentration.
   

View larger version (31K): [in this window] [in a new window]   FIGURE 7  Subseconds calcium release at a single channel site, according to Scenario II. During the 120 s lapsed after 3 µM IP3 stimulation, the channel displays two openings and the corresponding effects are visualized as evolutions of the ionic current (A), local cytosolic (B), or luminal (C) Ca2+ concentration, as well as the global cytosolic calcium signal and wave propagation within the cytosol (D). The buffering level is given by 200 µM endogenous fix buffer and 40 µM EGTA.

View larger version (40K): [in this window] [in a new window]   FIGURE 8  Fast calcium release dynamics simulated with nine-channel release site. Other conditions as in Fig. 7.

View larger version (40K): [in this window] [in a new window]   FIGURE 9  Fast calcium release at six-channel release site in conditions specified by Scenario I. The current through both Ca2+ and K+ channel in the central element on the ER membrane are shown in the upper panels (ICa and IK, respectively). Other notations as in Fig. 7. Over the simulated 250-ms period after IP3 stimulation, calcium is released through the IP3Rs cluster during the first 26 ms only. The stationary buffer level is 120 µM.

	DISCUSSION

TOP ABSTRACT INTRODUCTION INTRALUMINAL CALCIUM DYNAMICS MODEL FOR FAST GATING... SIMULATION OF CALCIUM SIGNALS... DISCUSSION APPENDIX 1: STEADY-STATE OPEN... APPENDIX 2. DWELL-TIME... APPENDIX 3: SIMULATION PROCEDURE REFERENCES

Initially, most models (De Young and Keizer, 1992; Atri et al., 1993; Bezprozvanny and Erlich, 1994) considered three sites responsible for channel activity regulation: one site for IP₃ and two sites for Ca²⁺, one activating and one inhibitory. Within this scheme, the IP₃ effect on P_o could partially be described by the De Young and Keizer model only. More recently, a low affinity IP₃ binding site was included (Kaftan et al., 1997) to explain how the P_o([Ca²⁺]) curve is rightward-shifted by increasing IP₃ and inhibition by up to 50 µM Ca²⁺ relieved. Later, this model has been reformulated (Moraru et al., 1999) to give better fits to the P_o and _o data in dependence on both Ca²⁺ and IP₃. To this end, it was assumed that each receptor monomer contains one Ca²⁺ regulatory site and two IP₃ sites, of medium and low affinity, respectively. Another combination has been used in agreement to the P_o data at 2 µM IP₃ (Bezprozvanny et al., 1991), by considering one IP₃ site, two Ca²⁺ activating sites, and two Ca²⁺ desensitizing sites, but the model (Swillens et al., 1998) is based on a limited set of data. Essentially, the various outcomes of these models point out to the need for continuous refinement in modeling the complex IP₃ regulation as more data are provided and the specific regulatory binding sites identified.

In addition, the different behavior of the channel in different membranes makes the results more confusing. It is possible that reconstitution of the receptor into lipid bilayers alters the protein-protein and protein-lipid interactions and basically reduces by one order the magnitude of P_o. It would be interesting to know whether P_o at [Ca²⁺]_cyt 50 µM show inactivation at high IP₃ levels in bilayer experiments (see Fig. 6 in Moraru et al., 1999). If this would be the case, the dependence of P_o on Ca²⁺ and IP₃ could be explained by the same mechanism as proposed for the channel studied in the ER membrane. Here the data (Mak et al., 1998) indicate that two antagonist classes of Ca²⁺ binding sites exist and that Ca²⁺ binds to two monomers to open and to four monomers to close the channel. Moreover, IP₃ binding to the medium affinity site only seems to be implicated in regulation of channel activity, and it does so by attenuating the inhibitory effect of Ca²⁺. Channel activity resumes to fast kinetics in excised membrane patches but spreads on a slow timescale in integral nuclear membrane, clearly showing loss of the slow component in the former case.

To describe the fast gating component of the ER Ca²⁺ channel, we assume it to be controlled by Ca²⁺ binding to an activating site, with a calculated Hill coefficient by Ca²⁺ binding to an inhibitory site, with a calculated Hill coefficient and by a hypothetical ligand (f_i) to account for the subunitary P_o,max at the plateau domain of the P_o(Ca²⁺) curve (Mak et al., 1998). We show that our model supports the ATP effect on channel activation by allosteric regulation with K_d = 215 µM and h = 1.3, similar to those deduced from a nonstandard Hill formula (Mak et al., 1999). It would be worth seeing whether the ATP binding site could have such low affinity for free ATP in the absence of Ca²⁺ as calculated here Among the possible candidates as f_i ligands, only comes in agreement with various data as analyzed with our model, but certainly more data are needed to define the exact mechanism of regulation.

Based on the results obtained here, we propose a number of possible mechanisms involved in calcium release regulation that fairly well describes many experimental data and observations, as follows:

1. The ER of cilliated cells is a huge reservoir of 4 M Ca²⁺ binding sites, most of which are masked by K⁺, and a Ca²⁺ ion must compete with six K⁺ ions inside the ER and with 11 K⁺ ions inside the mucin granules of goblet cells.
   
2. Ca²⁺ has a high (K_dis = 12 nM) or a 50 times lower affinity for its activating site on the cytosolic side of the IP₃R, depending on whether ATP is bound or not.
   
3. Luminal Ca²⁺ might inhibit IP₃R activity with a low, possibly subunitary Hill coefficient.
   
4. Both comparison of data on IP₃R in excised versus intact nuclear membrane and calcium release simulations lend support for a three-module decomposition of channel gating, as follows: one fast, activating, controlled by cytosolic Ca²⁺ and ATP; one fast, inhibitory, controlled by cytosolic and luminal Ca²⁺; and one slow, inactivating, controlled by cytosolic Ca²⁺ and IP₃.
   
5. Local calcium oscillations observed on a seconds scale in permeabilized cells might be driven by the slow inactivation of channels organized in discrete clusters composed of between six and 15 channels.
   
6. Out-of-phase calcium variations appear at rather constant luminal K⁺ whereas sustained K⁺ uptake by the lumen appears as a key factor in inducing parallel oscillations.
   
7. Different patterns of calcium signals, as seen in many experiments, can be induced depending on: Ca²⁺ and K⁺ fluxes between ER and cytosol, intrinsic rapidity of Ca²⁺ and K⁺ dissociation from the luminal matrix, and differences in channel distribution between cell types.
   

The model we propose here is the first to our knowledge that integrates in a quantitative manner the kinetic properties of the IP₃R with both the cytosolic and luminal calcium/potassium dynamics. On this basis we are able to reproduce various and puzzling aspects of calcium signaling observed in many experiments, as, for example, modest or large decreases (Barrero et al., 1997; Golovina and Blaustein, 1997) or increases (Golovina and Blaustein, 1997; Nguyen et al., 1998) in store filling after stimulation. So far, models for electrically excitable cells only could predict enhancement of calcium load during release.

Most important, our model can explain local calcium oscillations even in permeabilized cells, where intracellular calcium is not conserved. The nature of the allosteric transitions involved in the slow inactivation of the IP₃R (see Fig. 3 B) remains undefined, but may involve calmodulin (Patel et al., 1997; Missiaen et al., 1999) or IP₃ binding to its high affinity site (Kaftan et al., 1997; Moraru et al., 1999). An important prediction of the model is that [Ca²⁺]_cyt and [Ca²⁺]_lum variations result in- or out-of-phase if the ER transmembrane K⁺ gradient is high or low, respectively (here if the chemical equilibrium in the store is attained at tens or hundreds of µM calcium, respectively). This could have interesting implications in producing distinct local calcium signals arising from distinct ER compartments (Golovina and Blaustein, 1997) or in calcium signaling in electrically excitable cells, where parallel evolutions are generally thought to develop as, for example, in pancreatic ß-cells (Gilon et al., 1999). Nevertheless, many aspects of in vivo calcium release need further investigation. At present the data are not abundant and differences between cell types may lead to apparent discrepancies between various experimental results. However, the analysis here offers a more detailed description of the complex interplay between cytosolic and luminal control of calcium release, which undoubtedly has important effects on diverse calcium dependent processes, such as trigger of calcium influx through store-operated calcium channels (Hofer et al., 1998; Huang and Putney, 1998) birth and propagation of calcium waves (Dupont and Goldbeter, 1994; Lukyanenko et al., 1999), excitation-contraction coupling (Stern et al., 1997; Haddock et al., 1999), gene expression (Negulescu et al., 1994), or cell cycle progression (Baran, 1996).

	APPENDIX 1: STEADY-STATE OPEN PROBABILITY

TOP ABSTRACT INTRODUCTION INTRALUMINAL CALCIUM DYNAMICS MODEL FOR FAST GATING... SIMULATION OF CALCIUM SIGNALS... DISCUSSION APPENDIX 1: STEADY-STATE OPEN... APPE