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.     

The relationship between periodic dinucleotides and the nucleosomal DNA deformation revealed by normal mode analysis

Nucleosomes, which contain DNA and proteins, are the basic unit of eukaryotic chromatins. Polymers such as DNA and proteins are dynamic, and their conformational changes can lead to functional changes. Periodic dinucleotide patterns exist in nucleosomal DNA chains and play an important role in the nucleosome structure. In this paper, we use normal mode analysis to detect significant structural deformations of nucleosomal DNA and investigate the relationship between periodic dinucleotides and DNA motions. We have found that periodic dinucleotides are usually located at the peaks or valleys of DNA and protein motions, revealing that they dominate the nucleosome dynamics. Also, a specific dinucleotide pattern CA/TG appears most frequently.     

Correlation between dinucleotide periodicities and nucleosome positioning on mouse satellite DNA

Previous studies demonstrated 16 well-defined nucleosome locations (A-P) on a tandemly repeated prototype 234 base pair (bp) mouse satellite repeat unit. We have aligned the A-P fragments to search for DNA sequence elements that might contribute to nucleosome placement at these positions. Our results demonstrate a strikingly regular, uninterrupted, periodic pattern for the AA dinucleotide occurrences along the entire length of the aligned fragments. The periodicity of the AA occurrences is about 9.7 bp. The pattern exhibits a local minimum at position 74, near the nucleosome dyad axis of symmetry. Other dinucleotides--including AC: GT, CA: TG, and CC: GG--are also placed periodically, but their patterns of occurrence are less regular and less frequent than AA. The calculated spacings between consecutive preferred nucleosome locations on mouse satellite DNA are nearly identical, corresponding to multiples of 9.7 bp. The correlation between the periodicity of dinucleotide occurrences and the average spacing of nucleosome positions suggests that the preferred nucleosome locations recur at intervals that may correspond to the DNA helical repeat in the mouse satellite nucleosomes, and that the histone octamers sample (or slip along) the duplex in steps of 9.7 bp during nucleosome formation on mouse satellite DNA.     

Nucleosomes Shape DNA Polymorphism and Divergence

An estimated 80% of genomic DNA in eukaryotes is packaged as nucleosomes, which, together with the remaining interstitial linker regions, generate higher order chromatin structures [1]. Nucleosome sequences isolated from diverse organisms exhibit ∼10 bp periodic variations in AA, TT and GC dinucleotide frequencies. These sequence elements generate intrinsically curved DNA and help establish the histone-DNA interface. We investigated an important unanswered question concerning the interplay between chromatin organization and genome evolution: do the DNA sequence preferences inherent to the highly conserved histone core exert detectable natural selection on genomic divergence and polymorphism? To address this hypothesis, we isolated nucleosomal DNA sequences from Drosophila melanogaster embryos and examined the underlying genomic variation within and between species. We found that divergence along the D. melanogaster lineage is periodic across nucleosome regions with base changes following preferred nucleotides, providing new evidence for systematic evolutionary forces in the generation and maintenance of nucleosome-associated dinucleotide periodicities. Further, Single Nucleotide Polymorphism (SNP) frequency spectra show striking periodicities across nucleosomal regions, paralleling divergence patterns. Preferred alleles occur at higher frequencies in natural populations, consistent with a central role for natural selection. These patterns are stronger for nucleosomes in introns than in intergenic regions, suggesting selection is stronger in transcribed regions where nucleosomes undergo more displacement, remodeling and functional modification. In addition, we observe a large-scale (∼180 bp) periodic enrichment of AA/TT dinucleotides associated with nucleosome occupancy, while GC dinucleotide frequency peaks in linker regions. Divergence and polymorphism data also support a role for natural selection in the generation and maintenance of these super-nucleosomal patterns. Our results demonstrate that nucleosome-associated sequence periodicities are under selective pressure, implying that structural interactions between nucleosomes and DNA sequence shape sequence evolution, particularly in introns.     

Sequence motifs and free energies of selected natural and non-natural nucleosome positioning DNA sequences.

