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 Google Scholar Google Scholar Articles by Roux, E. Articles by Marhl, M. Search for Related Content PubMed PubMed Citation Articles by Roux, E. Articles by Marhl, M.

Role of Sarcoplasmic Reticulum and Mitochondria in Ca²⁺ Removal in Airway Myocytes

Etienne Roux ^* and Marko Marhl 

^* Laboratoire de Physiologie Cellulaire Respiratoire, Institut National de la Santé et de Recherche Médicale E 356, Université Victor Segalen Bordeaux 2, 33076 Bordeaux, France; and University of Maribor, Department of Physics, Maribor, Slovenia

Correspondence: Address reprint requests to Etienne Roux, Laboratoire de Physiologie Cellulaire Respiratoire, INSERM E 356, Université Victor Segalen Bordeaux 2, 146 Rue Léo-Saignat, 33076 Bordeaux Cedex, France. Tel.: +33-5-57-57-11-31; Fax: +33-5-57-57-16-95; E-mail: etienne.roux{at}u-bordeaux2.fr.

	ABSTRACT

TOP ABSTRACT INTRODUCTION MATERIALS AND METHODS RESULTS DISCUSSION ACKNOWLEDGEMENTS REFERENCES

 
The aim of this study was to use both a theoretical and experimental approach to determine the influence of the sarco-endoplasmic Ca²⁺-ATPase (SERCA) activity and mitochondria Ca²⁺ uptake on Ca²⁺ homeostasis in airway myocytes. Experimental studies were performed on myocytes freshly isolated from rat trachea. [Ca²⁺]_i was measured by microspectrofluorimetry using indo-1. Stimulation by caffeine for 30 s induced a concentration-graded response characterized by a transient peak followed by a progressive decay to a plateau phase. The decay phase was accelerated for 1-s stimulation, indicating ryanodine receptor closure. In Na²⁺-Ca²⁺-free medium containing 0.5 mM La³⁺, the [Ca²⁺]_i response pattern was not modified, indicating no involvement of transplasmalemmal Ca²⁺ fluxes. The mathematical model describing the mechanism of Ca²⁺ handling upon RyR stimulation predicts that after Ca²⁺ release from the sarcoplasmic reticulum, the Ca²⁺ is first sequestrated by cytosolic proteins and mitochondria, and pumped back into the sarcoplasmic reticulum after a time delay. Experimentally, we showed that the [Ca²⁺]_i decay after Ca²⁺ increase was not altered by the SERCA inhibitor cyclopiazonic acid, but was slightly but significantly modified by the mitochondria uncoupler carbonyl cyanide 4-(trifluoromethoxy)phenylhydrazone. The experimental and theoretical results indicate that, although Ca²⁺ pumping back by SERCA is active, it is not primarily involved in [Ca²⁺]_i decrease that is due, in part, to mitochondrial Ca²⁺ uptake.

	INTRODUCTION

TOP ABSTRACT INTRODUCTION MATERIALS AND METHODS RESULTS DISCUSSION ACKNOWLEDGEMENTS REFERENCES

 
Cell calcium homeostasis and calcium signaling are the result of the dynamic interactions between the ON mechanisms that lead to cytosolic Ca²⁺ concentration increase via extracellular calcium influx and/or calcium release from intracellular Ca²⁺ stores and the OFF mechanisms that tend to decrease [Ca²⁺]_i (Berridge et al., 2000). Although it is known that Ca²⁺ pumping back by sarco-endoplasmic Ca²⁺-ATPase (SERCA) and pumping out by plasma membrane Ca²⁺-ATPase (PMCA) and Na⁺-Ca²⁺ exchanger are not the unique processes involved in the OFF mechanisms and that mitochondrial Ca²⁺ uptake and Ca²⁺ binding to cytosolic proteins may also participate in the maintenance of a low [Ca²⁺]_i, SERCA and PMCA are generally considered as the major and the first structures implicated in [Ca²⁺]_i decrease after [Ca²⁺]_i increase due to cell stimulation (Berridge et al., 2000; Sanders, 2001). However, although theoretical studies have pointed out that Ca²⁺ binding to cytosolic proteins and mitochondrial uptake may greatly influence the Ca²⁺ response pattern (Marhl et al., 2000), few studies have been performed to investigate the relative participation of SERCA activity and mitochondrial Ca²⁺ uptake in Ca²⁺ signaling of airway smooth muscle cells.

Caffeine is a cell-permeant agonist of ryanodine receptors (RyR) that is largely used to investigate calcium signaling. In particular, pharmacological stimulation of RyRs by caffeine can be used to trigger calcium release from the sarcoplasmic reticulum and investigate the mechanisms activated by the initial [Ca²⁺]_i increase (Guibert et al., 1996; Roux et al., 1998; Drummond and Tuft, 1999; Pacher et al., 2000; Vallot et al., 2001; Kamishima and Quayle, 2002). On the other hand, RyRs are expressed in a large variety of cells, including airway myocytes, and may be implicated in their physiological response (Kannan et al., 1997; Prakash et al., 1998). The understanding of the effect of pharmacological stimulation of ryanodine receptors by caffeine as well as their physiological involvement in cell stimulation would be greatly improved by use of a theoretical model of calcium handling upon RyR stimulation.

The aim of the present study was first to characterize the calcium response to caffeine stimulation in rat freshly isolated cells and to build a theoretical model of calcium handling upon caffeine stimulation based on these experimental results. The model included calcium release from the sarcoplasmic reticulum (SR) through RyR stimulation and pumping back by SERCAs. In vascular smooth muscle cells, several recent studies have shown that mitochondria Ca²⁺ uptake occurs after SR Ca²⁺ release and [Ca²⁺]_i increase (Drummond and Fay, 1996; Drummond and Tuft, 1999; Pacher et al., 2000; Vallot et al., 2001; Kamishima and Quayle, 2002; Szado et al., 2003). It is also known that Ca²⁺ may bind to two classes of protein binding sites in the cytosol. The first class represents the buffering proteins such as parvalbumin, calbindin, and also C-domains of calmodulin, which bind calcium relatively slowly but with a high affinity (Falke et al., 1994; Smith et al., 1996). The second class, which is referred to as signaling proteins, comprises binding sites like N-domains of calmodulin that have very high rate constants of binding and dissociation with respect to calcium, but low affinity. Hence the model also included, as additional OFF mechanisms, mitochondrial Ca²⁺ uptake and Ca²⁺ binding to buffering proteins.

