Water and monovalent ions in the minor groove of B-DNA oligonucleotides as seen by NMR

Solvent Vibrations as a Proxy of the Telomere G-Quadruplex Rearrangements across Thermal Unfolding

G-quadruplexes (G4s) are noncanonical forms of DNA involved in many key genome functions. Here, we exploited UV Resonance Raman scattering to simultaneously explore the vibrational behavior of a human telomeric G4 (Tel22) and its aqueous solvent as the biomolecule underwent thermal melting. We found that the OH stretching band, related to the local hydrogen-bonded network of a water molecule, was in strict relation with the vibrational features of the G4 structure as a function of temperature. In particular, the modifications to the tetrahedral ordering of the water network were strongly coupled to the DNA rearrangements, showing changes in temperature that mirrored the multi-step melting process of Tel22. The comparison between circular dichroism and Raman results supported this view. The present findings provide novel insights into the impact of the molecular environment on G4 conformation. Improving current knowledge on the solvent structural properties will also contribute to a better understanding of the role played by water arrangement in the complexation of G4s with ligands.

Long-range DNA-water interactions

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

Probing the Hydrogen-Bonding Environment of Individual Bases in DNA Duplexes with Isotope-Edited Infrared Spectroscopy

Measuring the strength of the hydrogen bonds between DNA base pairs is of vital importance for understanding how our genetic code is physically accessed and recognized in cells, particularly during replication and transcription. Therefore, it is important to develop probes for these key hydrogen bonds (H-bonds) that dictate events critical to cellular function, such as the localized melting of DNA. The vibrations of carbonyl bonds are well-known probes of their H-bonding environment, and their signals can be observed with infrared (IR) spectroscopy. Yet, pinpointing a single bond of interest in the complex IR spectrum of DNA is challenging due to the large number of carbonyl signals that overlap with each other. Here, we develop a method using isotope editing and infrared (IR) spectroscopy to isolate IR signals from the thymine (T) C2=O carbonyl. We use solvatochromatic studies to show that the TC2=O signal’s position in the IR spectrum is sensitive to the H-bonding capacity of the solvent. Our results indicate that C2=O of a single T base within DNA duplexes experiences weak H-bonding interactions. This finding is consistent with the existence of a third, noncanonical CH···O H-bond between adenine and thymine in both Watson–Crick and Hoogsteen base pairs in DNA.

Unravelling the microhydration frameworks of prototype PAH by infrared spectroscopy: naphthalene–(water)1–3

Microhydration structures of the prototypical PAH, naphthalene, are probed by IR spectroscopy in helium droplets. The sequential water addition produces an extended hydrogen-bonded hydration network bound via π hydrogen bond to the aromatic ring.

The Role of Correlation and Solvation in Ion Interactions with B-DNA

The ionic atmospheres around nucleic acids play important roles in biological function. Large-scale explicit solvent simulations coupled to experimental assays such as anomalous small-angle x-ray scattering can provide important insights into the structure and energetics of such atmospheres but are time- and resource intensive. In this article, we use classical density functional theory to explore the balance among ion-DNA, ion-water, and ion-ion interactions in ionic atmospheres of RbCl, SrCl2, and CoHexCl3 (cobalt hexamine chloride) around a B-form DNA molecule. The accuracy of the classical density functional theory calculations was assessed by comparison between simulated and experimental anomalous small-angle x-ray scattering curves, demonstrating that an accurate model should take into account ion-ion correlation and ion hydration forces, DNA topology, and the discrete distribution of charges on the DNA backbone. As expected, these calculations revealed significant differences among monovalent, divalent, and trivalent cation distributions around DNA. Approximately half of the DNA-bound Rb(+) ions penetrate into the minor groove of the DNA and half adsorb on the DNA backbone. The fraction of cations in the minor groove decreases for the larger Sr(2+) ions and becomes zero for CoHex(3+) ions, which all adsorb on the DNA backbone. The distribution of CoHex(3+) ions is mainly determined by Coulomb and steric interactions, while ion-correlation forces play a central role in the monovalent Rb(+) distribution and a combination of ion-correlation and hydration forces affect the Sr(2+) distribution around DNA. This does not imply that correlations in CoHex solutions are weaker or stronger than for other ions. Steric inaccessibility of the grooves to large CoHex ions leads to their binding at the DNA surface. In this binding mode, first-order electrostatic interactions (Coulomb) dominate the overall binding energy as evidenced by low sensitivity of ionic distribution to the presence or absence of second-order electrostatic correlation interactions.

