Distribution of Residence Time of Water around DNA Base Pairs: Governing Factors and the Origin of Heterogeneity

Molecular Origin of Distinct Hydration Dynamics in Double Helical DNA and RNA Sequences

Water molecules are essential to determine the structure of nucleic acids and mediate their interactions with other biomolecules. Here, we characterize the hydration dynamics of analogous DNA and RNA double helices with unprecedented resolution and elucidate the molecular origin of their differences: first, the localization of the slowest hydration water molecules─in the minor groove in DNA, next to phosphates in RNA─and second, the markedly distinct hydration dynamics of the two phosphate oxygen atoms OR and OS in RNA. Using our Extended Jump Model for water reorientation, we assess the relative importance of previously proposed factors, including the local topography, water bridges, and the presence of ions. We show that the slow hydration dynamics at RNA OR sites is not due to bridging water molecules but is caused by both the larger excluded volume and the stronger initial H-bond next to OR, due to the different phosphate orientations in A-form double helical RNA.

Molecular dynamics study of collective water vibrations in a DNA hydration shell

The structure of DNA double helix is stabilized by water molecules and metal counterions that form the ion-hydration shell around the macromolecule. Understanding the role of the ion-hydration shell in the physical mechanisms of the biological functioning of DNA requires detailed studies of its structure and dynamics at the atomistic level. In the present work, the study of collective vibrations of water molecules around the DNA double helix was performed within the framework of classical all-atom molecular dynamics methods. Calculating the vibrational density of states, the vibrations of water molecules in the low-frequency spectra ranged from [Formula: see text]30 to [Formula: see text]300 [Formula: see text] were analyzed for the case of different regions of the DNA double helix (minor groove, major groove, and phosphate groups). The analysis revealed significant differences in the collective vibrations behavior of water molecules in the DNA hydration shell, compared to the vibrations of bulk water. All low-frequency modes of the DNA ion-hydration shell are shifted by about 15-20 [Formula: see text] towards higher frequencies, which is more significant for water molecules in the minor groove region of the double helix. The interactions of water molecules with the atoms of the macromolecule induce intensity decrease of the mode of hydrogen-bond symmetrical stretching near 150 [Formula: see text], leading to the disappearance of this mode in the DNA spectra. The obtained results can provide an interpretation of the experimental data for DNA low-frequency spectra and may be important for the understanding of the processes of indirect protein-nucleic recognition.

Residence Times from Molecular Dynamics Simulations

In this work, efficient methods are proposed for the calculation, from molecular dynamics trajectories, of residence times (RTs) and related quantities. One of these was designed to obtain RT distributions, from which mean residence times (MRTs), residual time distributions, and mean residual times can be computed. This method does not require the assumptions and approximations made by the most commonly used methods. Its link to the most popular method in the literature is identified. It is shown how the much faster new method can be used as a replacement for the latter and the advantages of doing so. Also, a prescription for estimating the uncertainty in the MRTs obtained though the proposed method is provided. Another even faster method for the calculation of the MRTs, their uncertainties, and the mean residual times is also proposed. It yields exactly the same results as the first one but does not allow to obtain the mentioned distributions. Another very popular method, based on autocorrelation functions, for computing MRTs is analyzed in terms of its assumptions and approximations. An alternative, also based on autocorrelation functions, which is faster and requires fewer assumptions is presented. A prescription for the calculation of the uncertainty of the MRTs obtained with the latter method is also provided. In the literature, there are a few methods to discard short transient escapes. Here, an algorithm is suggested that is much faster than the most used one and allows a more integral treatment of the process. Also, it is more widely applicable because it is a preprocessing step that can be used in conjunction with any of the proposed methods mentioned above. The main disadvantage of these two approaches to discard brief escapes is that the maximum duration allowed for an escape to be considered transient appears as a parameter to be determined for the particular system under study. As an alternative, a parameter-free method of a similar character is also proposed to estimate the mean residence time of particles that reached a constant probability of leaving the region of interest.

Long-range DNA-water interactions

