Functional muscle synergies to support the knee against moment specific loads while weight bearing

ObjectiveExternally applied abduction and rotational loads are major contributors to the knee joint injury mechanism; yet, how muscles work together to stabilize the knee against these loads remains unclear. Our study sought to evaluate lower limb functional muscle synergies in healthy young adults such that muscle activation can be directly related to internal knee joint moments.MethodsConcatenated non-negative matrix factorization extracted muscle and moment synergies of 22 participants from electromyographic signals and joint moments elicited during a weight-bearing force matching protocol.ResultsTwo synergy sets were extracted: Set 1 included four synergies, each corresponding to a general anterior, posterior, medial, or lateral force direction. Frontal and transverse moments were coupled during medial and lateral force directions. Set 2 included six synergies, each corresponding to a moment type (extension/flexion, ab/adduction, internal/external rotation). Hamstrings and quadriceps dominated synergies associated with respective flexion and extension moments while quadriceps-hamstring co-activation was associated with knee abduction. Rotation moments were associated with notable contributions from hamstrings, quadriceps, gastrocnemius, and hip ab/adductors, corresponding to a general co-activation muscle synergy.ConclusionOur results highlight the importance of muscular co-activation of all muscles crossing the knee to support it during injury-inducing loading conditions such as externally applied knee abduction and rotation. Functional muscle synergies can provide new insight into the relationship between neuromuscular control and knee joint stability by directly associating biomechanical variables to muscle activation.     

Dopaminergic modulation of multi-muscle synergies in postural tasks performed by patients with Parkinson’s disease

BackgroundPostural instability is one of most disabling motor symptoms in Parkinson’s disease. Indices of multi-muscle synergies are new measurements of postural stability.ObjectivesWe explored the effects of dopamine-replacement drugs on multi-muscle synergies stabilizing center of pressure coordinate and their adjustments prior to a self-triggered perturbation in patients with Parkinson’s disease. We hypothesized that both synergy indices and synergy adjustments would be improved on dopaminergic drugs.MethodsPatients at Hoehn-Yahr stages II and III performed whole-body tasks both off- and on-drugs while standing. Muscle modes were identified as factors in the muscle activation space. Synergy indices stabilizing center of pressure in the anterior-posterior direction were quantified in the muscle mode space during a load-release task.ResultsDopamine-replacement drugs led to more consistent organization of muscles in stable groups (muscle modes). On-drugs patients showed larger indices of synergies and anticipatory synergy adjustments. In contrast, no medication effects were seen on anticipatory postural adjustments or other performance indices.ConclusionsDopamine-replacement drugs lead to significant changes in characteristics of multi-muscle synergies in Parkinson’s disease. Studies of synergies may provide a biomarker sensitive to problems with postural stability and agility and to efficacy of dopamine-replacement therapy.     

Subject-specific thumb muscle activity during functional tasks of daily life

BackgroundThe trapeziometacarpal joint is subjected to high compressive forces during powerful pinch and grasp tasks due to muscle loading. In addition, muscle contraction is important for stability of the joint. The aim of the present study is to explore if different muscle activation patterns can be found between three functional tasks.MethodsIsometric forces and fine-wire electromyographic (fEMG) activity produced by three intrinsic and four extrinsic thumb muscles were measured in 10 healthy female volunteers. The participants performed isometric contractions in a lateral key pinch, a power grasp and a jar twist task. The tasks were executed with and without EMG recording to verify if electrode placement influenced force production.ResultsA subject-specific muscle recruitment was found which remained largely unchanged across tasks. Extrinsic thumb muscles were significantly more active than intrinsic muscles in all tasks. Insertion of the fEMG electrodes decreased force production significantly in all tasks.ConclusionThe thumb muscles display a high variability in muscle activity during functional tasks of daily life. The results of this study suggest that to produce a substantial amount of force, a well-integrated, but subject-specific, co-contraction between the intrinsic and extrinsic thumb muscles is necessary.     

