Solution structure of [d(GGTATACC)]2: wrinkled D structure of the TATA moiety

Phase-sensitive two-dimensional nuclear Overhauser effect spectra of [d(GGTATACC)]2 in aqueous deuterium oxide solution at four mixing times were quantified to give all nonoverlapping cross-peak intensities. A structural model for [d(GGTATACC)]2 was built in which the GG- and -CC moieties were in the B-DNA form, while the middle -TATA- moiety was in the wrinkled-D form (BDB model). This model was subjected to energy refinement by molecular mechanics calculations with the program AMBER. Counterions (Na+) were added to neutralize the charges, and water molecules were placed bridging across the minor groove. A complete relaxation matrix analysis was used to calculate two-dimensional nuclear Overhauser effect spectra of [d(GGTATACC)]2 from the above models (before and after energy refinement) and from four other [d(GGTATACC)]2 structural models: regular A, crystalline A, regular B, and energy-minimized B. Among them, the energy-minimized BDB model yielded a set of theoretical spectra that gave the best fit to the experimental spectra. It was also the energetically most stable. Therefore, it is a good representation of the ensemble- and time-averaged structure of the octamer in solution. This model has backbone torsion angles similar to those of B-form DNA in the GG- and -CC moieties and torsion angles similar to those of wrinkled D form DNA in the -TATA- moiety. The base stacking and base pairing are not interrupted at the junctions between the two structural moieties. Its minor groove is narrower than that of B DNA, and the solvent-accessible surface of the minor groove forms a closed hydration tunnel in the middle -TATA- segment.     

Comparison of the conformation of an oligonucleotide containing a central G-T base pair with the non-mismatch sequence by proton NMR.

We have recorded NOESY spectra of two non-selfcomplementary undecanucleotide duplexes. From the observed NOEs we do not detect any significant distortion of the helix when a G-C pair is replaced by a G-T pair and the normal interresidue connectivities can be followed through the mismatch site. We conclude that the 2D spectra of the non-exchangeable protons do not allow differentiation between a wobble or rare tautomer form for the mismatch. NOE measurements in H2O, however, clearly show that the mismatch adopts a wobble structure and give information on the hydration in the minor groove for the G-T base pair which is embedded between two A-T base pairs in the sequence.     

NMR observation of individual molecules of hydration water bound to DNA duplexes: direct evidence for a spine of hydration water present in aqueous solution.

The residence times of individual hydration water molecules in the major and minor grooves of DNA were measured by nuclear magnetic resonance (NMR) spectroscopy in aqueous solutions of d-(CGCGAATTCGCG)2 and d-(AAAAATTTTT)2. The experimental observations were nuclear Overhauser effects (NOE) between water protons and the protons of the DNA. The positive sign of NOEs with the thymine methyl groups shows that the residence times of the hydration water molecules near these protons in the major groove of the DNA must be shorter than about 500 ps, which coincides with the behavior of surface hydration water in peptides and proteins. Negative NOEs were observed with the hydrogen atoms in position 2 of adenine in both duplexes studied. This indicates that a 'spine of hydration' in the minor groove, as observed by X-ray diffraction in DNA crystals, is present also in solution, with residence times significantly longer than 1 ns. Such residence times are reminiscent of 'interior' hydration water molecules in globular proteins, which are an integral part of the molecular architecture both in solution and in crystals.     

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.     

Hydration of DNA in aqueous solution: NMR evidence for a kinetic destabilization of the minor groove hydration of d-(TTAA)2 versus d-(AATT)2 segments.

