Molecular dynamics simulations of xDNA

xDNA is a modified DNA, which contains natural as well as expanded bases. Expanded bases are generated by the addition of a benzene spacer to the natural bases. A set of AMBER force‐field parameters were derived for the expanded bases and the structural dynamics of the xDNA decamer ( xT5 ′ G xT A xC xG C xA xG T3′ ) · ( xA5′ C T xG C G xT A xC A3′) was explored using a 22 ns molecular dynamics simulation in explicit solvent. During the simulation, the duplex retained its Watson‐Crick base‐pairing and double helical structure, with deviations from the starting B‐form geometry towards A‐form; the deviations are mainly in the backbone torsion angles and in the helical parameters. The sugar pucker of the residues were distributed among a variety of modes; C2′ endo, C1′ exo, O4′ endo, C4′ exo, C2′ exo, and C3′ endo. The enhanced stacking interactions on account of the modification in the bases could help to retain the duplex nature of the helix with minor deviations from the ideal geometry. In our simulation, the xDNA showed a reduced minor groove width and an enlarged major groove width in comparison with the NMR structure. Both the grooves are larger than that of standard B‐DNA, but major groove width is larger than that of A‐DNA with almost equal minor groove width. The enlarged groove widths and the possibility of additional hydration in the grooves makes xDNA a potential molecule for various applications. © 2009 Wiley Periodicals, Inc. Biopolymers 91: 351–360, 2009. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

Molecular dynamics simulations of excimer forming (+)-anti-BPDE-DNA adducts in aqueous solution

The chemical carcinogen (+)-anti BPDE preferentially binds covalently to the guanine base in the minor groove of DNA. Fluorescence spectroscopic studies have shown that the BPDE molecules bound to DNA can interact in their photo-excited state giving strong excimer fluorescence when bound to poly(dGdC) · poly(dGdC). It was suggested that the formation of such excited state complexes is most probable when the two (+)-anti-BPDE bind to guanines of adjacent base pairs on the two different strands of the DNA. In the present work a model for such an excimer forming DNA-BPDE double adduct system has been constructed and shown to be stable over a 300 ps molecular dynamics simulation in a water box. The model is a d(CG)₃ · d(CG)₃ molecule with two BPDE molecules bound to the guanines at the 4th position on each strand, located in the minor groove and each oriented towards the 5 end of the modified strand, respectively. The results of 300 ps MD simulation show that the two BPDE chromophores exhibited on the average a relative geometry favourable for excimer formation. The local structure at the adduct position was considerably distorted and the helix axis was bent. The modified bases were found to be paired through a stable single non-Watson Crick type of hydrogen bond. Correspondence to: A. Gräslund     

Influence of divalent magnesium ion on DNA: molecular dynamics simulation studies

A large amount of experimental evidence is available on the effect of magnesium ions on the structure and stability of DNA double helix. Less is known, however, on how these ions affect the stability and dynamics of the molecule. The static time average pictures from X-ray structures or the quantum chemical energy minimized structures lack understanding of the dynamic DNA–ion interaction. The present work addresses these questions by molecular dynamics simulation studies on two DNA duplexes and their interaction with magnesium ions. Results show typical B-DNA character with occasional excursions to deviated states. We detected expected stability of the duplexes in terms of backbone conformations and base pair parameter by the CHARMM-27 force field. Ion environment analysis shows that Mg²⁺ retains the coordination sphere throughout the simulation with a preference for major groove over minor. An extensive analysis of the influence of the Mg²⁺ ion shows no evidence of the popular predictions of groove width narrowing by dipositive metal ion. The major groove atoms show higher occupancy and residence time compared to minor groove for magnesium, where no such distinction is found for the charge neutralizing Na⁺ ions. The determining factor of Mg²⁺ ion’s choice in DNA binding site evolves as the steric hindrance faced by the bulky hexahydrated cation where wider major groove gets the preference. We have shown that in case of binding of Mg²⁺ to DNA non electrostatic contributions play a major role.

An animated Interactive 3D Complement (I3DC) is available in Proteopedia at http://proteopedia.org/w/Journal:JBSD:5     

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

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

