Premium
DFT studies on key intermediates in thiamin catalysis
Author(s) -
Friedemann Rudolf,
Tittmann Kai,
Golbik Ralph,
Hübner Gerhard
Publication year - 2004
Publication title -
international journal of quantum chemistry
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.484
H-Index - 105
eISSN - 1097-461X
pISSN - 0020-7608
DOI - 10.1002/qua.20132
Subject(s) - chemistry , protonation , transketolase , density functional theory , computational chemistry , cofactor , proton , active site , thiamine , transition state , stereochemistry , catalysis , enzyme , organic chemistry , ion , physics , quantum mechanics
The diphosphate ester (ThDP) of thiamin (vitamin B 1 ) is an important cofactor of enzymes within the carbohydrate metabolism. From experiments of site‐specific variants and nuclear magnetic resonance (NMR) studies, it is known that the protonation of the N1′ atom is a significant step in the coenzyme activation by the enzymatic environment. Therefore, we have performed density functional theory (DFT) calculations on the B3LYP/6‐31G* level of N1′H and N1′CH 3 thiamin as model systems to study the protonation and methylation effect on the structure and the electronic properties of the 4′‐amino group. The relaxed rotational barriers related to the C4′‐4′N bond are correlated with findings of 1 H NMR studies and proton/deuterium exchange experiments. Moreover, the effect of N1′ protonation was studied in more detail on the hydroxyethyl‐thiamin carbanion (HETh − ), a key intermediate during catalysis of some ThDP‐dependent enzymes. The relaxed rotational barriers related to the C2C2α bond and the reaction coordinates of the proton transfer 4′NH→C2α of HETh − and N1′H‐HETh − show that they are significantly determined by the protonation at N1′ of HETh − . The influence of the apoenzyme environment on the active coenzyme conformation is modeled in a very simple way. The characteristic torsion angles Φ T and Φ P are considered to be restricted in terms of their values in the corresponding enzyme as well as free optimization parameters. Frequency calculations were performed to characterize the minima and transition state structures, respectively. The applicability of the DFT method was checked by comparing calculations on the MP2‐HF‐SCF/6‐31G* level. © 2004 Wiley Periodicals, Inc. Int J Quantum Chem, 2004