A unified model for the origin of DNA sequence-directed curvature

The fine structure of the DNA double helix and a number of its physical properties depend upon nucleotide sequence. This includes minor groove width, the propensity to undergo the B-form to A-form transition, sequence-directed curvature, and cation localization. Despite the multitude of studies conducted on DNA, it is still difficult to appreciate how these fundamental properties are linked to each other at the level of nucleotide sequence. We demonstrate that several sequence-dependent properties of DNA can be attributed, at least in part, to the sequence-specific localization of cations in the major and minor grooves. We also show that effects of cation localization on DNA structure are easier to understand if we divide all DNA sequences into three principal groups: A-tracts, G-tracts, and generic DNA. The A-tract group of sequences has a peculiar helical structure (i.e., B*-form) with an unusually narrow minor groove and high base-pair propeller twist. Both experimental and theoretical studies have provided evidence that the B*-form helical structure of A-tracts requires cations to be localized in the minor groove. G-tracts, on the other hand, have a propensity to undergo the B-form to A-form transition with increasing ionic strength. This property of G-tracts is directly connected to the observation that cations are preferentially localized in the major groove of G-tract sequences. Generic DNA, which represents the vast majority of DNA sequences, has a more balanced occupation of the major and minor grooves by cations than A-tracts or G-tracts and is thereby stabilized in the canonical B-form helix. Thus, DNA secondary structure can be viewed as a tug of war between the major and minor grooves for cations, with A-tracts and G-tracts each having one groove that dominates the other for cation localization. Finally, the sequence-directed curvature caused by A-tracts and G-tracts can, in both cases, be explained by the cation-dependent mismatch of A-tract and G-tract helical structures with the canonical B-form helix of generic DNA (i.e., a cation-dependent junction model).     

Uranyl photoprobing of nonbent A/T- and bent A-tracts. A difference of flexibility?

In this study, we have systematically compared the uranyl photocleavage of a range of bent A-tracts and nonbent TA-tracts as well as interrupted A-tracts. We demonstrate that uranyl photocleavage of A-tracts and TA-tracts is almost identical, indicating a very similar minor groove conformation. Furthermore, a 10 base pair A-tract is divided into two independent tracts by an intervening TA or GC step. Uranyl probing also clearly distinguishes the bent A4T4 and the nonbent T4A4 sequences as adopting different structures, and our interpretation of the data is consistent with a structure for the bent A4T4 sequence that resembles a continuous A-tract, whereas the nonbent T4A4 sequences are closer to two independent and opposite A-tracts that cancel each other in terms of macroscopic bending. Finally, we also note that even single TA and TAT steps are highly sensitive to uranyl photocleavage and propose that in addition to average minor groove width, uranyl also senses DNA helix flexibility/deformability. Thus, the structural difference of TA-tracts and A-tracts may to a large extent reflect a difference in flexibility, and DNA curvature may consequently require a rigid narrow minor groove conformation that creates distinct A-tract-B-DNA junctions as the predominant cause of the bending.     

Intrinsic bending and deformability at the T-A step of CCTTTAAAGG: a comparative analysis of T-A and A-T steps within A-tracts

Introduction of a T-A or pyrimidine-purine step into a straight and rigid A-tract can cause a positive roll deformation that kinks the DNA helix at that step. In CCTTTAAAGG, the central T-A step has an 8.6 degrees bend toward the major groove. We report the structural analysis of CCTTTAAAGG and a comparison with 25 other representative crystal structures from the NDB containing at least four consecutive A or T bases. On average, more local bending occurs at the disruptive T-A step (8.21 degrees ) than at an A-T step (5.71 degrees ). In addition, A-tracts containing an A-T step are more bent than are pure A-tracts, and hence A-A and A-T steps are not equivalent. All T-A steps examined exhibit positive roll, bending towards the major groove, while A-T steps display negative roll and bend slightly towards the minor groove. This illustrates how inherent negative and positive roll are, respectively, at A-T and T-A steps within A-tracts. T-A steps are more deformable, showing larger and more variable deformations of minor groove width, rise, cup, twist, and buckle. Standard deviations of twist, rise, and cup for T-A steps are 6.66 degrees, 0.55 A, and 15.90 degrees, versus 2.28 degrees, 0.21 A, and 2.99 degrees for A-T steps. Packing constraints determine which local values of these helical parameters an individual T-A step will adopt. For instance, with CCTTTAAAGG and three isomorphous structures, CGATTAATCG, CGATATATCG, and CGATCGATCG, crystal packing forces lead to a series of correlated changes: widened minor groove, large slide, low twist, and large rise. The difference in helical parameters between A-T steps lying within A-tracts, versus A-T steps within alternating AT sequences, demonstrates the importance of neighboring steps on the conformation of a given dinucleotide step.     

