DNA deformability and hydration studied by molecular dynamics simulation

Local dynamics coupled to hydration water determines DNA-sequence-dependent deformability

Molecular dynamics (MD) simulations and quasielastic neutron scattering (QENS) experiments were conducted on two hydrated DNA dodecamers with distinct deformability: 5'CGCG[under AATT]̲CGCG3' and 5'CGCG[under TTAA]̲CGCG3'. The former is known to be rigid and the latter to be flexible. The mean-square displacements of DNA dodecamers exhibit so-called dynamical transition around 200-240 K for both sequences. To investigate the DNA-sequence-dependent dynamics, the dynamics of DNA and hydration water above the transition temperature were examined using both MD simulations and QENS experiments. The fluctuation amplitude of the AATT central tetramer is smaller, and its relaxation time is longer, than that observed in TTAA, suggesting that the AT step is kinetically more stable than TA. The sequence-dependent local base pair step dynamics correlates with the kinetics of breaking the hydrogen bond between DNA and hydration water. The sequence-dependent DNA base pair step fluctuations appear above the dynamical transition temperature. Together with these results, we conclude that DNA deformability is related to the local dynamics of the base pair steps, themselves coupled to hydration water in the minor groove.

Sequence dependencies of DNA deformability and hydration in the minor groove

DNA deformability and hydration are both sequence-dependent and are essential in specific DNA sequence recognition by proteins. However, the relationship between the two is not well understood. Here, systematic molecular dynamics simulations of 136 DNA sequences that differ from each other in their central tetramer revealed that sequence dependence of hydration is clearly correlated with that of deformability. We show that this correlation can be illustrated by four typical cases. Most rigid basepair steps are highly likely to form an ordered hydration pattern composed of one water molecule forming a bridge between the bases of distinct strands, but a few exceptions favor another ordered hydration composed of two water molecules forming such a bridge. Steps with medium deformability can display both of these hydration patterns with frequent transition. Highly flexible steps do not have any stable hydration pattern. A detailed picture of this correlation demonstrates that motions of hydration water molecules and DNA bases are tightly coupled with each other at the atomic level. These results contribute to our understanding of the entropic contribution from water molecules in protein or drug binding and could be applied for the purpose of predicting binding sites.

Comparison of DNA hydration patterns obtained using two distinct computational methods, molecular dynamics simulation and three-dimensional reference interaction site model theory

Because proteins and DNA interact with each other and with various small molecules in the presence of water molecules, we cannot ignore their hydration when discussing their structural and energetic properties. Although high-resolution crystal structure analyses have given us a view of tightly bound water molecules on their surface, the structural data are still insufficient to capture the detailed configurations of water molecules around the surface of these biomolecules. Thanks to the invention of various computational algorithms, computer simulations can now provide an atomic view of hydration. Here, we describe the apparent patterns of DNA hydration calculated by using two different computational methods: Molecular dynamics (MD) simulation and three-dimensional reference interaction site model (3D-RISM) theory. Both methods are promising for obtaining hydration properties, but until now there have been no thorough comparisons of the calculated three-dimensional distributions of hydrating water. This rigorous comparison showed that MD and 3D-RISM provide essentially similar hydration patterns when there is sufficient sampling time for MD and a sufficient number of conformations to describe molecular flexibility for 3D-RISM. This suggests that these two computational methods can be used to complement one another when evaluating the reliability of the calculated hydration patterns.

Nucleic acid solvation: from outside to insight

Nucleic acids are polyanionic molecules that were historically considered to be solely surrounded by a shell of water molecules and a neutralizing cloud of monovalent and divalent cations. In this respect, recent experimental and theoretical reports demonstrate that water molecules within complex nucleic acid structures can display very long residency times, and assist drug binding and catalytic reactions. Finally, anions can also bind to these polyanionic systems. Many of these recent insights are provided by state-of-the-art molecular dynamics simulations of nucleic acid systems, which will be described together with relevant methodological issues.