Neuro-musculoskeletal simulation of instrumented contracture and spasticity assessment in children with cerebral palsy

Experiment-guided tuning of muscle-tendon parameters to estimate muscle fiber lengths and passive forces

The workflow to simulate motion with recorded data usually starts with selecting a generic musculoskeletal model and scaling it to represent subject-specific characteristics. Simulating muscle dynamics with muscle-tendon parameters computed from existing scaling methods in literature, however, yields some inconsistencies compared to measurable outcomes. For instance, simulating fiber lengths and muscle excitations during walking with linearly scaled parameters does not resemble established patterns in the literature. This study presents a tool that leverages reported in vivo experimental observations to tune muscle-tendon parameters and evaluates their influence in estimating muscle excitations and metabolic costs during walking. From a scaled generic musculoskeletal model, we tuned optimal fiber length, tendon slack length, and tendon stiffness to match reported fiber lengths from ultrasound imaging and muscle passive force-length relationships to match reported in vivo joint moment-angle relationships. With tuned parameters, muscle contracted more isometrically, and soleus's operating range was better estimated than with linearly scaled parameters. Also, with tuned parameters, on/off timing of nearly all muscles' excitations in the model agreed with reported electromyographic signals, and metabolic rate trajectories varied significantly throughout the gait cycle compared to linearly scaled parameters. Our tool, freely available online, can customize muscle-tendon parameters easily and be adapted to incorporate more experimental data.

Muscle synergy-informed neuromusculoskeletal modelling to estimate knee contact forces in children with cerebral palsy

Cerebral palsy (CP) includes a group of neurological conditions caused by damage to the developing brain, resulting in maladaptive alterations of muscle coordination and movement. Estimates of joint moments and contact forces during locomotion are important to establish the trajectory of disease progression and plan appropriate surgical interventions in children with CP. Joint moments and contact forces can be estimated using electromyogram (EMG)-informed neuromusculoskeletal models, but a reduced number of EMG sensors would facilitate translation of these computational methods to clinics. This study developed and evaluated a muscle synergy-informed neuromusculoskeletal modelling approach using EMG recordings from three to four muscles to estimate joint moments and knee contact forces of children with CP and typically developing (TD) children during walking. Using only three to four experimental EMG sensors attached to a single leg and leveraging an EMG database of walking data of TD children, the synergy-informed approach estimated total knee contact forces comparable to those estimated by EMG-assisted approaches that used 13 EMG sensors (children with CP, n = 3, R2 = 0.95 ± 0.01, RMSE = 0.40 ± 0.14 BW; TD controls, n = 3, R2 = 0.93 ± 0.07, RMSE = 0.19 ± 0.05 BW). The proposed synergy-informed neuromusculoskeletal modelling approach could enable rapid evaluation of joint biomechanics in children with unimpaired and impaired motor control within a clinical environment.

Identification of Neural and Non-Neural Origins of Joint Hyper-Resistance Based on a Novel Neuromechanical Model

Joint hyper-resistance is a common symptom in neurological disorders. It has both neural and non-neural origins, but it has been challenging to distinguish different origins based on clinical tests alone. Combining instrumented tests with parameter identification based on a neuromechanical model may allow us to dissociate the different origins of joint hyper-resistance in individual patients. However, this requires that the model captures the underlying mechanisms. Here, we propose a neuromechanical model that, in contrast to previously proposed models, accounts for muscle short-range stiffness (SRS) and its interaction with muscle tone and reflex activity. We collected knee angle trajectories during the pendulum test in 15 children with cerebral palsy (CP) and 5 typically developing children. We did the test in two conditions - hold and pre-movement - that have been shown to alter knee movement. We modeled the lower leg as an inverted pendulum actuated by two antagonistic Hill-type muscles extended with SRS. Reflex activity was modeled as delayed, linear feedback from muscle force. We estimated neural and non-neural parameters by optimizing the fit between simulated and measured knee angle trajectories during the hold condition. The model could fit a wide range of knee angle trajectories in the hold condition. The model with personalized parameters predicted the effect of pre-movement demonstrating that the model captured the underlying mechanism and subject-specific deficits. Our model may help with the identification of neural and non-neural origins of joint hyper-resistance and thereby opens perspectives for improved diagnosis and treatment selection in children with spastic CP, but such applications require further studies to establish the method's reliability.

SpasticSim: a synthetic data generation method for upper limb spasticity modelling in neurorehabilitation

In neurorehabilitation, assessment of functional problems is essential to define optimal rehabilitation treatments. Usually, this assessment process requires distinguishing between impaired and non-impaired behavior of limbs. One of the common muscle motor disorders affecting limbs is spasticity, which is complicated to quantify objectively due to the complex nature of motor control. Thus, the lack of heterogeneous samples of patients constituting an acceptable amount of data is an obstacle which is relevant to understanding the behavior of spasticity and, consequently, quantifying it. In this article, we use the 3D creation suite Blender combined with the MBLab add-on to generate synthetic samples of human body models, aiming to be as sufficiently representative as possible to real human samples. Exporting these samples to OpenSim and performing four specific upper limb movements, we analyze the muscle behavior by simulating the six degrees of spasticity contemplated by the Modified Ashworth Scale (MAS). The complete dataset of patients and movements is open-source and available for future research. This approach advocates the potential to generate synthetic data for testing and validating musculoskeletal models.

A Velocity Stretch Reflex Threshold Based on Muscle-Tendon Unit Peak Acceleration to Detect Possible Occurrences of Spasticity during Gait in Children with Cerebral Palsy

Spasticity might affect gait in children with cerebral palsy. Quantifying its occurrence during locomotion is challenging. One approach is to determine kinematic stretch reflex thresholds, usually on the velocity, during passive assessment and to search for their exceedance during gait. These thresholds are determined through EMG-Onset detection algorithms, which are variable in performance and sensitive to noisy data, and can therefore lack consistency. This study aimed to evaluate the feasibility of determining the velocity stretch reflex threshold from maximal musculotendon acceleration. Eighteen children with CP were recruited and underwent clinical gait analysis and a full instrumented assessment of their soleus, gastrocnemius lateralis, semitendinosus, and rectus femoris spasticity, with EMG, kinematics, and applied forces being measured simultaneously. Using a subject-scaled musculoskeletal model, the acceleration-based stretch reflex velocity thresholds were determined and compared to those based on EMG-Onset determination. Their consistencies according to physiological criteria, i.e., if the timing of the threshold was between the beginning of the stretch and the spastic catch, were evaluated. Finally, two parameters designed to evaluate the occurrence of spasticity during gait, i.e., the proportion of the gait trial time with a gait velocity above the velocity threshold and the number of times the threshold was exceeded, were compared. The proposed method produces velocity stretch reflex thresholds close to the EMG-based ones. For all muscles, no statistical difference was found between the two parameters designed to evaluate the occurrence of spasticity during gait. Contrarily to the EMG-based methods, the proposed method always provides physiologically consistent values, with median electromechanical delays of between 50 and 130 ms. For all subjects, the semitendinosus velocity during gait usually exceeded its stretch reflex threshold, while it was less frequent for the three other muscles. We conclude that a velocity stretch reflex threshold, based on musculotendon acceleration, is a reliable substitute for EMG-based ones.

How mechanics of individual muscle-tendon units define knee and ankle joint function in health and cerebral palsy-a narrative review

This study reviews the relationship between muscle-tendon biomechanics and joint function, with a particular focus on how cerebral palsy (CP) affects this relationship. In healthy individuals, muscle size is a critical determinant of strength, with muscle volume, cross-sectional area, and moment arm correlating with knee and ankle joint torque for different isometric/isokinetic contractions. However, in CP, impaired muscle growth contributes to joint pathophysiology even though only a limited number of studies have investigated the impact of deficits in muscle size on pathological joint function. As muscles are the primary factors determining joint torque, in this review two main approaches used for muscle force quantification are discussed. The direct quantification of individual muscle forces from their relevant tendons through intraoperative approaches holds a high potential for characterizing healthy and diseased muscles but poses challenges due to the invasive nature of the technique. On the other hand, musculoskeletal models, using an inverse dynamic approach, can predict muscle forces, but rely on several assumptions and have inherent limitations. Neither technique has become established in routine clinical practice. Nevertheless, identifying the relative contribution of each muscle to the overall joint moment would be key for diagnosis and formulating efficient treatment strategies for patients with CP. This review emphasizes the necessity of implementing the intraoperative approach into general surgical practice, particularly for joint correction operations in diverse patient groups. Obtaining in vivo data directly would enhance musculoskeletal models, providing more accurate force estimations. This integrated approach can improve the clinicians' decision-making process and advance treatment strategies by predicting changes at the muscle and joint levels before interventions, thus, holding the potential to significantly enhance clinical outcomes.

Influence of different calibration methods on surface electromyography-informed musculoskeletal models with few input signals

BACKGROUND

