Hydrogen bond dynamics of histamine monocation in aqueous solution: Car-Parrinello molecular dynamics and vibrational...

Hydration of histamine was examined by infrared spectroscopy and Car-Parrinello molecular dynamics simulation. Histamine is a neurotransmitter and inflammation mediator, which at physiological pH conditions is present mainly in monocationic form. Our focus was on the part of vibrational spectra that corresponds to histamine N-H stretching, since these degrees of freedom are essential for its interactions with either water molecules or transporters and receptors. Assignment of the experimental spectra revealed a broad feature between 3350 and 2300 cm(-1), being centered at 2950 cm(-1), which includes a mixed contribution from the ring N-H and the aminoethyl N-H stretching vibrations. Computational analysis was performed in two ways: first, by making Fourier transformation on the autocorrelation function of all four N-H bond distances recorded during CPMD run, and second, and most importantly, by incorporating quantum effects through applying an a posteriori quantization of all N-H stretching motions utilizing our snapshot analysis of the fluctuating proton potential. The one-dimensional vibrational Schrödinger equation was solved numerically for each snapshot, and the N-H stretching envelopes were calculated as a superposition of the 0→1 transitions. The agreement with the experiment was much better in the case of the second approach. Our calculations clearly demonstrated that the ring amino group absorbs at higher frequencies than the remaining three amino N-H protons of the protonated aminoethyl group, implying that the chemical bonding in the former group is stronger than in the three amino N-H bonds, thus forming weaker hydrogen bonding with the surrounding solvent molecules. In this way the results of the simulation complemented the experimental spectrum that cannot distinguish between the two sets of protons. The effects of deuteration were also considered. The resulting N-D absorption is narrower and red-shifted. The presented methodology is of general applicability to strongly correlated systems, and it is particularly tuned to provide computational support to vibrational spectroscopy. Perspectives are given for its future applications in computational studies of tunneling in enzyme reactive centers and for receptor activation.