Mode‐coupling Smoluchowski dynamics of a double‐stranded DNA oligomer

The local dynamics of a double‐stranded DNA d(TpCpGpCpG)₂ is obtained to second order in the mode‐coupling expansion of the Smoluchowski diffusion theory. The time correlation functions of bond variables are derived and the ¹³C‐nmr spin–lattice relaxation times T₁ of different ¹³C along the chains are calculated and compared to experimental data from the literature at three frequencies. The DNA is considered as a fluctuating three‐dimensional structure undergoing rotational diffusion. The fluctuations are evaluated using molecular dynamics simulations, with the ensemble averages approximated by time averages along a trajectory of length 1 ns. Any technique for sampling the configurational space can be used as an alternative. For a fluctuating three‐dimensional (3D) structure using the three first‐order vector modes of lower rates, higher order basis sets of second‐rank tensor are built to give the required mode coupling dynamics. Second‐ and even first‐order theories are found to be in close agreement with the experimental results, especially at high frequency, where the differences in T₁ for ¹³C in the base pairs, sugar, and backbone are well described. These atomistic calculations are of general application for studying, on a molecular basis, the local dynamics of fluctuating 3D structures such as double‐helix DNA fragments, proteins, and protein–DNA complexes. © 1999 John Wiley & Sons, Inc. Biopoly 50: 613–629, 1999     

Probing the interactions of the solvated electron with DNA by molecular dynamics simulations: II. bromodeoxyuridine-thymidine mismatched DNA

The interaction of solvated electrons with DNA results in various types of DNA lesions. The in vitro and in vivo sensitisation of DNA to -induced damage is achieved by incorporation of the electron-affinity radiosensitiser bromodeoxyuridine (BUdR) in place of thymidine. However, in DNA duplexes containing single-stranded regions (bulged BUdR-DNA), the type of lesion is different and the efficiency of damage is enhanced. In particular, DNA interstrand crosslinks (ICL) form at high efficiency in bulged DNA but are not detectable in completely duplex DNA. Knowledge about the processes and interactions leading to these differences is obscure. Previously, we addressed the problem by applying molecular modelling and molecular dynamics (MD) simulations to a system of normal (BUdR·A)-DNA and a hydrated electron, where the excess electron was modelled as a localised eˉ(H₂O)₆ anionic cluster. The goal of the present study was to apply the same MD simulation to a wobble system, containing a pyrimidine–pyrimidine mismatched base pair, BUdR·T. The results show an overall dynamic pattern similar to that of the motion around normal DNA. However, the number of configuration states when was particularly close to DNA is different. Moreover, in the (BUdR·T)-wobble DNA system, the electron frequently approaches the brominated strand, including BUdR, which was not observed with the normal (BUdR·A)-DNA. The structure and exchange of water at the sites of immobilisation near DNA were also characterised. The structural dynamics of the wobble DNA is prone to more extensive perturbations, including frequent formation of cross-strand (cs) interatomic contacts. The structural deviations correlated with approaching DNA from the major groove side, with sodium ions trapped deep in the minor groove. Altogether, the obtained results confirm and/or throw light on dynamic-structure determinants possibly responsible for the enhanced radiation damage of wobble DNA. Figure The structure of the tightly bound single water-layer between the DNA and the electron (Site-8, five H₂O molecules, bold capped sticks); the rest of the “second” shell waters (lines, in atom type colour) surround the ˉ(H₂O)₆ cluster (yellow, space fill). Orange dashed lines H-bonds; only one of the five molecules from the single H₂O layer mediates a single-step H-bond bridge with N7(A8); the other four present a network of two(three)-step H-bond bridges between DNA/ partner atoms     

An electrostatic model for the dielectric effects,the adsorption of multivalent ions,and the bending of B-DNA

