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, and Charcot-Marie-Tooth disease. While these deficits likely contribute to observed gait pathologies, elucidating 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 developed a planar model containing 9 degrees of freedom and 18 musculotendon actuators, and an 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 validated that our model could generate gaits that reproduced experimentally observed kinematic, kinetic, and metabolic trends for two cases: 1) walking at prescribed speeds between 0.50 m/s and 2.00 m/s and 2) walking at self-selected speed. We then 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. Our model, simulation and optimization framework, and results are freely shared so that others can reproduce and build upon our work. Author summary Deficits in the muscles that extend the ankle are thought to contribute to abnormal walking patterns in conditions such as cerebral palsy, stroke, muscular dystrophy, and Charcot-Marie-Tooth disease. To study how deficits in these muscles contribute to abnormal walking patterns, we used computer simulations to systematically introduce muscle deficits into a biomechanically accurate model. We first showed that our model could discover realistic walking patterns over a wide range of speeds when we posed a simple objective: walking while consuming a minimum amount of energy per distance, maintaining head stability, and avoiding injury. We then used the model to study the effect of two commonly observed problems: muscle weakness and muscle tightness. We found that severe weakness of the ankle extensors caused the model to adopt a slower, “heel-walking” gait, and severe tightness caused the model to adopt a crouched, “toe-walking” gait. These results highlight how deficits in the ankle extensor muscles may contribute to abnormal walking patterns commonly seen in pathological populations.

[1]  Michael A Sherman,et al.  Simbody: multibody dynamics for biomedical research. , 2011, Procedia IUTAM.

[2]  G E Fulford,et al.  Surgical management of ankle and foot deformities in cerebral palsy. , 1990, Clinical orthopaedics and related research.

[3]  Samuel R Ward,et al.  High resolution muscle measurements provide insights into equinus contractures in patients with cerebral palsy , 2015, Journal of orthopaedic research : official publication of the Orthopaedic Research Society.

[4]  B. Bril,et al.  Head Coordination as a Means to Assist Sensory Integration in Learning to Walk , 1998, Neuroscience & Biobehavioral Reviews.

[5]  J Koreska,et al.  MUSCLE GROWTH IN NORMAL AND SPASTIC MICE , 1984, Developmental medicine and child neurology.

[6]  K. Desloovere,et al.  Neuro-musculoskeletal simulation of instrumented contracture and spasticity assessment in children with cerebral palsy , 2016, Journal of NeuroEngineering and Rehabilitation.

[7]  T. Wren,et al.  Cross-correlation as a method for comparing dynamic electromyography signals during gait. , 2006, Journal of biomechanics.

[8]  Scott L. Delp,et al.  Full-Body Musculoskeletal Model for Muscle-Driven Simulation of Human Gait , 2016, IEEE Transactions on Biomedical Engineering.

[9]  Scott L. Delp,et al.  Predictive Simulation Generates Human Adaptations during Loaded and Inclined Walking , 2015, PloS one.

[10]  S. Delp,et al.  Muscle contributions to support and progression during single-limb stance in crouch gait. , 2010, Journal of biomechanics.

[11]  Seungmoon Song,et al.  Predictive neuromechanical simulations indicate why walking performance declines with ageing , 2018, The Journal of physiology.

[12]  Joseph Hamill,et al.  Evaluation of the minimum energy hypothesis and other potential optimality criteria for human running , 2012, Proceedings of the Royal Society B: Biological Sciences.

[13]  C. A. Byrne,et al.  Is rectus femoris really a part of quadriceps? Assessment of rectus femoris function during gait in able-bodied adults. , 2004, Gait & posture.

[14]  M. Clancy,et al.  A Geometric Method of Calculating Tendo Achillis Lengthening , 1985, Journal of pediatric orthopedics.

[15]  Thomas Geijtenbeek,et al.  SCONE: Open Source Software for Predictive Simulation of Biological Motion , 2019, J. Open Source Softw..

[16]  S. Delp,et al.  Preserving plantar flexion strength after surgical treatment for contracture of the triceps surae: A computer simulation study , 1995, Journal of orthopaedic research : official publication of the Orthopaedic Research Society.

