The effect of sarcomere length on triad location in intact feline caudofeomoralis muscle fibres

Force variability is mostly not motor noise: Theoretical implications for motor control

Variability in muscle force is a hallmark of healthy and pathological human behavior. Predominant theories of sensorimotor control assume 'motor noise' leads to force variability and its 'signal dependence' (variability in muscle force whose amplitude increases with intensity of neural drive). Here, we demonstrate that the two proposed mechanisms for motor noise (i.e. the stochastic nature of motor unit discharge and unfused tetanic contraction) cannot account for the majority of force variability nor for its signal dependence. We do so by considering three previously underappreciated but physiologically important features of a population of motor units: 1) fusion of motor unit twitches, 2) coupling among motoneuron discharge rate, cross-bridge dynamics, and muscle mechanics, and 3) a series-elastic element to account for the aponeurosis and tendon. These results argue strongly against the idea that force variability and the resulting kinematic variability are generated primarily by 'motor noise.' Rather, they underscore the importance of variability arising from properties of control strategies embodied through distributed sensorimotor systems. As such, our study provides a critical path toward developing theories and models of sensorimotor control that provide a physiologically valid and clinically useful understanding of healthy and pathologic force variability.

Muscle and Limb Mechanics

Understanding of the musculoskeletal system has evolved from the collection of individual phenomena in highly selected experimental preparations under highly controlled and often unphysiological conditions. At the systems level, it is now possible to construct complete and reasonably accurate models of the kinetics and energetics of realistic muscles and to combine them to understand the dynamics of complete musculoskeletal systems performing natural behaviors. At the reductionist level, it is possible to relate most of the individual phenomena to the anatomical structures and biochemical processes that account for them. Two large challenges remain. At a systems level, neuroscience must now account for how the nervous system learns to exploit the many complex features that evolution has incorporated into muscle and limb mechanics. At a reductionist level, medicine must now account for the many forms of pathology and disability that arise from the many diseases and injuries to which this highly evolved system is inevitably prone. © 2017 American Physiological Society. Compr Physiol 7:429-462, 2017.

Comparison of myoplasmic calcium movements during excitation-contraction coupling in frog twitch and mouse fast-twitch muscle fibers

Single twitch fibers from frog leg muscles were isolated by dissection and micro-injected with furaptra, a rapidly responding fluorescent Ca(2+) indicator. Indicator resting fluorescence (FR) and the change evoked by an action potential (ΔF) were measured at long sarcomere length (16°C); ΔF/FR was scaled to units of ΔfCaD, the change in fraction of the indicator in the Ca(2+)-bound form. ΔfCaD was simulated with a multicompartment model of the underlying myoplasmic Ca(2+) movements, and the results were compared with previous measurements and analyses in mouse fast-twitch fibers. In frog fibers, sarcoplasmic reticulum (SR) Ca(2+) release evoked by an action potential appears to be the sum of two components. The time course of the first component is similar to that of the entire Ca(2+) release waveform in mouse fibers, whereas that of the second component is severalfold slower; the fractional release amounts are ~0.8 (first component) and ~0.2 (second component). Similar results were obtained in frog simulations with a modified model that permitted competition between Mg(2+) and Ca(2+) for occupancy of the regulatory sites on troponin. An anatomical basis for two release components in frog fibers is the presence of both junctional and parajunctional SR Ca(2+) release channels (ryanodine receptors [RyRs]), whereas mouse fibers (usually) have only junctional RyRs. Also, frog fibers have two RyR isoforms, RyRα and RyRβ, whereas the mouse fibers (usually) have only one, RyR1. Our simulations suggest that the second release component in frog fibers functions to supply extra Ca(2+) to activate troponin, which, in mouse fibers, is not needed because of the more favorable location of their triadic junctions (near the middle of the thin filament). We speculate that, in general, parajunctional RyRs permit increased myofilament activation in fibers whose triadic junctions are located at the z-line.

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²⁺) from the sarcoplasmic reticulum. Experiments on individual mouse muscle fibers micro-injected with a rapidly responding fluorescent Ca²⁺ indicator dye reveal that the amount of Ca²⁺ released is three- to fourfold larger in fast-twitch fibers than in slow-twitch fibers, and the proportion of the released Ca²⁺ that binds to troponin to activate contraction is substantially smaller.

