Structure of a B-DNA dodecamer. II. Influence of base sequence on helix structure

Detailed examination of the structure of the B-DNA dodecamer C-G-C-G-A-A-T-T-C-G-C-G, obtained by single-crystal X-ray analysis (Drew et al., 1981), reveals that the local helix parameters, twist, tilt and roll, are much more strongly influenced by base sequence than by crystal packing or any other external forces. The central EcoRI restriction endonuclease recognition site, G-A-A-T-T-C, is a B helix with an average of 9.8 base-pairs per turn. It is flanked on either side by single-base-pair steps having aspects of an A-like helix character. The dodecamer structure suggests several general principles, whose validity must be tested by other B-DNA analyses. (1) When an external bending moment is applied to a B-DNA double helix, it bends smoothly, without kinks or breaks, and with relatively little effect on local helix parameters. (2) Purine-3′,5′-pyrimidine steps open their base planes towards the major groove, pyrimidine-purine steps open toward the minor groove, and homopolymer (Pur-Pur, Pyr-Pyr) steps resist rolling in either direction. This behavior is related to the preference of pyrimidines for more negative glycosyl torsion angles. (3) CpG steps have smaller helical twist angles than do GpC, as though in compensation for their smaller intrinsic base overlap. Data on A-T steps are insufficient for generalization. (4) G.C base-pairs have smaller propellor twist than A · T, and this arises mainly from interstrand base overlap rather than the presence of the third hydrogen bond. (5) DNAase I cuts preferentially at positions of high helical twist, perhaps because of increased exposure of the backbone to attack. The correlation of the digestion patterns in solution and helical twist in the crystal argues for the essential identity of the helix structure in the two environments. (6) In the two places where the sequence TpCpG occurs, the C slips from under T in order to stack more efficiently over G. At the paired bases of this CpG step, the G and C are tilted so the angle between base planes is splayed out to the outside of the helix. This TpC is the most favored cutting site for DNAase I by a factor of 4.5 (Lomonossoff et al., 1981). (7) The EcoRI restriction endonuclease and methylase both appear to prefer a cutting site of the type purine-purine-A-T-T-pyrimidine, involving two adjacent homopolymer triplets, and this may be a consequence of the relative stiffness of homopolymer base-stacking observed in the dodecamer.     

Molecular dynamics of an in vacuo model of duplex d(CGCGAATTCGCG) in the B-form based on the amber 3.0 force field.

The characteristics of 100 ps of molecular dynamics (MD) on the DNA dodecamer d(CGCGAATTCGCG) at 300 K are described and investigated. The simulation is based on an in vacuo model of the oligomer and the AMBER 3.0 force field configured in the manner of Singh, U. C., S. J. Weiner, and P. A. Kollman, (1985, Proc. Natl. Acad. Sci. USA. 82:755-759). The analysis of the results was carried out using the "curves, dials, and windows" procedure (Ravishanker, G., S. Swaminathan, D. L. Beveridge, R. Lavery, and H. Sklenar. 1989. J. Biomol. Struct. Dyn. 6:669-699). The results indicate this dynamical model to be a provisionally stable double helix which lies at approximately 3.2 A rms deviation from the canonical B-form. There is, however, a persistent nonplanarity in the base pair orientations which resemble that observed in canonical A-DNA. The major groove width is seen to narrow during the course of the simulation and the minor groove expands, contravariant to the alterations in groove width seen in the crystal structure of the native dodecamer (Drew, H. R., R. M. Wing, T. Takano, C. Broka, S. Tanaka, I. Itakura, and R. E. Dickerson, 1981. Proc. Natl. Acad. Sci. USA. 78:2179-2183). The propeller twist in the bases, the sequence dependence of the base pair roll and aspects of bending in the helix axis are in some degree of agreement with the crystal structure. The patterns in DNA bending are observed to follow Zhurkin theory (Zhurkin, V. B. 1985. J. Biomol. Struct. Dyn. 2:785-804.). The relationship between the dynamical model and structure in solution is discussed.     

Crystal and molecular structure of a DNA fragment: d(CGTGAATTCACG)