[17]  Ayman Habib,et al.  OpenSim: Open-Source Software to Create and Analyze Dynamic Simulations of Movement , 2007, IEEE Transactions on Biomedical Engineering.

[18]  Nikolaus Hansen,et al.  Evaluating the CMA Evolution Strategy on Multimodal Test Functions , 2004, PPSN.

[19]  J. Perry,et al.  Gait Analysis , 2024 .

[20]  T. Wren,et al.  Prevalence of Specific Gait Abnormalities in Children With Cerebral Palsy: Influence of Cerebral Palsy Subtype, Age, and Previous Surgery , 2005, Journal of pediatric orthopedics.

[21]  D. Thelen Adjustment of muscle mechanics model parameters to simulate dynamic contractions in older adults. , 2003, Journal of biomechanical engineering.

[22]  Philip E. Martin,et al.  Effects of age and physical activity status on the speed-aerobic demand relationship of walking. , 1992, Journal of applied physiology.

[23]  Matthew S. DeMers,et al.  Changes in tibiofemoral forces due to variations in muscle activity during walking , 2014, Journal of orthopaedic research : official publication of the Orthopaedic Research Society.

[24]  Hartmut Geyer,et al.  A Muscle-Reflex Model That Encodes Principles of Legged Mechanics Produces Human Walking Dynamics and Muscle Activities , 2010, IEEE Transactions on Neural Systems and Rehabilitation Engineering.

[25]  C. Frijns,et al.  Normal values of patellar and ankle tendon reflex latencies , 1997, Clinical Neurology and Neurosurgery.

[26]  G. Taga Emergence of Locomotion , 1997 .

[27]  M. Abel,et al.  Strategies for increasing walking speed in diplegic cerebral palsy. , 1996, Journal of pediatric orthopedics.

[28]  Rory O'Sullivan,et al.  The characteristics of gait in Charcot-Marie-Tooth disease types I and II. , 2007, Gait & posture.

[29]  Lucas R. Smith,et al.  Hamstring contractures in children with spastic cerebral palsy result from a stiffer extracellular matrix and increased in vivo sarcomere length , 2011, The Journal of physiology.

[30]  Seungmoon Song,et al.  A neural circuitry that emphasizes spinal feedback generates diverse behaviours of human locomotion , 2015, The Journal of physiology.

[31]  H. Ralston Energy-speed relation and optimal speed during level walking , 1958, Internationale Zeitschrift für angewandte Physiologie einschließlich Arbeitsphysiologie.

[32]  R. Fitzpatrick,et al.  Acceleration patterns of the head and pelvis when walking on level and irregular surfaces. , 2003, Gait & posture.

[33]  S. Delp,et al.  Changes in sarcomere lengths of the human vastus lateralis muscle with knee flexion measured using in vivo microendoscopy. , 2016, Journal of biomechanics.

[34]  J. Eng,et al.  The relationship of lower-extremity muscle torque to locomotor performance in people with stroke. , 2003, Physical therapy.

[35]  Rod Barrett,et al.  Passive muscle mechanical properties of the medial gastrocnemius in young adults with spastic cerebral palsy. , 2011, Journal of biomechanics.

[36]  J. Gage,et al.  Surgical treatment of knee dysfunction in cerebral palsy. , 1990, Clinical orthopaedics and related research.

[37]  Ray-Yau Wang,et al.  The relation between ankle impairments and gait velocity and symmetry in people with stroke. , 2006, Archives of physical medicine and rehabilitation.

[38]  Timothy J Doherty,et al.  Sarcopenia: prevalence, mechanisms, and functional consequences. , 2010, Interdisciplinary topics in gerontology.

[39]  Jeremy D Wong,et al.  Control of position and movement is simplified by combined muscle spindle and Golgi tendon organ feedback. , 2013, Journal of neurophysiology.

[40]  Philip E. Martin,et al.  A Model of Human Muscle Energy Expenditure , 2003, Computer methods in biomechanics and biomedical engineering.

[41]  V. Scaioli,et al.  Clinical and electrophysiological aspects of Charcot-Marie-Tooth disease , 2006 .