The residence times of the hydration water molecules near the base protons of d-(GTGGAATTCCAC)2 and d-(GTGGTTAACCAC)2 were investigated by nuclear magnetic resonance (NMR) spectroscopy. Nuclear Overhauser effects (NOE) were observed between base protons of the DNA and hydration water in NOESY and ROESY experiments. Large positive NOESY cross peaks observed between the resonances of the water and the adenine 2H protons of the central d-(AATT)2 segment in the duplex d-(GTGGAATTCCAC)2 indicate the presence of a 'spine of hydration' with water molecules exhibiting residence times on the DNA longer than 1 nanosecond. In contrast, no positive intermolecular NOESY cross peaks were detected in the d-(TTAA)2 segment of the duplex d-(GTGGTTAACCAC)2, indicating that no water molecules bound with similarly long residence times occur in the minor groove of this fragment. These results can be correlated with the larger width of the minor groove in d-(TTAA)2 segments as compared to that in d-(AATT)2 segments, as observed previously in single crystal structures of related oligonucleotide duplexes in B type conformation. The present experiments confirm earlier experimental results from single crystal studies and theoretical predictions that a 5'-dTA-3' step in the nucleotide sequence interrupts the spine of hydration in the minor groove.     

DNA B to D transition can be explained in terms of hydration economy of the minor groove atoms

Adjacent phosphate oxygen atoms in A and Z-DNA are located much closer together than in the B form and can be hydrated more economically due to the formation of water bridges between them, whereas in the B form phosphates are hydrated individually. This principle of hydration economy of phosphate groups discovered by Saenger and colleagues could not be applied to the B-D transition, which, like the B-A and B-Z transitions, occurs in a situation of water deficiency, because the distances between adjacent phosphates of individual polynucleotide chains in the D form are not much different from B-DNA. It follows from our calculations of B and D-DNA accessibility to solvent performed by the method of Lee & Richards, and from a simulation of solvent structure near DNA, that there is an economy of hydration only for the minor groove atoms. This feature and some experimental data can explain why only a limited range of sequences consisting of A.T or I.C pairs undergo the transition to the D form. The conformational transition in DNAs with such sequences to a poly[d(A]).poly[d(T])-like conformation (Bh-DNA), which is accompanied by a narrowing of the minor groove, can be explained in the same way. Calculations suggest that in the D-form minor groove of different A-T or I-C DNAs there is a double-layer hydration spine similar to that observed by Drew & Dickerson in the A-T tract of the d(C-G-C-G-A-A-T-T-C-G-C-G) dodecamer. The B-D and B-Bh transitions in A + T-rich DNAs can have biological implications, e.g. they can facilitate DNA bending upon the interaction with proteins.     

Dynamics and binding mode of Hoechst 33258 to d(GTGGAATTCCAC)2 in the 1:1 solution complex as determined by two-dimensional 1H NMR.

We have investigated the interaction of the bisbenzimidazole derivative Hoechst 33258 with the self-complementary dodecadeoxynucleotide duplex d(GTGGAATTCCAC)2 using one-dimensional (1D) and two-dimensional (2D) proton nuclear magnetic resonance (1H NMR) spectroscopy. To monitor the extent of complex formation, we used the imino proton region of the 1D 1H NMR spectra acquired in H2O solution. These spectra show that the DNA duplex loses its inherent C2v symmetry upon addition of the drug, indicating that the two molecules form a kinetically stable complex on the NMR time scale (the lifetime of the complex has been measured to be around 450 ms). We obtained sequence-specific assignments for all protons of the ligand and most protons of each separate strand of the oligonucleotide duplex using a variety of homonuclear 2D 1H NMR experiments. The aromatic protons of the DNA strands, which are symmetrically related in the free duplex, exhibit exchange cross peaks in the complex. This indicates that the drug binds in two equivalent sites on the 12-mer, with an exchange rate constant of 2.2 +/- 0.2 s-1. Twenty-five intermolecular NOEs were identified, all involving adenine 2 and sugar 1' protons of the DNA and protons in all four residues of the ligand, indicating that Hoechst 33258 is located in the minor groove at the AATT site. Only protons along the same edge of the two benzimidazole moieties of the drug show NOEs to DNA protons at the bottom of the minor groove. Using molecular mechanics, we have generated a unique model of the complex using distance constraints derived from the intermolecular NOEs. We present, however, evidence that the piperazine group may adopt at least two locally different conformations when the drug is bound to this dodecanucleotide.     

