TA-periodic (“601”-like) centromeric nucleosomes of A. thaliana

The bulk of strong nucleosomes (SNs, with visibly periodic DNA sequences) is described by consensus pattern of 5 or 6 base runs of purines alternating with similar runs of pyrimidines – RR/YY SNs. Yet, the strongest known nucleosome positioning sequence, the 601 clone of Lowary and Widom, is rather periodic repetition of TA dinucleotides following one another every 10 bases. We located “601”-like TA-periodic sequences in the genome of A. thaliana. Several families of such sequences are discovered repeating almost exclusively in centromeres. Thus, while A. thaliana SNs of RR/YY type have strong affinity to pericentromeric regions, as it has been previously found, the SNs of TA periodic type concentrate rather in centromeres.     

Universal full-length nucleosome mapping sequence probe

For the computational sequence-directed mapping of the nucleosomes, the knowledge of the nucleosome positioning motifs – 10–11 base long sequences – and respective matrices of bendability, is not sufficient, since there is no justified way to fuse these motifs in one continuous nucleosome DNA sequence. Discovery of the strong nucleosome (SN) DNA sequences, with visible sequence periodicity allows derivation of the full-length nucleosome DNA bendability pattern as matrix or consensus sequence. The SN sequences of three species (A. thaliana, C. elegans, and H. sapiens) are aligned (512 sequences for each species), and long (115 dinucleotides) matrices of bendability derived for the species. The matrices have strong common property – alternation of runs of purine–purine (RR) and pyrimidine–pyrimidine (YY) dinucleotides, with average period 10.4 bases. On this basis the universal [R,Y] consensus of the nucleosome DNA sequence is derived, with exactly defined positions of respective penta- and hexamers RRRRR, RRRRRR, YYYYY, and YYYYYY.     

Structure-based Analysis of DNA Sequence Patterns Guiding Nucleosome Positioning in vitro

Abstract

Recent studies of genome-wide nucleosomal organization suggest that the DNA sequence is one of the major determinants of nucleosome positioning. Although the search for underlying patterns encoded in nucleosomal DNA has been going on for about 30 years, our knowledge of these patterns still remains limited. Based on our evaluations of DNA deformation energy, we developed new scoring functions to predict nucleosome positioning. There are three principal differences between our approach and earlier studies: (i) we assume that the length of nucleosomal DNA varies from 146 to 147 bp; (ii) we consider the anisotropic flexibility of pyrimidine-purine (YR) dimeric steps in the context of their neighbors (e.g., YYRR versus RYRY); (iii) we postulate that alternating AT-rich and GC-rich motifs reflect sequence-dependent interactions between histone arginines and DNA in the minor groove. Using these functions, we analyzed 20 nucleosome positions mapped in vitro at single nucleotide resolution (including clones 601, 603, 605, the pGUB plasmid, chicken β-globin and three 5S rDNA genes). We predicted 15 of the 20 positions with 1-bp precision, and two positions with 2-bp precision. The predicted position of the ‘601’ nucleosome (i.e., the optimum of the computed score) deviates from the experimentally determined unique position by no more than 1 bp—an accuracy exceeding that of earlier predictions.

Our analysis reveals a clear heterogeneity of the nucleosomal sequences which can be divided into two groups based on the positioning ‘rules’ they follow. The sequences of one group are enriched by highly deformable YR/YYRR motifs at the minor-groove bending sites SHL ±3.5 and ±5.5, which is similar to the α-satellite sequence used in most crystallized nucleosomes. Apparently, the positioning of these nucleosomes is determined by the interactions between histones H2A/H2B and the terminal parts of nucleosomal DNA. In the other group (that includes the ‘601’ clone) the same YR/YYRR motifs occur predominantly at the sites SHL ±1.5. The interaction between the H3/H4 tetramer and the central part of the nucleosomal DNA is likely to be responsible for the positioning of nucleosomes of this group, and the DNA trajectory in these nucleosomes may differ in detail from the published structures.

Thus, from the stereochemical perspective, the in vitro nucleosomes studied here follow either an X-ray-like pattern (with strong deformations in the terminal parts of nucleosomal DNA), or an alternative pattern (with the deformations occurring predominantly in the central part of the nucleosomal DNA). The results presented here may be useful for genome-wide classification of nucleosomes, linking together structural and thermodynamic characteristics of nucleosomes with the underlying DNA sequence patterns guiding their positions.     