Recommendation for the minimum number of steps to analyze when performing the uncontrolled manifold analysis on walking data

The uncontrolled manifold (UCM) analysis quantifies the extent to which co-variation among a set of variables facilitates consistent performance by partitioning variance in those variables into two components then calculating their normalized difference (i.e., the synergy index). Although UCM-derived measures are thought to depend on the number of data points analyzed, the minimum number needed to reasonably approximate true values of these measures is unknown. For each of two performance variables related to mechanical stability of gait, we evaluated changes in UCM-derived measures when increasing the number of analyzed points, here steps. Fourteen older adults walked on a treadmill while motion capture tracked movement. For each subject, n steps (where n = 2–99) were randomly sampled from the first 100, then used to calculate UCM-derived variables. For each subject, variables were expressed as a percent of the subject-specific value with n = 100 and averaged across 50 simulations. For each n, 95% confidence intervals (CIs) were calculated from group data. The minimum number of steps to “reasonably approximate” a variables was defined as the value of n for which the lower CI was >90% of the value with n = 100. Regardless of performance variable, reasonable approximations of the synergy index were attained with n = 16 steps, whereas n = 50 steps were needed for each of the variance components However, the differences between using 16 steps and 50 steps were small. Collecting 15–20 steps is recommended for a reasonable approximation of the synergy indices considered herein, particularly when data collection is constrained to a limited number of steps.     

Muscle Synergies Heavily Influence the Neural Control of Arm Endpoint Stiffness and Energy Consumption

Much debate has arisen from research on muscle synergies with respect to both limb impedance control and energy consumption. Studies of limb impedance control in the context of reaching movements and postural tasks have produced divergent findings, and this study explores whether the use of synergies by the central nervous system (CNS) can resolve these findings and also provide insights on mechanisms of energy consumption. In this study, we phrase these debates at the conceptual level of interactions between neural degrees of freedom and tasks constraints. This allows us to examine the ability of experimentally-observed synergies—correlated muscle activations—to control both energy consumption and the stiffness component of limb endpoint impedance. In our nominal 6-muscle planar arm model, muscle synergies and the desired size, shape, and orientation of endpoint stiffness ellipses, are expressed as linear constraints that define the set of feasible muscle activation patterns. Quadratic programming allows us to predict whether and how energy consumption can be minimized throughout the workspace of the limb given those linear constraints. We show that the presence of synergies drastically decreases the ability of the CNS to vary the properties of the endpoint stiffness and can even preclude the ability to minimize energy. Furthermore, the capacity to minimize energy consumption—when available—can be greatly affected by arm posture. Our computational approach helps reconcile divergent findings and conclusions about task-specific regulation of endpoint stiffness and energy consumption in the context of synergies. But more generally, these results provide further evidence that the benefits and disadvantages of muscle synergies go hand-in-hand with the structure of feasible muscle activation patterns afforded by the mechanics of the limb and task constraints. These insights will help design experiments to elucidate the interplay between synergies and the mechanisms of learning, plasticity, versatility and pathology in neuromuscular systems.     

Muscle synergies are similar when typically developing children walk on a treadmill at different speeds and slopes

BackgroundThe aim of this study was to determine whether changes in synergies relate to changes in gait while walking on a treadmill at multiple speeds and slopes. The hypothesis was that significant changes in movement pattern would not be accompanied by significant changes in synergies, suggesting that synergies are not dependent on the mechanical constraints but are instead neurological in origin.MethodsSixteen typically developing children walked on a treadmill for nine combinations (stages) of different speeds and slopes while simultaneously collecting kinematics, kinetics, and surface electromyography (EMG) data. The kinematics for each stride were summarized using a modified version of the Gait Deviation Index that only includes the sagittal plane. The kinetics for each stride were summarized using a modified version of the Gait Deviation Index – Kinetic which includes sagittal plane moments and powers. Within each synergy group, the correlations of the synergies were calculated between the treadmill stages.ResultsWhile kinematics and kinetics were significantly altered at the highest slope compared to level ground when walking on a treadmill, synergies were similar across stages.ConclusionsThe high correlations between synergies across stages indicate that neuromuscular control strategies do not change as children walk at different speeds and slopes on a treadmill. However, the multiple significant differences in kinematics and kinetics between stages indicate real differences in movement pattern. This supports the theory that synergies are neurological in origin and not simply a response to the biomechanical task constraints.     

