Muscles consume metabolic energy for active movement, particularly when performing mechanical work or producing force. Less appreciated is the cost for activating muscle quickly, which adds considerably to the overall cost of cyclic force production. However, the cost magnitude relative to the cost of mechanical work, which features in many movements, is unknown. We therefore tested whether fast activation is costly compared with performing work or producing isometric force. We hypothesized that metabolic cost would increase with a proposed measure termed force rate (rate of increase in muscle force) in cyclic tasks, separate from mechanical work or average force level. We tested humans (N=9) producing cyclic knee extension torque against an isometric dynamometer (torque 22 N m, cyclic waveform frequencies 0.5–2.5 Hz), while also quantifying quadriceps muscle force and work against series elasticity (with ultrasonography), along with metabolic rate through respirometry. Net metabolic rate increased by more than four-fold (10.5 to 46.8 W) with waveform frequency. At high frequencies, the hypothesized force-rate cost accounted for nearly half (40%) of energy expenditure. This exceeded the cost for average force (17%) and was comparable to the cost for shortening work (43%). The force-rate cost is explained by additional active calcium transport necessary for producing forces at increasing waveform frequencies, owing to rate-limiting dynamics of force production. The force-rate cost could contribute substantially to the overall cost of movements that require cyclic muscle activation, such as locomotion.

Humans often expend energy to perform movement tasks where muscles are only intermittently or cyclically active. Two notable contributions to the metabolic cost of such tasks are for the mechanical work performed by muscle fascicles (Abbott et al., 1952; Margaria, 1968), and for the force exerted when fascicles are isometric (Crow and Kushmerick, 1982). Less appreciated is the cost for muscle activation during intermittent or cyclic conditions. This cost increases with briefer contractions (Hogan et al., 1998) or higher rates of force production (Doke and Kuo, 2007), and can exceed the cost for producing continuous isometric force (Chasiotis et al., 1987). However, there is also substantial metabolic energy expended in tasks that also entail work; for example, locomotion with the lower extremity (Margaria, 1976) or reaching with the upper extremity (Huang et al., 2012). But the relative cost for activating muscle versus performing work in such tasks remains unknown. It is therefore helpful to determine whether the cost for activating muscle is comparable to the cost for work and force production.

There is clearly an energetic cost for activating muscle under intermittent, isometric conditions. This has been demonstrated with square-wave, on–off activation patterns, which yield a considerably higher metabolic cost than continuous activation at similar overall contraction duration (Chasiotis et al., 1987; Spriet et al., 1988). Intermittent contraction also requires more metabolic energy when contraction frequency is higher (Bergström and Hultman, 1988; Hogan et al., 1998) or when duration of active force production is shorter (Beck et al., 2020), even when accounting for the cost of maintaining isometric force, which is roughly proportional to the force–time integral (Crow and Kushmerick, 1982). For example, metabolic cost increases with the frequency of muscle activation (Bergström and Hultman, 1988), leading to a doubling of energetic cost per unit force compared to continuous isometric force (Chasiotis et al., 1987). The underlying mechanism is thought not to be related to crossbridge work, as the activation cost has been attributed to active calcium transport (Hogan et al., 1998).

There is also a similar metabolic cost for cyclic, non-isometric movements. For example, cyclic leg swinging (Doke et al., 2005; Doke and Kuo, 2007) and ankle bouncing (Dean and Kuo, 2011) have costs increasing markedly with movement frequency (e.g. three- to four-fold for a 45% increase in frequency; Doke et al., 2005). The costs are not explained by the amount of mechanical work, and have also been associated with calcium transport (Doke and Kuo, 2007). For cyclic movements, the cost has been related to the rate of force production, or force rate, which may be considered an analogue of on–off intermittency, but for continuous, non-isometric movements. The cost of leg swinging, including a force-rate component, could account for one-third of the net metabolic cost of walking (Doke et al., 2005). The hypothesized force-rate cost may therefore be relevant to typical human activities.

However, it remains unclear to what extent the force-rate cost contributes to overall metabolic cost of movement. Everyday movements may have substantial costs for work and force production. This makes it difficult to separate the force-rate cost and quantify its contribution relative to other costs. The force-rate cost seems to be present when muscle–tendon complex work is controlled for (e.g. Chasiotis et al., 1987; Doke and Kuo, 2007). However, muscle fascicles can perform work on series elasticity, not observable from the overall muscle–tendon complex. Even in apparently isometric conditions, such work may contribute substantially to overall cost. After all, muscle fascicle work is biochemically and thermodynamically constrained to cost metabolic energy (Barclay, 2015). The hypothesized cost of force rate therefore needs to be quantified alongside the actual shortening work performed by fascicles.

The present study is intended to address that gap, by testing for a metabolic cost of cyclic force production while accounting for and estimating the costs for mechanical work and force. We simultaneously quantified work (against series elasticity), force and force rate during a cyclic force production task, along with the overall metabolic energy expenditure. The task was to cyclically produce voluntary knee extension torque against an isometric dynamometer, at an amplitude and range of frequencies comparable to everyday human movements. We hypothesized that (1) work and force–time integral fall short in accounting for metabolic cost at higher frequencies, and (2) the surplus cost is related to muscle force rate. We tested this by parametrically relating each contribution as a function of waveform frequency, and testing whether the cost of force rate is separable from and comparable to the costs of work and force–time integral.

We estimated quadriceps muscle force and mechanical work, and net metabolic rate as healthy adults produced cyclic, isometric knee extension torque (Fig. 1). Increasing frequency of cyclic torque was expected to yield greater fluctuations in muscle force and thereby greater rate of force development. We conducted a set of energetics trials to test for an associated metabolic cost, and a briefer set of mechanics trials to quantify muscle fascicle mechanical work (against series elasticity) and force via ultrasound imaging.