Molecular dynamics simulations of DNA curvature and flexibility: helix phasing and premelting

Recent studies of DNA axis curvature and flexibility based on molecular dynamics (MD) simulations on DNA are reviewed. The MD simulations are on DNA sequences up to 25 base pairs in length, including explicit consideration of counterions and waters in the computational model. MD studies are described for ApA steps, A-tracts, for sequences of A-tracts with helix phasing. In MD modeling, ApA steps and A-tracts in aqueous solution are essentially straight, relatively rigid, and exhibit the characteristic features associated with the B'-form of DNA. The results of MD modeling of A-tract oligonucleotides are validated by close accord with corresponding crystal structure results and nuclear magnetic resonance (NMR) nuclear Overhauser effect (NOE) and residual dipolar coupling (RDC) structures of d(CGCGAATTCGCG) and d(GGCAAAAAACGG). MD simulation successfully accounts for enhanced axis curvature in a set of three sequences with phased A-tracts studied to date. The primary origin of the axis curvature in the MD model is found at those pyrimidine/purine YpR "flexible hinge points" in a high roll, open hinge conformational substate. In the MD model of axis curvature in a DNA sequence with both phased A-tracts and YpR steps, the A-tracts appear to act as positioning elements that make the helix phasing more precise, and key YpR steps in the open hinge state serve as curvature elements. Our simulations on a phased A-tract sequence as a function of temperature show that the MD simulations exhibit a premelting transition in close accord with experiment, and predict that the mechanism involves a B'-to-B transition within A-tracts coupled with the prediction of a transition in key YpR steps from the high roll, open hinge, to a low roll, closed hinge substate. Diverse experimental observations on DNA curvature phenomena are examined in light of the MD model with no serious discrepancies. The collected MD results provide independent support for the "non-A-tract model" of DNA curvature. The "junction model" is indicated to be a special case of the non-A-tract model when there is a Y base at the 5' end of an A-tract. In accord with crystallography, the "ApA wedge model" is not supported by MD.     

Atomic-resolution crystal structures of B-DNA reveal specific influences of divalent metal ions on conformation and packing.

Crystal structures of B-form DNA have provided insights into the global and local conformational properties of the double helix, the solvent environment, drug binding and DNA packing. For example, structures of the duplex with sequence CGCGAATTCGCG, the Dickerson-Drew dodecamer (DDD), established a unique geometry of the central A-tract and a hydration spine in the minor groove. However, our knowledge of the various interaction modes between metal ions and DNA is very limited and almost no information exists concerning the origins of the different effects on DNA conformation and packing exerted by individual metal ions.Crystallization of the DDD duplex in the presence of Mg(2+)and Ca(2+)yields different crystal forms. The structures of the new Ca(2+)-form and isomorphous structures of oligonucleotides with sequences GGCGAATTCGCG and GCGAATTCGCG were determined at a maximum resolution of 1.3 A. These and the 1.1 A structure of the DDD Mg(2+)-form have revealed the most detailed picture yet of the ionic environment of B-DNA. In the Mg(2+)and Ca(2+)-forms, duplexes in the crystal lattice are surrounded by 13 magnesium and 11 calcium ions, respectively.Mg(2+)and Ca(2+)generate different DNA crystal lattices and stabilize different end-to-end overlaps and lateral contacts between duplexes, thus using different strategies for reducing the effective repeat length of the helix to ten base-pairs. Mg(2+)crystals allow the two outermost base-pairs at either end to interact laterally via minor groove H-bonds, turning the 12-mer into an effective 10-mer. Ca(2+)crystals, in contrast, unpair the outermost base-pair at each end, converting the helix into a 10-mer that can stack along its axis. This reduction of a 12-mer into a functional 10-mer is followed no matter what the detailed nature of the 5'-end of the chain: C-G-C-G-A-ellipsis, G-G-C-G-A-ellipsis, or a truncated G-C-G-A-ellipsis Rather than merely mediating close contacts between phosphate groups, ions are at the origin of many well-known features of the DDD duplex structure. A Mg(2+)coordinates in the major groove, contributing to kinking of the duplex at one end. While Ca(2+)resides in the minor groove, coordinating to bases via its hydration shell, two magnesium ions are located at the periphery of the minor groove, bridging phosphate groups from opposite strands and contracting the groove at one border of the A-tract.     

