Kinetics of DNA hydration

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.

Dynamical Crossover of Interfacial Water upon Melting of a Lipid Bilayer: Influence of Different Parts of the Headgroups

All-atom molecular dynamics simulations of a 1,2-dimyristoyl-sn-glycero-3-phosphocholine bilayer in contact with liquid water were performed at different temperatures ranging from 285 to 320 K. We have investigated the heterogeneity and dynamical transitions in interfacial water as the lipid bilayer undergoes a melting transition. Results are obtained for water at the outer surface of the bilayer and for those buried more deeply in the lipid chains of the bilayer. It is found that lipid bilayer melting influences both the structure and dynamics of interfacial water. The number of interfacial water molecules shows a jump in the melting of the bilayer. The temperature dependence of the diffusivity and orientational relaxation of interfacial water molecules exhibits a dynamical crossover upon melting of the bilayer. The extent of dynamical crossover is found to be rather strong with significant changes in activation barriers for interfacial water around the carbonyl groups, which are deeply buried toward the lipid chains of the bilayer. The dynamical crossover gradually decreases as one moves further away from the outer surface, and it essentially vanishes for water in the region of 5-10 Å from the outer surface. It is found that the lipid melting-induced dynamical crossover of interfacial water is significant only for water that is in close proximity to the bilayer surface or deeply buried into it. The current results reveal that water molecules in different parts of the interface respond differently on melting of the bilayer. The current study also shows that the carbonyl-bound water molecules can play an important role in the phase transition of the bilayer as the temperature is raised through its melting point.

The relayed nuclear Overhauser effect in magnetization transfer and chemical exchange saturation transfer MRI

Magnetic resonance (MR) is a powerful technique for noninvasively probing molecular species in vivo but suffers from low signal sensitivity. Saturation transfer (ST) MRI approaches, including chemical exchange saturation transfer (CEST) and conventional magnetization transfer contrast (MTC), allow imaging of low-concentration molecular components with enhanced sensitivity using indirect detection via the abundant water proton pool. Several recent studies have shown the utility of chemical exchange relayed nuclear Overhauser effect (rNOE) contrast originating from nonexchangeable carbon-bound protons in mobile macromolecules in solution. In this review, we describe the mechanisms leading to the occurrence of rNOE-based signals in the water saturation spectrum (Z-spectrum), including those from mobile and immobile molecular sources and from molecular binding. While it is becoming clear that MTC is mainly an rNOE-based signal, we continue to use the classical MTC nomenclature to separate it from the rNOE signals of mobile macromolecules, which we will refer to as rNOEs. Some emerging applications of the use of rNOEs for probing macromolecular solution components such as proteins and carbohydrates in vivo or studying the binding of small substrates are discussed.

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.

Modulation of the Directionality of Hole Transfer between the Base and the Sugar-Phosphate Backbone in DNA with the Number of Sulfur Atoms in the Phosphate Group

This work shows that S atom substitution in phosphate controls the directionality of hole transfer processes between the base and sugar-phosphate backbone in DNA systems. The investigation combines synthesis, electron spin resonance (ESR) studies in supercooled homogeneous solution, pulse radiolysis in aqueous solution at ambient temperature, and density functional theory (DFT) calculations of in-house synthesized model compound dimethylphosphorothioate (DMTP(O^-)═S) and nucleotide (5'-O-methoxyphosphorothioyl-2'-deoxyguanosine (G-P(O^-)═S)). ESR investigations show that DMTP(O^-)═S reacts with Cl₂^•- to form the σ²σ*¹ adduct radical -P-S[Formula: see text]Cl, which subsequently reacts with DMTP(O^-)═S to produce [-P-S[Formula: see text]S-P-]^-. -P-S[Formula: see text]Cl in G-P(O^-)═S undergoes hole transfer to Gua, forming the cation radical (G^•+) via thermally activated hopping. However, pulse radiolysis measurements show that DMTP(O^-)═S forms the thiyl radical (-P-S^•) by one-electron oxidation, which did not produce [-P-S[Formula: see text]S-P-]^-. Gua in G-P(O^-)═S is oxidized unimolecularly by the -P-S^• intermediate in the sub-picosecond range. DFT thermochemical calculations explain the differences in ESR and pulse radiolysis results obtained at different temperatures.

