Hydration and solution structure of d(CGCAAATTTGCG)2 and its complex with propamidine from NMR and molecular modelling

Single-Base Lesions and Mismatches Alter the Backbone Conformational Dynamics in DNA

DNA damage has been implicated in numerous human diseases, particularly cancer, and the aging process. Single-base lesions and mismatches in DNA can be cytotoxic or mutagenic and are recognized by a DNA glycosylase during the process of base excision repair. Altered local dynamics and conformational properties in damaged DNAs have previously been suggested to assist in recognition and specificity. Herein, we use solution nuclear magnetic resonance to quantify changes in BI-BII backbone conformational dynamics due to the presence of single-base lesions in DNA, including uracil, dihydrouracil, 1,N6-ethenoadenine, and T:G mismatches. Stepwise changes to the %BII and ΔG of the BI-BII dynamic equilibrium compared to those of unmodified sequences were observed. Additionally, the equilibrium skews toward endothermicity for the phosphates nearest the lesion/mismatched base pair. Finally, the phosphates with the greatest alterations correlate with those most relevant to the repair of enzyme binding. All of these results suggest local conformational rearrangement of the DNA backbone may play a role in lesion recognition by repair enzymes.

Bound Compound, Interfacial Water, and Phenyl Ring Rotation Dynamics of a Compound in the DNA Minor Groove

DB2277, a heterocyclic diamidine, is a successful design for mixed base pair (bp) DNA sequence recognition. The compound has a central aza-benzimidazole group that forms two H-bonds with a GC bp that has flanking AT bps. The nuclear magnetic resonance structure of the DB2277-DNA complex with an AAGATA recognition site sequence was determined, and here we report extended molecular dynamics (MD) simulations of the structure. DB2277 has two terminal phenyl-amidine groups, one of which is directly linked to the DB2277 heterocyclic core and the other through a flexible -OCH₂- group. The flexibly linked phenyl is too far from the minor groove floor to make direct H-bonds but is linked to an AT bp through water-mediated H-bonds. The flexibly linked phenyl-amidine with water-mediated H-bonds to the bases at the floor of the minor groove suggested that it might rotate in time spans accessible in MD. To test this idea, we conducted multimicrosecond MD simulations to determine if these phenyl rotations could be observed for a bound compound. In a 3 μs simulation, highly dynamic torsional motions were observed for the -OCH₂-linked phenyl but not for the other phenyl. The dynamics periodically reached a level to allow 180° rotation of the phenyl while it was still bound in the minor groove. This is the first observation of rotation of a phenyl bound to DNA, and the results provide mechanistic details about how a rotation can occur as well as how mixed bp recognition can occur for monomer compounds bound to the minor groove.

Uncovering the thermodynamics of monomer binding for RNA replication

The nonenzymatic replication of primordial RNA is thought to have been a critical step in the origin of life. However, despite decades of effort, the poor rate and fidelity of model template copying reactions have thus far prevented an experimental demonstration of nonenzymatic RNA replication. The overall rate and fidelity of template copying depend, in part, on the affinity of free ribonucleotides to the RNA primer-template complex. We have now used (1)H NMR spectroscopy to directly measure the thermodynamic association constants, Kas, of the standard ribonucleotide monophosphates (rNMPs) to native RNA primer-template complexes. The binding affinities of rNMPs to duplexes with a complementary single-nucleotide overhang follow the order C > G > A > U. Notably, these monomers bind more strongly to RNA primer-template complexes than to the analogous DNA complexes. The relative binding affinities of the rNMPs for complementary RNA primer-template complexes are in good quantitative agreement with the predictions of a nearest-neighbor analysis. With respect to G:U wobble base-pairing, we find that the binding of rGMP to a primer-template complex with a 5'-U overhang is approximately 10-fold weaker than to the complementary 5'-C overhang. We also find that the binding of rGMP is only about 2-fold weaker than the binding of rAMP to 5'-U, consistent with the poor fidelity observed in the nonenzymatic copying of U residues in RNA templates. The accurate Ka measurements for ribonucleotides obtained in this study will be useful for designing higher fidelity, more effective RNA replication systems.

Z-formamidoximes in molecular folding and macrocycles

The formamidoxime configurational Z isomer coupled with the pyridylbiscarboxamide conformational codon were used to fold planar, curved structures. When embedded into macrocycles, this folded motif promotes dimerization through π-π stacking and hydrogen-bonding and the formation of tubules akin to molecular channels in the solid state.

Alkoxyamine-derived formamidines: configurational control and molecular folding