This model was further used, in combination with experimental investigations, to evaluate the involvement of these various Ca²⁺ intracellular compartments in Ca²⁺ dynamics upon Ca²⁺ release from the SR. In particular, we determined the influence of SERCA activity and mitochondrial Ca²⁺ uptake in comparison with the other buffering mechanisms on Ca²⁺ homeostasis and on [Ca²⁺]_i variations after caffeine-induced [Ca²⁺]_i increase.

	MATERIALS AND METHODS

TOP ABSTRACT INTRODUCTION MATERIALS AND METHODS RESULTS DISCUSSION ACKNOWLEDGEMENTS REFERENCES

 
Tissue preparation
Rat tracheae were obtained from male Wistar rats 10–15 weeks old, weighing 300–400 g. For each experiment, a rat was stunned and killed by cervical dissociation. Heart and lungs were removed en bloc, and the trachea was rapidly dissected out. The muscular strip located on the dorsal face of the trachea was further dissected under binocular control, as previously described (Roux et al., 2002). The epithelium-free muscular strip was cut into several pieces and the tissue was then incubated overnight (14 h) in low-Ca²⁺ (200 µM) physiological saline solution (PSS; composition given below) containing 0.5 mg ml^-1 collagenase, 0.35 mg ml^-1 pronase, 0.03 mg ml^-1 elastase, and 3 mg ml^-1 bovine serum albumin at 4°C. After this time, the muscle pieces were triturated in a fresh enzyme-free solution with a fire polished Pasteur pipette to release cells. Cells were stored for 1–3 h to attach on glass coverslips at 4°C in PSS containing 0.8 mM Ca²⁺ and used on the same day. In control experiments, immunocytochemistry was performed using monoclonal mouse anti-smooth muscle -actin antibodies and fluorescein-5'-isothiocyanate-conjugated anti-mouse immunoglobulin G antibodies to verify that the isolated cells obtained by dissociation were smooth muscle cells (data not shown).

Fluorescence measurement and estimation of [Ca²⁺]_i
Changes in [Ca²⁺]_i were monitored fluorimetrically using the Ca²⁺-sensitive probe indo-1 as previously described (Roux et al., 2002). Briefly, freshly isolated cells were loaded with indo-1 by incubation in PSS containing 1 µM indo-1 acetoxymethylester for 25 min at room temperature and then washed in PSS for 25 min. Coverslips were then mounted in a perfusion chamber and continuously superfused at room temperature. A single cell was illuminated at 360 ± 10 nm. Emitted light from that cell was counted simultaneously at 405 nm and 480 nm by two photomultipliers (P100, Nikon, Tokyo, Japan). [Ca²⁺]_i was estimated from the 405/480 ratio using a calibration for indo-1 determined within cells.

Caffeine was applied to the tested cell by a pressure ejection from a glass pipette located close to the cell. No changes in [Ca²⁺]_i were observed during test ejections of PSS (data not shown). Generally, each record of [Ca²⁺]_i response to caffeine was obtained from a different cell. Each type of experiment was repeated for the number of cells indicated in the text.

Solutions, chemicals, and drugs
Normal PSS contained (in mM): 130 NaCl, 5.6 KCl, 1 MgCl₂, 2 CaCl₂, 11 glucose, 10 Hepes, pH 7.4 with NaOH. In Ca²⁺-free solution, Ca²⁺ was removed and 0.4 mM EGTA was added. In "Ca²⁺-confining" solution, Ca²⁺ and Na⁺ were omitted and 130 mM N-methyl-D-glucamine, 0.4 mM EGTA, and 0.5 mM LaCl₃ were added to the solution. Such a solution is Ca²⁺ confining, i.e., inhibits any Ca²⁺ fluxes throughout the plasma membrane (Tribe et al., 1994), because removal of external Ca²⁺ inhibits any Ca²⁺ influx, removal of external Na⁺ inhibits the Na⁺/Ca²⁺ exchanger, and lanthanum is a nonspecific inhibitor of Ca²⁺-dependent transport activities, in particular the PMCA (Herscher and Rega, 1996) and the Na⁺/Ca²⁺ exchanger (Iwamoto and Shigekawa, 1998).

Collagenase (type CLS1) was from Worthington Biochemical Corp. (Freehold, NJ). Bovine serum albumin, cyclopiazonic acid (CPA), thapsigargin, carbonyl cyanide 4-(trifluoromethoxy)phenylhydrazone (FCCP), and caffeine (CAF) were purchased from Sigma (Saint Quentin Fallavier, France). Indo-1 AM was from Calbiochem (France Biochem, Meudon, France). Indo-1 AM, CPA, and FCCP were dissolved in dimethyl sulphoxide, for which the maximal concentration used in our experiments was <0.1% and had no effect on the resting value of the [Ca²⁺]_i nor on the variation of the [Ca²⁺]_i induced by caffeine (data not shown).

Data analysis and statistics
Cytosolic calcium concentrations ([Ca²⁺]_i) in control conditions are expressed as mean ± SE. For each protocol, experiments were repeated on two or more rats, control cells were obtained in each rat, and [Ca²⁺]_i values recorded for each protocol are expressed as a percentage of the corresponding values obtained in control cells from the rats used for that protocol. [Ca²⁺]_i values that characterize the Ca²⁺ response, i. e., baseline, peak, and plateau values, were compared between control and each specific experimental condition using unpaired student t-tests (peak and plateau values were calculated above baseline). Time-dependent decrease in [Ca²⁺]_i induced by SERCA inhibition was fitted by a sigmoidal equation using the Levenberg-Marquardt algorithm with ORIGIN software (Microcal Software Inc., Northhampton, MA). Parameters derived form curve fitting are given with SE. Comparison of the decay phase following the maximal [Ca²⁺]_i value in control versus specific conditions (i.e., Ca²⁺ confining solution, CPA, and FCCP) were performed as follows: for each cell tested, [Ca²⁺]_i measurements were repeated every 0.2 s for 10 s, time 0 corresponding to the maximal [Ca²⁺]_i. Comparison between the points obtained in control versus specific condition was done by analysis of variance (ANOVA) for repeated measurements, using SPSS software (SPSS, Chicago, IL).

