Intracellular calcium movements during excitation-contraction coupling in mammalian slow-twitch and fast-twitch muscle fibers.

In skeletal muscle fibers, action potentials elicit contractions by releasing calcium ions (Ca(2+)) from the sarcoplasmic reticulum. Experiments on individual mouse muscle fibers micro-injected with a rapidly responding fluorescent Ca(2+) indicator dye reveal that the amount of Ca(2+) released is three- to fourfold larger in fast-twitch fibers than in slow-twitch fibers, and the proportion of the released Ca(2+) that binds to troponin to activate contraction is substantially smaller.

INTRODUCTION

The skeletal muscles of mammals account for ∼40% of body weight and for ∼25% of the basal metabolic rate (Rolfe and Brown, 1997). Approximately half of this muscle mass consists of slow-twitch fibers, and the other half consists of fast-twitch fibers (also called type I and type II fibers, respectively). Small numbers of other fiber types, such as superfast fibers and true-slow fibers, are also found in specialized muscles such as eye and jaw muscles (Hess, 1970; Hoh, 2002; Andrade et al., 2005). The basic unit of activity in the two main fiber types is the “twitch,” in which a single action potential (AP) on the fiber’s surface membrane elicits a brief contractile response. This article concerns the transient increase in the myoplasmic free calcium concentration (Δ[Ca2+]), which is the key intermediate signaling event between excitation and contraction. As [Ca2+] rises, Ca2+ binds to troponin, a regulatory protein on the actin filament, thereby enabling myosin to interact with actin and generate the contractile response. When [Ca2+] falls, Ca2+ dissociates from troponin and contractile activity declines.

Slow-twitch fibers have many similarities and also some important differences when compared with fast-twitch fibers (Berchtold et al., 2000; Allen et al., 2008; Schiaffino and Reggiani, 2011). Slow-twitch fibers are more resistant to fatigue than fast-twitch fibers, partly because of their larger mitochondrial volume and higher activity of oxidative enzymes (Schwerzmann et al., 1989; Jackman and Willis, 1996). Also, as the naming system implies, the time course of the contractile response of slow-twitch fibers is substantially slower than that of fast-twitch fibers. For example, in rat motor units studied in the anesthetized animal, the time of peak of the twitch response and the time between peak and half-relaxation are three- to fivefold larger in slow-twitch units of the soleus muscle than in fast-twitch units of the extensor digitorum longus (EDL) muscle (Table 1, A). Similar differences in twitch kinetics are found in individual fibers and in small bundles of fibers from mouse soleus and EDL muscles studied after removal of the muscles from the animal (Table 1, B and C).

Table 1.

Twitch kinetics of motor units and muscle fibers from young adult rats and mice

Experimental preparation	n	Time of peak	Half-relaxation time	FDHM
		ms	ms	ms
A. Motor units in anesthetized ratsa (normal sarcomere length; 35-36°C)
1. Soleus muscle (slow-twitch units)	34	38 ± 4	57 ± 8	—
2. EDL muscle (fast-twitch units)	48	11 ± 1	11 ± 2	—
B. Singly dissected mouse fibersb (normal sarcomere length; 20-21°C)
1. Soleus fibers (slow-twitch)	12	96 ± 7.5	130 ± 14	—
2. EDL fibers (fast-twitch)	11	33 ± 6.1	62 ± 12	—
C. Small bundles of mouse fibersc(long sarcomere length; 16°C)
1. Soleus fibers (slow-twitch)	14	173 ± 9	708 ± 73	843 ± 80
2. EDL fibers (fast-twitch)	16	52 ± 2	162 ± 11	197 ± 12

n is the number of motor units (A) or experiments (B and C). Time of peak is the time between the external stimulus and the peak of the twitch tension response, half-relaxation time is the time between the peak and half-decay of tension, and FDHM is the time between the half-rise and half-decay of tension. In A and B, the sarcomere length was set to give the maximum tension response. In C, average sarcomere length was 3.7 µm (range, 3.3–3.9 µm).

Close, 1967.

González et al., 2000.

Hollingworth et al., 1996, 2008, 2012; and Baylor and Hollingworth, 2003.

Part of the difference in contractile speed between slow-twitch and fast-twitch fibers is attributable to a difference in myosin isoforms. For example, in rat muscle, the major myosin isoform found in soleus fibers has a cycling rate at saturating [Ca2+] that is about one third of that of the major isoform in EDL fibers (Unsworth et al., 1982; Bottinelli et al., 1994). Some of the fiber-type differences in twitch kinetics likely also arise because of differences in Δ[Ca2+] (see below) and in the troponin isoforms (van Eerd and Takahashi, 1976; Potter et al., 1977; Sweeney et al., 1990).

The Ca2+ cycle during a muscle twitch