The crystal structure of the dodecanucleotide d(CGTGAATTCACG) has been determined to a resolution of 2.7 A and refined to an R factor of 17.0% for 1532 reflections. The sequence crystallizes as a B-form double helix, with Watson-Crick base pairing. This sequence contains the EcoRI restriction endonuclease recognition site, GAATTC, and is flanked by CGT on the 5'-end and ACG on the 3'-end, in contrast to the CGC on the 5'-end and GCG on the 3'-end in the parent dodecamer d(CGCGAATTCGCG). A comparison with the isomorphous parent compound shows that any changes in the structure induced by the change in the sequence in the flanking region are highly localized. The global conformation of the duplex is conserved. The overall bend in the helix is 10 degrees. The average helical twist values for the present and the parent structures are 36.5 degrees and 36.4 degrees, respectively, corresponding to 10 base pairs per turn. The buckle at the substituted sites are significantly different from those seen at the corresponding positions in the parent dodecamer. Step 2 (GpT) is underwound with respect to the parent structure (27 degrees vs 36 degrees) and step 3 (TpG) is overwound (34 degrees vs 27 degrees). There is a spine of hydration in the narrow minor groove. The N3 atom of adenine on the substituted A10 and A22 bases are involved in the formation of hydrogen bonds with other duplexes or with water; the N3 atom of guanine on G10 and G22 bases in the parent structure does not form hydrogen bonds.     

Chiral phosphorothioate analogues of B-DNA. The crystal structure of Rp-d[Gp(S)CpGp(S)CpGp(S)C]

The compound Rp-d[Gp(S)CpGp(S)CpGp(S)C], an analogue of the deoxyoligomer d(G-C)3, crystallizes in space group P2(1)2(1)2(1) with a = 34.90 A, b = 39.15 A and c = 20.64 A. The structure, which is not isomorphous with any previously determined deoxyoligonucleotide, was refined to an R factor of 14.5% at a resolution of 2.17 A, with 72 solvent molecules located. The two strands of the asymmetric unit form a right-handed double helix, which is a new example of a B-DNA conformation and brings to light an important and overlooked component of flexibility of the double helix. This flexibility is manifest in the alternation of the backbone conformation between two states, defined by the adjacent torsion angles epsilon and zeta, trans . gauche-(BI) and gauche-. trans (BII). BI is characteristic of classical of B-DNA and has an average C(1') to C(1') separation of 4.5 A. The corresponding separation for BII is 5.3 A. Each state is associated with a distinct phosphate orientation where the plane of the PO2 (or POS) group is alternately near horizontal or vertical with respect to the helix axis. The BI and BII conformations are out of phase on the two strands. As a consequence, on one strand purine-pyrimidine stacking is better than pyrimidine-purine, while the converse holds for the other strand. At each base-pair step, good and bad stacking alternate across the helix axis. The pattern of alternation is regular in the context of a fundamental dinucleotide repeat. Re-examination of the B-DNA dodecamer d(C-G-C-G-A-A-T-T-C-G-C-G) shows that the C-G-C-G regions contain the BI and BII conformations, and the associated dual phosphate orientation and asymmetric base stacking. Different mechanisms are used in the two structures to avoid clashes between guanine residues on opposite strands, a combination of lateral slide, tilt and helical twist in the present structure, and base roll, tilt and longitudinal slide (Calladine rules) in the dodecamer. The flexibility of the phosphate orientations demonstrated in this structure is important, since it offers a structural basis for protein-nucleic acid recognition.     

Absence of minor groove monovalent cations in the crosslinked dodecamer C-G-C-G-A-A-T-T-C-G-C-G.

The structure of a crosslinked B -DNA dodecamer of sequence C-G-C-G-A-A-T-T-C-G-C-G has been solved to a resolution of 1.43 A. The dithiobis-propane crosslink, -CH2-CH2-CH2-S-S-CH2-CH2-CH2-, bridges N7 atoms of adenine bases 6 and 18 in the two central base-pairs within the major groove. The crosslink is sufficiently long that no bending is induced in the helix, which is essentially isostructural with the native unlinked dodecamer at 1.9 A. A constellation of solvent peaks tentatively fitted as a spermine molecule in that earlier analysis is now seen at higher resolution to be a well-defined octahedral magnesium hexahydrate complex in the major groove. One end of the duplex curves around that complex to produce a roll-bend near base-pairs 3-5, and an overall bend in helix axis, as has long been noted. Two other magnesium complexes connect the helices and help to knit the crystal lattice together. No evidence exists for partial sodium or potassium ion substitution for solvent water molecules within the minor groove spine of hydration, as had been suggested previously: not coordination geometry and environment, nor B values, nor calculated valence values, nor difference map analyses. Indeed, the very numbers that have been claimed in support of partial substitution by sodium or potassium ions are reproduced with the present crystals, which by chemical analysis contains only one trace sodium ion per 160 bp, and one potassium ion per 41 bp. In contrast, our crystals contain one Mg2+ per base-pair, meaning that phosphate group charge neutrality is accomplished by divalent cations, not monovalent ions. Three of these magnesium cations per duplex are localized and visible in the X-ray analysis, and nine are disordered and invisible. Hence although binding of monovalent cations within the minor groove of A -tracts on occasion may be a consequence of groove narrowing, it cannot be the cause of that narrowing. Cations, contrary to what has been claimed, are not in charge.     