Computational
Simulations were performed using the Runge-Kutta routine for numerical solving of differential equations. The software used was MADONNA (University of Berkeley, Berkeley, CA).

	RESULTS

TOP ABSTRACT INTRODUCTION MATERIALS AND METHODS RESULTS DISCUSSION ACKNOWLEDGEMENTS REFERENCES

 
Characterization of caffeine-induced [Ca²⁺]_i response
[Ca²⁺]_i response to caffeine stimulation
In this series of experiments, cells were stimulated with caffeine for 30 s in a concentration range from 0.01 to 5 mM. Stimulation with caffeine induced a [Ca²⁺]_i response characterized by a first transient peak followed by a [Ca²⁺]_i decay to a plateau slightly above the resting [Ca²⁺]_i value. At the end of the stimulation the [Ca²⁺]_i rapidly returned to baseline. Both the amplitude of the peak and the plateau depended on caffeine concentration. Typical traces of [Ca²⁺]_i response to 0.1 and 5 mM caffeine are shown in Fig. 1 A. Fig. 1 B shows the mean peak (left panel) and plateau (right panel) values for the various caffeine concentrations are given in (n = 10–35 cells for each concentration). Log IC₅₀, i.e., the logarithm of the half-maximal response concentration determined from curve fitting, was -3.51 ± 0.07 logM and -3.85 ± 0.19 logM for the peak and the plateau, respectively.

View larger version (18K): [in this window] [in a new window]   FIGURE 1  Effect of various caffeine concentrations on [Ca2+]i responses. (A) Original traces of calcium response to 0.1 mM (left trace) and 5 mM (right trace) CAF (30-s stimulation). (B) Amplitude of the calcium peak (•, left panel) and plateau (, right panel) for CAF concentrations from 0.01 mM to 5 mM. Each point is a mean value from 10 to 35 cells. Vertical bars are SE.

View larger version (21K): [in this window] [in a new window]   FIGURE 2  Effect of caffeine stimulation for various durations on [Ca2+]i. (A) Original traces of cells stimulated by 5 mM CAF for 30 s, 5 s, and 1 s. (B) Mean [Ca2+]i decay upon stimulation with 5 mM CAF for 30 s (), 5 s (), and 1 s (shaded squares). Time 0 corresponds to max [Ca2+]i increase (peak). Values are expressed as percent of the peak obtained when cells were stimulated for 30 s (reference peak), 0 corresponding to baseline. Each curve is the mean value from 10 to 12 cells. (C) Effect of CAF ejection on fluorescence intensities and ratio. CAF was ejected for 5 s on a cell preincubated for 11 min with 10 µM CPA. During CAF ejection, fluorescence quenching at both wavelengths is observed, without alteration in the fluorescence ratio. Fluorescence intensities and ratio are in arbitrary units.

Role of transmembrane Ca²⁺ fluxes in the response to caffeine stimulation
The possible implication of calcium fluxes through the plasma membrane was assessed by a series of experiments in Ca²⁺-free medium and in Ca²⁺-confining solution. In normal PSS, baseline value was 139.7 ± 4.2 nM (n = 42). Removal of external Ca²⁺ or incubation in Ca²⁺-confining solution did not significantly modify the resting [Ca²⁺]_i value (88.0 ± 3.2% control, n = 23 and 99.9 ± 3.9% control, n = 23, respectively). When cells in Ca²⁺-free or Ca²⁺-confining medium were stimulated with 5 mM caffeine for 30 s, neither the peak nor the plateau was significantly different from control. Fig. 3 A shows the mean results, expressed as percent to control conditions, for the peak (left panel) and the plateau (right panel).

View larger version (33K): [in this window] [in a new window]   FIGURE 3  Effect of Ca2+-free and Ca2+-confining solutions on response to 5 mM caffeine. (A) Peak and plateau values in response to 30-s CAF stimulation, in Ca2+-free medium (open columns, n = 21) and Ca2+-confining medium (shaded columns, n = 11). (B) Original trace of Ca2+ response to two successive stimulations by CAF for 1 s, at a 15-s interval in the presence of external Ca2+ (left), in Ca2+-free medium (middle), and in Ca2+-confining solution (right). (C) Left panel, mean values of the Ca2+ response to two successive stimulations by CAF for 1 s, at a 15-s interval in Ca2+-free medium (open columns, n = 13) and Ca2+-confining solution (shaded columns, n = 23). (C) Right panel, mean [Ca2+]i decay after the first CAF-induced peak in control (•, n = 22) and Ca2+-confining solution (, n = 23), as percentage of the control peak, 0 corresponding to baseline. Vertical bars are SE.

Mathematical modeling
Presentation of the model
In accordance with our experiments, the model takes into account the calcium exchange between the cytosol and the intracellular calcium stores, neglecting any transplasmalemmal Ca²⁺ fluxes. We focus on four different intracellular calcium stores: the SR, mitochondria, and two classes of calcium-binding proteins in the cytosol, i.e., signaling and buffering proteins. The signaling proteins characterize protein-binding sites with high affinity and low capacity, whereas the buffering proteins characterize protein-binding sites with low affinity and high capacity. Regarding the SR, three different calcium fluxes are included in the model: the ATP-dependent calcium uptake from the cytosol into the SR (J_SERCA), the Ca²⁺ efflux from the SR through ryanodine sensitive Ca²⁺ channels following the calcium-induced calcium release (CICR) mechanism (J_RyR), and an additional Ca²⁺ leak flux from the SR into the cytosol (J_leak). For the exchange of Ca²⁺ between the mitochondria and the cytosol we take into account active Ca²⁺ uptake by mitochondrial uniporters (J_in) and calcium release through Na⁺/Ca²⁺ and H⁺/Ca²⁺ exchangers (J_out). In the cytosol Ca²⁺ binding to signaling and buffering proteins is considered.