Use of muscle synergies and wavelet transforms to identify fatigue during squatting

The objective of this study was to supplement continuous wavelet transforms with muscle synergies in a fatigue analysis to better describe the combination of decreased firing frequency and altered activation profiles during dynamic muscle contractions. Nine healthy young individuals completed the dynamic tasks before and after they squatted with a standard Olympic bar until complete exhaustion. Electromyography (EMG) profiles were analyzed with a novel concatenated non-negative matrix factorization method that decomposed EMG signals into muscle synergies. Muscle synergy analysis provides the activation pattern of the muscles while continuous wavelet transforms output the temporal frequency content of the EMG signals. Synergy analysis revealed subtle changes in two-legged squatting after fatigue while differences in one-legged squatting were more pronounced and included the shift from a general co-activation of muscles in the pre-fatigue state to a knee extensor dominant weighting post-fatigue. Continuous wavelet transforms showed major frequency content decreases in two-legged squatting after fatigue while very few frequency changes occurred in one-legged squatting. It was observed that the combination of methods is an effective way of describing muscle fatigue and that muscle activation patterns play a very important role in maintaining the overall joint kinetics after fatigue.     

Functional muscle synergies constrain force production during postural tasks

We recently demonstrated that a set of five functional muscle synergies were sufficient to characterize both hindlimb muscle activity and active forces during automatic postural responses in cats standing at multiple postural configurations. This characterization depended critically upon the assumption that the endpoint force vector (synergy force vector) produced by the activation of each muscle synergy rotated with the limb axis as the hindlimb posture varied in the sagittal plane. Here, we used a detailed, 3D static model of the hindlimb to confirm that this assumption is biomechanically plausible: as we varied the model posture, simulated synergy force vectors rotated monotonically with the limb axis in the parasagittal plane (r2=0.94+/-0.08). We then tested whether a neural strategy of using these five functional muscle synergies provides the same force-generating capability as controlling each of the 31 muscles individually. We compared feasible force sets (FFSs) from the model with and without a muscle synergy organization. FFS volumes were significantly reduced with the muscle synergy organization (F=1556.01, p<0.01), and as posture varied, the synergy-limited FFSs changed in shape, consistent with changes in experimentally measured active forces. In contrast, nominal FFS shapes were invariant with posture, reinforcing prior findings that postural forces cannot be predicted by hindlimb biomechanics alone. We propose that an internal model for postural force generation may coordinate functional muscle synergies that are invariant in intrinsic limb coordinates, and this reduced-dimension control scheme reduces the set of forces available for postural control.     

A Novel Application of Musculoskeletal Ultrasound Imaging

Ultrasound is an attractive modality for imaging muscle and tendon motion during dynamic tasks and can provide a complementary methodological approach for biomechanical studies in a clinical or laboratory setting. Towards this goal, methods for quantification of muscle kinematics from ultrasound imagery are being developed based on image processing. The temporal resolution of these methods is typically not sufficient for highly dynamic tasks, such as drop-landing. We propose a new approach that utilizes a Doppler method for quantifying muscle kinematics. We have developed a novel vector tissue Doppler imaging (vTDI) technique that can be used to measure musculoskeletal contraction velocity, strain and strain rate with sub-millisecond temporal resolution during dynamic activities using ultrasound. The goal of this preliminary study was to investigate the repeatability and potential applicability of the vTDI technique in measuring musculoskeletal velocities during a drop-landing task, in healthy subjects. The vTDI measurements can be performed concurrently with other biomechanical techniques, such as 3D motion capture for joint kinematics and kinetics, electromyography for timing of muscle activation and force plates for ground reaction force. Integration of these complementary techniques could lead to a better understanding of dynamic muscle function and dysfunction underlying the pathogenesis and pathophysiology of musculoskeletal disorders.     