Sequence-dependent recognition of DNA duplexes. Netropsin complexation to the AATT site of the d(G-G-A-A-T-T-C-C) duplex in aqueous solution

We have investigated intermolecular interactions and conformational features of the netropsin X d(G-G-A-A-T-T-C-C) complex by one- and two-dimensional NMR studies in aqueous solution. Netropsin removes the 2-fold symmetry of the d(G-G-A-A-T-T-C-C) duplex at the AATT binding site and to a lesser extent at adjacent dG X dC base pairs resulting in doubling of resonances for specific positions in the spectrum of the complex at 25 degrees C. We have assigned the amide, pyrrole, and CH2 protons of netropsin, and the base and sugar H1' protons of the nucleic acid from an analysis of the nuclear Overhauser effect (NOESY) and correlated (COSY) spectra of the complex at 25 degrees C. We observe intermolecular nuclear Overhauser effects (NOE) between all three amide and both pyrrole protons on the concave face of the antibiotic and the minor groove adenosine H2 proton of the two central A4 X T5 base pairs of the d(G1-G2-A3-A4-T5-T6-C7-C8) duplex. Weaker intermolecular NOEs are also observed between the pyrrole concave face protons and the sugar H1' protons of residues T5 and T6 in the AATT minor groove of the duplex. We also detect intermolecular NOEs between the guanidino CH2 protons at one end of netropsin and adenosine H2 proton of the two flanking A3 X T6 base pairs of the octanucleotide duplex. These studies establish a set of intermolecular contacts between the concave face of the antibiotic and the minor groove AATT segment of the d(G-G-A-A-T-T-C-C) duplex in solution. The magnitude of the NOEs require that there be no intervening water molecules sandwiched between the antibiotic and the DNA so that release of the minor groove spine of hydration is a prerequisite for netropsin complex formation.     

Anomalous structure and properties of poly (dA).poly(dT). Computer simulation of the polynucleotide structure with the spine of hydration in the minor groove.

The results of the search for low-energy conformations of poly(dA).poly(dT) and of the poly(dA).poly(dT) "complex" with the spine of hydration similar to that found by Dickerson and co-workers (Kopka, M.L., Fratini, A.V., Drew, H.R. and Dickerson, R.E. (1983) J. Mol. Biol. 163, 129-146) in the minor groove of the CGCGAATTCGCG crystals are described. It is shown that the existence of such a spine in the minor groove of poly(dA).poly(dT) is energetically favourable. Moreover, the spine of hydration makes the polynucleotide conformation similar to the poly(dA).poly(dT) structure in fibers and to the conformation of the central part of CGCGAATTCGCG in crystals; it also acquires features characteristic of the structure of poly(dA).poly(dT) and DNA oligo(dA)-tracts in solution. It is shown that the existence of the TpA step in conformations characteristic of the poly(dA).poly(dT) complex with the spine of hydration is energetically unfavourable (in contrast to the ApT step) and therefore this step should result in destabilization of the spine of hydration in the DNA minor groove. Thus, it appears that the spine of hydration as described by Dickerson and co-workers is unlikely to exist in the poly d(A-T).poly d(A-T) structure. The data obtained permit us to interpret a large body of experimental facts concerning the unusual structure and properties of poly(dA).poly(dT) and oligo(dA)-tracts in DNA both in fibers and in solution. The results provide evidence of the existence of the minor groove spine of hydration both in fibers and in solution on A/T tracts of DNA which do not contain the TpA step. The spine plays an active role in the formation of the anomalous conformation of these tracts.     

Sequence-specific interaction of Hoechst 33258 with the minor groove of an adenine-tract DNA duplex studied in solution by 1H NMR spectroscopy.