Although model personalization is critical when assessing individuals with morphological or neurological abnormalities, or even non-disabled subjects, its translation into routine clinical settings is hampered by the cumbersomeness of experimental data acquisition and lack of resources, which are linked to high costs and long processing pipelines. Quantifying the impact of neglecting subject-specific information in simulations that aim to estimate muscle forces with surface electromyography informed modeling approaches, can address their potential in relevant clinical questions. The present study investigates how different methods to fine-tune subject-specific neuromuscular parameters, reducing the number of electromyography input data, could affect the estimation of the unmeasured excitations and the musculotendon forces.

METHODS

Three-dimensional motion analysis was performed on 8 non-disabled adult subjects and 13 electromyographic signals captured. Four neuromusculoskeletal models were created for 8 participants: a reference model driven by a large set of sEMG signals; two models informed by four electromyographic signals but calibrated in different fashions; a model based on static optimization.

FINDINGS

The electromyography-informed models better predicted experimental excitations, including the unmeasured ones. The model based on static optimization obtained less reliable predictions of the experimental data. When comparing the different reduced models, no major differences were observed, suggesting that the less complex model may suffice for predicting muscle forces with a small set of input in clinical gait analysis tasks.

INTERPRETATION

Quantitative model performance evaluation in different conditions provides an objective indication of which method yields the most accurate prediction when a small set of electromyographic recordings is available.

Modeling Human Suboptimal Control: A Review

This review paper provides an overview of the approaches to model neuromuscular control, focusing on methods to identify nonoptimal control strategies typical of populations with neuromuscular disorders or children. Where possible, the authors tightened the description of the methods to the mechanisms behind the underlying biomechanical and physiological rationale. They start by describing the first and most simplified approach, the reductionist approach, which splits the role of the nervous and musculoskeletal systems. Static optimization and dynamic optimization methods and electromyography-based approaches are summarized to highlight their limitations and understand (the need for) their developments over time. Then, the authors look at the more recent stochastic approach, introduced to explore the space of plausible neural solutions, thus implementing the uncontrolled manifold theory, according to which the central nervous system only controls specific motions and tasks to limit energy consumption while allowing for some degree of adaptability to perturbations. Finally, they explore the literature covering the explicit modeling of the coupling between the nervous system (acting as controller) and the musculoskeletal system (the actuator), which may be employed to overcome the split characterizing the reductionist approach.

Predicting Gait Patterns of Children With Spasticity by Simulating Hyperreflexia

Spasticity is a common impairment within pediatric neuromusculoskeletal disorders. How spasticity contributes to gait deviations is important for treatment selection. Our aim was to evaluate the pathophysiological mechanisms underlying gait deviations seen in children with spasticity, using predictive simulations. A cluster analysis was performed to extract distinct gait patterns from experimental gait data of 17 children with spasticity to be used as comparative validation data. A forward dynamic simulation framework was employed to predict gait with either velocity- or force-based hyperreflexia. This framework entailed a generic musculoskeletal model controlled by reflexes and supraspinal drive, governed by a multiobjective cost function. Hyperreflexia values were optimized to enable the simulated gait to best match experimental gait patterns. Three experimental gait patterns were extracted: (1) increased knee flexion, (2) increased ankle plantar flexion, and (3) increased knee flexion and ankle plantar flexion when compared with typical gait. Overall, velocity-based hyperreflexia outperformed force-based hyperreflexia. The first gait pattern could mostly be explained by rectus femoris and hamstrings velocity-based hyperreflexia, the second by gastrocnemius velocity-based hyperreflexia, and the third by gastrocnemius, soleus, and hamstrings velocity-based hyperreflexia. This study shows how velocity-based hyperreflexia from specific muscles contributes to different spastic gait patterns, which may help in providing targeted treatment.

The gait pattern and not the femoral morphology is the main contributor to asymmetric hip joint loading

Gait asymmetry and skeletal deformities are common in many children with cerebral palsy (CP). Changes of the hip joint loading, i.e. hip joint contact force (HJCF), can lead to pathological femoral growth. A child's gait pattern and femoral morphology affect HJCFs. The twofold aim of this study was to (1) evaluate if the asymmetry in HJCFs is higher in children with CP compared to typically developing (TD) children and (2) identify if the bony morphology or the subject-specific gait pattern is the main contributor to asymmetric HJCFs. Magnetic resonance images (MRI) and three-dimensional gait analysis data of twelve children with CP and fifteen TD children were used to create subject-specific musculoskeletal models and calculate HJCF using OpenSim. Root-mean-square-differences between left and right HJCF magnitude and orientation were computed and compared between participant groups (CP versus TD). Additionally, the influence on HJCF asymmetries solely due to the femoral morphology and solely due to the gait pattern was quantified. Our findings demonstrate that the gait pattern is the main contributor to asymmetric HJCFs in CP and TD children. Children with CP have higher HJCF asymmetries which is probably the result of larger asymmetries in their gait pattern compared to TD children. The gained insights from our study highlight that clinical interventions should focus on normalizing the gait pattern and therefore the hip joint loading to avoid the development of femoral deformities.

Dysfunctional neuro-muscular mechanisms explain gradual gait changes in prodromal spastic paraplegia

BACKGROUND

In Hereditary Spastic Paraplegia (HSP) type 4 (SPG4) a length-dependent axonal degeneration in the cortico-spinal tract leads to progressing symptoms of hyperreflexia, muscle weakness, and spasticity of lower extremities. Even before the manifestation of spastic gait, in the prodromal phase, axonal degeneration leads to subtle gait changes. These gait changes - depicted by digital gait recording - are related to disease severity in prodromal and early-to-moderate manifest SPG4 participants.

METHODS

We hypothesize that dysfunctional neuro-muscular mechanisms such as hyperreflexia and muscle weakness explain these disease severity-related gait changes of prodromal and early-to-moderate manifest SPG4 participants. We test our hypothesis in computer simulation with a neuro-muscular model of human walking. We introduce neuro-muscular dysfunction by gradually increasing sensory-motor reflex sensitivity based on increased velocity feedback and gradually increasing muscle weakness by reducing maximum isometric force.

RESULTS

By increasing hyperreflexia of plantarflexor and dorsiflexor muscles, we found gradual muscular and kinematic changes in neuro-musculoskeletal simulations that are comparable to subtle gait changes found in prodromal SPG4 participants.

CONCLUSIONS

Predicting kinematic changes of prodromal and early-to-moderate manifest SPG4 participants by gradual alterations of sensory-motor reflex sensitivity allows us to link gait as a directly accessible performance marker to emerging neuro-muscular changes for early therapeutic interventions.

Quantitative assessment of spasticity: a narrative review of novel approaches and technologies

Spasticity is a complex neurological disorder, causing significant physical disabilities and affecting patients' independence and quality of daily lives. Current spasticity assessment methods are questioned for their non-standardized measurement protocols, limited reliabilities, and capabilities in distinguishing neuron or non-neuron factors in upper motor neuron lesion. A series of new approaches are developed for improving the effectiveness of current clinical used spasticity assessment methods with the developing technology in biosensors, robotics, medical imaging, biomechanics, telemedicine, and artificial intelligence. We investigated the reliabilities and effectiveness of current spasticity measures employed in clinical environments and the newly developed approaches, published from 2016 to date, which have the potential to be used in clinical environments. The new spasticity scales, taking advantage of quantified information such as torque, or echo intensity, the velocity-dependent feature and patients' self-reported information, grade spasticity semi-quantitatively, have competitive or better reliability than previous spasticity scales. Medical imaging technologies, including near-infrared spectroscopy, magnetic resonance imaging, ultrasound and thermography, can measure muscle hemodynamics and metabolism, muscle tissue properties, or temperature of tissue. Medical imaging-based methods are feasible to provide quantitative information in assessing and monitoring muscle spasticity. Portable devices, robotic based equipment or myotonometry, using information from angular, inertial, torque or surface EMG sensors, can quantify spasticity with the help of machine learning algorithms. However, spasticity measures using those devices are normally not physiological sound. Repetitive peripheral magnetic stimulation can assess patients with severe spasticity, which lost voluntary contractions. Neuromusculoskeletal modeling evaluates the neural and non-neural properties and may gain insights into the underlying pathology of spasticity muscles. Telemedicine technology enables outpatient spasticity assessment. The newly developed spasticity methods aim to standardize experimental protocols and outcome measures and enable quantified, accurate, and intelligent assessment. However, more work is needed to investigate and improve the effectiveness and accuracy of spasticity assessment.

Intra- and inter-subject variability of femoral growth plate stresses in typically developing children and children with cerebral palsy