Fig. 1 diagrams the major intracellular movements of Ca2+ that take place during a twitch. An AP initiated on the surface membrane propagates into the transverse (T)-tubular system, where specialized voltage-sensing molecules (dihydropyridine receptors [DHPRs]) move to their activating position (Schneider and Chandler, 1973; Hollingworth and Marshall, 1981; Ríos and Brum, 1987; Tanabe et al., 1988). The DHPRs then activate the SR Ca2+ release channels, whose large cytoplasmic domain spans the gap between the T-tubular and SR membranes at the triadic junction (Franzini-Armstrong and Nunzi, 1983). As the release channels open, Ca2+ moves from the SR into the myoplasm under a large chemical driving force. A large local increase in [Ca2+] ensues, which drives the diffusion of Ca2+ throughout the sarcomere (Cannell and Allen, 1984; Escobar et al., 1994; Baylor and Hollingworth, 1998; Hollingworth et al., 2000; Gómez et al., 2006). As Ca2+ diffuses, it binds to sites on various myoplasmic constituents (termed Ca2+ “buffers” in this article), the most important of which are indicated in Fig. 1: (a) troponin, which, as noted above, has regulatory sites (Ca2+-specific sites) that mediate contractile activation; (b) ATP, which acts as a mobile buffer for Ca2+ and Mg2+; (c) more slowly reacting, less mobile buffers for Ca2+ and Mg2+ (“SlowB”), such as parvalbumin; and (d) the SR Ca2+ pump, which resequesters Ca2+ within the SR, translocating two Ca2+ ions for each ATP consumed (Weber et al., 1966). (e) If a Ca2+ indicator dye is introduced into the myoplasm, some Ca2+ will also bind to the dye. If the indicator is well chosen and other features of the experiment are favorably arranged, the optical signal from the indicator can be used to monitor the amplitude and time course of Δ[Ca2+] (described in next section). From Δ[Ca2+] and information from the literature about the concentrations and kinetic reactions of the major myoplasmic Ca2+ buffers, quantitative estimates can be made about the complexation of Ca2+ with the buffers and thus about the amount and time course of SR Ca2+ release and reuptake (compare calculations described below).

Figure 1.

Schematic of intra-sarcomeric Ca2+ movements during excitation–contraction coupling in mammalian skeletal muscle fibers. Depolarization of the transverse tubular membranes initiates Ca2+ release from the SR at the triadic junctions. The rise in myoplasmic free [Ca2+] drives Ca2+ complexation with the major myoplasmic Ca2+ buffers: (a) troponin (located on the thin filaments); (b) ATP; (c) other proteins (“SlowB”) such as parvalbumin, which, in rodent fibers, is found at a large concentration in fast-twitch glycolytic fibers but not in slow-twitch fibers (Table 2); (d) Ca2+ pump molecules (located in the SR membranes); and (e) the Ca2+ indicator dye (which can be introduced experimentally to monitor Δ[Ca2+]).

Table 2.

Concentrations of divalent ion–binding sites in the computational models

Type of site	Slow-twitch fibers	Fast-twitch fibers
	µM	µM
Ca2+-specific sites:
Troponin regulatory sitesa	120	240
Ca2+/Mg2+ sites:
ATPb (one per molecule)	5,000	8,000
SlowB (slow-buffer sites)
Parvalbuminc (two/molecule)	0	1,500
Troponin-nonspecific sitesd (two/molecule)	240	0
Ca2+ pumpe (two per molecule)	96	240
Furaptraf (one per molecule)	100	100

The site concentrations are spatially averaged and referred to the myoplasmic water volume (Baylor et al., 1983). The fast-twitch model ignores the 240-µM troponin-nonspecific sites (Ca2+/Mg2+ sites) because their buffering effect is small compared with the sites on parvalbumin, which are present in the fast-twitch model at a large concentration. Other myoplasmic Ca2+-binding sites, such as those on inorganic phosphate, calmodulin, calpain, calcineurin, sorcin, annexin, S100A, the plasmalemma Ca2+ pump, and the Na+/Ca2+ exchanger, have been ignored because their concentrations are small relative to those of the major buffers listed in the table. The kinetic parameters for the sites at 16°C are given in Baylor and Hollingworth (2007) and Hollingworth et al. (2012). Both models assume that free [Ca2+] at rest is 0.05 µM and that free [Mg2+] is 1 mM and constant.

van Eerd and Takahashi, 1976; Potter et al., 1977; and Robertson et al., 1981.

Kushmerick et al., 1992, and Racay et al., 2006.

Heizmann et al., 1982, and Leberer and Pette, 1986.

Robertson et al., 1981.

Ferguson and Franzini-Armstrong, 1988.

Baylor and Hollingworth, 2007, and Hollingworth et al., 2012.

Comparison of Δ[Ca2+] in slow-twitch and fast-twitch mouse fibers

Several different Ca2+ indicators have been used in skeletal muscle fibers in pursuit of a reliable method to measure the amplitude and time course of Δ[Ca2+] (Baylor and Hollingworth, 2011). Furaptra (also called mag-fura-2; Raju et al., 1989), which is a fluorescent Ca2+ indicator of the “tricarboxylate” family, has yielded some of the most useful quantitative information to date about Δ[Ca2+] (Konishi et al., 1991; Delbono and Stefani, 1993; Claflin et al., 1994; Delbono and Meissner, 1996; Hollingworth et al., 1996, 2008, 2012; Rome et al., 1996; Baylor and Hollingworth, 2003, 2007; Harwood et al., 2011). Furaptra’s chromophore, which is similar to that of fura-2 (Grynkiewicz et al., 1985), undergoes a large change in both absorbance and fluorescence upon binding Ca2+. In contrast with fura-2, furaptra has a low affinity for Ca2+, with a dissociation constant (KD,Ca) that is estimated to be 96 µM in the myoplasmic environment (Konishi et al., 1991; Baylor and Hollingworth, 2003). Because furaptra’s low affinity is achieved by virtue of a fast off-rate constant for Ca2+, which is probably at least several thousand per second in myoplasm at 20°C (Naraghi, 1997; Baylor and Hollingworth, 1998), furaptra should always respond rapidly to changes in [Ca2+]. These features make furaptra a very valuable indicator for measuring the large and brief Δ[Ca2+] that occurs in skeletal muscle fibers during AP-evoked activity. Some valuable kinetic information about Δ[Ca2+] in mammalian fibers has also been obtained with mag-fluo-4 (Capote et al., 2005; Calderón et al., 2009, 2010; Hollingworth et al., 2012), another low-affinity tricarboxylate fluorescent indicator that reacts rapidly with Ca2+ (Hollingworth et al., 2009).

