Correspondence to: Eduardo Ríos, Department of Molecular Biophysics and Physiology, Rush University, 1750 W. Harrison Street, Chicago, IL 60612., erios{at}rush.edu (E-mail), Fax: 312-942-8711; (fax)
![]() |
Abstract |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
An algorithm for the calculation of Ca2+ release flux underlying Ca2+ sparks (Blatter, L.A., J. Hüser, and E. Ríos. 1997. Proc. Natl. Acad. Sci. USA. 94:41764181) was modified and applied to sparks obtained by confocal microscopy in single frog skeletal muscle fibers, which were voltage clamped in a two-Vaseline gap chamber or permeabilized and immersed in fluo-3containing internal solution. The performance of the algorithm was characterized on sparks obtained by simulation of fluorescence due to release of Ca2+ from a spherical source, in a homogeneous three-dimensional space that contained components representing cytoplasmic molecules and Ca2+ removal processes. Total release current, as well as source diameter and noise level, was varied in the simulations. Derived release flux or current, calculated by volume integration of the derived flux density, estimated quite closely the current used in the simulation, while full width at half magnitude of the derived release flux was a good monitor of source size only at diameters >0.7 µm. On an average of 157 sparks of amplitude >2 U resting fluorescence, located automatically in a representative voltage clamp experiment, the algorithm reported a release current of 16.9 pA, coming from a source of 0.5 µm, with an open time of 6.3 ms. Fewer sparks were obtained in permeabilized fibers, so that the algorithm had to be applied to individual sparks or averages of few events, which degraded its performance in comparable tests. The average current reported for 19 large sparks obtained in permeabilized fibers was 14.4 pA. A minimum estimate, derived from the rate of change of dye-bound Ca2+ concentration, was 8 pA. Such a current would require simultaneous opening of between 8 and 60 release channels with unitary Ca2+ currents of the level recorded in bilayer experiments. Real sparks differ from simulated ones mainly in having greater width. Correspondingly, the algorithm reported greater spatial extent of the source for real sparks. This may again indicate a multichannel origin of sparks, or could reflect limitations in spatial resolution.
Key Words: sarcoplasmic reticulum, excitationcontraction coupling, confocal microscopy, computer algorithm
![]() |
Introduction |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Excitationcontraction coupling in striated muscle takes place by the release of stored calcium in response to depolarization of the sarcolemma, via ryanodine receptor channels of the sarcoplasmic reticulum. When studied in frog skeletal muscle by confocal microscopic imaging of the fluorescence of fluo-3, this release is seen to occur largely in the form of discrete events, first described in cardiac muscle (
By this combination of forward (i.e., simulation of sparks) and backward calculations (the release algorithm), we found that the method recovers adequately the current of release flux. We also uncovered limitations of the method, manifested in errors in the determination of release magnitude, source diameter, or release flux duration, and established circumstances that make these errors greater or reduce them.
When the release calculation algorithm was applied to sparks obtained experimentally in skeletal muscle under different conditions, the resulting flux was large by comparison with Ca2+ flux measured through single channels in bilayers under what are thought to be near-physiological conditions (
![]() |
Methods |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Determination of the Point Spread Function
The point spread function (PSF),1 a function P(x,y,z), was determined generally following procedures discussed by (x2 + y2)].
The fluorescence averaged in this way is represented versus r and z in Figure 1 A. The dependence could be fitted as the product of two gaussians: P = G(r) H(z), where G = k exp[-0.5 (r/sigma}">xy)2], H = exp[-0.5(z/
sigma}">z)2],
sigma}">xy = 0.20 µm, and
sigma}">z = 0.614 µm. The two-gaussian description, which is used later in the analysis, is documented in B and C. Figure 1 B contains data as a function of the axial coordinate z. The fluorescence values represented are averages within each section, over a 0.1425-µm circle centered on the x = 0, y = 0 line (), or in a ring region 0.285 < r < 0.4275 µm (
). The lines are best fit gaussians with
sigma}"> = 0.608 and 0.736 µm, respectively, indicating some deviation from the product of gaussians at high values of r. In Figure 1 C, the fluorescence is represented as a function of distance r at constant z. The three sets of data are from the central section (|z| < 0.18 µm; ), an intermediate region (0.18 < |z| < 0.54 µm; dotted symbols), and an outlying region (0.72 < |z| < 0.98 µm;
). The best fit gaussians shown have very similar spread (
sigma}"> = 0.196, 0.2, and 0.207 µm, respectively). In the analysis of data, we neglect the deviations and approximate the PSF as a product of a gaussian function of r, with
sigma}">xy = 0.2 µm (corresponding to a full width at half magnitude [FWHM] of 0.47 µm) and a gaussian function of z, with
sigma}">z = 0.615 µm (FWHM of 1.44 µm). Similar values were obtained with a second C-apochromat 1.2 N.A. objective in our microscope, and with either objective in a LSM 410 (Carl Zeiss, Inc.) confocal scanner.
|
The xy or focal spread measured was close to the expected value for this objective and wavelengths of excitation and emission (e.g.,
Imaging of Fluorescence Sparks
Imaging of fluorescence sparks was carried out in segments of singly dissected semitendinosus muscle fibers of Rana pipiens (which were anesthetized in a 15% ethanol solution, and then killed by pithing), using two conventional techniques for eliciting sparks. In one, fiber segments were voltage clamped in a two-vaseline gap chamber, while their cytoplasmic medium was equilibrated with an "internal" solution containing (mM): 125 Cs-glutamate, 10 Cs-HEPES, 0.5 MgCl2, 1 EGTA, nominal [Ca2+] set to 100 nM, 0.1 fluo-3, 5 creatine-phosphate, 5 Mg-ATP, 5 glucose, pH 7.0. The "external" solution in the middle pool of the chamber contained (mM): 131.5 tetraethylammonium (TEA)methanesulfonate, 10 TEA-HEPES, 10 Ca-methanesulfonate, 10-6 g/liter tetrodotoxin. In this case, sparks were elicited by voltage clamp depolarization to low voltages (between -70 and -50 mV). In the other technique, similar to that described by
Ca2+ sparks were imaged in line scan mode. Images are formed by juxtaposition of 768 line scans of 512 pixels obtained at 2-ms intervals and stored as arrays F(xj,ti), [sometimes referred to as F(x,t)], where j and i are integers varying between 1 and 768 and 1 and 512, respectively. Experimental images are presented normalized to the resting fluorescence F0(x), calculated by averaging F(x,t) over time. In the voltage-clamp experiments, this averaging is restricted to the interval before the depolarizing pulse. In the experiments with permeabilized fibers, it is done after removing the spark regions located automatically by a simple amplitude criterion applied to a filtered version of the image (
Analysis of Fluorescence Images
Once the events were located, all subsequent analysis was carried out on subsets of the line scan, 45 x 45 arrays, centered at the peak of the located spark. Four morphological parameters were determined first: amplitude (a) full width at half magnitude, full duration at half magnitude (FDHM), and rise time. a was determined as the difference between the peak value of the normalized fluorescence and its baseline value at the same spatial position immediately before the event (the average of 10 values between 40 and 20 ms before the peak). FWHM was determined on the spatial distribution of fluorescence at the time of maximum change. FDHM was determined on the time dependence at the spark center. Rise time as the interval between time to 10% increase from pre-event baseline and time to peak, determined on a spline interpolate of the time-dependent fluorescence, averaged over three spatial pixels at the center of the spark. Even though the event locator uses filtration steps, the morphometric measurements were carried out on the unfiltered image, normalized as described. In some cases, averages were made of the events located automatically.
Fluorescence line scan images were analyzed to derive [Ca2+](x,t) and release flux, following a procedure modified from that of
![]() |
(1) |
where K 'd is the dye's dissociation constant in the cytoplasmic medium, R' is the ratio between the fluorescence of Ca2+-bound and free dye, M is the adjustable gain of the scanner, B'min is a constant and dyeT is total dye concentration. (As follows from Equation 1, the minimum and maximum fluorescence were derived as Fmin = B'min M dyeT and Fmax = R' Fmin). The parameters of fluo-3 were given values within the range estimated by
DyeT, a function of position x along the scanned line, was calculated replacing F in Equation 1 by F0(x), the fluorescence averaged over time during the periods of rest before the sparks of interest (or voltage pulse), and assuming the resting concentration [Ca2+]0 to be equal to the concentration in the internal solution, 100 nM.
When [Ca2+]0 is much lower than K'd (for instance at 100 nM), a simplification applies
![]() |
(2) |
where B'max = R' B'min. For the calibration situation Equation 2 applies, with Bmax and Kd substituted for the corresponding values inside the cell. From this calibration, an approximate formula follows for the dye concentration: dyeT/100 µM = (F0/F100) ([Ca2+]0/100 nM) (K'd/Kd) (Bmax/B'max), where F100 is the fluorescence in cuvette at 100 µM [Ca2+] and 100 µM dyeT. Kd is 0.48 µM in our cuvette calibrations. Assuming that the resting [Ca2+] is equal to 100 nM and that Bmax/B'max 1 (
2 x 100 µM (F0/F100). Because F0 is in general a function of x, dyeT is also a function of x. If, instead of 1.03 µM, we used the upper estimate of
When [Ca2+] is changing steeply in space and time, we evaluate it as that needed to produce the observed distribution of dye:Ca2+, the concentration of which is proportional to the observed changes in fluorescence. Specifically, [dye:Ca2+] was derived from the spatially resolved fluorescence F(x,t), using [dye:Ca2+](x,t) = dyeT(x) [F(x,t) - Fmin(x)]/[Fmax(x) - Fmin(x)]. Then [Ca2+](x,t) was obtained numerically solving the diffusionreaction equation that governs the evolution of [dye:Ca2+]:
![]() |
(3) |
where kon and koff are the rate constants of the fluo-3:Ca2+ reaction, DdyeCa is the diffusion coefficient, and is the laplacian operator (
2/
x2 +
2/
y2 +
2/
z2). The line scan only provides the partial derivative in the x direction.
2/
x2 as an approximation to
. M. Cannell (University of Auckland, New Zealand, personal communication) suggested assuming that the fluorescence increase is spherically symmetric, a function of time and the distance (
) to its center. In that case, the dependence of F with x gives all the information needed to calculate the laplacian correctly, as
2/
x2 + 2 (
/
x)/x. This approximation does not work when the sources of release are spatially complex (as in
Calculation of Release Flux
Calculation of release flux was done as described by , satisfies the equation
![]() |
(4) |
where DCa and Ddye:Ca are the diffusion coefficients of Ca2+ and its complex with the dye. rem/
t is a sum of terms of the form
![]() |
(5) |
for each intrinsic buffer b of the cell (parvalbumin, ATP, nondiffusible sites in troponin and the sarcoplasmic reticulum, SR), plus the flux density of sinks (a positive function of space and time). Buffers are assumed to be homogeneously distributed.
The terms of Equation 5 are calculated from the function [Ca2+] (x,t), solving for [b:Ca2+] the diffusionreaction equation of form 3 for the corresponding buffer. In general, this requires simultaneously solving an equation for the free buffer (see for instance Equation 1Equation 2Equation 3Equation 4 in
![]() |
(6) |
assuming [Mg2+] to be constant.
Parameter values are given in Table 1. The SR pump contribution to removal was proportional to the fractional occupancy of pump sites of dissociation constant 1 µM, with maximum (saturation) pump rate given in the table. An exploration of the effects of changing the assumed values of the parameters is presented later.
|
Numerical processing of data was done in the IDL environment (Research Systems Inc.). When specifically stated, the records of fluorescence were first subjected to three-point smoothing (replacing each value by the "boxcar" average of the three-elementwide array that surrounds it). Differentiation with respect to spatial coordinates was carried out by convolution with a 17-element Kaiser window kernel, with a corner frequency of either 2.45 or 1.75 µm-1 (
Simulations (Calculations of Fluorescence Sparks)
Reactiondiffusion equations (), of radial distance to the center of the release sphere), and was generated as a two-dimensional array for suitable values of
and time. To make it comparable with the experimental results, these arrays were spaced at increments of 0.15 µm and 2 ms.
The simulation approach, including the assumption of a homogeneous medium, isotropic for removal and diffusion, was similar to that of
Blurring and Deblurring
Blurring, the effect of the limited resolution of the microscope when imaging fluorescence sparks, was simulated assuming a point spread function similar to that determined experimentally. Let F represent the detected fluorescence, or fluorescence in image space. The object fluorescence, f(x,y,z), is always assumed to be a spherically symmetric function (of radial distance to the center of the spark). In general (
![]() |
(7) |
Because in all cases we processed or simulated line scans, we only concerned ourselves with the dependence of F on X. Additionally, we only processed sparks assumed to be centered on the scanning line. Therefore, the function F above only had to be evaluated at Y = 0 and Z = 0. Representing F(X,0,0) as F(X), Equation 7 simplifies to
![]() |
(8) |
We further simplified Equation 8 because the PSF, as shown in RESULTS, is approximately separable as the product of two gaussian functions: G(x,y) = exp{-[(x2 + y2)/2sigma}">xy2]}/2
sigma}">xy2, a function of the radial distance to the center of the scanned spot in the focal plane, and H(z) = exp{-[(z2)/2
sigma}">z2]}/
sigma}">z
2
, a function of the axial distance to the focal plane, in the z axis direction. The simulated fluorescence spark is spherically symmetrical [f
(
)]. It is well fitted by a gaussian function of polar radius (as shown with an example in DISCUSSION), namely
![]() |
(9) |
With these approximations, the blurred function F(X), Equation 8, simplifies to:
![]() |
(10) |
where G(x - X) G(x - X, 0) and G(y)
G(0, y).
This is done for every point in time and presented to the release calculation algorithm as a line scan image of a spark. A simulated spark before and after blurring is illustrated (see Figure 2A and Figure B).
|
The analysis of fluorescence sparks, real or simulated, starts with a deblurring step to recover the object spark. We assume, in agreement with our analysis of simulated sparks, that the object spark can be put as the spherically symmetric function f(0) g(x)g(y)g(z). The form of the function g is deduced deriving from Equation 10 the expression
![]() |
(11) |
where F(X) is the line scan image, and J and L are constants, the integrals over y and z in Equation 10. Therefore, g(x) is the deconvolution of F and the spatial spread G, scaled so that g(0) = 1. Evaluating Equation 11 at X = 0, it follows that f(0) = F(0)/J 2L. In more intuitive terms, deblurring changes the spatial shape of the spark from F to g, making it sharper, and rescales it by the factor 1/J 2L, which "undoes" the three-dimensional averaging and consequent loss in amplitude, due to the spread of the microscope. J and L quantify these averaging effects and are therefore smaller than 1. The final expression of the deblurred line scan is F(0)g(x)/J 2L, a function that is calculated at every time t j in the experimental arrays, or at arbitrary increments in simulations.
The deblurring operation is tested in Figure 2A and Figure B. The simulated spark before blurring is within 5% of that recovered after successive blurring and deblurring. When applied to noisy data, deblurring may lead to uninterpretably noisy results. In such cases, we used either no or partial deblurring (assuming a sharper PSF). The effects of partial deblurring, evaluated for simulated and experimental data, are illustrated (see Figure 4).
|
|
![]() |
Results |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
The basic self-consistency of the method is established by first applying it to simulated sparks. Then tests are presented on simulated sparks plus noise, which give a more realistic picture of the method's value and limitations. The algorithm is then applied to experimental sparks. The study ends with a consideration of lower bounds of the release estimate, crucial to evaluate the possibility that a spark may be caused by a single channel.
Tests of the Algorithm on Simulated Sparks
Sparks were simulated as described in METHODS, as the increase in fluorescence determined by release of Ca2+ from a spherical source into a homogeneous medium whose removal properties copied those of the skeletal myoplasm. It includes components with the properties of the most important removal processes in the cell; namely, parvalbumin, troponin, the SR pump, the extrinsic buffer EGTA, which in the experiments analyzed later is present at 1 mM, and ATP, a rapidly diffusing low affinity buffer (
Simulated sparks were first blurred (according to Equation 10) with half widths FWHMxy = 0.4 µm and FWHMz = 1.4 µm, to produce the effect of detection by the confocal microscope in its usual configuration. The final result of the simulation, with or without blurring, is represented in Figure 2A and Figure B, as normalized fluorescence profiles that pass through the maximum, at constant t and x. In thin lines is the spark before blurring ("simulated"). In thick trace is the spark after imaging ("blurred"). It is broader and of approximately half the amplitude of the object spark. The release calculation algorithm starts with a deblurring step (Equation 11), which seeks to reconstruct the object spark. Its result is represented by the dashed curves ("deblurred"), and is similar to the original simulated spark.
In Figure 2C and Figure D, the release flux density used in the simulations is compared with that reported by the release algorithm. The release flux density used in the simulations was calculated as I/(2F)/vs (where F is the Faraday constant and vs is the volume of the source, a sphere of radius 0.2 µm) and is represented by the thin trace. The curves in the thick trace represent release flux density calculated by the algorithm ("derived"), in Figure 2 C as a function of position (x) at the maximum of release, and in D as a function of time at the center of the source.
As seen in both representations, the derived release is a reduced and spread out version of the release used in the simulation. The peak of the derived flux is 1,750 mM/s rather than the theoretical 3,100 mM/s. The full duration at half magnitude, 4.7 ms, is in excellent agreement with release duration. The half width of the derived release is somewhat greater (0.46 µm) than that used in the simulation (0.4 µm). Total release flux was computed by volume integration of the spherically symmetric flux density. Expressed as a current (plotted versus time in Figure 2 E) and compared with the current of 20 pA used in the simulation. The current initially lags behind the correct value and later overshoots it, but yields a reasonable estimate. This is shown in Figure 2 F, where the peak of the derived release current is plotted against the simulation current. A good correspondence is found in the range explored, with some evidence of saturation at the high end.
Figure 2 C illustrates that release typically undershoots, going negative at the edges of the releasing volume. This error is due in part to the use of a wide kernel for spatial differentiation. It is also due to the inherent imprecision of an analysis that samples a continuous physical phenomenon at spatial and temporal frequencies that are comparable with those in the sampled object. Though small in amplitude relative to the positive peak of release, the error becomes proportionally greater when total flux is calculated by integration over the volume of the spark. For this reason, when calculating total flux, the volume integration is stopped at the first negative value of the computed release waveform.
The method was also tested with simulations that specifically violated the assumption of spherical symmetry at the source, in which release originated from a cylinder of 1-µm length and 0.2-µm radius (not shown). The current derived by the algorithm in this case overestimated the simulation current when the simulated spark was scanned along the axis of the cylinder and underestimated it when scanning transversely. The errors, however, were within 50% in either direction.
The conclusion from this section is that the release calculation algorithm is self-consistentit recovers approximately the input current used in simulations. More realistic tests required adding noise to the simulated sparks.
Performance of the algorithm in the presence of noise Figure 3 and Figure 4 illustrate performance of the algorithm in the presence of noise. In Figure 3, the algorithm operates on (A) the simulated fluorescence spark of Figure 2 (corresponding to a release current of 20 pA and 5 ms), to which pseudo-random noise has been added, with normal distribution of standard deviation 3 in fluorescence units. Such noise amounts to ~10% of F0, which is comparable with that of averaged experimental records. The main problem resulting from the presence of noise at this level is that the deblurring correction, which increases the higher frequency components more, may give uninterpretably noisy results. In the example, we used a partial deblurring correction for a less severe microscope spread (FWHMz = 0.88 instead of 1.4 µm). The consequences of partial deblurring are explored later.
Figure 3B and Figure C, plots multidimensionally the resulting release flux waveform; quantitative aspects are in DF. Due to the partial correction, the resulting release flux is a worse underestimate than in the case of Figure 2, and the "footprint" of release is wider than the actual source. The half duration of release was 6 ms, and the half width 0.55 µm, or 35% greater than the diameter of the simulation source. The underestimation of flux density and overestimation of spatial width tend to compensate, so that the integrated current peaked at 18.4 pA in the example shown, close to the 20 pA used to generate the spark. In summary, if the noise level prevents application of the deblurring operation, then the release flux density is broadened and more severely underestimated, while the source current is still evaluated fairly well.
Because it was not possible in general to use the full deblurring correction with experimental records, we evaluated in detail the effect of deblurring under different assumptions for the system's PSF. Results are represented in Figure 4 A, where the peak release current derived for simulations analogous to those in Figure 2 and Figure 3 is plotted against "deblurring length," the FWHMz assumed for correction. The values of deblurring length spanned the range from 0 (i.e., no correction), to 1.4 µm ("full correction"). In every case, FWHMxy was equal to FWHz/3.5. (Figure 4) represents simulations without noise and records with noise (at levels comparable with those in Figure 3 and an average spark presented later).
It can be seen that increasing the blurring correction increases the estimate of release current, which tends to the correct values as the deblurring length increases. The changes are moderate and vary continuously with deblurring length. When noise is present, the dependence of the parameters on deblurring length changes abruptly at or near 1 µm, indicating that the deblurring correction does not work properly beyond this point (another evidence, not shown, is that the results become variable with each different realization of pseudo-random noise). This test indicates that the algorithm applied with partial deblurring still evaluates source current acceptably, to some extent underestimating it. This issue will be considered further when processing experimental sparks.
Additional tests evaluated the ability of the algorithm to monitor current when spatial aspects of the source were varied. Simulations were similar to those shown before (20 pA, 5 ms, noise at 10% of baseline fluorescence, deblurring at 0.88 µm), but the source diameter varied between 0.2 and 1.5 µm. in Figure 4B and Figure C, plots peak current and FWHM of the derived release against source diameter. It can be seen that the method evaluates current well at every diameter, but only follows the changes in diameter at values beyond 0.5 µm (C). The spatial dimension of the source is overestimated at diameters below 0.6 µm, an evidence of limited spatial resolution, due to noise and the optical spread of the microscope.
, in Figure 4B and Figure C, illustrates results of applying the algorithm to simulations with noise at 30% of F0, a level characteristic of individual experimental sparks. Under these conditions, the algorithm can be used provided that fluorescence images are smoothed first, that a more heavily filtering differentiator is used (see METHODS), and that no deblurring correction is applied. Given the high level of noise, the results change with different noise realizations. Figure 4,
, represents averages of six realizations, bars are ±SEM. In this form, which can be applied to individual experimental sparks, the algorithm still evaluates release current within 50% error (B), but fails to report changes in source diameter under 1 µm (C). In summary, the algorithm evaluates release current well under realistic conditions of noise, within a wide range of source parameters. It is able to detect changes in source diameter above 0.7 or 0.8 µm when used with averages, or above 1 µm when used with individual sparks.
The Release Current of Experimental Sparks
The characteristics of sparks obtained with two different experimental techniques were somewhat different.
Cells under voltage clamp.
Figure 5 A contains a line scan image obtained from a fiber under voltage clamp, which had been exposed for 45 min to an internal solution with a nominal free [Ca2+] of 100 nM and 610 µM [Mg2+]. In all experiments prepared this way, almost no sparks were observed at rest, but small pulse depolarizations elicited many, as illustrated. The plot at bottom represents the relative fluorescence averaged in a 0.43-µm wide region at the position marked by the line aa' in Figure 5 A. The large spark there peaked at 6.0x the resting fluorescence (F/F0 = 5.0). Under voltage clamp, repeated depolarizations can be applied and large numbers of sparks collected in a short period of time within a small region of the fiber. We illustrate below the application of the release algorithm to an average of such sparks.
|
The average was made with the events identified and measured by an automatic routine (Cheng et al., 1999) in Figure 5 A and in nine other images obtained during a brief period. A total of 301 sparks were identified whose amplitude was >1.0 U resting fluorescence. Their morphology is illustrated in Figure 5B and Figure C, which plot, respectively, FWHM and rise time versus spark amplitude. The events of amplitude >2.0 U are marked with circles in A and identified with in the scatter plots.
The graphs demonstrate several aspects of spark morphology that were similar in five experiments studied in the same detail (listed in Table 2). Events span a wide range of amplitudes, between seven and the limit of detection. There does not seem to be any discrete segregation in morphology among sparks; that is, no systematic difference appears between the morphological parameters of large versus small amplitude sparks. Accordingly, the correlation between amplitude and the other variables is weak or nonexistent. The lines represent first order regressions to the open symbols. The regression coefficient r2 was 0.102 for the spark width (range 0.0120.127 in six experiments) and 0.008 (0.0000.017) for the rise time. This low correlation suggests that the observed spread in amplitude was not a consequence of the existence of multiple spatially resolved sources, which could collectively generate a large spark (such as those described in cardiac muscle by
|
The amplitude with which a spark is recorded in a line scan depends on the actual (object) amplitude of the spark and on its position relative to the scanned line, reaching a maximum (image) amplitude when it is centered on the scanned line. Events that appear large in a line scan image are likely to be closer to the scanned line, hence more suitable to our release analysis (which assumes centered sparks). Of course, sparks that are larger as objects are also likely to appear larger in line scans. In the following, we present the results of the release analysis applied to the average of the 157 events with amplitude >2.0 identified in the experiment illustrated. For the considerations above, this average selects a subset of events that are larger as objects and/or originated closer to the scan line.
The application of the release algorithm to the average spark is illustrated in Figure 6. Parameter values are listed in Table 1; [Mg2+] and [Ca2+] were assumed to have the nominal internal solution values. The dye concentration was set at the average of the estimates for each event (which ranged between 38 and 42 µM). The average was centered at the peaks of the sparks and is shown in Figure 6 A. Because of the low level of noise in the average, no smoothing or other filtering was applied. Partial deblurring was used with deblurring length = 0.88 µm; results with other degrees of deblurring are represented in Figure 4 A, . The release flux density calculated by the algorithm is represented in Figure 6BE. The derived release current is plotted in F.
|
The algorithm applied to this average spark results in flux that arises from a region 0.5 µm wide, and has a half magnitude duration of 6.3 ms and an amplitude of 16.9 pA. It is demonstrated later that duration and width are quite insensitive to assumptions in the method. The calculated current, however, depends steeply on [Ca2+]0. When, instead of 100 nM (the concentration set in the cut-end solutions), [Ca2+]0 was assumed to be 50 nM, peak current was reduced by ~35% (Figure 6 F). The dependence is explicable because [fluo-3] is derived from resting fluorescence, which is approximately proportional to both [fluo-3] and [Ca2+]0, hence a change in assumed [Ca2+]0 alters [fluo-3] and the derivation of [Ca2+](x,t).
In every fiber prepared the same way, several hundred sparks were identified (100 or more of which were >2.0 in amplitude). The average morphological parameters, listed in Table 2, and the calculated dyeT, were very similar in all cases, hence the release must have been similar.
Permeabilized fibers.
It was possible to reduce the uncertainty in the derivations by using the permeabilized fiber preparation, described in METHODS, in which the fiber membrane is cut at multiple locations and the cell is immersed in the internal solution. Intracellular fluorescence increases with a time constant of 1020 min. For analysis of these experiments, it is assumed that free [Ca2+] has the value present in the internal solution (an assumption consistent with model simulations of
|
Unlike the situation with voltage-clamped fibers, the sparks obtained with permeabilized fibers were fewer and did not allow for massive averaging. Limited averaging could be done occasionally, as in the case of Figure 7. The image shows several sparks arising at the same triad in a field that was of very low activity elsewhere. Therefore, it is likely that all those sparks were generated at the same release unit (
Figure 8 illustrates the application of the release algorithm to an average of the four largest sparks in Figure 7. No deblurring was used. The parameter values were the same as in the previous case, except for [Mg2+], which was set to the concentration applied, 0.15 mM. Figure 8 A represents the spark average, after smoothing, and BE represent the calculated release flux density. The current (F) peaked at 15.8. Half width of the release source was 1.1 µm and half duration was 7.3 ms.
|
Release was calculated in the same way for large sparks in 16 images from 12 similarly prepared fibers. The results are listed in Table 3. The average parameters of the release source include a current of 14.4 pA, duration of 7.5 ms, and half width of 0.89 µm.
|
Uncertainties in the Calculation of Release Flux
As illustrated above, the algorithm produces adequate estimates of release current when applied to sparks that are generated by the same process modeled in the algorithm and using the same parameter values. Neither condition applies with real data; there may be differences in parameter values, or other differences (for instance lumped components, diffusion barriers, etc., not considered in the present homogeneously distributed model of removal).
As a way to gauge the likely range of the first class of errors, we imposed parameter changes in the calculation of release for the average of 157 sparks studied in Figure 6. The changes and the resulting morphological measures of derived release current are listed in Table 4. Changes in parameters were introduced one at a time. For those parameters expected to have a major effect, two alternative values were used, defining a reasonable range that should comprise the actual fiber values. This was done for the diffusion coefficients of Ca2+, the dye, and ATP, and for the concentrations of ATP, EGTA, and Mg2+. For maximum pump rate and the concentrations of troponin and parvalbumin sites, only one change was explored, sufficient to demonstrate a relative lack of consequences.
|
Changing properties of the dye:Ca2+ reaction had the greatest effects. When the ON rate constant was reduced by a factor of 2.5, to bring K'd to the upper value in the range estimated in vivo by
To elicit a "low estimate" of release current, the method was applied with all buffer concentrations reduced threefold. The derived release parameters are in the last row of Table 4. With such low values, the release current was 8.05 pA.
A minimum estimate of another sort was reached by simply integrating the dye-binding term ([dye:Ca2+]/
t) in the representation of flux density given by Equation 4. The meaning of this estimate is illustrated in Figure 9, where it is developed for both the average of four sparks in Figure 8 and a simulated spark copying the experimental conditions. Figure 9, top, shows the concentration of dye:Ca2+, which for the experimental spark is calculated from the fluorescence following Equation 1. Figure 9, middle, shows the rate of change of [dye:Ca2+], which in both cases reaches a pronounced maximum right before the peak of fluorescence. The bottom graphs plot this rate of change (
[dye:Ca2+]/
t) at the time of its maximum, versus x. Intuitively, the Ca2+ release flux at this time should be greater than the total flux of Ca2+ binding to the dye, which can be calculated by volume integration of the rate of change shown. The result of the volume integration, expressed as a current intensity, is plotted in the bottom graphs as a function of position, starting from zero at the center of the spark. The limiting value, reached far from the spark center, is ~15 pA (0.074 fmol/s) in the experimental case. The simulated spark reaches a similar maximum rate of change, but has lesser spatial spread. Accordingly, the total flux binding to the dye in the simulation corresponds to <4 pA. In every case tested, the integrated rate at early times during a spark was less than the maximum release flux, in agreement with the idea that the other removal and diffusion processes are largely positive at these early times. The value of this dye-related estimate is listed for all events analyzed in permeabilized fibers, in the last column of Table 3. Its average is 8.26 pA.
|
Different microscope objectives. To test the possibility that release evaluations were affected by the less than ideal axial resolution determined in our conditions, we compared the spark morphology and release flux with our preferred C-apochromat objective and the 63x, 1.4 N.A., oil immersion plan-apochromat (the PSFs of both objectives are compared in METHODS). Pooled averages and their standard errors, for 1,850 sparks obtained with the C-apochromat and 1,591 in five experiments in which the system was aligned and all images were acquired with the plan-apochromat were, respectively: amplitude, 1.09 (0.04) and 1.10 (0.05); FWHM (µm), 1.40 (0.07) and 1.57 (0.09); FDHM (ms), 9.10 (0.4) and 9.87 (0.4). None of these differences were significant. (The average amplitudes were smaller than in Table 2 due to the use of a lower detection criterion.) In one experiment, the system was aligned and initial images were acquired with the C-apochromat, and then the fiber was imaged with the plan-apochromat and additional images were acquired. Again, no significant differences were observed between averages of 239 and 198 events obtained with the two objectives. Because sparks were similar, the release algorithm applied to the largest sparks or to spark averages gave similar results in both types of experiments with either objective. The conclusion is that the better axial resolution provided by the oil-immersion objective does not result in a spark of greater amplitude, which may be an indication of differences in the PSF that require more than the FWHM for their description. We continue to prefer the water immersion C-apochromat because the morphology of events imaged with it does not appear to depend on depth within the cell, which allows one to obtain many more images from the same fiber.
![]() |
Discussion |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
An algorithm introduced for the evaluation of release flux in cardiac Ca2+ sparks (, a function of space and time), which can then be integrated to obtain total flux (or current).
Evaluation of the Algorithm
We first tested the consistency of the algorithm (or "backward" calculation) on sparks simulated as the result of Ca2+ diffusion from a finite source (the "forward" calculation). The backward calculation retrieved approximately the same total current used in the forward construction, over a wide range of release current. There were some deviations, namely the flux density in the backward calculation was ~5560% of that in the simulation, and the half width of the region where the back-calculated source resided was slightly greater than the simulation source. These errors appear to be due to the finite spatial and temporal resolution of the experimental methods, which requires discrete and necessarily coarse approximations to steeply space- and time-dependent variables. They approximately compensated each other, however, resulting in a nearly correct determination of total current.
A set of more realistic tests, which characterized the performance and limitations of the method when applied to experimental data, included noise in the records. Because Ca2+ removal and the monitoring dye reaction are mass action-driven processes that occur at finite rates, the release algorithm involves differentiation stages that increase noise. Additionally, it starts with a deblurring step, to undo the blurring introduced by imaging. Deblurring is basically a deconvolution resulting in a greater enhancement of the high (spatial) frequency components of the signal, and a disproportionate increase in noise when applied to real data.
To counter these problems, we altered the method in two ways, we used less than full deblurring (that is, deblurring for a sharper PSF than actually measured) and we used smoothing and filtering digital differentiation. The effects of changing the deblurring length (FWHM of the PSF assumed in the correction) was explored on simulated sparks. This study, illustrated in Figure 4, shows that increasing the deblurring length gradually increases the estimate of release current. The change is gradual, albeit sizable (in going from no deblurring to full deblurring the current estimate may increase by as much as 40%). In the same tests, noise, added to the simulations at the levels found in spark averages, results in an erratic morphology of derived release beyond a certain point. This translates to a limiting value of ~0.9 µm for the deblurring length.
Using this practical level of deblurring correction, we then asked how would the algorithm depend on the spatial aspects of the source, namely the diameter (because the source is assumed to be spheric). As shown in Figure 4, the algorithm underestimated source current by ~30%, regardless of source diameter. The source radius reported by the algorithm was rather insensitive to changes in actual source geometry below 0.7 µm. This limits the ability of the technique to distinguish between extended and point sources.
In summary, the algorithm slightly underestimates source current, slightly overestimates source duration, and reflects source diameter beyond ~0.7 µm. Its performance is degraded when used with digital low-pass filtering, but release current is never overestimated.
Applications to Experimental Sparks
The method was applied to an average including every spark of amplitude above 2 F0 obtained under voltage clamp in the same cell. Because the images were acquired over a short period, the dye concentration varied <10%, which makes it legitimate to average the events and then process the average. There was a poor correlation between amplitude and rise time for these sparks, as well as a low positive correlation between amplitude and half width. This suggests that the selected sparks were not "special", other than by surpassing the amplitude criterion. Their large amplitude indicates that the sparks originated near the scanning line, were large as objects (i.e., before imaging), or both.
The derived source current is represented in Figure 4 () for various values of deblurring length. Comparing with the values obtained with simulations, it appears that the average spark originated from a source of at least 20 pA (the release current in the simulation). The range of source widths derived for the average spark, 0.470.7 µm, does not allow one to tell whether the flux originated at a point or an extended source.
The algorithm was also applied to sparks from permeabilized fibers immersed in internal solution. This technique had two advantages. Dye concentrations in the few hundred micromolar were achieved rapidly inside the cells, leading to a more robust determination of release flux. An even more important advantage of the permeabilized preparation was the better determination of resting [Ca2+], which should have been close to that in the equilibrating internal solution. The method was applied to individual events or to averages of just a few, which typically required low-pass filtering of the images and did not allow for deblurring.
For 19 large sparks from 11 permeabilized cells, the average release current was 14.4 pA, and half width 0.89 µm. Considering the conditions under which the algorithm was applied, the release current is probably an underestimate, and therefore not inconsistent with the estimate of 20 pA obtained under voltage clamp.
The large source width derived under these conditions should not be overinterpreted because the algorithm applied to simulations with high noise levels was quite insensitive to source width in the range below 1 µm. Nevertheless, it is interesting that in several experiments source diameters >1 µm were derived. In all those cases, fluorescence sparks had large width. One way in which sparks of such width could be simulated in forward calculations was by assuming source diameters >1 µm. These observations therefore raise the possibility that some wide sparks of high amplitude may originate from sources of resolvable size. One problem with this interpretation is that scanning was perpendicular to the Z disks (and the long axis of the transverse tubules), hence it is difficult to imagine a structural substrate for such extended release sources. The large width of some sparks therefore remains without a convincing explanation or simulation.
Estimates of Release Current and Their Implications
The present estimates of release current underlying a spark should be compared with single channel currents recorded in bilayers, in conditions as close as possible to physiological. Mejía-Alvarez et al. (1988) studied heavy SR cardiac release channels carrying Ca2+ current in the lumenal-to-cytoplasmic direction, driven by gradients of 210 mM, and in the presence of high concentrations of Cs+ (as a substitute for the inhibitory effects of K+ on Ca2+ current). At 2 mM lumenal Ca2+ (similar to [Ca2+]SR values determined in cardiac muscle by
The comparison implies that ~60 fully open channels are necessary to account for release in the sparks studied here. This number, however, is subject to multiple possible errors in parameter values or in the structure of the model itself. For this reason, we provide two lower estimates of current. One is 8.09 pA, obtained when concentrations of all the removal molecules of the muscle cell model were reduced to one third their consensus values. Another low estimate, not requiring any assumptions about removal processes, was obtained by computing the rate of change of dye-bound Ca2+ before the fluorescence peak. This estimate is equivalent to the rate of change of "signal mass," defined by
The unitary current could be greater if [Ca2+]SR was higher than the 2 mM used in the bilayer experiments. The range of possible values of [Ca2+]SR is limited, however. In cardiac muscle, [Ca2+]SR has been found to vary within a narrow band (
Three recent reviews (
In a model of
The present conclusions apply to sparks that were selected for their large amplitude. Examination of the amplitude distribution of sparks in view of the artifacts introduced by the process of line scanning (
Ca2+ Release in Sparks
![]() |
Footnotes |
---|
Dr. Shirokova's permanent address is A.A. Bogomoletz Institute of Physiology, Bogomoletz St. 4, Kiev, Ukraine.
![]() |
Acknowledgements |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
We are grateful to Dr. Wolfgang G. Kirsch (Rush University) for suggestions on the manuscript.
This work was supported by grants from the National Institutes of Health (NIH) and the Muscular Dystrophy Association to E. Ríos, by intramural research programs of NIH to M.D. Stern, and a grant from NIH to N. Shirokova. A. González was the recipient of a Senior Fellowship from the American Heart Association of metropolitan Chicago.
Submitted: October 13, 1998; Revised: May 7, 1999; Accepted: May 10, 1999.
1used in this paper: FDHM, full duration at half magnitude; FWHM, full width at half magnitude; PSF, point spread function; SR, sarcoplasmic reticulum
![]() |
References |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Agard, D.A., Hiraoka, Y., Shaw, P., Sedat, J.W. (1989) Fluorescence microscopy in three dimensions. Methods Cell Biol. 30:353-377[Medline].
Baylor, S.M., Chandler, W.K., Marshall, M.W. (1983) Sarcoplasmic reticulum calcium release in frog skeletal muscle fibres estimated from Arsenazo III calcium transients. J. Physiol. 344:625-666[Abstract].
Baylor, S.M., Hollingworth, S. (1998) The transient binding of calcium to ATP and diffusion of Ca-ATP help shape the amplitude and time course of the myoplasmic free Ca transient. Biophys. J. 74:A235. (Abstr.).
Bhat, M.B., Zhao, J., Zang, W., Balke, C.W., Takeshima, H., Wier, W.G., Ma, J. (1987) Caffeine-induced release of intracellular Ca2+ from Chinese hamster ovary cells expressing skeletal muscle ryanodine receptor. Effects on full-length and carboxyl-terminal portion of Ca2+ release channels. J. Gen. Physiol. 110:749-762
Blatter, L.A., Hüser, J., Ríos, E. (1997) Sarcoplasmic reticulum Ca2+ release flux underlying Ca2+ sparks in cardiac muscle. Proc. Natl. Acad. Sci. USA. 94:4176-4181
Bootman, M., Niggli, E., Berridge, M., Lipp, P. (1997) Imaging the hierarchical Ca2+ signalling system in cells. J. Physiol. 499:307-314[Abstract].
Cannell, M.B., Soeller, C. (1999) Perspective: mechanisms underlying calcium sparks in cardiac muscle. J. Gen. Physiol. 113:373-376
Chen, W., Steenbergen, C., Levy, L.A., Vance, J., London, R.E., Murphy, E. (1996) Measurements of free Ca2+ on sarcoplasmic reticulum in perfused rabbit heart loaded with 1,2-bis(2-amino-5,6-difluorophenoxy)ethane-N,N,N',N'-tetraacetic acid by 19F NMR. J. Biol. Chem. 271:7398-7403
Chen, W., London, R.E., Murphy, E., Steenbergen, C. (1998) Regulation of the Ca2+ gradient across the sarcoplasmic reticulum in perfused rabbit heart. A 19F nuclear magnetic resonance study. Circ. Res. 83:898-907
Cheng, H., Lederer, W.J., Cannell, M.B. (1993) Calcium sparks: elementary events underlying excitationcontraction coupling in heart muscle. Science. 262:740-744[Medline].
Cheng, H., Song, L.S., Shirokova, N., González, A., Lakatta, E.G., Ríos, E., Stern, M.D. (1998) Amplitude distribution of calcium sparks in confocal images. Theory and studies with an automatic detection method. Biophys. J. 76:606-617
Franzini-Armstrong, C., Jorgensen, A.O. (1994) Structure and development of E-C coupling units in skeletal muscle. Annu. Rev. Physiol. 56:509-534[Medline].
González, A., Shiorkova, N., Kirsch, W.G., Pizarro, G., Stern, M.D., Cheng, H., Ríos, E. (1999) Size distribution of Ca2+ sparks in caffine-stimulated skeletal muscle. Biophys. J. 76:A385.
Hamming, R.W. 1989. Digital Filters. Third edition. Prentice Hall, Englewood Cliffs, CA. 189 pp.
Harkins, A.B., Kurebayashi, N., Baylor, S.M. (1993) Resting myoplasmic free calcium in frog skeletal muscle fibers estimated with fluo-3. Biophys. J. 65:865-881[Abstract].
Inesi, G. (1994) Teaching active transport at the turn of the 21st century. Biophys. J. 66:554-560[Abstract].
Klein, M.G., Cheng, H., Santana, L.F., Jiang, Y.H., Lederer, W.J., Schneider, M.F. (1996) Two mechanisms of quantized calcium release in skeletal muscle. Nature. 379:455-458[Medline].
Klein, M.G., Lacampagne, A., Schneider, M.F. (1999) A repetitive mode of activation of discrete Ca2+ release events (Ca2+ sparks) in frog skeletal muscle fibers. J. Physiol. 515:391-411
Kushmerick, M.J., Podolsky, R.J. (1969) Ionic mobility in muscle cells. Science. 166:1297-1298[Medline].
Lacampagne, A., Klein, M.G., Schneider, M.F. (1998) Modulation of the frequency of spontaneous sarcoplasmic reticulum Ca2+ release events (Ca2+ sparks) by myoplasmic [Mg2+] in frog skeletal muscle. J. Gen. Physiol. 111:207-224
Lacampagne, A., Ward, C.W., Klein, M.G., Schneider, M.F. (1999) Time course of individual voltage-activated Ca2+ sparks recorded at ultra-high time resolution in frog skeletal muscle. J. Gen. Physiol. 113:187-198
Lipp, P., Niggli, E. (1996) A hierarchical concept of cellular and subcellular Ca2+ signalling. Prog. Biophys. Mol. Biol. 65:265-296[Medline].
Lipp, P., Niggli, E. (1998) Fundamental calcium release events revealed by two-photon excitation photolysis of caged calcium in guinea-pig cardiac myocytes. J. Physiol. 508:801-809
Mejía-Alvarez, R., Kettlun, C., Ríos, E., Stern, M.D., Fill, M. (1998) Unitary calcium currents trough cardiac ryanodine receptors under physiological conditions. J. Gen. Physiol. 113:177-186
Melzer, W., Ríos, E., Schneider, M.F. (1984) Time course of calcium release and removal in skeletal muscle fibers. Biophys. J. 45:637-641[Abstract].
Melzer, W., Ríos, E., Schneider, M.F. (1987) A general procedure for determining the rate of calcium release from the sarcoplasmic reticulum in skeletal muscle fibers. Biophys. J. 51:849-863[Abstract].
Parker, I., Zang, W.J., Wier, W.G. (1996) Ca2+ sparks involving multiple Ca2+ release sites along Z-lines in rat heart cells. J. Physiol. 497:31-38[Abstract].
Pechère, J.F., Capony, J.P., Demaille, J. (1973) Evolutionary aspects of the structure of muscular parvalbumins. Syst. Zool. 22:533-548.
Pizarro, G., Csernoch, L., Uribe, I., Rodríguez, M., Ríos, E. (1991) The relationship between Q and Ca release from the sarcoplasmic reticulum in skeletal muscle. J. Gen. Physiol. 97:913-947[Abstract].
Pratusevich, V.R., Balke, C.W. (1996) Factors shaping the confocal image of the calcium spark in cardiac muscle cells. Biophys. J. 71:2942-2957[Abstract].
Ríos, E., Stern, M.D. (1997) Calcium in close quarters: microdomain feedback in excitationcontraction coupling and other cell biological phenomena. Annu. Rev. Biophys. Biomol. Struct. 26:47-82[Medline].
Sako, Y., Sekihata, A., Yanagisawa, Y., Yamamoto, M., Shimada, K., Ozaki, K., Kusumi, A. (1997) Comparison of two-photon excitation laser scanning microscopy with UV-confocal laser scanning microscopy in three-dimensional calcium imaging using the fluorescence indicator Indo-1. J. Microsc. 185:9-20[Medline].
Schneider, M.F. (1999) Perspective: Ca2+ sparks in frog skeletal muscle: generation by one, some, or many SR Ca2+ release channels? J. Gen. Physiol. 113:365-371
Schneider, M.F., Ríos, E., Melzer, W. (1987) Determining the rate of calcium release from the sarcoplasmic reticulum in muscle fibers. Biophys. J. 51:1005-1007[Medline].
Shirokova, N., García, J., Ríos, E. (1998) Local calcium release in mammalian skeletal muscle. J. Physiol. 512:377-384
Shirokova, N., Ríos, E. (1996) Activation of Ca2+ release by caffeine and voltage in frog skeletal muscle. J. Physiol. 493:317-339[Abstract].
Shirokova, N., Ríos, E. (1997) Small event Ca2+ release: a probable precursor of Ca2+ sparks in frog skeletal muscle. J. Physiol. 502:3-11[Abstract].
Shtifman, A., Ward, C.W., Valdivia, H.H., Schneider, M.F. (1999) Induction of long duration Ca2+ release events by Imperatoxin A in frog skeletal muscle. Biophys. J. 76:A465. (Abstr.).
Smith, G.D., Keizer, J.E., Stern, J.D., Lederer, W.J., Cheng, H. (1998) A simple numerical model of calcium spark formation and detection in cardiac myocytes. Biophys. J. 75:15-32
Smith, P.D., Liesegang, G.W., Berger, R.L., Czerlinski, G., Podolsky, R.J. (1984) A stopped-flow investigation of calcium ion binding by ethylene glycol bis(ß-aminoethylether)-N,N'-tetraacetic acid. Anal Biochem. 143:188-195[Medline].
Stern, M.D., Pizarro, G., Ríos, E. (1997) Local control model of excitationcontraction coupling in skeletal muscle. J. Gen. Physiol. 110:415-440
Sun, X.P., Callamaras, N., Marchant, J.S., Parker, I. (1998) A continuum of InsP3-mediated elementary Ca2+ signalling events in Xenopus oocytes. J. Physiol. 509:67-80
Tsugorka, A., Ríos, E., Blatter, L.A. (1995) Imaging elementary events of calcium release in skeletal muscle cells. Science. 269:1723-1726[Medline].
Uttenweiler, D., Weber, C., Fink, R.H. (1998) Mathematical modeling and fluorescence imaging to study the Ca2+ turnover in skinned muscle fibers. Biophys. J. 74:1640-1653
Wilson, T. (1995) The role of the pinhole in confocal imaging system. In Pawley J.B., ed. Handbook of Biological Confocal Microscopy. New York, Plenum Publishing Corp, 167-182.