N,N'-Disubstituted formamidines, and amidines in general, have very rich configurational, conformational, and tautomeric diversities. As part of an effort to incorporate alkoxyamine-derived formamidine units into foldamers, the first evidence for the isolation of the up-to-now unknown E isomer, the conditions for its exclusive formation, its stability and self-assembly properties, and its configurational isomerization to its much more common Z counterpart are reported. Considering the distinctly different H-bonding patterns displayed by both E and Z isomers, such configurational control may find applications in self-assembly, molecular recognition, and biomimetic systems.

Sequence dependent femtosecond-resolved hydration dynamics in the minor groove of DNA and histone-DNA complexes

Understanding the sequence dependent molecular recognition of DNA is crucial for the rational design of many drugs. Femtosecond resolved studies on the hydration dynamics of the dodecamer duplexes having sequences (CGCGAATTCGCG)2 and (CGCAAATTTGCG)2, and their complexes with the nucleic protein histone 1 (H1) reveal significant correlation of the molecular recognition of the DNA and DNA-protein complexes with the dynamics of hydration. The different molecular recognition of DNA and DNA-protein complexes is also borne out by circular dichroism (CD) and fluorescence detected CD measurements.

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 and ion binding around RNA and DNA (C,G) oligomers

The dynamics, hydration, and ion-binding features of two duplexes, the A(r(CG)(12)) and the B(d(CG)(12)), in a neutralizing aqueous environment with 0.25 M added KCl have been investigated by molecular dynamics (MD) simulations. The regular repeats of the same C=G base-pair motif have been exploited as a statistical alternative to long MD simulations in order to extend the sampling of the conformational space. The trajectories demonstrate the larger flexibility of DNA compared to RNA helices. This flexibility results in less well defined hydration patterns around the DNA than around the RNA backbone atoms. Yet, 22 hydration sites are clearly characterized around both nucleic acid structures. With additional results from MD simulations, the following hydration scale for C=G pairs can be deduced: A-DNA<RNA (+3 H(2)O) and B-DNA<RNA (+2 H(2)O). The calculated residence times of water molecules in the first hydration shell of the helices range from 0.5 to 1 ns, in good agreement with available experimental data. Such water molecules are essentially found in the vicinity of the phosphate groups and in the DNA minor groove. The calculated number of ions that break into the first hydration shell of the nucleic acids is close to 0.5 per base-pair for both RNA and DNA. These ions form contacts essentially with the oxygen atoms of the phosphate groups and with the guanine N7 and O6 atoms; they display residence times in the deep/major groove approaching 500 ps. Further, a significant sequence-dependent effect on ion binding has been noted. Despite slight structural differences, K(+) binds essentially to GpC and not to CpG steps. These results may be of importance for understanding various sequence-dependent binding affinities. Additionally, the data help to rationalize the experimentally observed differences in gel electrophoretic mobility between RNA and DNA as due to the difference in hydration (two water molecules in favor of RNA) rather than to strong ion-binding features, which are largely similar for both nucleic acid structures.

The solution structure of [d(CGC)r(aaa)d(TTTGCG)](2): hybrid junctions flanked by DNA duplexes

The solution structure and hydration of the chimeric duplex [d(CGC)r(aaa)d(TTTGCG)](2), in which the central hybrid segment is flanked by DNA duplexes at both ends, was determined using two-dimensional NMR, simulated annealing and restrained molecular dynamics. The solution structure of this chimeric duplex differs from the previously determined X-ray structure of the analogous B-DNA duplex [d(CGCAAATTTGCG)](2)as well as NMR structure of the analogous A-RNA duplex [r(cgcaaauuugcg)](2). Long-lived water molecules with correlation time tau(c)longer than 0.3 ns were found close to the RNA adenine H2 and H1' protons in the hybrid segment. A possible long-lived water molecule was also detected close to the methyl group of 7T in the RNA-DNA junction but not with the other two thymines (8T and 9T). This result correlates with the structural studies that only DNA residue 7T in the RNA-DNA junction adopts an O4'-endo sugar conformation, while the other DNA residues including 3C in the DNA-RNA junction, adopt C1'-exo or C2'-endo conformations. The exchange rates for RNA C2'-OH were found to be approximately 5-20 s(-1). This slow exchange rate may be due to the narrow minor groove width of [d(CGC)r(aaa)d(TTTGCG)](2), which may trap the water molecules and restrict the dynamic motion of hydroxyl protons. The minor groove width of [d(CGC)r(aaa)d(TTTGCG)](2)is wider than its B-DNA analog but narrower than that of the A-RNA analog. It was further confirmed by its titration with the minor groove binding drug distamycin. A possible 2:1 binding mode was found by the titration experiments, suggesting that this chimeric duplex contains a wider minor groove than its B-DNA analog but still narrow enough to hold two distamycin molecules. These distinct structural features and hydration patterns of this chimeric duplex provide a molecular basis for further understanding the structure and recognition of DNA. RNA hybrid and chimeric duplexes.