Analysis of local helix bending in crystal structures of DNA oligonucleotides and DNA-protein complexes.

Sequence-dependent bending of the helical axes in 112 oligonucleotide duplex crystal structures resident in the Nucleic Acid Database have been analyzed and compared with the use of bending dials, a computer graphics tool. Our analysis includes structures of both A and B forms of DNA and considers both uncomplexed forms of the double helix as well as those bound to drugs and proteins. The patterns in bending preferences in the crystal structures are analyzed by base pair steps, and emerging trends are noted. Analysis of the 66 B-form structures in the Nucleic Acid Database indicates that uniform trends within all pyrimidine-purine and purine-pyrimidine steps are not necessarily observed but are found particularly at CG and GC steps of dodecamers. The results support the idea that AA steps are relatively straight and that larger roll bends occur at or near the junctions of these A-tracts with their flanking sequences. The data on 16 available crystal structures of protein-DNA complexes indicate that the majority of the DNA bends induced via protein binding are sharp localized kinks. The analysis of the 30 available A-form DNA structures indicates that these structures are also bent and show a definitive preference for bending into the deep major groove over the shallow minor groove.     

A-Tract bending: insights into experimental structures by computational models

While solution structures of adenine tract (A-tract) oligomers have indicated a unique bend direction equivalent to negative global roll (commonly termed "minor-groove bending"), crystallographic data have not unambiguously characterized the bend direction; nevertheless, many features are shared by all A-tract crystal and solution structures (e.g. propeller twisting, narrow minor grooves, and localized water spines). To examine the origin of bending and to relate findings to the crystallographic and solution data, we analyze molecular dynamics trajectories of two solvated A-tract dodecamers: 1D89, d(CGCGA(6)CG), and 1D98, d(CGCA(6)GCG), using a new general global bending framework for analyzing bent DNA and DNA/protein complexes. It is significant that the crystallographically-based initial structures are converted from dissimilar to similar bend directions equivalent to negative global roll, with the average helical-axis bend ranging from 10.5 degrees to 14.1 degrees. The largest bend occurs as positive roll of 12 degrees on the 5' side of the A-tracts (supporting a junction model) and is reinforced by gradual curvature at each A-tract base-pair (bp) step (supporting a wedge model). The precise magnitude of the bend is subtly sequence dependent (consistent with a curved general sequence model). The conversion to negative global roll only requires small local changes at each bp, accumulated over flexible moieties both outside and inside the A-tract. In contrast, the control sequence 1BNA, d(CGCGA(2)TTCGCG), bends marginally (only 6.9 degrees ) with no preferred direction. The molecular features that stabilize the bend direction in the A-tract dodecamers include propeller twisting of AT base-pairs, puckering differences between A and T deoxyriboses, a narrow minor groove, and a stable water spine (that extends slightly beyond the A-tract, with lifetimes approaching 0.2 ns). The sugar conformations, in particular, are proposed as important factors that support bent DNA. It is significant that all these curvature-stabilizing features are also observed in the crystallographic structures, but yield overall different bending paths, largely due to the effects of sequences outside the A-tract. These results merge structural details reported for A-tract structures by experiment and theory and lead to structural and dynamic insights into sequence-dependent DNA flexibility, as highlighted by the effect of an A-tract variant of a TATA-box element on bending and flexibility required for TBP binding.     

Structure and axis curvature in two dA6 x dT6 DNA oligonucleotides: comparison of molecular dynamics simulations with results from crystallography and NMR spectroscopy

Molecular dynamics (MD) simulations have been performed on the A6 containing DNA dodecamers d(GGCAAAAAACGG) solved by NMR and d(CGCAAAAAAGCG) solved by crystallography. The experimental structures differ in the direction of axis bending and in other small but important aspects relevant to the DNA curvature problem. Five nanosecond MD simulations of each sequence have been performed, beginning with both the NMR and crystal forms as well as canonical B-form DNA. The results show that all simulations converge to a common form in close proximity to the observed NMR structure, indicating that the structure obtained in the crystal is likely a strained form due to packing effects. A-tracts in the MD model are essentially straight. The origin of axis curvature is found at pyrimidine-purine steps in the flanking sequences.     

