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

Relation between stretch and activation of the medial gastrocnemius muscle during gait in children with cerebral palsy compared to typically developing children

Stretch hyperreflexia is often a target for treatment to improve gait in children with spastic cerebral palsy (CP). However, the presence of stretch hyperreflexia during gait remains debated. Therefore, we assessed the relation between gastrocnemius medialis muscle-tendon stretch and muscle activation during gait in children with CP compared to typically developing (TD) children. 3D gait analysis including electromyography (EMG) and dynamic ultrasound was carried out to assess, respectively gastrocnemius medialis activation and fascicle, belly, and tendon stretch during treadmill walking. Musculotendon-unit stretch was also estimated using OpenSim. Ratios of EMG/peak lengthening velocities and accelerations were compared between CP and TD. Velocity and acceleration peaks prior to EMG peaks were qualitatively assessed. EMG/velocity and EMG/acceleration ratios were up to 500% higher for CP (n = 14) than TD (n = 15) for most structures. Increased late swing muscle activation in CP was often preceded by fascicle and musculotendon-unit peak lengthening velocity, and early stance muscle activation by peaks in multiple structures. Increased muscle activation in CP is associated with muscle-tendon stretch during gait. Concluding, late swing muscle activation in CP appears velocity-dependent, whereas early stance activation can be velocity- and acceleration-dependent. These insights into stretch reflex mechanisms during gait can assist development of targeted interventions.

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.

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.

How Connecting the Legs with a Spring Improves Human Running Economy

Connecting the legs with a spring attached to the shoelaces reduces the energy cost of running, but how the spring reduces the energy burden of individual muscles remains unknown. We generated muscle-driven simulations of seven individuals running with and without the spring to discern whether savings occurred during the stance phase or the swing phase, and to identify which muscles contributed to energy savings. We computed differences in muscle-level energy consumption, muscle activations, and changes in muscle-fiber velocity and force between running with and without the spring. Across participants, running with the spring reduced the measured rate of energy expenditure by 0.9 W/kg (8.3%). Simulations predicted a 1.4 W/kg (12.0%) reduction in the average rate of energy expenditure and correctly identified that the spring reduced rates of energy expenditure for all participants. Simulations showed most of the savings occurred during stance (1.5 W/kg), though the rate of energy expenditure was also reduced during swing (0.3 W/kg). The energetic savings were distributed across the quadriceps, hip flexor, hip abductor, hamstring, hip adductor, and hip extensor muscle groups, whereas no changes in the rate of energy expenditure were observed in the plantarflexor or dorsiflexor muscles. Energetic savings were facilitated by reductions in the rate of mechanical work performed by muscles and their estimated rate of heat production. The simulations provide insight into muscle-level changes that occur when utilizing an assistive device and the mechanisms by which a spring connecting the legs improves running economy.

Robust walking control of a lower limb rehabilitation exoskeleton coupled with a musculoskeletal model via deep reinforcement learning

BACKGROUND

Few studies have systematically investigated robust controllers for lower limb rehabilitation exoskeletons (LLREs) that can safely and effectively assist users with a variety of neuromuscular disorders to walk with full autonomy. One of the key challenges for developing such a robust controller is to handle different degrees of uncertain human-exoskeleton interaction forces from the patients. Consequently, conventional walking controllers either are patient-condition specific or involve tuning of many control parameters, which could behave unreliably and even fail to maintain balance.

METHODS

We present a novel, deep neural network, reinforcement learning-based robust controller for a LLRE based on a decoupled offline human-exoskeleton simulation training with three independent networks, which aims to provide reliable walking assistance against various and uncertain human-exoskeleton interaction forces. The exoskeleton controller is driven by a neural network control policy that acts on a stream of the LLRE's proprioceptive signals, including joint kinematic states, and subsequently predicts real-time position control targets for the actuated joints. To handle uncertain human interaction forces, the control policy is trained intentionally with an integrated human musculoskeletal model and realistic human-exoskeleton interaction forces. Two other neural networks are connected with the control policy network to predict the interaction forces and muscle coordination. To further increase the robustness of the control policy to different human conditions, we employ domain randomization during training that includes not only randomization of exoskeleton dynamics properties but, more importantly, randomization of human muscle strength to simulate the variability of the patient's disability. Through this decoupled deep reinforcement learning framework, the trained controller of LLREs is able to provide reliable walking assistance to patients with different degrees of neuromuscular disorders without any control parameter tuning.

