Ammonium ion representation in Monte Carlo simulations of biomolecular solutions

Monte Carlo computer simulation techniques may be used to predict structural properties of solvent networks in helical fragments of nucleic acids, provided that suitable potential functions are available to describe the interactions between nucleic acid atoms, water and counterions. Previous studies have shown that simple non-bonded and point charge parameters are adequate for mononuclear ions such as sodium and calcium. In this study a model interaction potential for the polynuclear ammonium ion is evaluated. The parameters used take account of the distribution of charge over the constituent atoms in the ion. Simulations are carried out on the ammonium salt of a small nucleic acid crystal hydrate and a comparison is made between the predicted and experimental results. It is shown that the simulated structure is in reasonable agreement with experiment. It is therefore feasible to use this potential in studies of ammonium-containing bimolecular systems.