Preferential counterion binding to A-tract DNA oligomers

The free solution mobility of four 20 bp DNA oligomers, with and without A-tracts, has been measured by capillary electrophoresis in Tris-acetate buffer, to test the hypothesis that site-specific binding of monovalent counterions can occur in the narrow minor groove of A-tract DNAs. Preferential counterion binding has been proposed to cause A-tract bending because of asymmetric charge neutralization and collapse of the helix backbone toward the minor groove. Preferential counterion binding in A-tract DNAs should be manifested by a decrease in the electrophoretic mobility observed in free solution, compared to that of non-A-tract DNAs of the same size. Of the four sequences studied here, the slowest absolute mobility, indicative of the greatest counterion binding, was observed for a 20 bp oligomer containing two runs of A3T3 in phase with the helix repeat. A 20-mer containing phased CACA sequences migrated with the fastest mobility; 20-mers containing phased A5 tracts or phased runs of T3A3 migrated with intermediate mobilities. Very similar mobility differences were observed when 1-20 mM NaCl was added to the buffer. The results suggest that preferential counterion binding occurs in A-tract DNAs, especially those containing the AnTn sequence motif.     

Crystallographic study of one turn of G/C-rich B-DNA

The DNA decamer d(CCAGGCCTGG) has been studied by X-ray crystallography. At a nominal resolution of 1.6 A, the structure was refined to R = 16.9% using stereochemical restraints. The oligodeoxyribonucleotide forms a straight B-DNA double helix with crystallographic dyad symmetry and ten base-pairs per turn. In the crystal lattice, DNA fragments stack end-to-end along the c-axis to form continuous double helices. The overall helical structure and, notably, the groove dimensions of the decamer are more similar to standard, fiber diffraction-determined B-DNA than A-tract DNA. A unique stacking geometry is observed at the CA/TG base-pair step, where an increased rotation about the helix axis and a sliding motion of the base-pairs along their long axes leads to a superposition of the base rings with neighboring carbonyl and amino functions. Three-center (bifurcated) hydrogen bonds are possible at the CC/GG base-pair steps of the decamer. In their common sequence elements, d(CCAGGCCTGG) and the related G.A mismatch decamer d(CCAAGATTGG) show very similar three-dimensional structures, except that d(CCAGGCCTGG) appears to have a less regularly hydrated minor groove. The paucity of minor groove hydration in the center of the decamer may be a general feature of G/C-rich DNA and explain its relative instability in the B-form of DNA.     

(+)-CC-1065 as a probe for intrinsic and protein-induced bending of DNA

(+)-CC -1065 is biologically potent DNA-reactive antitumor antibiotic produced by Streptomyces zelensis. This antibiotic covalently modifies DNA by alkylation of N-3 of a adenine in the minor groove. As a Structural consequence of covalent modification of DNA, the helix axis id bent into the minor groove. The drug-induced bending of DNA has similarities to intrinsic. A-tract bending and the 3′ adenine of A-tracts shows a unique reactivity to alkylation by (+) -CC-1065. Upon covalent modification of A-tracts, the magnitude of bending is increased and helix is stiffened. Using high-field NMR, hydroxyl-radical footprinting and gel electrophoresis, the molecular basis for the high reactivity of the bonding sequence 5′ - AGTTA* (an asterisk indicates the covalent modification site) to (+)-CC-1065 has been shown to involve the inherent conformational flexibility of this sequence. Furthermore, these studies also demonstrate that after alkylation the drug-induced bending is focused over the TT region. By analogy with the junction bend model for A-tracts, a ‘truncated junction bend model’ is proposed for this structure. Last, the application of (+)-CC-1065 entrapped/induced bending of DNA as a probe for the Sp1-induced bending of the 21-base-pair repeat an Mu transpose bending of the att L3 sequence is described.     

Structural insights into the effect of hydration and ions on A-tract DNA: a molecular dynamics study