High day-to-day repeatability of lower extremity muscle activation patterns and joint biomechanics of dual-belt treadmill gait: A reliability study in healthy young adults

PurposeThe reliability of lower extremity muscle activation patterns has not been clearly studied in a dual-belt instrumented treadmill environment. The primary study objective was to quantify the day-to-day reliability of quadriceps, hamstrings, gastrocnemius and gluteus medius activation patterns in healthy young adult gait. Secondarily, the reliability of spatiotemporal, and knee/hip motion and moment-based gait outcomes was assessed.Scope: 20 young adults were recruited and tested on two separate days. Using standardized procedures, participants were prepared for surface electromyography and lower extremity motion capture. All individuals walked on a dual-belt instrumented treadmill while muscle activation, segment motions and ground reaction forces were recorded. Sagittal plane motion and net external sagittal and frontal plane moments were calculated. Discrete biomechanical and muscle activation measures were calculated, and non-negative matrix factorization extracted amplitude and temporal muscle activation features. Intraclass Correlation Coefficients, Standard Error of Measurement and Minimum Detectable Change were calculated.ConclusionsHigh to excellent Intraclass correlation coefficients were found between visits for most primary and secondary outcomes. The absolute and relative reliability, including Minimum Detectable Change values, provided in this study support the use of dual-belt instrumented treadmill walking as an acceptable medium to collect biomechanical and lower extremity EMG outcomes for future studies.     

Algorithms for task allocation in ants. (A study of temporal polyethism: Theory)

We present an algorithm for allocating individual ants to tasks that relies solely on task change being caused by the unavailability of work. We prove that such an algorithm will allocate the correct number of individuals to each job. Furthermore, we can demonstrate that if such an algorithm is used then an age structure emerges over the ants performing the various tasks. This matches closely with the weak temporal structure over tasks that is observed in Sendova-Franks and Franks (1993. Division of labour in ants nests within highly variable environments. (A study of temporal polyethism: experimental).Bull. math. Biol. 55, 75–96).     

Patterns of muscle activity during different behaviors in chicks: implications for neural control

The large behavioral repertoire that spans the embryonic and postembryonic stages of development make chicks an ideal system for identifying patterns of muscle activity that are common to different behaviors and those that are behavior-specific. The main goal of this work was to identify the similar and dissimilar aspects of the recruitment patterns and the regulation of muscle activity during three distinct postembryonic behaviors: walking, swimming and airstepping. We identified two synergies that were common to each of these behaviors. The synergies were not disrupted by the absence of FT1 activity in airstepping. Within each synergy the recruitment time, recruitment order and duration of activity were not rigid, but varied according to the context-specific resistance that the leg encountered. Unlike the other muscles, FT2 activity was not recruited as part of the same synergy in each behavior. When weight-bearing contact with the substrate did not occur, as in swimming and airstepping, as well as in walking in chicks with deafferented legs, FT2 activity was not recruited as part of either synergy, but was recruited during the time between them. Although not identical, embryonic motility and hatching motor pattern both show the two synergies described for the postembryonic behaviors. Like the latter behaviors, the synergies tolerated the absence of activity from specific muscles. Thus, we suggest that the CNS produces different behaviors using many of the same muscles by organizing the patterned activity around two common synergies while permitting the different muscles that participate in a synergy to be modified in tandem or on an individual basis. Furthermore, the common synergies are established early during prenatal development in chicks.Abbreviations aIF anterior iliofibularius - CF caudilioflexorius - CNS central nervous system - EMG electromyogram - FT femorotibialis - GL gastrocnemius lateralis - P1 protraction 1 - P2 protraction 2 - pIF posterior iliofibularis - R retraction - SA sartorius - TA tibialis anterior - pIF posterior iliofibularis - ms millisecond     