Fig. 2 A shows examples of indicator-related Ca2+ signals (top two rows of traces) and tension responses (third row of traces) recorded in mouse soleus fibers that were activated by a single AP (dashed traces) or a brief high-frequency train of APs (five APs at 67 Hz; continuous traces). Fig. 2 B shows similar recordings in EDL fibers. The traces are averaged results from several similar experiments in which one fiber within a bundle of intact fibers was micro-injected with the membrane-impermeant (K+ salt) form of either furaptra or mag-fluo-4; thus, only one fiber contributed to the change in fluorescence (ΔF) in each experiment. Based on the markedly different speeds of the contractile responses in Fig. 2 (see legend), the fibers from soleus and EDL muscles have been classified as slow-twitch and fast-twitch, respectively; biochemical properties such as the myosin heavy chain isoform were not determined.

Figure 2.

Ca2+-related fluorescence signals and tension responses in slow-twitch and fast-twitch mouse fibers stimulated by APs (16°C). (A and B) Bundles of fibers from soleus muscles (A) and EDL muscles (B) of adult mice, age 7–14 wk, were isolated by dissection, and one fiber within the bundle was micro-injected with either furaptra or mag-fluo-4. ΔF/FR and tension responses were recorded in response to external stimuli that elicited either one AP (dashed traces) or five APs at 67 Hz (continuous traces); the responses to both types of stimulation are shown superimposed. All traces represent averaged responses from several experiments in which movement artifacts in the ΔF/FR traces were small. All fast-twitch experiments were performed on fibers injected with furaptra (seven fibers in the measurements with a single AP, four in the measurements with five APs); the fast-twitch fibers are from a region of the EDL muscle thought to be rich in parvalbumin (see Baylor and Hollingworth, 2011). The slow-twitch experiments with a single AP are from 11 experiments (five with furaptra and six with mag-fluo-4); those with five APs are from three experiments (one with furaptra and two with mag-fluo-4). Before inclusion in the average responses, the mag-fluo-4 traces were scaled so that their amplitude matched that expected with furaptra (Hollingworth et al., 2009, 2012). ΔF/FR, in furaptra fluorescence units, was converted to ΔfCaD and Δ[Ca2+] as described in the text. For both fiber types, a value of 1 on the tension calibration bar corresponds to the peak tension response with the five-AP stimulus. With one AP, the time of peak and FDHM of the slow-twitch tension response are 192 and 1,001 ms, respectively; the corresponding values for the fast-twitch tension response are 55 and 215 ms; with the five-AP stimulus, the values are 229 and 993 ms for the slow-twitch response and 107 and 236 ms for the fast-twitch response. (C and D) The top traces show simulated furaptra ΔfCaD waveforms (noise-free traces) and measured ΔfCaD waveforms (traces with noise). The simulated traces were obtained with the multi-compartment model described in Fig. 3 and the text. In the model, the reaction scheme between Ca2+ and furaptra includes two pools of furaptra, protein-free and protein-bound (Baylor and Hollingworth, 2007); the protein-bound pool corresponds to the indicator molecules that appear to be bound to myoplasmic constituents of low mobility (Konishi et al., 1991). The bottom traces show the SR Ca2+ fluxes used to drive the simulations. Even though the measurements were made with the fibers stretched to a long sarcomere length (mean value, 3.7 µm), the experimental ΔfCaD traces may still be slightly contaminated with a small movement artifact.

Figure 3.

Geometry of the multi-compartment reaction–diffusion model used to simulate myoplasmic Ca2+ movements. The myoplasm of a half-sarcomere of one myofibril is divided into 18 equal-volume radially symmetric compartments (six longitudinal by three radial) at a half-sarcomere length of 2 µm. Ca2+ release (large downward arrow) occurs into the compartment located at the periphery of the myofibril about ∼0.5 µm from the z-line, the location of the triadic junctions in mammalian fibers (Smith, 1966; Eisenberg, 1983; Brown et al., 1998; Gómez et al., 2006). Ca2+ pumping (small upward arrows) is restricted to the six outermost compartments at the periphery of the myofibril (the location of the SR membrane), troponin is restricted to the nine compartments located within 1 µm of the z-line (the location of the thin filaments), and free Ca2+ and the soluble Ca2+ buffers (e.g., ATP, parvalbumin, and the Ca2+ indicator) have access to all compartments via diffusion. The buffer concentrations and reaction rate constants vary with the fiber type and compartment location (see text and Table 2). A set of first-order differential equations is specified for each compartment to describe the Ca2+ concentration changes that take place in that compartment caused by release, binding, pumping, and/or diffusion across compartment boundaries. The full set of equations can be integrated simultaneously once the SR Ca2+ release flux is specified.

