Muscle Strength Identification Based on Isokinetic Testing and Spine Musculoskeletal Modeling

Subject-specific spinal musculoskeletal modeling can help understand the spinal loading mechanism during human locomotion. However, existing literature lacks methods to identify the maximum isometric strength of individual spinal muscles. In this study, a muscle strength identification method combining isokinetic testing and musculoskeletal simulations was proposed, and the influence of muscle synergy and intra-abdominal pressure (IAP) on identified spinal muscle strength was further discussed. A multibody dynamic model of the spinal musculoskeletal system was established and controlled by a feedback controller. Muscle strength parameters were adjusted based on the measured isokinetic moments, and muscle synergy vectors and the IAP piston model were further introduced. The results of five healthy subjects showed that the proposed method successfully identified the subject-specific spinal flexor/extensor strength. Considering the synergistic activations of antagonist muscles improved the correlation between the simulated and measured spinal moments, and the introduction of IAP slightly increased the identified spinal extensor strength. The established method is beneficial for understanding spinal loading distributions for athletes and patients with sarcopenia.

Machine learning-based prediction of joint moments based on kinematics in patients with cerebral palsy

Joint moments during gait provide valuable information for clinical decision-making in patients with cerebral palsy (CP). Joint moments are calculated based on ground reaction forces (GRF) using inverse dynamics models. Obtaining GRF from patients with CP is challenging. Typically developed (TD) individuals' joint moments were predicted from joint angles using machine learning, but no such study has been conducted on patients with CP. Accordingly, we aimed to predict the dorsi-plantar flexion, knee flexion-extension, hip flexion-extension, and hip adduction-abduction moments based on the trunk, pelvis, hip, knee, and ankle kinematics during gait in patients with CP and TD individuals using one-dimensional convolutional neural networks (CNN). The anonymized retrospective gait data of 329 TD (26 years ± 14, mass: 70 kg ± 15, height: 167 cm ± 89) and 917 CP (17 years ± 9, mass:47 kg ± 19, height:153 cm ± 36) individuals were evaluated and after applying inclusion-exclusion criteria, 132 TD and 622 CP patients with spastic diplegia were selected. We trained specific CNN models and evaluated their performance using isolated test subject groups based on normalized root mean square error (nRMSE) and Pearson correlation coefficient (PCC). Joint moments were predicted with nRMSE between 18.02% and 13.58% for the CP and between 12.55% and 8.58% for the TD groups, whereas with PCC between 0.85 and 0.93 for the CP and between 0.94 and 0.98 for the TD groups. Machine learning-based joint moment prediction from kinematics could replace conventional moment calculation in CP patients in the future, but the current level of prediction errors restricts its use for clinical decision-making today.

Personalisation of Plantarflexor Musculotendon Model Parameters in Children with Cerebral Palsy

Neuromusculoskeletal models can be used to evaluate aberrant muscle function in cerebral palsy (CP), for example by estimating muscle and joint contact forces during gait. However, to be accurate, models should include representative musculotendon parameters. We aimed to estimate personalised parameters that capture the mechanical behaviour of the plantarflexors in children with CP and typically developing (TD) children. Ankle angle (using motion capture), torque (using a load-cell), and medial gastrocnemius fascicle lengths (using ultrasound) were measured during slow passive ankle dorsiflexion rotation for thirteen children with spastic CP and thirteen TD children. Per subject, the measured rotation was input to a scaled OpenSim model to simulate the torque and fascicle length output. Musculotendon model parameters were personalised by the best match between simulated and experimental torque–angle and fascicle length-angle curves according to a least-squares fit. Personalised tendon slack lengths were significantly longer and optimal fibre lengths significantly shorter in CP than model defaults and than in TD. Personalised tendon compliance was substantially higher in both groups compared to the model default. The presented method to personalise musculotendon parameters will likely yield more accurate simulations of subject-specific muscle mechanics, to help us understand the effects of altered musculotendon properties in CP.