Competition among Li(+), Na(+), K(+), and Rb(+) monovalent ions for DNA in molecular dynamics simulations using the additive CHARMM36 and Drude polarizable force fields

In the present study we report on interactions of and competition between monovalent ions for two DNA sequences in MD simulations. Efforts included the development and validation of parameters for interactions among the first-group monovalent cations, Li(+), Na(+), K(+), and Rb(+), and DNA in the Drude polarizable and additive CHARMM36 force fields (FF). The optimization process targeted gas-phase QM interaction energies of various model compounds with ions and osmotic pressures of bulk electrolyte solutions of chemically relevant ions. The optimized ionic parameters are validated against counterion condensation theory and buffer exchange-atomic emission spectroscopy measurements providing quantitative data on the competitive association of different monovalent ions with DNA. Comparison between experimental and MD simulation results demonstrates that, compared to the additive CHARMM36 model, the Drude FF provides an improved description of the general features of the ionic atmosphere around DNA and leads to closer agreement with experiment on the ionic competition within the ion atmosphere. Results indicate the importance of extended simulation systems on the order of 25 Å beyond the DNA surface to obtain proper convergence of ion distributions.

Balancing the interactions of ions, water, and DNA in the Drude polarizable force field

Recently we presented a first-generation all-atom Drude polarizable force field for DNA based on the classical Drude oscillator model, focusing on optimization of key dihedral angles followed by extensive validation of the force field parameters. Presently, we describe the procedure for balancing the electrostatic interactions between ions, water, and DNA as required for development of the Drude force field for DNA. The proper balance of these interactions is shown to impact DNA stability and subtler conformational properties, including the conformational equilibrium between the BI and BII states, and the A and B forms of DNA. The parametrization efforts were simultaneously guided by gas-phase quantum mechanics (QM) data on small model compounds and condensed-phase experimental data on the hydration and osmotic properties of biologically relevant ions and their solutions, as well as theoretical predictions for ionic distribution around DNA oligomer. In addition, fine-tuning of the internal base parameters was performed to obtain the final DNA model. Notably, the Drude model is shown to more accurately reproduce counterion condensation theory predictions of DNA charge neutralization by the condensed ions as compared to the CHARMM36 additive DNA force field, indicating an improved physical description of the forces dictating the ionic solvation of DNA due to the explicit treatment of electronic polarizability. In combination with the polarizable DNA force field, the availability of Drude polarizable parameters for proteins, lipids, and carbohydrates will allow for simulation studies of heterogeneous biological systems.

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.

NMR cross-correlated relaxation rates reveal ion coordination sites in DNA

In this work, a novel NMR method for the identification of preferential coordination sites between physiologically relevant counterions and nucleic acid bases is demonstrated. In this approach, the NMR cross-correlated relaxation rates between the aromatic carbon chemical shift anisotropy and the proton-carbon dipolar interaction are monitored as a function of increasing Na(+), K(+), and Mg(2+) concentrations. Increasing the counterion concentration modulates the residence times of the counterions at specific sites around the nucleic acid bases. It is demonstrated that the modulation of the counterion concentration leads to sizable variations of the cross-correlated relaxation rates, which can be used to probe the site-specific counterion coordination. In parallel, the very same measurements report on the rotational tumbling of DNA, which, as shown here, depends on the nature of the ion and its concentration. This methodology is highly sensitive and easily implemented. The method can be used to cross-validate and/or complement direct but artifact-prone experimental techniques such as X-ray diffraction, NMR analysis with substitutionary ions, and molecular dynamics simulations. The feasibility of this technique is demonstrated on the extraordinarily stable DNA mini-hairpin d(GCGAAGC).

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.

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.

