Analysis of puff dynamics in oocytes: interdependence of puff amplitude and interpuff interval

Puffs are localized Ca(2+) signals that arise in oocytes in response to inositol 1,4,5-trisphosphate (IP(3)). They are analogous to the sparks of myocytes and are believed to be the result of the liberation of Ca(2+) from the endoplasmic reticulum through the coordinated opening of IP(3) receptor/channels clustered at a functional release site. In this article, we analyze sequences of puffs that occur at the same site to help elucidate the mechanisms underlying puff dynamics. In particular, we show a dependence of the interpuff time on the amplitude of the preceding puff, and of the amplitude of the following puff on the preceding interval. These relationships can be accounted for by an inhibitory role of the Ca(2+) that is liberated during puffs. We construct a stochastic model for a cluster of IP(3) receptor/channels that quantitatively replicates the observed behavior, and we determine that the characteristic time for a channel to escape from the inhibitory state is of the order of seconds.     

The probability of triggering calcium puffs is linearly related to the number of inositol trisphosphate receptors in a cluster

Puffs are local Ca(2+) signals that arise by Ca(2+) liberation from the endoplasmic reticulum through concerted opening of tightly clustered inositol trisphosphate receptor/channels (IP(3)R). They serve both local signaling functions and trigger global Ca(2+) waves. The numbers of functional IP(3)R within clusters differ appreciably between different puff sites, and we investigated how the probability of puff occurrence varies with cluster size. We imaged puffs in SH-SY5Y cells using total internal fluorescence microscopy, and estimated cluster sizes from the magnitude of the largest puff observed at each site relative to the signal from a single channel. We find that the initial triggering rate of puffs following photorelease of IP(3), and the average frequency of subsequent repetitive puffs, vary about linearly with cluster size. These data accord well with stochastic simulations in which opening of any individual IP(3)R channel within a cluster triggers a puff via Ca(2+)-induced Ca(2+) release. An important consequence is that the signaling power of a puff site (average amount of Ca(2+) released per puff × puff frequency) varies about the square of cluster size, implying that large clusters contribute disproportionately to cellular signaling and, because of their higher puff frequency, preferentially act as pacemakers to initiate Ca(2+) waves.     

The dynamics of luminal depletion and the stochastic gating of Ca2+-activated Ca2+ channels and release sites

Single channel models of intracellular calcium (Ca(2+)) channels such as the 1,4,5-trisphosphate receptor and ryanodine receptor often assume that Ca(2+)-dependent transitions are mediated by constant background cytosolic [Ca(2+)]. This assumption neglects the fact that Ca(2+) released by open channels may influence subsequent gating through the processes of Ca(2+)-activation or inactivation. Similarly, the influence of the dynamics of luminal depletion on the stochastic gating of intracellular Ca(2+) channels is often neglected, in spite of the fact that the sarco/endoplasmic reticulum [Ca(2+)] near the luminal face of intracellular Ca(2+) channels influences the driving force for Ca(2+), the rate of Ca(2+) release, and the magnitude and time course of the consequent increase in cytosolic domain [Ca(2+)]. Here we analyze how the steady-state open probability of several minimal Ca(2+)-regulated Ca(2+) channel models depends on the conductance of the channel and the time constants for the relaxation of elevated cytosolic [Ca(2+)] and depleted luminal [Ca(2+)] to the bulk [Ca(2+)] of both compartments. Our approach includes Monte Carlo simulation as well as numerical solution of a system of advection-reaction equations for the multivariate probability density of elevated cytosolic [Ca(2+)] and depleted luminal [Ca(2+)] conditioned on each state of the stochastically gating channel. Both methods are subsequently used to study the role of luminal depletion in the dynamics of Ca(2+) puff/spark termination in release sites composed of Ca(2+) channels that are activated, but not inactivated, by cytosolic Ca(2+). The probability density approach shows that such minimal Ca(2+) release site models may exhibit puff/spark-like dynamics in either of two distinct parameter regimes. In one case, puffs/spark termination is due to the process of stochastic attrition and facilitated by rapid Ca(2+) domain collapse [cf. DeRemigio, H., Smith, G., 2005. The dynamics of stochastic attrition viewed as an absorption time on a terminating Markov chain. Cell Calcium 38, 73-86]. In the second case, puff/spark termination is promoted by the local depletion of luminal Ca(2+).     