RESULTS AND CONCLUSION

A universal, RL-based walking controller is trained and virtually tested on a LLRE system to verify its effectiveness and robustness in assisting users with different disabilities such as passive muscles (quadriplegic), muscle weakness, or hemiplegic conditions without any control parameter tuning. Analysis of the RMSE for joint tracking, CoP-based stability, and gait symmetry shows the effectiveness of the controller. An ablation study also demonstrates the strong robustness of the control policy under large exoskeleton dynamic property ranges and various human-exoskeleton interaction forces. The decoupled network structure allows us to isolate the LLRE control policy network for testing and sim-to-real transfer since it uses only proprioception information of the LLRE (joint sensory state) as the input. Furthermore, the controller is shown to be able to handle different patient conditions without the need for patient-specific control parameter tuning.

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.

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.

Analysis of muscle activation in children and adolescents with severe cerebral palsy

Abstract Introduction: Children and adolescents with spastic quadriparesis have a worse selective motor control, and studies with this population are still very scarce. The same is true with scientific evidence of one of the methods most used as a physiotherapeutic treatment in this population, the Bobath Concept. Objective: To evaluate spine erector muscles activation, gluteus medius and gluteus maximus, through the handling of the Bobath Concept and the sustained kneeling posture in subjects with a diagnosis of severe cerebral palsy; and to compare muscle activation with a reference group, in order to increase the reliability of this study. Methods: A cross-sectional study was carried out with 38 children and adolescents with cerebral palsy, classified by GMFCS at levels IV and V, and 20 healthy participants, aged between 3 and 18 years. They were submitted to the handling of the Bobath Concept and to the sustained kneeling posture, with muscle activation obtained by electromyography. Results: We observed significant muscle activation during handling in side-sitting, with weight transfer and without the help of another therapist, and in the sustained kneeling posture, for the erector of the spine and gluteus medius. Conclusion: The evidence from this study suggests that both the handling in side-sitting and the sustained kneeling posture cause significant muscle activation in the erector of the spine and gluteus medius for severe quadriparesis subjects, GMFCS IV and V, which can contribute to the improvement of postural control and decision-making in physical therapy practice.

Disentangling acceleration-, velocity-, and duration-dependency of the short- and medium-latency stretch reflexes in the ankle plantarflexors

Motorized assessment of the stretch reflex is instrumental to gain understanding of the stretch reflex, its physiological origin and to differentiate effects of neurological disorders, like spasticity. Both short-latency (M1) and medium-latency (M2) stretch reflexes have been reported to depend on the velocity and acceleration of an applied ramp-and-hold perturbation. In the upper limb, M2 has also been reported to depend on stretch duration. However, wrong conclusions might have been drawn in previous studies as the interdependence of perturbation parameters (amplitude, duration, velocity, and acceleration) possibly created uncontrolled, confounding effects. We disentangled the duration-, velocity-, and acceleration-dependence and their interactions of the M1 and M2 stretch reflex in the ankle plantarflexors. To disentangle the parameter interdependence, 49 unique ramp-and-hold joint perturbations elicited reflexes in 10 healthy volunteers during a torque control task. Linear mixed model analysis showed that M1 depended on acceleration, not velocity or duration, whereas M2 depended on acceleration, velocity, and duration. Simulations of the muscle spindle Ia afferents coupled to a motoneuron pool corroborated these experimental findings. In addition, this simulation model did show a nonlinear M1 velocity- and duration-dependence for perturbation parameters outside the experimental scope. In conclusion, motorized assessment of the stretch reflex or spasticity using ramp-and-hold perturbations should be systematically executed and reported. Our systematic motorized and simulation assessments showed that M1 and M2 depend on acceleration, velocity, and duration of the applied perturbation. The simulation model suggested that these dependencies emerge from: muscle-tendon unit and muscle cross-bridge dynamics, Ia sensitivity to force and yank, and motoneuron synchronization.NEW & NOTEWORTHY Previous research and definitions of the stretch reflex and spasticity have focused on velocity-dependence. We showed that perturbation acceleration, velocity, and duration all shape the M1 and M2 response, often via nonlinear or interacting dependencies. Consequently, systematic execution and reporting of stretch reflex and spasticity studies, avoiding uncontrolled parameter interdependence, is essential for proper understanding of the reflex neurophysiology.

