The maximum shortening velocity of muscle should be scaled with activation

A Biomechanical Simulation of Forearm Flexion Using the Finite Element Approach

Upper limb movement is vital in daily life. A biomechanical simulation of the forearm with consideration of the physiological characteristics of the muscles is instrumental in gaining deeper insights into the upper limb motion mechanisms. In this study, we established a finite element model of the forearm, including the radius, biceps brachii, and tendons. We simulated the motion of the forearm resulting from the contraction of the biceps brachii by using a Hill-type transversely isotropic hyperelastic muscle model. We adjusted the contraction velocity of the biceps brachii muscle in the simulation and found that a slower muscle contraction velocity facilitated forearm flexion. Then, we changed the percentage of fast-twitch fibers, the maximum muscle strength, and the neural excitation values of the biceps brachii muscle to investigate the forearm flexion of elderly individuals. Our results indicated that reduced fast-twitch fiber percentage, maximum muscle strength, and neural excitation contributed to the decline in forearm motion capability in elderly individuals. Additionally, there is a threshold for neural excitation, below which, motion capability sharply declines. Our model aids in understanding the role of the biceps brachii in forearm flexion and identifying the causes of upper limb movement disorders, which is able to provide guidance for enhancing upper limb performance.

More than energy cost: multiple benefits of the long Achilles tendon in human walking and running

Elastic strain energy that is stored and released from long, distal tendons such as the Achilles during locomotion allows for muscle power amplification as well as for reduction of the locomotor energy cost: as distal tendons perform mechanical work during recoil, plantar flexor muscle fibres can work over smaller length ranges, at slower shortening speeds, and at lower activation levels. Scant evidence exists that long distal tendons evolved in humans (or were retained from our more distant Hominoidea ancestors) primarily to allow high muscle-tendon power outputs, and indeed we remain relatively powerless compared to many other species. Instead, the majority of evidence suggests that such tendons evolved to reduce total locomotor energy cost. However, numerous additional, often unrecognised, advantages of long tendons may speculatively be of greater evolutionary advantage, including the reduced limb inertia afforded by shorter and lighter muscles (reducing proximal muscle force requirement), reduced energy dissipation during the foot-ground collisions, capacity to store and reuse the muscle work done to dampen the vibrations triggered by foot-ground collisions, reduced muscle heat production (and thus core temperature), and attenuation of work-induced muscle damage. Cumulatively, these effects should reduce both neuromotor fatigue and sense of locomotor effort, allowing humans to choose to move at faster speeds for longer. As these benefits are greater at faster locomotor speeds, they are consistent with the hypothesis that running gaits used by our ancestors may have exerted substantial evolutionary pressure on Achilles tendon length. The long Achilles tendon may therefore be a singular adaptation that provided numerous physiological, biomechanical, and psychological benefits and thus influenced behaviour across multiple tasks, both including and additional to locomotion. While energy cost may be a variable of interest in locomotor studies, future research should consider the broader range of factors influencing our movement capacity, including our decision to move over given distances at specific speeds, in order to understand more fully the effects of Achilles tendon function as well as changes in this function in response to physical activity, inactivity, disuse and disease, on movement performance.

How do differences in Achilles' tendon moment arm lengths affect muscle-tendon dynamics and energy cost during running?

Introduction

The relationship between the Achilles tendon moment arm length (AT_MA) and the energy cost of running (E_run) has been disputed. Some studies suggest a short AT_MA reduces E_run while others claim a long AT_MA reduces E_run. For a given ankle joint moment, a short AT_MA permits a higher tendon strain energy storage, whereas a long AT_MA reduces muscle fascicle force and muscle energy cost but shortening velocity is increased, elevating the metabolic cost. These are all conflicting mechanisms to reduce E_run, since AT energy storage comes at a metabolic cost. Neither of these proposed mechanisms have been examined together.

Methods

We measured AT_MA using the tendon travel method in 17 males and 3 females (24 ± 3 years, 75 ± 11 kg, 177 ± 7 cm). They ran on a motorized treadmill for 10 min at 2.5 m · s^-1 while E_run was measured. AT strain energy storage, muscle lengths, velocities and muscle energy cost were calculated during time-normalized stance from force and ultrasound data. A short (SHORT n = 11, AT_MA = 29.5 ± 2.0 mm) and long (LONG, n = 9, AT_MA = 36.6 ± 2.5 mm) AT_MA group was considered based on a bimodal distribution of measured AT_MA.

Results

Mean E_run was 4.9 ± 0.4 J · kg^-1 · m^-1. The relationship between AT_MA and E_run was not significant (r ² = 0.13, p = 0.12). Maximum AT force during stance was significantly lower in LONG (5,819 ± 1,202 N) compared to SHORT (6,990 ± 920 N, p = 0.028). Neither AT stretch nor AT strain energy storage was different between groups (mean difference: 0.3 ± 1 J · step^-1, p = 0.84). Fascicle force was significantly higher in SHORT (508 ± 93 N) compared to LONG (468 ± 84 N. p = 0.02). Fascicle lengths and velocities were similar between groups (p > 0.72). Muscle energy cost was significantly lower in LONG (0.028 ± 0.08 J · kg · step^-1) compared to SHORT (0.045 ± 0.14 J · kg · step^-1 p = 0.004). There was a significant negative relationship between AT_MA and total muscle energy cost relative to body mass across the stance phase (r = -0.699, p < 0.001).

Discussion

Together these results suggest that a LONG AT_MA serves to potentially reduce E_run by reducing the muscle energy cost of the plantarflexors during stance. The relative importance of AT energy storage and return in reducing E_run should be re-considered.

