Sequence-Dependent Bending of DNA and Phasing of Nucleosomes

Abstract

Conformational analysis has revealed anisotropic flexibility of the B-DNA double helix: it bends most easily into the grooves, being the most rigid when bent in a perpendicular direction. This result implies that DNA in a nucleosome is curved by means of relatively sharp bends (“mini-kinks”) which are directed into the major and minor grooves alternatively and separated by 5–6 base pairs. The “mini-kink” model proved to be in keeping with the x-ray structure of the B-DNA dodecamer resolved later, which exhibits two “annealed kinks”, also directed into the grooves.

The anisotropy of B DNA is sequence-dependent: the pyrimidine-purine dimers (YR) favor bending into the minor groove, and the purine-pyrimidine dinucleotides (RY), into the minor one. The RR and YY dimers appear to be the most rigid dinucleotides. Thus, a DNA fragment consisting of the interchanging oligopurine and oligopyrmidine blocks 5–6 base pairs long should manifest a spectacular curvature in solution.

Similarly, a nucleotide sequence containing the RY and YR dimers separated by a half-pitch of the double helix is the most suitable for wrapping around the nucleosomal core. Analysis of the numerous examples demonstrating the specific alignment of nucleosomes on DNA confirms this concept. So, the sequence-dependent “mechanical” properties of the double helix influence the spatial arrangement of DNA in chromatin.     

Anisotropic flexibility of DNA depends on the base sequence. Conformation calculations of double-stranded tetranucleotides AAAA:TTTT, (AATT)2, (TTAA)2, GGGG:CCCC, (GGCC)2, (CCGG)2

The bending flexibility of six tetramers was studied in an assumption that they were extended in the both directions by regular double helices. The bends of B-DNA in different directions were considered. The stiffness of the B-DNA double helix when bent into the both grooves proved to be less pronounced than in the perpendicular direction by the order of magnitude. Such an anisotropy is a feature of the sugar-phosphate backbone structure. The calculated fluctuations of the DNA bending along the dyad axis, 5-7 degrees, are in agreement with the experimental value of DNA persistence length. Anisotropy of the double helix is sequence-dependent: most easily bent into the minor groove are the tetramers with purine-pyrimidine dimer (RY) in the middle. In contrast, YR dinucleotides prefer bending into the major groove, moreover, they have an equilibrium bend of 6-12 degrees into this groove. The above inequality is caused by the stacking interaction of the bases. The bend in the central dimers is distributed to some extent between the adjacent links, though the main fraction of the bend remains within the central link. Variation of the sugar-phosphate geometry in the bent helix is unessential, so that DNA remains within the limits of the B-family of forms: namely, when the helical axis is bent by 20 degrees the backbone dihedral angles vary by no more than 15 degrees. The obtained results are in accord with the X-ray structure of B-DNA dodecamer; they further substantiate our earlier model of DNA wrapping in the nucleosome by means of "mini-kinks" separated by a half-pitch of the double helix, i.e. by 5-6 b. p. Sequence-dependent anisotropy of DNA presumably dictates the three-dimensional structure of DNA in solution as well. We have found that nonrandom allocation of YR dimers leads to the systematic bends in the equilibrium structure of certain DNA fragments. To the four "Calladine rules" two more can be added: the minor-groove steric clash of purines in the YR sequences are avoided by: (1) bending of the helix into the major groove; (2) increasing the distance between the base pairs (stretching the double helix).     

Sequence-dependent anisotropic flexibility of B-DNA. A conformational study

