Towards a realistic biomechanical model of the thumb: the choice of kinematic description may be more critical than the solution method or the variability/uncertainty of musculoskeletal parameters.

A biomechanical model of the thumb can help researchers and clinicians understand the clinical problem of how anatomical variability contributes to the variability of outcomes of surgeries to restore thumb function. We lack a realistic biomechanical model of the thumb because of the variability/uncertainty of musculoskeletal parameters, the multiple proposed kinematic descriptions and methods to solve the muscle redundancy problem, and the paucity of data to validate the model with in vivo coordination patterns and force output. We performed a multi-stage validation of a biomechanical computer model against our measurements of maximal static thumbtip force and fine-wire electromyograms (EMG) from 8 thumb muscles in each of five orthogonal directions in key and opposition pinch postures. A low-friction point-contact at the thumbtip ensured that subjects did not produce thumbtip torques during force production. The 3-D, 8-muscle biomechanical thumb model uses a 5-axis kinematic description with orthogonal and intersecting axes of rotation at the carpometacarpal and metacarpophalangeal joints. We represented the 50 musculoskeletal parameters of the model as stochastic variables based on experimental data, and ran Monte Carlo simulations in the "inverse" and "forward" directions for 5000 random instantiations of the model. Two inverse simulations (predicting the distribution of maximal static thumbtip forces and the muscle activations that maximized force) showed that: the model reproduces at most 50% of the 80 EMG distributions recorded (eight muscle excitations in 5 force directions in two postures); and well-directed thumbtip forces of adequate magnitude are predicted only if accompanied by unrealistically large thumbtip torques (0.64+/-0.28Nm). The forward simulation (which fed the experimental distributions of EMG through random instantiations of the model) resulted in misdirected thumbtip force vectors (within 74.3+/-24.5 degrees from the desired direction) accompanied by doubly large thumbtip torques (1.32+/-0.95Nm). Taken together, our results suggest that the variability and uncertainty of musculoskeletal parameters and the choice of solution method are not the likely reason for the unrealistic predictions obtained. Rather, the kinematic description of the thumb we used is not representative of the transformation of net joint torques into thumbtip forces/torques in the human thumb. Future efforts should focus on validating alternative kinematic descriptions of the thumb.