Preferred positions of AA and TT dinucleotides in aligned nucleosomal DNA sequences.

Multiple alignment of 118 nucleosomal DNA sequences by maximizing simultaneously match of AA dinucleotides and match of TT dinucleotides results in a pattern of the dinucleotide distributions which is characteristic of the nucleosomal DNA sequences. The AA dinucleotides are found to be distributed symmetrically relative to the TT dinucleotide distribution, around the middle point of the nucleosomal DNA sequence. The distances between major peaks of the distributions are multiples of about 10.4 bases. The peaks of the TT distribution are shifted by 6 bases downstream from the peaks of the AA distribution.     

Characteristics of nucleosome core DNA and their applications in predicting nucleosome positions

By analyzing dinucleotide position-frequency data of yeast nucleosome-bound DNA sequences, dinucleotide periodicities of core DNA sequences were investigated. Within frequency domains, weakly bound dinucleotides (AA, AT, and the combinations AA-TT-TA and AA-TT-TA-AT) present doublet peaks in a periodicity range of 10-11 bp, and strongly bound dinucleotides present a single peak. A time-frequency analysis, based on wavelet transformation, indicated that weakly bound dinucleotides of core DNA sequences were spaced smaller (∼10.3 bp) at the two ends, with larger (∼11.1 bp) spacing in the middle section. The finding was supported by DNA curvature and was prevalent in all core DNA sequences. Therefore, three approaches were developed to predict nucleosome positions. After analyzing a 2200-bp DNA sequence, results indicated that the predictions were feasible; areas near protein-DNA binding sites resulted in periodicity profiles with irregular signals. The effects of five dinucleotide patterns were evaluated, indicating that the AA-TT pattern exhibited better performance. A chromosome-scale prediction demonstrated that periodicity profiles perform better than previously described, with up to 59% accuracy. Based on predictions, nucleosome distributions near the beginning and end of open reading frames were analyzed. Results indicated that the majority of open reading frames’ start and end sites were occupied by nucleosomes.     

The intrinsic curvature of DNA in solution

We propose a detailed quantitative scheme for explaining the anomalous electrophoretic mobility in polyacrylamide gels of repeating sequence DNA. We assume that such DNA adopts a superhelical configuration in these circumstances, and migrates less quickly than straight DNA of the same length because it can only pass through larger holes. The retardation is maximal when the length of the DNA reaches one superhelical turn, but is less for shorter pieces. We attribute the curvature of the superhelix to different angles of roll at each kind of dinucleotide step, i.e. an opening up of an angle by an increased separation on the minor-groove side. The main effect is due to a difference of about 3 degrees in roll values between AA/TT and other steps, together with a difference of about 1 degree in the angle of helical twist: we deduce these values explicitly from some of the available data on gel-running. The scheme involves a simple calculation of the superhelical parameters for any given repeating sequence, and it gives a good correlation with all of the available data. We argue that these same base-step angular parameters are also consistent with observations from X-ray diffraction of crystallized oligomers, and particularly with the recent data on CGCA6GCG from Nelson et al. We are concerned here with the intrinsic curvature of unconstrained DNA, as distinct from the curvature of DNA in association with protein molecules; and this paper represents a first attempt at an absolute determination.     

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.     

Improved alignment of nucleosome DNA sequences using a mixture model

DNA sequences that are present in nucleosomes have a preferential approximately 10 bp periodicity of certain dinucleotide signals, but the overall sequence similarity of the nucleosomal DNA is weak, and traditional multiple sequence alignment tools fail to yield meaningful alignments. We develop a mixture model that characterizes the known dinucleotide periodicity probabilistically to improve the alignment of nucleosomal DNAs. We assume that a periodic dinucleotide signal of any type emits according to a probability distribution around a series of 'hot spots' that are equally spaced along nucleosomal DNA with 10 bp period, but with a 1 bp phase shift across the middle of the nucleosome. We model the three statistically most significant dinucleotide signals, AA/TT, GC and TA, simultaneously, while allowing phase shifts between the signals. The alignment is obtained by maximizing the likelihood of both Watson and Crick strands simultaneously. The resulting alignment of 177 chicken nucleosomal DNA sequences revealed that all 10 distinct dinucleotides are periodic, however, with only two distinct phases and varying intensity. By Fourier analysis, we show that our new alignment has enhanced periodicity and sequence identity compared with center alignment. The significance of the nucleosomal DNA sequence alignment is evaluated by comparing it with that obtained using the same model on non-nucleosomal sequences.     

Experimental evaluation of the Liu–Beveridge dinucleotide step model of DNA structure