The crystal structure analysis of d(CGCGAASSCGCG)2, a synthetic DNA dodecamer duplex containing four 4'-thio-2'-deoxythymidine nucleotides.

The crystal structure refinement of the synthetic dodecamer d(CGCGAASSCGCG), where S = 4'-thio-2'-deoxythymidine, has converged at R=0.201 for 2605 reflections with F > 2sigma(F) in the resolution range 8.0-2.4 A for a model consisting of the dodecamer duplex and 66 water molecules. A comparison of its structure with that of the native dodecamer d(CGCGAATTCGCG) has revealed that the major differences between the two structures is a change in the conformation of the sugar-phosphate backbone in the regions at and adjacent to the positions of the modified nucleosides. Examination of the fine structural parameters for each of the structures reveals that the thiosugars adopt a C3'-exo conformation in d(CGCGAASSCGCG), rather than the approximate C1'-exo conformation found for the analogous sugars in the structure of d(CGCGAATTCGCG). The observed differences in structure between the two duplexes may help to explain the enhanced resistance to nuclease digestion of synthetic oligonucleotides containing 4'-thio-2'-deoxynucleotides.     

Engineering DNA topology with locked nucleosides: a structural study

DNA dodecamers modified with nucleotide building blocks based on a bicyclo[3. 1.0]hexane system that effectively locks the ribose template into an RNA-like or North (N) conformation were analyzed by various biophysical techniques including high field nuclear magnetic resonance (NMR). Replacement of either one or both of the center thymidines in the Dickerson Drew dodecamer (CGCGAAT*T*CGCG) caused a progressive shift in the bending propensity of the double helix as shown by a newly developed rapid technique that compares the residual dipolar coupling (RDC) values of the modified duplexes with those previously determined for the native DNA.     

The effect of crystal packing on oligonucleotide double helix structure

One of the questions that constantly is asked regarding x-ray crystal structure analyses of macromolecules is: To what extent is the observed crystal structure representative of the molecular conformation when free in solution, and to what degree is the structure perturbed by intermolecular crystal forces? This can be assessed with DNA oligomers because of an unusual aspect of crystallization self-complementary oligomers should possess a twofold symmetry axis normal to their helix axis, yet more often than not crystal of such oligomers do not use this internal symmetry. The two ends of the helix are crystallographically distinct though chemically identical. Complexes of DNA oligomers with intercalating drugs such as triostin A tend to use their twofold symmetry when they crystallize, whereas complexes with non-intercalating, groove-binding drugs ignore this symmetry unless the drug molecule is very small. A detailed examination of crystal packing in the dodecamer C-G-C-G-A-A-T-T-C-G-C-G provides an explanation of all of the foregoing behavior in terms of the mechanism of nucleation of DNA or DNA-drug complexes on the surface of a growing crystal. Asymmetry of the ends of the DNA helix is the price that is paid for efficient lateral packing of helices within the crystal. The actual end-for-end variation in standard helix parameters is compared with the experimental noise level as gauged by independent re-refinement of the same oligonucleotide structure where available, and with the observed extent of variation of these same parameters along the helix. Oligomers analyzed are the B-DNA dodecamer C-G-C-G-A-A-T-T-C-G-C-G, the A-DNA octamer G-G-T-A-T-A-C-C, and the phosphorothioate analogue of the B-DNA hexamer G-C-G-C-G-C. End-for-end variation, presumably the result of crystal packing is typically double the experimental noise level, and half the variation in the same parameter along the helix. Analysis of crystal packing in the phosphorothioate hexamer, which uses the same P212121 space group as the dodecamer, shows that the highly unsymmetrical B1 vs. BII backbone conformation probably is to be ascribed to crystal packing forces, and not to the sequence of the hexamer.     

Molecular structure of the netropsin-d(CGCGATATCGCG) complex: DNA conformation in an alternating AT segment