Muscle and tendon adaptations to moderate load eccentric vs. concentric resistance exercise in young and older males

Resistance exercise training (RET) is well-known to counteract negative age-related changes in both muscle and tendon tissue. Traditional RET consists of both concentric (CON) and eccentric (ECC) contractions; nevertheless, isolated ECC contractions are metabolically less demanding and, thus, may be more suitable for older populations. However, whether submaximal (60% 1RM) CON or ECC contractions differ in their effectiveness is relatively unknown. Further, whether the time course of muscle and tendon adaptations differs to the above is also unknown. Therefore, this study aimed to establish the time course of muscle and tendon adaptations to submaximal CON and ECC RET. Twenty healthy young (24.5 ± 5.1 years) and 17 older males (68.1 ± 2.4 years) were randomly allocated to either isolated CON or ECC RET which took place 3/week for 8 weeks. Tendon biomechanical properties, muscle architecture and maximal voluntary contraction were assessed every 2 weeks and quadriceps muscle volume every 4 weeks. Positive changes in tendon Young's modulus were observed after 4 weeks in all groups after which adaptations in young males plateaued but continued to increase in older males, suggesting a dampened rate of adaptation with age. However, both CON and ECC resulted in similar overall changes in tendon Young's modulus, in all groups. Muscle hypertrophy and strength increases were similar between CON and ECC in all groups. However, pennation angle increases were greater in CON, and fascicle length changes were greater in ECC. Notably, muscle and tendon adaptations appeared to occur in synergy, presumably to maintain the efficacy of the muscle-tendon unit.

Parameter estimation and experimental design for Hill-type muscles: Impulses from optimization-based modeling

The benefits of optimization-based modeling for parameter estimation of Hill-type muscle models are demonstrated. Therefore, we examined the model and data of Günther et al. (2007), who analyzed isometric, concentric, and quick-release contractions of a piglet calf muscle. We found that the isometric experiments are suitable for derivative-based parameter estimation while the others did not provide any additional value. During the estimation process, certain parameters had to be fixed. We give possible reasons and provide impulses for modelers. Subsequently, unnecessarily complex or deprecated model parts were exchanged and the new model was fitted to the data. In order to be able to provide a reliable estimation of the whole parameter set, we propose two isometric and two quick-release experiments, which are real-life feasible and together allow an identification of all parameters based on a local sensitivity analysis. These experiments can be used as qualitative guidelines for practitioners to reduce the experimental effort when estimating parameters for macroscopic muscle models.

Muscle activation during leg-press exercise with or without eccentric overload

PURPOSE

This study sought to investigate the electromyographic activity of the vastus lateralis (VL) muscle during concentric-eccentric exercise using a new concept leg press machine enabling a preset overloading in the eccentric phase.

METHODS

Ten young males familiar with resistive exercise were recruited for this study. Tests were performed on a Leg-press Biostrength® (Technogym S.p.A., Italy). The load was set to 70% and 80% of one-repetition maximum (1-RM). The participants performed 2 sets of 6 repetitions at each relative load with (ECC +) and without (ISOW) an eccentric overload equivalent to 150% of the concentric load. A metronome was employed to maintain the selected cadence. Sets were separated by a 5-min rest. Surface electromyography (EMG) of VL was recorded and integrated (iEMG).

RESULTS

Results showed a higher iEMG in ECC + with respect to ISOW at both intensities (+ 29% for 70% 1-RM, p < 0.01 and + 31% for 80% 1-RM, p < 0.001). No statistically significant differences were detected between concentric and eccentric phase in both ECC + conditions.

CONCLUSIONS

Training with a 150% eccentric overload provides a ~ 30% greater motor unit recruitment of the VL muscle in leg press exercise. Moreover, the results show that the eccentric overloading provided by the Biostrength® machine enables training at the same level of neural activation of the concentric phase. Hence, the derecruitment of motor units, normally observed during the eccentric phase when using conventional training machines, was overcome using the Biostrength® machine; this observation seems particularly important for maximizing neuromuscular responses to strength training.

Effects of a titin mutation on force enhancement and force depression in mouse soleus muscles

The active isometric force produced by muscles varies with muscle length in accordance with the force-length relationship. Compared with isometric contractions at the same final length, force increases after active lengthening (force enhancement) and decreases after active shortening (force depression). In addition to cross-bridges, titin has been suggested to contribute to force enhancement and depression. Although titin is too compliant in passive muscles to contribute to active tension at short sarcomere lengths on the ascending limb and plateau of the force-length relationship, recent evidence suggests that activation increases titin stiffness. To test the hypothesis that titin plays a role in force enhancement and depression, we investigated isovelocity stretching and shortening in active and passive wild-type and mdm (muscular dystrophy with myositis) soleus muscles. Skeletal muscles from mdm mice have a small deletion in the N2A region of titin and show no increase in titin stiffness during active stretch. We found that: (1) force enhancement and depression were reduced in mdm soleus compared with wild-type muscles relative to passive force after stretch or shortening to the same final length; (2) force enhancement and force depression increased with amplitude of stretch across all activation levels in wild-type muscles; and (3) maximum shortening velocity of wild-type and mdm muscles estimated from isovelocity experiments was similar, although active stress was reduced in mdm compared with wild-type muscles. The results of this study suggest a role for titin in force enhancement and depression, which contribute importantly to muscle force during natural movements.

Calcium-dependent titin-thin filament interactions in muscle: observations and theory