Effect of molecular crowding and ionic strength on the isothermal hybridization of oligonucleotides

The isothermal hybridization of complementary oligonucleotides, 15-mer, 25-mer, 35-mer, and a molecular beacon, was investigated under varying conditions of molecular crowding and ionic strength, using hypochromicity to follow strand pairing and polyethylene glycol as a crowding agent. Thermodynamic analysis of the results revealed the addition of counterions to the oligonucleotide backbones, DeltaPsi, to be dependent on the strand GC content and the molecular crowding. A decrease in DeltaPsi was observed, with both increasing GC% and solution PEG content. In contrast, the number of bound water molecules depended on the activity of Na(+), where two regimes were observed. At a(Na(+)) < 0.05 and increasing molecular crowding, water molecules were released into the DNA solutions, and oligonucleotide pairing was favored with both increasing hydrophobic forces, whereas at a(Na(+)) >or= 0.05, water molecules were bound to the strands, and the extent of double strand formation decreased with increasing PEG wt %.

Facilitated assignment of adenine H2 resonances in oligonucleotides using homonuclear long-range couplings

The small homonuclear coupling constant between H2 and H8 protons of adenine (0.5 +/- 0.2 Hz) was determined and exploited to assign the H2 proton of adenine residues in DNA oligonucleotides. The utility of this approach is demonstrated for two different DNA sequences in D(2)O, for which we show a rapid and unambiguous assignment of the A H2 resonances.

Monovalent cation binding in the minor groove of DNA A-tracts

The binding of five different monovalent cations to DNA oligomers containing A-tracts, runs of four or more contiguous adenine residues, has been assessed by capillary electrophoresis, using the Replacement Ion method. In this method, a nonbinding cation in the background electrolyte is gradually replaced by a binding cation, keeping the ionic strength of the solution constant. Monovalent cation binding reduces the effective charge of an A-tract-containing oligomer, decreasing its free solution mobility. The cations bind in the A-tract minor groove, because the binding site can be blocked by the minor groove binding drug netropsin. Li(+), NH(4)(+), and Tris(+) ions bind to A-tracts with similar affinities; the binding of Na(+) ions is weaker, and K(+) ion binding is highly variable. Each A-tract appears to bind one monovalent cation upon saturation of the binding site(s). For a given cation, the apparent dissociation constants depend on A-tract sequence and orientation, but not on the phasing of the A-tracts with respect to the helix repeat. Differences in the cooperativity of binding of the various cations to A-tracts with different sequences suggest that monovalent cation binding may be coupled with a conformational transition leading to the formation of the characteristic narrow minor groove A-tract structure.

Dynamics of water and ions near DNA: comparison of simulation to time-resolved stokes-shift experiments

Time-resolved Stokes-shift experiments measure the dynamics of biomolecules and of the perturbed solvent near them on subnanosecond time scales, but molecular dynamics simulations are needed to provide a clear interpretation of the results. Here we show that simulations using standard methods quantitatively reproduce the main features of TRSS experiments in DNA and provide a molecular assignment for the dynamics. The simulations reproduce the magnitude and unusual power-law dynamics of the Stokes shift seen in recent experiments [ Andreatta, D., et al. J. Am. Chem. Soc. 2005, 127, 7270 ]. A polarization model is introduced to eliminate cross-correlations between the different components contributing to the signal. Using this model, well-defined contributions of the DNA, water, and counterion to the experimental signal are extracted. Water is found to have the largest contribution and to be responsible for the power-law dynamics. The counterions have a smaller, but non-negligible, contribution with a time constant of 220 ps. The contribution to the signal of the DNA itself is minor and fits a 30 ps stretched exponential. Both time-averaged and dynamic distributions are calculated. They show a small subset of ions with a different coupling but no other evidence of substates or rate heterogeneity.

