Active Site Dynamics in Substrate Hydrolysis Catalyzed by DapE Enzyme and Its Mutants from Hybrid QM/MM-Molecular Dynamics Simulation

Abstract

The mechanism of the catalytic hydrolysis of N-succinyl diaminopimelic acid (SDAP) by the microbial enzyme DapE in its wild-type (wt) form as well as three of its mutants (E134D, H67A, and H349A) is investigated employing a hybrid quantum mechanics/molecular mechanics (QM/MM) method coupled with molecular dynamics (MD) simulations, wherein the time evolution of the atoms of the QM and MM regions are obtained from the forces acting on the individual atoms. The free-energy profiles along the reaction coordinates of this multistep hydrolysis reaction process are explored using a combination of equilibrium and nonequilibrium (umbrella sampling) QM/MM-MD simulation techniques. In the enzyme-substrate complexes of wt-DapE and the E134D mutant, nucleophilic attack is found to be the rate-determining step involving a barrier of 15.3 and 21.5 kcal/mol, respectively, which satisfactorily explains the free energy of activation obtained from kinetic experiments in wt-DapE-SDAP (15.2 kcal/mol) and the 3 orders of magnitude decrease in the catalytic activity due to E134D mutation. The catalysis is found to be quenched in the H67A and H349A mutants of DapE due to conformational rearrangement in the active site induced by the absence of the active site His residues that prohibits activation of the catalytic water molecule.