Metabolic cost underlies task-dependent variations in motor unit recruitment

Mammalian skeletal muscles are comprised of many motor units, each containing a group of muscle fibres that have common contractile properties: these can be broadly categorized as slow and fast twitch muscle fibres. Motor units are typically recruited in an orderly fashion following the 'size principle', in which slower motor units would be recruited for low intensity contraction; a metabolically cheap and fatigue-resistant strategy. However, this recruitment strategy poses a mechanical paradox for fast, low intensity contractions, in which the recruitment of slower fibres, as predicted by the size principle, would be metabolically more costly than the recruitment of faster fibres that are more efficient at higher contraction speeds. Hence, it would be mechanically and metabolically more effective for recruitment strategies to vary in response to contraction speed so that the intrinsic efficiencies and contraction speeds of the recruited muscle fibres are matched to the mechanical demands of the task. In this study, we evaluated the effectiveness of a novel, mixed cost function within a musculoskeletal simulation, which includes the metabolic cost of contraction, to predict the recruitment of different muscle fibre types across a range of loads and speeds. Our results show that a metabolically informed cost function predicts favoured recruitment of slower muscle fibres for slower and isometric tasks versus recruitment that favours faster muscles fibres for higher velocity contractions. This cost function predicts a change in recruitment patterns consistent with experimental observations, and also predicts a less expensive metabolic cost for these muscle contractions regardless of speed of the movement. Hence, our findings support the premise that varying motor recruitment strategies to match the mechanical demands of a movement task results in a mechanically and metabolically sensible way to deploy the different types of motor unit.

Movement Complexity and Neuromechanical Factors Affect the Entropic Half-Life of Myoelectric Signals

Appropriate neuromuscular functioning is essential for survival and features underpinning motor control are present in myoelectric signals recorded from skeletal muscles. One approach to quantify control processes related to function is to assess signal variability using measures such as Sample Entropy. Here we developed a theoretical framework to simulate the effect of variability in burst duration, activation duty cycle, and intensity on the Entropic Half-Life (EnHL) in myoelectric signals. EnHLs were predicted to be <40 ms, and to vary with fluctuations in myoelectric signal amplitude and activation duty cycle. Comparison with myoelectic data from rats walking and running at a range of speeds and inclines confirmed the range of EnHLs, however, the direction of EnHL change in response to altered locomotor demand was not correctly predicted. The discrepancy reflected different associations between the ratio of the standard deviation and mean signal intensity (Ist:It¯) and duty factor in simulated and physiological data, likely reflecting additional information in the signals from the physiological data (e.g., quiescent phase content; variation in action potential shapes). EnHL could have significant value as a novel marker of neuromuscular responses to alterations in perceived locomotor task complexity and intensity.

Validation of Hill-type muscle models in relation to neuromuscular recruitment and force-velocity properties: predicting patterns of in vivo muscle force

We review here the use and reliability of Hill-type muscle models to predict muscle performance under varying conditions, ranging from in situ production of isometric force to in vivo dynamics of muscle length change and force in response to activation. Muscle models are frequently used in musculoskeletal simulations of movement, particularly when applied to studies of human motor performance in which surgically implanted transducers have limited use. Musculoskeletal simulations of different animal species also are being developed to evaluate comparative and evolutionary aspects of locomotor performance. However, such models are rarely validated against direct measures of fascicle strain or recordings of muscle-tendon force. Historically, Hill-type models simplify properties of whole muscle by scaling salient properties of single fibers to whole muscles, typically accounting for a muscle's architecture and series elasticity. Activation of the model's single contractile element (assigned the properties of homogenous fibers) is also simplified and is often based on temporal features of myoelectric (EMG) activation recorded from the muscle. Comparison of standard one-element models with a novel two-element model and with in situ and in vivo measures of EMG, fascicle strain, and force recorded from the gastrocnemius muscles of goats shows that a two-element Hill-type model, which allows independent recruitment of slow and fast units, better predicts temporal patterns of in situ and in vivo force. Recruitment patterns of slow/fast units based on wavelet decomposition of EMG activity in frequency-time space are generally correlated with the intensity spectra of the EMG signals, the strain rates of the fascicles, and the muscle-tendon forces measured in vivo, with faster units linked to greater strain rates and to more rapid forces. Using direct measures of muscle performance to further test Hill-type models, whether traditional or more complex, remains critical for establishing their accuracy and essential for verifying their applicability to scientific and clinical studies of musculoskeletal function.