Measurement and simulation of myoplasmic calcium transients in mouse slow-twitch muscle fibres

Bundles of intact fibres from soleus muscles of adult mice were isolated by dissection and one fibre within a bundle was micro-injected with either furaptra or mag-fluo-4, two low-affinity rapidly responding Ca(2+) indicators. Fibres were activated by action potentials to elicit changes in indicator fluorescence (ΔF), a monitor of the myoplasmic free Ca(2+) transient ([Ca(2+)]), and changes in fibre tension. All injected fibres appeared to be slow-twitch (type I) fibres as inferred from the time course of their tension responses. The full-duration at half-maximum (FDHM) of ΔF was found to be essentially identical with the two indicators; the mean value was 8.4 ± 0.3 ms (±SEM) at 16°C and 5.1 ± 0.3 ms at 22°C. The value at 22°C is about one-third that reported previously in enzyme-dissociated slow-twitch fibres that had been AM-loaded with mag-fluo-4: 12.4 ± 0.8 ms and 17.2 ± 1.7 ms. We attribute the larger FDHM in enzyme-dissociated fibres either to an alteration of fibre properties due to the enzyme treatment or to some error in the measurement of ΔF associated with AM loading. ΔF in intact fibres was simulated with a multi-compartment reaction-diffusion model that permitted estimation of the amount and time course of Ca(2+) release from the sarcoplasmic reticulum (SR), the binding and diffusion of Ca(2+) in the myoplasm, the re-uptake of Ca(2+) by the SR Ca(2+) pump, and Δ[Ca(2+)] itself. In response to one action potential at 16°C, the following estimates were obtained: 107 μm for the amount of Ca(2+) release; 1.7 ms for the FDHM of the release flux; 7.6 μm and 4.9 ms for the peak and FDHM of spatially averaged Δ[Ca(2+)]. With five action potentials at 67 Hz, the estimated amount of Ca(2+) release is 186 μm. Two important unknown model parameters are the on- and off-rate constants of the reaction between Ca(2+) and the regulatory sites on troponin; values of 0.4 × 10(8) m(-1) s(-1) and 26 s(-1), respectively, were found to be consistent with the ΔF measurements.

Calcium indicators and calcium signalling in skeletal muscle fibres during excitation-contraction coupling

During excitation-contraction coupling in skeletal muscle, calcium ions are released into the myoplasm by the sarcoplasmic reticulum (SR) in response to depolarization of the fibre's exterior membranes. Ca(2+) then diffuses to the thin filaments, where Ca(2+) binds to the Ca(2+) regulatory sites on troponin to activate muscle contraction. Quantitative studies of these events in intact muscle preparations have relied heavily on Ca(2+)-indicator dyes to measure the change in the spatially-averaged myoplasmic free Ca(2+) concentration (Δ[Ca(2+)]) that results from the release of SR Ca(2+). In normal fibres stimulated by an action potential, Δ[Ca(2+)] is large and brief, requiring that an accurate measurement of Δ[Ca(2+)] be made with a low-affinity rapidly-responding indicator. Some low-affinity Ca(2+) indicators monitor Δ[Ca(2+)] much more accurately than others, however, as reviewed here in measurements in frog twitch fibres with sixteen low-affinity indicators. This article also examines measurements and simulations of Δ[Ca(2+)] in mouse fast-twitch fibres. The simulations use a multi-compartment model of the sarcomere that takes into account Ca(2+)'s release from the SR, its diffusion and binding within the myoplasm, and its re-sequestration by the SR Ca(2+) pump. The simulations are quantitatively consistent with the measurements and appear to provide a satisfactory picture of the underlying Ca(2+) movements.

Comparison of the myoplasmic calcium transient elicited by an action potential in intact fibres of mdx and normal mice