Nine healthy participants (6 male, 3 female, leg length=87.4±4.9 cm, body mass=70.6±13.1 kg, mean±s.d.) produced cyclic bilateral knee extension torque against an isometric dynamometer. The task was intended to require approximately 10% of maximal quadriceps muscle force and to impose contraction frequencies f within ranges encountered in daily life (mean 11 N m and 0.5–2.5 Hz, respectively). Real-time visual feedback was displayed to participants, showing measured and target waveforms (Fig. 1) for each frequency, all at a fixed torque range (0–22 N m) from combined legs. Torques were exerted at 20 deg knee flexion (relative to straight leg) and measured with a dynamometer chair (Biodex, Biodex Medical Systems, NY, USA), which was also used to strap and constrain movement of the trunk and legs. As a reference for maximum torque, participants performed maximum voluntary contractions (MVC) with knee extensor and knee flexor muscles at two knee angles (20 and 70 deg). In line with the knee extension torque–angle relationship (e.g. Kulig et al., 1984), extensor MVCs resulted in less torque at a 20 deg knee angle than at a 70 deg knee angle (97±20 versus 221±70 N m). For each condition, we quantified knee torque in terms of its time average , peak amplitude and peak rate (Fig. 1). Peak torque rate was estimated from the change in torque from 50% to 150% of average, divided by the time duration of that change. Prior to data collection, participants provided their written informed consent as approved by the University of Calgary Conjoint Health Research Ethics Board.

The energetics trials entailed cyclic torque production while respirometry data were recorded. Participants were first given practice with matching real-time torque targets for at least 5 min. They then performed the six experimental task conditions for 6 min each, with torque and respirometry (rates of oxygen consumption O2 and carbon dioxide production CO2; K5 system, Cosmed, Rome, Italy) averaged only for the final 3 min, to allow time for oxygen kinetics to stabilize (Xu and Rhodes, 1999). Metabolic rate was determined from O2 and the respiratory exchange ratio (RER, i.e. the ratio between CO2 and O2) (Brockway, 1987). We also recorded electromyography (EMG) to monitor muscle activity and co-contraction [vastus lateralis (VL), rectus femoris (RF) and biceps femoris (BF)]. VL and BF EMG were recorded bilaterally and were averaged between legs, RF EMG was recorded unilaterally. EMG amplitude was determined from the band-pass filtered recordings (cut-off frequency fc=30–500 Hz, Butterworth) and low-pass filtered to obtain the envelope (fc=5 Hz, Butterworth) (Hof, 1984). EMG amplitude was defined as the absolute of the Hilbert transformed signal (Huang et al., 1998), expressed relative to its maximum value during MVC and corrected for background noise using baseline subtraction (La Delfa et al., 2014). In addition to waveform targets, there was also a constant torque target of 11 N m to determine the cost of steady force production, with the same average torque as the waveform targets. To avoid potential ordering effects such as from fatigue, the order of task conditions was randomized, with at least 1 min of rest between conditions.

The separate mechanics trials employed ultrasound imaging to estimate muscle fascicle shortening and muscle work against series elasticity. These trials consisted of single-legged isometric torque production with a low-frequency triangle waveform torque target (f=0.02 Hz, range=0–16 N m) at three knee angles (θknee=15, 20 and 25 deg). The low-frequency mechanics trials were to simplify tracking of images across a range of torques encompassing those applied in the higher frequency energetics trials. They yielded a force–displacement relationship for fascicle shortening, which could then be used to estimate the shortening work performed during the energetics trials. Images were obtained for the VL using ultrasound (5 cm probe, 11 MHz basic-mode; Logiq E9, General Electric, Fairfield, USA), recorded at 30 Hz and approximately synchronized with torque recordings (via a sync pulse). Separate trials were performed for each of both legs and then averaged between legs. On one occasion, data from one of the legs were discarded owing to poor image quality.

VL fascicle lengths Lfas and pennation angles φ were estimated using a custom (MathWorks, Natick, MA, USA) ultrasound algorithm (van der Zee and Kuo, 2020 preprint), similar to a previously described algorithm (Ryan et al., 2019). Our algorithm first filtered the images to highlight line-like structures, namely, fascicles and aponeuroses, with thin and thick lines, respectively. The filtered images were then thresholded, and overall fascicle angle was determined using a line detection technique termed Hough transform (Zhou et al., 2015, 2012; Zhou and Zheng, 2008). Aponeuroses were identified using feature detection based on object length, orientation and position. Superficial and deep aponeurosis were identified as the longest object with near-horizontal orientation within the top and bottom image halves, respectively (for details, see van der Zee and Kuo, 2020 preprint). Identified aponeuroses were checked using visual inspection. φ was defined as the difference between fascicle angle and superficial aponeurosis angle. Lfas was estimated as the distance between deep and superficial aponeuroses, as measured along the fascicle orientation. The junction between fascicle and superficial aponeurosis was linearly extrapolated in the typical case when it was beyond the image frame (see Fig. 1). This procedure was performed for 17 equidistantly spaced locations across the image width to yield an average Lfas per image.

At fixed muscle–tendon complex (MTC) length, the change in series elastic element (SEE) length is equal and opposite to the change in the fascicle length along the MTC. Consequently, the SEE length changes ΔLSEE could be estimated from changes in Lfas and φ (both sampled as a function of Tknee at 450 uniformly spaced points) (Fukunaga et al., 2001). At fixed knee torque Tknee (and thus fixed FSEE), the change in MTC length is equal and opposite to the change in the fascicle length along the MTC. Consequently, the MTC length changes could be estimated from changes in Lfas and φ. The moment arm rMTC of the quadriceps MTC about the knee (at θknee=20 deg) was estimated from the average difference in MTC length between the 15 and the 25 deg trials, divided by the corresponding knee angle difference (i.e. Δθknee=10 deg), similar to typical tendon excursion methods (Ito et al., 2000). Quadriceps MTC force FMTC and its peak rate were estimated as knee torque Tknee and its peak rate divided by moment arm rMTC. SEE force FSEE was assumed to equal quadriceps MTC force FMTC.

