Molecular dynamics simulations of the interaction of glucose with imidazole in aqueous solution

Molecular dynamics simulations were carried out on a concentrated system of β-d-glucopyranose and imidazole molecules in a periodic box of water at 298K. The purpose of the simulations was to determine whether or not there was any tendency for these two solutes to associate in an aqueous environment, as has previously been observed for other planar functional groups from amino acid side chains, such as the indole group of tryptophan or the phenolic group of tyrosine. A weak stacking interaction between β-d-glucopyranose, as a model for cellulose, and imidazole was indeed observed, with an energy of ∼0.25kcal/mol per pair, less than kT. Somewhat surprisingly, considerable imidazole self-association into small aggregates (dimers and trimers) was also observed, with binding energies of ∼0.4kcal/mol per pair, although still less than kT. Considerable non-stacked interactions between glucose and imidazole through hydrogen bonding were also found. These hydrogen bonds primarily involved the N3 atom of imidazole, because the N-H group of N1 was insufficiently polar to compete for water hydrogen bond partners.