Our laboratories recently completed SELEX experiments to isolate DNA sequences that most-strongly favor or disfavor nucleosome formation and positioning, from the entire mouse genome or from even more diverse pools of chemically synthetic random sequence DNA. Here we directly compare these selected natural and non-natural sequences. We find that the strongest natural positioning sequences have affinities for histone binding and nucleosome formation that are sixfold or more lower than those possessed by many of the selected non-natural sequences. We conclude that even the highest-affinity sequence regions of eukaryotic genomes are not evolved for the highest affinity or nucleosome positioning power. Fourier transform calculations on the selected natural sequences reveal a special significance for nucleosome positioning of a motif consisting of approximately 10 bp periodic placement of TA dinucleotide steps. Contributions to histone binding and nucleosome formation from periodic TA steps are more significant than those from other periodic steps such as AA (=TT), CC (=GG) and more important than those from the other YR steps (CA (=TG) and CG), which are reported to have greater conformational flexibility in protein-DNA complexes even than TA. We report the development of improved procedures for measuring the free energies of even stronger positioning sequences that may be isolated in the future, and show that when the favorable free energy of histone-DNA interactions becomes sufficiently large, measurements based on the widely used exchange method become unreliable.     

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.     

Curved DNA

A priori considerations and the concept of the sequence-dependent local curving of the DNA axis. Experimental evidence: electric dichroism (relaxation time measurements); anomalous electrophoretic mobility and gel-filtration of some restriction fragments of DNA; one-sided binding of the nucleosomal DNA to the mica surface. Theoretical predictions concerning the nucleotide sequences of the curved DNA. Discovery of the dinucleotide periodicity in the chromatin DNA. The sequence periodicity as a tool for mapping of the nucleosomes along the sequences. Preferential binding of the histone octamers to the curved pieces of DNA--sequence analysis predictions and comparison with experiments: Theoretical and experimental estimates of the tilt and roll angles for different combinations of the neighboring base-pairs. Inherent sequence-dependent curvature and apparent persistence length of DNA.     

The Molecular Evolution of Nucleosome Positioning Through Sequence-Dependent Deformation of the DNA Polymer

Abstract

The computational prediction of nucleosome positioning from DNA sequence now allows for in silico investigation of the molecular evolution of biophysical properties of the DNA molecule responsible for primary chromatin organization in the genome. To discern what signal components driving nucleosome positioning in the yeast genome are potentially targeted by natural selection, we compare the performance of various models predictive of nucleosome positioning within the context of a simple statistical test, the repositioned mutation test. We demonstrate that while nucleosome occupancy is driven largely by translational exclusion in response to AT content, there is also a strong signature of evolutionary conservation of regular patterns within nucleosomal DNA sequence related to the structural organization of the nucleosome core (e.g., 10-bp dinucleotide periodicity). We also use computer simulations to investigate hypothetical coding and regulatory constraints on the ability of sequence properties affecting nucleosome formation to adaptively evolve. Our results demonstrate that natural selection may act independently on different DNA sequence properties responsible for local chromatin organization. Furthermore, at least with respect to the deformation energy of the DNA molecule in the nucleosome, the presence of the genetic code has greatly restricted the ability of sequences to evolve the dynamic nucleosome organization typically observed in promoter regions.     

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.     

Cheminformatics methods for novel nanopore analysis of HIV DNA termini

BACKGROUND: Channel current feature extraction methods, using Hidden Markov Models (HMMs) have been designed for tracking individual-molecule conformational changes. This information is derived from observation of changes in ionic channel current blockade "signal" upon that molecule's interaction with (and occlusion of) a single nanometer-scale channel in a "nanopore detector". In effect, a nanopore detector transduces single molecule events into channel current blockades. HMM analysis tools described are used to help systematically explore DNA dinucleotide flexibility, with particular focus on HIV's highly conserved (and highly flexible/reactive) viral DNA termini. One of the most critical stages in HIV's attack is the binding between viral DNA and the retroviral integrase, which is influenced by the dynamic-coupling induced high flexibility of a CA/TG dinucleotide positioned precisely two base-pairs from the blunt terminus of the duplex viral DNA. This suggests the study of a family of such CA/TG dinucleotide molecules via nanopore measurement and cheminformatics analysis. RESULTS: HMMs are used for level identification on the current blockades, HMM/EM with boosted variance emissions are used for level projection pre-processing, and time-domain FSAs are used to parse the level-projected waveform for kinetic information. The observed state kinetics of the DNA hairpins containing the CA/TG dinucleotide provides clear evidence for HIV's selection of a peculiarly flexible/interactive DNA terminus.     