In the recordings in Fig. 2, the indicator’s fluorescence was collected from a length of fiber of ∼300 µm; thus, the fluorescence signals represent a spatial average from many sarcomeres. To ensure that the fluorescence signals were minimally influenced by movement artifacts, the fibers were stretched to a long sarcomere length (average value of 3.7 µm), which greatly reduces the fiber’s contractile response. The raw indicator measurements (topmost traces in Fig. 2) are shown in units of furaptra’s ΔF/FR, the change in the indicator’s fluorescence divided by its resting fluorescence. The second row of traces shows the ΔF/FR waveforms transformed to units of Δ[Ca2+] based on a two-step procedure (Baylor and Hollingworth, 2003). In the first step, the ΔF/FR traces were scaled by the factor −1.07 (a calibration constant specific to furaptra fluorescence measured with excitation wavelengths of 410 ± 20 nm and emission wavelengths of 530 ± 60 nm), which converts the fluorescence change to the underlying change in the fraction of furaptra bound with Ca2+, ΔfCaD. The ΔfCaD traces are shown in Fig. 2 (C and D; traces with noise, which are described further below). In the second step, the ΔfCaD traces were converted to Δ[Ca2+] with the equation for 1:1 binding between Ca2+ and a low-affinity indicator: Δ[Ca2+] = KD,Ca · ΔfCaD/(1 − ΔfCaD). As described in the next section, a more accurate estimate of spatially averaged Δ[Ca2+] can be obtained from the ΔfCaD traces with a different analytical approach.

In Fig. 2 A, Δ[Ca2+] during a single twitch in slow-twitch fibers has a peak of ∼10 µM and a full duration at half-maximum (FDHM) of ∼8 ms; the corresponding values in fast-twitch fibers are ∼18 µM and ∼5 ms, respectively (Fig. 2 B). These differences between fiber types are highly reproducible and are diagnostic of the fiber type (Hollingworth et al., 2012). During the high-frequency stimulation (five APs at 67 Hz), peak Δ[Ca2+] in slow-twitch fibers rises ∼20% above that elicited by the first AP; in contrast, in fast-twitch fibers, the peaks of Δ[Ca2+] elicited by subsequent APs have about the same amplitude as that elicited by the first AP. In both fiber types, the rate of decay of Δ[Ca2+] from its fifth peak is much smaller than that from its first peak. As described below, the reason for the progressive decrease in this decay rate with later APs in the train is that there is a progressive increase in the occupancy of the myoplasmic Ca2+ buffers with Ca2+ and, consequently, a decrease in the rate at which the buffers bind Ca2+ (Melzer et al., 1987; Hollingworth et al., 1996, 2012; Carroll et al., 1997; Baylor and Hollingworth, 2007).

Simulation of myoplasmic Ca2+ movements with a reaction–diffusion model

Given the spatially averaged estimate of Δ[Ca2+] obtained from the furaptra ΔfCaD signal (Fig. 2), a spatially averaged estimate of the change in concentration of Ca2+ bound to the sites on any particular myoplasmic Ca2+ buffer can be calculated if the buffer concentration is known and a reaction scheme and reaction rate constants are available for the site (e.g., Robertson et al., 1981; Baylor et al., 1983; Cannell and Allen, 1984; Hollingworth et al., 2006; see also Fig. 5, described below). Table 2 lists the concentrations of the buffer sites considered in the calculations described here and in the next section; the reaction schemes and rate constants are given in Baylor and Hollingworth (2007) and Hollingworth et al. (2012). If the calculated changes in concentration of Ca2+ bound to these various buffers are summed and added to that of Δ[Ca2+] plus Δ[CaDye], one obtains a spatially averaged estimate of the total concentration of Ca2+ released by the SR. A more accurate estimate of Δ[Ca2+] and the other Ca2+ changes, however, can be obtained if the ΔfCaD signal from the indicator is simulated with a multi-compartment model of the myoplasmic Ca2+ movements based on a reaction–diffusion model (Cannell and Allen, 1984; Baylor and Hollingworth, 1998, 2007; Jiang et al., 1999; Novo et al., 2003; Hollingworth et al., 2012). Models of this type, in addition to using information about the reactions of Ca2+ with its buffers, incorporate information about (a) the dimensions of the sarcomere; (b) the locations of the Ca2+ release sites, the Ca2+ buffers, and Ca2+ pumping; and (c) the myoplasmic diffusion coefficients of the mobile myoplasmic constituents (free Ca2+, mobile buffers such as ATP, etc.). Calculations with these models reveal steep gradients in Δ[Ca2+] within the sarcomere. One consequence of these gradients is that the degree of saturation of the buffers with Ca2+ varies markedly at different times and sarcomeric locations. In the case of furaptra, this variation introduces error into the estimation of Δ[Ca2+] from the measured ΔfCaD signal and the 1:1 binding equation (Hirota et al., 1989; Baylor and Hollingworth, 1998, 2007). By taking into account the Δ[Ca2+] gradients, reaction–diffusion models are thought to yield more accurate estimates of spatially averaged Δ[Ca2+], of the amplitude and time course with which Ca2+ binds to the myoplasmic Ca2+ buffers, and of the amount of SR Ca2+ release.

Figure 5.