The concentration of free-Ca²⁺ binding sites on signaling proteins, SPr, can be calculated by applying a rapid-equilibrium approximation to the fast binding reactions (see, e.g., Wagner and Keizer, 1994; Heinrich and Schuster, 1996; Marhl et al., 1998a; Höfer et al., 2001). Together with the conservation relation for the total concentration of Ca²⁺ binding sites on signaling proteins, SPr_tot, the concentration of free-Ca²⁺ binding sites on signaling proteins, SPr, is given by the following equation:

(1)

where is the dissociation constant of binding sites on signaling proteins.

The rapid-equilibrium approximation for the fast Ca²⁺-binding to signaling proteins is justified in view of the very high values of the rate constants for signaling proteins (Wagner and Keizer, 1994; Smith et al., 1996). The rate constants for binding and dissociation of 90–500 µM^-1s^-1 and 300–500 s^-1, respectively (Falke et al., 1994; Smith et al., 1996), imply a time constant of <0.01 s. Even under consideration of diffusional resistance and competition with Mg²⁺, this time constant is much smaller than the duration of any calcium pulse.

Using the conservation relation for the total concentration of Ca²⁺ binding sites on signaling proteins, , the concentration of bound Ca²⁺ binding sites on signaling proteins can be calculated as:

(2)

Taking into account the conservation relations for the total cellular calcium, , and the total concentration of Ca²⁺ binding sites on buffering proteins, , we obtain equations for bound and free-Ca²⁺ binding sites on buffering proteins, CaBPr and BPr, respectively:

(3)

(4)

Here and represent the volume ratio between the SR and the cytosol and between the mitochondria and the cytosol, respectively. Assuming very fast unsaturated buffering of Ca²⁺ in the SR and mitochondrial compartments, we use factors and , which are the free/total calcium concentration ratios, in the SR and in the mitochondria, respectively (Marhl et al., 1998b; Haberichter et al., 2001). Applying the rapid-equilibrium approximation for relating the free-calcium concentration in each organelle, (j stands for SR and m), with the total concentration in the corresponding organelle, , leads to , where is the dissociation constant of the binding sites on buffering proteins in the organelle. Taking into account that , we can assume that nearly equals the total buffer concentration, . Therefore, we consider Because the values for and have not been experimentally determined for airway myocytes, we take values that are in the range of experimentally obtained values for other cell types. For values usually span from 0.0025 (Li et al., 1995) to 0.01 (Smith et al., 1996), whereas can be found in the range from 0.00001 to 0.01 (Babcock et al., 1997; Fall and Keizer, 2001; Chalmers and Nicholls, 2003).

The time dependence of the free cytosolic calcium concentration, , is determined by Ca²⁺ fluxes across the SR membrane, by Ca²⁺ exchange between the cytosol and mitochondria, and by Ca²⁺ binding to signaling and buffering proteins in the cytosol:

(5)

Here and denote the off and on rate constants, respectively, of the Ca²⁺ binding to the buffering proteins. The first factor in Eq. 5, characterizing Ca²⁺ binding to signaling proteins, is obtained by the method of eliminating fast reactions (see Heinrich and Schuster, 1996; Marhl et al., 1998a).

The equation for the free-calcium concentration in the SR, , is linked with the fluxes across the SR membrane as follows:

(6)

The equation for the free-Ca²⁺ concentration in mitochondria, , reads:

(7)

Here, relates fluxes, i.e., , like in previous papers by Wagner and Keizer (1994) and Smith et al. (1996). Under the condition of very fast unsaturated buffering of Ca²⁺ in the organelles used above, this is in agreement with using the constant for relating the concentrations in Eq. 3. Differentiating with respect to time gives and dividing by gives the above formula .

There are three genetically distinct isoforms of RyR in Mammals, RyR1, RyR2, and RyR3. In smooth muscle, several isoforms may be expressed (Marin et al., 1999; Ogawa et al., 2000). In airway smooth muscle, RT-PCR in porcine trachea has revealed that both RyR2 and RyR3 mRNA were present (Kannan et al., 1997). However, in human bronchial smooth muscle, we have previously showed using RT-PCR and RNase protection assay that RyR3 was the unique isoform expressed (Hyvelin et al., 2000). On this basis, our model was built on RyR3 activity. The Ca²⁺ flux through ryanodine channels, , is caffeine-activated. The open probability depends on caffeine concentration, Caff, and follows the Hill kinetics with . As driving force for the Ca²⁺ channel flux the concentration gradient across the SR membrane is taken into account and the equation for reads:

(8)

where is the rate constant of the ryanodine channels and represents the half-saturation constant of the ryanodine channels for caffeine. The kinetics for caffeine activation was obtained by fitting our experimental data (see Fig. 1 B). This is also fully in accordance with previous experimental studies in airway smooth muscle cells (Roux et al., 1998). On the basis of experimental investigations (Ogawa et al., 2000) on the Ca²⁺-induced Ca²⁺ release (CICR) mechanism of ryanodine channels, CICRMg, the effect of Mg²⁺ on the activity of ryanodine receptors is taken into account as follows:

(9)

where and represent the half-saturation constants of activate sites for Ca²⁺ and Mg²⁺, respectively, and and represent the half-saturation constants of inactivate sites for Ca²⁺ and Mg²⁺, respectively.

For the SERCA flux into the SR lumen, , the Hill kinetics with is taken:

(10)

where is the rate constant of the SERCAs and stands for the half-saturation constant of SERCAs.

The leak flux depends on the concentration gradient across the SR membrane and a simple relation is taken:

(11)

where is the rate constant for Ca²⁺ leak flux through the SR membrane and has been estimated on the basis of our experimental data.

There is experimental evidence of a very fast and effective calcium sequestration by mitochondria through a specific uniporter (Hehl et al., 1996; Applegate et al., 1997; Babcock et al., 1997). In some cases the Ca²⁺ uptake by mitochondria can be extremely fast due to a mechanism called the rapid mode (RaM) (Gunter et al., 2000; Rizzuto et al., 2000). Based on the experimental results that Ca²⁺ sequestration takes place at free cytosolic calcium levels of >0.5–1.0 µM (Jouaville et al., 1995; Bernardi and Petronilli, 1996; Hehl et al., 1996; Herrington et al., 1996; Hoth et al., 1997; Ricken et al., 1998) a step-like kinetics is considered for the mitochondrial Ca²⁺ uptake by uniporters, (Marhl et al., 1998a; Grubelnik et al., 2001; Haberichter et al., 2001).