DNA functions only in aqueous environments and adopts different conformations depending on the hydration level. The dynamics of hydration water and hydrated DNA leads to rotating and oscillating dipoles that, in turn, give rise to a strong megahertz to terahertz absorption. Investigating the impact of hydration on DNA dynamics and the spectral features of water molecules influenced by DNA, however, is extremely challenging because of the strong absorption of water in the megahertz to terahertz frequency range. In response, we have employed a high-precision megahertz to terahertz dielectric spectrometer, assisted by molecular dynamics simulations, to investigate the dynamics of water molecules within the hydration shells of DNA as well as the collective vibrational motions of hydrated DNA, which are vital to DNA conformation and functionality. Our results reveal that the dynamics of water molecules in a DNA solution is heterogeneous, exhibiting a hierarchy of four distinct relaxation times ranging from ∼8 ps to 1 ns, and the hydration structure of a DNA chain can extend to as far as ∼18 Å from its surface. The low-frequency collective vibrational modes of hydrated DNA have been identified and found to be sensitive to environmental conditions including temperature and hydration level. The results reveal critical information on hydrated DNA dynamics and DNA-water interfaces, which impact the biochemical functions and reactivity of DNA.

Low-frequency vibrations of water molecules in DNA minor groove

Water molecules around the DNA form the hydration shell having different structural and dynamical features in different regions of the double helix. In the DNA minor groove, water molecules are highly ordered and in the case of AT nucleotide sequence, the formation of a hydration spine is observed. In the present research, the vibrations of the hydration spine have been studied to establish the mode of translational vibrations of water molecules in the DNA low-frequency spectra (water-spine vibrations). Using the developed phenomenological model with the parameters determined for different nucleotides of the DNA fragment CGCGAATTCGCG, the frequencies of vibrations of the hydration spine have been obtained within 185 ± 20 cm[Formula: see text] depending on type of nucleotide. The obtained frequencies are in the same region as the translational vibrations of water molecules in the bulk. To select the mode of water-spine vibrations from those modes that are present in the bulk water, the dynamics of DNA with different nucleotide contents has been analyzed, and the possible influence of heavy water has been estimated. The determined features of the mode of water vibrations in the hydration spine of DNA minor groove indicate that this mode may be observed in the experimental spectra.

Origin of Slow Solvation Dynamics in DNA: DAPI in Minor Groove of Dickerson-Drew DNA

The measurement and understanding of collective solvation dynamics in DNA have vital biological implications, as protein and ligand binding to DNA can be directly controlled by complex electrostatic interactions of anionic DNA and surrounding dipolar water, and ions. Time-resolved fluorescence Stokes shift (TRFSS) experiments revealed anomalously slow solvation dynamics in DNA much beyond 100 ps that follow either power-law or slow multiexponential decay over several nanoseconds. The origin of such dispersed dynamics remains difficult to understand. Here we compare results of TRFSS experiments to molecular dynamics (MD) simulations of well-known 4',6-diamidino-2-phenylindole (DAPI)/Dickerson-Drew DNA complex over five decades of time from 100 fs to 10 ns to understand the origin of such dispersed dynamics. We show that the solvation time-correlation function (TCF) calculated from 200 ns simulation trajectory (total 800 ns) captures most features of slow dynamics as measured in TRFSS experiments. Decomposition of TCF into individual components unravels that slow dynamics originating from dynamically coupled DNA-water motion, although contribution from coupled water-Na⁺ motion is non-negligible. The analysis of residence time of water molecules around the probe (DAPI) reveals broad distribution from ∼6 ps to ∼3.5 ns: Several (49 nos.) water molecules show residences time greater than 500 ps, of which at least 14 water molecules show residence times of more than 1 ns in the first solvation shell of DAPI. Most of these slow water molecules are found to occupy two hydration sites in the minor groove near DAPI binding site. The residence time of Na⁺, however, is found to vary within ∼17-120 ps. Remarkably, we find that freezing the DNA fluctuations in simulation eliminates slower dynamics beyond ∼100 ps, where water and Na⁺ dynamics become faster, although strong anticorrelation exists between them. These results indicate that primary origin of slow dynamics lies within the slow fluctuations of DNA parts that couple with nearby slow water and ions to control the dispersed collective solvation dynamics in DNA minor groove.