Comparison of the modeled Ca2+–troponin reactions in slow-twitch and fast-twitch fibers. (A and B) Schemes for the reaction of Ca2+ with the regulatory sites on slow-twitch troponin (A) and fast-twitch troponin (B). The text gives the values of the reaction rate constants. (C and D) Theoretical single-compartment responses of the slow-twitch (C) and fast-twitch (D) troponin regulatory sites at 16°C to a pulse of [Ca2+] that starts at 50 nM, rises to 5 µM for 50 ms, and returns to 50 nM; the response in D considers only the doubly occupied troponin state in B. The total site concentrations are 120 µM (C) and 240 µM (D).

Fig. 3 diagrams the multi-compartment model that we have used in the analysis of the furaptra Ca2+ measurements in mouse fibers. The myoplasmic water volume of a half-sarcomere of one myofibril is considered to be the basic unit of interest, and this volume is divided into 18 compartments of equal size. The buffer concentrations, as well as the locations of SR Ca2+ release and pumping, are assigned appropriately to the various compartments (see legend of Fig. 3 and Table 2). To relate the calculations with the model to the furaptra measurements, the model is driven by an assumed waveform for the SR Ca2+ release flux. The time course of this flux is estimated from the single-compartment (spatially averaged) version of the model (Baylor et al., 1983; Hollingworth et al., 1996), for which the release flux is fully determined by Δ[Ca2+] and the specified properties of the Ca2+ buffers. The amplitude of the flux waveform in the multi-compartment model is then adjusted iteratively until the amplitude of the simulated, spatially averaged ΔfCaD waveform matches that of the measured ΔfCaD waveform. The overall success of this approach is evaluated by a comparison of the time courses of the simulated and measured ΔfCaD waveforms. If these time courses are in agreement, it is reasonable to believe that the simulations give an approximate description of the underlying Ca2+ movements. In Fig. 2 (C and D), the top superimposed traces show examples of such comparisons, and the bottom traces show the SR Ca2+ release fluxes used to drive the simulations. The good agreement between the simulated ΔfCaD waveforms (noise-free traces) and the measured ΔfCaD waveforms (traces with noise) indicates that this strategy is successful in both fiber types.

Comparison of myoplasmic Ca2+ movements in slow-twitch and fast-twitch fibers

Fig. 4 summarizes some results from these multi-compartment simulations. Responses in slow-twitch and fast-twitch fibers are shown as continuous and dashed traces, respectively. Fig. 4 A compares responses to a single AP, and B compares them to five APs at 67 Hz. For simplicity in the presentation, all waveforms in Fig. 4 have been spatially averaged; i.e., each trace was calculated as the average of the corresponding changes in the 18 compartments of the model. (Examples of waveforms calculated for the individual compartments are given in Baylor and Hollingworth, 2007). The traces in Fig. 4 are (from bottom to top): (a) Δ[Ca2+]; (b) the change in concentration of Ca2+ bound to the regulatory sites on troponin (Δ[CaTrop]); (c) the fractional occupancy of the troponin regulatory sites with Ca2+ (fCaTrop, which, for this variable, includes the resting occupancy); (d) the change in concentration of Ca2+ bound to the sites on the slowly reacting buffers (Δ[CaSlowB], which includes Ca2+ binding to parvalbumin in fast-twitch fibers and to the troponin-nonspecific sites in slow-twitch fibers; see Table 2); (e) the change in concentration of Ca2+ bound to the SR Ca2+ pumps (Δ[CaPump]); (f) the change in concentration of Ca2+ returned to the SR as a result of Ca2+ pumping (Δ[CaPumped]); and (g) the total concentration of released Ca2+ (Δ[CaT]). Not shown, but included in the calculation of Δ[CaT], are Δ[CaATP] and Δ[CaDye] (see legend of Fig. 4). Except for fCaTrop, these traces are shown in concentration units referred to the myoplasmic water volume.

Figure 4.

Comparison of Ca2+ movements elicited by APs in slow-twitch and fast-twitch fibers (16°C). (A). Simulated Ca2+ movements in slow-twitch and fast-twitch fibers (continuous and dashed traces, respectively) in response to a single AP initiated at zero time. fCaTrop, the fractional occupancy of the troponin regulatory sites with Ca2+, includes the resting occupancy (0.071 in slow-twitch fibers and 0.004 in fast-twitch fibers) plus the change in occupancy during activity. See text for definitions of the other traces. Δ[CaT] is the sum of the five lower traces, excluding fCaTrop, plus Δ[CaATP] and Δ[CaDye] (the change in the concentration of Ca2+ bound to the indicator). Δ[CaATP] is approximately equal to 2.2 times Δ[Ca2+] in slow-twitch fibers and 3.4 times Δ[Ca2+] in fast-twitch fibers. Δ[CaDye] is relatively small (= ΔfCaD × [DyeT], where [DyeT], the total indicator concentration in myoplasm, is 100 µM). The units of all traces except fCaTrop are spatially averaged concentrations referred to the myoplasmic water volume. (B). As in A, but for five APs at 67 Hz.