Gaps in our understanding of muscle mechanics demonstrate that the current model is incomplete. Increasingly, it appears that a role for titin in active muscle contraction might help to fill these gaps. While such a role for titin is increasingly accepted, the underlying molecular mechanisms remain unclear. The goals of this paper are to review recent studies demonstrating Ca²⁺-dependent interactions between N2A titin and actin in vitro, to explore theoretical predictions of muscle behavior based on this interaction, and to review experimental data related to the predictions. In a recent study, we demonstrated that Ca²⁺ increases the association constant between N2A titin and F-actin; that Ca²⁺ increases rupture forces between N2A titin and F-actin; and that Ca²⁺ and N2A titin reduce sliding velocity of F-actin and reconstituted thin filaments in motility assays. Preliminary data support a role for Ig83, but other Ig domains in the N2A region may also be involved. Two mechanical consequences are inescapable if N2A titin binds to thin filaments in active muscle sarcomeres: (1) the length of titin's freely extensible I-band should decrease upon muscle activation; and (2) binding between N2A titin and thin filaments should increase titin stiffness in active muscle. Experimental observations demonstrate that these properties characterize wild type muscles, but not muscles from mdm mice with a small deletion in N2A titin, including part of Ig83. Given the new in vitro evidence for Ca²⁺-dependent binding between N2A titin and actin, it is time for skepticism to give way to further investigation.

The slack test does not assess maximal shortening velocity of muscle fascicles in humans

The application of a series of extremely high accelerative motor-driven quick releases while muscles contract isometrically (i.e. slack test) has been proposed to assess unloaded velocity in human muscle. This study aimed to measure gastrocnemius medialis fascicle shortening velocity (V_F) and tendinous tissue shortening velocity during motor-driven quick releases performed at various activation levels to assess the applicability of the slack test in humans. Gastrocnemius medialis peak V_F and joint velocity recorded from 25 participants using high frame rate ultrasound during quick releases (at activation levels from 0% to 60% of maximal voluntary isometric torque) and during fast contractions without external load (ballistic condition) were compared. Unloaded joint velocity calculated using the slack test method increased whereas V_F decreased with muscle activation level (P≤0.03). Passive and low-level quick releases elicited higher V_F values (≥41.8±10.7 cm s^-1) compared with the ballistic condition (36.3±8.7 cm s^-1), while quick releases applied at 60% of maximal voluntary isometric torque produced the lowest V_F These findings suggest that initial fascicle length, complex fascicle-tendon interactions, unloading reflex and motor-driven movement pattern strongly influence and limit the shortening velocity achieved during the slack test. Furthermore, V_F elicited by quick releases is likely to reflect substantial contributions of passive processes. Therefore, the slack test is not appropriate to assess maximal muscle shortening velocity in vivo.

Improved Hill-type musculotendon models with activation-force-length coupling

BACKGROUND

Hill-type musculotendon models are most commonly used in biomechanical simulations for their computational efficacy and efficiency. But these models are generally built for maximally-activated muscles and linearly scale muscle properties when applied to submaximal conditions. However, the precondition of this scaling, which is muscle activation and properties are independent each other, has been proven unreal in many studies. Actually, the maximal activation condition is not ubiquitous for muscles in vivo, so it is necessary to adapt the linear scaling approach to improve the model practicability.

OBJECTIVE

This paper aimed at proposing two improved Hill-type musculotendon models that are better suited for submaximal conditions.

METHOD

These two models were built by including the activation-force-length coupling and their biological accuracy and computation speed were evaluated by a series of benchmark simulations.

RESULTS

Compared to experimental measurements, the percent root mean square errors of forces calculated by the two AFLC models were less than 13.98% and 13.81% respectively. However, the average running time of the second AFLC model was nearly 17 times that of the first one with only a little improvement in accuracy.

CONCLUSION

The two AFLC models were validated more accurate than the common Hill-type model in submaximally activated conditions and the first one was recommended in the construction of upper-layer musculoskeletal models.

Theoretical considerations for muscle-energy savings during distance running

We have recently demonstrated that the triceps surae muscles energy cost (EC_TS) represents a substantial portion of the total metabolic cost of running (E_run). Therefore, it seems relevant to evaluate the factors which dictate EC_TS, namely the amount and velocity of shortening, since it is likely these factors will dictate E_run. E_run and triceps surae morphological and AT mechanical properties were obtained in 46 trained and elite male and female distance runners using ultrasonography and dynamometry. EC_TS (J·stride^-1) at the speed of lactate threshold (sLT) was estimated from AT force and crossbridge mechanics and energetics. To estimate the relative impact of these factors on EC_TS, mean values for running speed, body mass, resting fascicle length (L_f), Achilles tendon stiffness and moment arm and maximum isometric plantarflexion torque were obtained. EC_TS was calculated across a range (mean ± 1 sd) of values for each independent factor. Average sLT was 233 m·min^-1. At this speed, EC_TS was 255 J·stride^-1. Estimated fascicle shortening velocity was 0.08 V_max and the level of muscle activation was 84.7% of maximum isometric torque. Compared to the EC_TS calculated from the lowest range of values obtained for each independent factor, higher AT stiffness was associated with a 39% reduction in EC_TS, 81% reduction in fascicle shortening velocity and a 31% reduction in muscle activation. Longer AT moment arms and elevated body masses were associated with an increase in EC_TS of 18% and 23%, respectively. These results demonstrate that a low EC_TS is achieved primarily from a high AT stiffness and low body mass, which is exemplified in elite distance runners.

Running Economy from a Muscle Energetics Perspective