Surface Topography Effects of Globular Biomolecules on Hydration Water

Hydration water serves as a microscopic manifestation of structural stability and functions of biomolecules. To develop bio-nanomaterials in applications, it is important to study how the surface topography and heterogeneity of biomolecules result in their diversity of the hydration dynamics and energetics. We here performed molecular dynamics simulations combined with the steered molecular dynamics and umbrella sampling to investigate the dynamics and escape process associated with the free energy change of water molecules close to a globular biomolecule, i.e., hemoglobin (Hb) and G-quadruplex DNA (GDNA). The residence time, power of long-time tail, and dipole relaxation time were found to display drastic changes within the averaged hydration shell of 3.0-5.0 Å. Compared with bulk water, in the inner hydration shell, the water dipole moment displays a slower relaxation process and is more oriented toward GDNA than toward Hb, forming a hedgehog-like structure when it surrounds GDNA. In particular, a spine water structure is observed in the GDNA narrow groove. The water isotope effect not only prolongs the dynamic time scales of libration motion in the inner hydration shell and the dipole relaxation processes in the bulk but also strengthens the DNA spine water structure. The potential of the mean force profile reflects the integrity of the hydration shell structure and enables us to obtain detailed insights into the structures formed by water, such as the caged H-bond network and the edge bridge structures; it also reveals that local hydration shell free energy (LHSFE) depends on H-bond rupture processes and ranges from 0.2 to 4.2 kcal/mol. Our results demonstrate that the surface topography of a biomolecule influences the integrity of the hydration shell structure and LHSFE. Our studies are able to identify various further applications in the areas of microfluid devices and nano-dewetting on bioinspired surfaces.

Anticancer drug impact on DNA – a study by neutron spectroscopy coupled with synchrotron-based FTIR and EXAFS

Complementary information on drug–DNA interplay has been achieved for Pt/Pd anticancer agents, by a combined QENS, SR-FTIR-ATR and EXAFS approach.

Microscopic picture of water-ethylene glycol interaction near a model DNA by computer simulation: Concentration dependence, structure, and localized thermodynamics

It is known that crowded molecular environment affects the structure, thermodynamics, and dynamics of macromolecules. Most of the previous works on molecular crowding have majorly focused on the behavior of the macromolecule with less emphasis on the behavior of the crowder and water molecules. In the current study, we have precisely focused on the behavior of the crowder, (ethylene glycol (EG)), salt ions, and water in the presence of a DNA with the increase of the EG concentration. We have probed the behavior of water and crowder using molecular dynamics (MD) simulation and by calculating localized thermodynamic properties. Our results show an interesting competition between EG and water molecules to make hydrogen bonds (H-bond) with DNA. Although the total number of H-bonds involving DNA with both EG and water remains essentially same irrespective of the increase in EG concentration, there is a proportional change in the H-bonding pattern between water-water, EG-EG, and EG-water near DNA and in bulk. At low concentrations of EG, the displacement of water molecules near DNA is relatively easy. However, the displacement of water becomes more difficult as the concentration of EG increases. The density of Na⁺ (Cl^-) near DNA increases (decreases) as the concentration of EG is increased. The density of Cl^- near Na⁺ increases with the increase in EG concentration. It was also found that the average free energy per water in the first solvation shell increases with the increase in EG concentration. Putting all these together, a microscopic picture of EG, water, salt interaction in the presence of DNA, as a function of EG concentration, has emerged.

Effect of water and ionic liquids on biomolecules

The remarkable progress in the field of ionic liquids (ILs) in the last two decades has involved investigations on different aspects of ILs in various conditions. The nontoxic and biocompatible nature of ILs makes them a suitable substance for the storage and application of biomolecules. In this regard, the aqueous IL solutions have attracted a large number of studies to comprehend the role of water in modulating various properties of biomolecules. Here, we review some of the recent studies on aqueous ILs that concern the role of water in altering the behavior of ILs in general and in case of biomolecules solvated in ILs. The different structural and dynamic effects caused by water have been highlighted. We discuss the different modes of IL interaction that are responsible for stabilization and destabilization of proteins and enzymes followed by examples of water effect on this. The role of water in the case of nucleic acid storage in ILs, an area which has mostly been underrated, also has been emphasized. Our discussions highlight the fact that the effects of water on IL behavior are not general and are highly dependent on the nature of the IL under consideration. Overall, we aim to draw attention to the significance of water dynamics in the aqueous IL solutions, a better understanding of which can help in developing superior storage materials for application purposes.

