A new approach to the detection and statistical classification of Ca2+ sparks

The availability of high-speed, two-dimensional (2-D) confocal microscopes and the expanding armamentarium of fluorescent probes presents unprecedented opportunities and new challenges for studying the spatial and temporal dynamics of cellular processes. The need to remove subjectivity from the detection process, the difficulty of the human eye to detect subtle changes in fluorescence in these 2-D images, and the large volume of data produced by these confocal microscopes call for the need to develop algorithms to automatically mark the changes in fluorescence. These fluorescence signal changes are often subtle, so the statistical estimate of the likelihood that the detected signal is not noise is an integral part of the detection algorithm. This statistical estimation is fundamental to our new approach to detection; in earlier Ca(2+) spark detectors, this statistical assessment was incidental to detection. Importantly, the use of the statistical properties of the signal local to the spark, instead of over the whole image, reduces the false positive and false negative rates. We developed an automatic spark detection algorithm based on these principles and used it to detect sparks on an inhomogeneous background of transverse tubule-labeled rat ventricular cells. Because of the large region of the cell surveyed by the confocal microscope, we can detect a large enough number of sparks to measure the dynamic changes in spark frequency in individual cells. We also found, in contrast to earlier results, that cardiac sparks are spatially symmetric. This new approach puts the detection of fluorescent signals on a firm statistical foundation.     

Multidimensional detection and analysis of Ca2+ sparks in cardiac myocytes

Examining calcium spark morphology and its relationship to the structure of the cardiac myocyte offers a direct means of understanding excitation-contraction coupling mechanisms. Traditional confocal line scanning achieves excellent temporal spark resolution but at the cost of spatial information in the perpendicular dimension. To address this, we developed a methodology to identify and analyze sparks obtained via two-dimensional confocal or charge-coupled device microscopy. The technique consists of nonlinearly subtracting the background fluorescence, thresholding the data on the basis of noise level, and then localizing the spark peaks via a generalized extrema test, while taking care to detect and separate adjacent peaks. In this article, we describe the algorithm, compare its performance to a previously validated spark detection algorithm, and demonstrate it by applying it to both a synthetic replica and an experimental preparation of a two-dimensional isotropic myocyte monolayer exhibiting sparks during a calcium transient. We find that our multidimensional algorithm provides better sensitivity than the conventional method under conditions of temporally heterogeneous background fluorescence, and the inclusion of peak segmentation reduces false negative rates when spark density is high. Our algorithm is robust and can be effectively used with different imaging modalities and allows spark identification and quantification in subcellular, cellular, and tissue preparations.     

Amplitude distribution of calcium sparks in confocal images: theory and studies with an automatic detection method.

Determination of the calcium spark amplitude distribution is of critical importance for understanding the nature of elementary calcium release events in striated muscle. In the present study we show, on general theoretical grounds, that calcium sparks, as observed in confocal line scan images, should have a nonmodal, monotonic decreasing amplitude distribution, regardless of whether the underlying events are stereotyped. To test this prediction we developed, implemented, and verified an automated computer algorithm for objective detection and measurement of calcium sparks in raw image data. When the sensitivity and reliability of the algorithm were set appropriately, we observed highly left-skewed or monotonic decreasing amplitude distributions in skeletal muscle cells and cardiomyocytes, confirming the theoretical predictions. The previously reported modal or Gaussian distributions of sparks detected by eye must therefore be the result of subjective detection bias against small amplitude events. In addition, we discuss possible situations when a modal distribution might be observed.     

Long-Lasting Sparks: Multi-Metastability and Release Competition in the Calcium Release Unit Network