The effects of knee torque Tknee on muscle fascicle length change ΔLfas and pennation angle change Δφ, as well as the effect of SEE force FSEE on SEE length change ΔLSEE were averaged across knee angles and fitted using exponentials. An exponential toe region was used for the SEE force–length curve, appropriate for the relatively low torque levels, less than 10% of maximal torque (Kulig et al., 1984; Lichtwark and Wilson, 2008). The mechanical work done on the SEE by the quadriceps muscle fascicles (Wfas) was quantified by the integral of the (fitted) SEE force FSEE with respect to the (fitted) SEE length LSEE, for the torques applied during the energetics trials. Fits were done for each subject individually; work estimates were averaged across subjects.

We hypothesized that overall metabolic cost should include contributions from muscle force rate , as well as from mechanical work and steady force production. The force-rate cost was hypothesized to be a consequence of additional active calcium transport by sarco-endoplasmic reticulum calcium ATPase (SERCA) owing to increased myoplasmic calcium concentration (see Fig. 2). Based on imaging and modeling studies (Baylor and Hollingworth, 2012), myoplasmic calcium concentration was expected to increase with force rate (Fig. 2B) owing to rate-limiting dynamics between calcium and muscle force production (Fig. 2A), which act like a low-pass filter. This means that greater amplitudes of calcium should be required to produce a given amount of force at faster rates. (Technically, active transport occurs with muscle deactivation, but only as a consequence of activation, hence our term ‘activation cost’.) This general concept is illustrated with an example time constant of 30 ms (Fig. 2B), typical for human force development (Maffiuletti et al., 2016).

The force-rate cost hypothesis may be expressed as a metabolic cost per contraction, and experimentally tested in terms of a metabolic rate (or cost per time). The cost per muscle contraction EFR should increase in proportion to how quickly the force F increases:
(1)
and particularly the peak force rate . For a sinusoidal force waveform of fixed amplitude, peak force rate is expected to increase with waveform frequency f:
(2)
Under steady-state conditions, the corresponding metabolic rate (Fig. 2C) depends on the cost per contraction EFR and on the frequency of contractions, yielding:
(3)
The work and force–time contributions should have separable dependencies on frequency f. Despite an isometric joint, muscle fascicles perform work against series elastic tissues and imperfectly rigid leg and dynamometer. For a fixed torque amplitude, work should be performed in a fixed amount per contraction, and the rate of mechanical work should therefore increase with the rate of contractions, or waveform frequency f. Assuming fixed biochemical costs for fascicle work, the corresponding metabolic rate for work production should increase as:
(4)
The force maintenance contribution is for energy expended even when muscle fascicles are isometric. Metabolic cost should be proportional to the force amplitude and duration of force (or force–time integral; Crow and Kushmerick, 1982). Owing to constant amplitude of the waveform, average knee torque and therefore average muscle force were expected to be independent of waveform frequency. The corresponding force–time metabolic rate should be independent of waveform frequency:
(5)
In our experiments, we tested for an overall metabolic cost with all three of these contributions. We hypothesized that net metabolic rate would be the sum of force rate, work and force–time integral terms:
(6)

where ‘net’ is defined as gross minus the cost of quiet sitting. We assumed differences in net metabolic rate to mainly reflect knee muscle energy consumption, as producing low forces should have negligible effects on metabolism in other muscles (e.g. cardiovascular) and organs. The expected cost of performing fascicle work is derived from the inverse product of the efficiencies of (1) obtaining ATP from food stuff (about 60%; van Ingen Schenau et al., 1997) and (2) crossbridge formation (about 50%; Barclay, 2015), yielding the linear coefficient c1 of 3.33. The other coefficients (c0 and c2) were determined using regression with f as the independent variable. Unless stated otherwise, standard deviations (s.d.) refer to between-subjects variability.

Prior to examining the metabolic rate of cyclic force production, we first examined the experimental conditions, which showed that participants matched target torques well. As expected, participants' knee torque Tknee (t) resembled a sinusoidal waveform, with constant average and variable peak rate . There was no significant dependency of waveform frequency f on average torque (slope 0.2±0.4 N m Hz−1, mean±95% CI, P=0.5, linear regression). Average knee torque , torque minima and maxima were 10.7±1.0, 1.0±1.0 and 21.8±1.2 N m, respectively (means±s.d. across subjects and conditions; Fig. 3A). Average knee torque was 11.6±0.7 N m (mean±s.d.) in the constant torque condition. Despite constant average torque , quadriceps EMG amplitude increased with waveform frequency by about 76% (VL: 2.1±0.5% MVC at 0.5 Hz to 3.7±0.6% MVC at 2.5 Hz, mean±s.d.). The increase was approximately linear in frequency f, at rates of 0.8±0.3, 1.0±0.3 and 0.2±0.1% MVC Hz−1 for VL, RF and BF, respectively (means±95% CI, repeated-measures linear regression). Averaged across conditions, VL EMG amplitude was 1.0±2.4% MVC higher in the left leg than in the right leg. There was no difference in BF EMG amplitude between legs (i.e. 0.0±0.3% MVC difference). There was some co-contraction, with BF EMG amplitude less than 1% MVC in all conditions (Fig. 3C), on average approximately one-eighth of VL EMG amplitudes. Peak knee torque rate increased with waveform frequency f (slope 68.3±3.4 N m, mean±95% CI, linear regression without offset), and closely resembled the increase expected from the torque targets (slope 69.1 N m, R2=0.84, Fig. 3B).