(12)

where represents the maximal permeability of the uniporters in the mitochondrial membrane, and represents the half-saturation for Ca²⁺ (Marhl et al., 2000). As in our previous publication (Marhl et al., 1998a), in factor the constant value of the mitochondrial transmembrane potential, , is implicitly included. The potential difference is usually strongly changed only by fast release of calcium through the permeability transient pore (PTP). Note that under normal physiological conditions just a slow release of calcium from the mitochondria takes place (Bernardi and Petronilli, 1996; Eriksson et al., 1999; Svichar et al., 1999).

For the mitochondrial Ca²⁺ release through Na⁺/Ca²⁺ and H⁺/Ca²⁺ exchangers we consider a simple linear dependency on (Marhl et al., 1998b; Grubelnik et al., 2001):

(13)

where is the maximal rate for calcium efflux from mitochondria.

Model parameters used in our calculations are given in figure captions. The parameter values for Ca²⁺ fluxes across the SR membrane were estimated in accordance with our own experimental results. Data regarding the RyR activity were taken from the experiments done by Ogawa et al. (2000). The parameter values regarding Ca²⁺ binding to signaling and buffering proteins were taken from the experiments in the literature and were extensively discussed in our previous paper (Marhl et al., 1998a). The parameter values regarding the mitochondrial Ca²⁺ handling were also taken from the experimental results published in the literature and were exhaustively discussed in our previous papers (Marhl et al., 1998b, 2000; Haberichter et al., 2001).

Modeling of caffeine-induced [Ca²⁺]_i response
Cell stimulation with different concentrations of caffeine was simulated by varying the model parameter Caff. The calculated traces for Caff = 0.1 mM and 5 mM for 30 s are shown in Fig. 4 A, and correspond to the experimental results presented in Fig. 1 A. Fig. 4 B represents the calculated concentration-dependent RyR activity. The sigmoidal relationship between caffeine concentration and RyR activity is similar to that between caffeine concentration and [Ca²⁺]_i response showed in Fig. 1 B. Model calculations of [Ca²⁺]_i responses to 5-mM caffeine exposures for 1 s, 5 s, and 30 s are shown in Fig. 4 C (dotted line, dashed line, and solid line, respectively). They correspond to the experimental traces shown in Fig. 2 A, and reflect fast closure of RyR at the end of caffeine exposure.

View larger version (19K): [in this window] [in a new window]   FIGURE 4  Simulation of Ca2+ response to caffeine. (A) Simulated response to 0.1 (left) and 5 (right) mM CAF for 30 s. (B) Concentration-activity curve of RyR. (C) Simulated responses to 30-s (solid line), 5-s (dashed line), and 1-s (dotted line) caffeine (5 mM). Parameter values are SR = 0.01; m = 0.01; ßSR = 0.0025; ßm = 0.0025; Catot = 50 µM; SPrtot = 90 µM; BPrtot = 120 µM; Mgi = 500 µM; kRyR = 2000 s-1; kSERCA = 1 µM s-1; kleak = 0.02 s-1; kin = 20 µM s-1; kout = 0.1 s-1; kon = 0.1 µM-1 s-1; koff = 0.01 s-1; KCaff = 250 µM; KA,Ca = 2.5 µM; KA,Mg = 75 µM; KI,Ca = 400 µM; KI,Mg = 300 µM; KSERCA = 0.1 µM; Km = 1 µM; and KSPr = 5 µM.

View larger version (16K): [in this window] [in a new window]   FIGURE 5  Model prediction of Ca2+ distribution after Ca2+ release from the SR. (A and B) Model predictions on Ca2+ handling upon 5-mM CAF stimulation. Stimulation begins at time 0 s and runs up to 1 s (double arrow). Parameter values are the same as in Fig. 4. Concentration in the cytosol (Cai) is indicated by the dashed line. (A) Concentration in the SR (CaSR) and mitochondria (Cam). (B) Ca2+ bound to fast-kinetics (CaSPr) and to slow-kinetics (CaBPr) cytosolic proteins. (C) Predicted long-time variations in SR Ca2+ concentration for 0.1 (dashed line) and 5 mM CAF (solid line). Stimulation begins at time 0 s and runs up to 1 s (arrow). Parameter values are the same as in Fig. 4.

View larger version (18K): [in this window] [in a new window]   FIGURE 6  Model prediction of the effect of SERCA blockade on the Ca2+ response to two successive stimulations by 5 mM caffeine (1 s). (A) Predicted [Ca+]i response to two successive 1-s stimulations by 5 mM CAF at a 15-s interval, with active (left trace) and inactivated (right trace) SERCA. (B) Predicted values of the decay phase after the first [Ca+]i rise, with active (•) and inactivated () SERCA. Parameter values are the same as in Fig. 4.

View larger version (23K): [in this window] [in a new window]   FIGURE 7  Effect of CPA (10 µM) on [Ca2+]i response to successive 1-s stimulations with 5 mM caffeine. (A) Original traces of cells stimulated in the absence (left trace) and in the presence of CPA (right trace). (B) Effect of CPA on the calcium response to two successive CAF stimulations. Values were obtained from 21 cells, and are expressed as a percentage of control. Vertical bars are SE. *p < 0.05. (C) Mean [Ca2+]i decay after the first CAF stimulation in the absence (•) and in the presence () of 10 µM CPA. Each curve in the mean value of 21 cells in each experimental conditions, expressed as percentage of the control peak, 0 corresponding to baseline.