The effect of fast and slow motor unit activation on whole-muscle mechanical performance: the size principle may not pose a mechanical paradox

The output of skeletal muscle can be varied by selectively recruiting different motor units. However, our knowledge of muscle function is largely derived from muscle in which all motor units are activated. This discrepancy may limit our understanding of in vivo muscle function. Hence, this study aimed to characterize the mechanical properties of muscle with different motor unit activation. We determined the isometric properties and isotonic force-velocity relationship of rat plantaris muscles in situ with all of the muscle active, 30% of the muscle containing predominately slower motor units active or 20% of the muscle containing predominately faster motor units active. There was a significant effect of active motor unit type on isometric force rise time (p < 0.001) and the force-velocity relationship (p < 0.001). Surprisingly, force rise time was longer and maximum shortening velocity higher when all motor units were active than when either fast or slow motor units were selectively activated. We propose this is due to the greater relative effects of factors such as series compliance and muscle resistance to shortening during sub-maximal contractions. The findings presented here suggest that recruitment according to the size principle, where slow motor units are activated first and faster ones recruited as demand increases, may not pose a mechanical paradox, as has been previously suggested.

Early deactivation of slower muscle fibres at high movement frequencies

Animals produce rapid movements using fast cyclical muscle contractions. These types of movements are better suited to faster muscle fibres within muscles of mixed fibre types as they can shorten at faster velocities and achieve higher activation-deactivation rates than their slower counterparts. Preferential recruitment of faster muscle fibres has previously been shown during high velocity contractions. Additionally, muscle deactivation takes longer than activation and therefore may pose a limitation to fast cyclical contractions. It has been speculated that slower fibres maybe deactivated before faster fibres to accommodate their longer deactivation time. This study aimed to test whether shifts in muscle fibre recruitment occur with derecruitment of slow fibres before the faster fibres at high cycle frequencies. Electromyographic (EMG) signals were collected from the medial gastrocnemius at an extreme range of cycle frequencies and workloads. Wavelets were used to resolve the EMG signals into time and frequency space and the primary sources of variability within the EMG frequency spectra were identified through principal component analysis. A general early derecruitment of slower fibres was evident at the end of muscle excitation for the higher cycle frequencies, and additional slower fibre recruitment was present at the highest cycle frequency. The duration of muscle excitation reached a minimum of about 150 ms and did not change for the three highest cycle frequencies suggesting a duration limit for the medial gastrocnemius. This study provides further evidence of modifications of muscle fibre recruitment strategies to meet the mechanical demands of movement.

Accuracy of gastrocnemius muscles forces in walking and running goats predicted by one-element and two-element Hill-type models

Hill-type models are commonly used to estimate muscle forces during human and animal movement-yet the accuracy of the forces estimated during walking, running, and other tasks remains largely unknown. Further, most Hill-type models assume a single contractile element, despite evidence that faster and slower motor units, which have different activation-deactivation dynamics, may be independently or collectively excited. This study evaluated a novel, two-element Hill-type model with "differential" activation of fast and slow contractile elements. Model performance was assessed using a comprehensive data set (including measures of EMG intensity, fascicle length, and tendon force) collected from the gastrocnemius muscles of goats during locomotor experiments. Muscle forces predicted by the new two-element model were compared to the forces estimated using traditional one-element models and to the forces measured in vivo using tendon buckle transducers. Overall, the two-element model resulted in the best predictions of in vivo gastrocnemius force. The coefficient of determination, r(2), was up to 26.9% higher and the root mean square error, RMSE, was up to 37.4% lower for the two-element model than for the one-element models tested. All models captured salient features of the measured muscle force during walking, trotting, and galloping (r(2)=0.26-0.51), and all exhibited some errors (RMSE=9.63-32.2% of the maximum in vivo force). These comparisons provide important insight into the accuracy of Hill-type models. The results also show that incorporation of fast and slow contractile elements within muscle models can improve estimates of time-varying, whole muscle force during locomotor tasks.

Recruitment of faster motor units is associated with greater rates of fascicle strain and rapid changes in muscle force during locomotion