Little is known about the influence of mechanical loading on growth plate stresses and femoral growth. A multi-scale workflow based on musculoskeletal simulations and mechanobiological finite element (FE) analysis can be used to estimate growth plate loading and femoral growth trends. Personalizing the model in this workflow is time-consuming and therefore previous studies included small sample sizes (N < 4) or generic finite element models. The aim of this study was to develop a semi-automated toolbox to perform this workflow and to quantify intra-subject variability in growth plate stresses in 13 typically developing (TD) children and 12 children with cerebral palsy (CP). Additionally, we investigated the influence of the musculoskeletal model and the chosen material properties on the simulation results. Intra-subject variability in growth plate stresses was higher in cerebral palsy than in typically developing children. The highest osteogenic index (OI) was observed in the posterior region in 62% of the TD femurs while in children with CP the lateral region was the most common (50%). A representative reference osteogenic index distribution heatmap generated from data of 26 TD children's femurs showed a ring shape with low values in the center region and high values at the border of the growth plate. Our simulation results can be used as reference values for further investigations. Furthermore, the code of the developed GP-Tool ("Growth Prediction-Tool") is freely available on GitHub (https://github.com/WilliKoller/GP-Tool) to enable peers to conduct mechanobiological growth studies with larger sample sizes to improve our understanding of femoral growth and to support clinical decision making in the near future.

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.

Multi-level personalization of neuromusculoskeletal models to estimate physiologically plausible knee joint contact forces in children

Neuromusculoskeletal models are a powerful tool to investigate the internal biomechanics of an individual. However, commonly used neuromusculoskeletal models are generated via linear scaling of generic templates derived from elderly adult anatomies and poorly represent a child, let alone children with a neuromuscular disorder whose musculoskeletal structures and muscle activation patterns are profoundly altered. Model personalization can capture abnormalities and appropriately describe the underlying (altered) biomechanics of an individual. In this work, we explored the effect of six different levels of neuromusculoskeletal model personalization on estimates of muscle forces and knee joint contact forces to tease out the importance of model personalization for normal and abnormal musculoskeletal structures and muscle activation patterns. For six children, with and without cerebral palsy, generic scaled models were developed and progressively personalized by (1) tuning and calibrating musculotendon units' parameters, (2) implementing an electromyogram-assisted approach to synthesize muscle activations, and (3) replacing generic anatomies with image-based bony geometries, and physiologically and physically plausible muscle kinematics. Biomechanical simulations of gait were performed in the OpenSim and CEINMS software on ten overground walking trials per participant. A mixed-ANOVA test, with Bonferroni corrections, was conducted to compare all models' estimates. The model with the highest level of personalization produced the most physiologically plausible estimates. Model personalization is crucial to produce physiologically plausible estimates of internal biomechanical quantities. In particular, personalization of musculoskeletal anatomy and muscle activation patterns had the largest effect overall. Increased research efforts are needed to ease the creation of personalized neuromusculoskeletal models.

Evaluation of musculoskeletal models, scaling methods, and performance criteria for estimating muscle excitations and fiber lengths across walking speeds

Muscle-driven simulations have been widely adopted to study muscle-tendon behavior; several generic musculoskeletal models have been developed, and their biofidelity improved based on available experimental data and computational feasibility. It is, however, not clear which, if any, of these models accurately estimate muscle-tendon dynamics over a range of walking speeds. In addition, the interaction between model selection, performance criteria to solve muscle redundancy, and approaches for scaling muscle-tendon properties remain unclear. This study aims to compare estimated muscle excitations and muscle fiber lengths, qualitatively and quantitatively, from several model combinations to experimental observations. We tested three generic models proposed by Hamner et al., Rajagopal et al., and Lai-Arnold et al. in combination with performance criteria based on minimization of muscle effort to the power of 2, 3, 5, and 10, and four approaches to scale the muscle-tendon unit properties of maximum isometric force, optimal fiber length, and tendon slack length. We collected motion analysis and electromyography data in eight able-bodied subjects walking at seven speeds and compared agreement between estimated/modelled muscle excitations and observed muscle excitations from electromyography and computed normalized fiber lengths to values reported in the literature. We found that best agreement in on/off timing in vastus lateralis, vastus medialis, tibialis anterior, gastrocnemius lateralis, gastrocnemius medialis, and soleus was estimated with minimum squared muscle effort than to higher exponents, regardless of model and scaling approach. Also, minimum squared or cubed muscle effort with only a subset of muscle-tendon unit scaling approaches produced the best time-series agreement and best estimates of the increment of muscle excitation magnitude across walking speeds. There were discrepancies in estimated fiber lengths and muscle excitations among the models, with the largest discrepancy in the Hamner et al. model. The model proposed by Lai-Arnold et al. best estimated muscle excitation estimates overall, but failed to estimate realistic muscle fiber lengths, which were better estimated with the model proposed by Rajagopal et al. No single model combination estimated the most accurate muscle excitations for all muscles; commonly observed disagreements include onset delay, underestimated co-activation, and failure to estimate muscle excitation increments across walking speeds.

Design of a control framework for lower limb exoskeleton rehabilitation robot based on predictive assessment

BACKGROUND

Patients suffering from lower limb dyskinesia, especially in early stages of rehabilitation, have weak residual muscle strength in affected limb and require passive training by the lower limb rehabilitation robot. Anatomy indicates that the biceps femoris short head muscle has a strong influence on knee motion at the swing phase of walking. We sought to explore how it would influence on gait cycle in optimization framework. However, the training trajectory of conventional rehabilitation robots performing passive training usually follows gait planning based on general human gait data, which cannot simultaneously ensure both effective rehabilitation of affected limbs with varying severity pathological gait and comfort of the wearer within a safe motion trajectory.

METHODS

To elucidate the effects of weakness and contracture, we systematically introduced isolated defects into the musculoskeletal model and generated walking simulations to predict the gait adaptation due to these defects. An impedance control model of the rehabilitation robot is developed. Knee joint parameters optimized by predictive forward dynamics simulation are adopted as the expected values for the robot controller to achieve customized adjustment of the robot motion trajectory.

FINDINGS

Severe muscle contracture leads to severe knee flexion; severe muscle weakness induces a significant posterior tilt of the upper trunk, which hinders walking speed.

INTERPRETATION

Our simulation results attempt to reveal pathological gait features, which may help to reproduce the simulation of pathological gait. Furthermore, the robot simulation results show that the robot system achieves a speedy tracking by setting a larger stiffness value. The model also allows the implementation of different levels of damping or elasticity effects.

TRIAL REGISTRATION

The method proposed in this paper is an initial basic study that did not reach clinical trials and therefore retains retrospectively registered.

Investigation of neural and biomechanical impairments leading to pathological toe and heel gaits using neuromusculoskeletal modelling

KEY POINTS

Pathological toe and heel gaits are commonly present in various conditions such as spinal cord injury, stroke or cerebral palsy. These conditions present various neural and biomechanical impairments and the cause-effect relationships between these impairments and pathological gaits are hard to establish clinically. Based on neuromechanical simulation, this study focuses on the plantarflexor muscles and builds a new reflex circuit controller to model and evaluate the potential effect of both neural and biomechanical impairments on gait. Our results suggest an important contribution of active reflex mechanisms in pathological toe gait. This "what if" based on neuromechanical modelling is thus deemed of great interest to target potential pathological gait causes.

ABSTRACT

This study investigates the pathological toe and heel gaits in human locomotion using neuromusculoskeletal modelling and simulation. In particular, it aims at investigating potential cause-effect relationships between biomechanical or neural impairments and pathological gaits. Toe and heel gaits are commonly present in spinal cord injury, stroke or cerebral palsy. Toe walking is mainly attributed to spasticity and contracture at plantarflexor muscles, whereas heel walking can be attributed to muscle weakness from biomechanical or neural origin. To investigate the effect of these impairments on gait, this study focuses on the soleus and gastrocnemius muscles as they contribute to ankle plantarflexion. We built a reflex circuit model on top of Geyer and Herr's work (2010) with additional pathways affecting the plantarflexor muscles. The SCONE software, which provides optimisation tools for 2D neuromechanical simulation of human locomotion, is used to optimise the corresponding reflex parameters and simulate healthy gait. We then modelled various bilateral plantarflexors biomechanical and neural impairments, and individually introduced them in the healthy model. We characterised the resulting simulated gaits as pathological or not by comparing ankle kinematics and ankle moment with the healthy optimised gait based on metrics used in clinical studies. Our simulations suggest that toe walking can be generated by hyperreflexia, whereas muscle and neural weaknesses induce partially heel gait. Thus, this "what if" approach is deemed of great interest as it allows the investigation of the effect of various impairments on gait and suggests an important contribution of active reflex mechanisms in pathological toe gait. Abstract figure legend Various biomechanical and neural impairments are individually modelled at the level of the plantarflexor muscles in a musculoskeletal model and a complex reflex circuit-based gait controller. For instance, as shown on the left, the plantarflexors spindle reflex gain (KS) is increased to mimic hyperreflexia. The gait controller is then optimised for each of the impaired condition and the resulting gaits are characterised as pathological gait based on ankle kinematics and ankle moment metrics used in clinical studies. Thus, this "what if" approach allows the investigation of the effect of various impairments on gait presented in the table on the right. This article is protected by copyright. All rights reserved.

The Contributions of Extracellular Matrix and Sarcomere Properties to Passive Muscle Stiffness in Cerebral Palsy