View larger version (19K): [in this window] [in a new window]   FIGURE 8  Effect of SERCA blockade on SR Ca2+ load. (A) Typical traces of [Ca2+]i response to two successive 5-s CAF stimulations in cells preincubated with 10 µM CPA for 7 min (left trace), 10 min (middle trace), and 13 min (right trace). (B) Individual Ca2+ response to 5 mM CAF of cells incubated in Ca2+-free medium (crosses) and incubated in Ca2+-free medium with 10 µM CPA (•) or with 1 µM TG (). Abscissa: duration of incubation. Ordinate: first peak amplitude in percent of control, i.e., response of cells in normal PSS. [Ca2+]i in Ca2+-free medium are fitted by linear regression; [Ca2+]i responses in the presence of CPA and TG are fitted by a sigmoidal equation.

Model prediction of mitochondrial Ca²⁺uptake blockade
The theoretical and experimental results presented above show that Ca²⁺ pumping back by SERCA is not significantly responsible for the [Ca²⁺]_i decrease after Ca²⁺ release from the SR and, according to the model, Ca²⁺ uptake into mitochondria may be more important. To test this hypothesis, we compared model predictions of Ca²⁺ dynamics with active versus blocked mitochondrial Ca²⁺ uptake. We simulated the [Ca²⁺]_i response to 1-s stimulation by 5 mM caffeine, with and without active Ca²⁺ uniporter, and determined the respective [Ca²⁺]_i decays, as shown Fig. 9 A. According to the model, mitochondrial inhibition results in a slight change in the slope of the Ca²⁺ decay, indicating that mitochondrial Ca²⁺ uptake is actually involved in cytosolic Ca²⁺ clearance after Ca²⁺ release from the SR, but has a slight influence on the shape of the [Ca²⁺]_i decay.

View larger version (13K): [in this window] [in a new window]   FIGURE 9  Predicted and experimental effect of mitochondrial Ca2+ uptake inhibition on the Ca2+ response to 5 mM caffeine (1 s). (A) Predicted values of the decay phase after CAF stimulation, with active (solid line) and inactivated (dashed line) mitochondrial Ca2+ uptake. Parameter values are the same as in Fig. 4. (B) Mean [Ca2+]i decay after CAF-induced peak in the absence (n = 21, •) and in the presence (n = 25, ) of 5 µM FCCP, as a percentage of the control peak, 0 corresponding to baseline.

	DISCUSSION

TOP ABSTRACT INTRODUCTION MATERIALS AND METHODS RESULTS DISCUSSION ACKNOWLEDGEMENTS REFERENCES

 
Our experimental study has shown that calcium signaling upon short caffeine stimulation results in Ca²⁺ dynamics within the cell without significant involvement of Ca²⁺ fluxes through the plasma membrane. On the basis on these experimental data, we build a relevant theoretical model of RyR stimulation and subsequent Ca²⁺ handling. Using this model in combination with the experimental approach, we used short-time stimulation by caffeine to investigate the OFF mechanisms, i.e., the mechanisms implicated in the return of [Ca²⁺]_i to baseline after [Ca²⁺]_i increase. Our experimental results, in accordance with our theoretical model, indicate that although Ca²⁺ pumping back by SERCA is active after Ca²⁺ release from SR upon RyR stimulation, it is not primarily involved in [Ca²⁺]_i decrease. Ca²⁺ uptake by mitochondria slightly but significantly shapes the [Ca²⁺]_i decay that may also depend on additional intracellular buffering processes.

Our cell model of RyR stimulation is a "closed cell" one, i.e., in which Ca²⁺ fluxes through the plasma membrane are not considered and the total Ca²⁺ concentration within the cell is kept constant. Although transmembrane Ca²⁺ fluxes are likely to be present in these cells, our results showed that such an approximation is relevant in our experimental conditions. Indeed, Ca²⁺ influx was not triggered by caffeine-induced Ca²⁺ release, because the [Ca²⁺]_i response was not altered in Ca²⁺-free medium. Store-operated channels, activated by the depletion of intracellular Ca²⁺ stores, have been described in a variety of smooth muscle cells. Some authors have speculated that the main physiological role of store-operated channels may be the refilling of the SR (Marin et al., 1999; Ng and Gurney, 2001). The fact that the second Ca²⁺ response to two successive caffeine stimulations, which is actually due to SR refilling because it is abolished in the presence of the SERCA blocker CPA, was not modified in Ca²⁺-free medium indicates that such a mechanism is not involved in our experimental conditions. Usually, Ca²⁺ extrusion by PMCA and/or Na⁺/Ca²⁺ exchanger is considered as an important OFF mechanism (Berridge et al., 2000; Sanders, 2001). In Na⁺- and Ca²⁺-free solution containing 0.5 mM lanthanum, that inhibits any Ca²⁺ influx or efflux (Tribe et al., 1994), the decay phase of the Ca²⁺ response was not modified, indicating that Ca²⁺ extrusion is not implicated in the quick [Ca²⁺]_i return to baseline. This is in accordance with previous studies that showed that the [Ca²⁺]_i response to contractile agonists like acetylcholine is not altered by removal of extracellular Ca²⁺ (Roux et al., 1996, 1997; Prakash et al., 1998; Bergner and Sanderson, 2002). In the absence of compensatory Ca²⁺ influx, fast Ca²⁺ extrusion would rapidly induce a complete loss of intracellular Ca²⁺ pool and hence loss of Ca²⁺ response.

The predicted values for the resting Ca²⁺ levels, that are 40 nM, are lower than the experimental values, that were 130 nM, in accordance with our previous results (Roux et al., 1997, 1998). However, due to relative imprecision of the calculation of absolute [Ca²⁺]_i values and interindividual variations so that values measured in some cells from our experiments were similar to that predicted but the model, the predicted values, though lower than the mean estimated [Ca²⁺]_i, can be considered as in a physiological range.

We analyzed the OFF mechanisms after [Ca²⁺]_i increase using RyR stimulation. Single-channel activity measurement in lipid bilayer indicates that RyR2 and RyR3 differ in their sensitivity to Ca²⁺ inactivation and Mg²⁺ inhibition, so that for RyR3 Ca²⁺-induced Ca²⁺ release may be reduced in the presence of physiological Mg²⁺ concentration (Ogawa et al., 2000). In accordance with the simulated Ca²⁺ response, stimulation by various caffeine concentrations showed a sigmoidal relationship between the amplitude of the Ca²⁺ response and the logarithmic concentration of caffeine. The fact that the response is not an all-or-none phenomenon suggests that Ca²⁺-induced Ca²⁺ release is limited, which is consistent with a predominant RyR3 expression.