Internal sodium ions and water molecules in guanine quadruplexes: magnetic relaxation dispersion studies of [d(G3T4G3)]2 and [d(G4T4G4)]2

The structural stability of guanine quadruplexes depends critically on an unusual configuration of dehydrated Na (+) or K (+) ions, closely spaced along the central axis of the quadruplex. Crystallography and NMR spectroscopy indicate that these internal ions can be located between the G-quartet planes as well as in the thymine loops, but the precise ion coordination has been firmly established in only a few cases. Here, we examine the bimolecular diagonal-looped foldback quadruplexes [d(G 3T 4G 3)] 2 (Q3) and [d(G 4T 4G 4)] 2 (Q4) by (2)H, (17)O, and (23)Na magnetic relaxation dispersion (MRD). The MRD data indicate that both quadruplexes contain Na (+) ions between the T 4 loops and the terminal G-quartets and that these ions have one water ligand. These ions exchange with external ions on a time scale of 10-60 mus at 27 degrees C, while their highly ordered water ligands have residence times in the range 10 (-8)-10 (-6) s. The MRD data indicate that Q4 contains three Na (+) ions in the stem sites, in agreement with previous solid-state (23)Na NMR findings but contrary to the only crystal structure of this quadruplex. For Q3, the MRD data suggest a less symmetric coordination of the two stem ions. In both quadruplexes, the stem ions have residence times of 0.6-1.0 ms at 27 degrees C. The equilibrium constant for Na (+) --> K (+) exchange is approximately 4 for both loop and stem sites in Q3, in agreement with previous (1)H NMR findings.

Conductivity of natural and modified DNA measured by scanning tunneling microscopy. The effect of sequence, charge and stacking

The conductivity of DNA covalently bonded to a gold surface was studied by means of the STM technique. Various single- and double-stranded 32-nucleotide-long DNA sequences were measured under ambient conditions so as to provide a better understanding of the complex process of charge-carrier transport in natural as well as chemically modified DNA molecules. The investigations focused on the role of several features of DNA structure, namely the role of the negative charge at the backbone phosphate group and the related complex effects of counterions, and of the stacking interactions between the bases in Watson-Crick and other types of base pairs. The measurements have indicated that the best conductor is DNA in its biologically most relevant double-stranded form with Watson-Crick base pairs and charged phosphates equilibrated with counterions and water. All the studied modifications, including DNA with non-Watson-Crick base pairs, the abasic form, and especially the form with phosphate charges eliminated by chemical modifications, lower the conductivity of natural DNA.

Determination of alkali and halide monovalent ion parameters for use in explicitly solvated biomolecular simulations

