Quantum Mechanics/Molecular Mechanics Study of the Reaction Mechanism of Glyoxalase I.

Glyoxalase I (GlxI) is a member of the glyoxalase system, which is important in cell detoxification and converts hemithioacetals of methylglyoxal (a cytotoxic byproduct of sugar metabolism that may react with DNA or proteins and introduce nucleic acid strand breaks, elevated mutation frequencies, and structural or functional changes of the proteins) and glutathione into d-lactate. GlxI accepts both the S and R enantiomers of hemithioacetal, but converts them to only the S-d enantiomer of lactoylglutathione. Interestingly, the enzyme shows this unusual specificity with a rather symmetric active site (a Zn ion coordinated to two glutamate residues; Glu-99 and Glu-172), making the investigation of its reaction mechanism challenging. Herein, we have performed a series of combined quantum mechanics and molecular mechanics calculations to study the reaction mechanism of GlxI. The substrate can bind to the enzyme in two different modes, depending on the direction of its alcoholic proton (H2; toward Glu-99 or Glu-172). Our results show that the S substrate can react only if H2 is directed toward Glu-99 and the R substrate only if H2 is directed toward Glu-172. In both cases, the reactions lead to the experimentally observed S-d enantiomer of the product. In addition, the results do not show any low-energy paths to the wrong enantiomer of the product from neither the S nor the R substrate. Previous studies have presented several opposing mechanisms for the conversion of R and S enantiomers of the substrate to the correct enantiomer of the product. Our results confirm one of them for the S substrate, but propose a new one for the R substrate.

A novel mechanism of heme degradation to biliverdin studied by QM/MM and QM calculations.

Heme degradation by heme oxygenase enzymes is important for maintaining iron homeostasis and prevention of oxidative stress. Previous studies have reported that heme degradation proceeds through three consecutive steps of O2 activation: the regiospecific self-hydroxylation of heme, the conversion of hydroxyheme to verdoheme and CO, and the cleavage of the verdoheme macrocycle to release biliverdin and free ferrous iron. Our results indicate that in the second step of heme degradation, not only verdoheme is generated but ring opening and biliverdin production also occur. We have performed QM-cluster and QM/MM calculations, which show that calculations with H2O as the axial ligand of Fe give the lowest barrier. In the QM-cluster calculation, the reaction is exothermic by -85 kcal mol-1 and the rate-limiting barrier is 5 kcal mol-1, whereas the corresponding QM/MM calculations give a slightly lower barrier of 3 kcal mol-1, owing to strong hydrogen bonds and the protein environment.

QM/MM Study of the Conversion of Oxophlorin into Verdoheme by Heme Oxygenase.

Heme oxygenase is an enzyme that degrades heme, thereby recycling iron in most organisms, including humans. Pervious density functional theory (DFT) calculations have suggested that iron(III) hydroxyheme, an intermediate generated in the first step of heme degradation by heme oxygenase, is converted to iron(III) superoxo oxophlorin in the presence of dioxygen. In this article, we have studied the detailed mechanism of conversion of iron(III) superoxo oxophlorin to verdoheme by using combined quantum mechanics and molecular mechanics (QM/MM) calculations. The calculations employed the B3LYP method and the def2-QZVP basis set, considering dispersion effects with the DFT-D3 approach, obtaining accurate energies with large QM regions of almost 1000 atoms. The reaction was found to be exothermic by -35 kcal/mol, with a rate-determining barrier of 19 kcal/mol in the doublet state. The protein environment and especially water in the enzyme pocket significantly affects the reaction by decreasing the reaction activation energies and changing the structures by providing strategic hydrogen bonds.