DNA structure is known to be sensitive to hydration and ionic environment. To explore the dynamics, hydration, and ion binding features of A-tract sequences, a 7-ns Molecular dynamics (MD) study has been performed on the dodecamer d(CGCAAATTTGCG)(2). The results suggest that the intrusion of Na(+) ion into the minor groove is a rare event and the structure of this dodecamer is not very sensitive to the location of the sodium ions. The prolonged MD simulation successfully leads to the formation of sequence dependent hydration patterns in the minor groove, often called spine of hydration near the A-rich region and ribbon of hydration near the GC regions. Such sequence dependent differences in the hydration patterns have been seen earlier in the high resolution crystal structure of the Drew-Dickerson sequence, but not reported for the medium resolution structures (2.0 approximately 3.0 A). Several water molecules are also seen in the major groove of the MD simulated structure, though they are not highly ordered over the extended MD. The characteristic narrowing of the minor groove in the A-tract region is seen to precede the formation of the spine of hydration. Finally, the occurrence of cross-strand C2-H2.O2 hydrogen bonds in the minor groove of A-tract sequences is confirmed. These are found to occur even before the narrowing of the minor groove, indicating that such interactions are an intrinsic feature of A-tract sequences.     

Detection of an unusual distortion in A-tract DNA using KMnO4: effect of temperature and distamycin on the altered conformation.

The chemical probes potassium permanganate (KMnO4) and diethylpyrocarbonate (DEPC) can be used to study the conformational flexibility of short tracts of adenine (A-tracts) present in DNA. With these probes, we demonstrate that a novel distortion is induced in a 5 base pair A-tract at low temperature. Formation of this distorted A-tract structure, which occurs in a DNA fragment from the promoter region of the plasmid pBR322, is distinguished by a dramatic increase in the KMnO4 reactivity of the central thymines in this tract at 12 degrees C. This alteration occurs in the absence of any detectable rearrangement in the conformation of the adenines in the complementary strand. Induction of this low temperature A-tract structure is blocked by the minor groove binding drug distamycin. Hydroxyl radical footprinting of distamycin binding to the fragment containing the d(A)5 tract at 12 degrees C suggests that this drug has two different modes of binding to DNA in agreement with recent NMR data. These experiments show that short A-tracts are capable of forming more than one structural variant of B DNA in solution. The possible relationship between the intrinsic bending of DNA containing short phased A-tracts and the low temperature A-tract conformation is discussed.     

DNA Sequence Determinants Controlling Affinity,Stability and Shape of DNA Complexes Bound by the Nucleoid Protein Fis

The abundant Fis nucleoid protein selectively binds poorly related DNA sequences with high affinities to regulate diverse DNA reactions. Fis binds DNA primarily through DNA backbone contacts and selects target sites by reading conformational properties of DNA sequences, most prominently intrinsic minor groove widths. High-affinity binding requires Fis-stabilized DNA conformational changes that vary depending on DNA sequence. In order to better understand the molecular basis for high affinity site recognition, we analyzed the effects of DNA sequence within and flanking the core Fis binding site on binding affinity and DNA structure. X-ray crystal structures of Fis-DNA complexes containing variable sequences in the noncontacted center of the binding site or variations within the major groove interfaces show that the DNA can adapt to the Fis dimer surface asymmetrically. We show that the presence and position of pyrimidine-purine base steps within the major groove interfaces affect both local DNA bending and minor groove compression to modulate affinities and lifetimes of Fis-DNA complexes. Sequences flanking the core binding site also modulate complex affinities, lifetimes, and the degree of local and global Fis-induced DNA bending. In particular, a G immediately upstream of the 15 bp core sequence inhibits binding and bending, and A-tracts within the flanking base pairs increase both complex lifetimes and global DNA curvatures. Taken together, our observations support a revised DNA motif specifying high-affinity Fis binding and highlight the range of conformations that Fis-bound DNA can adopt. The affinities and DNA conformations of individual Fis-DNA complexes are likely to be tailored to their context-specific biological functions.     

Examination of the premelting transition of DNA A-tracts using a fluorescent adenosine analogue