Stochastic properties of Ca(2+) release of inositol 1,4,5-trisphosphate receptor clusters

Intracellular Ca(2+) release is controlled by inositol 1,4,5-trisphosphate (IP(3)) receptors or ryanodine receptors. These receptors are typically distributed in clusters with several or tens of channels. The random opening and closing of these channels introduces stochasticity into the elementary calcium release mechanism. Stochastic release events have been experimentally observed in a variety of cell types and have been termed sparks and puffs. We put forward a stochastic version of the Li-Rinzel model (the deactivation binding process is described by a Markovian scheme) and a computationally more efficient Langevin approach to model the stochastic Ca(2+) oscillation of single clusters. Statistical properties such as Ca(2+) puff amplitudes, lifetimes, and interpuff intervals are studied with both models and compared with experimental observations. For clusters with tens of channels, a simply decaying amplitude distribution is typically observed at low IP(3) concentration, while a single peak distribution appears at high IP(3) concentration.     

Reduction of calcium release site models via moment fitting of phase-type distributions

Models of calcium (Ca(2 +)) release sites derived from continuous-time Markov chain (CTMC) models of intracellular Ca(2 +) channels exhibit collective gating reminiscent of the experimentally observed phenomenon of Ca(2 +) puffs and sparks. In order to overcome the state-space explosion that occurs in compositionally defined Ca(2 +) release site models, we have implemented an automated procedure for model reduction that replaces aggregated states of the full release site model with much simpler CTMCs that have similar within-group phase-type sojourn times and inter-group transitions. Error analysis based on comparison of full and reduced models validates the method when applied to release site models composed of 20 three-state channels that are both activated and inactivated by Ca(2 +). Although inspired by existing techniques for fitting moments of phase-type distributions, the automated reduction method for compositional Ca(2 +) release site models is unique in several respects and novel in this biophysical context.     

Role of elementary Ca(2+) puffs in generating repetitive Ca(2+) oscillations

Inositol (1,4,5)-trisphosphate (IP(3)) liberates intracellular Ca(2+) both as localized 'puffs' and as repetitive waves that encode information in a frequency-dependent manner. Using video-rate confocal imaging, together with photorelease of IP(3) in Xenopus oocytes, we investigated the roles of puffs in determining the periodicity of global Ca(2+) waves. Wave frequency is not delimited solely by cyclical recovery of the cell's ability to support wave propagation, but further involves sensitization of Ca(2+)-induced Ca(2+) release by progressive increases in puff frequency and amplitude at numerous sites during the interwave period, and accumulation of pacemaker Ca(2+), allowing a puff at a 'focal' site to trigger a subsequent wave. These specific 'focal' sites, distinguished by their higher sensitivity to IP(3) and close apposition to neighboring puff sites, preferentially entrain both the temporal frequency and spatial directionality of Ca(2+) waves. Although summation of activity from many stochastic puff sites promotes the generation of regularly periodic global Ca(2+) signals, the properties of individual Ca(2+) puffs control the kinetics of Ca(2+) spiking and the (higher) frequency of subcellular spikes in their local microdomain.     

Markov chain models of coupled calcium channels: Kronecker representations and iterative solution methods

Mathematical models of calcium release sites derived from Markov chain models of intracellular calcium channels exhibit collective gating reminiscent of the experimentally observed phenomenon of stochastic calcium excitability (i.e., calcium puffs and sparks). Calcium release site models are stochastic automata networks that involve many functional transitions, that is, the transition probabilities of each channel depend on the local calcium concentration and thus the state of the other channels. We present a Kronecker-structured representation for calcium release site models and perform benchmark stationary distribution calculations using both exact and approximate iterative numerical solution techniques that leverage this structure. When it is possible to obtain an exact solution, response measures such as the number of channels in a particular state converge more quickly using the iterative numerical methods than occupation measures calculated via Monte Carlo simulation. In particular, multi-level methods provide excellent convergence with modest additional memory requirements for the Kronecker representation of calcium release site models. When an exact solution is not feasible, iterative approximate methods based on the power method may be used, with performance similar to Monte Carlo estimates. This suggests approximate methods with multi-level iterative engines as a promising avenue of future research for large-scale calcium release site models.     

Modeling the statistics of elementary calcium release events