The economy of running has traditionally been quantified from the mass-specific oxygen uptake; however, because fuel substrate usage varies with exercise intensity, it is more accurate to express running economy in units of metabolic energy. Fundamentally, the understanding of the major factors that influence the energy cost of running (Erun) can be obtained with this approach. Erun is determined by the energy needed for skeletal muscle contraction. Here, we approach the study of Erun from that perspective. The amount of energy needed for skeletal muscle contraction is dependent on the force, duration, shortening, shortening velocity, and length of the muscle. These factors therefore dictate the energy cost of running. It is understood that some determinants of the energy cost of running are not trainable: environmental factors, surface characteristics, and certain anthropometric features. Other factors affecting Erun are altered by training: other anthropometric features, muscle and tendon properties, and running mechanics. Here, the key features that dictate the energy cost during distance running are reviewed in the context of skeletal muscle energetics.

Early structural remodeling and deuterium oxide-derived protein metabolic responses to eccentric and concentric loading in human skeletal muscle

We recently reported that the greatest distinguishing feature between eccentric (ECC) and concentric (CON) muscle loading lays in architectural adaptations: ECC favors increases in fascicle length (Lf), associated with distal vastus lateralis muscle (VL) hypertrophy, and CON increases in pennation angle (PA). Here, we explored the interactions between structural and morphological remodeling, assessed by ultrasound and dual x-ray absorptiometry (DXA), and long-term muscle protein synthesis (MPS), evaluated by deuterium oxide (D₂O) tracing technique. Ten young males (23 ± 4 years) performed unilateral resistance exercise training (RET) three times/week for 4 weeks; thus, one-leg trained concentrically while the contralateral performed ECC exercise only at 80% of either CON or ECC one repetition maximum (1RM). Subjects consumed an initial bolus of D₂O (150 mL), while a 25-mL dose was thereafter provided every 8 days. Muscle biopsies from VL midbelly (MID) and distal myotendinous junction (MTJ) were collected at 0 and 4-weeks. MPS was then quantified via GC–pyrolysis–IRMS over the 4-week training period. Expectedly, ECC and CON RET resulted in similar increases in VL muscle thickness (MT) (7.5% vs. 8.4%, respectively) and thigh lean mass (DXA) (2.3% vs. 3%, respectively), albeit through distinct remodeling: Lf increasing more after ECC (5%) versus CON (2%) and PA increasing after CON (7% vs. 3%). MPS did not differ between contractile modes or biopsy sites (MID-ECC: 1.42 vs. MID-CON: 1.4% day⁻¹; MTJ-ECC: 1.38 vs. MTJ-CON: 1.39% day⁻¹). Muscle thickness at MID site increased similarly following ECC and CON RET, reflecting a tendency for a contractile mode-independent correlation between MPS and MT (P = 0.07; R² = 0.18). We conclude that, unlike MT, distinct structural remodeling responses to ECC or CON are not reflected in MPS; the molecular mechanisms of distinct protein deposition, and/or the role of protein breakdown in mediating these responses remain to be defined.

Spinal circuits can accommodate interaction torques during multijoint limb movements

The dynamic interaction of limb segments during movements that involve multiple joints creates torques in one joint due to motion about another. Evidence shows that such interaction torques are taken into account during the planning or control of movement in humans. Two alternative hypotheses could explain the compensation of these dynamic torques. One involves the use of internal models to centrally compute predicted interaction torques and their explicit compensation through anticipatory adjustment of descending motor commands. The alternative, based on the equilibrium-point hypothesis, claims that descending signals can be simple and related to the desired movement kinematics only, while spinal feedback mechanisms are responsible for the appropriate creation and coordination of dynamic muscle forces. Partial supporting evidence exists in each case. However, until now no model has explicitly shown, in the case of the second hypothesis, whether peripheral feedback is really sufficient on its own for coordinating the motion of several joints while at the same time accommodating intersegmental interaction torques. Here we propose a minimal computational model to examine this question. Using a biomechanics simulation of a two-joint arm controlled by spinal neural circuitry, we show for the first time that it is indeed possible for the neuromusculoskeletal system to transform simple descending control signals into muscle activation patterns that accommodate interaction forces depending on their direction and magnitude. This is achieved without the aid of any central predictive signal. Even though the model makes various simplifications and abstractions compared to the complexities involved in the control of human arm movements, the finding lends plausibility to the hypothesis that some multijoint movements can in principle be controlled even in the absence of internal models of intersegmental dynamics or learned compensatory motor signals.

Task-dependent control of the jaw during food splitting in humans

Although splitting of food items between the incisors often requires high bite forces, rarely do the teeth harmfully collide when the jaw quickly closes after split. Previous studies indicate that the force-velocity relationship of the jaw closing muscles principally explains the prompt dissipation of jaw closing force. Here, we asked whether people could regulate the dissipation of jaw closing force during food splitting. We hypothesized that such regulation might be implemented via differential recruitment of masseter muscle portions situated along the anteroposterior axis because these portions will experience a different shortening velocity during jaw closure. Study participants performed two different tasks when holding a peanut-half stacked on a chocolate piece between their incisors. In one task, they were asked to split the peanut-half only (single-split trials) and, in the other, to split both the peanut and the chocolate in one action (double-split trials). In double-split trials following the peanut split, the intensity of the tooth impact on the chocolate piece was on average 2.5 times greater than in single-split trials, indicating a substantially greater loss of jaw closing force in the single-split trials. We conclude that control of jaw closing force dissipation following food splitting depends on task demands. Consistent with our hypothesis, converging neurophysiological and morphometric data indicated that this control involved a differential activation of the jaw closing masseter muscle along the anteroposterior axis. These latter findings suggest that the regulation of jaw closing force after sudden unloading of the jaw exploits masseter muscle compartmentalization.