Hydration of [d(CGC)r(aaa)d(TTTGCG)](2)

We have studied the hydration and dynamics of RNA C2'-OH in a DNA. RNA hybrid chimeric duplex [d(CGC)r(aaa)d(TTTGCG)](2). Long-lived water molecules with correlation time tau(c) larger than 0.3 ns were found close to the RNA adenine H2 and H1' protons in the hybrid segment. A possible long-lived water molecule was also detected close to the methyl group of 7T in the RNA-DNA junction but not to the other two thymine bases (8T and 9T). This result correlates with the structural studies that only DNA residue 7T in the RNA-DNA junction adopts an O4'-endo sugar conformation (intermediate between B-form and A-form), while the other DNA residues including 3C in the DNA-RNA junction, adopt C1'-exo or C2'-endo conformations (in the B-form domain). Based on the NOE cross-peak patterns, we have found that RNA C2'-OH tends to orient toward the O3' direction, forming a possible hydrogen bond with the 3'-phosphate group. The exchange rates for RNA C2'-OH were found to be around 5-20 s(-1), compared to 26.7(+/-13.8) s(-1) reported previously for the other DNA.RNA hybrid duplex. This slow exchange rate may be due to the narrow minor groove width of [d(CGC)r(aaa)d(TTTGCG)](2), which may trap the water molecules and restrict the dynamic motion of hydroxyl protons. The distinct hydration patterns of the RNA adenine H2 and H1' protons and the DNA 7T methyl group in the hybrid segment, as well as the orientation and dynamics of the RNA C2'-OH protons, may provide a molecular basis for further understanding the structure and recognition of DNA.RNA hybrid and chimeric duplexes.

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.

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

The residence time of water in the minor groove of the d(CGCGAATTCGCG) duplex has been determined by a recent measurement combining nuclear Overhauser enhancements (NOE, ROE) and 17O relaxation dispersion. The time is in the range of nanoseconds, so that it may be measured by a rather simple method proposed here, namely the choice of conditions such that the NOE between the observed DNA proton and a nearby water proton is zero. This condition is realized when the residence time of the water molecule is 0.178 times the nuclear magnetic resonance period (e.g. 0.297 ns at 600 MHz). It may be achieved by varying the magnetic field and/or the temperature. The zero-NOE measurement may be performed by one-dimensional NMR, and has therefore good sensitivity. We have developed excitation sequences which suppress two spurious contributions to the NOE: from neighboring exchangeable protons and from H3' protons whose chemical shift is close to that of water. The method is applied here to the comparison of residence times of water next to B-DNA and next to B'-DNA, the latter corresponding to better stacked, propeller-twisted base-pairs and a correspondingly narrower minor groove. In the minor groove of [d(CGCGAATTCGCG)]2, a B'-DNA duplex, the residence time of the water molecule next to H2 of adenine(6) (underlined), is 0.6 ns at 10 degreesC, in good agreement with the value obtained previously. The residence time is slightly but distinctly shorter for the water next to A5, suggesting non-cooperative departure of these two molecules which are presumed to be part of the hydration spine. Near A5 and A4 of [d(AAAAATTTTT)]2, another B'-DNA duplex, the residence times are approximately twice as long, but the activation enthalpies are about the same, ca. 38 kJ/mol. The residence time in the minor groove of the regular B-DNA sequence d(CGCGATCGCG) was 0.3 ns at 10 degreesC, shorter than in the case of the B'-DNA sequences by factors of 2 and 4, respectively. The temperature dependence is less, with an activation enthalpy of 27 kJ/mol. The major groove residence times are comparable for the three sequences, and a few times shorter than those of minor groove water. A value of 0.36 ns, or even more in case of rotation of water, is obtained around -8 degreesC. The most striking aspect of these results is the relatively small difference in the residence times of reputedly fast and slow-exchanging water molecules bound to DNA in biological conditions. This suggests that the spine of hydration is perhaps not a major stabilizer of the B'-DNA structure as compared with B-DNA.

An NMR study of d(CTACTGCTTTAG).d(CTAAAGCAGTAG) showing hydration water molecules in the minor groove of a TpA step

