Phosphorylation of the myosin regulatory light chain (RLC) by skeletal myosin light chain kinase (skMLCK) potentiates rodent fast twitch muscle but is an ATP-requiring process. Our objective was to investigate the effect of skMLCK-catalyzed RLC phosphorylation on the energetic cost of contraction and the contractile economy (ratio of mechanical output to metabolic input) of mouse fast twitch muscle in vitro (25°C). To this end, extensor digitorum longus (EDL) muscles from wild-type (WT) and from skMLCK-devoid (skMLCK−/−) mice were subjected to repetitive low-frequency stimulation (10 Hz for 15 s) to produce staircase potentiation of isometric twitch force, after which muscles were quick frozen for determination of high-energy phosphate consumption (HEPC). During stimulation, WT muscles displayed significant potentiation of isometric twitch force while skMLCK−/− muscles did not (i.e. 23% versus 5% change, respectively). Consistent with this, RLC phosphorylation was increased ∼3.5-fold from the unstimulated control value in WT but not in skMLCK−/− muscles. Despite these differences, the HEPC of WT muscles was not greater than that of skMLCK−/− muscles. As a result of the increased contractile output relative to HEPC, the calculated contractile economy of WT muscles was greater than that of skMLCK−/− muscles. Thus, our results suggest that skMLCK-catalyzed phosphorylation of the myosin RLC increases the contractile economy of WT mouse EDL muscle compared with skMLCK−/− muscles without RLC phosphorylation.

The striated muscle cells that form skeletal muscle are biological motors that convert chemical energy into mechanical forces to power locomotion such as running, swimming and flight. Within each sarcomere, calcium (Ca2+)-regulated interactions between myriad myosin motor molecules and their designated docking sites on sarcomeric myofilaments transduce the energy provided by the hydrolysis of adenosine triphosphate (ATP) into force and work. The ability of cycling cross-bridges to generate the inter-filament forces required for muscle shortening and thus locomotion depends upon the ability of intracellular Ca2+ signals to activate cross-bridge formation in a feed-forward manner, while the availability of ATP and/or inhibitory effect of ATP hydrolysis product accumulation may limit cross-bridge formation cycling via metabolic feed-back (Lymn and Taylor, 1971).

List of abbreviations
     
  • BTS

    N-benzyl-p-toluene sulfonamide

  •  
  • Cr

    creatine

  •  
  • −dP/dt

    rate of force relaxation

  •  
  • +dP/dt

    rate of force development

  •  
  • EA

    activation energy

  •  
  • Etotal

    total energy turnover

  •  
  • EXB

    cross-bridge energy turnover

  •  
  • EDL

    extensor digitorum longus

  •  
  • FTI

    force–time integral

  •  
  • HEPC

    high-energy phosphate consumption

  •  
  • Lo

    optimal length (for peak muscle twitch force)

  •  
  • PCr

    phosphocreatine

  •  
  • Pt

    peak force

  •  
  • RLC

    regulatory light chain

  •  
  • RT1/2

    half-relaxation time

  •  
  • skMLCK−/−

    myosin light chain kinase knockout (skeletal muscle isoform)

  •  
  • TPT

    time to peak tension

  •  
  • Vmax

    maximal shortening velocity

The intracellular Ca2+ signal that regulates cross-bridge formation also activates skeletal myosin light chain kinase (skMLCK), a phosphotransferase enzyme responsible for introducing a negative charge to a serine residue on the myosin regulatory light chain (RLC) unit (Rayment et al., 1993; Sweeney et al., 1994). In turn, skMLCK-catalyzed phosphorylation of the RLC may induce structure–function changes to myosin that modulate the ability of Ca2+ to regulate force (see Vandenboom, 2016). As an example, skMLCK-catalyzed phosphorylation of the RLC increases the Ca2+ sensitivity of force in permeabilized skeletal muscle fibers (Metzger et al., 1989; Persechini et al., 1985; Sweeney and Kushmerick, 1985). Importantly, although steady-state force is enhanced, the force:ATPase relationship is not (Sweeney and Stull, 1990), suggesting that phosphorylation of the RLC promotes the ability of cycling cross-bridges to attain a force-generating state independent of any change to myosin ATPase rate (e.g. Morgan et al., 1976).

Myosin phosphorylation-mediated alterations to cross-bridge structure that increase the Ca2+ sensitivity of the contractile apparatus may also modulate force output of intact skeletal muscle, either in situ or in vitro (reviewed by MacIntosh, 2003). For example, strong associations between RLC phosphorylation and contraction-induced increases in muscle force and/or work, i.e. muscle potentiation, have been presented in a number of rodent skeletal muscle models (reviewed by Vandenboom et al., 2013). Moreover, fast twitch muscles from skMLCK knockout mice that do not display contraction-induced elevations in RLC phosphorylation display greatly reduced levels of potentiation compared with their wild-type (WT) counterparts (Bowslaugh et al., 2016; Gittings et al., 2011, 2017; Zhi et al., 2005). Although redundant mechanisms for potentiation may exist, an important question is whether skMLCK-catalyzed phosphorylation of the RLC significantly alters the contractile economy (i.e. mechanical output:metabolic input) of potentiated skeletal muscle. Crow and Kushmerick (1982) were the first to examine this relationship, and reported a positive correlation between RLC phosphorylation and reduced energy cost for isometric tetani of mouse extensor digitorum longus (EDL) muscle (in vitro, 22°C). However, subsequent work on rat gastrocnemius muscle (in situ, 35°C) showed that RLC phosphorylation had no consistent relationship with energy consumption during isometric tetani (Barsotti and Butler, 1984) or cross-bridge cycling rate during unloaded shortening (Butler et al., 1983). More recently, Abbate et al. (2001) found that potentiation had a greater effect on high-energy phosphate turnover than on work, thus reducing the economy of rat gastrocnemius muscle (in situ, 35°C). An important point was that it was the extra work caused by potentiation, and not the energetic cost of RLC phosphorylation per se, that decreased economy. Thus, no concordance in the literature exists regarding the influence of RLC phosphorylation-mediated potentiation on the contractile economy of rodent fast twitch muscle.

A limitation of previous studies examining the interaction of myosin RLC phosphorylation-mediated potentiation and contractile economy of rodent fast twitch muscle is the lack of a robust negative control for skMLCK activity and RLC phosphorylation. This is because even brief stimulation of most rodent fast twitch skeletal muscles results in fractional activation of skMLCK, making it difficult to parse out the influence of RLC phosphorylation from other processes. The purpose of this study was to compare the contractile economy of WT muscles and skMLCK−/− muscles lacking the enzymatic ability to phosphorylate myosin RLC. To do this, we used a paradigm that replicated the onset of low-frequency activity, a situation where the energetic cost of myosin phosphorylation, and its influence on economy, should be greatest. The use of different genotypes allowed us to parse out the influence of RLC phosphorylation on potentiation and high-energy phosphate consumption. We hypothesized that WT muscles with myosin phosphorylation would potentiate more than skMLCK−/− muscles, but that energy consumption would be increased in proportion to this increase in mechanical output, thus maintaining contractile economy compared with skMLCK−/− muscles without myosin phosphorylation. This may be a critical teleological question as a large decrease in economy may limit the usefulness of the potentiated state for working skeletal muscle (Brown and Loeb, 1998; Tsianos and Loeb, 2014).