Calcium (Ca) sparks are elementary events of biological Ca signaling. A normal Ca spark has a brief duration in the range of 10 to 100 ms, but long-lasting sparks with durations of several hundred milliseconds to seconds are also widely observed. Experiments have shown that the transition from normal to long-lasting sparks can occur when ryanodine receptor (RyR) open probability is either increased or decreased. Here, we demonstrate theoretically and computationally that long-lasting sparks emerge as a collective dynamical behavior of the network of diffusively coupled Ca release units (CRUs). We show that normal sparks occur when the CRU network is monostable and excitable, while long-lasting sparks occur when the network dynamics possesses multiple metastable attractors, each attractor corresponding to a different spatial firing pattern of sparks. We further highlight the mechanisms and conditions that produce long-lasting sparks, demonstrating the existence of an optimal range of RyR open probability favoring long-lasting sparks. We find that when CRU firings are sparse and sarcoplasmic reticulum (SR) Ca load is high, increasing RyR open probability promotes long-lasting sparks by potentiating Ca-induced Ca release (CICR). In contrast, when CICR is already strong enough to produce frequent firings, decreasing RyR open probability counter-intuitively promotes long-lasting sparks by decreasing spark frequency. The decrease in spark frequency promotes intra-SR Ca diffusion from neighboring non-firing CRUs to the firing CRUs, which helps to maintain the local SR Ca concentration of the firing CRUs above a critical level to sustain firing. In this setting, decreasing RyR open probability further suppresses long-lasting sparks by weakening CICR. Since a long-lasting spark terminates via the Kramers’ escape process over a potential barrier, its duration exhibits an exponential distribution determined by the barrier height and noise strength, which is modulated differently by different ways of altering the Ca release flux strength.     

Factors shaping the confocal image of the calcium spark in cardiac muscle cells.

The interpretation of confocal line-scan images of local [Ca2+]i transients (such as Ca2+ sparks in cardiac muscle) is complicated by uncertainties in the position of the origin of the Ca2+ spark (relative to the scan line) and by the dynamics of Ca(2+)-dye interactions. An investigation of the effects of these complications modeled the release, diffusion, binding, and uptake of Ca2+ in cardiac cells (producing a theoretical Ca2+ spark) and image formation in a confocal microscope (after measurement of its point-spread function) and simulated line-scan images of a theoretical Ca2+ spark (when it was viewed from all possible positions relative to the scan line). In line-scan images, Ca2+ sparks that arose in a different optical section or with the site of origin displaced laterally from the scan line appeared attenuated, whereas their rise times slowed down only slightly. These results indicate that even if all Ca2+ sparks are perfectly identical events, except for their site of origin, there will be an apparent variation in the amplitude and other characteristics of Ca2+ sparks as measured from confocal line-scan images. The frequency distributions of the kinetic parameters (i.e., peak amplitude, rise time, fall time) of Ca2+ sparks were calculated for repetitive registration of stereotyped Ca2+ sparks in two experimental situations: 1) random position of the scan line relative to possible SR Ca(2+)-release sites and 2) fixed position of the scan line going through a set of possible SR Ca(2+)-release sites. The effects of noise were incorporated into the model, and a visibility function was proposed to account for the subjective factors that may be involved in the evaluation of Ca(2+)-spark image parameters from noisy experimental recordings. The mean value of the resulting amplitude distributions underestimates the brightness of in-focus Ca2+ sparks because large numbers of out-of-focus Ca2+ sparks are detected (as small Ca2+ sparks). The distribution of peak amplitudes may split into more than one subpopulation even when one is viewing stereotyped Ca2+ sparks because of the discrete locations of possible SR Ca(2+)-release sites in mammalian ventricular heart cells.     

SparkMaster: automated calcium spark analysis with ImageJ