Neuromodulation Therapies for Spasticity Control: Now and Beyond

Spasticity is a major cause of disability following upper motor neuron (UMN) injury. The diagnosis and treatment of spasticity has been a focus of clinicians and researchers alike. In recent years, there have been significant advances both in strategies for spasticity assessment and in the development of novel treatments. Currently, several well-established spasticity management techniques fall into the major categories of physiotherapy, pharmacotherapy, and surgical management. The majority of recent developments in all of these broad categories have focused more on methods of neuromodulation instead of simple symptomatic treatment, attempting to address the underlying cause of spasticity more directly. The following narrative review briefly discusses the causes and clinical assessment of spasticity and also details the wide variety of current and developing treatment approaches for this often-debilitating condition.

Evaluation of force feedback in walking using joint torques as "naturalistic" stimuli

Control of adaptive walking requires the integration of sensory signals of muscle force and load. We have studied how mechanoreceptors (tibial campaniform sensilla) encode “naturalistic” stimuli derived from joint torques of stick insects walking on a horizontal substrate. Previous studies showed that forces applied to the legs using the mean torque profiles of a proximal joint were highly effective in eliciting motor activities. However, substantial variations in torque direction and magnitude occurred at the more distal femorotibial joint, which can generate braking or propulsive forces and provide lateral stability. To determine how these forces are encoded, we used torque waveforms of individual steps that had maximum values in stance in the directions of flexion or extension. Analysis of kinematic data showed that the torques in different directions tended to occur in different ranges of joint angles. Variations within stance were not accompanied by comparable changes in joint angle but often reflected vertical ground reaction forces and leg support of body load. Application of torque waveforms elicited sensory discharges with variations in firing frequency similar to those seen in freely walking insects. All sensilla directionally encoded the dynamics of force increases and showed hysteresis to transient force decreases. Smaller receptors exhibited more tonic firing. Our findings suggest that dynamic sensitivity in force feedback can modulate ongoing muscle activities to stabilize distal joints when large forces are generated at proximal joints. Furthermore, use of “naturalistic” stimuli can reproduce characteristics seen in freely moving animals that are absent in conventional restrained preparations.

NEW & NOTEWORTHY Sensory encoding of forces during walking by campaniform sensilla was characterized in stick insects using waveforms of joint torques calculated by inverse dynamics as mechanical stimuli. Tests using the mean joint torque and torques of individual steps showed the system is highly sensitive to force dynamics (dF/dt). Use of “naturalistic” stimuli can reproduce characteristics of sensory discharges seen in freely walking insects, such as load transfer among legs.

Perspective on musculoskeletal modelling and predictive simulations of human movement to assess the neuromechanics of gait

Locomotion results from complex interactions between the central nervous system and the musculoskeletal system with its many degrees of freedom and muscles. Gaining insight into how the properties of each subsystem shape human gait is challenging as experimental methods to manipulate and assess isolated subsystems are limited. Simulations that predict movement patterns based on a mathematical model of the neuro-musculoskeletal system without relying on experimental data can reveal principles of locomotion by elucidating cause-effect relationships. New computational approaches have enabled the use of such predictive simulations with complex neuro-musculoskeletal models. Here, we review recent advances in predictive simulations of human movement and how those simulations have been used to deepen our knowledge about the neuromechanics of gait. In addition, we give a perspective on challenges towards using predictive simulations to gain new fundamental insight into motor control of gait, and to help design personalized treatments in patients with neurological disorders and assistive devices that improve gait performance. Such applications will require more detailed neuro-musculoskeletal models and simulation approaches that take uncertainty into account, tools to efficiently personalize those models, and validation studies to demonstrate the ability of simulations to predict gait in novel circumstances.

