Connecting the wrist to the hand: A simulation study exploring changes in thumb-tip endpoint force following wrist surgery

Explainable AI Elucidates Musculoskeletal Biomechanics: A Case Study Using Wrist Surgeries

As datasets increase in size and complexity, biomechanists have turned to artificial intelligence (AI) to aid their analyses. This paper explores how explainable AI (XAI) can enhance the interpretability of biomechanics data derived from musculoskeletal simulations. We use machine learning to classify the simulated lateral pinch data as belonging to models with healthy or one of two types of surgically altered wrists. This simulation-based classification task is analogous to using biomechanical movement and force data to clinically diagnose a pathological state. The XAI describes which musculoskeletal features best explain the classifications and, in turn, the pathological states, at both the local (individual decision) level and global (entire algorithm) level. We demonstrate that these descriptions agree with assessments in the literature and additionally identify the blind spots that can be missed with traditional statistical techniques.

Sensitivity Analysis of Upper Limb Musculoskeletal Models During Isometric and Isokinetic Tasks

Sensitivity coefficients are used to understand how errors in subject-specific musculoskeletal model parameters influence model predictions. Previous sensitivity studies in the lower limb calculated sensitivity using perturbations that do not fully represent the diversity of the population. Hence, the present study performs sensitivity analysis in the upper limb using a large synthetic dataset to capture greater physiological diversity. The large dataset (n = 401 synthetic subjects) was created by adjusting maximum isometric force, optimal fiber length, pennation angle, and bone mass to induce atrophy, hypertrophy, osteoporosis, and osteopetrosis in two upper limb musculoskeletal models. Simulations of three isometric and two isokinetic upper limb tasks were performed using each synthetic subject to predict muscle activations. Sensitivity coefficients were calculated using three different methods (two point, linear regression, and sensitivity functions) to understand how changes in Hill-type parameters influenced predicted muscle activations. The sensitivity coefficient methods were then compared by evaluating how well the coefficients accounted for measurement uncertainty. This was done by using the sensitivity coefficients to predict the range of muscle activations given known errors in measuring musculoskeletal parameters from medical imaging. Sensitivity functions were found to best account for measurement uncertainty. Simulated muscle activations were most sensitive to optimal fiber length and maximum isometric force during upper limb tasks. Importantly, the level of sensitivity was muscle and task dependent. These findings provide a foundation for how large synthetic datasets can be applied to capture physiologically diverse populations and understand how model parameters influence predictions.

Predictions of thumb, hand, and arm muscle parameters derived using force measurements of varying complexity and neural networks

Subject-specific musculoskeletal models are a promising avenue for personalized healthcare. However, current methods for producing personalized models require dense, biomechanical datasets that include expensive and time-consuming physiological measurements. For personalized models to be clinically useful, we must be able to rapidly generate models from simple, easy to collect data. In this context, the objective of this paper is to evaluate if and how simple data, namely height/weight and pinch force data, can be used to achieve model personalization via machine learning. Using simulated lateral pinch force measurements from a synthetic population of 40,000 randomly generated subjects, we train neural networks to estimate four Hill-type muscle model parameters and bone density. We compare parameter estimates to the true parameters of 10,000 additional synthetic subjects. We also generate new personalized models using the parameter estimates and perform new lateral pinch simulations to compare predicted forces using these personalized models to those generated using a baseline model. We demonstrate that increasing force measurement complexity reduces the root-mean-square error in the majority of parameter estimates. Additionally, musculoskeletal models using neural network-based parameter estimates provide up to an 80% reduction in absolute error in simulated forces when compared to a generic model. Thus, easily obtained force measurements may be suitable for personalizing models of the thumb, although extending the method to more tasks and models involving other joints likely requires additional measurements.

Inverse distance weighting to rapidly generate large simulation datasets

Obtaining large biomechanical datasets for machine learning is an ongoing challenge. Physics-based simulations offer one approach for generating large datasets, but many simulation methods, such as computed muscle control (CMC), are computationally costly. In contrast, interpolation methods, such as inverse distance weighting (IDW), are computationally fast. We examined whether IDW is a low-cost and accurate approach for interpolating muscle activations from CMC.IDW was evaluated using lateral pinch simulations in OpenSim. Simulated pinch data were organized into grids of varying sparsity (high, medium, and low density), where each grid point represented the muscle activations associated with a unique combination of mass and height of a young adult. For each grid, muscle activations were calculated via CMC and IDW for 108 random mass-height pairs that were not coincident with simulation grid vertices. We evaluated the interpolation errors from IDW for each grid, as well as the sensitivity of lateral pinch force to these errors. The root mean square error (RMSE) associated with interpolated muscle activations decreased with increasing grid density and never exceeded 4%. While CMC received a target thumb-tip force of 40 N, errors from the interpolated muscle activations never impacted the simulated force magnitude by more than 0.1 N. Furthermore, the computation time for CMC simulations averaged 4.22 core-minutes, while IDW averaged 0.95 core-seconds per mass-height pair.These results indicate IDW is a practical approach for rapidly estimating muscle activations from sparse CMC datasets. Future works could adapt our IDW approach to evaluate other tasks, biomechanical features, and/or populations.