Ca sparks are elementary Ca-release events from intracellular Ca stores that are observed in virtually all types of muscle. Typically, Ca sparks are measured in the line-scan mode with confocal laser-scanning microscopes, yielding two-dimensional images (distance vs. time). The manual analysis of these images is time consuming and prone to errors as well as investigator bias. Therefore, we developed SparkMaster, an automated analysis program that allows rapid and reliable spark analysis. The underlying analysis algorithm is adapted from the threshold-based standard method of spark analysis developed by Cheng et al. (Biophys J 76: 606–617, 1999) and is implemented here in the freely available image-processing software ImageJ. SparkMaster offers a graphical user interface through which all analysis parameters and output options are selected. The analysis includes general image parameters (number of detected sparks, spark frequency) and individual spark parameters (amplitude, full width at half-maximum amplitude, full duration at half-maximum amplitude, full width, full duration, time to peak, maximum steepness of spark upstroke, time constant of spark decay). We validated the algorithm using images with synthetic sparks embedded into backgrounds with different signal-to-noise ratios to determine an analysis criteria at which a high sensitivity is combined with a low frequency of false-positive detections. Finally, we applied SparkMaster to analyze experimental data of sparks measured in intact and permeabilized ventricular cardiomyocytes, permeabilized mammalian skeletal muscle, and intact smooth muscle cells. We found that SparkMaster provides a reliable, easy to use, and fast way of analyzing Ca sparks in a wide variety of experimental conditions. myocytes; sarcoplasmic reticulum; confocal microscopy     

Increasing sensitivity of Ca2+ spark detection in noisy images by application of a matched-filter object detection algorithm

Microscopic calcium (Ca²⁺) events (such as Ca²⁺ sparks) are an important area for study, as they help clarify the mechanism(s) underlying intracellular signaling. In the heart, Ca²⁺ sparks occur as a result of Ca²⁺ release from the sarcoendoplasmic reticulum, via ryanodine receptor channels. Measurement of Ca²⁺ spark properties can provide valuable information about the control of ryanodine receptor channel gating in situ, but requires high spatiotemporal resolution imaging, which produces noisy datasets that are problematic for spark detection. Automated detection algorithms may overcome visual detection bias, but missed and false-positive events can distort the distribution of measured Ca²⁺ spark properties. We present a sensitive and reliable method for the automated detection of Ca²⁺ sparks in datasets obtained using confocal line-scanning or total internal reflection fluorescence microscopy. This matched-filter detection algorithm (MFDA) employs a user-defined object, chosen to mimic Ca²⁺ spark properties, and the experimental dataset is searched for instances of the object. Detection certainty is provided by nonparametric statistical testing. The supplied codes can also refine the search object on the basis of those detected to further increase detection sensitivity. In comparison to a commonly used, intensity-thresholding algorithm, the MFDA is more sensitive and reliable, particularly at low signal/noise ratios. The MFDA can also be easily adapted to other signal-detection problems in noisy datasets.     

Life and death of a cardiac calcium spark

Calcium sparks in cardiac myocytes are brief, localized calcium releases from the sarcoplasmic reticulum (SR) believed to be caused by locally regenerative calcium-induced calcium release (CICR) via couplons, clusters of ryanodine receptors (RyRs). How such regeneration is terminated is uncertain. We performed numerical simulations of an idealized stochastic model of spark production, assuming a RyR gating scheme with only two states (open and closed). Local depletion of calcium in the SR was inevitable during a spark, and this could terminate sparks by interrupting CICR, with or without assumed modulation of RyR gating by SR lumenal calcium. Spark termination by local SR depletion was not robust: under some conditions, sparks could be greatly and variably prolonged, terminating by stochastic attrition–a phenomenon we dub “spark metastability.” Spark fluorescence rise time was not a good surrogate for the duration of calcium release. Using a highly simplified, deterministic model of the dynamics of a couplon, we show that spark metastability depends on the kinetic relationship of RyR gating and junctional SR refilling rates. The conditions for spark metastability resemble those produced by known mutations of RyR2 and CASQ2 that cause life-threatening triggered arrhythmias, and spark metastability may be mitigated by altering the kinetics of the RyR in a manner similar to the effects of drugs known to prevent those arrhythmias. The model was unable to explain the distributions of spark amplitudes and rise times seen in chemically skinned cat atrial myocytes, suggesting that such sparks may be more complex events involving heterogeneity of couplons or local propagation among sub-clusters of RyRs.     

Theoretical analysis of the Ca2+ spark amplitude distribution.