Anomalous water dynamics at surfaces and interfaces: synergistic effects of confinement and surface interactions

In nature, water is often found in contact with surfaces that are extended on the scale of molecule size but small on a macroscopic scale. Examples include lipid bilayers and reverse micelles as well as biomolecules like proteins, DNA and zeolites, to name a few. While the presence of surfaces and interfaces interrupts the continuous hydrogen bond network of liquid water, confinement on a mesoscopic scale introduces new features. Even when extended on a molecular scale, natural and biological surfaces often have features (like charge, hydrophobicity) that vary on the scale of the molecular diameter of water. As a result, many new and exotic features, which are not seen in the bulk, appear in the dynamics of water close to the surface. These different behaviors bear the signature of both water-surface interactions and of confinement. In other words, the altered properties are the result of the synergistic effects of surface-water interactions and confinement. Ultrafast spectroscopy, theoretical modeling and computer simulations together form powerful synergistic approaches towards an understanding of the properties of confined water in such systems as nanocavities, reverse micelles (RMs), water inside and outside biomolecules like proteins and DNA, and also between two hydrophobic walls. We shall review the experimental results and place them in the context of theory and simulations. For water confined within RMs, we discuss the possible interference effects propagating from opposite surfaces. Similar interference is found to give rise to an effective attractive force between two hydrophobic surfaces immersed and kept fixed at a separation of d, with the force showing an exponential dependence on this distance. For protein and DNA hydration, we shall examine a multitude of timescales that arise from frustration effects due to the inherent heterogeneity of these surfaces. We pay particular attention to the role of orientational correlations and modification of the same due to interaction with the surfaces.

Complex Dynamics of Water in Protein Confinement

This paper studies single-molecule and collective dynamics of water confined in protein powders by means of molecular dynamics simulations. The single-particle dynamics show a modest retardation compared to the bulk but become highly stretched in the powder, with the stretching exponent of ≃0.2. The collective dynamics of the total water dipole are affected by intermolecular correlations inside water and by cross-correlations between the water and the protein. The dielectric spectrum of water in the powder has two nearly equal-amplitude peaks: a Debye peak with ≃16 ps relaxation time and a highly stretched peak with the relaxation time of ≃13 ns and a stretching exponent of ≃0.12. The slower relaxation component is not seen in the single-molecule correlation functions and can be assigned to elastic protein motions displacing water in the powder. The loss spectrum of the intermediate scattering function reported by neutron-scattering experiments is also highly stretched, with the high-frequency wing scaling according to a power law. Translational dynamics can become much slower in the powder than in the bulk but are overshadowed by the rotational loss in the overall loss spectrum of neutron scattering.

Water Dynamics in the Hydration Shells of Biomolecules

The structure and function of biomolecules are strongly influenced by their hydration shells. Structural fluctuations and molecular excitations of hydrating water molecules cover a broad range in space and time, from individual water molecules to larger pools and from femtosecond to microsecond time scales. Recent progress in theory and molecular dynamics simulations as well as in ultrafast vibrational spectroscopy has led to new and detailed insight into fluctuations of water structure, elementary water motions, electric fields at hydrated biointerfaces, and processes of vibrational relaxation and energy dissipation. Here, we review recent advances in both theory and experiment, focusing on hydrated DNA, proteins, and phospholipids, and compare dynamics in the hydration shells to bulk water.

[KCl] Dependence of B-DNA Groove Bending Anisotropy