Applying Stretch to Evoke Hyperreflexia in Spasticity Testing: Velocity vs. Acceleration

In neurological diseases, muscles often become hyper-resistant to stretch due to hyperreflexia, an exaggerated stretch reflex response that is considered to primarily depend on the muscle's stretch velocity. However, there is still limited understanding of how different biomechanical triggers applied during clinical tests evoke these reflex responses. We examined the effect of imposing a rotation with increasing velocity vs. increasing acceleration on triceps surae muscle repsonse in children with spastic paresis (SP) and compared the responses to those measured in typically developing (TD) children. A motor-operated ankle manipulator was used to apply different bell-shaped movement profiles, with three levels of maximum velocity (70, 110, and 150°/s) and three levels of maximum acceleration (500, 750, and 1,000°/s²). For each profile and both groups, we evaluated the amount of evoked triceps surae muscle activation. In SP, we evaluated two additional characteristics: the intensity of the response (peak EMG burst) and the time from movement initiation to onset of the EMG burst. As expected, the amount of evoked muscle activation was larger in SP compared to TD (all muscles: p < 0.001) and only sensitive to biomechanical triggers in SP. Further investigation of the responses in SP showed that peak EMG bursts increased in profiles with higher peak velocity (lateral gastrocnemius: p = 0.04), which was emphasized by fair correlations with increased velocity at EMG burst onset (all muscles: r > 0.33-0.36, p ≤ 0.008), but showed no significant effect for acceleration. However, the EMG burst was evoked faster with higher peak acceleration (all muscles p < 0.001) whereas it was delayed in profiles with higher peak velocity (medial gastrocnemius and soleus: p < 0.006). We conclude that while exaggerated response intensity (peak EMG burst) seems linked to stretch velocity, higher accelerations seem to evoke faster responses (time to EMG burst onset) in triceps surae muscles in SP. Understanding and controlling for the distinct effects of different biological triggers, including velocity, acceleration but also length and force of the applied movement, will contribute to the development of more precise clinical measurement tools. This is especially important when aiming to understand the role of hyperreflexia during functional movements where the biomechanical inputs are multiple and changing.

A Conceptual Blueprint for Making Neuromusculoskeletal Models Clinically Useful

The ultimate goal of most neuromusculoskeletal modeling research is to improve the treatment of movement impairments. However, even though neuromusculoskeletal models have become more realistic anatomically, physiologically, and neurologically over the past 25 years, they have yet to make a positive impact on the design of clinical treatments for movement impairments. Such impairments are caused by common conditions such as stroke, osteoarthritis, Parkinson’s disease, spinal cord injury, cerebral palsy, limb amputation, and even cancer. The lack of clinical impact is somewhat surprising given that comparable computational technology has transformed the design of airplanes, automobiles, and other commercial products over the same time period. This paper provides the author’s personal perspective for how neuromusculoskeletal models can become clinically useful. First, the paper motivates the potential value of neuromusculoskeletal models for clinical treatment design. Next, it highlights five challenges to achieving clinical utility and provides suggestions for how to overcome them. After that, it describes clinical, technical, collaboration, and practical needs that must be addressed for neuromusculoskeletal models to fulfill their clinical potential, along with recommendations for meeting them. Finally, it discusses how more complex modeling and experimental methods could enhance neuromusculoskeletal model fidelity, personalization, and utilization. The author hopes that these ideas will provide a conceptual blueprint that will help the neuromusculoskeletal modeling research community work toward clinical utility.

Neurophysiological validation of simultaneous intrinsic and reflexive joint impedance estimates

BACKGROUND