A difficulty of using confocal microscopy to study Ca2+ sparks is the uncertainty of the linescan position with respect to the source of Ca2+ release. Random placement of the linescan is expected to result in a broad distribution of measured Ca2+ spark amplitudes (a) even if all Ca2+ sparks were generated identically. Thus variations in Ca2+ spark amplitude due to positional differences between confocal linescans and Ca2+ release site are intertwined with variations due to intrinsic differences in Ca2+ release properties. To separate these two sources of variations on the Ca2+ spark amplitude, we determined the effect changes of channel current or channel open time--collectively called the source strength, alpha--had on the measured Ca2+ spark amplitude histogram, N(a). This was done by 1) simulating Ca2+ release, Ca2+ and fluo-3 diffusion, and Ca2+ binding reactions; 2) simulation of image formation of the Ca2+ spark by a confocal microscope; and 3) using a novel automatic Ca2+ spark detector. From these results we derived an integral equation relating the probability density function of source strengths, f alpha (alpha), to N(a), which takes into account random positional variations between the source and linescan. In the special, but important, case that the spatial distribution of Ca(2+)-bound fluo-3 is Gaussian, we show the following: 1) variations of Ca2+ spark amplitude due to positional or intrinsic differences can be separated, and 2) f alpha (alpha) can, in principle, be calculated from the Ca2+ spark amplitude histogram since N(a) is the sum of shifted hyperbolas, where the magnitudes of the shifts and weights depend on f alpha (alpha). In particular, if all Ca2+ sparks were generated identically, then the plot of 1/N(a) against a will be a straight line. Multiple populations of channels carrying distinct currents are revealed by discontinuities in the 1/N(a) plot. 3) Although the inverse relationship between Ca2+ spark amplitude and decay time might be used to distinguish Ca2+ sparks from different channel populations, noise can render the measured decay times meaningless for small amplitude Ca2+ sparks.     

CaMKIIδC slows [Ca]i decline in cardiac myocytes by promoting Ca sparks

Acute activation of calcium/calmodulin-dependent protein kinase (CaMKII) in permeabilized phospholamban knockout (PLN-KO) mouse myocytes phosphorylates ryanodine receptors (RyRs) and activates spontaneous local sarcoplasmic reticulum (SR) Ca release events (Ca sparks) even at constant SR Ca load. To assess how CaMKII regulates SR Ca release in intact myocytes (independent of SR Ca content changes or PLN effects), we compared Ca sparks in PLN-KO versus mice, which also have transgenic cardiac overexpression of CaMKIIδC in the PLN-KO background (KO/TG). Compared with PLN-KO mice, these KO/TG cardiomyocytes exhibited 1), increased twitch Ca transient and fractional release (both by ～35%), but unaltered SR Ca load; 2), increased resting Ca spark frequency (300%) despite a lower diastolic [Ca]i, which also slowed twitch [Ca]i decline (suggesting CaMKII-dependent RyR Ca sensitization); 3), elevated Ca spark amplitude and rate of Ca release (which might indicate that more RyR channels participate in a single spark); 4), prolonged Ca spark rise time (which implies that CaMKII either delays RyR closure or prolongs the time when openings can occur); 5), more frequent repetitive sparks at single release sites. Analysis of repetitive sparks from individual Ca release sites indicates that CaMKII enhanced RyR Ca sensitivity, but did not change the time course of SR Ca refilling. These results demonstrate that there are dramatic CaMKII-mediated effects on RyR Ca release that occur via regulation of both RyR activation and termination processes.     

Large currents generate cardiac Ca2+ sparks

