Neuromuscular disorders often lead to ankle plantar flexor muscle weakness, which impairs ankle push-off power and forward propulsion during gait. To improve walking speed and reduce metabolic cost of transport (mCoT), patients with plantar flexor weakness are provided dorsal-leaf spring ankle-foot orthoses (AFOs). It is widely believed that mCoT during gait depends on the AFO stiffness and an optimal AFO stiffness that minimizes mCoT exists. The biomechanics behind why and how an optimal stiffness exists and benefits individuals with plantar flexor weakness are not well understood. We hypothesized that the AFO would reduce the required support moment and, hence, metabolic cost contributions of the ankle plantar flexor and knee extensor muscles during stance, and reduce hip flexor metabolic cost to initiate swing. To test these hypotheses, we generated neuromusculoskeletal simulations to represent gait of an individual with bilateral plantar flexor weakness wearing an AFO with varying stiffness. Predictions were based on the objective of minimizing mCoT, loading rates at impact and head accelerations at each stiffness level, and the motor patterns were determined via dynamic optimization. The predictive gait simulation results were compared to experimental data from subjects with bilateral plantar flexor weakness walking with varying AFO stiffness. Our simulations demonstrated that reductions in mCoT with increasing stiffness were attributed to reductions in quadriceps metabolic cost during midstance. Increases in mCoT above optimum stiffness were attributed to the increasing metabolic cost of both hip flexor and hamstrings muscles. The insights gained from our predictive gait simulations could inform clinicians on the prescription of personalized AFOs. With further model individualization, simulations based on mCoT minimization may sufficiently predict adaptations to an AFO in individuals with plantar flexor weakness.