Editorial: Advances in Musculoskeletal Modeling and Their Application to Neurorehabilitation

NSF DARE-Transforming modeling in neurorehabilitation: Four threads for catalyzing progress

We present an overview of the Conference on Transformative Opportunities for Modeling in Neurorehabilitation held in March 2023. It was supported by the Disability and Rehabilitation Engineering (DARE) program from the National Science Foundation's Engineering Biology and Health Cluster. The conference brought together experts and trainees from around the world to discuss critical questions, challenges, and opportunities at the intersection of computational modeling and neurorehabilitation to understand, optimize, and improve clinical translation of neurorehabilitation. We organized the conference around four key, relevant, and promising Focus Areas for modeling: Adaptation & Plasticity, Personalization, Human-Device Interactions, and Modeling 'In-the-Wild'. We identified four common threads across the Focus Areas that, if addressed, can catalyze progress in the short, medium, and long terms. These were: (i) the need to capture and curate appropriate and useful data necessary to develop, validate, and deploy useful computational models (ii) the need to create multi-scale models that span the personalization spectrum from individuals to populations, and from cellular to behavioral levels (iii) the need for algorithms that extract as much information from available data, while requiring as little data as possible from each client (iv) the insistence on leveraging readily available sensors and data systems to push model-driven treatments from the lab, and into the clinic, home, workplace, and community. The conference archive can be found at (dare2023.usc.edu). These topics are also extended by three perspective papers prepared by trainees and junior faculty, clinician researchers, and federal funding agency representatives who attended the conference.

Scalable musculoskeletal model for dynamic simulations of upper body movement

A musculoskeletal (MSK) model is a valuable tool for assessing complex biomechanical problems, estimating joint torques during motion, optimizing motion in sports, and designing exoskeletons and prostheses. This study proposes an open-source upper body MSK model that supports biomechanical analysis of human motion. The MSK model of the upper body consists of 8 body segments (torso, head, left/right upper arm, left/right forearm, and left/right hand). The model has 20 degrees of freedom (DoFs) and 40 muscle torque generators (MTGs), which are constructed using experimental data. The model is adjustable for different anthropometric measurements and subject body characteristics: sex, age, body mass, height, dominant side, and physical activity. Joint limits are modeled using experimental dynamometer data within the proposed multi-DoF MTG model. The model equations are verified by simulating the joint range of motion (ROM) and torque; all simulation results have a good agreement with previously published research.

Does crouch alter the effects of neuromuscular impairments on gait? A simulation study

Cerebral palsy (CP) is a neurologic injury that impacts control of movement. Individuals with CP also often develop secondary impairments like weakness and contracture. Both altered motor control and secondary impairments influence how an individual walks after neurologic injury. However, understanding the complex interactions between and relative effects of these impairments makes analyzing and improving walking capacity in CP challenging. We used a sagittal-plane musculoskeletal model and neuromuscular control framework to simulate crouch and nondisabled gait. We perturbed each simulation by varying the number of synergies controlling each leg (altered control), and imposed weakness and contracture. A Bayesian Additive Regression Trees (BART) model was also used to parse the relative effects of each impairment on the muscle activations required for each gait pattern. By using these simulations to evaluate gait-pattern specific effects of neuromuscular impairments, we identified some advantages of crouch gait. For example, crouch tolerated 13 % and 22 % more plantarflexor weakness than nondisabled gait without and with altered control, respectively. Furthermore, BART demonstrated that plantarflexor weakness had twice the effect on total muscle activity required during nondisabled gait than crouch gait. However, crouch gait was also disadvantageous in the presence of vasti weakness: crouch gait increased the effects of vasti weakness on gait without and with altered control. These simulations highlight gait-pattern specific effects and interactions between neuromuscular impairments. Utilizing computational techniques to understand these effects can elicit advantages of gait deviations, providing insight into why individuals may select their gait pattern and possible interventions to improve energetics.

Scalable musculoskeletal model for dynamic simulations of lower body movement

A musculoskeletal (MSK) model is an important tool for analysing human motions, calculating joint torques during movement, enhancing sports activity, and developing exoskeletons and prostheses. To enable biomechanical investigation of human motion, this work presents an open-source lower body MSK model. The MSK model of the lower body consists of 7 body segments (pelvis, left/right thigh, left/right leg, and left/right foot). The model has 20 degrees of freedom (DoFs) and 28 muscle torque generators (MTGs), which are developed from experimental data. The model can be modified for different anthropometric measurements and subject body characteristics, including sex, age, body mass, height, physical activity, and skin temperature. The model is validated by simulating the torque within the range of motion (ROM) of isolated movements; all simulation findings exhibit a good level of agreement with the literature.

Modeling Musculoskeletal Dynamics during Gait: Evaluating the Best Personalization Strategy through Model Anatomical Consistency

No consensus exists on how to model human articulations within MSK models for the analysis of gait dynamics. We propose a method to evaluate joint models and we apply it to three models with different levels of personalization. The method evaluates the joint model’s adherence to the MSK hypothesis of negligible joint work by quantifying ligament and cartilage deformations resulting from joint motion; to be anatomically consistent, these deformations should be minimum. The contrary would require considerable external work to move the joint, violating a strong working hypothesis and raising concerns about the credibility of the MSK outputs. Gait analysis and medical resonance imaging (MRI) from ten participants were combined to build lower limb subject-specific MSK models. MRI-reconstructed anatomy enabled three levels of personalization using different ankle joint models, in which motion corresponded to different ligament elongation and cartilage co-penetration. To estimate the impact of anatomical inconsistency in MSK outputs, joint internal forces resulting from tissue deformations were computed for each joint model and MSK simulations were performed ignoring or considering their contribution. The three models differed considerably for maximum ligament elongation and cartilage co-penetration (between 5.94 and 50.69% and between −0.53 and −5.36 mm, respectively). However, the model dynamic output from the gait simulations were similar. When accounting for the internal forces associated with tissue deformation, outputs changed considerably, the higher the personalization level the smaller the changes. Anatomical consistency provides a solid method to compare different joint models. Results suggest that consistency grows with personalization, which should be tailored according to the research question. A high level of anatomical consistency is recommended when individual specificity and the behavior of articular structures is under investigation.