People with brain or neural injuries, such as cerebral palsy or spinal cord injury, commonly have joint hyper-resistance. Diagnosis and treatment of joint hyper-resistance is challenging due to a mix of tonic and phasic contributions. The parallel-cascade (PC) system identification technique offers a potential solution to disentangle the intrinsic (tonic) and reflexive (phasic) contributions to joint impedance, i.e. resistance. However, a simultaneous neurophysiological validation of both intrinsic and reflexive joint impedances is lacking. This simultaneous validation is important given the mix of tonic and phasic contributions to joint hyper-resistance. Therefore, the main goal of this paper is to perform a group-level neurophysiological validation of the PC system identification technique using electromyography (EMG) measurements.

METHODS

Ten healthy people participated in the study. Perturbations were applied to the ankle joint to elicit reflexes and allow for system identification. Participants completed 20 hold periods of 60 seconds, assumed to have constant joint impedance, with varying magnitudes of intrinsic and reflexive joint impedances across periods. Each hold period provided a paired data point between the PC-based estimates and neurophysiological measures, i.e. between intrinsic stiffness and background EMG, and between reflexive gain and reflex EMG.

RESULTS

The intrinsic paired data points, with all subjects combined, were strongly correlated, with a range of [Formula: see text] in both ankle plantarflexors and dorsiflexors. The reflexive paired data points were moderately correlated, with [Formula: see text] in the ankle plantarflexors only.

CONCLUSION

An agreement with the neurophysiological basis on which PC algorithms are built is necessary to support its clinical application in people with joint hyper-resistance. Our results show this agreement for the PC system identification technique on group-level. Consequently, these results show the validity of the use of the technique for the integrated assessment and training of people with joint hyper-resistance in clinical practice.

In Silico-Enhanced Treatment and Rehabilitation Planning for Patients with Musculoskeletal Disorders: Can Musculoskeletal Modelling and Dynamic Simulations Really Impact Current Clinical Practice?

Over the past decades, the use of computational physics-based models representative of the musculoskeletal (MSK) system has become increasingly popular in many fields of clinically driven research, locomotor rehabilitation in particular. These models have been applied to various functional impairments given their ability to estimate parameters which cannot be readily measured in vivo but are of interest to clinicians. The use of MSK modelling and simulations allows analysis of relevant MSK biomarkers such as muscle and joint contact loading at a number of different stages in the clinical treatment pathway in order to benefit patient functional outcome. Applications of these methods include optimisation of rehabilitation programs, patient stratification, disease characterisation, surgical pre-planning, and assistive device and exoskeleton design and optimisation. This review provides an overview of current approaches, the components of standard MSK models, applications, limitations, and assumptions of these modelling and simulation methods, and finally proposes a future direction.

Exploring the Contribution of Proprioceptive Reflexes to Balance Control in Perturbed Standing

Humans control balance using different feedback loops involving the vestibular system, the visual system, and proprioception. In this article, we focus on proprioception and explore the contribution of reflexes based on force and length feedback to standing balance. In particular, we address the questions of how much proprioception alone could explain balance control, and whether one modality, force or length feedback, is more important than the other. A sagittal plane neuro-musculoskeletal model was developed with six degrees of freedom and nine muscles in each leg. A controller was designed using proprioceptive reflexes and a dead zone. No feedback control was applied inside the dead zone. Reflexes were active once the center of mass moved outside the dead zone. Controller parameters were found by solving an optimization problem, where effort was minimized while the neuro-musculoskeletal model should remain standing upright on a perturbed platform. The ground was perturbed with random square pulses in the sagittal plane with different amplitudes and durations. The optimization was solved for three controllers: using force and length feedback (base model), using only force feedback, and using only length feedback. Simulations were compared to human data from previous work, where an experiment with the same perturbation signal was performed. The optimized controller yielded a similar posture, since average joint angles were within 5 degrees of the experimental average joint angles. The joint angles of the base model, the length only model, and the force only model correlated weakly (ankle) to moderately with the experimental joint angles. The ankle moment correlated weakly to moderately with the experimental ankle moment, while the hip and knee moment were only weakly correlated, or not at all. The time series of the joint angles showed that the length feedback model was better able to explain the experimental joint angles than the force feedback model. Changes in time delay affected the correlation of the joint angles and joint moments. The objective of effort minimization yielded lower joint moments than in the experiment, suggesting that other objectives are also important in balance control, which cause an increase in effort and thus larger joint moments.

