Determination of the residence time of water molecules hydrating B'- DNA and B-DNA, by one-dimensional zero-enhancement nuclear Overhauser effect spectroscopy

Terahertz spectroscopy as a method for investigation of hydration shells of biomolecules

The hydration of biomolecules is one of the fundamental processes underlying the construction of living matter. The formation of the native conformation of most biomolecules is possible only in an aqueous environment. At the same time, not only water affects the structure of biomolecules, but also biomolecules affect the structure of water, forming hydration shells. However, the study of the structure of biomolecules is given much more attention than their hydration shells. A real breakthrough in the study of hydration occurred with the development of the THz spectroscopy method, which showed that the hydration shell of biomolecules is not limited to 1-2 layers of strongly bound water, but also includes more distant areas of hydration with altered molecular dynamics. This review examines the fundamental features of the THz frequency range as a source of information about the structural and dynamic characteristics of water that change during hydration. The applied approaches to the study of hydration shells of biomolecules based on THz spectroscopy are described. The data on the hydration of biomolecules of all main types obtained from the beginning of the application of THz spectroscopy to the present are summarized. The emphasis is placed on the possible participation of extended hydration shells in the realization of the biological functions of biomolecules and at the same time on the insufficient knowledge of their structural and dynamic characteristics.

The Ubiquity of High-Energy Nanosecond Fluorescence in DNA Duplexes

During the past few years, several studies reported that a significant part of the intrinsic fluorescence of DNA duplexes decays with surprisingly long lifetimes (1-3 ns) at wavelengths shorter than the ππ* emission of their monomeric constituents. This high-energy nanosecond emission (HENE), hardly discernible in the steady-state fluorescence spectra of most duplexes, was investigated by time-correlated single-photon counting. The ubiquity of HENE contrasts with the paradigm that the longest-lived excited states correspond to low-energy excimers/exciplexes. Interestingly, the latter were found to decay faster than the HENE. So far, the excited states responsible for HENE remain elusive. In order to foster future studies for their characterization, this Perspective presents a critical summary of the experimental observations and the first theoretical approaches. Moreover, some new directions for further work are outlined. Finally, the obvious need for computations of the fluorescence anisotropy considering the dynamic conformational landscape of duplexes is stressed.

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.

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.

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.

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.

Structural study of the Fox-1 RRM protein hydration reveals a role for key water molecules in RRM-RNA recognition

The Fox-1 RNA recognition motif (RRM) domain is an important member of the RRM protein family. We report a 1.8 Å X-ray structure of the free Fox-1 containing six distinct monomers. We use this and the nuclear magnetic resonance (NMR) structure of the Fox-1 protein/RNA complex for molecular dynamics (MD) analyses of the structured hydration. The individual monomers of the X-ray structure show diverse hydration patterns, however, MD excellently reproduces the most occupied hydration sites. Simulations of the protein/RNA complex show hydration consistent with the isolated protein complemented by hydration sites specific to the protein/RNA interface. MD predicts intricate hydration sites with water-binding times extending up to hundreds of nanoseconds. We characterize two of them using NMR spectroscopy, RNA binding with switchSENSE and free-energy calculations of mutant proteins. Both hydration sites are experimentally confirmed and their abolishment reduces the binding free-energy. A quantitative agreement between theory and experiment is achieved for the S155A substitution but not for the S122A mutant. The S155 hydration site is evolutionarily conserved within the RRM domains. In conclusion, MD is an effective tool for predicting and interpreting the hydration patterns of protein/RNA complexes. Hydration is not easily detectable in NMR experiments but can affect stability of protein/RNA complexes.

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

Sequence dependency of canonical base pair opening in the DNA double helix

The flipping-out of a DNA base from the double helical structure is a key step of many cellular processes, such as DNA replication, modification and repair. Base pair opening is the first step of base flipping and the exact mechanism is still not well understood. We investigate sequence effects on base pair opening using extensive classical molecular dynamics simulations targeting the opening of 11 different canonical base pairs in two DNA sequences. Two popular biomolecular force fields are applied. To enhance sampling and calculate free energies, we bias the simulation along a simple distance coordinate using a newly developed adaptive sampling algorithm. The simulation is guided back and forth along the coordinate, allowing for multiple opening pathways. We compare the calculated free energies with those from an NMR study and check assumptions of the model used for interpreting the NMR data. Our results further show that the neighboring sequence is an important factor for the opening free energy, but also indicates that other sequence effects may play a role. All base pairs are observed to have a propensity for opening toward the major groove. The preferred opening base is cytosine for GC base pairs, while for AT there is sequence dependent competition between the two bases. For AT opening, we identify two non-canonical base pair interactions contributing to a local minimum in the free energy profile. For both AT and CG we observe long-lived interactions with water and with sodium ions at specific sites on the open base pair.