Interdependence of torque, joint angle, angular velocity and muscle action during human multi-joint leg extension

PURPOSE

Force and torque production of human muscles depends upon their lengths and contraction velocity. However, these factors are widely assumed to be independent of each other and the few studies that dealt with interactions of torque, angle and angular velocity are based on isolated single-joint movements. Thus, the purpose of this study was to determine force/torque-angle and force/torque-angular velocity properties for multi-joint leg extensions.

METHODS

Human leg extension was investigated (n = 18) on a motor-driven leg press dynamometer while measuring external reaction forces at the feet. Extensor torque in the knee joint was calculated using inverse dynamics. Isometric contractions were performed at eight joint angle configurations of the lower limb corresponding to increments of 10° at the knee from 30 to 100° of knee flexion. Concentric and eccentric contractions were performed over the same range of motion at mean angular velocities of the knee from 30 to 240° s(-1).

RESULTS

For contractions of increasing velocity, optimum knee angle shifted from 52 ± 7 to 64 ± 4° knee flexion. Furthermore, the curvature of the concentric force/torque-angular velocity relations varied with joint angles and maximum angular velocities increased from 866 ± 79 to 1,238 ± 132° s(-1) for 90-50° knee flexion. Normalised eccentric forces/torques ranged from 0.85 ± 0.12 to 1.32 ± 0.16 of their isometric reference, only showing significant increases above isometric and an effect of angular velocity for joint angles greater than optimum knee angle.

CONCLUSIONS

The findings reveal that force/torque production during multi-joint leg extension depends on the combined effects of angle and angular velocity. This finding should be accounted for in modelling and optimisation of human movement.

Architectural, functional and molecular responses to concentric and eccentric loading in human skeletal muscle

AIM

We investigated architectural, functional and molecular responses of human skeletal muscle to concentric (CON) or eccentric (ECC) resistance training (RT).

METHODS

Twelve young males performed 10 weeks of concentric (CON) or eccentric (ECC) resistance training (RT) (n = 6 CON, 6 ECC). An additional 14 males were recruited to evaluate acute muscle fascicle behaviour and molecular signalling in biopsies collected from vastus lateralis (VL) after 30 min of single bouts of CON or ECC exercise. VL volume was measured by magnetic resonance imaging. Muscle architecture (fascicle length, Lf; pennation angle, PA) was evaluated by ultrasonography. Muscle remodelling signals to CON or ECC loading [MAPK/AKT-mammalian target of rapamycin (mTOR) signalling] and inflammatory pathway (TNFαMurf-1-MAFbx) were evaluated by immunoblotting.

RESULTS

Despite the ~1.2-fold greater load of the ECC group, similar increases in muscle volume (+8% CON and +6% ECC) and in maximal voluntary isometric contraction (+9% CON and +11% ECC) were found after RT. However, increases in Lf were greater after ECC than CON (+12 vs. +5%) while increases in PA were greater in CON than ECC (+30 vs. +5%). Distinct architectural adaptations were associated with preferential growth in the distal regions of VL for ECC (+ECC +8% vs. +CON +2) and mid belly for CON (ECC +7 vs. CON +11%). While MAPK activation (p38MAPK, ERK1/2, p90RSK) was specific to ECC, neither mode affected AKT-mTOR or inflammatory signalling 30 min after exercise.

CONCLUSION

Muscle growth with CON and ECC RT occurs with different morphological adaptations reflecting distinct fibre fascicle behaviour and molecular responses.

A joint computational respiratory neural network-biomechanical model for breathing and airway defensive behaviors

Data-driven computational neural network models have been used to study mechanisms for generating the motor patterns for breathing and breathing related behaviors such as coughing. These models have commonly been evaluated in open loop conditions or with feedback of lung volume simply represented as a filtered version of phrenic motor output. Limitations of these approaches preclude assessment of the influence of mechanical properties of the musculoskeletal system and motivated development of a biomechanical model of the respiratory muscles, airway, and lungs using published measures from human subjects. Here we describe the model and some aspects of its behavior when linked to a computational brainstem respiratory network model for breathing and airway defensive behavior composed of discrete "integrate and fire" populations. The network incorporated multiple circuit paths and operations for tuning inspiratory drive suggested by prior work. Results from neuromechanical system simulations included generation of a eupneic-like breathing pattern and the observation that increased respiratory drive and operating volume result in higher peak flow rates during cough, even when the expiratory drive is unchanged, or when the expiratory abdominal pressure is unchanged. Sequential elimination of the model's sources of inspiratory drive during cough also suggested a role for disinhibitory regulation via tonic expiratory neurons, a result that was subsequently supported by an analysis of in vivo data. Comparisons with antecedent models, discrepancies with experimental results, and some model limitations are noted.

Force control of quadriceps muscle is bilaterally impaired in subacute stroke