The hydration properties of the non-palindromic duplex d(CTACTGCTTTAG). d(CTAAAGCAGTAG) were investigated by NMR spectroscopy. The oligonucleotide possesses a heterogeneous B-DNA structure. The H2(n)-H1'(m+1) distances reflect a minor groove narrowing within the TTT/AAA segment (approximately 3.9A) and a sudden widening at the T10:A15 base-pair (approximately 5.3A), the standard B-DNA distance being approximately 5A. The facing T10pA11 and T14pA15 steps at the end of the TTTA/AAAT segment have completely different behaviors. Only A15 ending the AAA run displays NMR features comparable to those shown by adenines of TpA steps occupying the central position of TnAn (n> or =2) segments. These involve particular chemical shifts and line broadening of the H2 and H8 protons. Positive NOESY cross-peaks were measured between the water protons and the H2 protons of A15, A16 and A17 reflecting the occurrence of hydration water molecules with residence times longer than 500 picoseconds along the minor groove of the TTT/AAA segment. In contrast no water molecules with long residence times were observed neither for A3, A20 and A23 nor for A11 ending the 5'TTTA run. We confirm thus that the binding of water molecules with long residence time to adenine residues correlates with the minor groove narrowing. In contrast, the widening of the minor groove at the A11:T14 base-pair ending the TTTA/TAAA segment, likely associated to a high negative propeller twist value at this base-pair, prevents the binding of a water molecule with long residence time to A11 but not to A15 of the preceding T10:A15 base-pair. Thus, in our non-palindromic oligonucleotide the water molecules bind differently to A11 and A15 although both adenines are part of a TpA step. The slower motions occurring at A15 compared to A11 are also well explained by the present results.

The orientation and dynamics of the C2'-OH and hydration of RNA and DNA.RNA hybrids

The stereochemical and dynamic properties of the C2' hydroxyl group in several DNA.RNA hybrids have been measured by NMR and compared with the homologous RNA duplex. The C2'-OH NMR signals of the RNA strands were identified, and numerous specific assignments were made. The rate constants for exchange of the hydroxyl protons with water were determined at 5 degrees C, and were found to depend on both the position within a particular sequence and the nature of the duplex. On average, the exchange rate constants were slowest for the hybrids of composition rR.dY, and fastest for the RNA duplex, with an overall range of approximately 10-50/s. In the DNA.RNA hybrids, strong NOEs and ROEs were observed between the OH and the H1' of the same sugar, unambiguously showing that the OH proton points toward the H1' most of the time, and not toward the O3' of the same sugar. Evidence for significant hydration in both grooves of the DNA.RNA hybrids and the DNA duplex was found in ROESY and NOESY experiments. On average, the minor groove of the DNA.RNA hybrids showed more kinetically significant hydration than the DNA, which can be attributed to the hydrophilic lining of hydroxyl groups in RNA.

Conformational and dynamic properties of a 14 residue antifreeze glycopeptide from Antarctic cod

The 1H and 13C NMR spectra of a 14-residue antifreeze glycopeptide from Antarctic cod (Tetramatomnus borchgrevinki) containing two proline residues have been assigned. 13C NMR relaxation experiments indicate motional anisotropy of the peptide, with a tumbling time in water at 5 degrees C of 3-4 ns. The relaxation data and lack of long-range NOEs are consistent with a linear peptide undergoing significant segmental motion. However, extreme values of some coupling constants and strong sequential NOEs indicate regions of local order, which are most evident at the two ATPA subsequences. Similar spectroscopic properties were observed in the 16-residue analogue containing an Arg-Ala dipeptide added to the C-terminus. Molecular modeling also showed no evidence of long-range order, but the two ATPA subsequences were relatively well determined by the experimental data. These motifs were quite distinct from helical structures or beta turns commonly found in proteins, but rather resemble sections of an extended polyproline helix. Thus, the NMR data provide a description of the local order, which is of relevance to the mechanism of action of the antifreeze activity of the antifreeze glycopeptides as well as their ability to protect cells during hypothermic storage.

Comparison of the electrophoretic and hydrodynamic properties of DNA and RNA oligonucleotide duplexes