Optimization of muscle activity for task-level goals predicts complex changes in limb forces across biomechanical contexts

Optimality principles have been proposed as a general framework for understanding motor control in animals and humans largely based on their ability to predict general features movement in idealized motor tasks. However, generalizing these concepts past proof-of-principle to understand the neuromechanical transformation from task-level control to detailed execution-level muscle activity and forces during behaviorally-relevant motor tasks has proved difficult. In an unrestrained balance task in cats, we demonstrate that achieving task-level constraints center of mass forces and moments while minimizing control effort predicts detailed patterns of muscle activity and ground reaction forces in an anatomically-realistic musculoskeletal model. Whereas optimization is typically used to resolve redundancy at a single level of the motor hierarchy, we simultaneously resolved redundancy across both muscles and limbs and directly compared predictions to experimental measures across multiple perturbation directions that elicit different intra- and interlimb coordination patterns. Further, although some candidate task-level variables and cost functions generated indistinguishable predictions in a single biomechanical context, we identified a common optimization framework that could predict up to 48 experimental conditions per animal (n = 3) across both perturbation directions and different biomechanical contexts created by altering animals' postural configuration. Predictions were further improved by imposing experimentally-derived muscle synergy constraints, suggesting additional task variables or costs that may be relevant to the neural control of balance. These results suggested that reduced-dimension neural control mechanisms such as muscle synergies can achieve similar kinetics to the optimal solution, but with increased control effort (≈2×) compared to individual muscle control. Our results are consistent with the idea that hierarchical, task-level neural control mechanisms previously associated with voluntary tasks may also be used in automatic brainstem-mediated pathways for balance.     

Unexpected Dual Task Benefits on Cycling in Parkinson Disease and Healthy Adults: A Neuro-Behavioral Model

BackgroundWhen performing two tasks at once, a dual task, performance on one or both tasks typically suffers. People with Parkinson’s disease (PD) usually experience larger dual task decrements on motor tasks than healthy older adults (HOA). Our objective was to investigate the decrements in cycling caused by performing cognitive tasks with a range of difficulty in people with PD and HOAs.MethodsTwenty-eight participants with Parkinson’s disease and 20 healthy older adults completed a baseline cycling task with no secondary tasks and then completed dual task cycling while performing 12 tasks from six cognitive domains representing a wide range of difficulty.ResultsCycling was faster during dual task conditions than at baseline, and was significantly faster for six tasks (all p<.02) across both groups. Cycling speed improved the most during the easiest cognitive tasks, and cognitive performance was largely unaffected. Cycling improvement was predicted by task difficulty (p<.001). People with Parkinson’s disease cycled slower (p<.03) and showed reduced dual task benefits (p<.01) than healthy older adults.ConclusionsUnexpectedly, participants’ motor performance improved during cognitive dual tasks, which cannot be explained in current models of dual task performance. To account for these findings, we propose a model integrating dual task and acute exercise approaches which posits that cognitive arousal during dual tasks increases resources to facilitate motor and cognitive performance, which is subsequently modulated by motor and cognitive task difficulty. This model can explain both the improvement observed on dual tasks in the current study and more typical dual task findings in other studies.     

Feed-forward control of a redundant motor system

We describe a model of feed-forward control of a redundant motor system and validate it using, as examples, tasks of multi-finger force production. The model assumes the existence of two input signals at an upper level of the control hierarchy, related and unrelated to a task variable. Knowledge of the Jacobian of the system is assumed at the level of generation of elemental variables (variables at the level of effectors). Variance at the level of elemental variables is considered as the sum of two components, related and unrelated to variability in the task variable. An index of stabilization of the task variable is similarly introduced as to how it was done in several studies using the framework of the uncontrolled manifold hypothesis. Several phenomena have been simulated including data point distributions corresponding to presence and absence of force-stabilizing synergies in two-finger tasks, changes in synergies with practice, and changes in synergy indices in preparation to a fast action. The model is discussed in comparison to other models of control of multi-element systems based on feedback processes. It shows that patterns of structured variability in the space of elemental variables can result from feed-forward processes. Relations of the model to the equilibrium-point hypothesis are also discussed.     