Study design

This study was approved by the Brock University Animal Care Committee and conforms to the standards of the Canadian Council for Animal Care. Adult skMLCK−/− mice (Mus musculus Linnaeus 1758, C57BL/6 background, 19.4±0.7 g, 8–12 weeks old, n=20) were obtained from our breeding colony at Brock University. WT C57BL/6 mice (Mus musculus, 19.4±0.3 g, 8–12 weeks old, n=20) were purchased from Charles River Laboratories (Saint Constant, QC, Canada). Identical experiments were conducted on the two genotypes to determine contractile economy with and without the ability to phosphorylate myosin RLC, and to investigate potential differences in total energy turnover (Etotal), as well as energy turnover for contractile (cross-bridge energy turnover, EXB) and non-contractile (activation energy, EA: i.e. SERCA and Na+/K+-ATPase) processes. A major limitation of rapid freezing followed by chemical analysis of extracts of the frozen muscle is that each muscle provides only a single measurement, at a single time point. Because it was impossible to adopt a repeated measures design, half of the muscles were used for determination of myosin RLC phosphorylation and the other half were used for measurement of high-energy phosphate consumption (HEPC) (i.e. multiple independent measurements).

Contractile experiments and analysis

On the day of an experiment, mice were anesthetized with an intraperitoneal injection of sodium pentobarbital (60 mg kg−1 body mass). Bilateral EDL muscles were surgically isolated from each hindlimb of WT and skMLCK−/− mice and immediately (<5 s) suspended in a vertical jacketed organ bath (Radnoti LLC, Monrovia, CA, USA) containing continuously oxygenated Tyrode's solution (Lännergren et al., 2000) maintained at 25.0±0.1°C using an Isotemp 4100 R20 circulator (Fisher Scientific Company, Ottawa, ON, Canada). Contractile experiments were conducted using a 1200a In Vitro Muscle Testing System and field stimulation was applied using flanking platinum electrodes driven by a Model 701B biphase stimulator (Aurora Scientific, Inc., Aurora, ON, Canada) with voltage set to 1.25 times the threshold required to activate all fibers and elicit maximal isometric twitch force. Muscle length was measured using digital vernier calipers with LCD display. All contractile data were sampled and collected at 2000 Hz from a 305B servomotor acquired through a 604C analog-to-digital interface, and controlled by a dual-mode lever system (Aurora Scientific, Inc.). Data acquisition and basic analysis were performed using Aurora Scientific, Inc. 600a software (v1.60) and further examined using Excel (Microsoft Canada Co., Mississauga, ON, Canada). After mounting, muscles were allowed a 30 min equilibration period, after which optimal length (Lo) for maximal isometric twitch force was determined. The muscle length was adjusted to 0.9 Lo to maximize the magnitude of isometric twitch potentiation (Rassier and MacIntosh, 2000), and following an 80 min period of rest, the muscle was repetitively stimulated at 10 Hz for 15 s for a total of 150 individual twitches (Fig. 1). In a subset of muscles, stimulation was performed in the presence of BTS (N-benzyl-p-toluene sulfonamide), a myosin ATPase inhibitor that does not interfere with intracellular Ca2+ handling (Pinniger et al., 2005; Shaw et al., 2003). In either case, relative twitch force was quantified as the ratio of final twitch force to initial twitch force, producing a value that reflects an increase or decrease in twitch force relative to baseline measures. Each experiment was terminated when the EDL was rapidly frozen (<2 s) using tongs pre-cooled in liquid nitrogen. Muscles were then stored at −80°C for future determination of either RLC phosphorylation or metabolite concentrations (with or without BTS).

Fig. 1.

Representative forcetraces during repetitive low-frequency stimulation (in vitro, 25°C). (A) Mice lacking skeletal muscle myosin light chain kinase (skMLCK−/−). (B) Wild-type (WT) muscles. Muscles were stimulated at 10 Hz for 15 s to induce staircase potentiation at 0.90 optimal length (Lo). In these examples, the final twitch was 1.04 and 1.60 of baseline values in the skMLCK−/− and WT muscle, respectively. Note also that the depression in isometric twitch force (i.e. negative staircase) evident over the first several pulses is considerably reduced in amplitude and duration in the WT muscle compared with the skMLCK−/− muscle.

Fig. 1.

Representative forcetraces during repetitive low-frequency stimulation (in vitro, 25°C). (A) Mice lacking skeletal muscle myosin light chain kinase (skMLCK−/−). (B) Wild-type (WT) muscles. Muscles were stimulated at 10 Hz for 15 s to induce staircase potentiation at 0.90 optimal length (Lo). In these examples, the final twitch was 1.04 and 1.60 of baseline values in the skMLCK−/− and WT muscle, respectively. Note also that the depression in isometric twitch force (i.e. negative staircase) evident over the first several pulses is considerably reduced in amplitude and duration in the WT muscle compared with the skMLCK−/− muscle.

Analysis of twitch characteristics

A rapid increase in force occurred over the initial 0.5 s of repetitive stimulation due to partial fusion of twitches and thus final twitch force was normalized to the twitch occurring at 0.5 s (i.e. the 6th twitch in the series) (Ryder et al., 2007). The initial, 25th, 50th and final twitch of repetitive stimulation were analyzed for peak force (Pt), time to peak tension (TPT), half-relaxation time (RT1/2), rate of force development (+dP/dt) and rate of force relaxation (−dP/dt). As all twitches were isometric (i.e. fixed length), we could not quantify work per se; instead, the force–time integral (FTI) was calculated and the area under the force–time curve was used as a surrogate for isometric work (de Haan et al., 1986).

Myosin RLC phosphorylation

In parallel experiments, a subset of muscles were removed from −80°C storage and myosin RLC phosphate content was quantified by urea–glycerol PAGE followed by immunoblotting. These methods were based on Zhi et al. (2005) and have been described in detail previously (see Gittings et al., 2016; Bowslaugh et al., 2016). The primary myosin RLC antibody used in this work was graciously donated by Dr Jim Stull and his lab at the University of Texas Southwestern Medical Center.

HEPC