The myoplasmic free [Ca2+] transient elicited by an action potential (Delta[Ca2+]) was compared in fast-twitch fibres of mdx (dystrophin null) and normal mice. Methods were used that maximized the likelihood that any detected differences apply in vivo. Small bundles of fibres were manually dissected from extensor digitorum longus muscles of 7- to 14-week-old mice. One fibre within a bundle was microinjected with furaptra, a low-affinity rapidly responding fluorescent calcium indicator. A fibre was accepted for study if it gave a stable, all-or-nothing fluorescence response to an external shock. In 18 normal fibres, the peak amplitude and the full-duration at half-maximum (FDHM) of Delta[Ca2+] were 18.4 +/- 0.5 microm and 4.9 +/- 0.2 ms, respectively (mean +/- s.e.m.; 16 degrees C). In 13 mdx fibres, the corresponding values were 14.5 +/- 0.6 microm and 4.7 +/- 0.2 ms. The difference in amplitude is statistically highly significant (P = 0.0001; two-tailed t test), whereas the difference in FDHM is not (P = 0.3). A multi-compartment computer model was used to estimate the amplitude and time course of the sarcoplasmic reticulum (SR) calcium release flux underlying Delta[Ca2+]. Estimates were made based on several differing assumptions: (i) that the resting myoplasmic free Ca2+ concentration ([Ca2+]R) and the total concentration of parvalbumin ([Parv(T)]) are the same in mdx and normal fibres, (ii) that [Ca2+](R) is larger in mdx fibres, (iii) that [Parv(T)] is smaller in mdx fibres, and (iv) that [Ca2+]R is larger and [Parv(T)] is smaller in mdx fibres. According to the simulations, the 21% smaller amplitude of Delta[Ca2+] in mdx fibres in combination with the unchanged FDHM of Delta[Ca2+] is consistent with mdx fibres having a approximately 25% smaller flux amplitude, a 6-23% larger FDHM of the flux, and a 9-20% smaller total amount of released Ca2+ than normal fibres. The changes in flux are probably due to a change in the gating of the SR Ca2+-release channels and/or in their single channel flux. The link between these changes and the absence of dystrophin remains to be elucidated.

Simulation of Ca2+ movements within the sarcomere of fast-twitch mouse fibers stimulated by action potentials

Ca²⁺ release from the sarcoplasmic reticulum (SR) of skeletal muscle takes place at the triadic junctions; following release, Ca²⁺ spreads within the sarcomere by diffusion. Here, we report multicompartment simulations of changes in sarcomeric Ca²⁺ evoked by action potentials (APs) in fast-twitch fibers of adult mice. The simulations include Ca²⁺ complexation reactions with ATP, troponin, parvalbumin, and the SR Ca²⁺ pump, as well as Ca²⁺ transport by the pump. Results are compared with spatially averaged Ca²⁺ transients measured in mouse fibers with furaptra, a low-affinity, rapidly responding Ca²⁺ indicator. The furaptra Δf_CaD signal (change in the fraction of the indicator in the Ca²⁺-bound form) evoked by one AP is well simulated under the assumption that SR Ca²⁺ release has a peak of 200–225 μM/ms and a FDHM of ∼1.6 ms (16°C). Δf_CaD elicited by a five-shock, 67-Hz train of APs is well simulated under the assumption that in response to APs 2–5, Ca²⁺ release decreases progressively from 0.25 to 0.15 times that elicited by the first AP, a reduction likely due to Ca²⁺ inactivation of Ca²⁺ release. Recovery from inactivation was studied with a two-AP protocol; the amplitude of the second release recovered to >0.9 times that of the first with a rate constant of 7 s⁻¹. An obvious feature of Δf_CaD during a five-shock train is a progressive decline in the rate of decay from the individual peaks of Δf_CaD. According to the simulations, this decline is due to a reduction in available Ca²⁺ binding sites on troponin and parvalbumin. The effects of sarcomere length, the location of the triadic junctions, resting [Ca²⁺], the parvalbumin concentration, and possible uptake of Ca²⁺ by mitochondria were also investigated. Overall, the simulations indicate that this reaction-diffusion model, which was originally developed for Ca²⁺ sparks in frog fibers, works well when adapted to mouse fast-twitch fibers stimulated by APs.

