Theoretical description of biomolecular hydration. Application to A-DNA.

The local density of water molecules around a biomolecule is constructed from calculated two- and three-points correlation functions of polar solvents in water using a Potential-of-Mean-Force (PMF) expansion. As a simple approximation, the hydration of all polar (including charged) groups in a biomolecule is represented by the hydration of water oxygen in bulk water, and the effect of non-polar groups on hydration are neglected, except for excluded volume effects. Pair and triplet correlation functions are calculated by molecular dynamics simulations. We present calculations of the structural hydration for ideal A-DNA molecules with sequences [d(CG)5]2 and [d(C5G5)]2. We find that this method can accurately reproduce the hydration patterns of A-DNA observed in neutron diffraction experiments on oriented DNA fibers (P. Langan et al. J. Biomol. Struct. Dyn., 10, 489 (1992)).