The molecular structure of the complex between a minor groove binding drug (netropsin) and the DNA dodecamer d(CGCGATATCGCG) has been solved and refined by single-crystal X-ray diffraction analysis to a final R factor of 20.0% to 2.4-A resolution. The crystal is similar to that of the other related dodecamers with unit cell dimensions of a = 25.48 A, b = 41.26 A, and c = 66.88 A in the space group P2(1)2(1)2(1). In the complex, netropsin binds to the central ATAT tetranucleotide segment in the narrow minor groove of the dodecamer B-DNA double helix as expected. However, in the structural refinement the drug is found to fit the electron density in two orientations equally well, suggesting the disordered model. This agrees with the results from solution studies (chemical footprinting and NMR) of the interactions between minor groove binding drugs (e.g., netropsin and distamycin A) and DNA. The stabilizing forces between drug and DNA are provided by a combination of ionic, van der Waals, and hydrogen-bonding interactions. No bifurcated hydrogen bond is found between netropsin and DNA in this complex due to the unique dispositions of the hydrogen-bond acceptors (N3 of adenine and O2 of thymine) on the floor of the DNA minor groove. Two of the four AT base pairs in the ATAT stretch have low propeller twist angles, even though the DNA has a narrow minor groove. Alternating helical twist angles are observed in the ATAT stretch with lower twist in the ApT steps than in the TpA step.     

ENGINEERING DNA TOPOLOGY WITH LOCKED NUCLEOSIDES: A STRUCTURAL STUDY

DNA dodecamers modified with nucleotide building blocks based on a bicyclo[3.1.0]hexane system that effectively “locks” the ribose template into an RNA-like or “North” (N) conformation were analyzed by various biophysical techniques including high field nuclear magnetic resonance (NMR). Replacement of either one or both of the center thymidines in the Dickerson Drew dodecamer (CGCGAAT*T*CGCG) caused a progressive shift in the bending propensity of the double helix as shown by a newly developed rapid technique that compares the residual dipolar coupling (RDC) values of the modified duplexes with those previously determined for the native DNA.     

Mediation of the A/B-DNA helix transition by G-tracts in the crystal structure of duplex CATGGGCCCATG

The crystal structure of the DNA dodecamer duplex CATGGGCCCATG lies on a structural continuum along the transition between A- and B-DNA. The dodecamer possesses the normal vector plot and inclination values typical of B-DNA, but has the crystal packing, helical twist, groove width, sugar pucker, slide and x-displacement values typical of A-DNA. The structure shows highly ordered water structures, such as a double spine of water molecules against each side of the major groove, stabilizing the GC base pairs in an A-like conformation. The different hydration of GC and AT base pairs provides a physical basis for solvent-dependent facilitation of the A↔B helix transition by GC base pairs. Crystal structures of CATGGGCCCATG and other A/B-DNA intermediates support a ‘slide first, roll later’ mechanism for the B→A helix transition. In the distribution of helical parameters in protein–DNA crystal structures, GpG base steps show A-like properties, reflecting their innate predisposition for the A conformation.     

Crystal and molecular structure of the A-DNA dodecamer d(CCGTACGTACGG). Choice of fragment helical axis.

The crystal structure of the dodecamer d(CCGTACGTACGG) has been determined at 2.5 A resolution. The crystals grow in the hexagonal space group P6(1)22, a = b = 46.2 A, c = 71.5 A with one strand as the asymmetric unit. Diffraction data were collected by the oscillation film method yielding 1664 unique reflections with an Rmerge of 0.04. The structure was solved by real-space rotational translational searches with idealized helical models of A, B and Z-DNA. The best agreement was given by an A-DNA model with its dyad axis along the diagonal crystallographic dyad axis, with an R-factor 0.43 and correlation coefficient of 0.59 for data between 10 and 5 A. Iterative map fitting and restrained least-squares refinement and addition of 40 solvent molecules brought the R-factor to 0.15 and the correlation coefficient to 0.97 for all data between 8.0 and 2.5 A. The stereochemistry of the atomic model is good, with a root-mean-square deviation in bond distances of 0.006 A. This is the first example of an A-DNA containing a full helical turn. The dodecamer displays a novel packing motif. In addition to the characteristic contacts between the terminal base-pairs and the minor grooves of symmetry-related molecules, there are also minor groove to minor groove interactions not previously observed. The packing leaves an approximately 25 A diameter solvent channel around the origin, along the c-axis. The presence of a prominent 3.4 A meridional reflection and other diffuse features in the diffraction pattern provided evidence for the presence of disordered B-DNA along the c-axis, which can be accommodated in these solvent channels. The molecular conformation of the dodecamer also displays novel features. The dyad-related halves of the molecule are bent at an angle of 20 degrees, and the helical parameters are affected by this bend. Unlike the shorter A-DNA octamers, the dimensions of the major groove can be directly measured. Novel correlations between local helical parameters and global conformational features are presented. Most of the solvent molecules are associated with the major groove and the sugar-phosphate backbone.     

