Abstract
In this DFT study, the substrate promiscuity of the binuclear [Fe(II)-Zn(II)] core containing glycerophosphodiesterase (GpdQ) from Enterobacter aerogenes has been investigated through the hydrolysis of three chemically diverse groups of substrates: i.e., phosphomono-, phosphodi-, and phosphotriesters. The hydrolysis of these substrates is studied by comparing stepwise, concerted, and substrate-assisted mechanisms. Both the stepwise and concerted mechanisms occur with similar barriers, while the energetics for the substrate-assisted mechanism are significantly less favorable. Irrespective of the mechanism, active site residue His217 plays a critical role, in agreement with structural, kinetics, and spectroscopic data, but the transition state of the reaction depends on the identity of the substrate (dissociative for the triester paraoxon, associative for the monoester 4-nitrophenyl phosphate (NPP), and in-between for the diesters glycerol-3-phosphoethanolamine (GPE) and bis(4-nitrophenyl)phosphate (BNPP)). In good agreement with available kinetic and spectrophotometric data, the calculations highlight the preference of GpdQ for diester substrates, followed by tri- and monoesters. For substrates with two different types of scissile bonds (paraoxon and GPE) a clear preference for the bond with the stronger electron withdrawing leaving group was observed. The extensive agreement between experimental data and DFT calculations enhances the understanding of the mechanism of GpdQ-catalyzed hydrolysis and paves the way for the rational design of optimized catalysts for the hydrolysis of different types of phosphoesters.