The interaction of Hoechst 33258 with the minor groove of the adenine-tract DNA duplex d(CTTTTGCAAAAG)2 has been studied in both D2O and H2O solutions by 1D and 2D 1H NMR spectroscopy. Thirty-one nuclear Overhauser effects between drug and nucleotide protons within the minor groove of the duplex, together with ring-current induced perturbations to the chemical shifts of basepair and deoxyribose protons, define the position and orientation of the bound dye molecules. Two drug molecules bind cooperatively and in symmetry related orientations at the centre of the 5'-TTTT and 5'-AAAA sequences with the binding interactions spanning only the four A-T basepairs. The positively charged N-methylpiperazine moieties point towards the centre of the duplex while the phenol groups are disposed towards the 3'-ends of the sequence. Resonance averaging is apparent for both the D2/D6 and D3/D5 phenol protons and D2"'/D6"' and D3"'/D5"' of the N-methylpiperazine ring and is consistent with these groups being involved in rapid rotation or ring-flipping motions in the bound state. Interstrand NOEs between adenine H2s and deoxyribose H1' are consistent with a high degree of propeller twisting of the A-T basepairs at the binding site of the aromatic benzimidazole and phenol rings of Hoechst. The data imply that the minor groove is particularly narrow with many contacts between the complementary curved surfaces of the drug and DNA indicating that strong van der Waals interactions, involving the floor and the walls of the minor groove, stabilize the complex. In our model the NH groups of the benzimidazole rings are positioned to make a pair of bifurcated hydrogen bonds with the adenine N3 and thymine O2 on the floor of the minor groove.     

Comparison of the conformation of poly(dI-dC) with poly(dI-dbr5C) and the B and Z forms of poly(dG-dC). One- and two-dimensional NMR studies

One- and two-dimensional nuclear Overhauser effects (2D NOE) have been used to compare the conformational properties of 60-80 base pair long duplexes of the synthetic DNA polymer poly(dI-dC) with those of poly(dI-dbr5C) and poly(dG-dC) in the B and Z conformations. Cross peaks in the 2D NOE spectra arising from proton-proton dipolar interactions which are more or less independent of the DNA conformation are used to assign the spectra of these molecules. Other cross peaks are sensitive to the conformational details, and these are used to make deductions about the average conformation in solution. The proton-proton interactions that give rise to the cross peaks in the 2D NOE spectrum of poly(dI-dC) are indicative of a B family conformation and rule out the possibility of some alternative conformations, including A, Z, alternating B, and left-handed B-DNA. The spectra are similar to those obtained from B-form poly(dI-dbr5C) and poly(dG-dC) but different from Z-form poly(dG-dC). Taken together, these results indicate that the solution conformation of poly(dI-dC) is not unusual but more closely resembles that of other B-form DNAs.     

Structure, dynamics and hydration of the nogalamycin-d(ATGCAT)2Complex determined by NMR and molecular dynamics simulations in solution.

The structure of the 1:1 nogalamycin:d(ATGCAT)2 complex has been determined in solution from high-resolution NMR data and restrained molecular dynamics (rMD) simulations using an explicit solvation model. The antibiotic intercalates at the 5'-TpG step with the nogalose lying along the minor groove towards the centre of the duplex. Many drug-DNA nuclear Overhauser enhancements (NOEs) in the minor groove are indicative of hydrophobic interactions over the TGCA sequence. Steric occlusion prevents a second nogalamycin molecule from binding at the symmetry-related 5'-CpA site, leading to the conclusion that the observed binding orientation in this complex is the preferred orientation free of the complication of end-effects (drug molecules occupy terminal intercalation sites in all X-ray structures) or steric interactions between drug molecules (other NMR structures have two drug molecules bound in close proximity), as previously suggested. Fluctuations in key structural parameters such as rise, helical twist, slide, shift, buckle and sugar pucker have been examined from an analysis of the final 500 ps of a 1 ns rMD simulation, and reveal that many sequence-dependent structural features previously identified by comparison of different X-ray structures lie within the range of dynamic fluctuations observed in the MD simulations. Water density calculations on MD simulation data reveal a time-averaged pattern of hydration in both the major and minor groove, in good agreement with the extensive hydration observed in two related X-ray structures in which nogalamycin is bound at terminal 5'-TpG sites. However, the pattern of hydration determined from the sign and magnitude of NOE and ROE cross-peaks to water identified in 2D NOESY and ROESY experiments identifies only a few "bound" water molecules with long residence times. These solvate the charged bicycloaminoglucose sugar ring, suggesting an important role for water molecules in mediating drug-DNA electrostatic interactions within the major groove. The high density of water molecules found in the minor groove in X-ray structures and MD simulations is found to be associated with only weakly bound solvent in solution.     

