A DFT study of solvation effects on the tautomeric equilibrium and catalytic ylide generation of thiamin models.

Thiamin diphosphate (ThDP) is the biologically active form of vitamin B1 and an essential cofactor for a number of enzymes. The effect of solvent polarity on the tautomeric equilibria of ThDP using three model systems of the 4'-aminopyrimidine ring is studied by density functional theory calculations (B3LYP/6-311+G(d,p)//B3LYP/6-31G(d)) in the gas phase and selected solvents (cyclohexane, ether, dichloroethane, and water). Solvation effects are investigated using three different schemes: implicit solvation by a continuum model, explicit solvation by inclusion of three water molecules mimicking the first solvation shell of the enzymatic environment, and by a mixed implicit/explicit solvation model. The 4'-aminopyrimidine tautomer is more stable than the 1',4'-iminopyrimidine tautomer in all solvation schemes employed; however, the trend for the stabilities of the 1',4'-iminopyrimidine tautomer in the solvents depends on the specific ThDP-model. Formation of the catalytic important ylide for ThDP-dependent enzymes by deprotonation of ThDP(C2) is also investigated by localization of transition states for two possible pathways. Only the less stable tautomer, 1',4'-iminopyrimidine ThDP, is able to form the catalytic active ylide. Generation of the ylide through a direct intramolecular proton transfer from ThDP(C2) to the ThDP(N4') nitrogen lone pair is favored by 6 kcal/mol in the gas phase, as compared to a water-mediated ylide generation. However, inclusion of a dielectric medium reduces this difference dramatically. Furthermore, inclusion of two water molecules to model the apoenzymatic environment lowers the activation energies of both direct and water-mediated ylide generation.