Strong nucleosomes of mouse genome including recovered centromeric sequences

Recently discovered strong nucleosomes (SNs) characterized by visibly periodical DNA sequences have been found to concentrate in centromeres of Arabidopsis thaliana and in transient meiotic centromeres of Caenorhabditis elegans. To find out whether such affiliation of SNs to centromeres is a more general phenomenon, we studied SNs of the Mus musculus. The publicly available genome sequences of mouse, as well as of practically all other eukaryotes do not include the centromere regions which are difficult to assemble because of a large amount of repeat sequences in the centromeres and pericentromeric regions. We recovered those missing sequences using the data from MNase-seq experiments in mouse embryonic stem cells, where the sequence of DNA inside nucleosomes, including missing regions, was determined by 100-bp paired-end sequencing. Those nucleosome sequences, which are not matching to the published genome sequence, would largely belong to the centromeres. By evaluating SN densities in centromeres and in non-centromeric regions, we conclude that mouse SNs concentrate in the centromeres of telocentric mouse chromosomes, with ~3.9 times excess compared to their density in the rest of the genome. The remaining non-centromeric SNs are harbored mainly by introns and intergenic regions, by retro-transposons, in particular. The centromeric involvement of the SNs opens new horizons for the chromosome and centromere structure studies.     

CENP-A Nucleosome is a Sensitive Allosteric Scaffold for DNA and Chromatin Factors

Centromeric loci of chromosomes are defined by nucleosomes containing the histone H3 variant CENP-A, which bind their DNA termini more permissively than their canonical counterpart, a feature that is critical for the mitotic fidelity. A recent cryo-EM study demonstrated that the DNA termini of CENP-A nucleosomes, reconstituted with the Widom 601 DNA sequence, are asymmetrically flexible, meaning one terminus is more clearly resolved than the other. However, an earlier work claimed that both ends could be resolved in the presence of two stabilizing single chain variable fragment (scFv) antibodies per nucleosome, and thus are likely permanently bound to the histone octamer. This suggests that the binding of scFv antibodies to the histone octamer surface would be associated with CENP-A nucleosome conformational changes, including stable binding of the DNA termini. Here, we present computational evidence that allows to explain at atomistic level the structural rearrangements of CENP-A nucleosomes resulting from the antibody binding. The antibodies, while they only bind the octamer façades, are capable of altering the dynamics of the nucleosomal core, and indirectly also the surrounding DNA. This effect has more drastic implications for the structure and the dynamics of the CENP-A nucleosome in comparison to its canonical counterpart. Furthermore, we find evidence that the antibodies bind the left and the right octamer façades at different affinities, another manifestation of the DNA sequence. We speculate that the cells could use induction of similar allosteric effects to control centromere function.     

Crystal Structures of Nucleosome Core Particles Containing the ‘601’ Strong Positioning Sequence

Nucleosome positioning plays a key role in genomic regulation by defining histone-DNA context and by modulating access to specific sites. Moreover, the histone-DNA register influences the double-helix structure, which in turn can affect the association of small molecules and protein factors. Analysis of genomic and synthetic DNA has revealed sequence motifs that direct nucleosome positioning in vitro; thus, establishing the basis for the DNA sequence dependence of positioning would shed light on the mechanics of the double helix and its contribution to chromatin structure in vivo. However, acquisition of well-diffracting nucleosome core particle (NCP) crystals is extremely dependent on the DNA fragment used for assembly, and all previous NCP crystal structures have been based on human α-satellite sequences. Here, we describe the crystal structures of Xenopus NCPs containing one of the strongest known histone octamer binding and positioning sequences, the so-called ‘601’ DNA.Two distinct 145-bp 601 crystal forms display the same histone-DNA register, which coincides with the occurrence of DNA stretching-overtwisting in both halves of the particle around five double-helical turns from the nucleosome center, giving the DNA an ‘effective length’ of 147 bp. As we have found previously with stretching around two turns from the nucleosome center for a centromere-based sequence, the terminal stretching observed in the 601 constructs is associated with extreme kinking into the minor groove at purine-purine (pyrimidine-pyrimidine) dinucleotide steps. In other contexts, these step types display an overall nonflexible behavior, which raises the possibility that DNA stretching in the nucleosome or extreme distortions in general have unique sequence dependency characteristics. Our findings indicate that DNA stretching is an intrinsically predisposed site-specific property of the nucleosome and suggest how NCP crystal structures with diverse DNA sequences can be obtained.     