Elementary Ca(2+) signals, such as "Ca(2+) puffs", which arise from the release of Ca(2+) from endoplasmic reticulum through small clusters of inositol 1,4,5-trisphosphate receptors, are the building blocks for intracellular Ca(2+) signaling. The small number of release channels involved during a Ca(2+) puff renders the puffs stochastic, with distributed amplitudes, durations, and frequency, well characterized experimentally. We present a stochastic model that accurately describes simultaneously the statistical properties of the duration, amplitudes, frequencies, and spatial spread with a single set of parameters.     

Calcium-dependent inactivation and the dynamics of calcium puffs and sparks

Localized intracellular Ca²⁺ elevations known as puffs and sparks arise from the cooperative activity of inositol 1,4,5-trisphosphate receptor Ca²⁺ channels (IP₃Rs) and ryanodine receptor Ca²⁺ channels (RyRs) clustered at Ca²⁺ release sites on the surface of the endoplasmic reticulum or sarcoplasmic reticulum. When Markov chain models of these intracellular Ca²⁺-regulated Ca²⁺ channels are coupled via a mathematical representation of a Ca²⁺ microdomain, simulated Ca²⁺ release sites may exhibit the phenomenon of “stochastic Ca²⁺ excitability” reminiscent of Ca²⁺ puffs and sparks where channels open and close in a concerted fashion. To clarify the role of Ca²⁺ inactivation of IP₃Rs and RyRs in the dynamics of puffs and sparks, we formulate and analyze Markov chain models of Ca²⁺ release sites composed of 10–40 three-state intracellular Ca²⁺ channels that are inactivated as well as activated by Ca²⁺. We study how the statistics of simulated puffs and sparks depend on the kinetics and dissociation constant of Ca²⁺ inactivation and find that puffs and sparks are often less sensitive to variations in the number of channels at release sites and strength of coupling via local [Ca²⁺] when the average fraction of inactivated channels is significant. Interestingly, we observe that the single channel kinetics of Ca²⁺ inactivation influences the thermodynamic entropy production rate of Markov chain models of puffs and sparks. While excessively fast Ca²⁺ inactivation can preclude puffs and sparks, moderately fast Ca²⁺ inactivation often leads to time-irreversible puffs and sparks whose termination is facilitated by the recruitment of inactivated channels throughout the duration of the puff/spark event. On the other hand, Ca²⁺ inactivation may be an important negative feedback mechanism even when its time constant is much greater than the duration of puffs and sparks. In fact, slow Ca²⁺ inactivation can lead to release sites with a substantial fraction of inactivated channels that exhibit puffs and sparks that are nearly time-reversible and terminate without additional recruitment of inactivated channels.     

Thermodynamically irreversible gating of ryanodine receptors in situ revealed by stereotyped duration of release in Ca(2+) sparks

For a single or a group of Markov channels gating reversibly, distributions of open and closed times should be the sum of positively weighted decaying exponentials. Violation of this microscopic reversibility has been demonstrated previously on a number of occasions at the single channel level, and has been attributed to possible channel coupling to external sources of free energy. Here we show that distribution of durations of Ca(2+) release underlying Ca(2+) sparks in intact cardiac myocytes exhibits a prominent mode at approximately 8 ms. Analysis of the cycle time for repetitive sparks at hyperactive sites revealed no intervals briefer than approximately 35 ms and a mode at approximately 90 ms. These results indicate that, regardless of whether Ca(2+) sparks are single-channel or multi-channel in origin, they are generated by thermodynamically irreversible stochastic processes. In contrast, data from planar lipid bilayer experiments were consistent with reversible gating of RyR under asymmetric cis (4 microM) and trans Ca(2+) (10 mM), suggesting that the irreversibility for Ca(2+) spark genesis may reside at a supramolecular level. Modeling suggests that Ca(2+)-induced Ca(2+) release among adjacent RyRs may couple the external energy derived from Ca(2+) gradients across the SR to RyR gating in situ, and drive the irreversible generation of Ca(2+) sparks.     

Invited review: mechanisms of calcium handling in smooth muscles.