Alkali (Li(+), Na(+), K(+), Rb(+), and Cs(+)) and halide (F(-), Cl(-), Br(-), and I(-)) ions play an important role in many biological phenomena, roles that range from stabilization of biomolecular structure, to influence on biomolecular dynamics, to key physiological influence on homeostasis and signaling. To properly model ionic interaction and stability in atomistic simulations of biomolecular structure, dynamics, folding, catalysis, and function, an accurate model or representation of the monovalent ions is critically necessary. A good model needs to simultaneously reproduce many properties of ions, including their structure, dynamics, solvation, and moreover both the interactions of these ions with each other in the crystal and in solution and the interactions of ions with other molecules. At present, the best force fields for biomolecules employ a simple additive, nonpolarizable, and pairwise potential for atomic interaction. In this work, we describe our efforts to build better models of the monovalent ions within the pairwise Coulombic and 6-12 Lennard-Jones framework, where the models are tuned to balance crystal and solution properties in Ewald simulations with specific choices of well-known water models. Although it has been clearly demonstrated that truly accurate treatments of ions will require inclusion of nonadditivity and polarizability (particularly with the anions) and ultimately even a quantum mechanical treatment, our goal was to simply push the limits of the additive treatments to see if a balanced model could be created. The applied methodology is general and can be extended to other ions and to polarizable force-field models. Our starting point centered on observations from long simulations of biomolecules in salt solution with the AMBER force fields where salt crystals formed well below their solubility limit. The likely cause of the artifact in the AMBER parameters relates to the naive mixing of the Smith and Dang chloride parameters with AMBER-adapted Aqvist cation parameters. To provide a more appropriate balance, we reoptimized the parameters of the Lennard-Jones potential for the ions and specific choices of water models. To validate and optimize the parameters, we calculated hydration free energies of the solvated ions and also lattice energies (LE) and lattice constants (LC) of alkali halide salt crystals. This is the first effort that systematically scans across the Lennard-Jones space (well depth and radius) while balancing ion properties like LE and LC across all pair combinations of the alkali ions and halide ions. The optimization across the entire monovalent series avoids systematic deviations. The ion parameters developed, optimized, and characterized were targeted for use with some of the most commonly used rigid and nonpolarizable water models, specifically TIP3P, TIP4P EW, and SPC/E. In addition to well reproducing the solution and crystal properties, the new ion parameters well reproduce binding energies of the ions to water and the radii of the first hydration shells.

Electrostatic, Steric, and Hydration Interactions Favor Na+ Condensation around DNA Compared with K+

Condensation of monovalent counterions around DNA influences polymer properties of the DNA chain. For example, the Na(+) ions show markedly stronger propensity to induce multiple DNA chains to assemble into compact structures compared with the K(+) ions. To investigate the similarities and differences in the sodium and potassium ion condensation around DNA, we carried out a number of extensive all-atom molecular dynamics simulations of a DNA oligomer consisting of 16 base pairs, [d(CGAGGTTTAAACCTCG)](2), in explicit water. We found that the Na(+) ions penetrate the DNA interior and condense around the DNA exterior to a significantly larger degree compared with the K(+) ions. We have provided a microscopic explanation for the larger Na(+) affinity toward DNA that is based on a combination of steric, electrostatic, and hydration effects. Unexpectedly, we found that the Cl(-) co-ions provide more efficient electrostatic screening for the K(+) ions than for the Na(+) ions, contributing to the larger Na(+) condensation around DNA. To examine the importance of the discrete nature of water and ions, we also computed the counterion distributions from the mean-field electrostatic theory, demonstrating significant disagreements with the all-atom simulations. Prior experimental results on the relative extent of the Na(+) and K(+) condensation around DNA were somewhat contradictory. Recent DNA compaction experiments may be interpreted to suggest stronger Na(+) condensation around DNA compared to K(+), which is consistent with our simulations. We also provide a simple interpretation for the experimentally observed increase in DNA electrophoretic mobility in the alkali metal series, Li(+) < Na(+) < K(+) < Rb(+). We compare the DNA segment conformational preferences in various buffers with the proposed NMR models.

Similarities and differences in interaction of K+ and Na+ with condensed ordered DNA. A molecular dynamics computer simulation study

Four 20 ns molecular dynamics simulations have been performed with two counterions, K⁺ or Na⁺, at two water contents, 15 or 20 H₂O per nucleotide. A hexagonal simulation cell comprised of three identical DNA decamers [d(5′-ATGCAGTCAG) × d(5′-TGACTGCATC)] with periodic boundary condition along the DNA helix was used. The simulation setup mimics the DNA state in oriented DNA fibers or in crystals of DNA oligomers. Variation of counterion nature and water content do not alter averaged DNA structure. K⁺ and Na⁺ binding to DNA are different. K⁺ binds to the electronegative sites of DNA bases in the major and the minor grooves, while Na⁺ interacts preferentially with the phosphate groups. Increase of water causes a shift of both K⁺ and Na⁺ from the first hydration shell of O1P/O2P and of the DNA bases in the minor groove with lesser influence for the cation binding to the bases in the major groove. Mobility of both water and cations in the K–DNA systems is faster than in the Na–DNA systems: Na⁺ organizes and immobilizes water structure around itself and near DNA while for K⁺ water is less organized and more dynamic.