Bending flexibility of the six tetrameric duplexes was investigated d(AAAA):d(TTTT), d(AATT)2, d(TTAA)2, d(GGGG):d(CCCC), d(GGCC)2 and d(CCGG)2,. The tetramers were extended in the both directions by regular double helices. The stiffness of the B-DNA double helix when bent into the both grooves proved to be less than that in the perpendicular direction by an order of magnitude. Such an anisotropy is a property of the sugar-phosphate backbone structure. The calculated fluctuations of the DNA bending along the dyad axis, 5-7 degree, are in agreement with experimental value of the DNA persistence length. Anisotropy of the double helix is sequence-dependent: most easily bent into the minor groove are the tetramers with purine-pyrimidine dimer (RY) in the middle. In contrast, YR dinucleotides prefer bending into the major groove. Moreover, they have an equilibrium bend of 6-12 degree into this groove. The above inequality is caused by stacking interaction of the bases. The bend in the central dimer is distributed to some extent between the adjacent links, though the main fraction of the bend remains within the central link. Variation of the sugar-phosphate geometry in the bent helix is inessential, so that DNA remains within the B-family of forms: namely, when the helical axis is bent by 20 degree. the backbone dihedral angles vary by no more than 15 degree. The obtained results are in accord with x-ray structure of the B-DNA dodecamer; they further substantiate our early model of DNA wrapping in the nucleosome by means of "mini-kinks" separated by a half-pitch of the double helix, i.e. by 5-6 b.p. Sequence-dependent anisotropy of DNA presumably dictates the three-dimensional structure of DNA in solution as well. We have found that nonrandom allocation of YR dimers leads to the systematic bends in equilibrium structure of certain DNA fragments.     

Sequence-Dependent Anisotropic Flexibility of B-DNA A Conformational Study

Abstract

Bending flexibility of the six tetrameric duplexes was investigated d(AAAA):d(TTTT), d(AATT)₂, d(TTAA) ₂, d(GGGG):d(CCCC), d(GGCC) ₂ and d(CCGG) ₂. The tetramers were extended in the both directions by regular double helices. The stiffness of the B-DNA double helix when bent into the both grooves proved to be less than that in the perpendicular direction by an order of magnitude. Such an anisotropy is a property of the sugar-phosphate backbone structure. The calculated fluctuations of the DNA bending along the dyad axis, 5–7°, are in agreement with experimental value of the DNA persistence length.

Anisotropy of the double helix is sequence-dependent: most easily bent into the minor groove are the tetramers with purine-pyrimidine dimer (RY) in the middle. In contrast, YR dinucleotides prefer bending into the major groove. Moreover, they have an equilibrium bend of 6–12° into this groove. The above inequality is caused by stacking interaction of the bases.

The bend in the central dimer is distributed to some extent between the adjacent links, though the main fraction of the bend remains within the central link. Variation of the sugar-phosphate geometry in the bent helix is inessential, so that DNA remains within the B-family of forms: namely, when the helical axis is bent by 20°, the backbone dihedral angles vary by no more than 15°.

The obtained results are in accord with x-ray structure of the B-DNA dodecamer; they further substantiate our early model of DNA wrapping in the nucleosome by means of “mini-kinks” separated by a half-pitch of the double helix, i.e. by 5–6 b.p. Sequence-dependent anisotropy of DNA presumably dictates the three-dimentional structure of DNA in solution as well. We have found that nonrandom allocation of YR dimers leads to the systematic bends in equilibrium structure of certain DNA fragments.     

Configurational statistics of the DNA duplex: Extended generator matrices to treat the rotations and translations of adjacent residues

The base-to-base virtual bond treatment of nucleic acids used in statistical mechanical calculations of polynucleotide chain properties has been refined by incorporating the six parameters that relate the positions and orientations of sequential rigid bodies. The scheme allows for the sequence-dependent bending, twisting, and displacement of base pairs as well as for asymmetry in the angular and translational fluctuations of individual residues. Expressions are developed for the generator matrices required for the computation, as a function of chain length, of various parameters measuring the overall mean extension and shape of the DNA. Quantities of interest include the end-to-end vector r , the square of the end-to-end distance r², the square radius of gyration s², the center-of-gravity vector g , the second moments of inertia S ^×2, and the higher moments of r and g . The matrix expressions introduced in the 1960s by Flory and co-workers for the determination of configuration-dependent polymer chain averages are decomposed into their translational and orientational contributions so that the methods can be extended to the rigid body analysis of chemical moieties. The new expressions permit, for the first time, examination of the effects of sequence-dependent translations, such as the lateral sliding of residues in A- and B-helices and the vertical opening of base pairs in drug–DNA complexes, on the average extension and shape of the long flexible double helix. The approach is illustrated in the following paper using conformational energy estimates of the base sequence-dependent flexibility of successive B-DNA base pairs. © 1994 John Wiley & Sons, Inc.     