We have studied electrostatic properties of DNA with a discrete charge model consisting of a cylindrical dielectric core with a radius of 8 Å and a dielectric constant D_i = 4, surrounded by two helical strings of phosphate point charges at 10 Å from the axis, immersed in an aqueous medium with dielectric constant D_w = 78.54. Eliminating the dielectric core makes potentials in the phosphate surface less negative by about 0.5 kT/e. Salt effects are evaluated for the model without a dielectric core, using the shielded Coulomb potential. Smearing the phosphate charges increases their potential by about 2.5 kT/e, due mostly to the self-potential of the smeared charge. Potentials in the center of the minor and major grooves vary less than 0.02 kT/e along their helical path. The potential in the center of the minor groove is from 1.0 to 1.7 kT/e, more negative than in the center of the major groove, depending on dielectric core and salt concentration. So multivalent cations and also larger cationic ligands, such as some antibiotics, are likely to adsorb in the minor groove, in agreement with earlier computations by A. and B. Pullman. Dielectric effects on the surface potential and the local potential variations are found to be relatively small. Bending of DNA is studied by placing a multivalent cation, M^Z+, in the center of the minor or major groove, curving DNA around it for a certain length, and calculating the free energy difference between the bent and the straight configuration. Boltzmann averaged bending angles, 〈β〉, are found to be maximal in 0.03M monovalent salt, for a length of about 50 or 25 Å of curved DNA when an M^Z+ ion is adsorbed in the minor or the major groove, respectively. When the dielectric constant of water is used throughout the calculation, we find maximal bends of 〈β〉 = 11° for M²⁺ and 〈β〉 = 16° for M³⁺ in the minor groove, 〈β〉 = 13° for M³⁺ in the major groove. The absence of bends in DNA adsorbed to mica in the presence of Mg salts supports the role of Mg²⁺ in “ion bridging” between DNA and mica. The treatment of the effective dielectric constant between two points outside a dielectric cylinder in water is appended. © 1998 John Wiley & Sons, Inc. Biopoly 46: 503–516, 1998     

Population dynamics of a thistle-feeding lady beetle Epilachna niponica (Coccinellidae: Epilachninae) in Kanazawa, Japan. 1. Adult demographic traits and population stability

 Population dynamics of a thistle-feeding univoltine lady beetle, Epilachna niponica Lewis (Coleoptera: Coccinellidae), was studied from 1996 to 1999 in Yuwaku, Kanazawa, Japan. The lady beetles often reached such a high density level that food was depleted. The Jolly–Seber method was used for adult marking, release, and recapture data to estimate population parameters of adult number, daily resident rate, longevity, reproductive rate (R, the number of new adults produced per overwintered adults), and survival rate of new adults to the reproductive seasons (S _w). These estimates were compared with those of the Asiu, Kutsuki (A and F), and Kyoto populations, which were previously studied with similar methods and have similar intensities. Asiu and Kutsuki F populations remained at a rather low density with a low R, while Kutsuki A and Kyoto populations reached a high density where food depletion occurred with a high R value. The Yuwaku population often reached a food-depleting level as in the Kutsuki A and Kyoto populations. It also shared the short life span of overwintered adults (13.5 days) of other high-density populations; however, it showed much shorter longevity of new adults (36.6 days), much lower R (1.0–2.5), and higher S _w (43%–53%). In some traits the Yuwaku population was similar to the Asiu population: low R, high S _w, and low population variability (SD of log densities; 0.103 and 0.115 for overwintered and new adults, respectively, which were lowest among the populations). Received: July 26, 2001 / Accepted: May 21, 2002     

DNA and its counterions: a molecular dynamics study

The behaviour of mobile counterions, Na⁺ and K⁺, was analysed around a B-DNA double helix with the sequence CCATGCGCTGAC in aqueous solution during two 50 ns long molecular dynamics trajectories. The movement of both monovalent ions remains diffusive in the presence of DNA. Ions sample the complete space available during the simulation time, although individual ions sample only about one-third of the simulation box. Ions preferentially sample electronegative sites around DNA, but direct binding to DNA bases remains a rather rare event, with highest site occupancy values of <13%. The location of direct binding sites depends greatly on the nature of the counterion. While Na⁺ binding in both grooves is strongly sequence-dependent with the preferred binding site in the minor groove, K⁺ mainly visits the major groove and binds close to the centre of the oligomer. The electrostatic potential of an average DNA structure therefore cannot account for the ability of a site to bind a given cation; other factors must also play a role. An extensive analysis of the influence of counterions on DNA conformation showed no evidence of minor groove narrowing upon ion binding. A significant difference between the conformations of the double helix in the different simulations can be attributed to extensive α/γ transitions in the phosphate backbone during the simulation with Na⁺. These transitions, with lifetimes over tens of nanoseconds, however, appear to be correlated with ion binding to phosphates. The ion-specific conformational properties of DNA, hitherto largely overlooked, may play an important role in DNA recognition and binding.     