Muscles were removed from −80°C storage and lyophilized (FreeZone 4.5, Labconco, Kansas City, MO, USA) for a minimum of 6 h to remove all water. After the silk sutures were removed, the muscles were dissected free of connective tissue and blood, then powdered and aliquoted into 1.5 ml microcentrifuge tubes. Metabolites were extracted from aliquots of lyophilized tissue using 0.5 mol l−1 perchloric acid (HClO4), and neutralized with 2.3 mol l−1 potassium carbonate (KHCO3). The concentrations of metabolites were analyzed in triplicate in opaque black polypropylene 96-well microplates (Greiner Bio-One America Inc., Monroe, NC, USA) using fluorometric techniques as previously described (Bergmeyer et al., 1983; Harris et al., 1974). Three specific assays were used to identify the concentration of the metabolites of interest, including ATP, phosphocreatine (PCr), creatine (Cr) and lactate. To adjust for variability in solid non-muscle constituents, all raw metabolite values were normalized to total Cr content (Zhang et al., 2006). HEPC, a measure of ATP turnover determined through anaerobic ATP production, was calculated for each group from changes of metabolite concentration using Eqn 1:
(1)
The change (Δ) signifies the differences in metabolite concentration between the experimental muscle and its contralateral resting control. The calculation for HEPC thus uses an equation that assumes that glycogen is the sole source for lactate production during contraction and any decrease in ATP within the muscle is associated with a stoichiometric increase in inosine monophosphate (IMP) (Zhang et al., 2006). Energy balance experiments show that the contribution to overall ATP production by oxidative metabolism is no more than 9% during brief contractile activity (<30 s) (Westra et al., 1988) and thus this contribution to energy turnover was ignored. Although we used a 15 s stimulation protocol in order to increase PCr breakdown and reduce uncertainty in its measurement, we cannot exclude entirely the possibility that errors in its measurement are propagated through calculations of HEPC.

Calculation of contractile economy

Contractile economy of both genotypes was calculated as the ratio of total FTI obtained during low-frequency stimulation to HEPC, as measured by NADH-linked fluorometric assays. Parallel experiments were conducted in which BTS was used to pharmacologically inhibit myosin II ATPase activity (Shaw et al., 2003). This procedure allowed us to partition Etotal into EXB and EA components after directly measured HEPC values were corrected for residual force production during the BTS incubation using Eqn 2 (Zhang et al., 2006):
(2)
Once this was done, new corrected HEPC values for cross-bridge and EA components were calculated. Contractile economy was defined as the ratio of mechanical output to metabolic input, i.e. FTI to HEPC, and expressed as millinewton-seconds per micromole of high-energy phosphate (mN s μmol−1 P) (de Haan et al., 1986). BTS was solubilized in dimethyl sulfoxide (DMSO) to create a 25 mmol l−1 stock solution. A small volume of BTS stock solution was added to the organ bath following preliminary procedures to yield a final concentration of 25 μmol l−1. Muscles were incubated in BTS for 80 min before the contractile experiment began.

Statistics

After checking for normal distribution and homogeneity of variances, a two-way repeated measures ANOVA was conducted to evaluate the effect of skMLCK expression (i.e. WT or skMLCK−/−) and repetitive stimulation (i.e. resting or stimulated) on isometric twitch force and work in the control condition. A two-way ANOVA was conducted to compare the effect of skMLCK expression and repetitive stimulation on RLC phosphorylation. Two-tailed, unpaired t-tests were conducted to examine the effect of skMLCK expression on total FTI, Etotal, EXB and EA and contractile economy. All statistical analyses were conducted using GraphPad Prism 6 (La Jolla, CA, USA) with significance set at P<0.05. All data are means±s.e.m. unless otherwise stated. Sample sizes (n) used to achieve statistical power were based upon previous work from our lab and others for both mechanical (Gittings et al., 2011, 2012, 2016; Abbate et al., 2001) and biochemical (Zhi et al., 2005;Zhang et al., 2006) analyses.

Effect of repetitive stimulation on isometric twitch force in WT and skMLCK−/− muscles

Representative twitch records obtained at select time points are shown in Fig. 2A,B and relative changes to WT and skMLCK−/− muscles during repetitive stimulation are summarized in Fig. 2C. Initial absolute isometric twitch force was similar between WT and skMLCK−/− muscles (45.8±3.5 and 40.2±2.8 mN, respectively, P=0.340); however, with repetitive stimulation, WT absolute isometric twitch force increased significantly more than skMLCK−/− values (56.3±2.4 and 42.2±3.0 mN, respectively, P<0.01). In addition to changes to peak force, the TPT and RT1/2 of twitches of each genotype decreased from the first to the final twitch in the series (Table 1). When genotypes are compared, WT TPT was significantly higher at each time point compared with that of skMLCK−/− muscles. By contrast, the RT1/2 of WT muscles was greater than that of skMLCK−/− muscles at the first and final twitches only. With respect to twitch kinetics, the initial +dP/dt and −dP/dt of each genotype was similar and each increased with repetitive stimulation; however, +dP/dt was increased more in the WT than in the skMLCK−/− genotype while the change to −dP/dt was similar for the two genotypes. WT muscles performed significantly more isometric work (i.e. FTI) than skMLCK−/− muscles (17.01±1.7 and 10.21±0.8 mN s, respectively, P<0.01) during the protocol.

Fig. 2.

Contractile data from WT and skMLCK−/− muscles subjected to repetitive low-frequency stimulation (in vitro, 25°C). (A,B) Representative isometric twitch force records of extensor digitorum longus (EDL) muscles from skMLCK−/− (A) and WT (B) mice at selected times during repetitive low-frequency stimulation. (C,D) Summary of relative change in mean isometric twitch force (C; means±s.e.m., n=8, **P<0.01 for one-tailed Mann–Whitney rank sum test), and force–time integral (FTI; D; means±s.e.m., n=8, **P<0.01 for one-tailed unpaired t-test) during repetitive stimulation of skMLCK−/− and WT muscles.

Fig. 2.

Contractile data from WT and skMLCK−/− muscles subjected to repetitive low-frequency stimulation (in vitro, 25°C). (A,B) Representative isometric twitch force records of extensor digitorum longus (EDL) muscles from skMLCK−/− (A) and WT (B) mice at selected times during repetitive low-frequency stimulation. (C,D) Summary of relative change in mean isometric twitch force (C; means±s.e.m., n=8, **P<0.01 for one-tailed Mann–Whitney rank sum test), and force–time integral (FTI; D; means±s.e.m., n=8, **P<0.01 for one-tailed unpaired t-test) during repetitive stimulation of skMLCK−/− and WT muscles.

Table 1.

Isometric twitch characteristics of WT and skMLCK−/− muscles during repetitive stimulation

Isometric twitch characteristics of WT and skMLCK−/− muscles during repetitive stimulation
Isometric twitch characteristics of WT and skMLCK−/− muscles during repetitive stimulation

Effects of BTS during repetitive stimulation in WT and skMLCK−/− muscles

An 80 min incubation in 25 µmol l−1 BTS abolished isometric twitch force of WT and skMLCK−/− muscles similarly (by 94±1% and 95±1%, respectively). Changes to work (as FTI) of WT and skMLCK−/− muscles were also similar (by 88±2% and 90±1%, respectively) (Fig. 3).