Revealing DNA Structure at Liquid/Solid Interfaces by AFM-Based High-Resolution Imaging and Molecular Spectroscopy

DNA covers the genetic information in all living organisms. Numerous intrinsic and extrinsic factors may influence the local structure of the DNA molecule or compromise its integrity. Detailed understanding of structural modifications of DNA resulting from interactions with other molecules and surrounding environment is of central importance for the future development of medicine and pharmacology. In this paper, we review the recent achievements in research on DNA structure at nanoscale. In particular, we focused on the molecular structure of DNA revealed by high-resolution AFM (Atomic Force Microscopy) imaging at liquid/solid interfaces. Such detailed structural studies were driven by the technical developments made in SPM (Scanning Probe Microscopy) techniques. Therefore, we describe here the working principles of AFM modes allowing high-resolution visualization of DNA structure under native (liquid) environment. While AFM provides well-resolved structure of molecules at nanoscale, it does not reveal the chemical structure and composition of studied samples. The simultaneous information combining the structural and chemical details of studied analyte allows achieve a comprehensive picture of investigated phenomenon. Therefore, we also summarize recent molecular spectroscopy studies, including Tip-Enhanced Raman Spectroscopy (TERS), on the DNA structure and its structural rearrangements.

Deprotonation Dynamics of Guanine Radical Cations†

This review is dedicated to guanine radical cations (G⁺ )^· that are precursors to oxidatively generated damage to DNA. (G⁺ )^· are unstable in neutral aqueous solution and tend to lose a proton. The deprotonation process has been studied by time-resolved absorption experiments in which (G⁺ )^· radicals are produced either by an electron abstraction reaction, using an external oxidant, or by low-energy/low-intensity photoionization of DNA. Both the position of the released proton and the dynamics of the process depend on the secondary DNA structure. While deprotonation in duplex DNA leads to (G-H1)^· radicals, in guanine quadruplexes the (G-H2)^· analogs are observed. Deprotonation in monomeric guanosine proceeds with a time constant of ˜60 ns; in genomic DNA, it is completed within 2 µs; and in guanine quadruplexes, it spans from at least 30 ns to over 50 µs. Such a deprotonation dynamics in four-stranded structures, extended over more than three decades of times, is correlated with the anisotropic structure of DNA and the mobility of its hydration shell. In this case, commonly used second-order reaction models are inappropriate for its description.

The Time Scale of Electronic Resonance in Oxidized DNA as Modulated by Solvent Response: An MD/QM-MM Study

The time needed to establish electronic resonant conditions for charge transfer in oxidized DNA has been evaluated by molecular dynamics simulations followed by QM/MM computations which include counterions and a realistic solvation shell. The solvent response is predicted to take ca. 800-1000 ps to bring two guanine sites into resonance, a range of values in reasonable agreement with the estimate previously obtained by a kinetic model able to correctly reproduce the observed yield ratios of oxidative damage for several sequences of oxidized 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.

Same fold, different properties: polarizable molecular dynamics simulations of telomeric and TERRA G-quadruplexes