Spectroscopic studies of 9-hydroxyellipticine binding to DNA

The binding of 9-hydroxyellipticine to calf thymus DNA, poly[d(A-T)]₂, and poly-[d(G-C)]₂ has been studied in detail by means of CD, linear dichroism, resonance light scattering, and molecular dynamics. The transition moment polarizations of 9-hydroxyelliptiycine were determined in polyvinyl alcohol stretched film. Spectroscopic solution studies of the DNA/drug complex are combined with theoretical CD calculations using the final 50 ps of a series of molecular dynamics simulations as input. The spectroscopic data shows 9-hydroxyellipticine to adopt two main binding modes, one intercalative and the other a stacked binding mode involving the formation of drug oligomers in the DNA major groove. Analysis of the intercalated binding mode in poly[d(A-T)]₂ suggests the 9-hydroxyellipticine hydroxyl group lies in the minor groove and hydrogen bonds to water with the pyridine ring protruding into the major groove. The stacked binding mode was examined using resonance light scattering and it was concluded that the drug was forming small oligomer stacks rather than extended aggregates. Reduced linear dichroism measurements suggested a binding geometry that precluded a minor groove binding mode where the plane of the drug makes a 45° angle with the plane of the bases. Thus it was concluded that the drug stacks in the major groove. No obvious differences in the mode of binding of 9-hydroxyellipticine were observed between different DNA sequences; however, the stacked binding mode appeared to be more favorable for calf thymus DNA and poly[d(G-C)]₂ than for poly[d(A-T)]₂, an observation that could be explained by the slightly greater steric hindrance of the poly[d(A-T)]₂ major groove. A strong concentration dependence was observed for the two binding modes where intercalation is favored at very low drug load, with stacking interactions becoming more prominent as the drug concentration is increased. Even at DNA : drug mixing ratios of 70:1 the stacked binding mode was still important for GC-rich DNAs. © 1998 John Wiley & Sons, Inc. Biopoly 46: 127–143, 1998     

DNA groove-binding and acid-base properties of a Ru(II) complex containing anthryl moieties

Abstract

DNA groove binders have been poorly studied as compared to the intercalators. A novel Ru(II) complex of [Ru(aeip)₂(Haip)](PF₆)₂ {Haip?=?2-(9-anthryl)-1H-imidazo[4,5-f][1,10]phenanthroline and aeip = 2-(anthracen-9-yl)-1-ethyl-imidazo[4,5-f][1, 10]phenanthroline} is synthesized and characterized by elemental analysis, ¹H NMR spectroscopy and mass spectrometry. The complex is evidenced to be a calf-thymus DNA groove binder with a large intrinsic binding constant of 10⁶ M^?1 order of magnitude as supported by UV–visible absorption spectral titrations, salt effects, DNA competitive binding with ethidium bromide, DNA melting experiment, DNA viscosity measurements and density functional theory calculations. The acid-base properties of the complex studied by UV–Vis spectrophotometric titrations are reported as well.     

Molecular dynamics simulations and binding free energy analysis of DNA minor groove complexes of curcumin