The electrophoretic behavior of defined DNA and RNA oligonucleotide duplexes from 10 to 20 bp in length has been investigated as a function of salt conditions, gel concentration, and temperature. The RNA oligomers migrated much more slowly than the DNA oligomers of the same sequence under all conditions. From sedimentation equilibrium and velocity measurements, the apparent partial specific volume in 0.1 M KCI, 20 mM NaPi, pH 7, was determined as 0.56 +/- 0.015 ml g(-1) for DNA and 0.508 ml g(-1) for RNA. The translational friction coefficients were determined and compared with the values calculated for cylinders. Taking into account the shape factors, the solution density, and partial specific volumes, the effective degree of hydration was estimated as 0.8-1 g g(-1) DNA. There was no significant difference in the frictional coefficients of the DNA and RNA oligomers, indicating that the effective sizes of DNA and RNA are very similar in solution. The differential electrophoretic mobility of DNA and RNA must arise from the differences in interaction with counterions, which is probably a global property of the oligonucleotides.

Specific binding of hoechst 33258 to the d(CGCAAATTTGCG)2 duplex: calorimetric and spectroscopic studies

Fluorescence spectroscopy and high-sensitivity isothermal titration calorimetry (ITC) techniques have been used to examine the binding characteristics of Hoechst 33258 with the extended AT-tract DNA duplex d(CGCAAATTTGCG)2 in aqueous solution. The method of continuous variation reveals a 1:1 binding stoichiometry. Fluorescence equilibrium studies carried out at three different, but fixed, ligand concentrations show that the binding isotherm shifts towards higher [DNA] as the concentration of ligand is increased. The data show tight binding with Kb=3.2(+/-0.6)x10(8) M(duplex)-1 at 25 degrees C in solutions containing 200 mM Na+. Based on UV studies of duplex melting, which show that strand separation starts at approximately 35 degrees C and has a Tm at 54 degrees C in 300 mM NaCl, binding enthalpies were determined by ITC in the 10 to 30 degrees C range. Binding is endothermic at all temperatures examined, with DeltaH values ranging from +10.24(+/-0.18) to +4.2(+/-0.10) kcal mol(duplex)-1 at 9.4 degrees C and 30.1 degrees C, indicating that the interaction is entropically driven. The temperature dependence of DeltaH shows a binding-induced change in heat capacity (DeltaCp) of -330(+/-50) cal mol-1 K-1. This value is similar to that predicted from a consideration of the effects of hydrophobic and hydrophilic solvent-accessible surface burial on complexation. This result, almost entirely dictated by a removal from exposure of the non-polar reactant surfaces, represents the first demonstration of such behavior in a DNA-drug system. The salt dependence of the binding constant was examined using reverse-salt fluorescence titrations, with a value of 0.99 determined for the deltalnK/deltaln[Na+] parameter. These data provide a detailed thermodynamic profile for the interaction that enables a dissection of DeltaGobs into the component free energy terms. Analysis of data obtained at 25 degrees C reveals that DeltaGobs is dominated by the free energy for hydrophobic transfer of ligand from solution to the DNA binding site. Molecular interactions, including H-bonding and van der Waals contacts, are found to play only a minor role in stabilizing the resulting complex, a somewhat surprising finding given the emphasis placed on such interactions from structural studies.

AT selectivity and DNA minor groove binding: modelling, NMR and structural studies of the interactions of propamidine and pentamidine with d(CGCGAATTCGCG)2

A molecular modelling strategy has been developed to identify potential binding sites for bis(amidine) ligands in the minor groove of duplex DNA. Calculations of interaction energy for propamidine and pentamidine with d(CGCGAAT TCGCG)2 show that this duplex contains two symmetrically equivalent binding sites of identical affinity, each displaced by 0.3-0.4 bp from the centre of the AT segment. The ligands occupy groove sites spanning approximately 4 and 4-5 bp, respectively with asymmetric binding to the 5'-AATT sequence. The DNA-bis(amidine) interactions have been examined by high-resolution 1H-NMR. The patterns of induced changes in DNA proton chemical shift and the DNA-ligand NOEs confirm that both agents bind in the AT minor groove in a non-centrosymmetric fashion. Detailed structures were determined for each complex using a NOE-restrained simulated annealing procedure, showing that the B-type DNA conformation is not significantly altered upon complexation with either ligand. The free DNA duplex has previously been shown to be extensively hydrated in the minor groove [Kubinec, M.G. and Wemmer, D.E. (1992) J. Am, Chem. Soc. 114, 8739-8740 Liepinsh, E. Otting, G. and Wüthrich, K. (1992) Nucleic Acids Res. 20. 6549-6553]. We detect hydration water close to the A(H2) protons in the presence of propamidine, which may stabilise certain waters against exchange. This conclusion supports recent crystallographic analyses, suggesting that such ligands may use water molecules to bridge between amidinium protons and host DNA bases Details of the ligand interactions with AT-tract DNA duplexes can now be compared for the subsequences 5'-AAT, 5'-AATT and 5'-AAATTT.