Experiment acknowledged the Watson-Crick hypothesis: a review of B-DNA double helix structural features at atomic resolution

The dodecamer d(CGCGAATTCGCG) forms a right-handed B-DNA double helix of a Watson-Crick type both in crystal and solution. It is the first piece of DNA longer than one helix turn whose molecular structure has become known at the atomic resolution. The article reviews qualitative aspects of its structure with a special emphasis on local variations in the disposition of base pairs in the double helix.     

Change in aggregation state of insulin upon conjugation with 5-dimethylaminonaphthalene-1-sulfonyl group

Change in aggregation state of insulin upon conjugation with 5-dimethylaminonaphthalene-1-sulfonyl (DNS) group was investigated at neutral pH. DNS group was introduced exclusively into B1 phenylalanine, the N-terminus of the B-chain of insulin. The association state of insulin shifted toward a more highly aggregated one upon conjugation, depending on the mole fraction (d) of DNS group to insulin monomer; at d equal 0.3 the equilibrium between dimer and hexamer was dominant over the range of 1-600 microM, while at d equal 1.0-1.5 DNS-insulin formed a larger aggregate (dodecamer) which is stable over the range of 67-600 microM. The dissociation constant of dimer-hexamer equilibrium at d=0.3 was evaluated to be 2.5 x 10(-10) M2 from the fluorescence anisotropy of the DNS group, which was about one order of magnitude smaller than that of the dimer-hexamer equilibrium in native insulin. Spectroscopic data and fluorescence decay analyses indicated that there exist at least two different environments surrounding the dye bound to B1 phenylalanine and that they are both relatively hydrophilic. It is considered that the major part of DNS group has excitation and emission maxima at longer wavelength with relatively low quantum yield, while the minor part has excitation and emission maxima at shorter wavelengths with relatively high quantum yield. The fluorescence lifetime of the dye was modified by the change in quaternary structure of DNS-lifetime of the dye was modified by the change in quaternary structure of DNS-insulin. Remarkable depolarization of DNS fluorescence was observed at d equal 1.0 and d equal 1.5 due to energy transfer between DNS groups conjugated to B1 phenylalanine in the hexamer or the dodecamer. Critical transfer distance for inter-DNS energy transfer was evaluated to be 15 A. From the molecular model of the insulin crystal, this energy transfer is ascribed to the close proximity, within about 15 A, between DNS groups in dimer units of the hexamer or the dodecamer.     

Crystal and solution structures of the B-DNA dodecamer d(CGCAAATTTGCG) probed by Raman spectroscopy: heterogeneity in the crystal structure does not persist in the solution structure

The self-complementary dodecamer d(CGCAAATTTGCG) crystallizes as a double helix of the B form and manifests a Raman spectrum with features not observed in Raman spectra of either DNA solutions or wet DNA fibers. A number of Raman bands are assigned to specific nucleoside sugar and phosphodiester conformations associated with this model B-DNA crystal structure. The Raman bands proposed as markers of the crystalline B-DNA structure are compared and contrasted with previously proposed markers of Z-DNA and A-DNA crystals. The results indicate that the three canonical forms of DNA can be readily distinguished by Raman spectroscopy. However, unlike Z-DNA and A-DNA, which retain their characteristic Raman fingerprints in aqueous solution, the B-DNA Raman spectrum is not completely conserved between crystal and solution states. The Raman spectra reveal greater heterogeneity of nucleoside conformations (sugar puckers) in the DNA molecules of the crystal structure than in those of the solution structure. The results are consistent with conversion of one-third of the dG residues from the C2'-endo/anti conformation in the solution structure to another conformation, deduced to be C1'-exo/anti, in the crystal. The dodecamer crystal also exhibits unusually broad Raman bands at 790 and 820 cm-1, associated with the geometry of the phosphodiester backbone and indicating a wider range of (alpha, zeta) backbone torsion angles in the crystal than in the solution structure. The results suggest that backbone torsion angles in the CGC and GCG sequences, which flank the central AAATTT sequence, are significantly different for crystal and solution structures, the former containing the greater diversity.(ABSTRACT TRUNCATED AT 250 WORDS)     