The fluorescent adenosine analogue 4-amino-8-(2-deoxy-beta-d-ribofuranosyl)-5'-O-dimethoxytrityl-6-methyl-7(8H)-pteridone (6MAP) has been used to perform residue specific analyses of DNA A-tracts during the premelting transition. DNA A-tracts, which exhibit sequence-induced curvature, adopt a B-DNA conformation as a function of increasing temperature. Fluorescence melting curves indicate that 6MAP is a more sensitive reporter of the premelting transition than UV absorption spectroscopy. Further, residue specific fluorescence analyses of A-tract and control duplexes reveal that some of the conformational changes associated with the premelting transition occur within A-tract regions. Analyses of the energetics of the premelting transition indicate that ApA steps make a larger enthalpic contribution to the premelting transition than ApT steps. To explore the effect of cations on the premelting transition, fluorescence melts were performed in the presence of NH(4)(+), Mg(2+), and low (0.05 M) and high (0.5 M) concentrations of Na(+). These studies show that the fluorescence intensity changes associated with the premelting transition are sensitive to cation type and concentration and are larger and more pronounced in the presence of 0.5 M Na(+), NH(4)(+), and Mg(2+). Incorporation of 6MAP into longer duplexes containing phased A-tracts shows that the local environment of adenosines in phased A-tracts is similar to that of individual A-tracts. Fluorescence quenching results indicate that ApA and ApT steps within A-tracts are less solvent exposed than their counterparts in control sequence isomers, possibly because of the narrowed minor groove of A-tract sequences.     

Flanking A·T Basepairs Destabilize the B? Conformation of DNA A-Tracts

Capillary electrophoresis has been used to characterize the interaction of monovalent cations with 26-basepair DNA oligomers containing A-tracts embedded in flanking sequences with different basepair compositions. A 26-basepair random-sequence oligomer was used as the reference; lithium and tetrabutylammonium (TBA⁺) ions were used as the probe ions. The free solution mobilities of the A-tract and random-sequence oligomers were identical in solutions containing <∼100 mM cation. At higher cation concentrations, the A-tract oligomers migrated faster than the reference oligomer in TBA⁺ and slower than the reference in Li⁺. Hence, cations of different sizes can interact very differently with DNA A-tracts. The increased mobilities observed in TBA⁺ suggest that the large hydrophobic TBA⁺ ions are preferentially excluded from the vicinity of the A-tract minor groove, increasing the effective net charge of the A-tract oligomers and increasing the mobility. By contrast, Li⁺ ions decrease the mobility of A-tract oligomers because of the preferential localization of Li⁺ ions in the narrow A-tract minor groove. Embedding the A-tracts in AT-rich flanking sequences markedly alters preferential interactions of monovalent cations with the B^∗ conformation. Hence, A-tracts embedded in genomic DNA may or may not interact preferentially with monovalent cations, depending on the relative number of A·T basepairs in the flanking sequences.     

DNA sequence-dependent positioning of the linker histone in a nucleosome: A single-pair FRET study

Linker histones (LHs) bind to nucleosomes with their globular domain (gH) positioned in either an on- or an off-dyad binding mode. Here, we study the effect of the linker DNA (L-DNA) sequence on the binding of a full-length LH, Xenopus laevis H1.0b, to a Widom 601 nucleosome core particle (NCP) flanked by two 40 bp long L-DNA arms, by single-pair FRET spectroscopy. We varied the sequence of the 11 bp of L-DNA adjoining the NCP on either side, making the sequence either A-tract, purely GC, or mixed with 64% AT. The labeled gH consistently exhibited higher FRET efficiency with the labeled L-DNA containing the A-tract than that with the pure-GC stretch, even when the stretches were swapped. However, it did not exhibit higher FRET efficiency with the L-DNA containing 64% AT-rich mixed DNA when compared to the pure-GC stretch. We explain our observations with a model that shows that the gH binds on dyad and that two arginines mediate recognition of the A-tract via its characteristically narrow minor groove. To investigate whether this on-dyad minor groove-based recognition was distinct from previously identified off-dyad major groove-based recognition, a nucleosome was designed with A-tracts on both the L-DNA arms. One A-tract was complementary to thymine and the other to deoxyuridine. The major groove of the thymine-tract was lined with methyl groups that were absent from the major groove of the deoxyuridine tract. The gH exhibited similar FRET for both these A-tracts, suggesting that it does not interact with the thymine methyl groups exposed on the major groove. Our observations thus complement previous studies that suggest that different LH isoforms may employ different ways of recognizing AT-rich DNA and A-tracts. This adaptability may enable the LH to universally compact scaffold-associated regions and constitutive heterochromatin, which are rich in such sequences.     

Studies of base pair sequence effects on DNA solvation based on all-atom molecular dynamics simulations