Animals modulate the power output needed for different locomotor tasks by changing muscle forces and fascicle strain rates. To generate the necessary forces, appropriate motor units must be recruited. Faster motor units have faster activation-deactivation rates than slower motor units, and they contract at higher strain rates; therefore, recruitment of faster motor units may be advantageous for tasks that involve rapid movements or high rates of work. This study identified motor unit recruitment patterns in the gastrocnemii muscles of goats and examined whether faster motor units are recruited when locomotor speed is increased. The study also examined whether locomotor tasks that elicit faster (or slower) motor units are associated with increased (or decreased) in vivo tendon forces, force rise and relaxation rates, fascicle strains and/or strain rates. Electromyography (EMG), sonomicrometry and muscle-tendon force data were collected from the lateral and medial gastrocnemius muscles of goats during level walking, trotting and galloping and during inclined walking and trotting. EMG signals were analyzed using wavelet and principal component analyses to quantify changes in the EMG frequency spectra across the different locomotor conditions. Fascicle strain and strain rate were calculated from the sonomicrometric data, and force rise and relaxation rates were determined from the tendon force data. The results of this study showed that faster motor units were recruited as goats increased their locomotor speeds from level walking to galloping. Slow inclined walking elicited EMG intensities similar to those of fast level galloping but different EMG frequency spectra, indicating that recruitment of the different motor unit types depended, in part, on characteristics of the task. For the locomotor tasks and muscles analyzed here, recruitment patterns were generally associated with in vivo fascicle strain rates, EMG intensity and tendon force. Together, these data provide new evidence that changes in motor unit recruitment have an underlying mechanical basis, at least for certain locomotor tasks.

Task-dependent activity of motor unit populations in feline ankle extensor muscles

Understanding the functional significance of the morphological diversity of mammalian skeletal muscles is limited by technical difficulties of estimating the contribution of motor units with different properties to unconstrained motor behaviours. Recently developed wavelet and principal components analysis of intramuscular myoelectric signals has linked signals with lower and higher frequency contents to the use of slower and faster motor unit populations. In this study we estimated the relative contributions of lower and higher frequency signals of cat ankle extensors (soleus, medial and lateral gastrocnemii, plantaris) during level, downslope and upslope walking and the paw-shake response. This was done using the first two myoelectric signal principal components (PCI, PCII), explaining over 90% of the signal, and an angle θ, a function of PCI/PCII, indicating the relative contribution of slower and faster motor unit populations. Mean myoelectric frequencies in all walking conditions were lowest for slow soleus (234 Hz) and highest for fast gastrocnemii (307 and 330 Hz) muscles. Motor unit populations within and across the studied muscles that demonstrated lower myoelectric frequency (suggesting slower populations) were recruited during tasks and movement phases with lower mechanical demands on the ankle extensors--during downslope and level walking and in early walking stance and paw-shake phases. With increasing mechanical demands (upslope walking, mid-phase of paw-shake cycles), motor unit populations generating higher frequency signals (suggesting faster populations) contributed progressively more. We conclude that the myoelectric frequency contents within and between feline ankle extensors vary across studied motor behaviours, with patterns that are generally consistent with muscle fibre-type composition.

A muscle's force depends on the recruitment patterns of its fibers

Biomechanical models of whole muscles commonly used in simulations of musculoskeletal function and movement typically assume that the muscle generates force as a scaled-up muscle fiber. However, muscles are comprised of motor units that have different intrinsic properties and that can be activated at different times. This study tested whether a muscle model comprised of motor units that could be independently activated resulted in more accurate predictions of force than traditional Hill-type models. Forces predicted by the models were evaluated by direct comparison with the muscle forces measured in situ from the gastrocnemii in goats. The muscle was stimulated tetanically at a range of frequencies, muscle fiber strains were measured using sonomicrometry, and the activation patterns of the different types of motor unit were calculated from electromyographic recordings. Activation patterns were input into five different muscle models. Four models were traditional Hill-type models with different intrinsic speeds and fiber-type properties. The fifth model incorporated differential groups of fast and slow motor units. For all goats, muscles and stimulation frequencies the differential model resulted in the best predictions of muscle force. The in situ muscle output was shown to depend on the recruitment of different motor units within the muscle.

Movement mechanics as a determinate of muscle structure, recruitment and coordination

During muscle contractions, the muscle fascicles may shorten at a rate different from the muscle-tendon unit, and the ratio of these velocities is its gearing. Appropriate gearing allows fascicles to reduce their shortening velocities and allows them to operate at effective shortening velocities across a range of movements. Gearing of the muscle fascicles within the muscle belly is the result of rotations of the fascicles and bulging of the belly. Variable gearing can also occur as a result of tendon length changes that can be caused by changes in the relative timing of muscle activity for different mechanical tasks. Recruitment patterns of slow and fast fibres are crucial for achieving optimal muscle performance, and coordination between muscles is related to whole limb performance. Poor coordination leads to inefficiencies and loss of power, and optimal coordination is required for high power outputs and high mechanical efficiencies from the limb. This paper summarizes key studies in these areas of neuromuscular mechanics and results from studies where we have tested these phenomena on a cycle ergometer are presented to highlight novel insights. The studies show how muscle structure and neural activation interact to generate smooth and effective motion of the body.