In Fig. 4 A, peak Δ[CaT] in response to one AP in a slow-twitch fiber is 107 µM, which is achieved with a peak SR Ca2+ flux of 57 µM/ms and a flux FDHM of 1.7 ms (see also Fig. 2 C); the corresponding values for a fast-twitch fiber are 349 µM, 205 µM/ms, and 1.6 ms (see also Fig. 2 D). The modeled spatially averaged Δ[Ca2+] in slow-twitch fibers has a peak amplitude of 7.6 µM and an FDHM of 4.9 ms. As expected from Fig. 2, the corresponding values in fast-twitch fibers are larger and briefer, 16.1 µM and 3.7 ms. These properties of spatially averaged Δ[Ca2+] obtained with the multi-compartment models are thought to be more accurate than those obtained above with the single-compartment analysis (Fig. 2, A and B). The multi-compartment estimates of the FDHM of Δ[Ca2+] are substantially briefer: 4.9 versus 8 ms for slow-twitch fibers and 3.7 versus 5 ms for fast-twitch fibers. In both fiber types, peak Δ[Ca2+] is a small fraction of peak Δ[CaT] (<0.1); i.e., the great majority of the released Ca2+ is rapidly bound by the Ca2+ buffers, with only a small fraction remaining free.

In Fig. 4 A, peak Δ[CaTrop] shows a large difference between the two fiber types: 86 µM for slow-twitch fibers and 227 µM for fast-twitch fibers. This difference arises primarily because fast-twitch troponin has two Ca2+ regulatory sites per molecule, whereas slow-twitch troponin has only one (van Eerd and Takahashi, 1976; Potter et al., 1977; Robertson et al., 1981; Sweeney et al., 1990). If, however, the concentration of Ca2+ bound to the troponin regulatory sites is normalized by the total concentration of these sites (120 µM in slow-twitch fibers and 240 µM in fast-twitch fibers; Table 2), peak fCaTrop is similar in slow-twitch and fast-twitch fibers: 0.79 and 0.95, respectively. Δ[CaSlowB] is also very different in the two fiber types, with values of 9 and 146 µM, respectively, at t = 50 ms. This difference reflects the large concentration of parvalbumin that is found in mouse fast-twitch glycolytic fibers compared with slow-twitch fibers, which have little or no parvalbumin (compare Table 2). Δ[CaPump] is also quite different in the two fiber types, with 10 µM in slow-twitch fibers and 28 µM in fast-twitch fibers at t = 50 ms. The primary reason is that the concentration of Ca2+ pump molecules is larger in fast-twitch fibers: 120 versus 48 µM (Table 2). A second reason is that, because Δ[Ca2+] is larger in fast-twitch fibers, the fractional occupancy of the pump sites with Ca2+ is larger. (Note: fast-twitch and slow-twitch fibers have different proportions of the two main isoforms of the Ca2+ pump molecules found in skeletal muscle, denoted SERCA1 and SERCA2a, but the functional properties of the two isoforms are similar; Lytton et al., 1992; Wu and Lytton, 1993.) Because of the larger amplitude of Δ[CaPump] in fast-twitch fibers, Δ[CaPumped] is also larger, with values at t = 50 ms of 3.7 µM for fast-twitch fibers and 1.1 µM for slow-twitch fibers. Ca2+ pumping is a relatively slow process, with an estimated turnover rate at saturating [Ca2+] of 3–4 s−1 (16°C; Inesi et al., 1982; Young et al., 2003); thus, the resequestration of the released Ca2+ into the SR is not expected to be completed until several seconds after stimulation.

Fig. 4 B, which compares responses for the case of five APs at 67 Hz, reveals several new features. First, in both fiber types, the increments in Δ[CaT] caused by the second through fifth APs are much smaller than that caused by the first. Considered as a fraction of the first release, the values vary from 0.30 to 0.13 in slow-twitch fibers and from 0.26 to 0.15 in fast-twitch fibers (see also bottommost traces in Fig. 2). These marked reductions in release are thought to be caused primarily by the process of Ca2+ inactivation of Ca2+ release, in which a rise in [Ca2+] is sensed by the release system to inhibit release (Baylor et al., 1983; Schneider and Simon, 1988; Baylor and Hollingworth, 1988, 2007; Jong et al., 1995). The molecule(s) that mediates this negative feedback remains to be fully identified, although an interaction between calmodulin and an identified site on the SR release channel may be involved (Yamaguchi et al., 2011). This feedback has an obvious physiological function, as it prevents [Ca2+] from rising to higher levels than is required to activate the myofilaments, thereby avoiding: (a) an unnecessary delay in fiber relaxation after cessation of APs, (b) extra expenditure of ATP by the Ca2+ pumps to resequester Ca2+, and (c) potential damage to the mitochondria caused by excessive uptake of Ca2+.

A second noteworthy feature in Fig. 4 B is that peak fCaTrop is close to 1 in both fiber types: 0.93 in slow-twitch fibers and 0.96 in fast-twitch fibers. Thus, in both fiber types, the models predict that Δ[Ca2+] during repetitive stimulation rises to a level at which activation of the myofilaments is not limited by [Ca2+]. The conclusion that the troponin regulatory sites are rapidly and nearly completely occupied with Ca2+ in both fiber types while the rate of tension development is much slower in slow-twitch than in fast-twitch fibers (see tension traces in Fig. 2, A and B) is consistent with observations in chemically skinned rodent fibers that the rate at which myosin produces force at saturating [Ca2+] is inherently slower with slow-twitch than with fast-twitch myosin (Metzger and Moss, 1990). Because the troponin regulatory sites and the sites on the other buffers (e.g., “SlowB”) become increasingly occupied with Ca2+ during the high-frequency stimulus, the ability of these sites to bind newly released Ca2+ decreases; as a result, there is a progressive decrease in the rate of decay of Δ[Ca2+] from its individual peaks during the train. Consistent with the longer lasting Δ[Ca2+] waveform, Δ[CaSlowB], Δ[CaPump], and Δ[CaPumped] all rise to substantially higher levels than those that occur with a single AP.