Detailed analyses of the sequence-dependent solvation and ion atmosphere of DNA are presented based on molecular dynamics (MD) simulations on all the 136 unique tetranucleotide steps obtained by the ABC consortium using the AMBER suite of programs. Significant sequence effects on solvation and ion localization were observed in these simulations. The results were compared to essentially all known experimental data on the subject. Proximity analysis was employed to highlight the sequence dependent differences in solvation and ion localization properties in the grooves of DNA. Comparison of the MD-calculated DNA structure with canonical A- and B-forms supports the idea that the G/C-rich sequences are closer to canonical A- than B-form structures, while the reverse is true for the poly A sequences, with the exception of the alternating ATAT sequence. Analysis of hydration density maps reveals that the flexibility of solute molecule has a significant effect on the nature of observed hydration. Energetic analysis of solute-solvent interactions based on proximity analysis of solvent reveals that the GC or CG base pairs interact more strongly with water molecules in the minor groove of DNA that the AT or TA base pairs, while the interactions of the AT or TA pairs in the major groove are stronger than those of the GC or CG pairs. Computation of solvent-accessible surface area of the nucleotide units in the simulated trajectories reveals that the similarity with results derived from analysis of a database of crystallographic structures is excellent. The MD trajectories tend to follow Manning's counterion condensation theory, presenting a region of condensed counterions within a radius of about 17 A from the DNA surface independent of sequence. The GC and CG pairs tend to associate with cations in the major groove of the DNA structure to a greater extent than the AT and TA pairs. Cation association is more frequent in the minor groove of AT than the GC pairs. In general, the observed water and ion atmosphere around the DNA sequences is the MD simulation is in good agreement with experimental observations.     

Sequence-dependent variation in DNA minor groove width dictates orientational preference of Hoechst 33258 in A-tract recognition: solution NMR structure of the 2:1 complex with d(CTTTTGCAAAAG)(2)

The solution structure of the dodecamer duplex d(CTTTTGCAAAAG)₂ and its 2:1 complex with the bis-benzimidazole Hoechst 33258 has been investigated by NMR and NOE-restrained molecular dynamics (rMD) simulations. Drug molecules are bound in each of the two A-tracts with the bulky N-methylpiperazine ring of each drug located close to the central TG (CA) step, binding essentially to the narrow minor groove of each A-tract. MD simulations over 1 ns, using an explicit solvation model, reveal time-averaged sequence-dependent narrowing of the minor groove from the 3′-end towards the 5′-end of each TTTT sequence. Distinct junctions at the TpG (CpA) steps, characterised by large positive roll, low helical and propeller twists and rapid AT base pair opening rates, add to the widening of the groove at these sites and appear to account for the bound orientation of the two drug molecules with the N-methylpiperazine ring binding in the wider part of the groove close to the junctions. Comparisons between the free DNA structure and the 2:1 complex (heavy atom RMSD 1.55 Å) reveal that these sequence-dependent features persist in both structures. NMR studies of the sequence d(GAAAAGCTTTTC)₂, in which the A-tracts have been inverted with the elimination of the TpG junctions, results in loss of orientational specificity of Hoechst 33258 and formation of multiple bound species in solution, consistent with the drug binding in a number of different orientations.     

MPD and DNA Bending in Crystals and in Solution

Bending of 15 to 24° is observed within crystal structures ofB-DNA duplexes, is strongly sequence-dependent, and exhibits no correlation with the concentration of MPD (2-methyl-2,4-pentanediol) in the crystallizing solution. Two types of bends are observed: facultative bends or flexible hinges at junctions between regions of G·C and A·T base-pairs, and a persistent and almost obligatory bend at the center of the sequence R-G-C-Y. Only A-tracts are characteristically straight and unbent in every crystal structure examined to date. A detailed examination of normal vector plots for individual strands of a double helix provides an explanation, in terms of the stacking properties of guanine and adenine bases. The effect of high MPD concentrations, in both solution and crystal, is to decrease local bending somewhat without removing it altogether. MPD gel retardation experiments provide no basis for choosing among the three models that seek to explain macroscopic curvature of DNA by means of microscopic bending: junction bending, bent A-tracts, or bent general- sequence DNA. Crystallographic data on the straightness of A-tracts, the bendability of non-A sequences, and the identity of inclination angles in A-tract and non-A-tractB-DNA support only the general-sequence bending model. The pre-melting transition observed in A-tract DNA probably represents a relaxation of stiff adenine stacks to a flexible conformation more typical of general-sequence DNA.