Methods for predicting DNA curvature have many possible applications. Dinucleotide step models describe DNA shape by characterization of helical twist, deflection angles and the direction of deflection for nearest neighbor base pairs. Liu and Beveridge have extended previous applications of dinucleotide step models with the development and qualitative validation of a predictive method for sequence-dependent DNA curvature (the LB model). We tested whether the LB model accurately predicts experimentally deduced curvature angles and helical repeat parameters for DNA sequences not in its training set, particularly when challenged with quantitative data and subtle sequence phasings. We examined a series of 17 well-characterized DNA sequences to compare electrophoretic and computational results. The LB model is superior to two other models in the prediction of helical repeat parameters. We observed a strong linear correlation between curvature magnitudes predicted using the LB model and those determined by electrophoretic ligation ladder experiments, although the LB model somewhat underestimated apparent curvature. With longer electrophoretic phasing probes the LB model slightly overestimated gel mobility anomalies, with modest deviations in predicted helical repeat parameters. Overall, our analyses suggest that the LB model provides reasonably accurate predictions for the electrophoretic behavior of DNA.     

Structural characterization of intrinsically curved AT-rich DNA sequences.

AT-rich DNA sequences other than AnTm tracts (n + m > or = 4) are known to be intrinsically curved. The AATAT-element constitutes one known example of these sequences. In this paper, the elucidation of the structural basis of the curvature induced by this sequence element was addressed. As judged by the patterns of cleavage by the hydroxyl radical and DNase I, the AATAT sequence shows a narrow minor groove. Furthermore, the 5' adenine residue of the AA dinucleotide contained within the sequence is hyperreactive to diethylpyrocarbonate. Similar structural properties are shown by several sequences inducing intrinsic DNA curvature, such as an A5-tract or the closely related ATAAT, AATATA and TAATAT sequences, which are also shown here to induce curvature. On the other hand, other related sequences, such as TATAA and ATATA, that do not induce curvature, show different structural characteristics.     

Sequence-dependent contribution of distal binding domains to CAP protein-DNA binding affinity.

We report measurements of the relative binding affinity of CAP for DNA sequences which have been systematically mutated in the region flanking the consensus binding site. Our experiments focus on the locus one helical turn from the dyad axis where DNA bending toward the minor groove is induced upon C-AP binding. The binding free energy and extent of bending are moderately well correlated for the set of 56 sequences. Changes in binding affinity spanning a factor of about 50 could be accounted for by additive contributions of dinucleotides; with a few exceptions, the relative ranking of dinucleotide contributions to binding and bending are similar. We conclude that dinucleotides are the smallest independent unit required for quantitative interpretation of CAP-induced DNA bending and binding in the distal domains of the CAP consensus binding site. The imperfect correlation between binding strength and extent of bending implies that sequence changes affect protein binding strength not only by altering the DNA deformation energy required to form the complex, but also by affecting directly the free energy of interaction between protein and DNA.     

Dinucleosome DNA of human K562 cells: experimental and computational characterizations

Dinucleosome formation is the first step in the organization of the higher order chromatin structure. With the ultimate aim of elucidating the dinucleosome structure, we constructed a library of human dinucleosome DNA. The library consists of PCR-amplifiable DNA fragments obtained by treatment of nuclei of erythroid K562 cells with micrococcal nuclease followed by extraction of DNA and adaptor ligation to the blunt-ended DNA fragments. The library was then cloned using a plasmid vector and the sequences of the clones were determined. The dominating clones containing the Alu elements were removed. A total of 1002 clones, which comprised a dinucleosome database, contained 84 and 918 clones from the clones before and after removing Alu elements, respectively. Approximately 70% of the clones were between 300 and 400 bp in size and they were distributed to various locations of all chromosomes except the Y chromosome. The clones containing A(2)N(8)A(2)N(8)A(2) or T(2)N(8)T(2)N(8)T(2) sequences were classified into three types, Type I (N shape), Type II (V shape) and Type III (M shape) according to DNA curvature plots. The locations of experimentally determined curved DNA segments matched well with the calculated ones though the clones of Types I and III showed additional curved DNA segments as revealed by the curvature plots. The distributions of complementary dinucleotides in the nucleosome DNA, at the ends of the dinucleosome DNA clones, allowed us to predict the positions of the nucleosome dyad axis, and estimate the size of the nucleosome core DNA, 125nt. The distributions of AA and TT dinucleotides, as well as other RR and YY dinucleotides, showed a periodicity with an average period of 10.4 bases, close to the values observed before. Mapping of nucleosome positions in the dinucleosome database based on the observed periodicity revealed that the nucleosomes were separated by a linker of 7.5+ approximately 10 x n nt. This indicates that the nucleosome-nucleosome orientations are, typically, halfway between parallel and antiparallel. Also an important finding is that the distributions of AA/TT and other RR/YY dinucleotides, apparently, reflect both DNA curvature and DNA bendability, cooperatively contributing to the nucleosome formation.     