Competitive Na(+) and Rb(+) binding in the minor groove of DNA

Sequence-dependent coordination of alkali ions to the nucleotide bases in the minor groove of AT-tract B-DNA has recently been inferred from X-ray crystallography, solution NMR and computer simulations. Here, we present new (23)Na and (87)Rb magnetic relaxation dispersion (MRD) data that demonstrate competitive and long-lived binding of Na(+) and Rb(+) ions in the minor groove of the B-DNA duplex [d(CGCGAATTCGCG)](2). The Na(+)/Rb(+) selectivity of the minor groove is found to be weak, consistent with local structural flexibility. The ion occupancies derived from the MRD data are substantially higher than previously reported, suggesting that groove-bound ions significantly influence the energetics and structural polymorphism of DNA in vivo. For example, in the presence of 0.20 M Na(+) and 0.56 M Rb(+) at 4 degrees C, the ApT site in the minor groove is occupied by a Rb(+) ion, a Na(+) ion, or a water molecule 40, 10, and 50% of the time, respectively. In the absence of Rb(+), the Na(+) occupancy increases to 50%. At 4 degrees C, the mean residence time of groove-bound ions is 0.2 +/- 0.1 micros for Rb(+) and 10 ns to 100 micros for Na(+). A shorter correlation time of 2 ns is attributed to counterions bridging cross-strand phosphate groups.

Dynamics of protein and peptide hydration

Biological processes often involve the surfaces of proteins, where the structural and dynamic properties of the aqueous solvent are modified. Information about the dynamics of protein hydration can be obtained by measuring the magnetic relaxation dispersion (MRD) of the water (2)H and (17)O nuclei or by recording the nuclear Overhauser effect (NOE) between water and protein protons. Here, we use the MRD method to study the hydration of the cyclic peptide oxytocin and the globular protein BPTI in deeply supercooled solutions. The results provide a detailed characterization of water dynamics in the hydration layer at the surface of these biomolecules. More than 95% of the water molecules in contact with the biomolecular surface are found to be no more than two-fold motionally retarded as compared to bulk water. In contrast to small nonpolar molecules, the retardation factor for BPTI showed little or no temperature dependence, suggesting that the exposed nonpolar residues do not induce clathrate-like hydrophobic hydration structures. New NOE data for oxytocin and published NOE data for BPTI were analyzed, and a mutually consistent interpretation of MRD and NOE results was achieved with the aid of a new theory of intermolecular dipolar relaxation that accounts explicitly for the dynamic perturbation at the biomolecular surface. The analysis indicates that water-protein NOEs are dominated by long-range dipolar couplings to bulk water, unless the monitored protein proton is near a partly or fully buried hydration site where the water molecule has a long residence time.

Site- and sequence-selective ultrafast hydration of DNA

Water molecules in the DNA grooves are critical for maintaining structural integrity, conformational changes, and molecular recognition. Here we report studies of site- and sequence-specific hydration dynamics, using 2-aminopurine (Ap) as the intrinsic fluorescence probe and with femtosecond resolution. The dodecamer d[CGCA(Ap)ATTTGCG]2 was investigated, and we also examined the effect of a specific minor groove-binding drug, pentamidine, on hydration dynamics. Two time scales were observed: approximately 1 ps (bulk-like) and 10-12 ps (weakly bound type), consistent with layer hydration observed in proteins and DNA. However, for denatured DNA, the cosolvent condition of 40% formamide hydration is very different: it becomes that of bulk (in the presence of formamide). Well known electron transfer between Ap and nearby bases in stacked assemblies becomes inefficient in the single-stranded state. The rigidity of Ap in the single strands is significantly higher than that in bulk water and that attached to deoxyribose, suggesting a unique role for the dynamics of the phosphate-sugar-base in helix formation. The disparity in minor and major groove hydration is evident because of the site selection of Ap and in the time scale observed here (in the presence and absence of the drug), which is different by a factor of 2 from that observed in the minor groove-drug recognition.