Modified curved DNA that could allow local DNA underwinding at the nucleosomal pseudodyad fails to position a nucleosome in vivo.

In competitive in vitro reconstitution experiments synthetic DNA composed of tandem repeats of the repetitive sequence (A/T)3NN(G/C)3NN, specifically the 20 bp 'TG sequence' (5'-TCGGTGTTAGAGCCTGTAAC-3'), was reported to associate with the histone octamer with an affinity higher than that of nucleosomally derived DNA. However, at least two groups have independently shown that tandem repeats of the TG sequence do not accommodate a stably positioned nucleosome in vivo. It was suggested that the anisotropic flexibility of the TG sequence, governed by a 10 bp sequence periodicity, is incompatible with the required underwinding of the DNA helix at the nucleosome pseudodyad while maintaining a bending preference that can be accommodated in the remainder of the nucleosome. Here we test this hypothesis directly by studying the in vivo nucleosomal structure of modified TG sequences designed to accommodate underwinding at the pseudodyad. We show that these modifications are not sufficient to allow stable incorporation of the TG sequence repeat into a nucleosome in vivo, but do note invasion from one end of the TG heptamer of a translationally random but rotationally constrained nucleosome. We discuss possible reasons for the absence of nucleosomes from the TG sequence in vivo.     

Oligonucleotide Sequence Motifs as Nucleosome Positioning Signals

To gain a better understanding of the sequence patterns that characterize positioned nucleosomes, we first performed an analysis of the periodicities of the 256 tetranucleotides in a yeast genome-wide library of nucleosomal DNA sequences that was prepared by in vitro reconstitution. The approach entailed the identification and analysis of 24 unique tetranucleotides that were defined by 8 consensus sequences. These consensus sequences were shown to be responsible for most if not all of the tetranucleotide and dinucleotide periodicities displayed by the entire library, demonstrating that the periodicities of dinucleotides that characterize the yeast genome are, in actuality, due primarily to the 8 consensus sequences. A novel combination of experimental and bioinformatic approaches was then used to show that these tetranucleotides are important for preferred formation of nucleosomes at specific sites along DNA in vitro. These results were then compared to tetranucleotide patterns in genome-wide in vivo libraries from yeast and C. elegans in order to assess the contributions of DNA sequence in the control of nucleosome residency in the cell. These comparisons revealed striking similarities in the tetranucleotide occurrence profiles that are likely to be involved in nucleosome positioning in both in vitro and in vivo libraries, suggesting that DNA sequence is an important factor in the control of nucleosome placement in vivo. However, the strengths of the tetranucleotide periodicities were 3–4 fold higher in the in vitro as compared to the in vivo libraries, which implies that DNA sequence plays less of a role in dictating nucleosome positions in vivo. The results of this study have important implications for models of sequence-dependent positioning since they suggest that a defined subset of tetranucleotides is involved in preferred nucleosome occupancy and that these tetranucleotides are the major source of the dinucleotide periodicities that are characteristic of positioned nucleosomes.     

A nonlinear dynamic model of DNA with a sequence-dependent stacking term

No simple model exists that accurately describes the melting behavior and breathing dynamics of double-stranded DNA as a function of nucleotide sequence. This is especially true for homogenous and periodic DNA sequences, which exhibit large deviations in melting temperature from predictions made by additive thermodynamic contributions. Currently, no method exists for analysis of the DNA breathing dynamics of repeats and of highly G/C- or A/T-rich regions, even though such sequences are widespread in vertebrate genomes. Here, we extend the nonlinear Peyrard–Bishop–Dauxois (PBD) model of DNA to include a sequence-dependent stacking term, resulting in a model that can accurately describe the melting behavior of homogenous and periodic sequences. We collect melting data for several DNA oligos, and apply Monte Carlo simulations to establish force constants for the 10 dinucleotide steps (CG, CA, GC, AT, AG, AA, AC, TA, GG, TC). The experiments and numerical simulations confirm that the GG/CC dinucleotide stacking is remarkably unstable, compared with the stacking in GC/CG and CG/GC dinucleotide steps. The extended PBD model will facilitate thermodynamic and dynamic simulations of important genomic regions such as CpG islands and disease-related repeats.     