Fig. 3.

Summary of effect of 25 µmol l−1 BTS on isometric twitch force of WT and skMLCK−/− muscles (in vitro, 25°C). (A) Time course of relative force inhibition for WT and skMLCK−/− muscles over an 80 min incubation period. Values are means±s.e.m. for each group (n=7–8). (B) Representative traces obtained 5 min apart during N-benzyl-p-toluene sulfonamide (BTS) incubation depicting the effect on twitch amplitude and time course (only WT shown). (C) Representative WT record during repetitive low-frequency stimulation in the absence of BTS; staircase potentiation was ∼20% in this example. (D) Representative WT record during repetitive low-frequency stimulation in the presence of BTS; note that despite a force depression of ∼95%, staircase potentiation was still observed.

Fig. 3.

Summary of effect of 25 µmol l−1 BTS on isometric twitch force of WT and skMLCK−/− muscles (in vitro, 25°C). (A) Time course of relative force inhibition for WT and skMLCK−/− muscles over an 80 min incubation period. Values are means±s.e.m. for each group (n=7–8). (B) Representative traces obtained 5 min apart during N-benzyl-p-toluene sulfonamide (BTS) incubation depicting the effect on twitch amplitude and time course (only WT shown). (C) Representative WT record during repetitive low-frequency stimulation in the absence of BTS; staircase potentiation was ∼20% in this example. (D) Representative WT record during repetitive low-frequency stimulation in the presence of BTS; note that despite a force depression of ∼95%, staircase potentiation was still observed.

Myosin RLC phosphorylation in WT and skMLCK−/− muscles

At rest, WT muscles had 11-fold more RLC phosphorylation than skMLCK−/− muscles (i.e. 0.22±0.02 and 0.02±0.01 mol phosphate mol−1 RLC, respectively) (Fig. 4). In WT muscles frozen immediately after repetitive stimulation, RLC phosphorylation was increased ∼3.5-fold but was not changed in skMLCK−/− muscles (to 0.76±0.02 and 0.02±0.01 mol phosphate mol−1 RLC, respectively). Data for muscles from control and matched BTS conditions were similar and thus were pooled for this analysis.

Fig. 4.

Effect of skMLCK expression on myosin regulatory light chain (RLC) phosphorylation of mouse EDL muscle (in vitro, 25°C). (A) Representative urea–glycerol PAGE blot of non-phosphorylated and monophosphorylated myosin RLC of isolated EDL muscles from WT and skMLCK−/− mice. One muscle was frozen after an 80 min equilibration period (Resting) while the contralateral muscle from the same animal was frozen after 15 s of repetitive low-frequency stimulation (Stimulated). Fractional myosin phosphorylation was calculated as the ratio of the phosphorylated band to total light chain content (RLC phosphorylation=[RLC-P]/[RLC]+[RLC-P]). Note that the top bands in each lane are unidentified in each genotype. (B) Summary of effect of genotype (skMLCK−/− versus WT) and stimulation (resting versus stimulated) on RLC phosphorylation of mouse EDL muscle. RLC phosphorylation is expressed as mol phosphate mol−1 RLC. Data from BTS-untreated and -treated muscles were pooled for this analysis. Values are means±s.e.m.; n=6, ****P<0.0001 for two-way ANOVA with post hoc analyses conducted using Tukey's HSD.

Fig. 4.

Effect of skMLCK expression on myosin regulatory light chain (RLC) phosphorylation of mouse EDL muscle (in vitro, 25°C). (A) Representative urea–glycerol PAGE blot of non-phosphorylated and monophosphorylated myosin RLC of isolated EDL muscles from WT and skMLCK−/− mice. One muscle was frozen after an 80 min equilibration period (Resting) while the contralateral muscle from the same animal was frozen after 15 s of repetitive low-frequency stimulation (Stimulated). Fractional myosin phosphorylation was calculated as the ratio of the phosphorylated band to total light chain content (RLC phosphorylation=[RLC-P]/[RLC]+[RLC-P]). Note that the top bands in each lane are unidentified in each genotype. (B) Summary of effect of genotype (skMLCK−/− versus WT) and stimulation (resting versus stimulated) on RLC phosphorylation of mouse EDL muscle. RLC phosphorylation is expressed as mol phosphate mol−1 RLC. Data from BTS-untreated and -treated muscles were pooled for this analysis. Values are means±s.e.m.; n=6, ****P<0.0001 for two-way ANOVA with post hoc analyses conducted using Tukey's HSD.

Effects of repetitive stimulation on metabolites in WT and skMLCK−/− muscles without BTS

Resting metabolite concentrations did not differ between genotypes or incubation conditions in unstimulated muscles (Table 2). In addition, no significant differences existed between genotypes or incubation conditions in concentrations of PCr, Cr or ATP following repetitive stimulation; however, WT control muscles had significantly greater lactate concentrations after stimulation than WT BTS and skMLCK−/− BTS muscles (Table 2).

Table 2.

Effect of BTS on metabolites of unstimulated and stimulated WT and skMLCK−/− muscles

Effect of BTS on metabolites of unstimulated and stimulated WT and skMLCK−/− muscles
Effect of BTS on metabolites of unstimulated and stimulated WT and skMLCK−/− muscles

No differences in PCr, Cr, ATP or lactate were detected between genotypes, with or without BTS, in unstimulated muscles (Table 2). Similar to the basal condition, no significant differences in PCr, Cr or ATP were found between genotypes with or without BTS following stimulation; in both genotypes, control muscles did have significantly greater lactate concentrations than BTS-treated muscles, however (Table 2).

Energy consumption and contractile economy of WT and skMLCK−/− muscles

Total energy consumption, partitioned into cross-bridge and activation energy consumption, is shown for both genotypes in Fig. 5. No effect of genotype was found as the Etotal of WT and skMLCK−/− muscles was similar (Table 3). When Etotal was partitioned by BTS treatment, EA was also similar between genotypes, but EXB was greater for skMLCK−/− than for WT muscles. Because work was increased but Etotal was similar, the contractile economy of WT muscles was significantly greater than for skMLCK−/− muscles (Fig. 6).

Fig. 5.

High-energy phosphate consumption (HEPC) of WT and skMLCK−/− muscles during repetitive low-frequency stimulation (in vitro, 25°C). Bars show total energy (ETOTAL), cross-bridge (EXB) and activation (EA) energy components for each genotype when corrected for residual force production in the presence of BTS. All values are means±s.e.m. (n=7, contralateral pairs). *P<0.05 for two-tailed unpaired t-test. See Materials and methods for details.

Fig. 5.