Curcumin is a natural phytochemical that exhibits a wide range of pharmacological properties, including antitumor and anticancer activities. The similarity in the shape of curcumin to DNA minor groove binding drugs is the motivation for exploring its binding affinity in the minor grooves of DNA sequences. Interactions of curcumin with DNA have not been extensively examined, while its pharmacological activities have been studied and documented in depth. Curcumin was docked with two DNA duplexes, d(GTATATAC)₂ and d(CGCGATATCGCG)₂, and molecular dynamics simulations of the complexes were performed in explicit solvent to determine the stability of the binding. In all systems, the curcumin is positioned in the minor groove in the A·T region, and was stably bound throughout the simulation, causing only minor modifications to the structural parameters of DNA. Water molecules were found to contribute to the stability of the binding of the ligand. Free energy analyses of the complexes were performed with MM-PBSA, and the binding affinities that were calculated are comparable to the values reported for other similar nucleic acid–ligand systems, indicating that curcumin is a suitable natural molecule for the development of minor groove binding drugs.     

Peculiar features of kink dynamics in inhomogeneous DNA

A mathematical model was constructed for studies of the peculiar features of the dynamics of local conformational distortions (kinks) in inhomogeneous DNA. General formulas were obtained that determine the dependence of the main dynamic characteristics of kinks: size, energy, energy density and velocity, on DNA composition. Quantitative estimates were made for the characteristics of kinks activated in polynucleotide chains analogous to promoter sequences A₁, A₂, and A₃ of bacteriophage T7.     

DrRecQ regulates guanine quadruplex DNA structure dynamics and its impact on radioresistance in Deinococcus radiodurans

The GC‐rich genome of Deinococcus radiodurans contains a very high density of putative guanine quadruplex (G4) DNA motifs and its RecQ (drRecQ) was earlier characterized as a 3′→5′ dsDNA helicase. We saw that N‐Methyl mesoporphyrin IX (NMM), a G4 DNA binding drug affected normal growth as well as the gamma radiation resistance of the wild‐type bacterium. Interestingly, NMM treatment and recQ deletion showed additive effect on normal growth but there was no effect of NMM on gamma radiation resistance of recQ mutant. The recombinant drRecQ showed ~400 times higher affinity to G4 DNA (K_d = 11.74 ± 1.77 nM) as compared to dsDNA (K_d = 4.88 ± 1.30 µM). drRecQ showed ATP independent helicase function on G4 DNA, which was higher than ATP‐dependent helicase activity on dsDNA. Unlike wild‐type cells that sparingly stained for G4 structure with Thioflavin T (ThT), recQ mutant showed very high‐density of ThT fluorescence foci on DNA indicating an important role of drRecQ in regulation of G4 DNA structure dynamics in vivo. These results together suggested that drRecQ is an ATP independent G4 DNA helicase that plays an important role in the regulation of G4 DNA structure dynamics and its impact on radioresistance in D. radiodurans.     

1H NMR studies: dynamics of water in gelatin

 Proton magnetic resonance was used to characterize the dynamics of water in gelatin. Both sol and gel states were investigated. Transverse relaxation rates (R ₂) were dependent on the proton frequency measurement. (R ₂) measured with the Carr-Purcell-Meiboom-Gill pulse sequence was dependent on pulse spacing. These observations were interpreted in terms of chemical exchanges between water protons and those of the macromolecules in the sol state, whereas in the gel state the contribution of diffusion through microheterogeneities in the sample seems to provide an additional transverse relaxation mechanism. Received: 10 May 1999 / Revised version: 13 December 1999 / Accepted: 25 January 2000     

Influence of Spermine on DNA Conformation in a Molecular Dynamics Trajectory of d(CGCGAATTCGCG)2: Major Groove Binding by One Spermine Molecule Delays the A→B Transition

Abstract

The effect of spermine on the A-DNA to B-DNA transition in d(CGCGAATTCGCG)₂ has been investigated by five A-start molecular dynamics simulations, using the Cornell et al. potential. In the absence of spermine an A→B transition is initiated immediately and the DNA becomes equidistant from the A- and B-forms at 200ps. In three DNA-spermine simulations, when a spermine is located across the major groove of A-DNA in one of three different initial locations, the time taken to reach equidistance from the A- and B-forms is delayed until 800, 950 or 1000ps. In each case the A-form appears to be temporarily stabilized by spermine's electrostatic interactions with phosphates on both sides of the major groove. The onset of the A→B transition can be correlated with the spermine losing contact with phosphates on one side of the groove and with A-like → B-like sugar pucker transitions in the vicinity of the spermine bridge. However in the fifth trajectory, in which the spermine initially threads from the major groove via the backbone into the minor groove, the B→A transition occurs rapidly once again and the DNA is equidistant between the A- and B-forms within 300ps. This indicates that the mere presence of spermine is insufficient to delay the transition and that major groove binding stabilizes A-DNA.     