Cerebral palsy results from an upper motor neuron lesion and significantly affects skeletal muscle stiffness. The increased stiffness that occurs is partly a result of changes in the microstructural components of muscle. In particular, alterations in extracellular matrix, sarcomere length, fibre diameter, and fat content have been reported; however, experimental studies have shown wide variability in the degree of alteration. Many studies have reported changes in the extracellular matrix, while others have reported no differences. A consistent finding is increased sarcomere length in cerebral palsy affected muscle. Often many components are altered simultaneously, making it difficult to determine the individual effects on muscle stiffness. In this study, we use a three dimensional modelling approach to isolate individual effects of microstructural alterations typically occurring due to cerebral palsy on whole muscle behaviour; in particular, the effects of extracellular matrix volume fraction, stiffness, and sarcomere length. Causation between the changes to the microstructure and the overall muscle response is difficult to determine experimentally, since components of muscle cannot be manipulated individually; however, utilising a modelling approach allows greater control over each factor. We find that extracellular matrix volume fraction has the largest effect on whole muscle stiffness and mitigates effects from sarcomere length.

Alterations in the Functional Knee Alignment Are Not an Effective Strategy to Modify the Mediolateral Distribution of Knee Forces During Closed Kinetic Chain Exercises

Pain felt while performing rehabilitation exercises could be a reason for the low adherence of knee osteoarthritis patients to physical rehabilitation. Reducing compressive forces on the most affected knee regions may help to mitigate the pain. Knee frontal plane positioning with respect to pelvis and foot (functional knee alignment) has been shown to modify the mediolateral distribution of the tibiofemoral joint contact force in walking. Hence, different functional knee alignments could be potentially used to modify joint loading during rehabilitation exercises. The aim was to understand whether utilizing different alignments is an effective strategy to unload specific knee areas while performing rehabilitation exercises. Eight healthy volunteers performed 5 exercises with neutral, medial, and lateral knee alignment. A musculoskeletal model was modified for improved prediction of tibiofemoral contact forces and used to evaluate knee joint kinematics, moments, and contact forces. Functional knee alignment had only a small and inconsistent effect on the mediolateral distribution joint contact force. Moreover, the magnitude of tibiofemoral and patellofemoral contact forces, knee moments, and measured muscle activities was not significantly affected by the alignment. Our results suggest that altering the functional knee alignment is not an effective strategy to unload specific knee regions in physical rehabilitation.

In vivo assessment of the passive stretching response of the bi-compartmental human semitendinosus muscle using shear wave elastography

The semitendinosus muscle contains distinct proximal and distal compartments arranged anatomically in-series but separated by a tendinous inscription, with each compartment innervated by separate nerve branches. Although extensively investigated in other mammals, compartment-specific mechanical properties within the human semitendinosus have scarcely been assessed in vivo. Experimental data obtained during muscle-tendon unit stretching (e.g., slack angle) can also be used to validate and/or improve musculoskeletal model estimates of semitendinosus muscle force. The purpose of this study was to investigate the passive stretching response of proximal and distal humans semitendinosus compartments to distal joint extension. Using two-dimensional shear wave elastography, we bilaterally obtained shear moduli of both semitendinosus compartments from 14 prone-positioned individuals at ten knee flexion angles (from 90° to 0° [full extension] at 10° intervals). Passive muscle mechanical characteristics (slack angle, slack shear modulus, and the slope of the increase in shear modulus) were determined for each semitendinosus compartment by fitting a piecewise exponential model to the shear modulus-joint angle curves. We found no differences between compartments or legs for slack angle, slack shear modulus, or the slope of the increase in shear modulus. We also found the experimentally determined slack angle occurred at ~15-80° higher knee flexion angles compared to estimates from two commonly used musculoskeletal models, depending on participant and model used. Overall, these findings demonstrate that passive shear modulus-joint angle curves do not differ between proximal and distal human semitendinosus compartments, and provide experimental data to improve semitendinosus force estimates derived from musculoskeletal models.

Muscle actions on crossed and non-crossed joints during upright standing and gait: A comprehensive description based on induced acceleration analysis

The multibody nature of the musculoskeletal system makes each applied force potentially accelerate all body segments. Hence, muscles' actions on the kinematics of crossed and non-crossed joints should be estimated based on multibody dynamics. The objective of this study was to systematically investigate the actions of main lower limb muscles on the sagittal-plane angular kinematics of the hip, knee, and ankle joints, during upright standing and gait. Subject-specific simulations were performed to compute the muscle-tendon forces based on three-dimensional kinematic data collected from 10 able-bodied subjects during walking at preferred speed and during relaxed standing posture. A subject-scaled model consisting of the lower limb segments, 19 degrees of freedom and 92 Hill-type muscle-tendon units was used. Muscle-induced joint angular accelerations were estimated by Induced Acceleration Analysis in OpenSim. A comprehensive description of the estimated joint accelerations induced by lower limb muscles was presented, for upright standing and for the whole gait cycle. The observed muscle actions on crossed and non-crossed joints were phase- and task-specific. The main flexors and extensors for each joint were reported. Particular biarticular muscles presented actions opposite to their anatomical classification for specific joints. Antagonist muscle actions were revealed, such as the hitherto unknown opposite actions of the soleus and gastrocnemius at the ankle, and of the iliopsoas and soleus at the knee and ankle, during upright standing. Agonist actions among remote muscles were also identified. The presented muscle actions and their roles in joint kinematics of bipedal standing and walking contribute to understanding task-specific coordination.

The reliability of measuring medial gastrocnemius muscle-tendon unit lengths during gait

BACKGROUND

Ultrasound imaging combined with 3D motion analysis allows for in-vivo assessment of muscle-tendon unit lengths during gait. The clinical relevance of analysing the medial gastrocnemius (MG) and Achilles muscle-tendon junction (MTJ), MG mid-muscle belly fascicles (FAS) and muscle thickness was shown. However, their reliability error estimations over the gait cycle is unknown.

RESEARCH QUESTION

What are the intra- and inter-session errors associated with extracting MG belly, thickness, FAS and tendon lengths using ultrasound during gait in healthy participants?

METHODS

3D gait analysis was carried out in ten healthy adults as they walked on an instrumented treadmill at a comfortable walking speed. An ultrasound probe was secured on the leg and tracked by 3D motion analysis. Images were collected during two walking trials with the probe on the MG muscle-belly to estimate FAS lengths and muscle thickness, and during two trials with the probe on the MTJ to estimate muscle-belly and tendon lengths. A second session was performed after 5 ± 4 days where a different operator placed the ultrasound probe. The standard deviation (SD) of absolute and relative lengths changes during the gait cycle over different trials were calculated per participant. SD values averaged over participants represented intra- and inter-session errors.

RESULTS

For all assessed variables, the intra-session errors were <2.2 mm, except for the FAS lengths (3.1 mm). The inter-session errors were larger than the intra-session, with the highest values found for the absolute muscle-tendon unit lengths (5.6 mm). Relative length errors were smaller than absolute length errors.

SIGNIFICANCE

Intra-session errors, which may reflect natural variability and data processing errors, seem more critical when extracting absolute FAS than muscle-tendon lengths. Standardized probe positioning on the MTJ between sessions may improve the inter-session reliability. Expressing the lengths relative to their lengths as the beginning of the gait cycle reduces the inter-session errors.

Quantifying Muscle Forces and Joint Loading During Hip Exercises Performed With and Without an Elastic Resistance Band

An increase in hip joint contact forces (HJCFs) is one of the main contributing mechanical causes of hip joint pathologies, such as hip osteoarthritis, and its progression. The strengthening of the surrounding muscles of the joint is a way to increase joint stability, which results in the reduction of HJCF. Most of the exercise recommendations are based on expert opinions instead of evidence-based facts. This study aimed to quantify muscle forces and joint loading during rehabilitative exercises using an elastic resistance band (ERB). Hip exercise movements of 16 healthy volunteers were recorded using a three-dimensional motion capture system and two force plates. All exercises were performed without and with an ERB and two execution velocities. Hip joint kinematics, kinetics, muscle forces, and HJCF were calculated based on the musculoskeletal simulations in OpenSim. Time-normalized waveforms of the different exercise modalities were compared with each other and with reference values found during walking. The results showed that training with an ERB increases both target muscle forces and HJCF. Furthermore, the ERB reduced the hip joint range of motion during the exercises. The type of ERB used (soft vs. stiff ERB) and the execution velocity of the exercise had a minor impact on the peak muscle forces and HJCF. The velocity of exercise execution, however, had an influence on the total required muscle force. Performing hip exercises without an ERB resulted in similar or lower peak HJCF and lower muscle forces than those found during walking. Adding an ERB during hip exercises increased the peak muscle and HJCF but the values remained below those found during walking. Our workflow and findings can be used in conjunction with future studies to support exercise design.

Assessment of functional balance in children with sensory impairments undergoing hippotherapy