High-energy phosphate consumption (HEPC) of WT and skMLCK−/− muscles during repetitive low-frequency stimulation (in vitro, 25°C). Bars show total energy (ETOTAL), cross-bridge (EXB) and activation (EA) energy components for each genotype when corrected for residual force production in the presence of BTS. All values are means±s.e.m. (n=7, contralateral pairs). *P<0.05 for two-tailed unpaired t-test. See Materials and methods for details.

Fig. 6.

Contractile economy of WT and skMLCK−/− muscles during repetitive low-frequency stimulation (in vitro, 25°C). Values are expressed as means±s.e.m. (n=7). *P<0.05 for two-tailed unpaired t-test with Welch's correction.

Fig. 6.

Contractile economy of WT and skMLCK−/− muscles during repetitive low-frequency stimulation (in vitro, 25°C). Values are expressed as means±s.e.m. (n=7). *P<0.05 for two-tailed unpaired t-test with Welch's correction.

Table 3.

Total high-energy phosphate consumption of WT and skMLCK−/− muscles during repetitive low-frequency stimulation

Total high-energy phosphate consumption of WT and skMLCK−/− muscles during repetitive low-frequency stimulation
Total high-energy phosphate consumption of WT and skMLCK−/− muscles during repetitive low-frequency stimulation

The purpose of this study was to investigate the effects of RLC phosphorylation-mediated force potentiation on the energetic cost of contraction during low-frequency stimulation of fast twitch skeletal muscle isolated from WT and skMLCK−/− mice. A novel aspect of this work is that it is the first to directly examine the effects of RLC phosphorylation on muscle energetics and contractile economy by using BTS to partition energy consumption into contractile and non-contractile components. Our results in indicate that skMLCK-catalyzed increases in RLC phosphorylation augment peak isometric twitch force and work without increasing energy consumption, thus increasing low-frequency contractile economy of mouse fast twitch muscle in vitro.

Isometric force and work during staircase potentiation

The progressive increase in isometric force during repetitive stimulation displayed by rodent fast twitch skeletal muscle, referred to as staircase potentiation, is well documented (Close and Hoh, 1968; Krarup, 1981a; MacIntosh et al., 1993; MacIntosh and Kupsh, 1987; Smith et al., 2014). Our findings of increased force and work by WT but not skMLCK−/− muscles is consistent with the idea that RLC phosphorylation causes staircase potentiation. Evidence that RLC phosphorylation is the main, but not only, mechanism for this response is that muscles from skMLCK−/− mice display much reduced staircase potentiation compared with muscles from WT mice (Zhi et al., 2005). In addition, rat disuse or denervation models show that muscles with reduced RLC phosphorylation also display reduced staircase potentiation (MacIntosh et al., 2008; Rassier et al., 1999). The fact that the initial decrease in force observed during stimulation (i.e. negative staircase) was greater in skMLCK−/− than in WT muscles also accords with the idea that, without RLC phosphorylation, low-frequency forces are compromised. By contrast, changes to twitch TPT or RT1/2 noted during staircase may be unrelated to RLC phosphorylation (Smith et al., 2014).

Myosin RLC phosphorylation

The levels of RLC phosphorylation we observed following low-frequency stimulation producing staircase potentiation are similar to those reported previously using similar stimulus protocols (MacIntosh et al., 1993) or similar muscle models (Moore et al., 1990). Consistent with previous work comparing the two genotypes, we found that RLC phosphorylation was increased by stimulation in WT but not in skMLCK−/− muscles (Bowslaugh et al., 2016; Gittings et al., 2011, 2016; Zhi et al., 2005). These data confirm the critical influence of skMLCK for stimulation-induced RLC phosphorylation levels. Interestingly, because skMLCK activity is Ca2+ dependent (Stull et al., 2011), the fact that stimulated values for WT RLC phosphorylation were similar with or without BTS is consistent with the idea that this drug does not influence Ca2+ handling (Pinniger et al., 2005).

Energy consumption and economy of WT versus skMLCK−/− muscle

The ATP cost of skMLCK-catalyzed RLC phosphorylation has been estimated to be ∼5% of total energy used during a 5 s tetanus (at 0°C) (Homsher, 1987). After correcting for residual force development during BTS incubation, no differences in Etotal or EA were found between genotypes in the present study (Fig. 5). This suggests that the energetic cost of skMLCK-catalyzed RLC phosphorylation was negligible and below the detection limit of our assays. This may be due to the very slow kinetics of the PP1 phosphatase enzyme which greatly reduces the skMLCK activity required to maintain RLC phosphorylation levels during stimulation (Stull et al., 2011).

An important assumption in the present study is that the genotype-dependent difference in contractile economy we observed was due to the presence of skMLCK-catalyzed phosphorylation of the RLC in WT but not skMLCK−/− muscles. Although the precise mechanism for this putative effect of RLC phosphorylation cannot be identified by our data, work on single rabbit myosin molecules indicates that RLC phosphorylation does not directly influence myosin ATPase rate (Morgan et al., 1976). In accord with this, genotype-dependent differences in potentiation due to phosphorylation-mediated increases in the number of attached cross-bridges should produce a commensurate increase in EXB (see fig. 2 of Sweeney and Stull, 1990). Thus, a working hypothesis for our results is that, in addition to a myosin phosphorylation-mediated increase in the portion of cycling cross-bridges able to participate in isometric force generation, there may have also been an increase in the force generated per attached cross-bridge. Results from studies on single myosin motor molecules from rabbit psoas muscle may provide clues as to how this unitary property may have been altered by RLC phosphorylation. For example, Greenberg et al. (2009) showed that phosphorylated cross-bridges detach more slowly than unphosphorylated cross-bridges, thus increasing cross-bridge duty cycle (cf. Duggal et al., 2014). Alternatively, Greenberg et al. (2010) showed that phosphorylation may increase the rigidity of the light chain-binding domain, thus increasing cross-bridge stiffness. In principle, either of these molecular outcomes could increase the thermodynamic efficiency of the cross-bridge force-generating cycle by increasing the life time of the force-generating state or by increasing the force exerted per force-generating cross-bridge, respectively. However, a phosphorylation-mediated slowing of cross-bridge detachment rate is incompatible with results from intact muscle showing no change in unloaded shortening velocity, a mechanical surrogate for cross-bridge detachment kinetics, between WT and skMLCK−/− or between phosphorylated and unphosphorylated muscles (Barsotti and Butler, 1984;Butler et al., 1983; Gittings et al., 2011; Palmer and Moore, 1989). Thus, we favor the idea that the disproportionate increase in cross-bridge force and work to EXB we found was due to stiffness-mediated increases in force per attached cross-bridge.

Relationship with pertinent work