The solution structure of full‐length dodecameric MCM by SANS and molecular modeling

The solution structure of the full‐length DNA helicase minichromosome maintenance protein from Methanothermobacter thermautotrophicus was determined by small‐angle neutron scattering (SANS) data together with all‐atom molecular modeling. The data were fit best with a dodecamer (dimer of hexamers). The 12 monomers were linked together by the B/C domains, and the adenosine triphosphatase (AAA+) catalytic regions were found to be freely movable in the full‐length dodecamer both in the presence and absence of Mg²⁺ and 50‐meric single‐stranded DNA (ssDNA). In particular, the SANS data and molecular modeling indicate that all 12 AAA+ domains in the dodecamer lie approximately the same distance from the axis of the molecule, but the positions of the helix–turn–helix region at the C‐terminus of each monomer differ. In addition, the A domain at the N‐terminus of each monomer is tucked up next to the AAA+ domain for all 12 monomers of the dodecamer. Finally, binding of ssDNA does not lock the AAA+ domains in any specific position, which leaves them with the flexibility to move both for helicase function and for binding along the ssDNA. Proteins 2014; 82:2364–2374. © 2014 Wiley Periodicals, Inc.     

Secondary structure of the r(CUUCGG) tetraloop.

The oligonucleotide r(GGACUUCGGUCC) has been observed to adopt a hairpin conformation by solution NMR and a double helical conformation by X-ray diffraction. In order to understand this apparent conflict, we used time-resolved fluorescence depolarization and 19fluorine NMR to follow the secondary structure of this dodecamer as the solution composition was changed stepwise from the NMR experimental conditions to those used for crystallization. Calculation of the dodecamer concentration in the crystal (180 mM strands) and the cation concentration needed for neutrality (>2 M) prompted investigation of a tethered species, in which two dodecamers are connected by a string of 4 nt, geometrically equivalent to approximately 100 mM strands, in 2.5 M NaCl. The RNA tetraloop and its DNA analog maintain a single-strand hairpin conformation in solution, even under the conditions used to grow the crystal. Under high salt conditions, the tethered RNA and DNA analogs of this sequence yield secondary components which could be the double helical conformation. Crystal contacts in addition to solvent changes and high RNA concentrations are needed to obtain the double helix as the predominant species.     

Structure of d(CACGTG), a Z-DNA hexamer containing AT base pairs.

The left-handed Z-DNA conformation has been observed in crystals made from the self-complementary DNA hexamer d(CACGTG). This is the first time that a non disordered Z form is found in the crystal structure of an alternating sequence containing AT base pairs without methylated or brominated cytosines. The structure has been determined and refined to an agreement factor R = 22.9% using 746 reflections in the resolution in the resolution shell 7 to 2.5 A. The overall shape of the molecule is very similar to the Z-structure of the related hexamer d(CG)3 confirming the rigidity of the Z form. No solvent molecules were detected in the minor groove of the helix near the A bases. The disruption of the spine of hydration in the AT step appears to be a general fact in the Z form in contrast with the B form. The biological relevance of the structure in relation to the CA genome repeats is discussed.     

Binding of Hoechst 33258 to the minor groove of B-DNA

An X-ray crystallographic structure analysis has been carried out on the complex between the antibiotic and DNA fluorochrome Hoechst 33258 and a synthetic B-DNA dodecamer of sequence C-G-C-G-A-A-T-T-C-G-C-G. The drug molecule, which can be schematized as: phenol-benzimidazole-benzimidazole-piperazine, sits within the minor groove in the A-T-T-C region of the DNA double helix, displacing the spine of hydration that is found in drug-free DNA. The NH groups of the benzimidazoles make bridging three-center hydrogen bonds between adenine N-3 and thymine O-2 atoms on the edges of base-pairs, in a manner both mimicking the spine of hydration and calling to mind the binding of the auti-tumor drug netropsin. Two conformers of Hoechst are seen in roughly equal populations, related by 180 degrees rotation about the central benzimidazole-benzimidazole bond: one form in which the piperazine ring extends out from the surface of the double helix, and another in which it is buried deep within the minor groove. Steric clash between the drug and DNA dictates that the phenol-benzimidazole-benzimidazole portion of Hoechst 33258 binds only to A.T regions of DNA, whereas the piperazine ring demands the wider groove characteristic of G.C regions. Hence, the piperazine ring suggests a possible G.C-reading element for synthetic DNA sequence-reading drug analogs.