Objective: To evaluate the functional balance of children with sensory impairments submitted to hippotherapy. Methods: 24 children of both genres equally participated in three groups: cerebral palsy (CP), Down syndrome (DS), and intellectual disability (ID) with the respective age groups (10.71 ± 2.69 years), (12.83 ± 2.64 years), and (11 ± 1.69 years).There were 15 attendances in hippotherapy with riding materials specific to each group. The Pediatric Balance Scale (EEP) was used before the 1st and after the 15th hippotherapy session (pre and post moments) to assess functional balance. The data were analyzed using the Shapiro-Wilk tests (normality), Bartlett test (homogeneity), and, between the pre- and post-care times, the paired t-test (intra-groups) and ANOVA with Bonferroni's multiple comparison tests (between groups), with statistical significance for p < 0.05. Results: in the post-attendance, there was an increase in the EEP score for the three groups (intra groups) with significance for children with ID (p = 0.003) and DS (p = 0.033); the CP group had a lower score (inter groups) in both times, pre (p = 0.003) and post (p = 0.002) attendance. Conclusion: hippotherapy contributed to the functional balance of children with distinct sensory impairments, according to the clinical diagnosis and riding material specific to the group, thus being able to be considered a therapeutic method with relevant benefits regarding the sensory aspects of the population.

Impact of scaling errors of the thigh and shank segments on musculoskeletal simulation results

BACKGROUND

Musculoskeletal simulations are widely used in the research community. The locations of surface markers are mostly used to scale a generic model to the participant's anthropometry. Marker-based scaling approaches include errors due to inaccuracies in marker placements.

RESEARCH QUESTION

How do scaling errors of the thigh and shank segments influence simulation results?

METHODS

Motion capture data and magnetic resonance images from a child with cerebral palsy and a typically developing child were used to create a subject-specific reference model for each child. These reference models were modified to mimic scaling errors due to inaccurately placed lateral epicondyle markers, which are frequently used to scale the thigh and shank segments. The thigh length was altered in 1 % steps from the original length and the shank length was accordingly adjusted to keep the total leg length constant. Thirty additional models were created, which included models with an altered thigh length of ±15 %. Subsequently, musculoskeletal simulations with OpenSim were performed with all models. Joint kinematics, joint kinetics, muscle forces and joint contact forces (JCF) were compared between the reference and altered models.

RESULTS

The investigated scaling error influenced joint kinematics and joint kinetics by up to 9.4° (hip flexion angle) and 0.15 Nm/kg (knee flexion moment), respectively. Maximum muscle and JCF differences of 46 % (medial gastrocnemius) and 72 % (hip JCF) bodyweight, respectively, were observed between the reference and altered models. Scaling errors mainly changed the magnitude but not the shape of most analyzed parameters. The influence of scaling errors on simulation results were similar in both participants.

SIGNIFICANCE

Scaling errors of the thigh segment influence simulation results at all joints due to the global optimization approach used in musculoskeletal simulations. Our findings can be used to estimate potential errors due to marker-based scaling approaches in previous and future studies.

The Effects of Mechanical Scale on Neural Control and the Regulation of Joint Stability

Recent work has demonstrated how the size of an animal can affect neural control strategies, showing that passive viscoelastic limb properties have a significant role in determining limb movements in small animals but are less important in large animals. We extend that work to consider effects of mechanical scaling on the maintenance of joint integrity; i.e., the prevention of aberrant contact forces within joints that might lead to joint dislocation or cartilage degradation. We first performed a literature review to evaluate how properties of ligaments responsible for joint integrity scale with animal size. Although we found that the cross-sectional area of the anterior cruciate ligament generally scaled with animal size, as expected, the effects of scale on the ligament’s mechanical properties were less clear, suggesting potential adaptations in passive contributions to the maintenance of joint integrity across species. We then analyzed how the neural control of joint stability is altered by body scale. We show how neural control strategies change across mechanical scales, how this scaling is affected by passive muscle properties and the cost function used to specify muscle activations, and the consequences of scaling on internal joint contact forces. This work provides insights into how scale affects the regulation of joint integrity by both passive and active processes and provides directions for studies examining how this regulation might be accomplished by neural systems.

SENSITIVITY OF MODEL-PREDICTED MUSCLE FORCES OF PATIENTS WITH CEREBRAL PALSY TO VARIATIONS IN MUSCLE-TENDON PARAMETERS

Computational musculoskeletal modeling and simulation platforms are efficient tools to gain insight into the muscular coordination of patients with motor disabilities such as cerebral palsy (CP). Muscle force predictions from simulation programs are influenced by the architectural and contractile properties of muscle-tendon units. In this study, we aimed to evaluate the sensitivity of major lower limb muscle forces in patients with CP to changes in muscle-tendon parameters. Open-access datasets of children with CP ([Formula: see text]) and healthy children ([Formula: see text]) were considered. Monte Carlo analysis was executed to specify how sensitive the muscle forces to perturbations between [Formula: see text]% and [Formula: see text]% of the nominal value of the maximum isometric muscle force, optimal muscle fiber length, muscle pennation angle, tendon slack length, and maximum contraction velocity of muscle. The sensitivity analysis revealed that muscle forces of CP patients and healthy individuals were most sensitive to perturbations in the tendon slack length ([Formula: see text]), while forces of CP patients were more sensitive to tendon slack length when compared to the healthy group ([Formula: see text]). Muscle forces of patients and healthy individuals were insensitive to the other four parameters ([Formula: see text]), except for the gracilis and sartorius muscles in which the proportion of optimal muscle fiber length to tendon slack length is higher than 1; forces of these two muscles were also sensitive to the optimal muscle fiber length. The results of this study are expected to contribute to our understanding of which parameters should be personalized when conducting musculoskeletal modeling and simulation of patients with CP.

The effect of hip muscle weakness and femoral bony deformities on gait performance

BACKGROUND

Children with cerebral palsy (CP) present with a pathological gait pattern due to musculoskeletal impairments, such as muscle weakness and altered bony geometry. However, the effect of these impairments on gait performance is still unknown. Research aim:This study aimed to explore the effect of hip muscle weakness and femoral deformities on the gait performance of CP and typical developing (TD) subjects.

METHODS

6400 musculoskeletal models were created by weakening the hip extensors, abductors, adductors and flexors from 0% to 75 % and increasing the femoral anteversion angle (FAA) and neck shaft angle (NSA) from 20° to 60° and 120° to 160°, respectively. One TD and five CP gait patterns were imposed to each model and muscle forces were calculated. The effect of weakness and bony deformities on the capability gap (CG) at the hip, i.e. the lack in hip moment generating capacity to perform the gait pattern, was investigated using regression analysis.

RESULTS

The CG of apparent equinus, stiff knee gait, TD gait, jump gait and true equinus increased with 0.080, 0.038, 0.015, 0.023 and 0.005 Nm/kg per 10 percent hip abductor weakness increase, with 0.211, 0.130, 0.056, 0.045 and 0.011 Nm/kg per 10 degrees FAA increase and with 0.163, 0.080, 0.036, 0.043 and 0.011 Nm/kg per 10 degrees NSA increase, respectively. Combined weakness and bony deformities explained 96 %, 85 %, 82 %, 65 %, 40 % and 35 % of the variance in the CG of apparent equinus, TD gait, stiff knee gait, jump gait, true equinus and crouch gait, respectively.

SIGNIFICANCE

The results suggest that surgical correction of femoral deformities is more likely to be effective than strength training of hip muscles in enhancing CP gait performance. Jump gait, true equinus and especially crouch were more robust, while apparent equinus and stiff knee gait were limited by hip weakness and femoral deformities.

A multi-scale modelling framework combining musculoskeletal rigid-body simulations with adaptive finite element analyses, to evaluate the impact of femoral geometry on hip joint contact forces and femoral bone growth

Multi-scale simulations, combining muscle and joint contact force (JCF) from musculoskeletal simulations with adaptive mechanobiological finite element analysis, allow to estimate musculoskeletal loading and predict femoral growth in children. Generic linearly scaled musculoskeletal models are commonly used. This approach, however, neglects subject- and age-specific musculoskeletal geometry, e.g. femoral neck-shaft angle (NSA) and anteversion angle (AVA). This study aimed to evaluate the impact of proximal femoral geometry, i.e. altered NSA and AVA, on hip JCF and femoral growth simulations. Musculoskeletal models with NSA ranging from 120° to 150° and AVA ranging from 20° to 50° were created and used to calculate muscle and hip JCF based on the gait analysis data of a typically developing child. A finite element model of a paediatric femur was created from magnetic resonance images. The finite element model was morphed to the geometries of the different musculoskeletal models and used for mechanobiological finite element analysis to predict femoral growth trends. Our findings showed that hip JCF increase with increasing NSA and AVA. Furthermore, the orientation of the hip JCF followed the orientation of the femoral neck axis. Consequently, the osteogenic index, which is a function of cartilage stresses and defines the growth rate, barely changed with altered NSA and AVA. Nevertheless, growth predictions were sensitive to the femoral geometry due to changes in the predicted growth directions. Altered NSA had a bigger impact on the growth results than altered AVA. Growth simulations based on mechanobiological principles were in agreement with reported changes in paediatric populations.

Subject-specific muscle properties from diffusion tensor imaging significantly improve the accuracy of musculoskeletal models