Our data are consistent with results from mouse EDL muscle showing stimulation duration-dependent reductions in the average rate of energy consumption (Barsotti and Butler, 1984). They also accord with the results of Crow and Kushmerick (1982), who reported a correlation between RLC phosphorylation and reduced energy cost for isometric tetani in mouse EDL muscle. However, our findings are in stark contrast to those from rat skeletal muscle in situ showing that RLC phosphorylation increased the energetic cost of high-frequency, isovelocity contractions to a higher extent than it increased work output, thus decreasing economy (Abbate et al., 2001). Although isometric twitch force potentiation magnitude increases with temperature in the range 25–35°C (Krarup, 1981b; Moore et al., 1990), cross-bridge efficiency of mammalian muscle has been shown to be constant between 25 and 30°C (Barclay et al., 2010). Thus, excluding the possibility of a steep change at temperatures just above 30°C, we suggest that differences between our study and that of Abbate et al. (2001) are due to stimulation frequency and/or contraction-type dependencies for the influence of RLC phosphorylation on force (e.g. Gittings et al., 2012). For example, low-frequency stimulation resulting in low levels of thin filament activation may result in lower work performed per cycling cross-bridge compared with the effects of stimulation at high frequencies (Lewis and Barclay, 2014). Thus, our staircase protocol may have preferentially enhanced contractile economy by amplifying the structure–function influence of RLC phosphorylation on cross-bridge stiffness and/or duty cycle compared with the tetanic protocol used by Abbate et al. (2001). Alternatively, the increase in cross-bridge efficiency displayed by rodent skeletal muscle at moderate compared with slow or isometric speeds of shortening (see fig. 11 of Barclay, 2015) may have negatively biased the ability of the RLC phosphorylation mechanism to further enhance efficiency. More work is needed to distinguish between these possibilities.

BTS experiments

Consistent with previous work on rat skeletal muscle, our BTS experiments revealed a relatively low contribution of EXB to Etotal in either genotype (∼40%) (Siegman et al., 1994; Zhang et al., 2006). Although BTS is a selective inhibitor of myosin II ATPase activity (Bruton et al., 2006; Pinniger et al., 2005), this may be accounted for on the basis of temperature as SERCA ATPase activity has been shown to comprise ∼50% of total ATP consumption at 20°C in rat muscle (Rall, 1982; Stienen et al., 1995). Alternatively, because mouse EDL muscle consists predominately of IIX and IIB fibers (Smith et al., 2013), this profile may have contributed to a relatively high EA versus EXB compared with other muscle models with more heterogeneous fiber-type composition. Despite our finding that EA accounts for a high portion of Etotal during repetitive low-frequency stimulation of mouse EDL (Barclay and Loiselle, 2007), the EXB is still large enough for a myosin phosphorylation-mediated influence on muscle energetics.

Summary

This study provides evidence that RLC phosphorylation-mediated force potentiation improves the contractile economy of fast twitch mammalian muscle during low-frequency, repetitive stimulation of mouse EDL muscle in vitro. Use of skMLCK−/− muscles provided a negative control allowing us to show that the energetic overhead for skMLCK-catalyzed phosphorylation of the RLC is insignificant compared with the overall energy cost, and that RLC phosphorylation enhances the thermodynamic efficiency of the cross-bridge force-generating cycle. Finally, isometric twitch potentiation is length dependent (Rassier and MacIntosh, 2002) and our experiments were performed at a length (0.90 Lo) intended to maximize this outcome. It thus remains to be shown whether experiments performed at longer lengths or those involving length excursions display the same changes to contractile economy as observed in the present experiments.

The authors would like to acknowledge the support of Dr J. Stull (Department of Physiology, University of Texas Southwestern Medical Center) in development of the skMLCK−/− mouse colony and with methods for quantification of myosin RLC phosphate content.

Author contributions

Conceptualization: J.B., W.G., R.V.; Methodology: J.B., W.G., R.V.; Formal analysis: J.B.; Investigation: J.B., W.G.; Resources: R.V.; Writing - original draft: J.B.; Writing - review & editing: R.V., J.B., W.G.; Supervision: R.V.; Project administration: R.V.; Funding acquisition: R.V.

Funding

This research was made possible by funding from the Natural Sciences and Engineering Council of Canada (RV Discovery Grant: 2014-05122).

Data availability

The datasets used and/or analyzed during the current study are available from the corresponding author on request.

