Insights into muscle metabolic energetics: Modelling muscle-tendon mechanics and metabolic rates during walking across speeds
Israel Luis,Maarten Afschrift,Friedl De Groote,Elena M. Gutierrez-Farewik
DOI: https://doi.org/10.1371/journal.pcbi.1012411
2024-09-14
PLoS Computational Biology
Abstract:The metabolic energy rate of individual muscles is impossible to measure without invasive procedures. Prior studies have produced models to predict metabolic rates based on experimental observations of isolated muscle contraction from various species. Such models can provide reliable predictions of metabolic rates in humans if muscle properties and control are accurately modeled. This study aimed to examine how muscle-tendon model individualization and metabolic energy models influenced estimation of muscle-tendon states and time-series metabolic rates, to evaluate the agreement with empirical data, and to provide predictions of the metabolic rate of muscle groups and gait phases across walking speeds. Three-dimensional musculoskeletal simulations with prescribed kinematics and dynamics were performed. An optimal control formulation was used to compute muscle-tendon states with four levels of individualization, ranging from a scaled generic model and muscle controls based on minimal activations, inclusion of calibrated muscle passive forces, personalization of Achilles and quadriceps tendon stiffnesses, to finally informing muscle controls with electromyography. We computed metabolic rates based on existing models. Simulations with calibrated passive forces and personalized tendon stiffness most accurately estimate muscle excitations and fiber lengths. Interestingly, the inclusion of electromyography did not improve our estimates. The whole-body average metabolic cost was better estimated with a subset of metabolic energy models. We estimated metabolic rate peaks near early stance, pre-swing, and initial swing at all walking speeds. Plantarflexors accounted for the highest cost among muscle groups at the preferred speed and were similar to the cost of hip adductors and abductors combined. Also, the swing phase accounted for slightly more than one-quarter of the total cost in a gait cycle, and its relative cost decreased with walking speed. Our prediction might inform the design of assistive devices and rehabilitation treatment. The code and experimental data are available online. Simulations of the musculoskeletal system can provide insights into motion mechanics and energetics. To increase confidence in the approach, it is important to establish how predictions are affected by modelling assumptions. In this study, we evaluated how simulations with increasing levels of sophistication and metabolic energy models from the literature agreed with recorded muscle excitations, fiber lengths, and whole-body metabolic rates during walking at different speeds. We found that calibrating muscle passive components: muscle passive forces, and tendon stiffness, estimated better muscle excitations and fiber lengths than using original musculoskeletal model values. The metabolic energy cost of walking was closely estimated with only a subset of metabolic energy models. We estimated that the metabolic energy demands are the highest during three periods of the gait cycle—early stance, pre-swing, and initial swing, at all walking speeds. In all walking speeds, muscles that support plantarflexion accounted for the highest metabolic cost among all the muscle groups. As walking speed increases, the metabolic cost of the stance phase also increases. Our detailed description of metabolic costs at muscle groups and gait phases might help to design assistive devices and rehabilitation treatments.
biochemical research methods,mathematical & computational biology