The concentration of cytoplasmic Ca(2+) regulates the contractile state of smooth muscle cells and tissues. Elevations in global cytoplasmic Ca(2+) resulting in contraction are accomplished by Ca(2+) entry and release from intracellular stores. Pathways for Ca(2+) entry include dihydropyridine-sensitive and -insensitive Ca(2+) channels and receptor and store-operated nonselective channels permeable to Ca(2+). Intracellular release from the sarcoplasmic reticulum (SR) is accomplished by ryanodine and inositol trisphosphate receptors. The impact of Ca(2+) entry and release on cytoplasmic concentration is modulated by Ca(2+) reuptake into the SR, uptake into mitochondria, and extrusion into the extracellular solution. Highly localized Ca(2+) transients (i.e., sparks and puffs) regulate ionic conductances in the plasma membrane, which can provide feedback to cell excitability and affect Ca(2+) entry. This short review describes the major transport mechanisms and compartments that are utilized for Ca(2+) handling in smooth muscles.     

The number and spatial distribution of IP3 receptors underlying calcium puffs in Xenopus oocytes

Calcium puffs are local Ca(2+) release events that arise from a cluster of inositol 1,4,5-trisphosphate receptor channels (IP(3)Rs) and serve as a basic "building block" from which global Ca(2+) waves are generated. Important questions remain as to the number of IP(3)Rs that open during a puff, their spatial distribution within a cluster, and how much Ca(2+) current flows through each channel. The recent discovery of "trigger" events-small Ca(2+) signals that immediately precede puffs and are interpreted to arise through opening of single IP(3)R channels-now provides a useful yardstick by which to calibrate the Ca(2+) flux underlying puffs. Here, we describe a deterministic numerical model to simulate puffs and trigger events. Based on confocal linescan imaging in Xenopus oocytes, we simulated Ca(2+) release in two sequential stages; representing the trigger by the opening of a single IP(3)R in the center of a cluster for 12 ms, followed by the concerted opening of some number of IP(3)Rs for 19 ms, representing the rising phase of the puff. The diffusion of Ca(2+) and Ca(2+)-bound indicator dye were modeled in a three-dimensional cytosolic volume in the presence of immobile and mobile Ca(2+) buffers, and were used to predict the observed fluorescence signal after blurring by the microscope point-spread function. Optimal correspondence with experimental measurements of puff spatial width and puff/trigger amplitude ratio was obtained assuming that puffs arise from the synchronous opening of 25-35 IP(3)Rs, each carrying a Ca(2+) current of approximately 0.4 pA, with the channels distributed through a cluster 300-800 nm in diameter.     

Factors Determining the Recruitment of Inositol Trisphosphate Receptor Channels During Calcium Puffs

Puffs are localized, transient elevations in cytosolic Ca²⁺ that serve both as the building blocks of global cellular Ca²⁺ signals and as local signals in their own right. They arise from clustered inositol 1,4,5-trisphosphate receptor/channels (IP₃Rs), whose openings are coordinated by Ca²⁺-induced Ca²⁺ release (CICR). We utilized total internal reflection fluorescence imaging of Ca²⁺ signals in neuroblastoma cells with single-channel resolution to elucidate the mechanisms determining the triggering, amplitudes, kinetics, and spatial spread of puffs. We find that any given channel in a cluster has a mean probability of ∼66% of opening following opening of an initial “trigger” channel, and the probability of puff triggering thus increases steeply with increasing number of channels in a cluster (cluster size). Mean puff amplitudes scale with cluster size, but individual amplitudes vary widely, even at sites of similar cluster size, displaying similar proportions of events involving any given number of the channels in the cluster. Stochastic variation in numbers of Ca²⁺-inhibited IP₃Rs likely contributes to the variability of amplitudes of repeated puffs at a site but the amplitudes of successive puffs were uncorrelated, even though we observed statistical correlations between interpuff intervals and puff amplitudes. Initial puffs evoked following photorelease of IP₃—which would not be subject to earlier Ca²⁺-inhibition—also showed wide variability, indicating that mechanisms such as stochastic variation in IP₃ binding and channel recruitment by CICR further determine puff amplitudes. The mean termination time of puffs lengthened with increasing puff amplitude size, consistent with independent closings of channels after a given mean open time, but we found no correlation of termination time with cluster size independent of puff amplitude. The spatial extent of puffs increased with their amplitude, and puffs of similar size were of similar width, independent of cluster size.     

Factors Determining the Recruitment of Inositol Trisphosphate Receptor Channels During Calcium Puffs