A Musculoskeletal Model of the Hand and Wrist Capable of Simulating Functional Tasks

OBJECTIVE

The purpose of this work was to develop an open-source musculoskeletal model of the hand and wrist and to evaluate its performance during simulations of functional tasks.

METHODS

The current model was developed by adapting and expanding upon existing models. An optimal control theory framework that combines forward-dynamics simulations with a simulated-annealing optimization was used to simulate maximum grip and pinch force. Active and passive hand opening were simulated to evaluate coordinated kinematic hand movements.

RESULTS

The model's maximum grip force production matched experimental measures of grip force, force distribution amongst the digits, and displayed sensitivity to wrist flexion. Simulated lateral pinch strength replicated in vivo palmar pinch strength data. Additionally, predicted activations for 7 of 8 muscles fell within variability of EMG data during palmar pinch. The active and passive hand opening simulations predicted reasonable activations and demonstrated passive motion mimicking tenodesis, respectively.

CONCLUSION

This work advances simulation capabilities of hand and wrist models and provides a foundation for future work to build upon.

SIGNIFICANCE

This is the first open-source musculoskeletal model of the hand and wrist to be implemented during both functional kinetic and kinematic tasks. We provide a novel simulation framework to predict maximal grip and pinch force which can be used to evaluate how potential surgical and rehabilitation interventions influence these functional outcomes while requiring minimal experimental data.

Musculoskeletal Modeling of the Wrist via a Multi Body Simulation

In this study, three different musculoskeletal modeling approaches were compared to each other. The objective was to show the possibilities in the case of a simple mechanical model of the wrist, using a simple multi-body-simulation (MBS) model, and using a more complex and patient-specific adaptable wrist joint MBS model. Musculoskeletal modeling could be a useful alternative, which can be practiced as a non-invasive approach to investigate body motion and internal loads in a wide range of conditions. The goal of this study was the introduction of computer-based modelling of the physiological wrist with (MBS-) models focused on the muscle and joint forces acting on the wrist.

Classifying muscle parameters with artificial neural networks and simulated lateral pinch data

Objective

Hill-type muscle models are widely employed in simulations of human movement. Yet, the parameters underlying these models are difficult or impossible to measure in vivo. Prior studies demonstrate that Hill-type muscle parameters are encoded within dynamometric data. But, a generalizable approach for estimating these parameters from dynamometric data has not been realized. We aimed to leverage musculoskeletal models and artificial neural networks to classify one Hill-type muscle parameter (maximum isometric force) from easily measurable dynamometric data (simulated lateral pinch force). We tested two neural networks (feedforward and long short-term memory) to identify if accounting for dynamic behavior improved accuracy.

Methods

We generated four datasets via forward dynamics, each with increasing complexity from adjustments to more muscles. Simulations were grouped and evaluated to show how varying the maximum isometric force of thumb muscles affects lateral pinch force. Both neural networks classified these groups from lateral pinch force alone.

Results

Both neural networks achieved accuracies above 80% for datasets which varied only the flexor pollicis longus and/or the abductor pollicis longus. The inclusion of muscles with redundant functions dropped model accuracies to below 30%. While both neural networks were consistently more accurate than random guess, the long short-term memory model was not consistently more accurate than the feedforward model.

Conclusion

Our investigations demonstrate that artificial neural networks provide an inexpensive, data-driven approach for approximating Hill-type muscle-tendon parameters from easily measurable data. However, muscles of redundant function or of little impact to force production make parameter classification more challenging.

Anthropometric scaling of musculoskeletal models of the hand captures age-dependent differences in lateral pinch force

Musculoskeletal models and computer simulations enable non-invasive study of muscle function and contact forces. Hand models are useful for understanding the complexities of hand strength, precision movement, and the dexterity required during daily activities. Yet, generic models fail to accurately represent the entire scope of the population, while subject-specific models are labor-intensive to create. The objective of this study was to assess the efficacy of scaled generic models to represent the broad spectrum of strength profiles across the lifespan. We examined one hundred lateral pinch simulations using a generic model of the wrist and thumb anthropometrically scaled to represent the full range of heights reported for four ages across childhood, puberty, older adolescence, and adulthood. We evaluated maximum lateral pinch force produced, muscle control strategies, and the effect of linearly scaling the maximum isometric force. Our simulations demonstrated three main concepts. First, anthropometric scaling could capture age-dependent differences in pinch strength. Second, a generic muscle control strategy is not representative of all populations. Lastly, simulations do not employ optimal fiber length to complete a lateral pinch task. These results demonstrate the potential of anthropometrically-scaled models to study hand strength across the lifespan, while also highlighting that muscle control strategies may adapt as we age. The results also provide insight to the force-length relationship of thumb muscles during lateral pinch. We conclude that anthropometric scaling can accurately represent age characteristics of the population, but subject-specific models are still necessary to represent individuals.