Journal of Physical Chemistry B, Vol.117, No.46, 14261-14272, 2013
Conformational Compression and Barrier Height Heterogeneity in the N-Acetylglutamate Kinase
The role of motions in the catalytic cycle of an enzyme is the subject of much debate. Crystallographic results for the enzyme N-acetyl-L-glutamate kinase (NAGK), which is a suitable target for antibacterial drugs, suggest that a conformational compression of the active site favors catalysis. We have used a QM/MM scheme to compute energy profiles of the phosphoryl transfer reaction for 20 conformations of NAGK, starting from four crystal structures that represent different stages of the catalytic process. All paths show a common associative mechanism but with a wide range of barrier heights. The position of several active site residues and water molecules are found to determine the energetic barrier of each conformation, as revealed by principal component and partial least-squares chemometric analyses. In particular, conformations in which the two substrates have a shorter distance separation and a more linear mutual orientation tend to have lower energetic barriers, thus supporting the putative role of conformational compressive motions in catalysis. Interestingly, these motions are the same that lead to opening of the active site, which molecular dynamics simulations indicate is a fast process when the enzyme is free of substrates. Despite the lack of extended sampling, the energy barrier we calculate for the chemical step lies significantly below the apparent energetic barrier derived from experiment. Although not conclusive, this result supports a previous hypothesis, also derived from experiment, that conformational motions, rather than the chemical step, are rate limiting.