DNA and RNA sequences rich in guanine can fold into noncanonical structures called G-quadruplexes (GQs), which exhibit a common stem structure of Hoogsteen hydrogen-bonded guanine tetrads and diverse loop structures. GQ sequence motifs are overrepresented in promoters, origins of replication, telomeres, and untranslated regions in mRNA, suggesting roles in modulating gene expression and preserving genomic integrity. Given these roles and unique aspects of different structures, GQs are attractive targets for drug design, but greater insight into GQ folding pathways and the interactions stabilizing them is required. Here, we performed molecular dynamics simulations to study two bimolecular GQs, a telomeric DNA GQ and the analogous telomeric repeat-containing RNA (TERRA) GQ. We applied the Drude polarizable force field, which we show outperforms the additive CHARMM36 force field in both ion retention and maintenance of the GQ folds. The polarizable simulations reveal that the GQs bind bulk K+ ions differently, and that the TERRA GQ accumulates more K+ ions, suggesting different ion interactions stabilize these structures. Nucleobase dipole moments vary as a function of position and also contribute to ion binding. Finally, we show that the TERRA GQ is more sensitive than the telomeric DNA GQ to water-mediated modulation of ion-induced dipole-dipole interactions.

A New Look into the Mode of Action of Metal-Based Anticancer Drugs

The mode of action of Pt- and Pd-based anticancer agents (cisplatin and Pd₂Spm) was studied by characterising their impact on DNA. Changes in conformation and mobility at the molecular level in hydrated DNA were analysed by quasi-elastic and inelastic neutron scattering techniques (QENS and INS), coupled to Fourier transform infrared (FTIR) and microRaman spectroscopies. Although INS, FTIR and Raman revealed drug-triggered changes in the phosphate groups and the double helix base pairing, QENS allowed access to the nanosecond motions of the biomolecule’s backbone and confined hydration water within the minor groove. Distinct effects were observed for cisplatin and Pd₂Spm, the former having a predominant effect on DNA’s spine of hydration, whereas the latter had a higher influence on the backbone dynamics. This is an innovative way of tackling a drug’s mode of action, mediated by the hydration waters within its pharmacological target (DNA).

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.

Overhauser Dynamic Nuclear Polarization for the Study of Hydration Dynamics, Explained

We outline the physical properties of hydration water that are captured by Overhauser Dynamic Nuclear Polarization (ODNP) relaxometry and explore the insights that ODNP yields about the water and the surface that this water is coupled to. As ODNP relies on the pairwise cross-relaxation between the electron spin of a spin probe and a proton nuclear spin of water, it captures the dynamics of single-particle diffusion of an ensemble of water molecules moving near the spin probe. ODNP principally utilizes the same physics as other nuclear magnetic resonance (NMR) relaxometry (i.e., relaxation measurement) techniques. However, in ODNP, electron paramagnetic resonance (EPR) excites the electron spins probes and their high net polarization acts as a signal amplifier. Furthermore, it renders ODNP parameters highly sensitive to water moving at rates commensurate with the EPR frequency of the spin probe (typically 10GHz). Also, ODNP selectively enhances the NMR signal contributions of water moving within close proximity to the spin label. As a result, ODNP can capture ps-ns movements of hydration waters with high sensitivity and locality, even in samples with protein concentrations as dilute as 10 µM. To date, the utility of the ODNP technique has been demonstrated for two major applications: the characterization of the spatial variation in the properties of the hydration layer of proteins or other surfaces displaying topological diversity, and the identification of structural properties emerging from highly disordered proteins and protein domains. The former has been shown to correlate well with the properties of hydration water predicted by MD simulations and has been shown capable of evaluating the hydrophilicity or hydrophobicity of a surface. The latter has been demonstrated for studies of an interhelical loop of proteorhodopsin, the partial structure of α-synuclein embedded at the lipid membrane surface, incipient structures adopted by tau proteins en route to fibrils, and the structure and hydration profile of a transmembrane peptide. This chapter focuses on offering a mechanistic understanding of the ODNP measurement and the molecular dynamics encoded in the ODNP parameters. In particular, it clarifies how the electron-nuclear dipolar coupling encodes information about the molecular dynamics in the nuclear spin self-relaxation and, more importantly, the electron-nuclear spin cross-relaxation rates. The clarification of the molecular dynamics underlying ODNP should assist in establishing a connection to theory and computer simulation that will offer far richer interpretations of ODNP results in future studies.

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.

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.

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.

DNA's Chiral Spine of Hydration