Quadriceps muscle fascicles performed work during isometric torque production, as evidenced by VL fascicle shortening (Fig. 4A) and pennation angle changes (Fig. 4B) observed during mechanics trials. All results reported here are averaged between legs, unless stated otherwise. As a baseline, VL fascicle length was estimated at 9.5±0.8 cm (mean±s.d.) during rest (20 deg knee flexion). Fascicle lengths decreased with greater torque, as described by a fitted exponential (R2=0.36; see Fig. 4A). At 11 N m, muscle fascicles had shortened by 1.4±0.6 cm. Pennation angle increased with greater torque, also described by a fitted exponential (R2=0.36; Fig. 4B). Pennation angle was 15.9±1.2 deg (mean±s.d.) at rest, and increased by 3.4±1.1 deg at 11 N m. Combining fascicle length and pennation angle changes with knee angle, quadriceps MTC moment arm rMTC was estimated at 5.0±1.9 cm (mean±s.d.). These data also yielded an estimate of length change of SEEs, fitted by an exponential toe region (R2=0.37; Fig. 4C). From that fit, quadriceps muscle fascicle work Wfas was estimated at 1.2±0.6 J (mean±s.d.) for an 11 N m contraction. The same estimate was obtained when exclusively using left leg data (1.2±0.8 J) or right leg data (1.2±0.7 J). Using the measured cyclic torque maxima (see Fig. 3A) instead of assuming 11 N m contractions resulted in similar work estimates (see Fig. 5A). The increase in mechanical work rate with waveform frequency was therefore approximately 2.4±1.2 J.

Net metabolic rate increased with waveform frequency f (P=4×10−6, repeated-measures regression; Fig. 5B), in agreement with predictions. The RER was below 0.8 in all conditions, indicating aerobic metabolism. Net metabolic rate increased from a value of 10.5±13.2 W (mean±s.d.) at the lowest frequency (f=0.5 Hz) to 46.8±8.3 W at the highest frequency (f=2.5 Hz). The overall cost appeared to agree with hypothesized contributions from force rate , work and force–time . At the highest frequency, these three costs were estimated at 18.9, 20.1 and 7.8 W, respectively, or approximately 40.4%, 43.0% and 16.6% of net metabolic rate . Costs of force rate and force–time had coefficients of 3.2±0.6 J s and 7.8±2.2 W, respectively (mean±95% CI, P=9×10−6 and P=0.001, respectively, linear mixed-effects regression). The force–time cost estimate was comparable to the metabolic cost of the constant torque condition (9.8±12.7 W, mean±s.d.). Examining the force-rate cost as a cost per contraction EFR (Fig. 5C), it increased linearly with measured peak force rate (2.6±1.9×10−3 m s, mean±95% CI, linear mixed-effects regression without offset, P=0.008; Fig. 5C).

The present study aimed at investigating the metabolic cost related to muscle force rate, separately from both force and mechanical work costs. Net metabolic rate increased with the frequency of isometric torque production, faster than would be expected for force production alone. We also estimated muscle fascicle shortening, and observed non-negligible work performed by fascicles, but also insufficient to explain metabolic cost. A considerable portion of the overall metabolic cost can be ascribed to the proposed force-rate cost. We next examined the hypothesis in relation to the data, the potential mechanisms for the force-rate cost, and implications for movement.

We first considered why force and work production alone do not explain the observed metabolic cost. The cyclic, isometric task resulted in nearly the same average force across conditions. The muscle force–time integral was therefore unchanged across frequencies, and likely to have the same metabolic cost in all conditions. The force–time cost estimate accounted for 74% of metabolic cost in the 0.5 Hz condition, as opposed to only 17% in the 2.5 Hz condition. Expressed as a cost per unit force, the force–time cost estimate (i.e. 0.038 W N−1) is well within the range of values reported previously (i.e. 0.020–0.058 W N–1, Ortega et al., 2015; van der Zee et al., 2019). Muscle co-contraction could have affected the force–time cost, as antagonist force needs to be offset by additional agonist force. We did observe antagonist EMG of relatively small and linearly increasing amplitude, but consider it unlikely to explain the quadratically increasing energetic cost with waveform frequency. The cost of steady isometric contraction has in part been ascribed to muscle crossbridges, which cycle even when net work is zero (Barclay et al., 2007), with energy dissipated as heat. The force–time coefficient c0 in part reflects such dissipative crossbridge cycling in isometric conditions. As shortening velocity increases, a smaller fraction of metabolic energy is dissipated as heat, while a greater fraction is converted into mechanical work (Barclay et al., 2010). This suggests that we may have overestimated the force–time coefficient c0 at higher waveform frequencies. Although the force–time cost did not increase with waveform frequency, the mechanical work rate did increase. We assumed a proportional work cost, implying a constant crossbridge efficiency. However, crossbridge efficiency is thought to increase with shortening velocity, peaking at approximately 15–25% of maximal shortening velocity (Smith et al., 2005). This increase is thought to reflect more work done per crossbridge cycle owing to an increased working stroke length with crossbridge strain similar to isometric values (Barclay et al., 2010), implying fewer active crossbridges per unit work. Assuming previously reported muscle optimal length and maximal shortening velocity (De Ruiter et al., 2000; van Soest et al., 1993), shortening velocities in the present study were always well below 15% of maximal shortening velocity. Consequently, if anything, crossbridge efficiency would have increased with waveform frequency. Therefore, we may have overestimated the work cost at higher waveform frequencies by assuming a fixed mechanical efficiency. Altogether, the assumptions we have applied would lead to overestimation of work and/or force–time costs, so that their sum might actually explain less than the 60% of net metabolic cost we have estimated.

