Myosin regulates muscle function through a complex cycle of conformational rearrangements coupled with the hydrolysis of adenosine triphosphate (ATP). The recovery stroke reorganizes the myosin active site to hydrolyze ATP and cross bridge with the thin filament to produce muscle contraction. Engineered mutations K84M and R704E in Dictyostelium myosin have been designed to specifically inhibit the recovery stroke and have been shown to indirectly affect the ATPase activity of myosin. We investigated these mutagenic perturbations to the recovery stroke and generated thermodynamically correct and unbiased trajectories for native ATP hydrolysis with computationally enhanced sampling methods. Our methodology was able to resolve experimentally observed changes to kinetic and equilibrium dynamics for the recovery stroke with the correct prediction in the severity of these changes. For ATP hydrolysis, the sequential nature along with the stabilization of a metaphosphate intermediate was observed in agreement with previous studies. However, we observed glutamate 459 being utilized as a proton abstractor to prime the attacking water instead of a lytic water, a phenomenon not well categorized in myosin but has in other ATPases. Both rare event methodologies can be extended to human myosin to investigate isoformic differences from Dictyostelium and scan cardiomyopathic mutations to see differential perturbations to kinetics of other conformational changes in myosin such as the power stroke.