Puffs are localized, transient elevations in cytosolic Ca²⁺ that serve both as the building blocks of global cellular Ca²⁺ signals and as local signals in their own right. They arise from clustered inositol 1,4,5-trisphosphate receptor/channels (IP₃Rs), whose openings are coordinated by Ca²⁺-induced Ca²⁺ release (CICR). We utilized total internal reflection fluorescence imaging of Ca²⁺ signals in neuroblastoma cells with single-channel resolution to elucidate the mechanisms determining the triggering, amplitudes, kinetics, and spatial spread of puffs. We find that any given channel in a cluster has a mean probability of ∼66% of opening following opening of an initial “trigger” channel, and the probability of puff triggering thus increases steeply with increasing number of channels in a cluster (cluster size). Mean puff amplitudes scale with cluster size, but individual amplitudes vary widely, even at sites of similar cluster size, displaying similar proportions of events involving any given number of the channels in the cluster. Stochastic variation in numbers of Ca²⁺-inhibited IP₃Rs likely contributes to the variability of amplitudes of repeated puffs at a site but the amplitudes of successive puffs were uncorrelated, even though we observed statistical correlations between interpuff intervals and puff amplitudes. Initial puffs evoked following photorelease of IP₃—which would not be subject to earlier Ca²⁺-inhibition—also showed wide variability, indicating that mechanisms such as stochastic variation in IP₃ binding and channel recruitment by CICR further determine puff amplitudes. The mean termination time of puffs lengthened with increasing puff amplitude size, consistent with independent closings of channels after a given mean open time, but we found no correlation of termination time with cluster size independent of puff amplitude. The spatial extent of puffs increased with their amplitude, and puffs of similar size were of similar width, independent of cluster size.     

Inositol (1,4,5)-trisphosphate receptor microarchitecture shapes Ca2+ puff kinetics

Inositol (1,4,5)-trisphosphate receptors (IP(3)Rs) release intracellular Ca(2+) as localized Ca(2+) signals (Ca(2+) puffs) that represent the activity of small numbers of clustered IP(3)Rs spaced throughout the endoplasmic reticulum. Although much emphasis has been placed on estimating the number of active Ca(2+) release channels supporting Ca(2+) puffs, less attention has been placed on understanding the role of cluster microarchitecture. This is important as recent data underscores the dynamic nature of IP(3)R transitions between heterogeneous cellular architectures and the differential behavior of IP(3)Rs socialized into clusters. Here, we applied a high-resolution model incorporating stochastically gating IP(3)Rs within a three-dimensional cytoplasmic space to demonstrate: 1), Ca(2+) puffs are supported by a broad range of clustered IP(3)R microarchitectures; 2), cluster ultrastructure shapes Ca(2+) puff characteristics; and 3), loosely corralled IP(3)R clusters (>200 nm interchannel separation) fail to coordinate Ca(2+) puffs, owing to inefficient triggering and impaired coupling due to reduced Ca(2+)-induced Ca(2+) release microwave velocity (<10 nm/s) throughout the channel array. Dynamic microarchitectural considerations may therefore influence Ca(2+) puff occurrence/properties in intact cells, contrasting with a more minimal role for channel number over the same simulated conditions in shaping local Ca(2+) dynamics.     

Activation and propagation of Ca(2+) release during excitation-contraction coupling in atrial myocytes

Fast two-dimensional confocal microscopy and the Ca(2+) indicator fluo-4 were used to study excitation-contraction (E-C) coupling in cat atrial myocytes which lack transverse tubules and contain both subsarcolemmal junctional (j-SR) and central nonjunctional (nj-SR) sarcoplasmic reticulum. Action potentials elicited by field stimulation induced transient increases of intracellular Ca(2+) concentration ([Ca(2+)](i)) that were highly inhomogeneous. Increases started at distinct subsarcolemmal release sites spaced approximately 2 microm apart. The amplitude and the latency of Ca(2+) release from these sites varied from beat to beat. Subsarcolemmal release fused to build a peripheral ring of elevated [Ca(2+)](i), which actively propagated to the center of the cells via Ca(2+)-induced Ca(2+) release. Resting myocytes exhibited spontaneous Ca(2+) release events, including Ca(2+) sparks and local (microscopic) or global (macroscopic) [Ca(2+)](i) waves. The microscopic [Ca(2+)](i) waves propagated in a saltatory fashion along the sarcolemma ("coupled" Ca(2+) sparks) revealing the sequential activation of Ca(2+) release sites of the j-SR. Moreover, during global [Ca(2+)](i) waves, Ca(2+) release was evident from individual nj-SR sites. Ca(2+) release sites were arranged in a regular three-dimensional grid as deduced from the functional data and shown by immunostaining of ryanodine receptor Ca(2+) release channels. The longitudinal and transverse distances between individual Ca(2+) release sites were both approximately 2 microm. Furthermore, electron microscopy revealed a continuous sarcotubular network and one peripheral coupling of j-SR with the sarcolemma per sarcomere. The results demonstrate directly that, in cat atrial myocytes, the action potential-induced whole-cell [Ca(2+)](i) transient is the spatio-temporal summation of Ca(2+) release from subsarcolemmal and central sites. First, j-SR sites are activated in a stochastic fashion by the opening of voltage-dependent sarcolemmal Ca(2+) channels. Subsequently, nj-SR sites are activated by Ca(2+)-induced Ca(2+) release propagating from the periphery.     

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.     