This leaves a substantial cost explained by the hypothesized force-rate cost. The sharp increase in metabolic cost with muscle contraction frequency (Fig. 5) was consistent with the expected force-rate cost increasing with the square of contraction frequency. The hypothesis is based on the rate-limited dynamics between muscle excitation and force production, which acts as a low-pass filter of excitation and associated calcium release into the myoplasm. In the present experiment, increasing contraction frequencies called for higher force rates (, Fig. 3), which are expected to require more calcium release per muscle contraction. Active calcium pumping returns the calcium to the sarcoplasmic reticulum, with an associated metabolic cost (SERCA; Inoue et al., 2019). This contrasts with the energetic costs for force maintenance or work production, both associated with cycling of actomyosin crossbridges (actomyosin ATPase). Our hypothesis is that the combination of rate-limiting dynamics and calcium pumping imply a separate force-rate cost. This cost is in addition to the traditional force- and work-related costs, and could explain 40% (or more, if work and force–time costs are overestimated) of metabolic cost observed at the higher frequencies of the present experiment. In movement conditions that require cyclic muscle contraction, the force-rate cost could be quite substantial.

A potential source of uncertainty is imaging of VL fascicle length and pennation angle. Fascicle angle estimates obtained from the Hough transform have been found to correlate well with traditional (manual) estimates (Zhou et al., 2015, 2012; Zhou and Zheng, 2008). Different in the current algorithm was the (semi-)automated muscle thickness estimation from aponeurosis objects identified using feature detection. To avoid erroneous muscle thickness estimates, identified aponeurosis objects were visually inspected and corrected where necessary. Owing to the VL fascicles extending beyond the ultrasound image frame, we linearly extrapolated fascicles to an aponeurosis point outside the frame (Brennan et al., 2017). To avoid excess linear extrapolation, we selected conditions (20 deg knee angle) to keep fascicles at relatively short length (Marzilger et al., 2018). To ensure clear images, we used a slow waveform frequency for torque production. Our quadriceps muscle work estimate was the same when exclusively using left leg versus right leg data, suggesting good consistency across legs. The resulting estimates of fascicle shortening (approximately 1.4 cm) were comparable to manually obtained estimates under comparable conditions (Fukunaga et al., 1997; Ichinose et al., 1997).

We used knee torque, VL fascicle length and VL pennation angle to estimate quadriceps muscle work against SEEs during isometric torque production. We assumed that SEE length changes, as estimated from VL fascicle shortening, were shared across the quadriceps heads, which all have a common distal tendon. Each head also has its own aponeurosis, which we assumed would not drastically affect our results. Similar activation patterns have been observed between VL and vastus medialis by others (Wakeling and Horn, 2009). While activation can differ between rectus femoris and vastii in some conditions (Wakeling and Horn, 2009), we observed quite comparable EMGs for VL and rectus femoris for the conditions in the present study. We also assumed equal contributions of left and right leg to knee torque. While left and right leg EMG increased to similar extents with waveform frequency, their average values were different. However, even if there was a difference in contribution between legs, this would not have affected the work estimate because the work–force curve was approximately linear around the peak forces during the cyclic force experiment (see Fig. 4C, inset). Our ultrasound estimates can also be compared against two alternatives. One is to assume that muscle fascicles do not shorten against an isometric dynamometer, thus yielding zero work. Another is to consult the literature for SEE stress–strain curves. Assuming 4–5% SEE strain at maximal muscle force (e.g. Millard et al., 2013), a SEE slack length of 20 cm (based on VL and RF values in van Soest and Bobbert, 1993) and a typical toe region (Thelen, 2003), a contraction at 10% of maximal muscle force would perform 0.24 J of work against SEE, 5 times less than we reported. This discrepancy may in part reflect an overestimation of work in the present study. Altogether, potential inaccuracies in our quadriceps muscle work estimate would lead to an overestimate of work and its energetic cost. This would suggest that the hypothesized force-rate cost is greater than estimated here.

Evidence for the proposed force-rate cost has previously been demonstrated in both muscle and whole-body experiments. For cyclic isometric contraction of mouse muscle in situ (1:3 on:off duty factor), rate of oxygen consumption is higher when cycle frequency is higher (Hogan et al., 1998). For prolonged isometric contraction of human VL muscle in vivo, rate of ATP utilization is higher when stimulation is cyclic (1:1 on:off duty factor) compared with constant (Chasiotis et al., 1987; Spriet et al., 1988). For cyclic isometric contraction of human calf muscle at constant frequency, metabolic cost decreases with increasing duty factor (Beck et al., 2020). While these studies have indicated a metabolic cost for cyclic muscle contraction, this cost was not parametrically related to muscle force-rate. We believe that force rate, rather than work or force–time, could explain many such findings on metabolic cost of cyclic isometric contraction. In contrast to previous studies on isometric contraction varying on–off durations, we parametrically related such metabolic cost to muscle force rate, which we find more suitable for the continuously varying forces of most daily movements. Furthermore, metabolic cost of cyclic movement has been observed to increase with movement frequency in leg swinging (Doke et al., 2005; Doke and Kuo, 2007), ankle bouncing (Dean and Kuo, 2011) and arm reaching (Wong et al., 2020 preprint). Altogether, the force-rate cost seems to be observable during a variety of movements and contraction conditions. Here, we provide evidence that the force-rate cost exists during voluntary isometric contraction, independently of work and force–time costs.

We have separated total metabolic cost into three terms: work, force time and force rate. Such phenomenological division is potentially explained mechanistically as costs due to individual ATPases, particularly actomyosin (for work and force time) and SERCA (for force rate). Whether SERCA is primarily responsible for the force-rate cost could be tested in future experimental studies on isolated muscle, for example using crossbridge (and therefore actomyosin ATPase) blocking agents (Barclay et al., 2007). We expect SERCA ATP consumption rate to increase with contraction frequency because of the increase in calcium release amplitudes, owing to rate-limiting force production dynamics. Higher calcium amplitudes would require higher muscle activation levels, which agrees with the observed increase in quadriceps EMG amplitude with waveform frequency. Such properties are not included in Hill-type muscle models (Zajac, 1989), which do have activation dynamics attributed to calcium release (sometimes called ‘active state’) and muscle force–length and –velocity dependency, but not rate-limiting dynamics between calcium release and muscle force. As demonstration, we simulated the present experiment with a commonly used metabolic muscle model (OpenSim software, Uchida et al., 2016). It predicted a 2.4 W increase in metabolic rate with waveform frequency, only approximately 7% of the amount experimentally measured. Current models therefore appear to underestimate the energetic cost of cyclic force production. Improved estimates might be obtained by including rate-limiting dynamics, crudely modeled as a low-pass filter here, such as calcium release and binding (Baylor and Hollingworth, 2012) and Huxley-type crossbridge binding (van Soest et al., 2019). Future modeling studies could develop such mechanistic models that explain the force-rate cost.