Yeast nucleosome DNA pattern: deconvolution from genome sequences of S. cerevisiae

Positional correlation analysis for the complete genome of Saccharomyces cerevisiae is performed with the aim to reveal possible chromatin-related sequence features. A strong periodicity with the period 10.4 bases is detected in the distance histograms for the dinucleotides AA and TT, with the characteristic decay distance of approximately 50 base pairs. The oscillations are observed as well in the distributions of other dinucleotides. However, the respective amplitudes are small, consistent with secondary effects, due to dominant periodicity of AA and TT. The observations are in accord with earlier data on the chromatin sequence periodicities and nucleosome DNA sequence patterns. The autocorrelations of AA and TT dinucleotides in yeast include also a counter-phase component. A tentative DNA sequence pattern for the yeast nucleosomes is suggested and verified by comparison of its autocorrelation plots with the respective natural autocorrelations. The nucleosome mapping guided by the pattern is in accord with experimental data on the linker length distribution in yeast.     

Propeller-Twisting of Base-pairs and the Conformational Mobility of Dinucleotide Steps in DNA

When DNA is bent around a protein, it must distort. The distortion occurs by changes in the conformation of successive dinucleotide steps. Bending does not necessarily occur uniformly: some steps might remain particularly rigid, i.e. they might deform relatively little, while others might take more than their proportional share of deformation. We investigate here the deformational capacity of specific dinucleotide steps by examining a database of crystallized oligomers. Dividing the steps into ten types by sequence (AA(=TT), AC(=GT), AG(=CT), AT, CA(=TG), CG, GA(=TC), GC, GG(=CC) and TA), we find that some step types are practically rigid, while others have considerable internal mobility or conformational flexibility. Now in general base-pairs are not planar, but have Propeller-Twist. We find a clear empirical correlation between the level of Propeller-Twist in the base-pairs and the flexibility of the dinucleotide step which they constitute. Propeller-Twist in the base-pairs makes stacking into a dinucleotide step more awkward than in plane base-pairs. In particular, it provides a stereochemical “locking” effect which can make steps with highly Propeller-Twisted base-pairs rigid. Although the origins of Propeller-Twist are not yet clearly understood, this result provides a key to understanding the flexibility of DNA in bending around proteins.     

Yeast Nucleosome DNA Pattern: Deconvolution from Genome Sequences of S. cerevisiae

Abstract

Positional correlation analysis for the complete genome of Saccharomyces cerevisiae is performed with the aim to reveal possible chromatin-related sequence features. A strong periodicity with the period 10.4 bases is detected in the distance histograms for the dinucleotides AA and TT, with the characteristic decay distance of approximately 50 base pairs. The oscillations are observed as well in the distributions of other dinucleotides. However, the respective amplitudes are small, consistent with secondary effects, due to dominant periodicity of AA and TT. The observations are in accord with earlier data on the chromatin sequence periodicities and nucleosome DNA sequence patterns. The autocorrelations of AA and TT dinucleotides in yeast include also a counter-phase component. A tentative DNA sequence pattern for the yeast nucleosomes is suggested and verified by comparison of its autocorrelation plots with the respective natural autocorrelations. The nucleosome mapping guided by the pattern is in accord with experimental data on the linker length distribution in yeast.     

The conserved CA/TG motif at Mu termini: T specifies stable transpososome assembly

The dinucleotide CA/TG found at the termini of transposable phage Mu occurs also at the termini of a large class of transposable elements, including HIV, all retroviruses and many retrotransposons. It was shown recently that mutations of this sequence block transpososome assembly, that A/T is more critical for activity than C/G, and that the hierarchy of reactivity of mutant termini follows closely the reported hierarchy of flexibility of their dinucleotide steps. In order to test the hypothesis that the terminal dinucleotide plays an essential structural role during "open termini" formation accompanying assembly, we have examined the activity of substrates carrying 100 different pairs of mismatched termini. Consistent with the flexibility hypothesis, we find that mismatched substrates are extremely efficient at assembly. A wild-type T residue on the bottom strand is essential for stable assembly, but the identity of the dinucleotide on the top strand is irrelevant for transposition chemistry. In addition, we have found a new rule for suppression of terminal defects by MuB protein, as well as a role for metal ions in DNA opening at the termini.