Distamycin-stabilized Antiparallel-Parallel Combination (APC) DNA

Abstract

The formation of Antiparallel-Parallel-Combination (APC) DNA, a liner duplex with a segment of parallel-stranded (ps) helix flanked by conventional B-DNA, was tested with a number of synthetic oligonucleotides. The groove-binding ligand distamycin A (DstA) was used to stabilize the ps segment comprising five A·T base pairs. Two drug molecules bound per APC, one in each of the two equivalent grooves characteristic of ps-DNA. APC-DNA, reference molecules and their complexes with DstA were analysed by several methods: circular dichroism and absorption spectroscopy, thermal denaturation, chemical modification, and molecular modeling. The dye binding stoichiometry differed significantly due to inherent structural differences in the groove geometries of ps-DNA (trans base pairs, similar grooves) and conventional antiparallel-stranded (aps) B-DNA (cis base pairs, distinct major and minor grooves). The data support the existence of APC folding in solution.     

Base pair buckling can eliminate the interstrand purine clash at the CpG steps in B-DNA caused by the base pair propeller twisting

Results of calculations using various empirical potentials suggest that base pair buckling, which commonly occurs in DNA crystal structures, is sufficient to eliminate the steric clash at CpG steps in B-DNA, originating from the base pair propeller twisting. The buckling is formed by an inclination of cytosines while deviations of guanines from a plane perpendicular to the double helix axis are unfavorable. The buckling is accompanied by an increased vertical separation of the base pair centers but the buckled arrangement of base pairs is at least as stable as when the vertical separation is normal and buckle zero. In addition, room is created by the increased vertical separation for the bases to propeller twist as is observed in DNA crystal structures. Further stabilization of base stacking is introduced into the buckled base pair arrangement by roll opening the base pairs into the double helix minor groove. The roll may lead to the double helix bending and liberation of guanines from the strictly perpendicular orientation to the double helix axis. The liberated guanines further contribute to the base pair buckling and stacking improvement. This work also suggests a characteristic very stable DNA structure promoted by nucleotide sequences in which runs of purines follow runs of pyrimidine bases.     

Antitumor drug nogalamycin binds DNA in both grooves simultaneously: molecular structure of nogalamycin-DNA complex

The three-dimensional molecular structures of the complexes between an interesting antitumor drug, nogalamycin, and two DNA hexamers, d[CGT(pS)ACG] and d[m5CGT(pS)Am5CG], were determined at high resolution by X-ray diffraction analyses. Two nogalamycins bind to the DNA double helix in a 2:1 ratio with the aglycon chromophore intercalated between the CpG steps at both ends of the helix. The nogalose and aminoglucose sugars lie in the minor and major grooves, respectively, of the distorted B-DNA double helix. The binding of nogalamycin to DNA requires that the base pairs in DNA open up transiently to allow the bulky sugars to go through. Specific hydrogen bonds are found in the complex between the drug and guanine bases. We suggest that nogalamycin may prefer GC sequences embedded in a stretch of AT sequences.     

Analyzing ion distributions around DNA: sequence-dependence of potassium ion distributions from microsecond molecular dynamics

Microsecond molecular dynamics simulations of B-DNA oligomers carried out in an aqueous environment with a physiological salt concentration enable us to perform a detailed analysis of how potassium ions interact with the double helix. The oligomers studied contain all 136 distinct tetranucleotides and we are thus able to make a comprehensive analysis of base sequence effects. Using a recently developed curvilinear helicoidal coordinate method we are able to analyze the details of ion populations and densities within the major and minor grooves and in the space surrounding DNA. The results show higher ion populations than have typically been observed in earlier studies and sequence effects that go beyond the nature of individual base pairs or base pair steps. We also show that, in some special cases, ion distributions converge very slowly and, on a microsecond timescale, do not reflect the symmetry of the corresponding base sequence.     