Fast contraction of muscle may require high calcium release and pumping rates, at the cost of metabolic energy spend by SERCA. Such energetic penalty for rapid increase in muscle force may explain why (ground reaction) forces during human walking are smoother than expected based on work considerations (Rebula and Kuo, 2015). Nevertheless, there may be a considerable force-rate cost during human locomotion, as VL fascicles produce relatively short bursts of high force without shortening much (0–30% of stride; Bohm et al., 2018). The force-rate cost may also apply to a variety of other animals. In locust flight muscle, oxygen consumption during cyclic isometric contraction is 87% of the consumption during maximal power output, indicating a high cost for cyclic muscle activation (Josephson and Stevenson, 1991). Locusts use synchronous muscle with one contraction per wing beat, whereas other insects such as beetles use asynchronous muscle with fewer contractions per beat (Josephson et al., 2001), which could reduce the activation cost as a potentially advantageous adaptation (Syme and Josephson, 2002). A force-rate cost could be an important factor in energy expenditure in a wide range of movements and species.

Conclusions

We observed an increase in average net metabolic rate with frequency of cyclic torque production, which could not be explained by force–time nor by muscle fascicle mechanical work. Average net metabolic rate was related to force rate, suggesting that increasing muscle force abruptly requires more metabolic energy than when done slowly. We propose that this metabolic cost of force rate may be explained by an increase in the amount of required active calcium transport, and may be relevant for a wide range of movements.

The authors would like to thank Ashna Subramanium for assistance with data collection.

Author contributions

Conceptualization: T.J.v.d.Z., A.D.K.; Methodology: T.J.v.d.Z., A.D.K.; Software: T.J.v.d.Z.; Validation: T.J.v.d.Z., A.D.K.; Formal analysis: T.J.v.d.Z.; Investigation: T.J.v.d.Z.; Resources: A.D.K.; Data curation: T.J.v.d.Z.; Writing - original draft: T.J.v.d.Z.; Writing - review & editing: T.J.v.d.Z., A.D.K.; Visualization: T.J.v.d.Z., A.D.K.; Supervision: A.D.K.; Project administration: A.D.K.; Funding acquisition: A.D.K.

Funding

This work was supported in part by the Natural Sciences and Engineering Research Council of Canada (NSERC Discovery and Canada Research Chair, Tier 1) and the Dr Benno Nigg Research Chair in Biomechanics.

Data availability

Ultrasound algorithm and typical example of ultrasound data are available from github at: https://github.com/timvanderzee/ultrasound-automated-algorithm