Musculoskeletal modelling is an important platform on which to study the biomechanics of morphological structures in vertebrates and is widely used in clinical, zoological and palaeontological fields. The popularity of this approach stems from the potential to non-invasively quantify biologically important but difficult-to-measure functional parameters. However, while it is known that model predictions are highly sensitive to input values, it is standard practice to build models by combining musculoskeletal data from different sources resulting in 'generic' models for a given species. At present, there are little quantitative data on how merging disparate anatomical data in models impacts the accuracy of these functional predictions. This issue is addressed herein by quantifying the accuracy of both subject-specific human limb models containing individualised muscle force-generating properties and models built using generic properties from both elderly and young individuals, relative to experimental muscle torques obtained from an isokinetic dynamometer. The results show that subject-specific models predict isokinetic muscle torques to a greater degree of accuracy than generic models at the ankle (root-mean-squared error - 7.9% vs. 49.3% in elderly anatomy-based models), knee (13.2% vs. 57.3%) and hip (21.9% vs. 32.8%). These results have important implications for the choice of musculoskeletal properties in future modelling studies, and the relatively high level of accuracy achieved in the subject-specific models suggests that such models can potentially address questions about inter-subject variations of muscle functions. However, despite relatively high levels of overall accuracy, models built using averaged generic muscle architecture data from young, healthy individuals may lack the resolution and accuracy required to study such differences between individuals, at least in certain circumstances. The results do not wholly discourage the continued use of averaged generic data in musculoskeletal modelling studies but do emphasise the need for to maximise the accuracy of input values if studying intra-species form-function relationships in the musculoskeletal system.

Pre-operative gastrocnemius lengths in gait predict outcomes following gastrocnemius lengthening surgery in children with cerebral palsy

Equinus deformity is one of the most common gait deformities in children with cerebral palsy. We examined whether estimates of gastrocnemius length in gait could identify limbs likely to have short-term and long-term improvements in ankle kinematics following gastrocnemius lengthening surgery to correct equinus. We retrospectively analyzed data of 891 limbs that underwent a single-event multi-level surgery (SEMLS), and categorized outcomes based on the normalcy of ankle kinematics. Limbs with short gastrocnemius lengths that received a gastrocnemius lengthening surgery as part of a SEMLS (case limbs) were 2.2 times more likely than overtreated limbs (i.e., limbs who did not have short lengths, but still received a lengthening surgery) to have a good surgical outcome at the follow-up gait visit (good outcome rate of 71% vs. 33%). Case limbs were 1.2 times more likely than control limbs (i.e., limbs that had short gastrocnemius lengths but no lengthening surgery) to have a good outcome (71% vs. 59%). Three-fourths of the case limbs with a good outcome at the follow-up gait visit maintained this outcome over time, compared to only one-half of the overtreated limbs. Our results caution against over-prescription of gastrocnemius lengthening surgery and suggest gastrocnemius lengths can be used to identify good surgical candidates.

Single-event multilevel surgery, but not botulinum toxin injections normalize joint loading in cerebral palsy patients

BACKGROUND

Many patients with cerebral palsy present a pathologic gait pattern, which presumably induces aberrant musculoskeletal loading that interferes with natural bone growth, causing bone deformations on the long term. Botulinum toxin interventions and single-event multilevel surgeries are used to restore the gait pattern, assuming that a normal gait pattern restores musculoskeletal loading and thus prevents further bone deformation. However, it is unknown if these interventions are able to restore musculoskeletal loading. Hence, we investigated the impact of botulinum toxin injections and single-event multilevel surgery on musculoskeletal loading.

METHODS

Gait data collected in 93 children with bilateral cerebral palsy, which included pre- and post multi-level botulinum toxin (49 children) and single-event multilevel surgery (44 children) assessments, and 15 typically developing children were retrospectively processed using a musculoskeletal modelling workflow to calculate joint angles, moments, muscle and joint contact force magnitudes and orientations. Differences from the typically developing waveform were expressed by a root-mean square difference were compared using paired t-tests for each intervention separately (alpha <0.05).

FINDINGS

Botulinum toxin induced significant changes in the joint angles, but did not improve the muscle and joint contact forces. Single-event multilevel surgery induced significant kinematic and kinetic changes, which were associated with improved muscle and joint contact forces.

INTERPRETATION

The present results indicate that botulinum toxin injections were not able to restore normal gait kinematics nor musculoskeletal loading, whereas single-event multilevel surgery did successfully restore both. Therefore, single-event multilevel surgery might be protective against the re-occurrence of bone deformation on the longer term.

Physics-Based Simulations to Predict the Differential Effects of Motor Control and Musculoskeletal Deficits on Gait Dysfunction in Cerebral Palsy: A Retrospective Case Study

Physics-based simulations of walking have the theoretical potential to support clinical decision-making by predicting the functional outcome of treatments in terms of walking performance. Yet before using such simulations in clinical practice, their ability to identify the main treatment targets in specific patients needs to be demonstrated. In this study, we generated predictive simulations of walking with a medical imaging based neuro-musculoskeletal model of a child with cerebral palsy presenting crouch gait. We explored the influence of altered muscle-tendon properties, reduced neuromuscular control complexity, and spasticity on gait dysfunction in terms of joint kinematics, kinetics, muscle activity, and metabolic cost of transport. We modeled altered muscle-tendon properties by personalizing Hill-type muscle-tendon parameters based on data collected during functional movements, simpler neuromuscular control by reducing the number of independent muscle synergies, and spasticity through delayed muscle activity feedback from muscle force and force rate. Our simulations revealed that, in the presence of aberrant musculoskeletal geometries, altered muscle-tendon properties rather than reduced neuromuscular control complexity and spasticity were the primary cause of the crouch gait pattern observed for this child, which is in agreement with the clinical examination. These results suggest that muscle-tendon properties should be the primary target of interventions aiming to restore an upright gait pattern for this child. This suggestion is in line with the gait analysis following muscle-tendon property and bone deformity corrections. Future work should extend this single case analysis to more patients in order to validate the ability of our physics-based simulations to capture the gait patterns of individual patients pre- and post-treatment. Such validation would open the door for identifying targeted treatment strategies with the aim of designing optimized interventions for neuro-musculoskeletal disorders.

Muscle Synergy Constraints Do Not Improve Estimates of Muscle Activity From Static Optimization During Gait for Unimpaired Children or Children With Cerebral Palsy

Neuromusculoskeletal simulation provides a promising platform to inform the design of assistive devices or inform rehabilitation. For these applications, a simulation must be able to accurately represent the person of interest, such as an individual with a neurologic injury. If a simulation fails to predict how an individual recruits and coordinates their muscles during movement, it will have limited utility for informing design or rehabilitation. While inverse dynamic simulations have previously been used to evaluate anticipated responses from interventions, like orthopedic surgery or orthoses, they frequently struggle to accurately estimate muscle activations, even for tasks like walking. The simulated muscle activity often fails to represent experimentally measured muscle activity from electromyographic (EMG) recordings. Research has theorized that the nervous system may simplify the range of possible activations used during dynamic tasks, by constraining activations to weighted groups of muscles, referred to as muscle synergies. Synergies are altered after neurological injury, such as stroke or cerebral palsy (CP), and may provide a method for improving subject-specific models of neuromuscular control. The aim of this study was to test whether constraining simulation to synergies could improve estimated muscle activations compared to EMG data. We evaluated modeled muscle activations during gait for six typically developing (TD) children and six children with CP. Muscle activations were estimated with: (1) static optimization (SO), minimizing muscle activations squared, and (2) synergy SO (SynSO), minimizing synergy activations squared using the weights identified from EMG data for two to five synergies. While SynSO caused changes in estimated activations compared to SO, the correlation to EMG data was not higher in SynSO than SO for either TD or CP groups. The correlations to EMG were higher in CP than TD for both SO (CP: 0.48, TD: 0.36) and SynSO (CP: 0.46, TD: 0.26 for five synergies). Constraining activations to SynSO caused the simulated muscle stress to increase compared to SO for all individuals, causing a 157% increase with two synergies. These results suggest that constraining simulated activations in inverse dynamic simulations to subject-specific synergies alone may not improve estimation of muscle activations during gait for generic musculoskeletal models.

Comprimento de isquiotibiais, função motora grossa e marcha em crianças e adolescentes com paralisia cerebral

RESUMO O objetivo desse estudo é verificar se há relação entre o comprimento dos isquiotibiais, função motora grossa e marcha em crianças e adolescentes com paralisia cerebral (PC). Os participantes, entre 6 e 18 anos, foram classificados pelo Sistema de Classificação da Função Motora Grossa (GMFCS) nos níveis I, II e III através da escala Tardieu modificada, para avaliar o comprimento dos músculos isquiotibiais, sendo identificados em: R1 (primeira resistência da extensão da perna), R2 (segunda resistência da extensão da perna) e R2-R1 (diferença entre R1 e R2) do membro inferior esquerdo e direito. Para avaliar a função motora grossa, foi utilizado o teste gross motor function measure (GMFM-88), e a escala physicians rating scale modificada (PRS) foi utilizada para avaliar a marcha. Vinte e três participantes foram incluídos e os resultados evidenciaram correlações moderadas entre R1, R2 e PRS do membro inferior esquerdo e GMFM. As demais variáveis apresentaram uma correlação fraca.