The iconic helical structure of DNA is stabilized by the solvation environment, where a change in the hydration state can lead to dramatic changes to the DNA structure. X-ray diffraction experiments at cryogenic temperatures have shown crystallographic water molecules in the minor groove of DNA, which has led to the notion of a spine of hydration of DNA. Here, chiral nonlinear vibrational spectroscopy of two DNA sequences shows that not only do such structural water molecules exist in solution at ambient conditions but that they form a chiral superstructure: a chiral spine of hydration. This is the first observation of a chiral water superstructure templated by a biomolecule. While the biological relevance of a chiral spine of hydration is unknown, the method provides a direct way to interrogate the properties of the hydration environment of DNA and water in biological settings without the use of labels.

[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.

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.

Effect of temperature on the structure and hydration layer of TATA-box DNA: A molecular dynamics simulation study

DNA within the living cells experiences a diverse range of temperature, ranging from freezing condition to hot spring water. How the structure, the mechanical properties of DNA, and the solvation dynamics around DNA changes with the temperature is important to understand the functionality of DNA under those acute temperature conditions. In that notion, we have carried out molecular dynamics simulations of a DNA oligomer, containing TATA-box sequence for three different temperatures (250K, 300K and 350K). We observed that the structure of the DNA, in terms of backbone torsion angles, sugar pucker, base pair parameters, and base pair step parameters, did not show any unusual properties within the studied range of temperatures, but significant structural alteration was noticed between BI and BII forms at higher temperature. As expected, the flexibility of the DNA, in terms of the torsional rigidity and the bending rigidity is highly temperature dependent, confirming that flexibility increases with increase in temperature. Additionally, the groove widths of the studied DNA showed temperature sensitivity, specifically, the major groove width decreases and the minor groove width increases, respectively, with the increase in temperature. We observed that at higher temperature, water around both the major and the minor groove of the DNA is less structured. However, the water dynamics around the minor groove of the DNA is more restricted as compared to the water around the major groove throughout the studied range of temperatures, without any anomalous behavior.

Are Waters around RNA More than Just a Solvent? - An Insight from Molecular Dynamics Simulations

Hydrating water molecules are believed to be an inherent part of the RNA structure and have a considerable impact on RNA conformation. However, the magnitude and mechanism of the interplay between water molecules and the RNA structure are still poorly understood. In principle, such hydration effects can be studied by molecular dynamics (MD) simulations. In our recent MD studies, we observed that the choice of water model has a visible impact on the predicted structure and structural dynamics of RNA and, in particular, has a larger effect than type, parametrization, and concentration of the ions. Furthermore, the water model effect is sequence dependent and modulates the sequence dependence of A-RNA helical parameters. Clearly, the sensitivity of A-RNA structural dynamics to the water model parametrization is a rather spurious effect that complicates MD studies of RNA molecules. These results nevertheless suggest that the sequence dependence of the A-RNA structure, usually attributed to base stacking, might be driven by the structural dynamics of specific hydration. Here, we present a systematic MD study that aimed to (i) clarify the atomistic mechanism of the water model sensitivity and (ii) discover whether and to what extent specific hydration modulates the A-RNA structural variability. We carried out an extended set of MD simulations of canonical A-RNA duplexes with TIP3P, TIP4P/2005, TIP5P, and SPC/E explicit water models and found that different water models provided a different extent of water bridging between 2'-OH groups across the minor groove, which in turn influences their distance and consequently also inclination, roll, and slide parameters. Minor groove hydration is also responsible for the sequence dependence of these helical parameters. Our simulations suggest that TIP5P is not optimal for RNA simulations.

Thermodynamic properties of water molecules in the presence of cosolute depend on DNA structure: a study using grid inhomogeneous solvation theory

In conditions that mimic those of the living cell, where various biomolecules and other components are present, DNA strands can adopt many structures in addition to the canonical B-form duplex. Previous studies in the presence of cosolutes that induce molecular crowding showed that thermal stabilities of DNA structures are associated with the properties of the water molecules around the DNAs. To understand how cosolutes, such as ethylene glycol, affect the thermal stability of DNA structures, we investigated the thermodynamic properties of water molecules around a hairpin duplex and a G-quadruplex using grid inhomogeneous solvation theory (GIST) with or without cosolutes. Our analysis indicated that (i) cosolutes increased the free energy of water molecules around DNA by disrupting water–water interactions, (ii) ethylene glycol more effectively disrupted water–water interactions around Watson–Crick base pairs than those around G-quartets or non-paired bases, (iii) due to the negative electrostatic potential there was a thicker hydration shell around G-quartets than around Watson–Crick-paired bases. Our findings suggest that the thermal stability of the hydration shell around DNAs is one factor that affects the thermal stabilities of DNA structures under the crowding conditions.

Anomalously Rapid Hydration Water Diffusion Dynamics Near DNA Surfaces

The emerging Overhauser effect dynamic nuclear polarization (ODNP) technique measures the translational mobility of water within the vicinity (5-15 Å) of preselected sites. The work presented here expands the capabilities of the ODNP technique and illuminates an important, previously unseen, property of the translational diffusion dynamics of water at the surface of DNA duplexes. We attach nitroxide radicals (i.e., spin labels) to multiple phosphate backbone positions of DNA duplexes, allowing ODNP to measure the hydration dynamics at select positions along the DNA surface. With a novel approach to ODNP analysis, we isolate the contributions of water molecules at these sites that undergo free translational diffusion from water molecules that either loosely bind to or exchange protons with the DNA. The results reveal that a significant population of water in a localized volume adjacent to the DNA surface exhibits fast, bulk-like characteristics and moves unusually rapidly compared to water found in similar probe volumes near protein and membrane surfaces. Control studies show that the observation of these characteristics are upheld even when the DNA duplex is tethered to streptavidin or the mobility of the nitroxides is altered. This implies that, as compared to protein or lipid surfaces, it is an intrinsic feature of the DNA duplex surface that it interacts only weakly with a significant fraction of the surface hydration water network. The displacement of this translationally mobile water is energetically less costly than that of more strongly bound water by up to several kBT and thus can lower the activation barrier for interactions involving the DNA surface.

Thermodynamics-hydration relationships within loops that affect G-quadruplexes under molecular crowding conditions

We systematically investigated the effects of loop length on the conformation, thermodynamic stability, and hydration of DNA G-quadruplexes under dilute and molecular crowding conditions in the presence of Na(+). Structural analysis showed that molecular crowding induced conformational switches of oligonucleotides with the longer guanine stretch and the shorter thymine loop. Thermodynamic parameters further demonstrated that the thermodynamic stability of G-quadruplexes increased by increasing the loop length from two to four, whereas it decreased by increasing the loop length from four to six. Interestingly, we found by osmotic pressure analysis that the number of water molecules released from the G-quadruplex decreased with increasing thermodynamic stability. We assumed that base-stacking interactions within the loops not only stabilized the whole G-quadruplex structure but also created hydration sites by accumulating nucleotide functional groups. The molecular crowding effects on the stability of G-quadruplexes composed of abasic sites, which reduce the stacking interactions at the loops, further demonstrated that G-quadruplexes with fewer stacking interactions within the loops released a larger number of water molecules upon folding. These results showed that the stacking interactions within the loops determined the thermodynamic stability and hydration of the whole G-quadruplex.

Characterization of metal ion-nucleic acid interactions in solution

Metal ions are inextricably involved with nucleic acids due to their polyanionic nature. In order to understand the structure and function of RNAs and DNAs, one needs to have detailed pictures on the structural, thermodynamic, and kinetic properties of metal ion interactions with these biomacromolecules. In this review we first compile the physicochemical properties of metal ions found and used in combination with nucleic acids in solution. The main part then describes the various methods developed over the past decades to investigate metal ion binding by nucleic acids in solution. This includes for example hydrolytic and radical cleavage experiments, mutational approaches, as well as kinetic isotope effects. In addition, spectroscopic techniques like EPR, lanthanide(III) luminescence, IR and Raman as well as various NMR methods are summarized. Aside from gaining knowledge about the thermodynamic properties on the metal ion-nucleic acid interactions, especially NMR can be used to extract information on the kinetics of ligand exchange rates of the metal ions applied. The final section deals with the influence of anions, buffers, and the solvent permittivity on the binding equilibria between metal ions and nucleic acids. Little is known on some of these aspects, but it is clear that these three factors have a large influence on the interaction between metal ions and nucleic acids.

Observation of water molecules within the bimolecular d(G₃CT₄G₃C)₂G-quadruplex

G-Rich oligonucleotides with cytosine residues in their sequences can form G-quadruplexes where G-quartets are flanked by G·C Watson-Crick base pairs. In an attempt to probe the role of cations in stabilization of a structural element with two G·C base pairs stacked on a G-quartet, we utilized solution state nuclear magnetic resonance to study the folding of the d(G(3)CT(4)G(3)C) oligonucleotide into a G-quadruplex upon addition of (15)NH(4)(+) ions. Its bimolecular structure exhibits antiparallel strands with edge-type loops. Two G-quartets in the core of the structure are flanked by a couple of Watson-Crick G·C base pairs in a sheared arrangement. The topology is equivalent to the solution state structure of the same oligonucleotide in the presence of Na(+) and K(+) ions [Kettani, A., et al. (1998) J. Mol. Biol.282, 619, and Bouaziz, S., et al. (1998) J. Mol. Biol.282, 637). A single ammonium ion binding site was identified between adjacent G-quartets, but three sites were expected. The remaining potential cation binding sites between G-quartets and G·C base pairs are occupied by water molecules. This is the first observation of long-lived water molecules within a G-quadruplex structure. The flanking G·C base pairs adopt a coplanar arrangement and apparently do not require cations to neutralize unfavorable electrostatic interactions among proximal carbonyl groups. A relatively fast movement of ammonium ions from the inner binding site to bulk with the rate constants of 21 s(-1) was attributed to the lack of hydrogen bonds between adjacent G·C base pairs and the flexibility of the T(4) loops.

Effects of an unnatural base pair replacement on the structure and dynamics of DNA and neighboring water and ions

Incorporating small molecule probes into biomolecular systems to report on local structure and dynamics is a powerful strategy that underlies a wide variety of experimental techniques, including fluorescence, electron paramagnetic resonance (EPR), and Forster resonance energy transfer (FRET) measurements. When an unnatural probe is inserted into a protein or DNA, the degree to which the presence of the probe has perturbed the local structure and dynamics it was intended to study is always an important concern. Here, molecular dynamics (MD) simulations are used to systematically study the effect of replacing a DNA base pair with a fluorescent probe, coumarin 102 deoxyriboside, at six unique sites along an A-tract DNA dodecamer. While the overall structure of the DNA oligonucleotide remains intact, replacement of A*T base pairs leads to widespread structural and dynamic perturbations up to four base pairs away from the probe site, including widening of the minor groove and increased DNA flexibility. New DNA conformations, not observed in the native sequence, are sometimes found in the vicinity of the probe and its partner abasic site analog. Strong correlations are demonstrated between DNA surface topology and water mobility.

Interpreting protein/DNA interactions: distinguishing specific from non-specific and electrostatic from non-electrostatic components

We discuss the effectiveness of existing methods for understanding the forces driving the formation of specific protein-DNA complexes. Theoretical approaches using the Poisson-Boltzmann (PB) equation to analyse interactions between these highly charged macromolecules to form known structures are contrasted with an empirical approach that analyses the effects of salt on the stability of these complexes and assumes that release of counter-ions associated with the free DNA plays the dominant role in their formation. According to this counter-ion condensation (CC) concept, the salt-dependent part of the Gibbs energy of binding, which is defined as the electrostatic component, is fully entropic and its dependence on the salt concentration represents the number of ionic contacts present in the complex. It is shown that although this electrostatic component provides the majority of the Gibbs energy of complex formation and does not depend on the DNA sequence, the salt-independent part of the Gibbs energy--usually regarded as non-electrostatic--is sequence specific. The CC approach thus has considerable practical value for studying protein/DNA complexes, while practical applications of PB analysis have yet to demonstrate their merit.

A Coarse Grained Model for Atomic-Detailed DNA Simulations with Explicit Electrostatics

Coarse-grain (CG) techniques allow considerable extension of the accessible size and time scales in simulations of biological systems. Although many CG representations are available for the most common biomacromolecules, very few have been reported for nucleic acids. Here, we present a CG model for molecular dynamics simulations of DNA on the multi-microsecond time scale. Our model maps the complexity of each nucleotide onto six effective superatoms keeping the "chemical sense" of specific Watson-Crick recognition. Molecular interactions are evaluated using a classical Hamiltonian with explicit electrostatics calculated under the framework of the generalized Born approach. This CG representation is able to accurately reproduce experimental structures, breathing dynamics, and conformational transitions from the A to the B form in double helical fragments. The model achieves a good qualitative reproduction of temperature-driven melting and its dependence on size, ionic strength, and sequence specificity. Reconstruction of atomistic models from CG trajectories give remarkable agreement with structural, dynamic, and energetic features obtained from fully atomistic simulation, opening the possibility to acquire nearly atomic detail data from CG trajectories.

The unique structure of A-tracts and intrinsic DNA bending

Short runs of adenines are a ubiquitous DNA element in regulatory regions of many organisms. When runs of 4–6 adenine base pairs (‘A-tracts’) are repeated with the helical periodicity, they give rise to global curvature of the DNA double helix, which can be macroscopically characterized by anomalously slow migration on polyacrylamide gels. The molecular structure of these DNA tracts is unusual and distinct from that of canonical B-DNA. We review here our current knowledge about the molecular details of A-tract structure and its interaction with sequences flanking them of either side and with the environment. Various molecular models were proposed to describe A-tract structure and how it causes global deflection of the DNA helical axis. We review old and recent findings that enable us to amalgamate the various findings to one model that conforms to the experimental data. Sequences containing phased repeats of A-tracts have from the very beginning been synonymous with global intrinsic DNA bending. In this review, we show that very often it is the unique structure of A-tracts that is at the basis of their widespread occurrence in regulatory regions of many organisms. Thus, the biological importance of A-tracts may often be residing in their distinct structure rather than in the global curvature that they induce on sequences containing them.

Circular dichroism of polynucleotides: Interactions of NiCl2 with poly(dA-dT).poly(dA-dT) and poly(dG-dC).poly(dG-dC) in a water-in-oil microemulsion

The thermal behavior of the synthetic, high molecular weight, double stranded polynucleotides poly(dA-dT).poly(dA-dT) [polyAT] and poly(dG-dC).poly(dG-dC) [polyGC] solubilized in the aqueous core of the quaternary water-in-oil cationic microemulsion CTAB|n-pentanol|n-hexane|water in the presence of increasing amounts of NiCl(2) at several constant ionic strength values (NaCl) has been studied by means of circular dichroism and electronic absorption spectroscopies. In the microemulsive medium, both polynucleotides show temperature-induced modifications that markedly vary with both Ni(II) concentration and ionic strength. An increase of temperature causes denaturation of the polyAT duplex at low nickel concentrations, while more complex CD spectral modifications are observed at higher nickel concentrations and ionic strengths. By contrast, thermal denaturation is never observed for polyGC. At low Ni(II) concentrations, the increase of temperature induces conformational transitions from B-DNA to Z-DNA form, or, more precisely, to left-handed helical structures. In some cases, at higher nickel concentrations, the CD spectra suggest the presence of Z'-type forms of the polynucleotide.