The DNA double helix, a molecule that stores biological information, has become an iconic image of biomedical research. In order to use or repair the information it carries, the bases that are stacked in the helix need to be chemically exposed. This can happen either by separating the two strands in the helix or by flipping out individual bases. Here, we focus on the latter process. Usually proteins are involved in interactions with bases, but it is still unclear if bases are pulled out actively by proteins or if they act on spontaneously flipped bases. Although experiments can detect base pair opening, it is difficult to detect which base moves in which direction. Here, we present results from molecular dynamics simulations using a recently developed sampling method which improves the statistics in the simulations by enhancing the probability of the base pair opening event. We observe differences in probability, modes and mechanism of opening that depend not only on the types of the bases in the pair, but also strongly on their neighbors. This provides essential information for understanding how DNA functions.

Rotational dynamics of water molecules near biological surfaces with implications for nuclear quadrupole relaxation

Based on Molecular Dynamics simulations of two different systems, the protein ubiquitin dissolved in water and an AOT reverse micelle, we present a broad analysis of the single particle rotational dynamics of water. A comprehensive connection to NQR, which is a prominent experimental method in this field, is developed, based on a reformulation of its theoretical framework. Interpretation of experimental NQR results requires a model which usually assumes that the NQR experiences retardation only in the first hydration shell. Indeed, the present study shows that this first-shell model is correct. Moreover, previous experimental retardation factors are quantitatively reproduced. All of this is seemingly contradicted by results of other methods, e.g., dielectric spectroscopy, responsible for a long-standing debate in this field. Our detailed analysis shows that NQR omits important information contained in overall water dynamics, most notably, the retardation of the water dipole axis in the electric field exerted by a biological surface.

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.

Decelerated water dynamics and vibrational couplings of hydrated DNA mapped by two-dimensional infrared spectroscopy

Double-stranded DNA oligomers containing 23 alternating adenine-thymine base pairs are studied at different hydration levels by femtosecond two-dimensional (2D) infrared spectrosopy. Coupled NH stretching modes of the A-T pairs and OH stretching excitations of the water shell are discerned in the 2D spectra. Limited changes of NH stretching frequencies and line shapes with increasing hydration suggest spectral dynamics governed by DNA rather than water fluctuations. In contrast, OH stretching excitations of the water shell around fully hydrated DNA undergo spectral diffusion on a ~500 fs time scale. The center line slopes of the 2D spectra of hydrated DNA demonstrate a slower decay of the frequency-time correlation function (TCF) than that in neat water, as is evident from a comparison with 2D spectra of neat H(2)O and theoretical TCFs. We attribute this behavior to reduced structural fluctuations of the water shell and a reduced rate of resonant OH stretching energy transfer.

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.

Ultrafast energy exchange via water-phosphate interactions in hydrated DNA

The ionic phosphate groups in the DNA backbone play a key role for DNA hydration. We study ultrafast vibrational dynamics and local interactions of phosphate groups and water by femtosecond two-color pump-probe spectroscopy. The asymmetric (PO(2))(-) stretching vibration nu(AS)(PO(2))(-) of artificial DNA oligomers containing 23 alternating adenine-thymine base pairs displays a lifetime of 340 fs, independent of the hydration level. For DNA at zero relative humidity, excess energy from the decay of the phosphate excitation is transferred within DNA on a 20 ps time scale. For fully hydrated DNA, the water shells around the phosphates serve as a primary heat sink accepting vibrational excess energy from DNA on a femtosecond time scale. OH stretching excitation of water molecules around fully hydrated DNA induces an ultrafast nu(AS)(PO(2))(-) response which includes rearrangements of the hydration shell and a reduction of the average number of phosphate-water hydrogen bonds.

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.