Investigation of the solution structure of a DNA octamer [d(GGAATTCC)]2 using two-dimensional nuclear Overhauser enhancement spectroscopy

Proton two-dimensional nuclear Overhauser enhancement (2D NOE) spectra in the pure absorption phase were obtained at 500 MHz for [d(GGAATTCC)]2 in aqueous solution at a series of mixing times. The experimental data were analyzed by comparison with theoretical spectra calculated using the complete 70 X 70 relaxation matrix including all proton dipole-dipole interactions and spin diffusion [Keepers, J. W. & James, T. L. (1984) J. Magn. Reson. 57, 404-426]. The theoretical spectra at each mixing time were calculated using two structures: a standard B-form DNA structure and an energy-minimized structure based on the similarity of the six internal residues of the title octamer with those of the dodecamer [d(CGCGAATTCGCG)]2, for which the crystal structure has been determined. Neither the standard B-form nor the energy-minimized structure will yield theoretical 2D NOE spectra which accurately reproduce all peak intensities in the experimental spectra. However, many features of the experimental spectra can be represented by both the B-form and the energy-minimized structure. Sequence-dependent structural characteristics are manifest in the 2D NOE spectra, in particular at the purine-pyrimidine junction as noted previously in the crystal structure. On the whole, the energy-minimized structure appears to yield theoretical 2D NOE spectra which mimic many, if not all, aspects of the experimental spectra. All 2D NOE data were consistent with nanosecond correction times as implied by proton spin-lattice relaxation time measurements. But better fits of some of the 2D NOE data using small variations in an effective isotropic correlation time suggest that there may be some local variations in mobility within the octamer duplex structure in solution.     

Solution structure of [d(GTATATAC)]2 via restrained molecular dynamics simulations with nuclear magnetic resonance constraints derived from relaxation matrix analysis of two-dimensional nuclear Overhauser effect experiments

Two-dimensional nuclear Overhauser effect (2D NOE) spectra have been used as the experimental basis for determining the solution structure of the duplex [d(GTATATAC)]2 employing restrained molecular dynamics (rMD) simulations. The MARDIGRAS algorithm has been employed to construct a set of 233 interproton distance constraints via iterative complete relaxation matrix analysis utilizing the peak intensities from the 2D NOE spectra obtained for different mixing times and model structures. The upper and lower bounds for each of the constraints, defining size of a flat-well potential function term used in the rMD simulations, were conservatively chosen as the largest or smallest value calculated by MARDIGRAS. Three different starting models were utilized in several rMD calculations: energy-minimized A-DNA, B-DNA, and a structure containing wrinkled D-DNA in the interior. Considerable effort was made to define the appropriate force constants to be employed with the NOE terms in the AMBER force field, using as criteria the average constraints deviation, the constraints violation energy and the total energy. Of the 233 constraints, one was generated indirectly, but proved to be crucial in defining the structure: the cross-strand A5-H2 A5-H2 distance. As those two protons resonate isochronously for the self-complementary duplex, the distance cannot be determined directly. However, the general pattern of 2D NOE peak intensities, spin-lattice relaxation time (T1) values, and 31P nuclear magnetic resonance spectra lead to use of the A3-H2 A7-H2 distance for A5-H2 A5-H2 as well. Five rMD runs, with different random number seeds, were made for each of the three starting structures with the full distance constraint set. The average structure from all 15 runs and the five-structure averages from each starting structure were all quite similar. Two rMD runs for each starting structure were made with the A5-H2 A5-H2 constraint missing. The average of these six rMD runs revealed differences in structure, compared to that with the full set of constraints, primarily for the middle two base-pairs involving the missing cross-strand constraint but global deviations also were found. Conformational analysis of the resulting structures revealed that the inner four to six base-pairs differed in structure from the termini. Furthermore, an alternating structure was suggested with features alternating for the A-T and T-A steps.     