Previous models of cardiac Ca2+ sparks have assumed that Ca2+ currents through the Ca2+ release units (CRUs) were approximately 1-2 pA, producing sparks with peak fluorescence ratio (F/F(0)) of approximately 2.0 and a full-width at half maximum (FWHM) of approximately 1 microm. Here, we present actual Ca2+ sparks with peak F/F(0) of >6 and a FWHM of approximately 2 microm, and a mathematical model of such sparks, the main feature of which is a much larger underlying Ca2+ current. Assuming infinite reaction rates and no endogenous buffers, we obtain a lower bound of approximately 11 pA needed to generate a Ca2+ spark with FWHM of 2 microm. Under realistic conditions, the CRU current must be approximately 20 pA to generate a 2- microm Ca2+)spark. For currents > or =5 pA, the computed spark amplitudes (F/F(0)) are large (approximately 6-12 depending on buffer model). We considered several factors that might produce sparks with FWHM approximately 2 microm without using large currents. Possible protein-dye interactions increased the FWHM slightly. Hypothetical Ca2+ "quarks" had little effect, as did blurring of sparks by the confocal microscope. A clusters of CRUs, each producing 10 pA simultaneously, can produce sparks with FWHM approximately 2 microm. We conclude that cardiac Ca2+ sparks are significantly larger in peak amplitude than previously thought, that such large Ca2+ sparks are consistent with the measured FWHM of approximately 2 microm, and that the underlying Ca2+ current is in the range of 10-20 pA.     

Ryanodine receptor allosteric coupling and the dynamics of calcium sparks

Puffs and sparks are localized intracellular Ca²⁺ elevations that arise from the cooperative activity of Ca²⁺-regulated inositol 1,4,5-trisphosphate receptors and ryanodine receptors clustered at Ca²⁺ release sites on the surface of the endoplasmic reticulum or the sarcoplasmic reticulum. While the synchronous gating of Ca²⁺-regulated Ca²⁺ channels can be mediated entirely though the buffered diffusion of intracellular Ca²⁺, interprotein allosteric interactions also contribute to the dynamics of ryanodine receptor (RyR) gating and Ca²⁺ sparks. In this article, Markov chain models of Ca²⁺ release sites are used to investigate how the statistics of Ca²⁺ spark generation and termination are related to the coupling of RyRs via local [Ca²⁺] changes and allosteric interactions. Allosteric interactions are included in a manner that promotes the synchronous gating of channels by stabilizing neighboring closed-closed and/or open-open channel pairs. When the strength of Ca²⁺-mediated channel coupling is systematically varied (e.g., by changing the Ca²⁺ buffer concentration), simulations that include synchronizing allosteric interactions often exhibit more robust Ca²⁺ sparks; however, for some Ca²⁺ coupling strengths the sparks are less robust. We find no evidence that the distribution of spark durations can be used to distinguish between allosteric interactions that stabilize closed channel pairs, open channel pairs, or both in a balanced fashion. On the other hand, the changes in spark duration, interspark interval, and frequency observed when allosteric interactions that stabilize closed channel pairs are gradually removed from simulations are qualitatively different than the changes observed when open or both closed and open channel pairs are stabilized. Thus, our simulations clarify how changes in spark statistics due to pharmacological washout of the accessory proteins mediating allosteric coupling may indicate the type of synchronizing allosteric interactions exhibited by physically coupled RyRs. We also investigate the validity of a mean-field reduction applicable to the dynamics of a ryanodine receptor cluster coupled via local [Ca²⁺] and allosteric interactions. In addition to facilitating parameter studies of the effect of allosteric coupling on spark statistics, the derivation of the mean-field model establishes the correct functional form for cooperativity factors representing the coupled gating of RyRs. This mean-field formulation is well suited for use in computationally efficient whole cell simulations of excitation-contraction coupling.     

A mathematical analysis of the generation and termination of calcium sparks

Calcium sparks are local regenerative releases of Ca(2+) from a cluster of ryanodine receptors on the sarcoplasmic reticulum. During excitation-contraction coupling in cardiac cells, Ca(2+) sparks are triggered by Ca(2+) entering the cell via the T-tubules (Ca(2+)-induced Ca(2+) release). However under conditions of calcium overload, Ca(2+) sparks can be triggered spontaneously. The exact process by which Ca(2+) sparks terminate is still an open question, although both deterministic and stochastic processes are likely to be important. In this article, asymptotic methods are used to analyze a single Ca(2+) spark model, which includes both deterministic and stochastic biophysical mechanisms. The analysis calculates both spark frequencies and spark duration distributions, and shows under what circumstances stochastic transitions are important. Additionally, a model of the coupling of the release channels via the FK-binding protein is analyzed.     