Catapult effect in pole vaulting: Is muscle coordination determinant?

This study focused on the phase between the time of straightened pole and the maximum height (HP) of vaulter and aimed at determining the catapult effect in pole vaulting on HP. Seven experienced vaulters performed 5-10 vaults recorded by two video cameras, while the surface electromyography (sEMG) activity of 10 upper limbs muscles was recorded. HP was compared with an estimated maximum height (HP_est) allowing the computation of a push-off index. Muscle synergies were extracted from the sEMG activity profiles using a non-negative matrix factorization algorithm. No significant difference (p > 0.47) was found between HP_est (4.64 ± 0.21 m) and HP (4.69 ± 0.23 m). Despite a high inter-individual variability in sEMG profiles, two muscle synergies were extracted for all the subjects which accounted for 96.1 ± 2.9% of the total variance. While, the synergy activation coefficients were very similar across subjects, a higher variability was found in the muscle synergy vectors. Consequently, whatever the push-off index among the pole vaulters, the athletes used different muscle groupings (i.e., muscle synergy vectors) which were activated in a similar fashion (i.e., synergy activation coefficients). Overall, these results suggested that muscle coordination adopted between the time of straightened pole and the maximum height does not have a major influence on HP.     

Understanding muscle coordination of the human leg with dynamical simulations

Muscles coordinate multijoint motion by generating forces that cause reaction forces throughout the body. Thus, a muscle can redistribute existing segmental energy by accelerating some segments and decelerating others. In the process, a muscle may also produce or absorb energy, in which case its summed energetic effect on the segments is positive or negative, respectively. This Borelli Lecture shows how dynamical simulations derived from musculoskeletal models reveal muscle-induced segmental energy redistribution and muscle co-functions and synergies. Synergy occurs when co-excited muscles distribute segmental energy differently to execute the motor task. In maximum height jumping, high vertical velocity at lift-off occurs desirably at full body extension because biarticular leg muscles redistribute the energy produced by the uniarticular leg muscles. In pedaling, synergistic ankle plantarflexor force generation during leg extension allows the high energy produced by the uniarticular hip and knee extensors to be delivered to the crank. An analogous less-powerful flexor synergy exists during leg flexion. Hamstrings reduce crank deceleration during the leg extension-to-flexion transition by not only producing energy but delivering it to the crank through its contribution to the tangential (accelerating) crank force, though this hamstrings function occurs at the opposite (flexion-extension) transition when pedaling backwards. In walking, the eccentric quadriceps activity in early stance not only decelerates the leg but also accelerates the trunk. In mid-stance, the uni- and biarticular plantarflexors, by having opposite segmental energetic effects, act in synergy to support the whole body, so segmental potential and kinetic energy exchange can occur. To conclude, the extraction of unmeasurable variables from dynamical simulations emulating task kinematics, kinetics, and EMGs shows how the production of force and energy by individual muscles contribute to the energy flow among the individual segments during task execution.     

The development of rhythm regularity, neuromuscular strategies, and movement smoothness during repetitive reaching in typically developing children