[42]  F. Zajac,et al.  Contributions of the individual ankle plantar flexors to support, forward progression and swing initiation during walking. , 2001, Journal of biomechanics.

[43]  Scott L. Delp,et al.  A Model of the Lower Limb for Analysis of Human Movement , 2010, Annals of Biomedical Engineering.

[44]  C. Meyer,et al.  Relationships of 35 lower limb muscles to height and body mass quantified using MRI. , 2014, Journal of biomechanics.

[45]  P. Mermelstein,et al.  Opposite Effects of mGluR1a and mGluR5 Activation on Nucleus Accumbens Medium Spiny Neuron Dendritic Spine Density , 2016, PloS one.

[46]  J. Gage,et al.  An update on the treatment of gait problems in cerebral palsy. , 2001, Journal of pediatric orthopedics. Part B.

[47]  T. Wren,et al.  Gastrocnemius and soleus lengths in cerebral palsy equinus gait--differences between children with and without static contracture and effects of gastrocnemius recession. , 2004, Journal of biomechanics.

[48]  C. Meyer,et al.  Heterogeneity of muscle sizes in the lower limbs of children with cerebral palsy , 2016, Muscle & nerve.

[49]  J. Crosbie,et al.  Rest length and compliance of non-immobilised and immobilised rabbit soleus muscle and tendon , 1997, European Journal of Applied Physiology and Occupational Physiology.

[50]  T. Hortobágyi,et al.  The biomechanical mechanism of how strength and power training improves walking speed in old adults remains unknown , 2013, Ageing Research Reviews.

[51]  Fan Gao,et al.  Changes in passive mechanical properties of the gastrocnemius muscle at the muscle fascicle and joint levels in stroke survivors. , 2009, Archives of physical medicine and rehabilitation.

[52]  Mariano Serrao,et al.  Foot drop and plantar flexion failure determine different gait strategies in Charcot-Marie-Tooth patients. , 2007, Clinical biomechanics.

[53]  S. Delp,et al.  Rectus femoris and vastus intermedius fiber excursions predicted by three-dimensional muscle models. , 2006, Journal of biomechanics.

[54]  Scott L. Delp,et al.  Simulation-Based Design for Wearable Robotic Systems: An Optimization Framework for Enhancing a Standing Long Jump , 2016, IEEE Transactions on Biomedical Engineering.

[55]  Vladlen Koltun,et al.  Optimizing locomotion controllers using biologically-based actuators and objectives , 2012, ACM Trans. Graph..

[56]  Matthew Millard,et al.  Flexing computational muscle: modeling and simulation of musculotendon dynamics. , 2013, Journal of biomechanical engineering.

[57]  Sebastian I Wolf,et al.  Long-term results after gastrocnemius-soleus intramuscular aponeurotic recession as a part of multilevel surgery in spastic diplegic cerebral palsy. , 2012, The Journal of bone and joint surgery. American volume.

[58]  Marko Ackermann,et al.  Optimality principles for model-based prediction of human gait. , 2010, Journal of biomechanics.

[59]  D. Sutherland,et al.  The role of the ankle plantar flexors in normal walking. , 1980, The Journal of bone and joint surgery. American volume.

[60]  E Bakker,et al.  Signs and symptoms of Duchenne muscular dystrophy and Becker muscular dystrophy among carriers in the Netherlands: a cohort study , 1999, The Lancet.

[61]  May Q. Liu,et al.  Muscle contributions to support and progression over a range of walking speeds. , 2008, Journal of biomechanics.

[62]  Jessica C. Selinger,et al.  Humans Can Continuously Optimize Energetic Cost during Walking , 2015, Current Biology.

[63]  D Lyttle,et al.  Pathogenesis of pes cavus in Charcot-Marie-Tooth disease. , 1983, Clinical orthopaedics and related research.

[64]  R. Kram,et al.  Effects of obesity and sex on the energetic cost and preferred speed of walking. , 2006, Journal of applied physiology.

[65]  F.E. Zajac,et al.  An interactive graphics-based model of the lower extremity to study orthopaedic surgical procedures , 1990, IEEE Transactions on Biomedical Engineering.