Solution structure of a trans-opened (10S)-dA adduct of (+)-(7S,8R,9S,10R)-7,8-dihydroxy-9,10-epoxy-7,8,9,10-tetrahydrobenzo[a]pyrene in a fully complementary DNA duplex: evidence for a major syn conformation

Two-dimensional NMR was used to determine the solution structure of an undecanucleotide duplex, d(CGGTCACGAGG).d(CCTCGTGACCG), in which (+)-(7S,8R,9S,10R)-7,8-dihydroxy-9,10-epoxy-7,8,9,10-tetrahydrobenzo[a]pyrene is covalently bonded to the exocyclic N(6)() amino group of the central deoxyadenosine, dA(6), through trans addition at C10 of the epoxide (to give a 10S adduct). The present study represents the first NMR structure of a benzo[a]pyrene (10S)-dA adduct in DNA with a complementary T opposite the modified dA. Exchangeable and nonexchangeable protons of the modified duplex were assigned by the use of TOCSY (in D(2)O) and NOESY spectra (in H(2)O and D(2)O). Sequential NOEs expected for a B-type DNA conformation with typical Watson-Crick base pairing are observed along the duplex, except at the lesion site. We observed a strong intraresidue NOE cross-peak between H1' and H8 of the modified dA(6). The sugar H2' and H2' ' of dC(5) lacked NOE cross-peaks with H8 of dA(6) but showed weak interactions with H2 of dA(6) instead. In addition, the chemical shift of the H8 proton (7.51 ppm) of dA(6) appears at a higher field than that of H2 (8.48 ppm). These NOE and chemical shift data for the dA(6) base protons are typical of a syn glycosidic bond at the modified base. Restrained molecular dynamics/energy minimization calculations show that the hydrocarbon is intercalated from the major groove on the 3'-side of the modified base between base pairs A(6)-T(17) and C(7)-G(16) and confirm the syn glycosidic angle (58 degrees ) of the modified dA(6). In the syn structure, a weak A-T hydrogen bond is possible between the N3-H proton of T(17) and N7 of dA(6) (at a distance of 3.11 A), whereas N1, the usual hydrogen bonding partner for N3-H of T when dA is in the anti conformation, is 6.31 A away from this proton. The 10(S)-dA modified DNA duplex remains in a right-handed helix, which bends in the direction of the aliphatic ring of BaP at about 42 degrees from the helical axis. ROESY experiments provided evidence for interconversion between the major, syn conformer and a minor, possibly anti, conformer.     

1H NMR studies on the interaction between distamycin A and a symmetrical DNA dodecamer