We tested the hypothesis that force variability and error during maintenance of submaximal isometric knee extension are greater in subacute stroke patients than in controls and are related to motor impairments. Contralesional (more-affected) and ipsilesional (less-affected) legs of 33 stroke patients with sufficiently high motor abilities (62 ± 13 yr, 16 ± 2 days postinjury) and the dominant leg of 20 controls (62 ± 10 yr) were tested in sitting position. After peak knee extension torque [maximum voluntary contraction (MVC)] was established, subjects maintained 10, 20, 30, and 50% of MVC as steady and accurate as possible for 10 s by matching voluntary force to the target level displayed on a monitor. Coefficient of variation (CV) and root-mean-square error (RMSE) were used to quantify force variability and error, respectively. The MVC was significantly smaller in the more-affected than less-affected leg, and both were significantly lower than in controls. The CV was significantly larger in the more-affected than less-affected leg at 20 and 50% MVC, whereas both were significantly larger compared with controls across all force levels. Both more-affected and less-affected legs of patients showed significantly greater RMSE than controls at 30 and 50% MVC. The CV and RMSE were not related to the Fugl-Meyer motor score or to the Rivermead Mobility Index. The CV negatively correlated with MVC in controls but only in the less-affected leg of patients. It is concluded that isometric knee extension strength and force control are bilaterally impaired soon after stroke but more so in the more-affected leg. Future studies should examine possible mechanisms and the evolution of these changes.

Unloaded shortening velocity of voluntarily and electrically activated human dorsiflexor muscles in vivo

We have previously shown that unloaded shortening velocity (V(0)) of human plantar flexors can be determined in vivo, by applying the "slack test" to submaximal voluntary contractions (J Physiol 567:1047-1056, 2005). In the present study, to investigate the effect of motor unit recruitment pattern on V(0) of human muscle, we modified the slack test and applied this method to both voluntary and electrically elicited contractions of dorsiflexors. A series of quick releases (i.e., rapid ankle joint rotation driven by an electrical dynamometer) was applied to voluntarily activated dorsiflexor muscles at three different contraction intensities (15, 50, and 85% of maximal voluntary contraction; MVC). The quick-release trials were also performed on electrically activated dorsiflexor muscles, in which three stimulus conditions were used: submaximal (equal to 15%MVC) 50-Hz stimulation, supramaximal 50-Hz stimulation, and supramaximal 20-Hz stimulation. Modification of the slack test in vivo resulted in good reproducibility of V(0), with an intraclass correlation coefficient of 0.87 (95% confidence interval: 0.68-0.95). Regression analysis showed that V(0) of voluntarily activated dorsiflexor muscles significantly increased with increasing contraction intensity (R(2) = 0.52, P<0.001). By contrast, V(0) of electrically activated dorsiflexor muscles remained unchanged (R(2)<0.001, P = 0.98) among three different stimulus conditions showing a large variation of tetanic torque. These results suggest that the recruitment pattern of motor units, which is quite different between voluntary and electrically elicited contractions, plays an important role in determining shortening velocity of human skeletal muscle in vivo.

Changes in tendon stiffness and running economy in highly trained distance runners

The purpose of this study was to determine if changes in triceps-surae tendon stiffness (TST K) could affect running economy (RE) in highly trained distance runners. The intent was to induce increased TST K in a subgroup of runners by an added isometric training program. If TST K is a primary determinant of RE, then the energy cost of running (EC) should decrease in the trained subjects. EC was measured via open-circuit spirometry in 12 highly trained male distance runners, and TST K was measured using ultrasonography and dynamometry. Runners were randomly assigned to either a training or control group. The training group performed 4 × 20 s isometric contractions at 80% of maximum voluntary plantarflexion moment three times per week for 8 weeks. All subjects (V(O)₂(max)) = 67.4 ± 4.6 ml kg(-1) min(-1)) continued their usual training for running. TST K was measured every 2 weeks. EC was measured in both training and control groups before and after the 8 weeks at three submaximal velocities, corresponding to 75, 85 and 95% of the speed at lactate threshold (sLT). Isometric training did neither result in a mean increase in TST K (0.9 ± 25.8%) nor a mean improvement in RE (0.1 ± 3.6%); however, there was a significant relationship (r(2) = 0.43, p = 0.02) between the change in TST K and change in EC, regardless of the assigned group. It was concluded that TST K and EC are somewhat labile and change together.

High-frequency oscillations as a consequence of neglected serial damping in Hill-type muscle models

High-frequency vibrations e.g., induced by legs impacting with the ground during terrestrial locomotion can provoke damage within tendons even leading to ruptures. So far, macroscopic Hill-type muscle models do not account for the observed high-frequency damping at low-amplitudes. Therefore, former studies proposed that protective damping might be explained by modelling the contractile machinery of the muscles in more detail, i.e., taking the microscopic processes of the actin-myosin coupling into account. In contrast, this study formulates an alternative hypothesis: low but significant damping of the passive material in series to the contractile machinery--e.g., tendons, aponeuroses, titin--may well suffice to damp these hazardous vibrations. Thereto, we measured the contraction dynamics of a piglet muscle-tendon complex (MTC) in three contraction modes at varying loads and muscle-tendon lengths. We simulated all three respective load situations on a computer: a Hill-type muscle model including a contractile element (CE) and each an elastic element in parallel (PEE) and in series (SEE) to the CE pulled on a loading mass. By comparing the model to the measured output of the MTC, we extracted a consistent set of muscle parameters. We varied the model by introducing either linear damping in parallel or in series to the CE leading to accordant re-formulations of the contraction dynamics of the CE. The comparison of the three cases (no additional damping, parallel damping, serial damping) revealed that serial damping at a physiological magnitude suffices to explain damping of high-frequency vibrations of low amplitudes. The simulation demonstrates that any undamped serial structure within the MTC enforces SEE-load eigenoscillations. Consequently, damping must be spread all over the MTC, i.e., rather has to be de-localised than localised within just the active muscle material. Additionally, due to suppressed eigenoscillations Hill-type muscle models taking into account serial damping are numerically more efficient when used in macroscopic biomechanical neuro-musculo-skeletal models.

An improved method to determine neuromuscular properties using force laws - From single muscle to applications in human movements