IntroductionThis study examined the development of paced coordinated reaching characterized by the successful entrainment of the movement to an external pacer, synchronous muscle activations and movement smoothness.MethodsThirty children, 5–10 years of age, and ten adults were instructed to repeatedly reach for and move an object from a lower shelf to an upper shelf in time to a metronome. Surface electromyography data were recorded. Amplitude and cross-correlations were calculated on three muscle pairs crossing the shoulder and elbow. A motion capture system captured the space curve accelerations of hand, forearm and upper arm segments to quantify movement smoothness.ResultsThe 5–6 year old children showed the greatest amount of temporal variability, followed by 7–10 year olds and then the adults. Correlations between muscle pairs stabilizing the shoulder girdle were higher in each group as compared to the other two muscle pairs but the correlations for all pairs were consistently higher for adults. Movement smoothness for children 9–10 years of age was closer to an adult-like pattern with respect to control of the upper arm, but the hand segment had the greatest variability across groups.ConclusionsThe increased temporal variability and decreased movement smoothness of the hand and forearm segments suggest that control of more distal musculature may be more difficult in children. The neuromuscular strategies adopted by adults were more optimal than those adopted by children as reflected by smoother and more consistent reaching.     

No evidence of expertise-related changes in muscle synergies during rowing

The purpose of the present study was to determine whether expertise in rowing is driven by a specific structure in muscular coordination. We compared seven experienced rowers and eight untrained (i.e., inexperienced) subjects during rowing on an ergometer. Both surface electromyography activity and mechanical patterns (forces exerted at the handle and the foot-stretcher) were recorded during a high intensity rowing exercise. A non-negative matrix factorization was applied to 23 electromyographic patterns to differentiate muscle synergies. Results showed that expertise was not associated with different dimensionality in the electromyographic data and that three muscle synergies were sufficient to explain the majority of the variance accounted for (i.e., >90% of the total variance) in the two populations. The synergies extracted were similar in the two populations, with identical functional roles. While the temporal organization of the propulsive synergies was very similar, slight differences were found in the composition of the muscle synergies (muscle synergy vectors) between the two populations. The results suggests that rowing expertise would not require the development of novel muscle synergies but would imply intrinsic synergies already used in different behaviors. Performance in rowing is more probably linked to adjustments in the mechanical output of the muscle synergies rather than to differences in the shape and timing of their activations.     

Subject-specific tendon-aponeurosis definition in hill-type model predicts higher muscle forces in dynamic tasks

Neuromusculoskeletal models are a common method to estimate muscle forces. Developing accurate neuromusculoskeletal models is a challenging task due to the complexity of the system and large inter-subject variability. The estimation of muscles force is based on the mechanical properties of tendon-aponeurosis complex. Most neuromusculoskeletal models use a generic definition of the tendon-aponeurosis complex based on in vitro test, perhaps limiting their validity. Ultrasonography allows subject-specific estimates of the tendon-aponeurosis complex's mechanical properties. The aim of this study was to investigate the influence of subject-specific mechanical properties of the tendon-aponeurosis complex on a neuromusculoskeletal model of the ankle joint. Seven subjects performed isometric contractions from which the tendon-aponeurosis force-strain relationship was estimated. Hopping and running tasks were performed and muscle forces were estimated using subject-specific tendon-aponeurosis and generic tendon properties. Two ultrasound probes positioned over the muscle-tendon junction and the mid-belly were combined with motion capture to estimate the in vivo tendon and aponeurosis strain of the medial head of gastrocnemius muscle. The tendon-aponeurosis force-strain relationship was scaled for the other ankle muscles based on tendon and aponeurosis length of each muscle measured by ultrasonography. The EMG-driven model was calibrated twice - using the generic tendon definition and a subject-specific tendon-aponeurosis force-strain definition. The use of subject-specific tendon-aponeurosis definition leads to a higher muscle force estimate for the soleus muscle and the plantar-flexor group, and to a better model prediction of the ankle joint moment compared to the model estimate which used a generic definition. Furthermore, the subject-specific tendon-aponeurosis definition leads to a decoupling behaviour between the muscle fibre and muscle-tendon unit in agreement with previous experiments using ultrasonography. These results indicate the use of subject-specific tendon-aponeurosis definitions in a neuromusculoskeletal model produce better agreement with measured external loads and more physiological model behaviour.