[66]  J. Eng,et al.  Magnitude and pattern of 3D kinematic and kinetic gait profiles in persons with stroke: relationship to walking speed. , 2004, Gait & posture.

[67]  Thomas M Best,et al.  Simulation of biceps femoris musculotendon mechanics during the swing phase of sprinting. , 2005, Medicine and science in sports and exercise.

[68]  Michael H Schwartz,et al.  The effect of walking speed on the gait of typically developing children. , 2008, Journal of biomechanics.

[69]  Manoj Srinivasan,et al.  Robotic lower limb prosthesis design through simultaneous computer optimizations of human and prosthesis costs , 2016, Scientific reports.

[70]  S. Thomas,et al.  Calcaneal Gait in spastic Diplegia After Heel Cord Lemgthening: A Study with Gait Analysis , 1989, Journal of pediatric orthopedics.

[71]  Laura H. Smallwood,et al.  Are Current Measurements of Lower Extremity Muscle Architecture Accurate? , 2009, Clinical orthopaedics and related research.

[72]  D. Thelen,et al.  The contribution of passive-elastic mechanisms to lower extremity joint kinetics during human walking. , 2008, Gait & posture.

[73]  Antonie J van den Bogert,et al.  Implicit methods for efficient musculoskeletal simulation and optimal control. , 2011, Procedia IUTAM.

[74]  D. Newham,et al.  Knee muscle isometric strength, voluntary activation and antagonist co-contraction in the first six months after stroke , 2001, Disability and rehabilitation.

[75]  Christopher L. Dembia,et al.  Stretching Your Energetic Budget: How Tendon Compliance Affects the Metabolic Cost of Running , 2016, PloS one.

[76]  F. Dietz,et al.  Medium-term follow-up of Achilles tendon lengthening in the treatment of ankle equinus in cerebral palsy. , 2006, The Iowa orthopaedic journal.

[77]  R. Finkel,et al.  Diagnosis and management of Duchenne muscular dystrophy, part 1: diagnosis, and pharmacological and psychosocial management , 2010, The Lancet Neurology.

[78]  J. Duysens,et al.  Gating and reversal of reflexes in ankle muscles during human walking , 2004, Experimental Brain Research.

[79]  Richard Baker,et al.  Medial gastrocnemius muscle volume and fascicle length in children aged 2 to 5 years with cerebral palsy , 2011, Developmental medicine and child neurology.

[80]  Samuel R. Hamner,et al.  How muscle fiber lengths and velocities affect muscle force generation as humans walk and run at different speeds , 2013, Journal of Experimental Biology.

[81]  Francesco Muntoni,et al.  Muscle MRI findings in patients with limb girdle muscular dystrophy with calpain 3 deficiency (LGMD2A) and early contractures , 2005, Neuromuscular Disorders.

[82]  J C Wall,et al.  A comparison of the range of walking speeds between normal and hemiplegic subjects. , 1995, Scandinavian journal of rehabilitation medicine.

[83]  K. H. Hunt,et al.  Coefficient of Restitution Interpreted as Damping in Vibroimpact , 1975 .

[84]  S Khodadadeh,et al.  Variations of Gait Parameters in Duchenne Muscular Dystrophy , 1990, Proceedings of the Institution of Mechanical Engineers. Part H, Journal of engineering in medicine.

[85]  Craig McDonald,et al.  Diagnosis and management of Duchenne muscular dystrophy, part 2: implementation of multidisciplinary care , 2010, The Lancet Neurology.

[86]  Katherine M Steele,et al.  How much muscle strength is required to walk in a crouch gait? , 2012, Journal of biomechanics.

[87]  Sergey Levine,et al.  Learning to Run challenge: Synthesizing physiologically accurate motion using deep reinforcement learning , 2018, ArXiv.

[88]  Ayman Habib,et al.  OpenSim: Simulating musculoskeletal dynamics and neuromuscular control to study human and animal movement , 2018, PLoS Comput. Biol..

[89]  Stephen F. Keevil,et al.  Lower limb muscle volumes in bilateral spastic cerebral palsy , 2014, Brain and Development.

[90]  Yuval Tassa,et al.  Emergence of Locomotion Behaviours in Rich Environments , 2017, ArXiv.