We evaluate an improved method for individually determining neuromuscular properties in vivo. The method is based on Hill's equation used as a force law combined with Newton's equation of motion. To ensure the range of validity of Hill's equation, we first perform detailed investigations on in vitro single muscles. The force-velocity relation determined with the model coincides well with results obtained by standard methods (r=.99) above 20% of the isometric force. In addition, the model-predicted force curves during work loop contractions very well agree with measurements (mean difference: 2-3%). Subsequently, we deduce theoretically under which conditions it is possible to combine several muscles of the human body to model muscles. This leads to a model equation for human leg extension movements containing parameters for the muscle properties and for the activation. To numerically determine these invariant neuromuscular properties we devise an experimental method based on concentric and isometric leg extensions. With this method we determine individual muscle parameters from experiments such that the simulated curves agree well with experiments (r=.99). A reliability test with 12 participants revealed correlations r=.72-.91 for the neuromuscular parameters (p<.01). Predictions of similar movements under different conditions show mean errors of about 5%. In addition, we present applications in sports practise and theory.

Steady-state force–velocity relation in human multi-joint movement determined with force clamp analysis

To study the force-velocity characteristics of human knee-hip extension movement, a dynamometer, in which force was controlled by a servo system, was developed. Seated subjects pressed either bilaterally or unilaterally a force plate, a horizontal position of which was servo-controlled so as to equalize the measured force and a force command generated by a computer at a time resolution of 2 ms (force clamp). The force command was based on the relation between maximum isometric force and foot position within the range between 70% and 90% of "leg length" (LL: longitudinal distance between the sole of the foot and the hip joint), so that the same force relative to the maximum isometric force was consistently applied regardless of the foot position. By regulating the force according to this function, the force-velocity relation was determined. The force-velocity relation obtained was described by a linear function (n=17, r=-0.986 for 80% LL, r=-0.968 for 85% LL) within a range of force between 0.1 and 0.8F(0) (maximum isometric force). The maximum force extrapolated from the linear regression (F(max)) coincided with F(0) (n=17, F(0)/F(max)=1.00+/-0.09 for 80% LL and 1.00+/-0.20 for 85% LL). Also, the velocity at zero force (V(max)) was obtained from the extrapolation. When compared to the bilateral movements, unilateral movements gave rise to a smaller F(max) but the same V(max), suggesting that V(max) is independent of force and therefore represents the proper unloaded velocity. It is suggested that some neural mechanisms may be involved in the force-velocity relation of the knee-hip extension movement, and make it exhibit a linear appearance rather than a hyperbola.

Are the maximum shortening velocity and the shape parameter in a Hill-type model of whole muscle related to activation?

Mathematical models of the inter-relationship of muscle force, velocity, and activation are useful in forward dynamic simulations of human movement tasks. The objective of this work was to determine whether the parameters (maximum shortening velocity V(max) and shape parameter k) of a Hill-type muscle model, interrelating muscle force, velocity, and activation, are themselves dependent on the activation. To fulfill this objective, surface EMG signals from four muscles, as well as the kinematics and kinetics of the arm, were recorded from 14 subjects who performed rapid-release elbow extension tasks at 25%, 50%, 75%, and 100% activation (MVC). The experimental elbow flexion angle was tracked by a forward dynamic simulation of the task in which V(max) and k of the triceps brachii were varied at each activation level to minimize the difference between the simulated and experimental elbow flexion angle. Because a preliminary analysis demonstrated no dependency of k on activation, additional simulations were performed with constant k values of 0.15, 0.20, and 0.25. The optimized values of V(max) normalized to the average value within a subject were then regressed onto the activation. Normalized V(max) depended significantly on the activation (p<0.001) for all values of k. Furthermore, the estimated V(max) values were not sensitive to the selected k value. The results support the use of Hill-type models in which V(max) depends on activation in forward dynamic simulations modeling muscles with mixed fiber-type composition recruited in the range of 25-100% activation. The use of more accurate models will lend greater confidence to the results of forward dynamic simulations.

Shortening velocity of human triceps surae muscle measured with the slack test in vivo

Unloaded shortening velocity (V(0)) of human triceps surae muscle was measured in vivo by applying the 'slack test', originally developed for determining V(0) of single muscle fibres, to voluntary contractions at varied activation levels (ALs). V(0) was measured from 10 subjects at five different ALs defined as a fraction (5, 10, 20, 40 and 60%) of the maximum voluntary contraction (MVC) torque. Although individual variability was apparent, V(0) tended to increase with AL (R(2) = 0.089; P = 0.035) up to 60%MVC (8.6 +/- 2.6 rad s(-1)). This value of V(0) at 60%MVC was comparable to the maximum shortening velocity of plantar flexors reported in the previous studies. Electromyographic analysis showed that the activities of soleus, medial gastrocnemius and lateral gastrocnemius muscles increased with AL during isometric contraction and after the application of quick release in a similar manner. Also, it showed that the activity of an antagonist, tibialis anterior muscle, was negligible, even though a slight increase took place after the quick release of agonist. Correlation analysis showed that there were no significant correlations between V(0) and MVC torque normalized with respect to body mass, although the correlation coefficient was relatively high at low ALs. The results suggest that in human muscle, V(0) represents the unloaded velocity of the fastest muscle fibres recruited, and increases with AL possibly because of progressive recruitment of faster fibres. Individual variability may be explained, at least partially, by the difference in fibre-type composition.

ISOFIT: a model-based method to measure muscle–tendon properties simultaneously