A simple numerical model of calcium spark formation and detection in cardiac myocytes.

The elementary events of excitation-contraction coupling in heart muscle are Ca2+ sparks, which arise from one or more ryanodine receptors in the sarcoplasmic reticulum (SR). Here a simple numerical model is constructed to explore Ca2+ spark formation, detection, and interpretation in cardiac myocytes. This model includes Ca2+ release, cytosolic diffusion, resequestration by SR Ca2+-ATPases, and the association and dissociation of Ca2+ with endogenous Ca2+-binding sites and a diffusible indicator dye (fluo-3). Simulations in a homogeneous, isotropic cytosol reproduce the brightness and the time course of a typical cardiac Ca2+ spark, but underestimate its spatial size (approximately 1.1 micron vs. approximately 2.0 micron). Back-calculating [Ca2+]i by assuming equilibrium with indicator fails to provide a good estimate of the free Ca2+ concentration even when using blur-free fluorescence data. A parameter sensitivity study reveals that the mobility, kinetics, and concentration of the indicator are essential determinants of the shape of Ca2+ sparks, whereas the stationary buffers and pumps are less influential. Using a geometrically more complex version of the model, we show that the asymmetric shape of Ca2+ sparks is better explained by anisotropic diffusion of Ca2+ ions and indicator dye rather than by subsarcomeric inhomogeneities of the Ca2+ buffer and transport system. In addition, we examine the contribution of off-center confocal sampling to the variance of spark statistics.     

Automated Detection and Analysis of Ca Sparks in x-y Image Stacks Using a Thresholding Algorithm Implemented within the Open-Source Image Analysis Platform ImageJ

Previous studies have used analysis of Ca²⁺ sparks extensively to investigate both normal and pathological Ca²⁺ regulation in cardiac myocytes. The great majority of these studies used line-scan confocal imaging. In part, this is because the development of open-source software for automatic detection of Ca²⁺ sparks in line-scan images has greatly simplified data analysis. A disadvantage of line-scan imaging is that data are collected from a single row of pixels, representing only a small fraction of the cell, and in many instances x-y confocal imaging is preferable. However, the limited availability of software for Ca²⁺ spark analysis in two-dimensional x-y image stacks presents an obstacle to its wider application. This study describes the development and characterization of software to enable automatic detection and analysis of Ca²⁺ sparks within x-y image stacks, implemented as a plugin within the open-source image analysis platform ImageJ. The program includes methods to enable precise identification of cells within confocal fluorescence images, compensation for changes in background fluorescence, and options that allow exclusion of events based on spatial characteristics.     

Automated Detection and Analysis of Ca2+ Sparks in x-y Image Stacks Using a Thresholding Algorithm Implemented within the Open-Source Image Analysis Platform ImageJ

Previous studies have used analysis of Ca²⁺ sparks extensively to investigate both normal and pathological Ca²⁺ regulation in cardiac myocytes. The great majority of these studies used line-scan confocal imaging. In part, this is because the development of open-source software for automatic detection of Ca²⁺ sparks in line-scan images has greatly simplified data analysis. A disadvantage of line-scan imaging is that data are collected from a single row of pixels, representing only a small fraction of the cell, and in many instances x-y confocal imaging is preferable. However, the limited availability of software for Ca²⁺ spark analysis in two-dimensional x-y image stacks presents an obstacle to its wider application. This study describes the development and characterization of software to enable automatic detection and analysis of Ca²⁺ sparks within x-y image stacks, implemented as a plugin within the open-source image analysis platform ImageJ. The program includes methods to enable precise identification of cells within confocal fluorescence images, compensation for changes in background fluorescence, and options that allow exclusion of events based on spatial characteristics.     

A close association of RyRs with highly dense clusters of Ca2+-activated Cl- channels underlies the activation of STICs by Ca2+ sparks in mouse airway smooth muscle

