Microhydration of Protonated Nα-Acetylhistidine: A Theoretical Approach

Extensive exploration of the potential energy surfaces of protonated Nα-acetylhistidine hydrated by 0-3 molecules of water was performed. The methodology combined hierarchical and genealogical (Darwin family tree) approaches using polarizable AMOEBA force field and M06 functional. It is demonstrated that this mixed approach allows recovering a larger number of conformers than the number recovered by using any one of the two methods alone. Hydration enthalpies of protonated Nα-acetylhistidine and of model compounds have been computed using higher theoretical methods, up to the G4MP2 procedure. Excellent agreement with experiment is observed for successive hydration of methylamonium and imidazolium cations using MP2/6-311++G(2d,2p)//M06/6-311++G(d,p) and G4MP2 methods, thereby validating the theory levels used for hydrated protonated Nα-acetylhistidine. It is found that the first hydration enthalpy of protonated Nα-acetylhistidine is ca. 10 kJ mol^-1 lower than that of imidazolium, a result explained by the local environment of the positively charged imidazolium moiety.