To investigate the mechanisms responsible for Ca²⁺ removal in smooth muscle cells, some authors have used voltage pulse to induce Ca²⁺ influx through voltage-dependent L-type Ca²⁺ channels (Kamishima and McCarron, 1998; Shmigol et al., 1999). In such a protocol, initial [Ca²⁺]_i increase is due to additional Ca²⁺ influx into the cell without SR Ca²⁺ depletion. However, in airway smooth muscle, although L-type Ca²⁺ channels may be activated during stimulation (Tomasic et al., 1992), contractile agonists such as acetylcholine act primarily via Ca²⁺ release from the SR (Roux et al., 1996, 1997; Prakash et al., 1998; Bergner and Sanderson, 2002). Hence, in this type of cells, it seems physiologically relevant to induce Ca²⁺ release from the SR rather than extracellular Ca²⁺ influx to investigate the subsequent OFF mechanisms. Analysis of fluorescence quenching and comparison of the decay phase after 1-, 5-, and 30-s caffeine stimulation indicate that caffeine washout and RyR closure occurred 2–3 s after the beginning of caffeine ejection. Because the maximal [Ca²⁺]_i increase occurred approximately in the same time (see Fig. 2), we assumed that during the decay phase, the duration of which was 10 s after maximal [Ca²⁺]_i increase, the RyRs were closed. Hence the mechanisms involved in cytosolic clearance in our experiments do not seem to be specific to RyR stimulation and are likely to be also implicated upon airway smooth muscle stimulation by physiological agonists that induce Ca²⁺ release from the SR.

Usually, Ca²⁺ uptake by SERCAs is considered as the main OFF mechanism (Marin et al., 1999; Berridge et al., 2000; Sanders, 2001). The analysis of the decay phase after short caffeine stimulation in vitro and in silico demonstrated that, although Ca²⁺ pumping back is actually active and does reload the SR, SERCA activity does not significantly modulate the shape of the [Ca²⁺]_i decay. Several authors have used CPA to inhibit the SERCA activity, with significant effect on Ca²⁺ decrease, including in airway myocytes (Sims et al., 1996, 1997; Yoshikawa et al., 1996; Shmigol et al., 1999). Global comparison of these different studies is difficult because the experimental protocols used to trigger [Ca²⁺]_i increase as well as to analyze [Ca²⁺]_i decay differed between studies. However, it appears in these studies that even if the [Ca²⁺]_i decay rate was altered in the presence of CPA, the general pattern was not deeply modified. Hence, [Ca²⁺]_i return to baseline should be attributed to mechanisms other than SERCA activity.

The model prediction about Ca²⁺ sequestration in mitochondria is in agreement with experimental results in some other systems, showing that mitochondria indeed sequester a large amount of Ca²⁺ released from the sarco/endoplasmic reticulum (for review see Schuster et al., 2002). For example, in chromaffin cells, around 80% of the Ca²⁺ released from the ER is cleared first into mitochondria (Babcock and Hille, 1998). In vascular smooth muscle cells, several recent studies have shown that mitochondria Ca²⁺ uptake occurs after SR Ca²⁺ release and [Ca²⁺]_i increase (Drummond and Fay, 1996; Drummond and Tuft, 1999; Pacher et al., 2000; Vallot et al., 2001; Kamishima and Quayle, 2002; Szado et al., 2003) but, to the best of our knowledge, this issue has not been previously addressed in airway smooth muscle cells. It should be noted that although our model predicts that mitochondrial Ca²⁺ uptake is effective after [Ca²⁺]_i increase, simulated inhibition of mitochondria uniporter does not greatly alter the shape of the [Ca²⁺]_i decay after 1-s caffeine stimulation, and this was confirmed by our experimental results. In particular, the [Ca²⁺]_i decay phase itself was slowed in the presence of FCCP, but the subsequent resting value did not differ from control.

In conclusion, we have built a relevant theoretical model of Ca²⁺ handling upon RyR stimulation in airway myocytes. We have experimentally confirmed the model predictions that Ca²⁺ pumping back into the SR by SERCA, though effective, is not primarily involved in [Ca²⁺]_i decay upon Ca²⁺ release from the SR, and that Ca²⁺ uptake by mitochondria slightly but significantly shapes the [Ca²⁺]_i decay. The model provides a possible scenario for Ca²⁺ handling in which cytosolic Ca²⁺-binding proteins play a key role in cytosolic Ca²⁺ clearance. Further studies are needed to check the model predictions about the role of these buffering processes in [Ca²⁺]_i modulation.

	ACKNOWLEDGEMENTS

TOP ABSTRACT INTRODUCTION MATERIALS AND METHODS RESULTS DISCUSSION ACKNOWLEDGEMENTS REFERENCES

 
The authors thank Professors J. P. Mazat and R. Marthan from Bordeaux for fruitful discussions, Dr. C. Guibert from Bordeaux for reading the manuscript, and Professor S. Schuster from Jena and Professors R. Heinrich and T. Höfer from Berlin for checking the mathematical expressions used in the model.

This work was financially supported by the French Ministère des Affaires Étrangères and the Slovenian Ministry of Education, Science and Sport (bilateral project PROTEUS, F-No. 03811NH/SI-No. FR-2002-04), and INSERM.

Submitted on July 2, 2003; accepted for publication December 24, 2003.

	REFERENCES

TOP ABSTRACT INTRODUCTION MATERIALS AND METHODS RESULTS DISCUSSION ACKNOWLEDGEMENTS REFERENCES

 
Albrecht, M. A., S. L. Colegrove, and D. D. Friel. 2002. Differential regulation of ER Ca²⁺ uptake and release rates accounts for multiple modes of Ca²⁺-induced Ca²⁺ release. J. Gen. Physiol. 119:211–233.[Abstract/Free Full Text]

Applegate, T. L., A. Karjalainen, and F. L. Bygrave. 1997. Rapid Ca²⁺ influx induced by the action of dibutylhydroquinone and glucagon in the perfused rat liver. Biochem. J. 323:463–467.[Medline]