Identification of sequence elements contributing to the intrinsic curvature of the mouse satellite DNA repeat.

In this paper, the contribution of different sequence elements to the intrisic curvature of the mouse satellite DNA repeat was investigated. This DNA fragment contains nineteen groups of three or more consecutive adenines which are only poorly phased with respect to the helical repeat. The mouse satellite DNA repeat shows a sinusoidal pattern of cleavage by the hydroxyl radical; the waves of reactivity are phased with respect to the A-tracts. Some interesting observations arise from a detailed analysis of these cleavage patterns: a) the maxima of hydroxyl radical cleavage are more periodically spaced along the DNA sequence than the A-tracts themselves. As a consequence, the position of each maximum with respect to the A-tract is variable; b) the sequence 5' TGGAATATG/AA 3' shows a sinusoidal pattern of hydroxyl radical cleavage. This sequence shows a retarded migration in polyacrylamide gels indicating that it is actually intrinsically curved. These results are discussed in view of the current models for DNA curvature.     

Electron microscopy mapping of pBR322 DNA curvature. Comparison with theoretical models.

A map of local curvature of the pBR322 DNA has been established by electron microscopy analysis of linearized plasmid molecules. To determine their polarity these molecules are one end labelled with an avidin-ferritin-biotin complex and the images are digitized. Local curvature is calculated from two mathematical treatments of the DNA trajectory and expressed in term of a mean dinucleotide wedge angle. Eight regions of curvature are distinguished. The four main regions of curvature have a high content of phased AA runs. The experimental curvature map is compared to theoretical maps of curvature obtained from four available models for DNA curvature.     

In vitro selection of DNAs with an increased propensity to form small circles

A protocol was devised to select for DNA molecules that efficiently form circles from a library of 126 base pair DNAs containing 90 randomized base pairs. After six rounds of selection, individual molecules from the library showed 20‐ to 100‐fold greater j‐factors compared with the starting library, validating the selection protocol. High‐throughput sequencing revealed a sinusoidal pattern of enrichment and de‐enrichment of A/T dinucleotides in the random region with a 10.4 base pair period associated with the helicity of DNA. A similar, but more moderate pattern of C/G dinucleotides was offset by precisely half a helical turn. While C/G dinucleotide enrichments were evenly distributed, A/T dinucleotide enrichments displayed a preference to cluster in individual DNA molecules. The most highly enriched 10 base pair sequences in the random region contained adjacent blocks of A/T and C/G trinucleotides present in some, but not all, rapidly cyclizing molecules. The phased dinucleotide enrichments closely match those present in accurately mapped yeast nucleosomes, confirming the importance of DNA bending in nucleosome formation. However, at certain sites the nucleosomal DNAs show dinucleotide enrichments that differ substantially from the cyclization data. These discrepancies can often be correlated with sequence specific contacts that form between histones and DNA. © 2015 Wiley Periodicals, Inc. Biopolymers 103: 303–320, 2015.     

Sequence structure of human nucleosome DNA

Positional distributions of various dinucleotides in experimentally derived human nucleosome DNA sequences are analyzed. Nucleosome positioning in this species is found to depend largely on GG and CC dinucleotides periodically distributed along the nucleosome DNA sequence, with the period of 10.4 bases. The GG and CC dinucleotides oscillate counterphase, i.e., their respective preferred positions are shifted about a half-period from one another, as it was observed earlier for AA and TT dinucleotides. Other purine-purine and pyrimidine-pyrimidine dinucleotides (RR and YY) display the same periodical and counterphase pattern. The dominance of oscillating GG and CC dinucleotides in human nucleosomes and the contribution of AG(CT), GA(TC), and AA(TT) suggest a general nucleosome DNA sequence pattern - counterphase oscillation of RR and YY dinucleotides. AA and TT dinucleotides, commonly accepted as major players, are only weak contributors in the case of human nucleosomes.     

Repertoires of the Nucleosome-Positioning Dinucleotides

It is generally accepted that the organization of eukaryotic DNA into chromatin is strongly governed by a code inherent in the genomic DNA sequence. This code, as well as other codes, is superposed on the triplets coding for amino acids. The history of the chromatin code started three decades ago with the discovery of the periodic appearance of certain dinucleotides, with AA/TT and RR/YY giving the strongest signals, all with a period of 10.4 bases. Every base-pair stack in the DNA duplex has specific deformation properties, thus favoring DNA bending in a specific direction. The appearance of the corresponding dinucleotide at the distance 10.4 xn bases will facilitate DNA bending in that direction, which corresponds to the minimum energy of DNA folding in the nucleosome. We have analyzed the periodic appearances of all 16 dinucleotides in the genomes of thirteen different eukaryotic organisms. Our data show that a large variety of dinucleotides (if not all) are, apparently, contributing to the nucleosome positioning code. The choice of the periodical dinucleotides differs considerably from one organism to another. Among other 10.4 base periodicities, a strong and very regular 10.4 base signal was observed for CG dinucleotides in the genome of the honey bee A. mellifera. Also, the dinucleotide CG appears as the only periodical component in the human genome. This observation seems especially relevant since CpG methylation is well known to modulate chromatin packing and regularity. Thus, the selection of the dinucleotides contributing to the chromatin code is species specific, and may differ from region to region, depending on the sequence context.     