Predicting gait adaptations due to ankle plantarflexor muscle weakness and contracture using physics-based musculoskeletal simulations

Deficits in the ankle plantarflexor muscles, such as weakness and contracture, occur commonly in conditions such as cerebral palsy, stroke, muscular dystrophy, Charcot-Marie-Tooth disease, and sarcopenia. While these deficits likely contribute to observed gait pathologies, determining cause-effect relationships is difficult due to the often co-occurring biomechanical and neural deficits. To elucidate the effects of weakness and contracture, we systematically introduced isolated deficits into a musculoskeletal model and generated simulations of walking to predict gait adaptations due to these deficits. We trained a planar model containing 9 degrees of freedom and 18 musculotendon actuators to walk using a custom optimization framework through which we imposed simple objectives, such as minimizing cost of transport while avoiding falling and injury, and maintaining head stability. We first generated gaits at prescribed speeds between 0.50 m/s and 2.00 m/s that reproduced experimentally observed kinematic, kinetic, and metabolic trends for walking. We then generated a gait at self-selected walking speed; quantitative comparisons between our simulation and experimental data for joint angles, joint moments, and ground reaction forces showed root-mean-squared errors of less than 1.6 standard deviations and normalized cross-correlations above 0.8 except for knee joint moment trajectories. Finally, we applied mild, moderate, and severe levels of muscle weakness or contracture to either the soleus (SOL) or gastrocnemius (GAS) or both of these major plantarflexors (PF) and retrained the model to walk at a self-selected speed. The model was robust to all deficits, finding a stable gait in all cases. Severe PF weakness caused the model to adopt a slower, "heel-walking" gait. Severe contracture of only SOL or both PF yielded similar results: the model adopted a "toe-walking" gait with excessive hip and knee flexion during stance. These results highlight how plantarflexor weakness and contracture may contribute to observed gait patterns.

SimCP: A Simulation Platform to Predict Gait Performance Following Orthopedic Intervention in Children With Cerebral Palsy

Gait deficits in cerebral palsy (CP) are often treated with a single-event multi-level surgery (SEMLS). Selecting the treatment options (combination of bony and soft tissue corrections) for a specific patient is a complex endeavor and very often treatment outcome is not satisfying. A deterioration in 22.8% of the parameters describing gait performance has been reported and there is need for additional surgery in 11% of the patients. Computational simulations based on musculoskeletal models that allow clinicians to test the effects of different treatment options before surgery have the potential to drastically improve treatment outcome. However, to date, no such simulation and modeling method is available. Two important challenges are the development of methods to include patient-specific neuromechanical impairments into the models and to simulate the effect of different surgical procedures on post-operative gait performance. Therefore, we developed the SimCP framework that allows the evaluation of the effect of different simulated surgeries on gait performance of a specific patient and includes a graphical user interface (GUI) that enables performing virtual surgery on the models. We demonstrated the potential of our framework for two case studies. Models reflecting the patient-specific musculoskeletal geometry and muscle properties are generated based solely on data collected before the treatment. The patient's motor control is described based on muscle synergies derived from pre-operative EMG. The GUI is then used to modify the musculoskeletal properties according to the surgical plan. Since SEMLS does not affect motor control, the same motor control model is used to define gait performance pre- and post-operative. We use the capability gap (CG), i.e., the difference between the joint moments needed to perform healthy walking and the joint moments the personalized model can generate, to quantify gait performance. In both cases, the CG was smaller post- then pre-operative and this was in accordance with the measured change in gait kinematics after treatment.

Computational modeling of neuromuscular response to swing-phase robotic knee extension assistance in cerebral palsy

Predicting subject-specific responses to exoskeleton assistance may aid in maximizing functional gait outcomes, such as achieving full knee-extension at foot contact in individuals with crouch gait from cerebral palsy (CP). The purpose of this study was to investigate the role of volitional and non-volitional muscle activity in subject-specific responses to knee extension assistance during walking with an exoskeleton. We developed a simulation framework to predict responses to exoskeleton torque by applying a stretch-reflex spasticity model with muscle excitations computed during unassisted walking. The framework was validated with data collected from six individuals with CP. Framework-predicted knee angle at terminal swing was within 4 ± 4° (mean ± sd) of the knee angle measured experimentally without the addition of spasticity. Kinematic responses in two-thirds of the participants could be accurately modeled using only underlying muscle activity and the applied exoskeleton torque; incorporating hamstring spasticity was necessary to recreate the measured kinematics to within 1 ± 1° in the remaining participants. We observed strong positive linear relationships between knee extension and exoskeleton assistance, and strong negative quadratic relationships between knee extension and spasticity. We utilized our framework to identify optimal torque profiles necessary to achieve full knee-extension at foot contact. An angular impulse of 0.061 ± 0.025 Nm·s·kg^-1·deg^-1 with 0.013 ± 0.002 Nm·kg^-1·deg^-1 of peak torque and 4.1 ± 1.9 W·kg^-1·deg^-1 peak mechanical power was required to achieve full knee extension (values normalized by knee excursion). This framework may aid the prescription of exoskeleton control strategies in pathologies with muscle spasticity. https://simtk.org/projects/knee-exo-pred/.

Neuromusculoskeletal Simulation Reveals Abnormal Rectus Femoris-Gluteus Medius Coupling in Post-stroke Gait

Post-stroke gait is often accompanied by muscle impairments that result in adaptations such as hip circumduction to compensate for lack of knee flexion. Our previous work robotically enhanced knee flexion in individuals post-stroke with Stiff-Knee Gait (SKG), however, this resulted in greater circumduction, suggesting the existence of abnormal coordination in SKG. The purpose of this work is to investigate two possible mechanisms of the abnormal coordination: (1) a reflex coupling between stretched quadriceps and abductors, and (2) a coupling between volitionally activated knee flexors and abductors. We used previously collected kinematic, kinetic and EMG measures from nine participants with chronic stroke and five healthy controls during walking with and without the applied knee flexion torque perturbations in the pre-swing phase of gait in the neuromusculoskeletal simulation. The measured muscle activity was supplemented by simulated muscle activations to estimate the muscle states of the quadriceps, hamstrings and hip abductors. We used linear mixed models to investigate two hypotheses: (H1) association between quadriceps and abductor activation during an involuntary period (reflex latency) following the perturbation and (H2) association between hamstrings and abductor activation after the perturbation was removed. We observed significantly higher rectus femoris (RF) activation in stroke participants compared to healthy controls within the involuntary response period following the perturbation based on both measured (H1, p < 0.001) and simulated (H1, p = 0.022) activity. Simulated RF and gluteus medius (GMed) activations were correlated only in those with SKG, which was significantly higher compared to healthy controls (H1, p = 0.030). There was no evidence of synergistic coupling between any combination of hamstrings and hip abductors (H2, p > 0.05) when the perturbation was removed. The RF-GMed coupling suggests an underlying abnormal coordination pattern in post-stroke SKG, likely reflexive in origin. These results challenge earlier assumptions that hip circumduction in stroke is simply a kinematic adaptation due to reduced toe clearance. Instead, abnormal coordination may underlie circumduction, illustrating the deleterious role of abnormal coordination in post-stroke gait.

A spasticity model based on feedback from muscle force explains muscle activity during passive stretches and gait in children with cerebral palsy

Muscle spasticity is characterized by exaggerated stretch reflexes and affects about 85% of the children with cerebral palsy. However, the mechanisms underlying spasticity and its influence on gait are not well understood. Here, we first aimed to model the response of spastic hamstrings and gastrocnemii in children with cerebral palsy to fast passive stretches. Then, we evaluated how the model applied to gait. We developed three models based on exaggerated proprioceptive feedback. The first model relied on feedback from muscle fiber length and velocity (velocity-related model), the second model relied on feedback from muscle fiber length, velocity, and acceleration (acceleration-related model), and the third model relied on feedback from muscle force and its first time derivative (force-related model). The force-related model better reproduced measured hamstrings and gastrocnemii activity during fast passive stretches (coefficients of determination (R2): 0.73 ± 0.10 and 0.60 ± 0.13, respectively, and root mean square errors (RMSE): 0.034 ± 0.031 and 0.009 ± 0.007, respectively) than the velocity-related model (R2: 0.46 ± 0.15 and 0.07 ± 0.13, and RMSE: 0.053 ± 0.051 and 0.015 ± 0.009), and the acceleration-related model (R2: 0.47 ± 0.15 and 0.09 ± 0.14, and RMSE: 0.052 ± 0.050 and 0.015 ± 0.008). Additionally, the force-related model predicted hamstrings and gastrocnemii activity that better correlated with measured activity during gait (cross correlations: 0.82 ± 0.09 and 0.85 ± 0.06, respectively) than the activity predicted by the velocity-related model (cross correlations: 0.49 ± 0.17 and 0.71 ± 0.22) and the acceleration-related model (cross correlations: 0.51 ± 0.16 and 0.67 ± 0.20). Our results therefore suggest that force encoding in muscle spindles in combination with altered feedback gains and thresholds underlie activity of spastic muscles during passive stretches and gait. Our model of spasticity opens new perspectives for studying movement impairments due to spasticity through simulation.

