Musculotendon Parameters in Lower Limb Models: Simplifications, Uncertainties, and Muscle Force Estimation Sensitivity

Muscle Moment Arm-Joint Angle Relations in the Hip, Knee, and Ankle: A Visualization of Datasets

Muscle moment arm is a property that associates muscle force with joint moment and is crucial to biomechanical analysis. In musculoskeletal simulations, the accuracy of moment arm is as important as that of muscle force, and calibrating moment arms in a musculoskeletal model requires data from anatomical measurements. Nonetheless, such data are elusive, and the complex relation between moment arm and joint angle can be unclear. Using common techniques in systematic review, we collected a total of 300 moment arm datasets from literature and visualized the muscle moment arm-joint angle relations in the human hip, knee, and ankle. The findings contribute to the analysis of musculoskeletal mechanics and providing reference regarding the experimental design for future moment arm measurements.

A registration strategy to characterize DTI-observed changes in skeletal muscle architecture due to passive shortening

Skeletal muscle architecture is a key determinant of muscle function. Architectural properties such as fascicle length, pennation angle, and curvature can be characterized using Diffusion Tensor Imaging (DTI), but acquiring these data during a contraction is not currently feasible. However, an image registration-based strategy may be able to convert muscle architectural properties observed at rest to their contracted state. As an initial step toward this long-term objective, the aim of this study was to determine if an image registration strategy could be used to convert the whole-muscle average architectural properties observed in the extended joint position to those of a flexed position, following passive rotation. DTI and high-resolution fat/water scans were acquired in the lower leg of seven healthy participants on a 3T MR system in + 20° and -10° ankle positions. The diffusion and anatomical images from the two positions were used to propagate DTI fiber-tracts from seed points along a mesh representation of the aponeurosis of fiber insertion. The -10° and + 20° anatomical images were registered and the displacement fields were used to transform the mesh and fiber-tracts from the + 20° to the -10° position. Student's paired t-tests were used to compare the mean architectural parameters between the original and transformed fiber-tracts. The whole-muscle average fiber-tract length, pennation angle, curvature, and physiological cross-sectional areas estimates did not differ significantly. DTI fiber-tracts in plantarflexion can be transformed to dorsiflexion position without significantly affecting the average architectural characteristics of the fiber-tracts. In the future, a similar approach could be used to evaluate muscle architecture in a contracted state.

Predicting prosthetic gait and the effects of induced stiff-knee gait

Prosthetic gait differs considerably from the unimpaired gait. Studying alterations in the gait patterns could help to understand different adaptation mechanisms adopted by these populations. This study investigated the effects of induced stiff-knee gait (SKG) on prosthetic and healthy gait patterns and the capabilities of predictive simulation. Self-selected speed gait of two participants was measured: one healthy subject and one knee disarticulation subject using a variable-damping microprocessor controlled knee prosthesis. Both performed unperturbed gait and gait with restricted knee flexion. Experimental joint angles and moments were computed using OpenSim and muscle activity was measured using surface electromyography (EMG). The differences between the conditions were analyzed using statistical parametric mapping (SPM). Predictive models based on optimal control were created to represent the participants. Additionally, a hypothetical unimpaired predictive model with the same anthropometric characteristics as the amputee was created. Some patterns observed in the experimental prosthetic gait were predicted by the models, including increased knee flexion moment on the contralateral side caused by SKG in both participants, which was statistically significant according to SPM. With the exception of the rectus femoris muscle, we also found overall good agreement between measured EMG and predicted muscle activation. We predicted more alterations in activation of the hip flexors than other muscle groups due to the amputation and in the activation of the biceps femoris short head, quadratus femoris, and tibialis anterior due to SKG. In summary, we demonstrated that the method applied in this study could predict gait alterations due to amputation of the lower limb or due to imposed SKG.

Influence of muscle packing on the three-dimensional architecture of rabbit M. plantaris

In their physiological condition, muscles are surrounded by connective tissue, other muscles and bone. These tissues exert transverse forces that change the three-dimensional shape of the muscle compared to its isolated condition, in which all surrounding tissues are removed. A change in shape affects the architecture of a muscle and therefore its mechanical properties. The rabbit M. plantaris is a multi-pennate calf muscle consisting of two compartments. A smaller, bi-pennate inner muscle compartment is embedded in a larger, uni-pennate outer compartment (Böl et al., 2015). As part of the calf, the plantaris is tightly packed between other muscles. It is unclear how packing affects the shape and architecture of the plantaris. Therefore, we examined the isolated and packed plantaris of the contralateral legs of three rabbits to determine the influence of the surrounding muscles on its shape and architectural properties using photogrammetric reconstruction and manual digitization, respectively. In the packed condition, the plantaris showed a 27% increase in fascicle pennation and a 54% increase in fascicle curvature compared to the isolated condition. Fascicle length was not affected by muscle packing. The change in muscle architecture occurred mainly in the outer compartment of the plantaris. Furthermore, the isolated plantaris showed a more circular shape and a reduced width of its muscle belly. It can be concluded that the packed plantaris is flattened by the forces exerted by the surrounding muscles, causing a complex architectural change. The data provided improve our understanding of muscle packages in general and can be used to develop and validate realistic three-dimensional muscle models.

