Non-empirical study of the phosphorylation reaction catalyzed by 4-methyl-5-β-hydroxyethylthiazole kinase:: relevance of the theory of intermolecular interactions

The subject of this study was an analysis of the role of active site residues in the phosphoryl transfer reaction catalyzed by 4-methyl-5-beta-hydroxyethylthiazole kinase (ThiK). The ThiK-catalyzed reaction is of special interest due to the lack of a highly conserved aspartate residue serving as a catalytic base. ONIOM(B3LYP:PM3) models of stationary points along the reaction pathway consisted of reactants, two magnesium ions and several highly conserved ThiK active site residues. The results indicate that an S(N)2-like mechanism of ThiK, with gamma-phosphate acting as an alcohol-activating base is reasonable. Geometries of substrates, transition state and products were utilized in the non-empirical analysis of the physical nature of catalytic interactions taking place in the ThiK active site. The role of particular residues was investigated in terms of their ability to preferentially stabilize the transition state relative to substrates (differential transition state stabilization, DTSS) or products (differential product stabilization, DPS). It seems that Mg2, Glu126 and Cys198 play a major catalytic role, whereas Mg1 and the same Cys198 are responsible for product release. It is remarkable that no dominant role of an electrostatic term in the interactions involved in catalytic activity is observed for product release. Determination of catalytic fields expressing differential electrostatic potential of the transition state with respect to substrates revealed the optimal electrostatic features of an ideal catalyst for the studied reaction. The predicted catalytic environment is in agreement with experimental data showing increased catalytic activity of ThiK upon mutation of Cys198 to aspartate.