Methylglyoxal (MG) is a reactive and toxic compound produced in carbohydrate, lipid, and amino acid metabolism. The glyoxalase system is the main detoxifying route for MG and consists of two enzymes, glyoxalase I (GlxI) and glyoxalase II (GlxII). GlxI catalyzes the formation of S-d-lactoylglutathione from hemithioacetal, and GlxII converts this intermediate to d-lactate. A relationship between the glyoxalase system and some diseases like diabetes has been shown, and inhibiting enzymes of this system may be an effective means of controlling certain diseases. A detailed understanding of the reaction mechanism of an enzyme is essential to the rational design of competitive inhibitors. In this work, we use quantum mechanics/molecular mechanics (QM/MM) calculations and energy refinement utilizing the big-QM and QM/MM thermodynamic cycle perturbation methods to propose a mechanism for the GlxII reaction that starts with a nucleophilic attack of the bridging OH– group on the substrate. The coordination of the substrate to the Zn ions places its electrophilic center close to the hydroxide group, enabling the reaction to proceed. Our estimated reaction energies are in excellent agreement with experimental data, thus demonstrating the reliability of our approach and the proposed mechanism. Additionally, we examined alternative protonation states of Asp-29, Asp-58, Asp-134, and the bridging hydroxide ion in the catalytic process. However, these give less favorable reactions, a poorer reproduction of the crystal structure geometry of the active site, and higher root-mean-squared deviations of the active site residues in molecular dynamics simulations.