Voltage dependence of the coupling of Ca(2+) sparks to BK(Ca) channels in urinary bladder smooth muscle

Large-conductance Ca(2+)-dependent K(+) (BK(Ca)) channels play a critical role in regulating urinary bladder smooth muscle (UBSM) excitability and contractility. Measurements of BK(Ca) currents and intracellular Ca(2+) revealed that BK(Ca) currents are activated by Ca(2+) release events (Ca(2+) sparks) from ryanodine receptors (RyRs) in the sarcoplasmic reticulum. The goals of this project were to characterize Ca(2+) sparks and BK(Ca) currents and to determine the voltage dependence of the coupling of RyRs (Ca(2+) sparks) to BK(Ca) channels in UBSM. Ca(2+) sparks in UBSM had properties similar to those described in arterial smooth muscle. Most Ca(2+) sparks caused BK(Ca) currents at all voltages tested, consistent with the BK(Ca) channels sensing approximately 10 microM Ca(2+). Membrane potential depolarization from -50 to -20 mV increased Ca(2+) spark and BK(Ca) current frequency threefold. However, membrane depolarization over this range had a differential effect on spark and current amplitude, with Ca(2+) spark amplitude increasing by only 30% and BK(Ca) current amplitude increasing 16-fold. A major component of the amplitude modulation of spark-activated BK(Ca) current was quantitatively explained by the known voltage dependence of the Ca(2+) sensitivity of BK(Ca) channels. We, therefore, propose that membrane potential, or any other agent that modulates the Ca(2+) sensitivity of BK(Ca) channels, profoundly alters the coupling strength of Ca(2+) sparks to BK(Ca) channels.     

Ca(2+) spark sites in smooth muscle cells are numerous and differ in number of ryanodine receptors, large-conductance K(+) channels, and coupling ratio between them

Ca(2+) sparks are highly localized Ca(2+) transients caused by Ca(2+) release from sarcoplasmic reticulum through ryanodine receptors (RyR). In smooth muscle, Ca(2+) sparks activate nearby large-conductance, Ca(2+)-sensitive K(+) (BK) channels to generate spontaneous transient outward currents (STOC). The properties of individual sites that give rise to Ca(2+) sparks have not been examined systematically. We have characterized individual sites in amphibian gastric smooth muscle cells with simultaneous high-speed imaging of Ca(2+) sparks using wide-field digital microscopy and patch-clamp recording of STOC in whole cell mode. We used a signal mass approach to measure the total Ca(2+) released at a site and to estimate the Ca(2+) current flowing through RyR [I(Ca(spark))]. The variance between spark sites was significantly greater than the intrasite variance for the following parameters: Ca(2+) signal mass, I(Ca(spark)), STOC amplitude, and 5-ms isochronic STOC amplitude. Sites that failed to generate STOC did so consistently, while those at the remaining sites generated STOC without failure, allowing the sites to be divided into STOC-generating and STOC-less sites. We also determined the average number of spark sites, which was 42/cell at a minimum and more likely on the order of at least 400/cell. We conclude that 1) spark sites differ in the number of RyR, BK channels, and coupling ratio of RyR-BK channels, and 2) there are numerous Ca(2+) spark-generating sites in smooth muscle cells. The implications of these findings for the organization of the spark microdomain are explored.