The conversion of 4-chlorobenzoyi-CoA to 4-hydroxybenzoyl-CoA catalyzed by 4-chlorobenzoyl-CoA dehalogenase is investigated using combined QM/MM approaches. The calculated potential of mean force at the PM3/CHARMM level supports the proposed nucleophilic aromatic substitution mechanism. In particular, a Meisenheimer intermediate was found, stabilized by hydrogen bonds between the benzoyl carbonyl of the ligand and two backbone amide NHs at positions 64 and 114. Mutation of Glyl 13 to Ala significantly increases the barrier by disrupting the hydrogen bond with the Glyl 14 backbone. The formation of the Meisenheimer complex is accompanied by significant charge redistribution and structural changes in the substrate benzoyl moiety, consistent with experimental observations. Theoretical results suggest that the reaction rate is limited by the formation of the Meisenheimer complex, rather than by its decomposition. A kinetic model based on the calculated free energy profile is found to be consistent with the experimental time course data.