The Biomechanics of Musculoskeletal Tissues during Activities of Daily Living: Dynamic Assessment Using Quantitative Transmission-Mode Ultrasound Techniques

The measurement of musculoskeletal tissue properties and loading patterns during physical activity is important for understanding the adaptation mechanisms of tissues such as bone, tendon, and muscle tissues, particularly with injury and repair. Although the properties and loading of these connective tissues have been quantified using direct measurement techniques, these methods are highly invasive and often prevent or interfere with normal activity patterns. Indirect biomechanical methods, such as estimates based on electromyography, ultrasound, and inverse dynamics, are used more widely but are known to yield different parameter values than direct measurements. Through a series of literature searches of electronic databases, including Pubmed, Embase, Web of Science, and IEEE Explore, this paper reviews current methods used for the in vivo measurement of human musculoskeletal tissue and describes the operating principals, application, and emerging research findings gained from the use of quantitative transmission-mode ultrasound measurement techniques to non-invasively characterize human bone, tendon, and muscle properties at rest and during activities of daily living. In contrast to standard ultrasound imaging approaches, these techniques assess the interaction between ultrasound compression waves and connective tissues to provide quantifiable parameters associated with the structure, instantaneous elastic modulus, and density of tissues. By taking advantage of the physical relationship between the axial velocity of ultrasound compression waves and the instantaneous modulus of the propagation material, these techniques can also be used to estimate the in vivo loading environment of relatively superficial soft connective tissues during sports and activities of daily living. This paper highlights key findings from clinical studies in which quantitative transmission-mode ultrasound has been used to measure the properties and loading of bone, tendon, and muscle tissue during common physical activities in healthy and pathological populations.

EMG-Constrained and Ultrasound-Informed Muscle-Tendon Parameter Estimation in Post-Stroke Hemiparesis

Secondary morphological and mechanical property changes in the muscle-tendon unit at the ankle joint are often observed in post-stroke individuals. These changes may alter the force generation capacity and affect daily activities such as locomotion. This work aimed to estimate subject-specific muscle-tendon parameters in individuals after stroke by solving the muscle redundancy problem using direct collocation optimal control methods based on experimental electromyography (EMG) signals and measured muscle fiber length. Subject-specific muscle-tendon parameters of the gastrocnemius, soleus, and tibialis anterior were estimated in seven post-stroke individuals and seven healthy controls. We found that the maximum isometric force, tendon stiffness and optimal fiber length in the post-stroke group were considerably lower than in the control group. We also computed the root mean square error between estimated and experimental values of muscle excitation and fiber length. The musculoskeletal model with estimated subject-specific muscle tendon parameters (from the muscle redundancy solver), yielded better muscle excitation and fiber length estimations than did scaled generic parameters. Our findings also showed that the muscle redundancy solver can estimate muscle-tendon parameters that produce force behavior in better accordance with the experimentally-measured value. These muscle-tendon parameters in the post-stroke individuals were physiologically meaningful and may shed light on treatment and/or rehabilitation planning.

A registration strategy to characterize DTI-observed changes in skeletal muscle architecture due to passive shortening

Skeletal muscle architecture is a key determinant of muscle function. Architectural properties such as fascicle length, pennation angle, and curvature can be characterized using Diffusion Tensor Imaging (DTI), but acquiring these data during a contraction is not currently feasible. However, an image registration-based strategy may be able to convert muscle architectural properties observed at rest to their contracted state. As an initial step toward this long-term objective, the aim of this study was to determine if an image registration strategy could be used to convert the whole-muscle average architectural properties observed in the extended joint position to those of a flexed position, following passive rotation. DTI and high-resolution fat/water scans were acquired in the lower leg of seven healthy participants on a 3T MR system in +20° (plantarflexion) and -10° (dorsiflexion) foot positions. The diffusion and anatomical images from the two positions were used to propagate DTI fiber-tracts from seed points along a mesh representation of the aponeurosis of fiber insertion. The -10° and +20° anatomical images were registered and the displacement fields were used to transform the mesh and fiber-tracts from the +20° to the -10° position. Student's paired t-tests were used to compare the mean architectural parameters between the original and transformed fiber-tracts. The whole-muscle average fiber-tract length, pennation angle, curvature, and physiological cross-sectional areas estimates did not differ significantly. DTI fiber-tracts in plantarflexion can be transformed to dorsiflexion position without significantly affecting the average architectural characteristics of the fiber-tracts. In the future, a similar approach could be used to evaluate muscle architecture in a contracted state.