Ca(2+) sparks are highly localized, transient releases of Ca(2+) from sarcoplasmic reticulum through ryanodine receptors (RyRs). In smooth muscle, Ca(2+) sparks trigger spontaneous transient outward currents (STOCs) by opening nearby clusters of large-conductance Ca(2+)-activated K(+) channels, and also gate Ca(2+)-activated Cl(-) (Cl((Ca))) channels to induce spontaneous transient inward currents (STICs). While the molecular mechanisms underlying the activation of STOCs by Ca(2+) sparks is well understood, little information is available on how Ca(2+) sparks activate STICs. In the present study, we investigated the spatial organization of RyRs and Cl((Ca)) channels in spark sites in airway myocytes from mouse. Ca(2+) sparks and STICs were simultaneously recorded, respectively, with high-speed, widefield digital microscopy and whole-cell patch-clamp. An image-based approach was applied to measure the Ca(2+) current underlying a Ca(2+) spark (I(Ca(spark))), with an appropriate correction for endogenous fixed Ca(2+) buffer, which was characterized by flash photolysis of NPEGTA. We found that I(Ca(spark)) rises to a peak in 9 ms and decays with a single exponential with a time constant of 12 ms, suggesting that Ca(2+) sparks result from the nonsimultaneous opening and closure of multiple RyRs. The onset of the STIC lags the onset of the I(Ca(spark)) by less than 3 ms, and its rising phase matches the duration of the I(Ca(spark)). We further determined that Cl((Ca)) channels on average are exposed to a [Ca(2+)] of 2.4 microM or greater during Ca(2+) sparks. The area of the plasma membrane reaching this level is <600 nm in radius, as revealed by the spatiotemporal profile of [Ca(2+)] produced by a reaction-diffusion simulation with measured I(Ca(spark)). Finally we estimated that the number of Cl((Ca)) channels localized in Ca(2+) spark sites could account for all the Cl((Ca)) channels in the entire cell. Taken together these results lead us to propose a model in which RyRs and Cl((Ca)) channels in Ca(2+) spark sites localize near to each other, and, moreover, Cl((Ca)) channels concentrate in an area with a radius of approximately 600 nm, where their density reaches as high as 300 channels/microm(2). This model reveals that Cl((Ca)) channels are tightly controlled by Ca(2+) sparks via local Ca(2+) signaling.     

Termination of cardiac Ca(2+) sparks: an investigative mathematical model of calcium-induced calcium release

A Ca(2+) spark arises when a cluster of sarcoplasmic reticulum (SR) channels (ryanodine receptors or RyRs) opens to release calcium in a locally regenerative manner. Normally triggered by Ca(2+) influx across the sarcolemmal or transverse tubule membrane neighboring the cluster, the Ca(2+) spark has been shown to be the elementary Ca(2+) signaling event of excitation-contraction coupling in heart muscle. However, the question of how the Ca(2+) spark terminates remains a central, unresolved issue. Here we present a new model, "sticky cluster," of SR Ca(2+) release that simulates Ca(2+) spark behavior and enables robust Ca(2+) spark termination. Two newly documented features of RyR behavior have been incorporated in this otherwise simple model: "coupled gating" and an opening rate that depends on SR lumenal [Ca(2+)]. Using a Monte Carlo method, local Ca(2+)-induced Ca(2+) release from clusters containing between 10 and 100 RyRs is modeled. After release is triggered, Ca(2+) flux from RyRs diffuses into the cytosol and binds to intracellular buffers and the fluorescent Ca(2+) indicator fluo-3 to produce the model Ca(2+) spark. Ca(2+) sparks generated by the sticky cluster model resemble those observed experimentally, and Ca(2+) spark duration and amplitude are largely insensitive to the number of RyRs in a cluster. As expected from heart cell investigation, the spontaneous Ca(2+) spark rate in the model increases with elevated cytosolic or SR lumenal [Ca(2+)]. Furthermore, reduction of RyR coupling leads to prolonged model Ca(2+) sparks just as treatment with FK506 lengthens Ca(2+) sparks in heart cells. This new model of Ca(2+) spark behavior provides a "proof of principle" test of a new hypothesis for Ca(2+) spark termination and reproduces critical features of Ca(2+) sparks observed experimentally.     