Abbate
,
F.
,
Velden
,
J. V. D.
,
Stienen
,
G. J. M.
and
De Haan
,
A.
(
2001
).
Post-tetanic potentiation increases energy cost to a higher extent than work in rat fast skeletal muscle
.
J. Muscle Res. Cell M.
22
,
703
-
710
.
Barclay
,
C. J.
(
2015
).
Energetics of contraction
.
Compr Physiol.
5
,
961
-
995
.
Barclay
,
C. J.
and
Loiselle
,
D. S.
(
2007
).
Can activation account for 80% of skeletal muscle energy use during isometric contraction?
Am. J. Physiol. Cell Physiol.
292
,
C612
.
Barclay
,
C. J.
,
Woledge
,
R. C.
and
Curtin
,
N. A.
(
2010
).
Is the efficiency of mammalian (mouse) skeletal muscle temperature dependent?
J. Physiol.
588
,
3819
-
3831
.
Barsotti
,
R. J.
and
Butler
,
T. M.
(
1984
).
Chemical energy usage and myosin light chain phosphorylation in mammalian skeletal muscle
.
J. Muscle Res. Cell M.
5
,
45
-
64
.
Bergmeyer
,
H. U.
,
Bergmeyer
,
J.
and
Grassl
,
M.
(
1983
).
Methods of Enzymatic Analysis
.
Weinheim
,
Germany
:
Verlag Chemie
.
Bowslaugh
,
J.
,
Gittings
,
W.
and
Vandenboom
,
R.
(
2016
).
Myosin light chain phosphorylation is required for peak power output of mouse fast skeletal muscle in vitro
.
Pflugers Arch. – Eur. J. Physiol.
468
,
2007
-
2016
.
Brown
,
I. E.
and
Loeb
,
G. E.
(
1998
).
Post-activation potentiation: A clue for simplifying models of muscle dynamics
.
Am. Zool.
38
,
743
-
754
.
Bruton
,
J.
,
Pinniger
,
G. J.
,
Lännergren
,
J.
and
Westerblad
,
H.
(
2006
).
The effects of the myosin-II inhibitor N-benzyl-p-toluene sulphonamide on fatigue in mouse single intact toe muscle fibres
.
Acta Physiol.
186
,
59
-
66
.
Butler
,
T. M.
,
Seigman
,
M. J.
,
Mooers
,
S. V.
and
Barsotti
,
R. J.
(
1983
).
Myosin light chain phosphorylation does not modulate cross-bridge cycling rate in mouse skeletal muscle
.
Science
220
,
1167
-
1169
.
Close
,
R.
and
Hoh
,
J. F. Y.
(
1968
).
The after-effects of repetitive stimulation on the isometric twitch contraction of rat fast skeletal muscle
.
J. Physiol.
197
,
461
-
477
.
Crow
,
M. T.
and
Kushmerick
,
M. J.
(
1982
).
Myosin light chain phosphorylation is associated with a decrease in the energy cost for contraction in fast twitch mouse muscle
.
J. Biol. Chem.
257
,
2121
-
2124
.
de Haan
,
A.
,
de Jong
,
J.
,
van Doorn
,
J. E.
,
Huijing
,
P. A.
,
Woittiez
,
R. D.
and
Westra
,
H. G.
(
1986
).
Muscle economy of isometric contractions as a function of stimulation time and relative muscle length
.
Pflügers Arch. Eur. J. Physiol.
407
,
445
-
450
.
Duggal
,
D.
,
Nagwekar
,
J.
,
Rich
,
R.
,
Midde
,
K.
,
Fudala
,
R.
,
Gryczynski
,
I.
and
Borejdo
,
J.
(
2014
).
Phosphorylation of myosin regulatory light chain has minimal effect on kinetics and distribution of orientations of cross bridges of rabbit skeletal muscle
.
Am. J. Physiol. Regul. Integr. Comp. Physiol.
306
,
R222
-
R233
.
Gittings
,
W.
,
Huang
,
J.
,
Smith
,
I. C.
,
Quadrilatero
,
J.
and
Vandenboom
,
R.
(
2011
).
The effect of skeletal myosin light chain kinase gene ablation on the fatigability of mouse fast muscle
.
J. Muscle Res. Cell M.
31
,
337
-
348
.
Gittings
,
W.
,
Huang
,
J.
and
Vandenboom
,
R.
(
2012
).
Tetanic force potentiation of mouse fast muscle is shortening speed dependent
.
J. Muscle Res. Cell M.
33
,
359
-
368
.
Gittings
,
W.
,
Bunda
,
J.
,
Stull
,
J. T.
and
Vandenboom
,
R.
(
2016
).
Interaction of posttetanic potentiation and the catchlike property in mouse skeletal muscle
.
Muscle Nerve
54
,
308
-
316
.
Gittings
,
W.
,
Bunda
,
J.
and
Vandenboom
,
R.
(
2017
).
Shortening speed dependent force potentiation is attenuated but not eliminated in skeletal muscles without myosin phosphorylation
.
J. Muscle Res. Cell Motil.
38
,
157
-
162
.
Greenberg
,
M. J.
,
Mealy
,
T. R.
,
Watt
,
J. D.
,
Jones
,
M.
,
Szczesna-Cordary
,
D.
and
Moore
,
J. R.
(
2009
).
The molecular effects of skeletal muscle myosin regulatory light chain phosphorylation
.
Am. J. Physiol. Regul. Integr. Comp. Physiol.
297
,
R265
-
R274
.
Greenberg
,
M. J.
,
Mealy
,
T. R.
,
Jones
,
M.
,
Szczesna-Cordary
,
D.
and
Moore
,
J. R.
(
2010
).
The direct molecular effects of fatigue and myosin regulatory light chain phosphorylation on the actomyosin contractile apparatus
.
Am. J. Physiol. Regul. Integr. Comp. Physiol.
298
,
989
-
996
.
Harris
,
R. C.
,
Hultman
,
E.
and
Nordesjö
,
L.-O.
(
1974
).
Glycogen, glycolytic intermediates and high-energy phosphates determined in biopsy samples of musculus quadriceps femoris of man at rest. Methods and variance of values
.
Scan. J. Clin. Lab. Invest.
33
,
109
-
120
.
Homsher
,
E.
(
1987
).
Muscle enthalpy production and its relationship to actomyosin ATPase
.
Annu. Rev. Physiol.
49
,
673
-
690
.
Krarup
,
C.
(
1981a
).
Enhancement and diminution of mechanical tension evoked by staircase and by tetanus in rat muscle
.
J. Physiol.
311
,
355
-
372
.
Krarup
,
C.
(
1981b
).
Temperature dependence of enhancement and diminution of tension evoked by staircase and by tetanus in rat muscle
.
J. Physiol.
311
,
373
-
387
.
Lännergren
,
J.
,
Bruton
,
J. D.
and
Westerblad
,
H.
(
2000
).
Vacuole formation in fatigued skeletal muscle fibres from frog and mouse: effects of extracellular lactate
.
J. Physiol.
526
,
597
-
611
.
Lewis
,
D. B.
and
Barclay
,
C. J.
(
2014
).
Efficiency and cross-bridge work output of skeletal muscle is decreased at low levels of activation
.
Pflügers Arch. Eur. J. Physiol.
466
,
599
-
609
.
Lymn
,
R. W.
and
Taylor
,
E. W.
(
1971
).
Mechanism of adenosine triphosphate hydrolysis by actomyosin
.
Biochemistry
10
,
4617
-
4624
.
MacIntosh
,
B. R.
(
2003
).
Role of calcium sensitivity modulation in skeletal muscle performance
.
News Physiol. Sci.
18
,
222
-
225
.
MacIntosh
,
B. R.
and
Kupsh
,
C. C.
(
1987
).
Staircase, fatigue, and caffeine in skeletal muscle in situ
.
Muscle Nerve
10
,
717
-
722
.
MacIntosh
,
B. R.
,
Grange
,
R. W.
,
Cory
,
C. R.
and
Houston
,
M. E.
(
1993
).
Myosin light chain phosphorylation during staircase in fatigued skeletal muscle
.
Pflügers Arch.
425
,
9
-
15
.
MacIntosh
,
B. R.
,
Smith
,
M. J.
and
Rassier
,
D. E.
(
2008
).
Staircase but not posttetanic potentiation in rat muscle after spinal cord hemisection
.
Muscle Nerve
38
,
1455
-
1465
.
Metzger
,
J. M.
,
Greaser
,
M. L.
and
Moss
,
R. L.
(
1989
).
Variations in cross-bridge attachment rate and tension with phosphorylation of myosin in mammalian skinned skeletal muscle fibres
.
J. Gen. Physiol.
93
,
855
-
883
.
Moore
,
R. L.
,
Palmer
,
B. L.
,
Williams
,
S. L.
,
Tanabe
,
H.
,
Grange
,
R. W.
and
Houston
,
M. E.
(
1990
).
Effect of temperature on myosin phosphorylation in mouse skeletal muscle
.
Am. J. Physiol.
259
,
C432
-
C438
.
Morgan
,
M.
,
Perry
,
S. V.
and
Ottaway
,
J.
(
1976
).
Myosin light-chain phosphatase
.
Biochem. J.
157
,
687
-
697
.
Palmer
,
B. M.
and
Moore
,
R. L.
(
1989
).
Myosin light chain phosphorylation and tension potentiation in mouse skeletal muscle
.
Am. J. Physiol.
257
,
C1012
-
C1019
.
Persechini
,
A.
,
Stull
,
J. T.
and
Cooke
,
R.
(
1985
).
The effect of myosin phosphorylation on the contractile properties of skinned rabbit skeletal muscle fibers
.
J. Biol. Chem.
260
,
7951
-
7954
.
Pinniger
,
G. J.
,
Bruton
,
J. D.
,
Westerblad
,
H.
and
Ranatunga
,
K. W.
(
2005
).
Effects of a myosin-II inhibitor (N-benzyl-p-toluene sulphonamide, BTS) on contractile characteristics of intact fast-twitch mammalian muscle fibres
.
J. Muscle Res. Cell M.
26
,
135
-
141
.
Rall
,
J. A.
(
1982
).
Energetics of Ca2+ cycling during skeletal muscle contraction
.
Fed. Proc.
41
,
155
-
160
.
Rassier
,
D. E.
and
MacIntosh
,
B. R.
(
2000
).
Length dependence of staircase potentiation: interactions with caffeine and dantrolene sodium
.
Can. J. Physiol. Pharmacol.
78
,
350
-
357
.
Rassier
,
D. E.
and
MacIntosh
,
B. R.
(
2002
).
Sarcomere length-dependence of activity-dependent twitch potentiation in mouse skeletal muscle
.
BMC Physiol.
2
,
19
.
Rassier
,
D. E.
,
Tubman
,
L. A.
and
MacIntosh
,
B. R.
(
1999
).
Staircase in mammalian muscle without light chain phosphorylation
.
Braz. J. Med. Biol. Res.
32
,
121
-
1299
.
Rayment
,
I.
,
Rypniewski
,
W. R.
,
Schmidt-Bäse
,
K.
,
Smith
,
R.
,
Tomchick
,
D. R.
,
Benning
,
M. M.
and
Winkelmann
,
D. A.
(
1993
).
Three-dimensional structure of myosin subfragment- 1: A molecular motor
.
Science
261
,
50
-
58
.
Ryder
,
J. W.
,
Lau
,
K. S.
,
Kamm
,
K. E.
and
Stull
,
J. T.
(
2007
).
Enhanced skeletal muscle contraction with myosin light chain phosphorylation by a calmodulin-sensing kinase
.
Journal of Biol. Chem.
282
,
20447
-
20454
.
Shaw
,
M. A.
,
Ostap
,
E. M.
and
Goldman
,
Y. E.
(
2003
).
Mechanism of inhibition of skeletal muscle actomyosin by N-benzyl-p-toluenesulfonamide
.
Biochem.
42
,
6128
-
6135
.
Siegman
,
M. J.
,
Mooers
,
S. U.
,
Warren
,
T. B.
,
Warshaw
,
D. M.
,
Ikebe
,
M.
and
Butler
,
T. M.
(
1994
).
Comparison of the effects of 2,3-butanedione monoxime on force production, myosin light chain phosphorylation and chemical energy usage in intact and permeabilized smooth and skeletal muscles
.
J. Muscle Res. Cell M.
15
,
457
-
472
.
Smith
,
I. C.
,
Gittings
,
W.
,
Huang
,
J.
,
McMillan
,
E. M.
,
Quadrilatero
,
J.
,
Tupling
,
A. R.
and
Vandenboom
,
R.
(
2013
).
Potentiation in mouse lumbrical muscle without myosin light chain phosphorylation: Is resting calcium responsible?
J. Gen. Physiol.
141
,
297
-
308
.
Smith
,
I. C.
,
Vandenboom
,
R.
and
Tupling
,
A. R.
(
2014
).
Juxtaposition of the changes in intracellular calcium and force during staircase potentiation at 30 and 37°C
.
J. Gen. Physiol.
144
,
561
-
570
.
Stienen
,
G. J.
,
Zaremba
,
R.
and
Elzinga
,
G.
(
1995
).
ATP utilization for calcium uptake and force production in skinned muscle fibres of Xenopus laevis
.
J. Physiol.
482
,
109
.
Stull
,
J. T.
,
Kamm
,
K. E.
and
Vandenboom
,
R.
(
2011
).
Myosin light chain kinase and the role of myosin light chain phosphorylation in skeletal muscle
.
Arch. Biochem. Biophys.
510
,
120
-
128
.
Sweeney
,
H. L.
and
Kushmerick
,
M. J.
(
1985
).
Myosin phosphorylation in permeabilized rabbit psoas fibers
.
Am. J. Physiol.
249
,
C362
-
C365
.
Sweeney
,
H. L.
and
Stull
,
J. T.
(
1990
).
Alteration of cross-bridge kinetics by myosin light chain phosphorylation in rabbit skeletal muscle: Implications for regulation of actin-myosin interaction
.
Proc. Natl. Acad. Sci. USA
87
,
414
-
418
.
Sweeney
,
H. L.
,
Yang
,
Z.
,
Zhi
,
G.
,
Stull
,
J. T.
and
Trybus
,
K. M.
(
1994
).
Charge replacement near the phosphorylatable serine of the myosin regulatory light chain mimics aspects of phosphorylation
.
Proc. Natl. Acad. Sci. USA
91
,
1490
-
1494
.
Tsianos
,
G. A.
and
Loeb
,
G. E.
(
2014
).
Physiology and computational principles of muscle force generation
. In
Encyclopedia of Computational Neuroscience
.
New York
:
Springer
.
Vandenboom
,
R.
(
2016
).
Modulation of skeletal muscle contraction by myosin phosphorylation
.
Compr. Phyiol.
7
,
171
-
212
.
Vandenboom
,
R.
,
Gittings
,
W.
,
Smith
,
I. C.
,
Grange
,
R. W.
and
Stull
,
J. T.
(
2013
).
Myosin phosphorylation and force potentiation in skeletal muscle: Evidence from animal models
.
J. Muscle Res. Cell M.
34
,
317
-
332
.
Westra
,
H. G.
,
de Haan
,
A.
,
van Doorn
,
J. E.
and
de Haan
,
E. J.
(
1988
).
Anaerobic chemical changes and mechanical output during isometric tetani of rat skeletal muscle in situ
.
Pflügers Arch. Eur. J. Physiol.
412
,
121
-
127
.
Zhang
,
S.-J.
,
Andersson
,
D. C.
,
Sandström
,
M. E.
,
Westerblad
,
H.
and
Katz
,
A.
(
2006
).
Cross bridges account for only 20% of total ATP consumption during submaximal isometric contraction in mouse fast-twitch skeletal muscle
.
Am. J. Physiol. Cell Physiol.
291
,
C147
-
C154
.
Zhi
,
G.
,
Ryder
,
J. W.
,
Huang
,
J.
,
Ding
,
P.
,
Chen
,
Y.
,
Zhao
,
Y.
,
Kamm
,
K. E.
and
Stull
,
J. T.
(
2005
).
Myosin light chain kinase and myosin phosphorylation effect frequency-dependent potentiation of skeletal muscle contraction
.
Proc. Natl. Acad. Sci. USA
102
,
17519
-
17524
.

Competing interests

The authors declare no competing or financial interests.