Quantitative Modeling of Spasticity for Clinical Assessment, Treatment and Rehabilitation

Spasticity, a common symptom in patients with upper motor neuron lesions, reduces the ability of a person to freely move their limbs by generating unwanted reflexes. Spasticity can interfere with rehabilitation programs and cause pain, muscle atrophy and musculoskeletal deformities. Despite its prevalence, it is not commonly understood. Widely used clinical scores are neither accurate nor reliable for spasticity assessment and follow up of treatments. Advancement of wearable sensors, signal processing and robotic platforms have enabled new developments and modeling approaches to better quantify spasticity. In this paper, we review quantitative modeling techniques that have been used for evaluating spasticity. These models generate objective measures to assess spasticity and use different approaches, such as purely mechanical modeling, musculoskeletal and neurological modeling, and threshold control-based modeling. We compare their advantages and limitations and discuss the recommendations for future studies. Finally, we discuss the focus on treatment and rehabilitation and the need for further investigation in those directions.

Movement History Influences Pendulum Test Kinematics in Children With Spastic Cerebral Palsy

The pendulum test assesses quadriceps spasticity by dropping the lower leg of a relaxed patient from the horizontal position and observing limb movement. The first swing excursion (FS) decreases with increasing spasticity severity. Our recent simulation study suggests that the reduced initial swing results from muscle short-range stiffness and its interaction with reflex hyper-excitability. Short-range stiffness emerges from the thixotropic behavior of muscles where fiber stiffness upon stretch increases when the muscle is held isometric. Fiber stiffness might thus be higher during the first swing of the pendulum test than during consecutive swings. In addition, it has recently been suggested that muscle spindle firing reflects fiber force rather than velocity and therefore, reflex activity might depend on fiber stiffness. If this hypothesized mechanism is true, we expect to observe larger first swing excursions and reduced reflex muscle activity when the leg is moved rather than kept isometric before release, especially in patients with increased reflex activity. We performed the pendulum test in 15 children with cerebral palsy (CP) and 15 age-matched typically developing (TD) children in two conditions. In the hold condition, the leg was kept isometric in the extended position before release. In the movement condition, the leg was moved up and down before release to reduce the contribution of short-range stiffness. Knee kinematics and muscle activity were recorded. Moving the leg before release increased first swing excursion (p < 0.001) and this increase was larger in children with CP (21°) than in TD children (8°) (p < 0.005). In addition, pre-movement delayed reflex onset by 87 ms (p < 0.05) and reduced reflex activity as assessed through the area under the curve of rectus femoris electromyography (p < 0.05) in children with CP. The movement history dependence of pendulum kinematics and reflex activity supports our hypothesis that muscle short-range stiffness and its interaction with reflex hyper-excitability contribute to joint hyper-resistance in spastic CP. Our results have implications for standardizing movement history in clinical tests of spasticity and for understanding the role of spasticity in functional movements, where movement history differs from movement history in clinical tests.

The Interzonal Distribution of the Load on the Plantar Surface of the Foot During Walking in the Patients with Cerebral Palsy as an Objective Criterion of Functional Impairment Severity