The kinetic landscape of nucleosome assembly: A coarse-grained molecular dynamics study

The organization of nucleosomes along the Eukaryotic genome is maintained over time despite disruptive events such as replication. During this complex process, histones and DNA can form a variety of non-canonical nucleosome conformations, but their precise molecular details and roles during nucleosome assembly remain unclear. In this study, employing coarse-grained molecular dynamics simulations and Markov state modeling, we characterized the complete kinetics of nucleosome assembly. On the nucleosome-positioning 601 DNA sequence, we observe a rich transition network among various canonical and non-canonical tetrasome, hexasome, and nucleosome conformations. A low salt environment makes nucleosomes stable, but the kinetic landscape becomes more rugged, so that the system is more likely to be trapped in off-pathway partially assembled intermediates. Finally, we find that the co-operativity between DNA bending and histone association enables positioning sequence motifs to direct the assembly process, with potential implications for the dynamic organization of nucleosomes on real genomic sequences.     

Sequence Signatures of Nucleosome Positioning in Caenorhabditis elegans

Our recent investigation in the protist Trichomonas vaginalis suggested a DNA sequence periodicity with a unit length of 120.9 nt, which represents a sequence signature for nucleosome positioning. We now extended our observation in higher eukaryotes and identified a similar periodicity of 175 nt in length in Caenorhabditis elegans. In the process of defining the sequence compositional characteristics, we found that the 10.5-nt periodicity, the sequence signature of DNA double helix, may not be sufficient for cross-nucleosome positioning but provides essential guiding rails to facilitate positioning. We further dissected nucleosome-protected sequences and identified a strong positive purine (AG) gradient from the 5′-end to the 3′-end, and also learnt that the nucleosome-enriched regions are GC-rich as compared to the nucleosome-free sequences as purine content is positively correlated with GC content. Sequence characterization allowed us to develop a hidden Markov model (HMM) algorithm for decoding nucleosome positioning computationally, and based on a set of training data from the fifth chromosome of C. elegans, our algorithm predicted 60%-70% of the well-positioned nucleosomes, which is 15%-20% higher than random positioning. We concluded that nucleosomes are not randomly positioned on DNA sequences and yet bind to different genome regions with variable stability, well-positioned nucleosomes leave sequence signatures on DNA, and statistical positioning of nucleosomes across genome can be decoded computationally based on these sequence signatures.     

Nucleosome positioning on yeast plasmids is determined only by the internal signal of DNA sequence occupied by the nucleosome

The possible role of border factors in determining the nucleosome positioning on a DNA sequence was investigated. To this end a family of recombinant plasmids based on Gal10Cyc1 promoter and neomycin phosphotransferase gene NPTII were created. A DNA sequence adjoining the GalCyc promoter was varied in these plasmids. Three nearly equally represented nucleosome positions on the GalCyc promoter were found. In the basal plasmid an FRT sequence adjoins the GalCyc promoter at the right. It contains an internal signal of multiple positioning. Its replacement with different DNA sequences does not affect nucleosome positioning on the GalCyc promoter. The nucleosome positioning on the GalCyc promoter does not depend on nucleosome positioning (or its absence) on adjoining sequences. The same is true for nucleosome positioning on FRT sequence. It was found also that nucleosomes' positioning on the NPTII gene and their mutual disposition, namely the spacing between neighboring nucleosomes (linker length) are determined by the location of positioning signals only. Generally the nucleosome positioning in our experimental model is determined solely by internal DNA sequence occupied by nucleosome. On the other hand, the action of this internal positioning signal does not extend to neighboring DNA sequences.     

Nucleosomes structure and dynamics: effect of CHAPS

