Background: The stiffness of a dorsal leaf AFO that minimizes walking energy cost in people with plantarflexor weakness varies between individuals. Using predictive simulations, we studied the effects of plantarflexor weakness, passive plantarflexor stiffness, body mass, and walking speed on the optimal AFO stiffness for energy cost reduction. Methods: We employed a planar, nine degrees-of-freedom musculoskeletal model, in which for validation maximal strength of the plantar flexors was reduced by 80%. Walking simulations, driven by minimizing a comprehensive cost function of which energy cost was the main contributor, were generated using a reflex-based controller. Simulations of walking without and with an AFO with stiffnesses between 0.9 and 8.7 Nm/degree were generated. After validation against experimental data of 11 people with plantarflexor weakness using the Root-mean-square error (RMSE), we systematically changed plantarflexor weakness (range 40–90% weakness), passive plantarflexor stiffness (range: 20–200% of normal), body mass (+ 30%) and walking speed (range: 0.8–1.2 m/s) in our baseline model to evaluate their effect on the optimal AFO stiffness for energy cost minimization. Results: Our simulations had a RMSE < 2 for all lower limb joint kinetics and kinematics except the knee and hip power for walking without AFO. When systematically varying model parameters, more severe plantarflexor weakness, lower passive plantarflexor stiffness, higher body mass and walking speed increased the optimal AFO stiffness for energy cost minimization, with the largest effects for severity of plantarflexor weakness. Conclusions: Our forward simulations demonstrate that in individuals with bilateral plantarflexor the necessary AFO stiffness for walking energy cost minimization is largely affected by severity of plantarflexor weakness, while variation in walking speed, passive muscle stiffness and body mass influence the optimal stiffness to a lesser extent. That gait deviations without AFO are overestimated may have exaggerated the required support of the AFO to minimize walking energy cost. Future research should focus on improving predictive simulations in order to implement personalized predictions in usual care. Trial Registration Nederlands Trial Register 5170. Registration date: May 7th 2015. http://www.trialregister.nl/trialreg/admin/rctview.asp?TC=5170.

@en

To maximize effects of dorsal leaf ankle foot orthoses (AFOs) on gait in people with bilateral plantarflexor weakness, the AFO properties should be matched to the individual. However, how AFO properties interact regarding their effect on gait function is unknown. We studied the interaction of AFO bending stiffness with neutral angle and footplate stiffness on the effect of bending stiffness on walking energy cost, gait kinematics and kinetics in people with plantarflexor weakness by employing predictive simulations. Our simulation framework consisted of a planar 11 degrees of freedom model, containing 11 muscles activated by a reflex-based neuromuscular controller. The controller was optimized by a comprehensive cost function, predominantly minimizing walking energy cost. The AFO bending and footplate stiffness were modelled as torsional springs around the ankle and metatarsal joint. The neutral angle of the AFO was defined as the angle in the sagittal plane at which the moment of the ankle torsional spring was zero. Simulations without AFO and with AFO for 9 bending stiffnesses (0–14 Nm/degree), 3 neutral angles (0–3-6 degrees dorsiflexion) and 3 footplate stiffnesses (0–0.5–2.0 Nm/degree) were performed. When changing neutral angle towards dorsiflexion, a higher AFO bending stiffness minimized energy cost of walking and normalized joint kinematics and kinetics. Footplate stiffness mainly affected MTP joint kinematics and kinetics, while no systematic and only marginal effects on energy cost were found. In conclusion, the interaction of the AFO bending stiffness and neutral angle in bilateral plantarflexor weakness, suggests that these should both be considered together when matching AFO properties to the individual patient.

@en

Computational modelling of the neuromusculoskeletal system (NMSS) can potentially provide detailed insight into muscle function to optimize treatment planning and evaluation in cerebral palsy (CP). Commonly, static optimization is used to solve the redundancy problem in estimating muscle forces, assuming, for example, minimization of the muscle activation squared. However, since the primary problem in children with CP is an aberrant motor control [1], it is questionable if using this criterion is applicable in this clinical population. Electromyography (EMG) might be used to further inform optimization procedures and improve model performance.@en