Structure of the nucleosome core particle at 8 A resolution

The x-ray crystallographic structure of the nucleosome core particle has been determined using 8 A resolution diffraction data. The particle has a mean diameter of 106 A and a maximum thickness of 65 A in the superhelical axis direction. The longest chord through the histone core measures 85 A and is in a non-axial direction. The 1.87 turn superhelix consists of B-DNA with about 78 base pairs or 7.6 helical repeats per superhelical turn. The mean DNA helical repeat contains 10.2 +/- 0.05 base pairs and spans 35 A, slightly more than standard B-DNA. The superhelix varies several Angstroms in radius and pitch, and has three distinct domains of curvature (with radii of curvature of 60, 45 and 51 A). These regions are separated by localized sharper bends +/- 10 and +/- 40 base pairs from the center of the particle, resulting in an overall radius of curvature about 43 A. Compression of superhelical DNA grooves on the inner surface and expansion on the outer surface can be seen throughout the DNA electron density. This density has been fit with a double helical ribbon model providing groove width estimates of 12 +/- 1 A inside vs. 19 +/- 1 A outside for the major groove, and 8 +/- 1 A inside vs. 13 +/- 1 A outside for the minor groove. The histone core is primarily contained within the bounds defined by the superhelical DNA, contacting the DNA where the phosphate backbone faces in toward the core. Possible extensions of density between the gyres have been located, but these are below the significance level of the electron density map. In cross-section, a tripartite organization of the histone octamer is apparent, with the tetramer occupying the central region and the dimers at the extremes. Several extensions of histone density are present which form contacts between nucleosomes in the crystal, perhaps representing flexible or "tail" histone regions. The radius of gyration of the histone portion of the electron density is calculated to be 30.4 A (in reasonable agreement with solution scattering values), and the histone core volume in the map is 93% of its theoretical volume.     

Bends in SV40 DNA: use of mutagenesis to identify the critical bases involved.

Five fragments of DNA exhibiting sequence directed bends were isolated from the Simian Virus 40 genome using a two-dimensional polyacrylamide gel fractionation. The bend sites were mapped for each fragment using the circular permutation test. All five sites have multiple, short runs of A residues with helical spacing typical of other bent fragments. Base pairs important for the bends were determined for one fragment by utilizing a random, single base pair mutagenesis. Of 28 mutants with decreased or increased bends, 14 had alterations that could be interpreted to affect the spaced runs of A residues, supporting their role in bends as predicted by the ApA wedge model. One major mutation was not explainable by existing models. The remaining minor mutations may only be due to small, local DNA conformational changes in the surrounding B-DNA.     

Structural principles of B-DNA grooves hydration in fibers as revealed by Monte Carlo simulations and X-ray diffraction

Being interested in possible effects of sequence-dependent hydration of B-DNA with mixed sequence in fibers, we performed a series of Monte Carlo calculations of hydration of polydeoxyribonucleotides in B form, considering all sequences with dinucleotide repeat. The computational results allow the ten base-stacking types to be classified in accordance with their primary hydration in the minor groove. As a rule, the minor groove is occupied by two water molecules per base pair in the depth of the groove, which are located nearly midway between the planes of successive base pairs and symmetrically according to the dyad there. The primary hydration of the major groove depends on the type of the given base pair. The coordinates of 3 water molecules per base pair in the depth of the major groove are determined by the type of this pair together with its position and orientation in the helix, and are practically independent on the adjacent base pairs. A/T-homopolymer tracts do not fit into this hydration pattern; the base pair edges are hydrated autonomously in both grooves. Analysis of the Li-B-DNA x-ray diffraction intensities reveals those two water positions in the minor groove. In the major groove, no electronic density peaks in sufficient distance from the base edges were found, thus confirming the absence of any helical invariance of primary hydration in this region. With the help of the rules proposed in this paper it is possible to position the water molecules of the first hydration shell in the grooves of canonical B-DNA for any given sequence.     

CAP binding sites reveal pyrimidine-purine pattern characteristic of DNA bending