Dynamics of nucleosomes and spontaneous unwrapping of DNA are fundamental property of the chromatin enabling access to nucleosomal DNA for regulatory proteins. Probing of such dynamics of nucleosomes performed by single molecule techniques revealed a large scale dynamics of nucleosomes including their spontaneous unwrapping. Dissociation of nucleosomes at low concentrations is a complicating issue for studies with single molecule techniques. In this paper, we tested the ability of 3-[(3-Cholamidopropyl)dimethylammonio]-l-propanesulfonate (CHAPS) to prevent dissociation of nucleosomes. The study was performed with mononucleosome system assembled with human histones H2A, H2B, H3 and H4 on the DNA substrate containing sequence 601 that provides the sequencespecific assembly of nucleosomes. We used Atomic Force Microscopy (AFM) to directly identify nucleosomes and analyze their structure at the nanometer level. These studies showed that in the presence of CHAPS at millimolar concentrations, nucleosomes, even at sub-nanomolar concentrations, remain intact over days compared to a complete dissociation of the same nucleosome sample over 10 min in the absence of CHAPS. Importantly, CHAPS does not change the conformation of nucleosomes as confirmed by the AFM analysis. Moreover, 16 µM CHAPS stabilizes nucleosomes in over one hour incubation in the solution containing as low as 0.4 nM in nucleosomes. The stability of nucleosomes is slightly reduced at physiological conditions (150 mM NaCl), although the nucleosomes dissociate rapidly at 300 mM NaCl. The sequence specificity of the nucleosome in the presence of CHAPS decreased suggesting that the histone core translocates along the DNA substrate utilizing sliding mechanism.     

盐离子作用下组蛋白变体H2A.Z增强核小体结构的稳定性

真核生物染色质的基本结构组成单元是核小体,基因组DNA被压缩在染色质中,核小体的存在通常会抑制转录、复制、修复和重组等发生在DNA模板上的生物学过程。组蛋白变体H2A.Z可以调控染色质结构进而影响基因的转录过程,但其详细的调控机制仍未研究清楚。为了比较含有组蛋白变体H2A.Z的核小体和常规核小体在盐离子作用下的稳定性差异,本文采用Förster共振能量转移的方法检测氯化钠、氯化钾、氯化锰、氯化钙、氯化镁等离子对核小体的解聚影响。实验对Widom 601 DNA序列进行双荧光Cy3和Cy5标记,通过荧光信号值的变化来反映核小体的解聚变化。Förster共振能量转移检测结果显示:在氯化钠、氯化钾、氯化锰、氯化钙和氯化镁作用下,含有组蛋白变体H2A.Z的核小体解聚速度相比于常规核小体要慢,且氯化钙、氯化锰和氯化镁的影响更明显。电泳分析结果表明,在75℃条件下含有组蛋白变体H2A.Z的核小体的解聚速率明显低于常规核小体。采用荧光热漂移检测(fluorescence thermal shift analysis , FTS)进一步分析含有组蛋白变体H2A.Z核小体的稳定性,发现两类核小体的荧光信号均呈现2个明显的增长期,含有组蛋白变体H2A.Z核小体的第1个荧光信号增速期所对应的温度明显高于常规核小体,表明核小体中H2A.Z/H2B二聚体的解聚变性温度要高于常规的H2A/H2B二聚体,含有组蛋白变体H2A.Z核小体的热稳定性高。研究结果均表明,含有组蛋白变体H2A.Z的核小体的结构比常规核小体的结构稳定。     

Biochemical screening of stable dinucleosomes using DNA fragments from a dinucleosome DNA library

The dinucleosome is an informative unit for analysis of the higher-order chromatin structure. DNA fragments forming stable dinucleosomes were screened from a dinucleosome DNA library after the reconstitution of nucleosomes in vitro and digestion with micrococcal nuclease. Reconstituted dinucleosomes showed a diversity of sensitivity to micrococcal nuclease, suggesting that the biochemical stability of a dinucleosome depends, in part, on the DNA fragments. The DNA fragments after the screening were classified into three groups represented by clones bf10, af14 and af32 according to the sensitivity to micrococcal nuclease. Mapping of the nucleosome boundaries by Southern blotting of the DNA after restriction digestion and by primer extension analysis showed that each nucleosome position of clone af32 was fixed. Analysis of reconstituted dinucleosomes using mutant DNA fragments of clone af32 revealed a unique property characteristic of a key nucleosome, given that the replacement of a DNA fragment corresponding to the right nucleosome position resulted in marked sensitivity to micrococcal nuclease, whereas the replacement of the other nucleosome fragment had almost no effect on sensitivity as compared to the original af32 construct. The mutant construct in which the right nucleosome was removed showed multiple nucleosome phases, suggesting that the right nucleosome stabilized first each mononucleosome and then the dinucleosome. An oligonucleotide bending assay revealed that the DNA fragment in the right nucleosome included curved DNA, suggesting that the positioning activity of the nucleosome was attributed to its DNA structure. These results suggest that information for forming stable dinucleosome is embedded in the genomic DNA and that a further characterization of the key nucleosome is useful for understanding the building up of the chromatin structure.     