Estimation of muscle parameters specifying force-length and force-velocity behavior requires in general a large number of sophisticated experiments often including a combination of isometric, isokinetic, isotonic, and quick-release experiments. This study validates a simpler method (ISOFIT) to determine muscle properties by fitting a Hill-type muscle model to a set of isovelocity data. Muscle properties resulting from the ISOFIT method agreed well with muscle properties determined separately in in vitro measurements using frog semitendinosus muscles. The force-length curve was described well by the results of the model. The force-velocity curve resulting from the model coincided with the experimentally determined curve above approximately 20% of maximum isometric force (correlation coefficient R>0.99). At lower forces and thus higher velocities the predicted curve underestimated velocity. The stiffness of the series elastic component determined with direct experiments was approximately 10% lower than that determined by the ISOFIT method. Use of the ISOFIT method can decrease experimental time up to 80% and reduce potential changes in muscle parameters due to fatigue.

Development and Validation of a 3-D Model to Predict Knee Joint Loading During Dynamic Movement

The purpose of this study was to develop a subject-specific 3-D model of the lower extremity to predict neuromuscular control effects on 3-D knee joint loading during movements that can potentially cause injury to the anterior cruciate ligament (ACL) in the knee. The simulation consisted of a forward dynamic 3-D musculoskeletal model of the lower extremity, scaled to represent a specific subject. Inputs of the model were the initial position and velocity of the skeletal elements, and the muscle stimulation patterns. Outputs of the model were movement and ground reaction forces, as well as resultant 3-D forces and moments acting across the knee joint. An optimization method was established to find muscle stimulation patterns that best reproduced the subject’s movement and ground reaction forces during a sidestepping task. The optimized model produced movements and forces that were generally within one standard deviation of the measured subject data. Resultant knee joint loading variables extracted from the optimized model were comparable to those reported in the literature. The ability of the model to successfully predict the subject’s response to altered initial conditions was quantified and found acceptable for use of the model to investigate the effect of altered neuromuscular control on knee joint loading during sidestepping. Monte Carlo simulations (N=100,000) using randomly perturbed initial kinematic conditions, based on the subject’s variability, resulted in peak anterior force, valgus torque and internal torque values of 378 N, 94 Nm and 71 Nm, respectively, large enough to cause ACL rupture. We conclude that the procedures described in this paper were successful in creating valid simulations of normal movement, and in simulating injuries that are caused by perturbed neuromuscular control.

A computer-graphics model of muscle activation and contraction dynamics

An interactive computer-graphics model of the mechanical interaction between the structural components of a Hill-type muscle model is presented. The model allows the length, velocity, and force in the contractile component, parallel elastic component, and series elastic of a generic (normalised) muscle-tendon complex to be computed from input defining instantaneous muscle stimulation and muscle-tendon length. By altering model inputs and model parameter values the user can observe corresponding model behaviour in schematic and graphical form. In this paper a general description of the model is provided together with examples of how the model can be used to simulate different contractile conditions. It is intended that the model will serve as a useful tool for demonstrating the behaviour of Hill-type muscle models.

Diaphragm contractile dysfunction in MyoD gene-inactivated mice

MyoD is one of four myogenic regulatory factors found exclusively in skeletal muscle. In an effort to better understand the role that MyoD plays in determining muscle contractile properties, we examined the effects of MyoD deletion on both diaphragmatic contractile properties and myosin heavy chain (MHC) phenotype. Regions of the costal diaphragm from wild-type and MyoD knockout [MyoD (-/-)] adult male BALB/c mice (n = 8/group) were removed, and in vitro diaphragmatic contractile properties were measured. Diaphragmatic contractile measurements revealed that MyoD (-/-) animals exhibited a significant (P < 0.05) downward shift in the force-frequency relationship, a decrement in maximal specific tension (P(o); -33%), a decline in maximal shortening velocity (V(max); -37%), and concomitant decrease in peak power output (-47%). Determination of MHC isoforms in the diaphragm via gel electrophoresis revealed that MyoD elimination resulted in a fast-to-slow shift (P < 0.05) in the MHC phenotype toward MHC types IIA and IIX in MyoD (-/-) animals. These data indicate that MyoD deletion results in a decrease in diaphragmatic submaximal force generation and P(o), along with decrements in both V(max) and peak power output. Hence, MyoD plays an important role in determining diaphragmatic contractile properties.

Cadence, power, and muscle activation in cycle ergometry

PURPOSE

Based on the resistance-rpm relationship for cycling, which is not unlike the force-velocity relationship of muscle, it is hypothesized that the cadence which requires the minimal muscle activation will be progressively higher as power output increases.

METHODS

To test this hypothesis, subjects were instrumented with surface electrodes placed over seven muscles that were considered to be important during cycling. Measurements were made while subjects cycled at 100, 200, 300, and 400 W at each cadence: 50, 60, 80, 100, and 120 rpm. These power outputs represented effort which was up to 32% of peak power output for these subjects.

RESULTS

When all seven muscles were averaged together, there was a proportional increase in EMG amplitude each cadence as power increased. A second-order polynomial equation fit the EMG:cadence results very well (r2 = 0.87- 0.996) for each power output. Optimal cadence (cadence with lowest amplitude of EMG for a given power output) increased with increases in power output: 57 +/- 3.1, 70 +/- 3.7, 86 +/- 7.6, and 99 +/- 4.0 rpm for 100, 200, 300, and 400 W, respectively.

CONCLUSION

The results confirm that the level of muscle activation varies with cadence at a given power output. The minimum EMG amplitude occurs at a progressively higher cadence as power output increases. These results have implications for the sense of effort and preferential use of higher cadences as power output is increased.