Clustering of glycine molecules in aqueous solution studied by molecular dynamics simulation

The nature of glycine--glycine interactions in aqueous solution has been studied using molecular dynamics simulations at four different concentrations and, in each case, four different temperatures. Although evidence is found for formation of small, transient hydrogen-bonded clusters of glycine molecules, the main type of interaction between glycine molecules is found to be single NH...OC hydrogen bonds. Double-hydrogen-bonded "dimers", which have often been cited as a significant species present in aqueous solutions of glycine, are only observed infrequently. When double-hydrogen-bonded dimers are formed, they dissociate quickly (typically within less than ca. 4 ps), although the broken hydrogen bonds have a higher than average probability of reforming. Several aspects of the clustering of glycine molecules are investigated as a function of both temperature and concentration, including the size distribution of glycine clusters, the radii of gyration of the clusters, and aspects of the lifetimes of glycine-glycine hydrogen bonding by means of hydrogen-bond correlation functions. Diffusion coefficients for the glycine clusters and water molecules are also investigated and provide results in realistic agreement with experimental results.