Strong nucleosomes reside in meiotic centromeres of C. elegans

Recently discovered strong nucleosomes (SNs) are characterized by strongly periodical DNA sequence, with visible rather than hidden sequence periodicity. In a quest for possible functions of the SNs, it has been found that the SNs concentrate within centromere regions of A. thaliana chromosomes . They, however, have been detected in Caenorhabditis elegans as well, although the holocentric chromosomes of this species do not have centromeres. Scrutinizing the SNs of C. elegans and their distributions along the DNA sequences of the chromosomes, we have discovered that the SNs are located mainly at the ends of the chromosomes of C. elegans. This suggests that, perhaps, the ends of the chromosomes fulfill some function(s) of centromeres in this species, as also indicated by the cytogenetic studies on meiotic chromosomes in spermatocytes of C. elegans, where the end-to-end association is observed. The centromeric involvement of the SNs, also found in A. thaliana, opens new horizons for the chromosome and centromere structure studies.     

Activity of FEN1 endonuclease on nucleosome substrates is dependent upon DNA sequence but not flap orientation

We demonstrated previously that human FEN1 endonuclease, an enzyme involved in excising single-stranded DNA flaps that arise during Okazaki fragment processing and base excision repair, cleaves model flap substrates assembled into nucleosomes. Here we explore the effect of flap orientation with respect to the surface of the histone octamer on nucleosome structure and FEN1 activity in vitro. We find that orienting the flap substrate toward the histone octamer does not significantly alter the rotational orientation of two different nucleosome positioning sequences on the surface of the histone octamer but does cause minor perturbation of nucleosome structure. Surprisingly, flaps oriented toward the nucleosome surface are accessible to FEN1 cleavage in nucleosomes containing the Xenopus 5S positioning sequence. In contrast, neither flaps oriented toward nor away from the nucleosome surface are cleaved by the enzyme in nucleosomes containing the high-affinity 601 nucleosome positioning sequence. The data are consistent with a model in which sequence-dependent motility of DNA on the nucleosome is a major determinant of FEN1 activity. The implications of these findings for the activity of FEN1 in vivo are discussed.     

Nucleosome positioning in the <Emphasis Type="Italic">NPTII</Emphasis> neomycin phosphotransferase gene on a yeast plasmid in the repressed and active states

A yeast plasmid was constructed to contain a hybrid GAL-CYC promoter, the NPTII neomycin phosphotransferase gene, and the FRT sequence between them. The CYC part of the GAL-CYC promoter harbored four upstream activating sequences (UASs) and two close TATA boxes. NPTII was efficiently expressed upon induction with galactose, conferring G418 resistance on yeast cells. Nucleosome positioning was studied in repressed and induced NPTII in transformed cells. A stable positioning of three nucleosomes was detected under repressive conditions (growth on glucose). Two nucleosomes were on the CYC part of the promoter, one including both of the TATA boxes. The third nucleosome overlapped the FRT sequence and the start of the NPTII coding region. Each of the three nucleosomes displayed multiple positions, suggesting their sliding along DNA. After induction of NPTII expression with galactose, a sliding of two nucleosomes was detected, exposing the TATA box and a long promoter segment. The 5′-distal nucleosome moved closer to the UASs, bringing them closer to the TATA box, which was assumed to facilitate the assembly of the preinitiation complex. The two nucleosomes slid independently of each other. The second nucleosome moved towards the FRT sequence and repositioned at its nucleosome positioning signal. Galactose-induced expression did not affect the nucleosome positioning in the coding region of NPTII. Unidirectional sliding and repositioning were detected without induction after deacetylase inhibition with trichostatin A. Basal NPTII expression was observed without activation of the GAL-CYC promoter and after a spatial uncoupling of the coding sequence and promoter via gene inversion and was probably driven by the FRT TATA-like element, which is in the region permanently exposed in vivo.     