To investigate the intrinsic bending of DNA at sites where proteins bind, we analyzed catabolite gene activator protein (CAP) binding sites and various operators from the viewpoint of DNA bending flexibility. Theoretical conformational analysis. DNase I digestion and x-ray crystallography data indicate that bending of B-DNA is highly anisotropic and sequence-dependent. Certain dimers prefer to bend into the major groove ("major-philic") and others prefer to bend into the minor groove ("minor-philic" dimers). From these data we considered TA, CG, CA:TG and GG:CC as major-philic dimers and AT,AA:TT and GT:AC as minor-philic ones. Analysis of 31 CAP binding sites has identified strong major-philic tendencies 5-7 base pairs (bp) away from the center. In addition, we found minor-philic poly-A tracts extending 4-5 bp away from the proposed major-philic bends. Finally, to analyze the central regions we followed the lead of Shumilov and classified the DNA sites by their spacer lengths [V.Y. Shumilov, Mol. Biol. (Mosk) 21, 168-187 (1987)]. In this way, we identified two subsets of CAP binding sites: one with 6 bp between the TGTGA:TCACA consensus boxes (N6-set) and one with 8 central bp (N8-set). We discovered that the dimer at the center of an N6-set site was usually major-philic, whereas at the center of an N8-set site more often minor-philic. Analysis of phages 434, P22 lambda and trp operators revealed similar results. In conclusion, our data show that CAP binding sites have major-philic and minor-philic dimers at specific positions; the location of these dimers may facilitate wrapping of DNA around CAP. A similar pattern is seen in nucleosomes.     

A Possible Role of DNA Superstructures in Genome Evolution

The concept of DNA as a simple repository of the gene information has changed in that of a polymorphic macromolecule, which plays a relevant part in the management of the complex biochemical transformations in living matter. As a consequence of the slight stereochemical differences between base pairs, the direction of the DNA double helix axis undergoes deterministic writhing. A useful representation of such sequence-dependent structural distortions is the curvature diagram. Here, it is reported as an evolution simulation obtained by extensive point mutations along a biologically important DNA tract. The curvature changes, consequence of the point mutations. were compared to the related experimental gel electrophoresis mobility. The curvature of most mutants decreases and the mobility increases accordingly, suggesting the curvature of that tract is genetically selected. Moreover, DNA images by scanning force microscopy, show evidence of a sequence-dependent adhesion of curved DNA tracts to inorganic crystal surfaces. In particular, mica shows a large affinity towards the TT-rich dinucleotide sequences. This suggests a possible mechanism of selection of curved DNA regions, characterized by AA.TT dinucleotides in phase with double-helical periodicity, in the very early evolution steps.     

Spatial translational motions of base pairs in DNA molecules: Application of the extended matrix generator method

We have used the elementary generator matrices outlined in the preceding paper to examine the conformational plasticity of the nucleic acid double helix. Here we investigate kinked DNA structures made up of alternating B- and A-type helices and intrinsically curved duplexes perturbed by the intercalation of ligands. We model the B-to-A transition by the lateral translation of adjacent base pairs, and the intercalation of ligands by the vertical displacement of neighboring residues. We report a complete set of average configuration-dependent parameters, ranging from scalars (i.e., persistence lengths) to first- and second-order tensor parameters (i.e., average second moments of inertia), as well as approximations of the associated spatial distributions of the DNA and their angular correlations. The average structures of short chains (of lengths less than 100 base pairs) with local kinks or intrinsically curved sequences are essentially rigid rods. At the smallest chain lengths (10 base pairs), the kinked and curved chains exhibit similar average properties, although they are structurally perturbed compared to the standard B-DNA duplex. In contrast, at lengths of 200 base pairs, the curved and kinked chains are more compact on average and are located in a different space from the standard B- or A-DNA helix. While A-DNA is shorter and thicker than B-DNA in x-ray models, the long flexible A-DNA helix is thinner and more extended on average than its B-DNA counterpart because of more limited fluctuations in local structure. Curved polymers of 50 base pairs or longer also show significantly greater asymmetry than other DNAs (in terms of the distribution of base pairs with respect to the center of gravity of the chain). The intercalation of drugs in the curved DNA straightens and extends the smoothly deformed template. The dimensions of the average ellipsoidal boundaries defining the configurations of the intercalated polymers are roughly double those of the intrinsically curved chain. The altered proportions and orientations of these density functions reflect the changing shape and flexibility of the double helix. The calculations shed new light on the possible structural role of short A-DNA fragments in long B-type duplexes and also offer a model for understanding how GC-specific intercalative ligands can straighten naturally curved DNA. The mechanism is not immediately obvious from current models of DNA curvature, which attribute the bending of the chain to a perturbed structure in repeating tracts of A · T base pairs. © 1994 John Wiley & Sons, Inc.     