Role of 6-Mercaptopurine in the potential therapeutic targets DNA base pairs and G-quadruplex DNA: insights from quantum chemical and molecular dynamics simulations

The theoretical studies on DNA with the anticancer drug 6-Mercaptopurine (6-MP) are investigated using theoretical methods to shed light on drug designing. Among the DNA base pairs considered, 6-MP is stacked with GC with the highest interaction energy of –46.19 kcal/mol. Structural parameters revealed that structure of the DNA base pairs is deviated from the planarity of the equilibrium position due to the formation of hydrogen bonds and stacking interactions with 6-MP. These deviations are verified through the systematic comparison between X–H bond contraction and elongation and the associated blue shift and red shift values by both NBO analysis and vibrational analysis. Bent’s rule is verified for the C–H bond contraction in the 6-MP interacted base pairs. The AIM results disclose that the higher values of electron density (ρ) and Laplacian of electron density (?²ρ) indicate the increased overlap between the orbitals that represent the strong interaction and positive values of the total electron density show the closed-shell interaction. The relative sensitivity of the chemical shift values for the DNA base pairs with 6-MP is investigated to confirm the hydrogen bond strength. Molecular dynamics simulation studies of G-quadruplex DNA d(TGGGGT)₄ with 6-MP revealed that the incorporation of 6-MP appears to cause local distortions and destabilize the G-quadruplex DNA.     

Induced fit DNA recognition by a minor groove binding analogue of Hoechst 33258: fluctuations in DNA A tract structure investigated by NMR and molecular dynamics simulations

NMR analysis and molecular dynamics simulations of d(GGTAATTACC)₂ and its complex with a tetrahydropyrimidinium analogue of Hoechst 33258 suggest that DNA minor groove recognition in solution involves a combination of conformational selection and induced fit, rather than binding to a preorganised site. Analysis of structural fluctuations in the bound and unbound states suggests that the degree of induced fit observed is primarily a consequence of optimising van der Waals contacts with the walls of the minor groove resulting in groove narrowing through: (i) changes in base step parameters, including increased helical twist and propeller twist; (ii) changes to the sugar–phosphate backbone conformation to engulf the bound ligand; (iii) suppression of bending modes at the TpA steps. In contrast, the geometrical arrangement of hydrogen bond acceptors on the groove floor appears to be relatively insensitive to DNA conformation (helical twist and propeller twist). We suggest that effective recognition of DNA sequences (in this case an A tract structure) appears to depend to a significant extent on the sequence being flexible enough to be able to adopt the geometrically optimal conformation compatible with the various binding interactions, rather than involving ‘lock and key’ recognition.     

Enhanced uranyl photocleavage across the minor groove of all (A/T)4 sequences indicates a similar narrow minor groove conformation

The uranyl(VI)-mediated photocleavage of a Drew–Dickerson sequence oligonucleotide (5′-dGATCACGCGAATTCGCGT) either as the (self-complementary) duplex or cloned into the BamH1 site of pUC19 has been studied. At pH 6.5 in acetate buffer relatively enhanced photocleavage is observed at the 3′-end of the AATT sequence corresponding to maximum cleavage across the minor groove in the A/T tract. Thus maximum cleavage correlates with minimum minor groove width in the crystal structure and also with the largest electronegative potential according to computations. Using plasmid constructs with cloned inserts of the type [CGCG(A/T₄)]_n, we also analysed all possible sequence combinations of the (A/T)₄ tract and in all cases we observed maximum uranyl-mediated photocleavage across the minor groove in the (A/T)₄ tract without any significant differences between the various sequences. From these results we infer that DNA double helices of all (A/T)₄ sequences share the same narrow minor groove helix conformation.