Functional differences between Ca2+ movements in slow-twitch and fast-twitch fibers

As noted above, the calculations indicate that the troponin regulatory sites rapidly bind a large fraction of the released Ca2+ in both slow-twitch and fast-twitch fibers. The fraction is smaller in fast-twitch fibers, however, because fast-twitch fibers have relatively larger concentrations of sites on the other buffers (Table 2), which compete with troponin for Ca2+. Thus, to achieve full occupancy of their regulatory sites, fast-twitch fibers must release substantially more Ca2+ than is expected from their twofold larger concentration of regulatory sites. The ability of fast-twitch fibers to release three- to fourfold larger amounts of Ca2+ than slow-twitch fibers (topmost traces in Fig. 4) at three- to fourfold larger rates (bottommost traces in Fig. 2) is achieved by means of comparable increases in the fractional volume of their SR (Eisenberg, 1983) and, concomitantly, in the number of the SR Ca2+ release channels and their controlling DHPR molecules (Hollingworth and Marshall, 1981; Lamb and Walsh, 1987; Franzini-Armstrong et al., 1988; Delbono and Meissner, 1996). The functional advantage of the large concentration of the slow-buffer sites in fast-twitch fibers is that a substantial fraction of the Ca2+ that is initially bound by troponin is subsequently bound by these sites (Gillis et al., 1982; Baylor et al., 1983; see also Fig. 4). By serving as a temporary myoplasmic storage site for released Ca2+, the slow buffers contribute importantly to the more rapid decline of Δ[Ca2+] that occurs in fast-twitch fibers, which helps to speed their contractile relaxation. An example of this slow-buffer function taken to an extreme occurs in the superfast swimbladder fibers of the toadfish (Opsanau tau), where the concentration of parvalbumin is nearly fourfold larger than that in mouse fast-twitch fibers (2.7 vs. 0.75 mM; Hamoir et al., 1980; Appelt et al., 1991; Tikunov and Rome, 2009), and the unusually brief Δ[Ca2+] elicited by an AP in swimbladder fibers can be maintained for extended periods at stimulation frequencies as high as 100 Hz (16°C; Rome et al., 1996; Harwood et al., 2011).

Fig. 5 illustrates several relevant differences between the modeled troponin reactions in slow-twitch and fast-twitch fibers. Ca2+ binds to the single regulatory site on each slow-twitch troponin molecule with a simple 1:1 reaction (Fig. 5 A), whereas Ca2+ binds to the two sites on each fast-twitch troponin molecule in a two-step reaction (Fig. 5 B). The rate constants for the slow-twitch reaction, which were estimated from information in Davis et al. (2007) in combination with comparisons of the type shown in Fig. 2 C (Hollingworth et al., 2012), are 0.4 × 108 M−1 s−1 for the on-rate (k+1) and 26 s−1 for the off-rate (k−1) (16°C); the corresponding value of KD (= k−1/k+1) is 0.65 µM. In contrast, with fast-twitch troponin, Ca2+ is thought to react with the two regulatory sites with positive cooperativity (Fuchs and Bayuk, 1976; Grabarek et al., 1983; see also Baylor et al., 2002). To achieve cooperativity in the modeled troponin reaction, the first Ca2+ ion binds with low affinity and the second with high affinity. The reaction rate constants are k+1 = 1.77 × 108 M−1 s−1 and k−1 = 1,542 s−1 (KD,1 = 8.72 µM), and k+2 = 0.885 × 108 M−1 s−1 and k-2 = 17.1 s−1 (KD,2 = 0.193 µM), respectively (Baylor et al., 2002; Davis, J.P., S.B. Tikunova, D.R. Swartz, and J.A. Rall. 2004. Biophysical Society 48th Annual Meeting. Abstr. 1135; Hollingworth et al., 2006). The steady-state level of free [Ca2+] that gives half-occupancy of these sites (Ca2+50) is 1.3 µM (= sqrt (KD,1 × KD,2)), which is twice Ca2+50 for the slow-twitch troponin reaction. These features are consistent with observations on mammalian skinned fibers that the tension–pCa curve in fast-twitch fibers is right-shifted with respect to that in slow-twitch fibers and has a steeper dependence on [Ca2+] (Kerrick et al., 1976; Stephenson and Williams, 1981; Ruff, 1989; Danieli-Betto et al., 1990). For example, in rat fibers, the [Ca2+] levels for half-activation of tension are ∼1.2 and ∼0.6 µM, respectively, for fast-twitch and slow-twitch fibers (22°C; sarcomere length of 2.7–2.8 µm; Stephenson and Williams, 1981).