Binding of the antitumor drug nogalamycin and its derivatives to DNA: structural comparison

The three-dimensional molecular structures of the complexes between a novel antitumor drug nogalamycin and its derivative U-58872 with a modified DNA hexamer d[m5CGT(pS)Am5CG] have been determined at 1.7- and 1.8-A resolution, respectively, by X-ray diffraction analyses. Both structures (in space group P6(1)) have been refined with constrained refinement procedure to final R factors of 0.208 (3386 reflections) and 0.196 (2143 reflections). In both complexes, two nogalamycins bind to the DNA hexamer double helix in a 2:1 ratio with the elongated aglycon chromophore intercalated between the CpG steps at both ends of the helix. The aglycon chromophore spans across the GC Watson-Crick base pairs with its nogalose lying in the minor groove and the aminoglucose lying in the major groove of the distorted B-DNA double helix. Most of the sugars remain in the C2'-endo pucker family, except three deoxycytidine residues (terminal C1, C7, and internal C5). All nucleotides are in the anti conformation. Specific hydrogen bonds are found in the complex between the drug and guanine-cytosine bases in both grooves of the helix. One hydroxyl group of the aminoglucose donates a hydrogen bond to the N7 of guanine, while the other receives a hydrogen bond from the N4 amino group of cytosine. The orientation of these two hydrogen bonds suggests that nogalamycin prefers a GC base pair with its aglycon chromophore intercalating at the 5'-side of a guanine (between NpG), or at the 3'-side of a cytosine (between CpN) with the sugars pointing toward the GC base pair. The binding of nogalamycin to DNA requires that the base pairs in DNA open up transiently to allow the bulky sugars to go through, suggesting that nogalamycin prefers GC sequences embedded in a stretch of AT sequences.     

CAP Binding Sites Reveal Pyrimidine-Purine Pattern Characteristic of DNA Bending

Abstract

To investigate the intrinsic bending of DNA at sites where proteins bind, we analyzed catabolite gene activator protein (CAP) binding sites and various operators from the viewpoint of DNA bending flexibility. Theoretical conformational analysis, DNase I digestion and x-ray crystallography data indicate that bending of B-DNA is highly anisotropic and sequence-dependent. Certain dimers prefer to bend into the major groove (“major-philic”) and others prefer to bend into the minor groove (“minor-philic” dimers). From these data we considered TA CG, CA:TG and GG:CC as major-philic dimers and AT, AA:TT and GT:AC as minor-philic ones.

Analysis of 31 CAP binding sites has identified strong major-philic tendencies 5–7 base pairs (bp) away from the center. In addition, we found minor-philic poly-A tracts extending 4–5 bp away from the proposed major-philic bends. Finally, to analyze the central regions we followed the lead of Shumilov and classified the DNA sites by their spacer lengths [V.Y. Shumilov, Mol. Biol. (Mosk) 21, 168–187 (1987)]. In this way, we identified two subsets of CAP binding sites: one with 6 bp between the TGTGA:TCACA consensus boxes (N6-set) and one with 8 central bp (N8-set). We discovered that the dimer at the center of an N6-set site was usually major-philic, whereas at the center of an N8-set site more often minor-philic. Analysis of phages 434, P22, λ and trp operators revealed similar results.

In conclusion, our data show that CAP binding sites have major-philic and minor-philic dimers at specific positions; the location of these dimers may facilitate wrapping of DNA around CAP. A similar pattern is seen in nucleosomes.     

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.