Hydration pattern of A4T4 and T4A4 DNA: A molecular dynamics study

Hydration pattern and energetics of 'A-tract' containing duplexes have been studied using molecular dynamics on 12-mer self-complementary sequences 5'-d(GCA4T4GC)-3' and 5'-d(CGT4A4CG)-3'. The structural features for the simulated duplexes showed correlation with the corresponding experimental structures. Analysis of the hydration pattern confirmed that water network around the simulated duplexes is more conformation specific rather than sequence specific. The calculated heat capacity change upon duplex formation showed that the process is entropically driven for both the sequences. Furthermore, the theoretical free energy estimates calculated using MMPBSA approach showed a higher net electrostatic contribution for A4T4 duplex formation than for T4A4, however, energetically both the duplexes are observed to be equally stable.

Sensitivity of hydrogen bonds of DNA and RNA to hydration, as gauged by 1JNH measurements in ethanol-water mixtures

Hydrogen-bond lengths of nucleic acids are (1) longer in DNA than in RNA, and (2) sequence dependent. The physicochemical basis for these variations in hydrogen-bond lengths is unknown, however. Here, the notion that hydration plays a significant role in nucleic acid hydrogen-bond lengths is tested. Watson-Crick N1...N3 hydrogen-bond lengths of several DNA and RNA duplexes are gauged using imino 1J(NH) measurements, and ethanol is used as a cosolvent to lower water activity. We find that 1J(NH) values of DNA and RNA become less negative with added ethanol, which suggests that mild dehydration reduces hydrogen-bond lengths even as the overall thermal stabilities of these duplexes decrease. The 1J(NH) of DNA are increased in 8 mol% ethanol to those of RNA in water, which suggests that the greater hydration of DNA plays a significant role in its longer hydrogen bonds. The data also suggest that ethanol-induced dehydration is greater for the more hydrated G:C base pairs and thereby results in greater hydrogen-bond shortening than for the less hydrated A:T/U base pairs of DNA and RNA.

High-resolution protein hydration NMR experiments: Probing how protein surfaces interact with water and other non-covalent ligands