Prediction of nucleosome positions in the yeast genome based on matched mirror position filtering

Nucleosome positioning can affect the accessibility of the underlying DNA to the nuclear environment and as such plays an essential role in the regulation of cellular processes. Specific patterns have been found in the underlying DNA sequences of the nucleosome, and one of the most important patterns includes dinucleotides distributed every 10 to 11 base pairs. Based on this property, we propose to match each dinucleotide in the sequence against its mirror occurrences for 10 to 11 base pairs on both left-hand and right­hand sides. A large number of matches in a local region will then signify the existence of a nucleosome. In this paper, we propose the matched mirror position filters for efficient matching of periodic dinucleotide patterns and computationally predict the nucleosome positions. Experimental results on the Saccharomyces cerevisiae (yeast) genome show that the proposed algorithm can predict nucleosome positions effectively. More than 50% of our predicted nucleosomes are within 35 base pairs of those detected by biological experiments.     

Sequence Structure of Hidden 10.4-base Repeat in the Nucleosomes of C. elegans

Abstract

By measuring prevailing distances between YY, YR, RR, and RY dinucleotides in the large database of the nucleosome DNA fragments from C. elegans, the consensus sequence structure of the nucleosome DNA repeat of C. elegans was reconstructed: (YYYYYRRRRR)_n. An actual period was estimated to be 10.4 bases. The pattern is fully consistent with the nucleosome DNA patterns of other eukaryotes, as established earlier, and, thus, the YYYYYRRRRR repeat can be considered as consensus nucleosome DNA sequence repeat across eukaryotic species. Similar distance analysis for [A, T] dinucleotides suggested the related pattern (TTTYTARAAA)_n where the TT and AA dinucleotides display rather out of phase behavior, contrary to the “AA or TT” in-phase periodicity, considered in some publications. A weak 5-base periodicity in the distribution of TA dinucleotides was detected.     

Positioning of a Nucleosome on Mouse Satellite DNA Inserted into a Yeast Plasmid Is Determined by Its DNA Sequence and an Adjacent Nucleosome

It has earlier been shown that multiple positioning of nucleosomes on mouse satellite DNA is determined by its nucleotide sequence. To clarify whether other factors, such as boundary ones, can affect the positionings, we modified the environment of satellite DNA monomer by inserting it into a yeast plasmid between inducible GalCyc promoter and a structural region of the yeast FLP gene. We have revealed that the positions of nucleosomes on satellite DNA are identical to those detected upon reconstruction in vitro. The positioning signal (GAAAAA sequence) of satellite DNA governs nucleosome location at the adjacent nucleotide sequence as well. Upon promoter induction the nucleosome, translationally positioned on the GalCyc promoter, transfers to the satellite DNA and its location follows the positioning signal of the latter. Thus, the alternatives of positioning of a nucleosome on satellite DNA are controlled by its nucleotide sequence, though the choice of one of them is determined by the adjacent nucleosome.     

Isw1a does not have strict limitations on the length of extranucleosomal DNAs for mobilization of nucleosomes assembled with HeLa cell histones

The Saccharomyces cerevisiae Isw1a and Isw2 ATP-dependent chromatin-remodeling complexes have important roles in vivo in the regulation of nucleosome positioning and modulation of gene activity. We studied the ability of the Isw1a- and Isw2-remodeling enzymes to reposition nucleosomes in mono- and dinucleosomes templates with variably positioned histone octamers (in the center or at the ends of the DNA fragment). To compare the Isw1a and Isw2 nucleosome-mobilizing activities, we utilized mono- and dinucleosome templates reconstituted with purified HeLa cell histones and DNA containing one or two copies of the “601” nucleosome high-affinity sequence used to specifically position nucleosomes on the DNA. The obtained data suggest that Isw1a is able to mobilize HeLa cell histone-assembled mononucleosomes with long (more than 30?bp) extranucleosomal DNAs protruding from both sides, which contrasts to the previously reported inability of Isw1 to mobilize similar nucleosomes assembled with recombinant yeast histones. The results also suggest that Isw1a and Isw2 can mobilize nucleosomes with unfavorably short linker DNA lengths, and the presence of internucleosomal interactions promotes mobilization of nucleosomes even when the linkers are short.