Molecular recognition of aldehydes by aldehyde dehydrogenase and mechanism of nucleophile activation

Experimental structural data on the state of substrates bound to class 3 Aldehyde Dehydrogenases (ALDH3A1) is currently unknown. We have utilized molecular mechanics (MM) simulations, in conjunction with new force field parameters for aldehydes, to study the atomic details of benzaldehyde binding to ALDH3A1. Our results indicate that while the nucleophilic Cys243 must be in the neutral state to form what are commonly called near‐attack conformers (NACs), these structures do not correlate with increased complexation energy calculated with the MM‐Generalized Born Molecular Volume (GBMV) method. The negatively charged Cys243 (thiolate form) of ALDH3A1 also binds benzaldehyde in a stable conformation but in this complex the sulfur of Cys243 is oriented away from benzaldehyde yet yields the most favorable MM‐GBMV complexation energy. The identity of the general base, Glu209 or Glu333, in ALDHs remains uncertain. The MM simulations reveal structural and possible functional roles for both Glu209 and Glu333. Structures from the MM simulations that would support either glutamate residue as the general base were further examined with Hybrid Quantum Mechanical (QM)/MM simulations. These simulations show that, with the PM3/OPLS potential, Glu209 must go through a step‐wise mechanism to activate Cys243 through an intervening water molecule while Glu333 can go through a more favorable concerted mechanism for the same activation process. Proteins 2004. © 2004 Wiley‐Liss, Inc.