The energetics of B-DNA bending toward the major and minor grooves were quantified by free energy simulations at four different KCl concentrations. Increased [KCl] led to more flexible DNA, with persistence lengths that agreed well with experimental values. At all salt concentrations, major groove bending was preferred, although preferences for major and minor groove bending were similar for the A-tract containing sequence. Since the phosphate repulsions and DNA internal energy favored minor groove bending, the preference for major groove bending was thought to originate from differences in solvation. Water in the minor groove was tighter bound than water in the major groove, and harder to displace than major groove water, which favored the compression of the major groove upon bending. Higher [KCl] decreased the persistence length for both major and minor groove bending but did not greatly affect the free energy spacing between the minor and major groove bending curves. For sequences without A-tracts, salt affected major and minor bending to nearly the same degree, and did not change the preference for major groove bending. For the A-tract containing sequence, an increase in salt concentration decreased the already small energetic difference between major and minor groove bending. Since salts did not significantly affect the relative differences in bending energetics and hydration, it is likely that the increased bending flexibilities upon salt increase are simply due to screening.

Anomalous Dynamics of Water Confined in Protein-Protein and Protein-DNA Interfaces

Confined water often exhibits anomalous properties not observable in the bulk phase. Although water in hydrophobic confinement has been the focus of intense investigation, the behavior of water confined between hydrophilic surfaces, which are more frequently found in biological systems, has not been fully explored. Here, we investigate using molecular dynamics simulations dynamical properties of the water confined in hydrophilic protein-protein and protein-DNA interfaces. We find that the interfacial water exhibits glassy slow relaxations even at 300 K. In particular, the rotational dynamics show a logarithmic decay that was observed in glass-forming liquids at deeply supercooled states. We argue that such slow water dynamics are indeed induced by the hydrophilic binding surfaces, which is in opposition to the picture that the hydration water slaves protein motions. Our results will significantly impact the view on the role of water in biomolecular interactions.

Dispersed dynamics of solvation in G-quadruplex DNA: comparison of dynamic Stokes shifts of probes in parallel and antiparallel quadruplex structures

G-quadruplex DNA (GqDNA) structures play an important role in many specific cellular functions and are promising anti-tumor targets for small molecules (ligands). Here, we measured the dynamic Stokes shift of a ligand (Hoechst) bound to parallel c-Myc (mPu22) GqDNA over five decades of time from 100 fs to 10 ns, and compared it with the previously reported dynamics of DAPI bound to antiparallel human telomeric (hTelo22) GqDNA (Pal et al 2015 J. Phys. Chem. Lett. 6 1754). Stokes shift data from fluorescence up-conversion and time-correlated single photon counting experiments was combined to cover the broad dynamic range. The results show that the solvation dynamics of Hoechst in parallel mPu22 GqDNA follow a power law relaxation, added to fast 2 ps exponential relaxation, from 100 fs to 10 ns, with only a subtle difference of power law exponents in the two ligand-GqDNA systems (0.06 in Hoechst-mPu22 compared to 0.16 in DAPI-hTelo22). We measured steady-state fluorescence spectra and time-resolved anisotropy decays which confirm the tight binding of Hoechst to parallel mPu22 with a binding constant of ~1  ×  10⁵ M^-1. The molecular docking of Hoechst in parallel GqDNA followed by a 50 ns molecular dynamics (MD) simulation on a Hoechst-GqDNA complex reveals that Hoechst binds to one of the outer G-tetrads by end-stacking near G13 and G4, which is different from the binding site of DAPI inside a groove of antiparallel hTelo22 GqDNA. Reconciling previous experimental and simulation results, we assign the 2 ps component to the hydration dynamics of only weakly perturbed water near mPu22 and the power law relaxation to the coupled motion of water and DNA (i.e. DNA backbone, unpaired bases and loops connecting G-tetrads) which come near the Hoechst inside parallel GqDNA.

How Mg2+ ion and water network affect the stability and structure of non-Watson-Crick base pairs in E. coli loop E of 5S rRNA: a molecular dynamics and reference interaction site model (RISM) study