Triadins are not triad-specific proteins: two new skeletal muscle triadins possibly involved in the architecture of sarcoplasmic reticulum

We have cloned two new triadin isoforms from rat skeletal muscle, Trisk 49 and Trisk 32, which were named according to their theoretical molecular masses (49 and 32 kDa, respectively). Specific antibodies directed against each protein were produced to characterize both new triadins. Both are expressed in adult rat skeletal muscle, and their expression in slow twitch muscle is lower than that in fast twitch muscle. Using double immunofluorescent labeling, the localization of these two triadins was studied in comparison to well-characterized proteins such as ryanodine receptor, calsequestrin, desmin, Ca(2+)-ATPase, and titin. None of these two triadins are localized within the rat skeletal muscle triad. Both are instead found in different parts of the longitudinal sarcoplasmic reticulum. We attempted to identify partners for each isoform: neither is associated with ryanodine receptor; Trisk 49 could be associated with titin or another sarcomeric protein; and Trisk 32 could be associated with IP(3) receptor. These results open further fields of research concerning the functions of these two proteins; in particular, they could be involved in the set up and maintenance of a precise sarcoplasmic reticulum structure.

Reversible vacuolation of T-tubules in skeletal muscle: mechanisms and implications for cell biology

The majority of investigations of the transverse tubules (T-system) of skeletal muscle have been devoted to their role in excitation-contraction coupling, with particular reference to contact with the sarcoplasmic reticulum and the mechanism of Ca2- release. By contrast, this review is concerned with structural and functional aspects of the vacuolation of T-tubules. It covers experimental procedures used in reversible vacuolation induced by the efflux-influx of glycerol and other small nonelectrolytes, sugars, and ions. The characteristics of the phenomenon, associated alterations in muscle function, and the swelling of analogous structures in nonmuscle cells are considered. Possible functions of reversible vacuolation in water balance, transport, membrane repair, muscle pathology, and fatigue are considered, and the potential application of reversible vacuolation in the transfection of skeletal muscle is discussed. In relation to the possible mechanisms involved in reversible vacuolation, particular attention is given to the dynamic and structural aspects of the opening and closing of T-tubules, the origin of vacuolar membranes, and the localized character of tubular swelling.

Measured and modeled properties of mammalian skeletal muscle. II. The effects of stimulus frequency on force-length and force-velocity relationships

Interactions between physiological stimulus frequencies, fascicle lengths and velocities were analyzed in feline caudofemoralis (CF), a hindlimb skeletal muscle composed exclusively of fast-twitch fibers. Split ventral roots were stimulated asynchronously to produce smooth contractions at sub-tetanic stimulus frequencies. As described previously, the peak of the sub-tetanic force-length relationship was found to shift to longer lengths with decreases in stimulus frequency, indicating a length dependence for activation that is independent of filament overlap. The sub-tetanic force-velocity (FV) relationship was affected strongly both by stimulus frequency and by length; decreases in either decreased the slope of the FV relationship around isometric. The shapes of the force transients following stretch or shortening revealed that these effects were not due to a change in the instantaneous FV relationship; the relative shape of the force transients following stretch or shortening was independent of stimulus frequency and hardly affected by length. The effects of stimulus frequency and length on the sub-tetanic FV relationship instead appear to be caused by a time delay in the length-dependent changes of activation. In contrast to feline soleus muscle, which is composed exclusively of slow-twitch fibers, CF did not yield at sub-tetanic stimulus frequencies for the range of stretch velocities tested (up to 2 L0/s). The data presented here were used to build a model of muscle that accounted well for all of the effects described. We extended our model to account for slow twitch muscle by comparing our fast-twitch model with previously published data and then changing the necessary parameters to fit the data. Our slow-twitch model accounts well for all previous findings including that of yielding.