Water at DNA surfaces: ultrafast dynamics in minor groove recognition

Water molecules at the surface of DNA are critical to its equilibrium structure, DNA-protein function, and DNA-ligand recognition. Here we report direct probing of the dynamics of hydration, with femtosecond resolution, at the surface of a DNA dodecamer duplex whose native structure remains unperturbed on recognition in minor groove binding with the bisbenzimide drug (Hoechst 33258). By following the temporal evolution of fluorescence, we observed two well separated hydration times, 1.4 and 19 ps, whereas in bulk water the same drug is hydrated with time constants of 0.2 and 1.2 ps. For comparison, we also studied calf thymus DNA for which the hydration exhibits similar time scales to that of dodecamer DNA. However, the time-resolved polarization anisotropy is very different for the two types of DNA and clearly elucidates the rigidity in drug binding and difference in DNA rotational motions. These results demonstrate that hydration at the surface of the groove is a dynamical process with two general types of trajectories; the slowest of them (approximately 20 ps) are those describing dynamically ordered water. Because of their ultrafast time scale, the "ordered" water molecules are the most weakly bound and are accordingly involved in the entropic (hydration/dehydration) process of recognition.

Water structure theory and some implications for drug design

The development of theories of water structure has been hindered in the past by the difficulty of experimental measurement. Both measurement and computer modelling studies have now reached the stage where theoretical treatments of water structure are converging to a broadly acceptable model. In current understanding, water is a mixture of randomly hydrogen-bonded molecules and larger structures comprised of tetrahedral oxygen centres which, when hydrogen-bonded to each other, lead to five-membered and other rings which can aggregate to form three-dimensional structures. Evidence is taken from studies of the ices, from clathrates and other solid solutions, as well as from liquid solutions, that certain motifs occur very frequently and have relatively high stability, such as the (H2O)20 cavity-forming structure known from studies on clathrates. The implications of recent models of water structure for an understanding of biological events, including the interactions of drugs with receptors, are profound. It is becoming clear that modelling of aqueous solutions of any molecule must consider the explicit interactions with water molecules, which should not be regarded as a continuum: water itself is not a continuum. Solute molecules which possess hydrogen-bonding groups will provoke the formation of further hydrogen-bonding chains of water molecules: if these can form rings, such ringswilltend to persist longerthan chains, givingthesolute a secondary identity of associated water which may play a role in molecular recognition. Solutes that do not have hydrogen-bonding capability, or regions of solutes which are non-polar, may also produce partial cage-like water structures that are characteristic of the solute. The classification of many solutes as structure makers or structure breakers has relevance to the interactions between ligands and large biomolecules such as proteins. While it is generally accepted that sulfate and urea, respectively structure maker and breaker, may alter protein conformation through effects on water, it has not been recognised that bioactive ligands, which also change the conformation of proteins, may do so by a related, but more selective, mechanism. Very early studies of cell contents suggested that the associated water might be different from bulk water, a concept that lost support in the mid-20th century. Current theories of water structure may invite a reappraisal of this position, given the observation that structuring may extend for many molecular diameters from an ordered surface.

DNA-cation interactions: quo vadis?

The interactions between double helical DNA and cations, specifically mono- and divalent metal ions, have recently received increased attention. Molecular dynamics simulations, solution NMR, and X-ray crystallography have all shed light on the coordination of ions in the major and minor grooves of DNA. Metal ion interactions may play key roles in the control of DNA conformation and topology, but despite progress in locating the ions and determining their precise binding modes, it remains difficult to figure out just how important ions really are. What have we learned and what remains to be done?