Skeletal Muscle Adaptations and Passive Muscle Stiffness in Cerebral Palsy: A Literature Review and Conceptual Model

This literature review focuses on the primary morphological and structural characteristics, and mechanical properties identified in muscles affected by spastic cerebral palsy (CP). CP is a non-progressive neurological disorder caused by brain damage and is commonly diagnosed at birth. Although the brain damage is not progressive, subsequent neuro-physiological developmental adaptations may initiate changes in muscle structure, function, and composition, causing abnormal muscle activity and coordination. The symptoms of CP vary among patients. However, muscle spasticity is commonly present and is one of the most debilitating effects of CP. Here, we present the current knowledge regarding the mechanical properties of skeletal tissue affected by spastic CP. An increase in sarcomere length, collagen content, and fascicle diameter, and a reduction in the number of satellite cells within spastic CP muscle were consistent findings in the literature. Studies differed, however, in changes in fascicle lengths and fiber diameters. We also present a conceptual mechanical model of fascicle force transmission that incorporates mechanisms that impact both serial and lateral force production, highlighting the connections between the macro and micro structures of muscle to assist in deducing specific mechanisms for property changes and reduced force production.

The influence of maximum isometric muscle force scaling on estimated muscle forces from musculoskeletal models of children with cerebral palsy

BACKGROUND

Musculoskeletal models do not include patient-specific muscle forces but rely on a scaled generic model, with muscle forces left unscaled in most cases. However, to use musculoskeletal simulations to inform clinical decision-making in children with cerebral palsy (CP), inclusion of subject-specific muscle forces is of utmost importance in order to represent each child's compensation mechanisms introduced through muscle weakness.

RESEARCH AIM

The aims of this study were to (i) evaluate if maximum isometric muscle forces (MIMF) in musculoskeletal models of children with CP can be scaled based on strength measurements obtained with a hand-held-dynamometer (HHD), (ii) evaluate the impact of the HHD based scaling approach and previously published MIMF scaling methods on computed muscle forces during gait, and (iii) compare maximum muscle forces during gait between CP and typically developing (TD) children.

METHODS

Strength and motion capture data of six CP and motion capture data of six TD children were collected. The HHD measurements to obtain hip, knee and ankle muscle strength were simulated in OpenSim and used to modify MIMF of the 2392-OpenSim model. These muscle forces were compared to the MIMF scaled on the child's body mass and a scaling approach, which included the body mass and muscle-tendon lengths. OpenSim was used to calculate peak muscle forces during gait.

RESULTS

Ankle muscle strength was insufficient to reproduce joint moments during walking when MIMF were scaled based on HHD. During gait, peak hip and knee extensor muscle forces were higher and peak ankle dorsi-flexor forces were lower in CP compared to TD participants.

SIGNIFICANCE

HHD measurements can be used to scale MIMF for the hip and knee muscle groups but underestimate the force capacity of the ankle muscle groups during walking. Muscle-tendon-length and mass based scaling methods affected muscle activations but had little influence on peak muscle forces during gait.

OpenSim: Simulating musculoskeletal dynamics and neuromuscular control to study human and animal movement

Movement is fundamental to human and animal life, emerging through interaction of complex neural, muscular, and skeletal systems. Study of movement draws from and contributes to diverse fields, including biology, neuroscience, mechanics, and robotics. OpenSim unites methods from these fields to create fast and accurate simulations of movement, enabling two fundamental tasks. First, the software can calculate variables that are difficult to measure experimentally, such as the forces generated by muscles and the stretch and recoil of tendons during movement. Second, OpenSim can predict novel movements from models of motor control, such as kinematic adaptations of human gait during loaded or inclined walking. Changes in musculoskeletal dynamics following surgery or due to human-device interaction can also be simulated; these simulations have played a vital role in several applications, including the design of implantable mechanical devices to improve human grasping in individuals with paralysis. OpenSim is an extensible and user-friendly software package built on decades of knowledge about computational modeling and simulation of biomechanical systems. OpenSim's design enables computational scientists to create new state-of-the-art software tools and empowers others to use these tools in research and clinical applications. OpenSim supports a large and growing community of biomechanics and rehabilitation researchers, facilitating exchange of models and simulations for reproducing and extending discoveries. Examples, tutorials, documentation, and an active user forum support this community. The OpenSim software is covered by the Apache License 2.0, which permits its use for any purpose including both nonprofit and commercial applications. The source code is freely and anonymously accessible on GitHub, where the community is welcomed to make contributions. Platform-specific installers of OpenSim include a GUI and are available on simtk.org.

Credibility, Replicability, and Reproducibility in Simulation for Biomedicine and Clinical Applications in Neuroscience

Modeling and simulation in computational neuroscience is currently a research enterprise to better understand neural systems. It is not yet directly applicable to the problems of patients with brain disease. To be used for clinical applications, there must not only be considerable progress in the field but also a concerted effort to use best practices in order to demonstrate model credibility to regulatory bodies, to clinics and hospitals, to doctors, and to patients. In doing this for neuroscience, we can learn lessons from long-standing practices in other areas of simulation (aircraft, computer chips), from software engineering, and from other biomedical disciplines. In this manuscript, we introduce some basic concepts that will be important in the development of credible clinical neuroscience models: reproducibility and replicability; verification and validation; model configuration; and procedures and processes for credible mechanistic multiscale modeling. We also discuss how garnering strong community involvement can promote model credibility. Finally, in addition to direct usage with patients, we note the potential for simulation usage in the area of Simulation-Based Medical Education, an area which to date has been primarily reliant on physical models (mannequins) and scenario-based simulations rather than on numerical simulations.

Spasticity Measurement Based on Tonic Stretch Reflex Threshold in Children with Cerebral Palsy Using the PediAnklebot

Nowadays, objective measures are becoming prominent in spasticity assessment, to overcome limitations of clinical scales. Among others, Tonic Stretch Reflex Threshold (TSRT) showed promising results. Previous studies demonstrated the validity and reliability of TSRT in spasticity assessment at elbow and ankle joints in adults. Purposes of the present study were to assess: (i) the feasibility of measuring TSRT to evaluate spasticity at the ankle joint in children with Cerebral Palsy (CP), and (ii) the correlation between objective measures and clinical scores. A mechatronic device, the pediAnklebot, was used to impose 50 passive stretches to the ankle of 10 children with CP and 3 healthy children, to elicit muscles response at 5 different velocities. Surface electromyography, angles, and angular velocities were recorded to compute dynamic stretch reflex threshold; TSRT was computed with a linear regression through angles and angular velocities. TSRTs for the most affected side of children with CP resulted into the biomechanical range (95.7 ± 12.9° and 86.7 ± 17.4° for Medial and Lateral Gastrocnemius, and 75.9 ± 12.5° for Tibialis Anterior). In three patients, the stretch reflex was not elicited in the less affected side. TSRTs were outside the biomechanical range in healthy children. However, no correlation was found between clinical scores and TSRT values. Here, we demonstrated the capability of TSRT to discriminate between spastic and non-spastic muscles, while no significant outcomes were found for the dorsiflexor muscle.

Robust Real-Time Musculoskeletal Modeling Driven by Electromyograms

OBJECTIVE

Current clinical biomechanics involves lengthy data acquisition and time-consuming offline analyses with biomechanical models not operating in real-time for man-machine interfacing. We developed a method that enables online analysis of neuromusculoskeletal function in vivo in the intact human.

METHODS

We used electromyography (EMG)-driven musculoskeletal modeling to simulate all transformations from muscle excitation onset (EMGs) to mechanical moment production around multiple lower-limb degrees of freedom (DOFs). We developed a calibration algorithm that enables adjusting musculoskeletal model parameters specifically to an individual's anthropometry and force-generating capacity. We incorporated the modeling paradigm into a computationally efficient, generic framework that can be interfaced in real-time with any movement data collection system.

RESULTS

The framework demonstrated the ability of computing forces in 13 lower-limb muscle-tendon units and resulting moments about three joint DOFs simultaneously in real-time. Remarkably, it was capable of extrapolating beyond calibration conditions, i.e., predicting accurate joint moments during six unseen tasks and one unseen DOF.

CONCLUSION

The proposed framework can dramatically reduce evaluation latency in current clinical biomechanics and open up new avenues for establishing prompt and personalized treatments, as well as for establishing natural interfaces between patients and rehabilitation systems.

SIGNIFICANCE

The integration of EMG with numerical modeling will enable simulating realistic neuromuscular strategies in conditions including muscular/orthopedic deficit, which could not be robustly simulated via pure modeling formulations. This will enable translation to clinical settings and development of healthcare technologies including real-time bio-feedback of internal mechanical forces and direct patient-machine interfacing.