Time Course of Individual Ca2+ Sparks in Frog Skeletal Muscle Recorded at High Time Resolution

Discrete Ca²⁺ release events (Ca²⁺ “sparks”) were recorded in cut segments of single frog skeletal muscle fibers using a video-rate laser-scanning confocal microscope operating in line-scan mode (63 μs per line). Fibers loaded with the Ca²⁺ indicator fluo-3 were voltage clamped at a holding potential of 0 mV, briefly reprimed at −90 mV, and then strongly depolarized with a large test pulse to activate any reprimed voltage sensors. Using this high time resolution system, it was possible to record individual Ca²⁺ sparks at ∼30-fold higher time resolution than previously attained. The resulting new experimental data provides a means of characterizing the time course of fluorescence during the brief (a few milliseconds) rising phase of a spark, which was not possible with the previously used 1.5–2 ms per line confocal systems. Analysis of the time course of individual identified events indicates that fluorescence begins to rise rather abruptly at the start of the spark, continues to rise at a slightly decreasing rate to a relatively sharp peak, and then declines along a quasi-exponential time course. The mean rise time of 198 sparks was 4.7 ± 0.1 ms, and there was no correlation between rise time and peak amplitude. Average sparks constructed by temporally and spatially superimposing and summing groups of individual sparks having similar rise times gave a lower noise representation of the sparks, consistent with the time course of individual events. In theory, the rising phase of a spark provides a lower bound estimation of the time that Ca²⁺ ions are being released by the sarcoplasmic reticulum Ca²⁺ channel(s) generating the spark. The observed time course of fluorescence suggests that the Ca²⁺ release underlying a spark could continue at a fairly constant rate throughout the rising phase of the spark, and then stop rather abruptly at the time of the peak.     

ET-1 activates Ca2+ sparks in PASMC: local Ca2+ signaling between inositol trisphosphate and ryanodine receptors

Ca+ sparks originating from ryanodine receptors (RyRs) are known to cause membrane hyperpolarization and vasorelaxation in systemic arterial myocytes. By contrast, we have found that Ca2+ sparks of pulmonary arterial smooth muscle cells (PASMCs) are associated with membrane depolarization and activated by endothelin-1 (ET-1), a potent vasoconstrictor that mediates/modulates acute and chronic hypoxic pulmonary vasoconstriction. In this study, we characterized the effects of ET-1 on the physical properties of Ca2+ sparks and probed the signal transduction mechanism for spark activation in rat intralobar PASMCs. Application of ET-1 at 0.1-10 nM caused concentration-dependent increases in frequency, duration, and amplitude of Ca2+ sparks. The ET-1-induced increase in spark frequency was inhibited by BQ-123, an ETA-receptor antagonist; by U-73122, a PLC inhibitor; and by xestospongin C and 2-aminoethyl diphenylborate, antagonists of inositol trisphosphate (IP3) receptors (IP3Rs). However, it was unrelated to sarcoplasmic reticulum Ca2+ content, activation of L-type Ca2+ channels, PKC, or cADP ribose. Photorelease of caged-IP3 indicated that Ca2+ release from IP3R could cross-activate RyRs to generate Ca2+ sparks. Immunocytochemistry showed that the distributions of IP3Rs and RyRs were similar in PASMCs. Moreover, inhibition of Ca2+ sparks with ryanodine caused a significant rightward shift in the ET-1 concentration-tension relationship in pulmonary arteries. These results suggest that ET-1 activation of Ca2+ sparks is mediated via the ETA receptor-PLC-IP3 pathway and local Ca2+ cross-signaling between IP3Rs and RyRs; in addition, this novel signaling mechanism contributes significantly to the ET-1-induced vasoconstriction in pulmonary arteries.