High-resolution NMR techniques have been used to examine the structural and dynamical features of the interaction between distamycin A and the self-complementary DNA dodecamer duplex d-(CGCGAATTCGCG)2. The proton resonances of d(CGCGAATTCGCG)2 have been completely assigned by previous two-dimensional NMR studies [Hare, D. R., Wemmer, D. E., Chou, S. H., Drobny, G., & Reid, B. R. (1983) J. Mol. Biol. 171, 319-336]. Addition of the asymmetric drug molecule to the symmetric dodecamer leads to the formation of an asymmetric complex as evidenced by a doubling of DNA resonances over much of the spectrum. In two-dimensional exchange experiments, strong cross-peaks were observed between uncomplexed DNA and drug-bound DNA resonances, permitting direct assignment of many drug-bound DNA resonances from previously assigned free DNA resonances. Weaker exchange cross-peaks between formerly symmetry related DNA resonances indicate that the drug molecule flips head-to-tail on one duplex with half the frequency at which it leaves the DNA molecule completely. In experiments performed in H2O, nuclear Overhauser effects (NOEs) were observed from each drug amide proton to an adenine C2H and a pyrrole H3 ring proton. In two-dimensional nuclear Overhauser experiments performed on D2O solutions, strong intermolecular NOEs were observed between each of the three pyrrole H3 resonances of the drug and an adenine C2H resonance, with weaker NOEs observed between the drug H3 resonances and C1'H resonances. The combined NOE data allow us to position the distamycin A unambiguously on the DNA dodecamer, with the drug spanning the central AATT segment in the minor groove.     

Complicated water orientations in the minor groove of the B-DNA decamer d(CCATTAATGG)2 observed by neutron diffraction measurements

It has long been suspected that the structure and function of a DNA duplex can be strongly dependent on its degree of hydration. By neutron diffraction experiments, we have succeeded in determining most of the hydrogen (H) and deuterium (D) atomic positions in the decameric d(CCATTAATGG)₂ duplex. Moreover, the D positions in 27 D₂O molecules have been determined. In particular, the complex water network in the minor groove has been observed in detail. By a combined structural analysis using 2.0 Å resolution X-ray and 3.0 Å resolution neutron data, it is clear that the spine of hydration is built up, not only by a simple hexagonal hydration pattern (as reported in earlier X-ray studies), but also by many other water bridges hydrogen-bonded to the DNA strands. The complexity of the hydration pattern in the minor groove is derived from an extraordinary variety of orientations displayed by the water molecules.     

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.     

Conversions of poly(2-aminodeoxyadenylate-5-halodeoxyuridylate) from B to A forms in high salt. An NMR and circular dichroism study

Proton one- and two-dimensional nuclear Overhauser enhancement (1D and 2D NOE) spectroscopy has been used to demonstrate that poly(d2NH2A-d5IU) and poly(d2NH2A-d5BrU) are converted from the B to the A conformation in high salt, as found previously for poly(d2NH2A-dT) [Borah, B., Cohen, J. S., Howard, F. B., & Miles, H. T. (1985) Biochemistry 24, 7456-7462]. The 2D NOE and 1D NOE spectra exhibit strong base proton (H8,H6)-H3' cross relaxation, suggesting short interproton distances. These results are indicative of a C3'-endo sugar pucker for both purine and pyrimidine residues in an A or closely related structure. The circular dichroism and UV spectra are consistent with the interpretation of an A conformation in high salt.     

1H NMR and optical spectroscopic investigation of the sequence-dependent dimerization of a symmetrical cyanine dye in the DNA minor groove

A symmetrical cyanine dye was previously shown to bind as a cofacial dimer to alternating A-T sequences of duplex DNA. Indirect evidence suggested that dimerization of the dye occurred in the minor groove. 1H NMR experiments reported here verify this model based on broadening and shifting of signals due to protons on carbon 2 of adenine and imino protons at the central five A-T pairs of the 11 base pair duplex: 5'-GCGTATATGCG-3'/3'-CGCATATACGC-5'. This binding mode is similar to that of distamycin A, even though the dye lacks the hydrogen-bonding groups used by distamycin for sequence-specific recognition. Surprisingly, the third base pair (G-C) was also implicated in the binding site. UV-vis experiments were used to compare the extent of dimerization of the dye for 11 different sequence variants. These experiments verified the importance of a G-C pair at the third position: replacing this pair with A-T suppressed dimerization. These results indicate that the dye binding site spans six base pairs: 5'-GTATAT-3'. The initial G-C pair seems to be important for widening the minor groove rather than for making important contacts with the dye molecules since inverting its orientation to C-G or replacing it with I-C still led to favorable dimerization of the dye.