Abbott
,
B. C.
,
Bigland
,
B.
and
Ritchie
,
J. M.
(
1952
).
The physiological cost of negative work
.
J. Physiol.
117
,
380
-
390
.
Barclay
,
C. J.
(
2015
).
Energetics of contraction
.
Compr. Physiol.
5
,
961
-
995
.
Barclay
,
C. J.
,
Woledge
,
R. C.
and
Curtin
,
N. A.
(
2007
).
Energy turnover for Ca2+ cycling in skeletal muscle
.
J. Muscle Res. Cell Motil.
28
,
259
-
274
.
Barclay
,
C. J.
,
Woledge
,
R. C.
and
Curtin
,
N. A.
(
2010
).
Inferring crossbridge properties from skeletal muscle energetics
.
Prog. Biophys. Mol. Biol.
102
,
53
-
71
.
Baylor
,
S. M.
and
Hollingworth
,
S.
(
2012
).
Intracellular calcium movements during excitation–contraction coupling in mammalian slow-twitch and fast-twitch muscle fibers
.
J. Gen. Physiol.
139
,
261
-
272
.
Beck
,
O. N.
,
Gosyne
,
J.
,
Franz
,
J. R.
and
Sawicki
,
G. S.
(
2020
).
Cyclically producing the same average muscle-tendon force with a smaller duty increases metabolic rate
.
Proc. R. Soc. B
287
,
20200431
.
Bergström
,
M.
and
Hultman
,
E.
(
1988
).
Energy cost and fatigue during intermittent electrical stimulation of human skeletal muscle
.
J. Appl. Physiol.
65
,
1500
-
1505
.
Bohm
,
S.
,
Marzilger
,
R.
,
Mersmann
,
F.
,
Santuz
,
A.
and
Arampatzis
,
A.
(
2018
).
Operating length and velocity of human vastus lateralis muscle during walking and running
.
Sci. Rep.
8
,
5066
.
Brennan
,
S. F.
,
Cresswell
,
A. G.
,
Farris
,
D. J.
and
Lichtwark
,
G. A.
(
2017
).
In vivo fascicle length measurements via B-mode ultrasound imaging with single vs dual transducer arrangements
.
J. Biomech.
64
,
240
-
244
.
Brockway
,
J. M.
(
1987
).
Derivation of formulae used to calculate energy expenditure in man
.
Hum. Nutr. Clin. Nutr.
41
,
463
-
471
.
Chasiotis
,
D.
,
Bergström
,
M.
and
Hultman
,
E.
(
1987
).
ATP utilization and force during intermittent and continuous muscle contractions
.
J. Appl. Physiol.
63
,
167
-
174
.
Crow
,
M. T.
and
Kushmerick
,
M. J.
(
1982
).
Chemical energetics of slow- and fast-twitch muscles of the mouse
.
J. Gen. Physiol.
79
,
147
-
166
.
De Ruiter
,
C. J.
,
Didden
,
W. J. M.
,
Jones
,
D. A.
and
De Haan
,
A.
(
2000
).
The force–velocity relationship of human adductor pollicis muscle during stretch and the effects of fatigue
.
J. Physiol.
526
,
671
-
681
.
Dean
,
J. C.
and
Kuo
,
A. D.
(
2011
).
Energetic costs of producing muscle work and force in a cyclical human bouncing task
.
J. Appl. Physiol.
110
,
873
-
880
.
Doke
,
J.
and
Kuo
,
A. D.
(
2007
).
Energetic cost of producing cyclic muscle force, rather than work, to swing the human leg
.
J. Exp. Biol.
210
,
2390
-
2398
.
Doke
,
J.
,
Donelan
,
J. M.
and
Kuo
,
A. D.
(
2005
).
Mechanics and energetics of swinging the human leg
.
J. Exp. Biol.
208
,
439
-
445
.
Fukunaga
,
T.
,
Ichinose
,
Y.
,
Ito
,
M.
,
Kawakami
,
Y.
and
Fukashiro
,
S.
(
1997
).
Determination of fascicle length and pennation in a contracting human muscle in vivo
.
J. Appl. Physiol.
82
,
354
-
358
.
Fukunaga
,
T.
,
Kubo
,
K.
,
Kawakami
,
Y.
,
Fukashiro
,
S.
,
Kanehisa
,
H.
and
Maganaris
,
C. N.
(
2001
).
In vivo behaviour of human muscle tendon during walking
.
Proc. R. Soc. B Biol. Sci.
268
,
229
-
233
.
Hof
,
A. L.
(
1984
).
EMG and muscle force: an introduction
.
Hum. Mov. Sci.
3
,
119
-
153
.
Hogan
,
M. C.
,
Ingham
,
E.
and
Kurdak
,
S. S.
(
1998
).
Contraction duration affects metabolic energy cost and fatigue in skeletal muscle
.
Am. J. Physiol.
274
,
E397
-
E402
.
Huang
,
N. E.
,
Shen
,
Z.
,
Long
,
S. R.
,
Wu
,
M. C.
,
Shih
,
H. H.
,
Zheng
,
Q.
,
Yen
,
N.-C.
,
Tung
,
C. C.
and
Liu
,
H. H.
(
1998
).
The empirical mode decomposition and the Hilbert spectrum for nonlinear and non-stationary time series analysis
.
Proc. Math. Phys. Eng. Sci.
454
,
903
-
995
.
Huang
,
H. J.
,
Kram
,
R.
and
Ahmed
,
A. A.
(
2012
).
Reduction of metabolic cost during motor learning of arm reaching dynamics
.
J. Neurosci.
32
,
2182
-
2190
.
Ichinose
,
Y.
,
Kawakami
,
Y.
,
Ito
,
M.
and
Fukunaga
,
T.
(
1997
).
Estimation of active force–length characteristics of human vastus lateralis muscle
.
Cells Tissues Organs
159
,
78
-
83
.
Inoue
,
M.
,
Sakuta
,
N.
,
Watanabe
,
S.
,
Zhang
,
Y.
,
Yoshikaie
,
K.
,
Tanaka
,
Y.
,
Ushioda
,
R.
,
Kato
,
Y.
,
Takagi
,
J.
,
Tsukazaki
,
T.
et al. 
(
2019
).
Structural basis of Sarco/Endoplasmic Reticulum Ca2+-ATPase 2b regulation via transmembrane helix interplay
.
Cell Rep.
27
,
1221
-
1230.e3
.
Ito
,
M.
,
Akima
,
H.
and
Fukunaga
,
T.
(
2000
).
In vivo moment arm determination using B-mode ultrasonography
.
J. Biomech.
33
,
215
-
218
.
Josephson
,
R. K.
and
Stevenson
,
R. D.
(
1991
).
The efficiency of a flight muscle from the locust Schistocerca americana
.
J. Physiol.
442
,
413
-
429
.
Josephson
,
R. K.
,
Malamud
,
J. G.
and
Stokes
,
D. R.
(
2001
).
The efficiency of an asynchronous flight muscle from a beetle
.
J. Exp. Biol.
204
,
4125
-
4139
.
Kulig
,
K.
,
Andrews
,
J. G.
and
Hay
,
J. G.
(
1984
).
Human strength curves
.
Exerc. Sport Sci. Rev.
12
,
417
-
466
.
La Delfa
,
N. J.
,
Sutherland
,
C. A.
and
Potvin
,
J. R.
(
2014
).
EMG processing to interpret a neural tension-limiting mechanism with fatigue
.
Muscle Nerve
50
,
384
-
392
.
Lichtwark
,
G. A.
and
Wilson
,
A. M.
(
2008
).
Optimal muscle fascicle length and tendon stiffness for maximising gastrocnemius efficiency during human walking and running
.
J. Theor. Biol.
252
,
662
-
673
.
Maffiuletti
,
N. A.
,
Aagaard
,
P.
,
Blazevich
,
A. J.
,
Folland
,
J.
,
Tillin
,
N.
and
Duchateau
,
J.
(
2016
).
Rate of force development: physiological and methodological considerations
.
Eur. J. Appl. Physiol.
116
,
1091
-
1116
.
Margaria
,
R.
(
1968
).
Positive and negative work performances and their efficiencies in human locomotion
.
Int. Z. Angew. Physiol. Einschl. Arbeitsphysiol.
25
,
339
-
351
.
Margaria
,
R.
(
1976
).
Biomechanics and Energetics of Muscular Exercise
.
Oxford: Clarendon Press
.
Marzilger
,
R.
,
Legerlotz
,
K.
,
Panteli
,
C.
,
Bohm
,
S.
and
Arampatzis
,
A.
(
2018
).
Reliability of a semi-automated algorithm for the vastus lateralis muscle architecture measurement based on ultrasound images
.
Eur. J. Appl. Physiol.
118
,
291
-
301
.
Millard
,
M.
,
Uchida
,
T.
,
Seth
,
A.
and
Delp
,
S. L.
(
2013
).
Flexing computational muscle: modeling and simulation of musculotendon dynamics
.
J. Biomech. Eng.
135
,
021005
.
Ortega
,
J. O.
,
Lindstedt
,
S. L.
,
Nelson
,
F. E.
,
Jubrias
,
S. A.
,
Kushmerick
,
M. J.
and
Conley
,
K. E.
(
2015
).
Muscle force, work and cost: a novel technique to revisit the Fenn effect
.
J. Exp. Biol.
218
,
2075
-
2082
.
Rebula
,
J. R.
and
Kuo
,
A. D.
(
2015
).
The cost of leg forces in bipedal locomotion: a simple optimization study
.
PLoS ONE
10
,
e0117384
.
Ryan
,
D. S.
,
Stutzig
,
N.
,
Siebert
,
T.
and
Wakeling
,
J. M.
(
2019
).
Passive and dynamic muscle architecture during transverse loading for gastrocnemius medialis in man
.
J. Biomech.
86
,
160
-
166
.
Smith
,
N. P.
,
Barclay
,
C. J.
and
Loiselle
,
D. S.
(
2005
).
The efficiency of muscle contraction
.
Prog. Biophys. Mol. Biol.
88
,
1
-
58
.
Spriet
,
L. L.
,
Soderlund
,
K.
and
Hultman
,
E.
(
1988
).
Energy cost and metabolic regulation during intermittent and continuous tetanic contractions in human skeletal muscle
.
Can. J. Physiol. Pharmacol.
66
,
134
-
139
.
Syme
,
D. A.
and
Josephson
,
R. K.
(
2002
).
How to build fast muscles: synchronous and asynchronous designs
.
Integr. Comp. Biol.
42
,
762
-
770
.
Thelen
,
D. G.
(
2003
).
Adjustment of muscle mechanics model parameters to simulate dynamic contractions in older adults
.
J. Biomech. Eng.
125
,
70
-
77
.
Uchida
,
T. K.
,
Hicks
,
J. L.
,
Dembia
,
C. L.
and
Delp
,
S. L.
(
2016
).
Stretching your energetic budget: how tendon compliance affects the metabolic cost of running
.
PLoS ONE
11
,
e0150378
.
van der Zee
,
T. J.
and
Kuo
,
A. D.
(
2020
).
Fully automated algorithm estimates muscle fascicle length from ultrasound image
.
BioRxiv
,
2020.08.23.263574
.
van der Zee
,
T. J.
,
Lemaire
,
K. K.
and
van Soest
,
A. J.
(
2019
).
The metabolic cost of in vivo constant muscle force production at zero net mechanical work
.
J. Exp. Biol.
222
,
jeb199158
.
van Ingen Schenau
,
G. J.
,
Bobbert
,
M. F.
and
de Haan
,
A.
(
1997
).
Does elastic energy enhance work and efficiency in the stretch–shortening cycle?
J. Appl. Biomech.
13
,
389
-
415
.
van Soest
,
A. J.
and
Bobbert
,
M. F.
(
1993
).
The contribution of muscle properties in the control of explosive movements
.
Biol. Cybern.
69
,
195
-
204
.
van Soest
,
A. J.
,
Schwab
,
A. L.
,
Bobbert
,
M. F.
and
van Ingen Schenau
,
G. J.
(
1993
).
The influence of the biarticularity of the gastrocnemius muscle on vertical-jumping achievement
.
J. Biomech.
26
,
1
-
8
.
van Soest
,
A. J. K.
,
Casius
,
L. J. R.
and
Lemaire
,
K. K.
(
2019
).
Huxley-type cross-bridge models in largeish-scale musculoskeletal models; an evaluation of computational cost
.
J. Biomech.
83
,
43
-
48
.
Wakeling
,
J. M.
and
Horn
,
T.
(
2009
).
Neuromechanics of muscle synergies during cycling
.
J. Neurophysiol. 101, 843-854
.
Wong
,
J. D.
,
Cluff
,
T.
and
Kuo
,
A. D.
(
2020
).
The energetic basis for smooth human arm movements
.
BioRxiv
.
Xu
,
F.
and
Rhodes
,
E. C.
(
1999
).
Oxygen uptake kinetics during exercise
.
Sports Med.
27
,
313
-
327
.
Zajac
,
F. E.
(
1989
).
Muscle and tendon: properties, models, scaling, and application to biomechanics and motor control
.
Crit. Rev. Biomed. Eng.
17
,
359
-
411
.
Zhou
,
Y.
and
Zheng
,
Y.-P.
(
2008
).
Estimation of muscle fiber orientation in ultrasound images using revoting Hough transform (RVHT)
.
Ultrasound Med. Biol.
34
,
1474
-
1481
.
Zhou
,
Y.
,
Li
,
J.-Z.
,
Zhou
,
G.
and
Zheng
,
Y.-P.
(
2012
).
Dynamic measurement of pennation angle of gastrocnemius muscles during contractions based on ultrasound imaging
.
Biomed. Eng. Online
11
,
63
.
Zhou
,
G.-Q.
,
Chan
,
P.
and
Zheng
,
Y.-P.
(
2015
).
Automatic measurement of pennation angle and fascicle length of gastrocnemius muscles using real-time ultrasound imaging
.
Ultrasonics
57
,
72
-
83
.

Competing interests

The authors declare no competing or financial interests.