PNNL

Abstract

The research described in this product was performed in part in the Environmental Molecular Sciences Laboratory, a national scientific user facility sponsored by the Department of Energy's Office of Biological and Environmental Research and located at Pacific Northwest National Laboratory. Dynamic protonation equilibria in water of one 4-methylimidazole molecule as well as for pairs and groups consisting of 4- methylimidazole, acetic acid and bridging water molecules are studied using Q-HOP molecular dynamics simulation. We find a qualitatively different protonation behavior of 4-methylimidazole compared to that of acetic acid. On one hand, deprotonated, neutral 4-methylimidazole cannot as easily attract a freely diffusing extra proton from solution. Once the proton is bound, however, it remains tightly bound on a time scale of tens of nanoseconds. In a linear chain composed of acetic acid, a separating water molecule and 4-methylimidazole, an excess proton is equally shared between 4-methylimidazole and water. When a water molecule is linearly placed between two acetic acid molecules, the excess proton is always found on the central water. On the other hand, an excess proton in a 4-methylimidazole-water- 4-methylimidazole chain is always localized on one of the two 4- methylimidazoles. These findings are of interest to the discussion of proton transfer along chains of amino acids and water molecules in biomolecules.