Relevance. The main direction of rehabilitation of children with cerebral palsy is the preservation and enhancement of the existing level of support and locomotion, as well as compensation of its impairment through various methods of rehabilitation. For an adequate prescription and reliable assessment of these measures effectiveness, it is necessary to use objective indicators of functional impairment characteristic of cerebral palsy. The purpose of this study was to substantiate objective biomechanical indicators of functional impairment in children with cerebral palsy based on the analysis of the interzonal distribution of the load on the foot during walking, taking into account the level of global motor functions impairment. Materials and Methods. 47 children with cerebral palsy at the GMFCS levels of impairment 1 to 3 were examined. The control group consisted of 14 children without anatomical and functional signs of support and locomotion system impairment. Biomechanical examination was performed on the complex «DiaSled-M-Scan» with matrix plantar pressure meters in the form of insoles. The statistical analysis of the data was carried out by nonparametric methods using the SPSS for Widows software. Results. The analysis of the anatomical and functional impairment of 94 feet of the children with cerebral palsy and 28 feet of the control group revealed differences in the interzonal distribution of the load under the feet in six variables (p from 0.001 to 0.003). The most typical were: an increase in the toe-to-heel load ratio (on average by 80%), an increase in the load on the arch (by 49%), and a decrease in the medio-lateral load ratio on the toe (by 37%). For GMFCS 1 patients, a significant indicator of impairment was an increase in the partial load on the arch, for GMFCS 2 and 3 patients — a decrease in the load on the heel and an increase it under the toe. This leads to an increase in the toe-to-heel load ratio. Conclusion. It is advisable to use the revealed indicators of roll-over-the-foot impairment in the functional diagnosis of the condition and in assessing the effectiveness of rehabilitation of children with cerebral palsy.

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.

Non-ambulatory children with cerebral palsy: effects of four months of static and dynamic standing exercise on passive range of motion and spasticity in the hip

Purpose

The aim of this study was to compare the effects of four months of two types of structured training regimes, static standing (StS) versus dynamic standing (DyS), on passive range of motion (PROM) and spasticity in the hip among non-ambulatory children with cerebral palsy.

Method

Twenty non-ambulatory children with cerebral palsy participated in an exercise intervention study with a crossover design. During StS, the Non-ambulatory children with cerebral palsy were encouraged to exercise according to standard care recommendations, including daily supported StS for 30-90 min. During DyS, daily exercise for at least 30 min at a speed between 30 and 50 rpm in an Innowalk (Made for movement, Norway) was recommended. We assessed adaptive effects from the exercise programs through PROM in the hip assessed with a handheld goniometer, and spasticity in the hip assessed with the Modified Ashworth Scale before and after 30 min of StS or DyS. A trained physiotherapist performed the assessments. The exercise test and exercise training were performed in the children's habitual environment. Non-parametric statistics were used and each leg was used as its own control.

Result

PROM increased in all directions after 30 min (p < 0.001), and after four months of exercise training (p < 0.001) of DyS. Thirty minutes of DyS lowered the spasticity in the muscles around the hip (p < 0.001) more than 30 min of StS (p < 0.001).

Conclusion

Thirty minutes of DyS increased PROM and decreased spasticity among non-ambulatory children with CP. Four months of DyS increased PROM but did not decrease spasticity. These results can help inform individualised standing recommendations.

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.

Pre-treatment EMG can be used to model post-treatment muscle coordination during walking in children with cerebral palsy

When treating children with Cerebral Palsy (CP), computational simulations based on musculoskeletal models have a great potential in assisting the clinical decision-making process towards the most promising treatments. In particular, predictive simulations could be used to predict and compare the functional outcome of a series of candidate interventions. In order to be able to benefit from these predictive simulations however, it is important to know how much information about the post-treatment patient’s motor control could be gathered from data available before the intervention. Within this paper, we quantified how much of the muscle activity measured after a treatment could be explained by subject-specific muscle synergies computed from EMG data collected before the intervention. We also investigated whether generic synergies could be used, in case no EMG data is available when running predictive simulations, to reproduce both pre- and post-treatment muscle activity in children with CP. Subject-specific synergies proved to be a good indicator of the patient’s post-treatment motor control, explaining on average more than 85% of the post-treatment muscle activity, compared to an average of 94% when applied to the original pre-treatment data. Generic synergies explained 84% of the pre-treatment and 83% of the post-treatment muscle activity on average, but performed relatively well for patients with low selective motor control and poorly in patients with more selectivity. Our results suggest that subject-specific muscle synergies computed from pre-treatment EMG data could be used with confidence to represent the post-treatment motor control of children with CP during walking. In addition, when performing simulations involving patients with a low selective motor control, generic synergies could be a valid alternative.

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.