High-resolution solution NMR experiments are extremely useful to characterize the location and the dynamics of hydrating water molecules at atomic resolution. However, these methods are severely limited by undesired incoherent transfer pathways such as those arising from exchange-relayed intra-molecular cross-relaxation. Here, we review several complementary exchange network editing methods that can be used in conjunction with other types of NMR hydration experiments such as magnetic relaxation dispersion and 1J(NC') measurements to circumvent these limitations. We also review several recent contributions illustrating how the original solution hydration NMR pulse sequence architecture has inspired new approaches to map other types of non-covalent interactions going well beyond the initial scope of hydration. Specifically, we will show how hydration NMR methods have evolved and have been adapted to binding site mapping, ligand screening, protein-peptide and peptide-lipid interaction profiling.

Large-amplitude fast motions in double-stranded DNA driven by solvent thermal fluctuations

The nature of the internal dynamics of double-stranded DNA in aqueous environment remains to be established. We consider the motions to stem from thermal fluctuations/dissipations of the harmonic modes of beads (bases and sugars) in a cylindrical geometry that are tracked through the stochastic Langevin trajectories; these are characterized by parameters obtained from published data. The present approach has allowed a comparative study of the dynamics for DNA lengths in the range of 20-600 base pairs. For this range, we find that rotational motions about directions parallel to the helix axis (opening, twist) and perpendicular to it (propeller-twist, roll) contribute significantly to the dynamics. For a 20-mer at a solvent viscosity of 1 cP, the calculated fluorescence anisotropy profile exhibits a fast decay in the subnanosecond range due to large-amplitude fluctuations at the mesoscopic level. This feature reproduces the experimental behavior well, and suggests a possible way for the initiation of biological processes: they may be suddenly triggered on this scale through the occurrence of favorable thermal fluctuations. This analysis also reveals that, as is the case for a 20-mer, the dynamics of longer N-mers are dominated by internal motions, and are modulated by the viscosity of the solvent, in agreement with our previous experimental observations. Moreover, the model indicates that occurrence of partially concerted rotations of the bases due to thermal fluctuations can possibly be sustained over a DNA length of the order of 100 A at 1 ns, suggesting a possible mechanism for action-at-a-distance in transcription.

Effects of small neutral osmolytes on the supercoiling free energy and intrinsic twist of p30? DNA

Both theory and experiments are employed to investigate the effects of small neutral osmolytes on the average intrinsic twist (l0), the torsion and bending elastic constants, and the twist energy parameter (ET) that governs the supercoiling free energy. The experimental data for ethylene glycol and acetamide at 37 degrees C suggest, and are interpreted in terms of, a model wherein the DNA exhibits an equilibrium between two distinct conformational states that possess different numbers of bound water molecules and exhibit different intrinsic twists and torsion and bending elastic constants. Expressions are derived to relate the effective ET and l0 to the equilibrium constant, water activity (aw), and number (n) of bound water molecules released per cooperative domain undergoing the two-state transition. The variations of l0 and ET with -ln(aw) are similar for acetamide and ethylene glycol at 37 degrees C. Fitting the theory to those data yields the range n = 103-125 for ethylene glycol and n = 71-113 for acetamide, depending on the assumed value of ET for the dehydrated state. The cooperative domain size of the two-state transition is estimated to exceed about 25-30 base pairs (bp). Between 0 and 19.4 w/v % ethylene glycol, the torsion elastic constant, measured by time-resolved fluorescence polarization anisotropy (FPA), increases by 1.37-fold, whereas the measured ET decreases by 1.15-fold over that same range. The implied decrease in bending rigidity over that range is by a factor of about 0.7. The variations of l0 and ET with increasing -ln(aw) due to added ethylene glycol at 37 degrees C are far smaller than the corresponding variations observed previously at 14 and 15 degrees C. However, at 21 degrees C, upon adding either ethylene glycol or acetamide, l0 and ET initially decline steeply with increasing -ln(aw), with slopes possibly comparable to those seen at 14 and 15 degrees C, but then flatten out and follow curves similar to those at 37 degrees C. Possible origins of such mixed behavior are discussed. The effects of betaine at both 37 and 21 degrees C differ qualitatively and quantitatively in various respects from those of ethylene glycol and acetamide. Upon adding sucrose, l0 initially jumps to higher plateaus at both 37 and 21 degrees C, but its effects on ET cannot be reliably assessed, due to the limited range of -ln(aw).

Water-mediated contacts in thetrp-repressor operator complex recognition process

Water-mediated contacts are known as an important recognition tool in trp-repressor operator systems. One of these contacts involves two conserved base pairs (G(6).C(-6) and A(5). T(-5)) and three amino acids (Lys 72, Ile 79, and Ala 80). To investigate the nature of these contacts, we analyzed the X-ray structure (PDB code: 1TRO) of the trp-repressor operator complex by means of molecular dynamics simulations. This X-ray structure contains two dimers that exhibit structural differences. From these two different starting structures, two 10 ns molecular dynamics simulations have been performed. Both of our simulations show an increase of water molecules in the major groove at one side of the dimer, while the other side remains unchanged compared to the X-ray structure. Though the maximum residence time of the concerned water molecules decreases with an increase of solvent at the interface, these water molecules continue to play an important role in mediating DNA-protein contacts. This is shown by new stable amino acids-DNA distances and a long water residence time compared to free DNA simulation. To maintain stability of the new contacts, the preferential water binding site on O6(G6) is extended. This extension agrees with mutation experiment data on A5 and G6, which shows different relative affinity due to mutation on these bases [A. Joachimiak, T. E. Haran, P. B. Sigler, EMBO Journal 1994, Vol. 13, No. (2) pp. 367-372]. Due to the rearrangements in the system, the phosphate of the base G6 is able to interconvert to the B(II) substate, which is not observed on the other half side of the complex. The decrease of the number of hydrogen bonds between protein and DNA backbone could be the initial step of the dissociation process of the complex, or in other words an intermediate complex conformation of the association process. Thus, we surmise that these features show the importance of water-mediated contacts in the trp-repressor operator recognition process.

Base pair opening within B-DNA: free energy pathways for GC and AT pairs from umbrella sampling simulations

The conformational pathways and the free energy variations for base opening into the major and minor grooves of a B-DNA duplex are studied using umbrella sampling molecular dynamics simulations. We compare both GC and AT base pair opening within a double-stranded d(GAGAGAGAGAGAG)* d(CTCTCTCTCTCTC) oligomer, and we are also able to study the impact of opening on the conformational and dynamic properties of DNA and on the surrounding solvent. The results indicate a two-stage opening process with an initial coupling of the movements of the bases within the perturbed base pair. Major and minor groove pathways are energetically comparable in the case of the pyrimidine bases, but the major groove pathway is favored for the larger purine bases. Base opening is coupled to changes in specific backbone dihedrals and certain helical distortions, including untwisting and bending, although all these effects are dependent on the particular base involved. Partial opening also leads to well defined water bridging sites, which may play a role in stabilizing the perturbed base pairs.

Melting of the solvent structure around a RNA duplex: a molecular dynamics simulation study

From three 2.4-ns molecular dynamics simulations of the r(CpG)(12) duplex conducted at 5, 25 and 37 degrees C, a strong temperature dependence of the dynamics of the water molecules and ions located in the first nucleic acid coordination shell is observed. At 5 degrees C, the highest residence times of bound water molecules exceed 1 ns while, at 37 degrees C, they decrease to 0.5 ns in agreement with available NMR data. Similar temperature dependencies are observed for the potassium ions bound to the duplex. In this temperature range, the structure of the RNA helix remains essentially unchanged. Thus, the observed alterations correspond to a 'premelting' of the solvent structure around the duplex. It is proposed that, before the nucleic acid structure melts, the entropy of the solvent increases to a point where it is no longer compensated by the enthalpic contribution of solute-solute and solute-solvent interactions. At this stage, the weakest structural elements start to melt. In other terms, the experimentally observed melting processes are preceded by a melting of the more labile solvent structure.

Water and ion binding around r(UpA)12 and d(TpA)12 oligomers--comparison with RNA and DNA (CpG)12 duplexes

The structural and dynamic properties of the water and ion first coordination shell of the r(A-U) and d(A-T) base-pairs embedded within the r(UpA)12 and d(TpA)12 duplexes are described on the basis of two 2.4 ns molecular dynamics simulations performed in a neutralizing aqueous environment with 0.25 M added KCl. The results are compared to previous molecular dynamics simulations of the r(CpG)12 and d(CpG)12 structures performed under similar conditions. It can be concluded that: (i) RNA helices are more rigid than DNA helices of identical sequence, as reflected by the fact that RNA duplexes keep their initial A-form shape while DNA duplexes adopt more sequence-specific shapes. (ii) Around these base-pairs, the water molecules occupy 21 to 22 well-defined hydration sites, some of which are partially occupied by potassium ions. (iii) These hydration sites are occupied by an average of 21.9, 21.0, 20.1, and 19.8 solvent molecules (water and ions) around the r(G=C), r(A-U), d(G=C), and d(A-T) pairs, respectively. (iv) From a dynamic point of view, the stability of the hydration shell is the strongest for the r(G=C) pairs and the weakest for the d(A-T) pairs. (v) For RNA, the observed long-lived hydration patterns are essentially non-sequence dependent and involve water bridges located in the deep groove and linking OR atoms of adjacent phosphate groups. Maximum lifetimes are close to 400 ps. (vi) In contrast, for DNA, long-lived hydration patterns are sequence dependent and located in the minor groove. For d(CpG)12, water bridges linking the (G)N3 and (C)O2 with the O4' atoms of adjacent nucleotides with 400 ps maximum lifetimes are characterized while no such bridges are observed for d(TpA)12. (vii) Potassium ions are observed to bind preferentially to deep/major groove atoms at RpY steps, essentially d(GpC), r(GpC), and r(ApU), by forming ion-bridges between electronegative atoms of adjacent base-pairs. On average, about half an ion is observed per base-pair. Positive ion-binding determinants are related to the proximity of two or more electronegative atoms. Negative binding determinants are associated with the electrostatic and steric hindrance due to the proximity of electropositive amino groups and neutral methyl groups. Potassium ions form only transient contacts with phosphate groups.

Molecular dynamics simulation reveals conformational switching of water-mediated uracil-cytosine base-pairs in an RNA duplex

A 4 ns molecular dynamics simulation of an RNA duplex (r-GGACUUCGGUCC)(2 )in solution with Na+ and Cl- as counterions was performed. The X-ray structure of this duplex includes two water-mediated uracil-cytosine pairs. In contrast to the other base-pairs in the duplex the water-mediated pairs switch between different conformations. One conformation corresponds to the geometry of the water-mediated UC pairs in the duplex X-ray structure with water acting both as hydrogen-bond donor and acceptor. Another conformation is close to that of a water-mediated UC base-pair found in the X-ray structure of the 23 S rRNA sarcin/ricin domain. In this case the oxygen of the water molecule is linked to two-base donor sites. For a very short time also a direct UC base-pair and a further conformation that is similar to the one found in the RNA duplex structure but exhibits an increased H3(U)...N3(C) distance is observed. Water molecules with unusually long residence times are involved in the water-mediated conformations. These results indicate that the dynamic behaviour of the water-mediated UC base-pairs differs from that of the duplex Watson-Crick and non-canonical guanine-uracil pairs with two or three direct hydrogen bonds. The conformational variability and increased flexibility has to be taken into account when considering these base-pairs as RNA building blocks and as recognition motifs.

Structure and hydration of BamHI DNA recognition site: a molecular dynamics investigation

The results of a 3-ns molecular dynamics simulation of the dodecamer duplex d(TATGGATCCATA)(2) recognized by the BamHI endonuclease are presented here. The DNA has been simulated as a flexible molecule using an AMBER force field and the Ewald summation method, which eliminates the undesired effects of truncation and permits evaluation of the full effects of electrostatic forces. The starting B conformation evolves toward a configuration quite close to that observed through x-ray diffraction in its complex with BamHI. This configuration is fairly stable and the Watson-Crick hydrogen bonds are well maintained over the simulation trajectory. Hydration analysis indicates a preferential hydration for the phosphate rather than for the ester oxygens. Hydration shells in both the major and minor groove were observed. In both grooves the C-G pairs were found to be more hydrated than A-T pairs. The "spine of hydration" in the minor groove was clear. Water residence times are longer in the minor groove than in the major groove, although relatively short in both cases. No special long values are observed for sites where water molecules were observed by x-ray diffraction, indicating that water molecules having a high probability of being located in a specific site are also fast-exchanging.

A-Tract bending: insights into experimental structures by computational models

While solution structures of adenine tract (A-tract) oligomers have indicated a unique bend direction equivalent to negative global roll (commonly termed "minor-groove bending"), crystallographic data have not unambiguously characterized the bend direction; nevertheless, many features are shared by all A-tract crystal and solution structures (e.g. propeller twisting, narrow minor grooves, and localized water spines). To examine the origin of bending and to relate findings to the crystallographic and solution data, we analyze molecular dynamics trajectories of two solvated A-tract dodecamers: 1D89, d(CGCGA(6)CG), and 1D98, d(CGCA(6)GCG), using a new general global bending framework for analyzing bent DNA and DNA/protein complexes. It is significant that the crystallographically-based initial structures are converted from dissimilar to similar bend directions equivalent to negative global roll, with the average helical-axis bend ranging from 10.5 degrees to 14.1 degrees. The largest bend occurs as positive roll of 12 degrees on the 5' side of the A-tracts (supporting a junction model) and is reinforced by gradual curvature at each A-tract base-pair (bp) step (supporting a wedge model). The precise magnitude of the bend is subtly sequence dependent (consistent with a curved general sequence model). The conversion to negative global roll only requires small local changes at each bp, accumulated over flexible moieties both outside and inside the A-tract. In contrast, the control sequence 1BNA, d(CGCGA(2)TTCGCG), bends marginally (only 6.9 degrees ) with no preferred direction. The molecular features that stabilize the bend direction in the A-tract dodecamers include propeller twisting of AT base-pairs, puckering differences between A and T deoxyriboses, a narrow minor groove, and a stable water spine (that extends slightly beyond the A-tract, with lifetimes approaching 0.2 ns). The sugar conformations, in particular, are proposed as important factors that support bent DNA. It is significant that all these curvature-stabilizing features are also observed in the crystallographic structures, but yield overall different bending paths, largely due to the effects of sequences outside the A-tract. These results merge structural details reported for A-tract structures by experiment and theory and lead to structural and dynamic insights into sequence-dependent DNA flexibility, as highlighted by the effect of an A-tract variant of a TATA-box element on bending and flexibility required for TBP binding.