The DNA sequence specificity of bleomycin cleavage in telomeric sequences in human cells

Bleomycin is an antibiotic drug that is widely used in cancer chemotherapy. Telomeres are located at the ends of chromosomes and comprise the tandemly repeated DNA sequence (GGGTTA) _n in humans. Since bleomycin cleaves DNA at 5??-GT dinucleotide sequences, telomeres are expected to be a major target for bleomycin cleavage. In this work, we determined the DNA sequence specificity of bleomycin cleavage in telomeric sequences in human cells. This was accomplished using a linear amplification procedure, a fluorescently labelled oligonucleotide primer and capillary gel electrophoresis with laser-induced fluorescence detection. This represents the first occasion that the DNA sequence specificity of bleomycin cleavage in telomeric DNA sequences in human cells has been reported. The bleomycin DNA sequence selectivity was mainly at 5??-GT dinucleotides, with lesser amounts at 5??-GG dinucleotides. The cellular bleomycin telomeric DNA damage was also compared with bleomycin telomeric damage in purified human genomic DNA and was found to be very similar. The implications of these results for the understanding of bleomycin??s mechanism of action in human cells are discussed.     

Epigenetic nucleosomes: Alu sequences and CG as nucleosome positioning element

Alu sequences carry periodical pattern with CG dinucleotides (CpG) repeating every 31-32 bases. Similar distances are observed in distribution of DNA curvature in crystallized nucleosomes, at positions +/-1.5 and +/-4.5 periods of DNA from nucleosome DNA dyad. Since CG elements are also found to impart to nucleosomes higher stability when positioned at +/-1.5 sites, it suggests that CG dinucleotides may play a role in modulation of the nucleosome strength when the CG elements are methylated. Thus, Alu sequences may harbor special epigenetic nucleosomes with methylation-dependent regulatory functions. Nucleosome DNA sequence probe is suggested to detect locations of such regulatory nucleosomes in the sequences.     

Evidence for opposite groove-directed curvature of GGGCCC and AAAAA sequence elements.

The repetitive sequence (AGGGCCCTAGAGGGGCCC-TAG)n was previously shown to be curved by gel mobility assays. Here we show, using hydroxy radical/DNase I digestion and differential helical phasing experiments that the curvature is directed towards the major groove and is located in the GGGCCC, but not the CTAGAG segments. The effect of the GC step in the context of the GGGCCC motif is apparently about as large as that of AA/TT, i.e. enough to cancel the macroscopic curvature of helically phased A-tracts. These data are in agreement with positive roll-like curvature of the GCC/GGC motif, predicted from nucleosome packing data and the 3D structure of the GGGGCCCC octamer, but they are not in agreement with the dinucleotide-based roll angle values predicted for AG/CT, TA, GG/CC and GC steps. Our results thus indicate the importance of interactions beyond the dinucleotide steps in predictive models of DNA curvature.     

Sequence-dependent kinks induced in curved DNA

In certain curved DNA fragments without AA dinucleotides, the gel retardation anomaly associated with curvature passes through a maximum with fragment length, indicating length (and electric field) dependent structural transitions in the DNA. We suggest that thermally induced stereochemical kinks in DNA are stabilized in the gel, thus relieving the effects of curvature. These kinks are shown to occur specifically at CA/TG and TA/TA stacks. Other physical and biological evidence points to frequent structural dislocations at CA and TA steps. These reversible sequence dependent kinks may therefore represent a novel class of structural protein-DNA recognition elements.     

Molecular modelling of (A4T4NN)n and (T4A4NN)n: sequence elements responsible for curvature.

The molecular modelling program JUMNA has been used to investigate the origins of the strikingly different curvature of the two sequences, (A4T4NN)n and (T4A4NN)n. Gel electrophoresis and cyclisation studies have shown that only the former of these two sequences is significantly curved. By developing novel superhelical symmetry constraints we were able to study the energetic and structural aspects of polymeric DNA having a controlled curvature. The results obtained (which do not take into account specific hydration effects) correlate well with the experimental data and offer a molecular level explanation of curvature. Although curvature is found to be initiated by specific dinucleotide junctions, deformations spread to surrounding dinucleotide steps and, moreover, sequence effects beyond the dinucleotide level are observed.