Babcock, D. F., J. Herrington, P. C. Goodwin, Y. B. Park, and B. Hille. 1997. Mitochondrial participation in the intracellular Ca²⁺ network. J. Cell Biol. 136:833–844.[Abstract/Free Full Text]

Babcock, D. F., and B. Hille. 1998. Mitochondrial oversight of cellular Ca²⁺ signaling. Curr. Opin. Neurobiol. 8:398–404.[Medline]

Bergner, A., and M. J. Sanderson. 2002. Acetylcholine-induced calcium signaling and contraction of airway smooth muscle cells in lung slices. J. Gen. Physiol. 119:187–198.[Abstract/Free Full Text]

Bernardi, P., and V. Petronilli. 1996. The permeability transition pore as a mitochondrial calcium release channel: a critical appraisal. J. Bioenerg. Biomembr. 28:131–138.[Medline]

Berridge, M. J., P. Lipp, and M. D. Bootman. 2000. The versatility and universality of calcium signalling. Nat. Rev. Mol. Cell Biol. 1:11–21.[Medline]

Chalmers, S., and D. G. Nicholls. 2003. The relationship between free and total calcium concentrations in the matrix of liver and brain mitochondria. J. Biol. Chem. 278:19062–19070.[Abstract/Free Full Text]

Drummond, R. M., and F. S. Fay. 1996. Mitochondria contribute to Ca²⁺ removal in smooth muscle cells. Pflugers Arch. 431:473–482.[Medline]

Drummond, R. M., and R. A. Tuft. 1999. Release of Ca²⁺ from the sarcoplasmic reticulum increases mitochondrial [Ca²⁺] in rat pulmonary artery smooth muscle cells. J. Physiol. 516:139–147.[Abstract/Free Full Text]

Eriksson, O., P. Pollesello, and E. Geimonen. 1999. Regulation of total mitochondrial Ca²⁺ in perfused liver is independent of the permeability transition pore. Am. J. Physiol. 276:C1297–C1302.[Medline]

Falke, J. J., S. K. Drake, A. L. Hazard, and O. B. Peersen. 1994. Molecular tuning of ion binding to calcium signaling proteins. Q. Rev. Biophys. 27:219–290.[Medline]

Fall, C. P., and J. E. Keizer. 2001. Mitochondrial modulation of intracellular Ca²⁺ signaling. J. Theor. Biol. 210:151–165.[Medline]

Grubelnik, V., A. Z. Larsen, U. Kummer, L. F. Olsen, and M. Marhl. 2001. Mitochondria regulate the amplitude of simple and complex calcium oscillations. Biophys. Chem. 94:59–74.[Medline]

Guibert, C., R. Marthan, and J. P. Savineau. 1996. Angiotensin II-induced Ca²⁺-oscillations in vascular myocytes from the rat pulmonary artery. Am. J. Physiol. 270:L637–L642.[Medline]

Gunter, T. E., L. Buntinas, G. Sparagna, R. Eliseev, and K. Gunter. 2000. Mitochondrial calcium transport: mechanisms and functions. Cell Calcium. 28:285–296.[Medline]

Gurney, A. M., R. M. Drummond, and F. S. Fay. 2000. Calcium signalling in sarcoplasmic reticulum, cytoplasm and mitochondria during activation of rabbit aorta myocytes. Cell Calcium. 27:339–351.[Medline]

Haberichter, T., M. Marhl, and R. Heinrich. 2001. Birhythmicity, trirhythmicity and chaos in bursting calcium oscillations. Biophys. Chem. 90:17–30.[Medline]

Hehl, S., A. Golard, and B. Hille. 1996. Involvement of mitochondria in intracellular calcium sequestration by rat gonadotropes. Cell Calcium. 20:515–524.[Medline]

Heinrich, R., and S. Schuster. 1996. The Regulation of Cellular Systems. Chapman Hall, New York.

Herrington, J., Y. B. Park, D. F. Babcock, and B. Hille. 1996. Dominant role of mitochondria in clearance of large Ca²⁺ loads from rat adrenal chromaffin cells. Neuron. 16:219–228.[Medline]

Herscher, C. J., and A. F. Rega. 1996. Pre-steady-state kinetic study of the mechanism of inhibition of the plasma membrane Ca²⁺-ATPase by lanthanum. Biochemistry. 35:14917–14922.[Medline]

Höfer, T., A. Politi, and R. Heinrich. 2001. Intercellular Ca²⁺ wave propagation through gap-junctional Ca²⁺ diffusion: a theoretical study. Biophys J. 80:75–87.[Abstract/Free Full Text]

Hoth, M., C. M. Fanger, and R. S. Lewis. 1997. Mitochondrial regulation of store-operated calcium signaling in T lymphocytes. J. Cell Biol. 137:633–648.[Abstract/Free Full Text]

Hyvelin, J. M., C. Martin, E. Roux, R. Marthan, and J. P. Savineau. 2000. Human isolated bronchial smooth muscle contains functional ryanodine/caffeine-sensitive Ca-release channels. Am. J. Respir. Crit. Care Med. 162:687–694.[Abstract/Free Full Text]

Iwamoto, T., and M. Shigekawa. 1998. Differential inhibition of Na⁺/Ca²⁺ exchanger isoforms by divalent cations and isothiourea derivative. Am. J. Physiol. 275:C423–C430.[Medline]

Jouaville, L. S., F. Ichas, E. L. Holmuhamedov, P. Camacho, and J. D. Lechleiter. 1995. Synchronization of calcium waves by mitochondrial substrates in Xenopus laevis oocytes. Nature. 377:438–441.[Medline]

Kamishima, T., and J. G. McCarron. 1998. Ca²⁺ removal mechanisms in rat cerebral resistance size arteries. Biophys. J. 75:1767–1773.[Abstract/Free Full Text]

Kamishima, T., and J. M. Quayle. 2002. Mitochondrial Ca²⁺ uptake is important over low [Ca²⁺]_i range in arterial smooth muscle. Am. J. Physiol. Heart Circ. Physiol. 283:H2431–H2439.