Fig. 5 (C and D) presents a theoretical single-compartment analysis of some of the kinetic consequences of the modeled troponin reactions in the two fiber types. As revealed by this analysis, the more rapid rate of dissociation of Ca2+ from the regulatory sites that is observed in Fig. 4 B in fast-twitch fibers compared with slow-twitch fibers is explained by the more rapid decline of Δ[Ca2+] in fast twitch fibers, not by an inherently faster rate of Ca2+ dissociation from the fast-twitch regulatory sites. In Fig. 5 C, the bottom trace shows a pulse of [Ca2+], which starts at 50 nM (= [Ca2+]R), rises in a square-wave fashion to 5 µM, where it remains for 50 ms, and then returns to the resting level, again in a square-wave fashion. The top trace shows the calculated Δ[CaTrop] response for slow-twitch troponin (compare Fig. 5 A), which rises exponentially from the resting level of 8.6 µM with a rate of 226 s−1 to 106 µM (= 88.5% of saturation), and then decays back to the resting level with a rate of 28 s−1. Fig. 5 D shows a similar calculation for the doubly occupied Ca2+ state of fast-twitch troponin (Ca2Trop in Fig. 5 B), which is the troponin state that likely controls activation of the myofilaments. Here, the concentration of Ca2+ that is bound to the active state (= 2 × [Ca2Trop]) is 0.4 µM at rest and rises to 217 µM at peak (= 90.5% of saturation). Because the fast-twitch troponin reaction involves two-steps, neither the “on” nor “off” response exactly follows a single-exponential time course. Nevertheless, both responses are very well fitted by single-exponential time courses (not depicted), with effective on- and off-rates of 153 and 17.2 s−1, respectively. The effective off-rate in Fig. 5 D is smaller than that in Fig. 5 C, although Ca2+50 is larger in fast-twitch than in slow-twitch fibers (1.3 vs. 0.65 µM). A significant advantage of the two-step cooperative troponin reaction (Fig. 5 B) is that it allows a more switchlike (all-or-nothing) control of the actomyosin system during activity. For example, the fraction of the troponin molecules in the activated state (the Ca2Trop state) is negligible at rest in fast-twitch fibers (0.001), whereas, in slow-twitch fibers, the resting fraction in the CaTrop state (0.071) is not entirely negligible. The latter value would be expected to lead to a slight activation of the actomyosin system in slow-twitch fibers at rest and could contribute to the higher rate of resting energy consumption that has been reported in mouse soleus muscle compared with EDL muscle (Decrouy et al., 1993; Norris et al., 2010).

Calcium movements in slow-twitch and fast-twitch fibers at higher temperatures

The results described in connection with Figs. 2, 4, and 5 were obtained at 16°C, an experimental temperature that facilitates survival of the fibers outside the animal. To assess the functioning of the fibers at more physiological temperatures, some measurements like those in Fig. 2 have also been performed at higher temperatures (22, 28, and 35°C; Hollingworth et al., 1996, 2012; Baylor and Hollingworth, 2003). For example, in fibers at 28°C stimulated by a single AP, the following parameter values were reported for slow-twitch (soleus) fibers versus fast-twitch (EDL) fibers: (a) time of peak and FDHM of twitch tension: 97 and 480 ms versus 23 and 97 ms, respectively; (b) peak amplitude and FDHM of Δ[Ca2+]: 9 µM and 4.4 ms versus 21 µM and 2.0 ms; and (c) the amount and peak rate of SR Ca2+ release: 120 µM and 105 µM/ms versus 358 µM and 473 µM/ms (Baylor and Hollingworth, 2003). These relative differences between fiber types are quite similar to those reported above at 16°C, and it is reasonable to believe that similar relative differences apply under fully physiological conditions.

Summary

The waveforms in Fig. 4 provide detailed estimates of the major myoplasmic Ca2+ movements that take place in mouse slow-twitch and fast-twitch fibers during twitches and brief tetani. Although it is too much to expect that all aspects of these comparisons are error free, we believe that the basic features shown in Fig. 4 provide a reasonable working model of Ca2+ signaling in the two fiber types. Hopefully, future measurements will permit corrections and refinements of the estimated Ca2+ movements. The analysis reveals that, for comparable stimuli, slow-twitch fibers release ∼0.3 times the amount of SR Ca2+ as fast-twitch fibers. A corollary prediction is that, to resequester the released Ca2+, the SR Ca2+ pumps in slow-twitch fibers consume ∼0.3 times the ATP as that in fast-twitch fibers. This factor is similar to that which characterizes the total ATP consumption during brief contractions of slow-twitch versus fast-twitch fibers (20°C; Crow and Kushmerick, 1982); in this case, however, the ATP consumption is mainly caused by myosin, with a minor fraction caused by the SR Ca2+ pumps (Barclay et al., 2007).

One interesting avenue for future investigation concerns the time course with which Δ[Ca2+] declines to zero after stimuli like those in Fig. 2. In both slow-twitch and fast-twitch fibers, our measurements with high-affinity fluorescent Ca2+ indicators reveal that this decline—and hence the time of pumping Ca2+ back into the SR—continues for some seconds (16°C). It remains to be determined how well the model predictions of Δ[Ca2+] on this slow time scale agree with these measurements. Another avenue concerns the properties of Δ[Ca2+] and the associated myoplasmic Ca2+ movements in fibers of the fast-twitch oxidative subtype (fatigue-resistant fast-twitch fibers; Berchtold et al., 2000; Schiaffino and Reggiani, 2011), which are not considered here. In mice, these fibers make up a significant fraction of the fibers in both soleus and EDL muscles (e.g., Crow and Kushmerick, 1982; Haida et al., 1989; Maréchal et al., 1995; Allen et al., 2001; Asmussen et al., 2003), and they remain to be studied with the methods described in this article.