EMG analysis tuned for determining the timing and level of activation in different motor units

Recruitment patterns and activation dynamics of different motor units greatly influence the temporal pattern and magnitude of muscle force development, yet these features are not often considered in muscle models. The purpose of this study was to characterize the recruitment and activation dynamics of slow and fast motor units from electromyographic (EMG) recordings and twitch force profiles recorded directly from animal muscles. EMG and force data from the gastrocnemius muscles of seven goats were recorded during in vivo tendon-tap reflex and in situ nerve stimulation experiments. These experiments elicited EMG signals with significant differences in frequency content (p<0.001). The frequency content was characterized using wavelet and principal components analysis, and optimized wavelets with centre frequencies, 149.94 Hz and 323.13 Hz, were obtained. The optimized wavelets were used to calculate the EMG intensities and, with the reconstructed twitch force profiles, to derive transfer functions for slow and fast motor units that estimate the activation state of the muscle from the EMG signal. The resulting activation-deactivation time constants gave r values of 0.98-0.99 between the activation state and the force profiles. This work establishes a framework for developing improved muscle models that consider the intrinsic properties of slow and fast fibres within a mixed muscle, and that can more accurately predict muscle force output from EMG.

The influence of strain and activation on the locomotor function of rat ankle extensor muscles

The ankle extensor muscles of the rat have different mechanical and physiological properties, providing a means of studying how changes in locomotor demands influence muscle fascicle behaviour, force and mechanical power output in different populations of muscle fibre types. Muscle fascicle strain, strain rate and activation patterns in the soleus, plantaris and medial gastrocnemius muscles of the rat were quantified from sonomicrometric and myoelectric data, collected during treadmill locomotion under nine velocity/incline conditions. Significant differences in peak-to-peak muscle fascicle strains and strain rates were identified between the three muscles (P&lt;0.001, all cases), with much smaller strains (&lt;0.1) and strain rates (&lt;0.5 s−1) occurring in soleus and plantaris compared with medial gastrocnemius (&gt;0.2 and &gt;1.0 s−1, respectively). The proportion of stride duration that each muscle was active (duty cycle) differed between locomotor conditions as did the timing of the activation and deactivation phases. A simple Hill-based muscle model was used to determine the influence of muscle activation relative to maximum fascicle strain and duty cycle on total force production and mechanical power output, from a slow and a fast muscle fibre, simulated through two peak-to-peak strain cycles (0.1 and 0.3). The predictions of the model did not complement conclusions that may be drawn from the observation of myoelectric timing and fascicle strain trajectories in each of the muscles. The model predicted that changes in mechanical power output were more sensitive to changes in activation parameters than to changes in strain trajectories, with subtle changes in activation phase and duty cycle significantly affecting predicted mechanical power output.

Motor unit recruitment patterns 1: responses to changes in locomotor velocity and incline

Mammalian skeletal muscles are composed of a mixture of motor unit types, which contribute a range of mechanical and physiological properties to the muscle. For a muscle to effectively contribute to smooth, co-ordinated movement it must activate an appropriate number and combination of motor units to generate the required force over a suitable time period. Much evidence exists indicating that motor units are activated in an orderly fashion, from the slowest through to the fastest. A growing body of evidence, however, indicates that such a recruitment strategy does not always hold true. Here we investigate how motor unit recruitment patterns were influenced by changes in locomotor velocity and incline. Kinematics data and myoelectric signals were collected from three rat ankle extensor muscles during running on a treadmill at nine velocity and incline combinations. Wavelet and principal component analysis were used to simultaneously decompose the signals into time and frequency space. The relative frequency components of the signals were quantified during 20 time windows of a stride from each locomotor condition. Differences in signal frequency components existed between muscles and locomotor conditions. Faster locomotor velocities led to a relative increase in high frequency components, whereas greater inclines led to a relative increase in the low frequency components. These data were interpreted as representing changes in motor unit recruitment patterns in response to changes in the locomotor demand. Motor units were not always recruited in an orderly manner, indicating that recruitment is a multi-factorial phenomenon that is not yet fully understood.