The non-Watson-Crick (non-WC) base pairs of Escherichia coli loop E of 5S rRNA are stabilized by Mg²⁺ ions through water-mediated interaction. It is important to know the synergic role of Mg²⁺ and the water network surrounding Mg²⁺ in stabilizing the non-WC base pairs of RNA. For this purpose, free energy change of the system is calculated using molecular dynamics (MD) simulation as Mg²⁺ is pulled from RNA, which causes disturbance of the water network. It was found that Mg²⁺ remains hexahydrated unless it is close to or far from RNA. In the pentahydrated form, Mg²⁺ interacts directly with RNA. Water network has been identified by two complimentary methods; MD followed by a density-based clustering algorithm and three-dimensional-reference interaction site model. These two methods gave similar results. Identification of water network around Mg²⁺ and non-WC base pairs gives a clue to the strong effect of water network on the stability of this RNA. Based on sequence analysis of all Eubacteria 5s rRNA, we propose that hexahydrated Mg²⁺ is an integral part of this RNA and geometry of base pairs surrounding it adjust to accommodate the [Formula: see text]. Overall the findings from this work can help in understanding the basis of the complex structure and stability of RNA with non-WC base pairs.

Anisotropy of B-DNA Groove Bending

DNA bending is critical for DNA packaging, recognition, and repair, and occurs toward either the major or the minor groove. The anisotropy of B-DNA groove bending was quantified for eight DNA sequences by free energy simulations employing a novel reaction coordinate. The simulations show that bending toward the major groove is preferred for non-A-tracts while the A-tract has a high tendency of bending toward the minor groove. Persistence lengths were generally larger for bending toward the minor groove, which is thought to originate from differences in groove hydration. While this difference in stiffness is one of the factors determining the overall preference of bending direction, the dominant contribution is shown to be a free energy offset between major and minor groove bending. The data suggests that, for the A-tract, this offset is largely determined by inherent structural properties, while differences in groove hydration play a large role for non-A-tracts. By quantifying the energetics of DNA groove bending and rationalizing the origins of the anisotropy, the calculations provide important new insights into a key biological process.

The Clusters-in-a-Liquid Approach for Solvation: New Insights from the Conformer Specific Gas Phase Spectroscopy and Vibrational Optical Activity Spectroscopy

Vibrational optical activity spectroscopies, namely vibrational circular dichroism (VCD) and Raman optical activity (ROA), have been emerged in the past decade as powerful spectroscopic tools for stereochemical information of a wide range of chiral compounds in solution directly. More recently, their applications in unveiling solvent effects, especially those associated with water solvent, have been explored. In this review article, we first select a few examples to demonstrate the unique sensitivity of VCD spectral signatures to both bulk solvent effects and explicit hydrogen-bonding interactions in solution. Second, we discuss the induced solvent chirality, or chiral transfer, VCD spectral features observed in the water bending band region in detail. From these chirality transfer spectral data, the related conformer specific gas phase spectroscopic studies of small chiral hydration clusters, and the associated matrix isolation VCD experiments of hydrogen-bonded complexes in cold rare gas matrices, a general picture of solvation in aqueous solution emerges. In such an aqueous solution, some small chiral hydration clusters, rather than the chiral solutes themselves, are the dominant species and are the ones that contribute mainly to the experimentally observed VCD features. We then review a series of VCD studies of amino acids and their derivatives in aqueous solution under different pHs to emphasize the importance of the inclusion of the bulk solvent effects. These experimental data and the associated theoretical analyses are the foundation for the proposed "clusters-in-a-liquid" approach to account for solvent effects effectively. We present several approaches to identify and build such representative chiral hydration clusters. Recent studies which applied molecular dynamics simulations and the subsequent snapshot averaging approach to generate the ROA, VCD, electronic CD, and optical rotatory dispersion spectra are also reviewed. Challenges associated with the molecular dynamics snapshot approach are discussed and the successes of the seemingly random "ad hoc explicit solvation" reported before are also explained. To further test and improve the "clusters-in-a-liquid" model in practice, future work in terms of conformer specific gas phase spectroscopy of sequential solvation of a chiral solute, matrix isolation VCD measurements of small chiral hydration clusters, and more sophisticated models for the bulk solvent effects would be highly valuable.

Water modulates the ultraslow dynamics of hydrated ionic liquids near CG rich DNA: consequences for DNA stability

The higher affinity of water towards CG-base pairs can alter the dynamics of ionic liquids in solution, causing entropic adversity in binding to the base pairs.