The mechanics behind DNA sequence-dependent properties of the nucleosome

Chromatin organization and composition impart sophisticated regulatory features critical to eukaryotic genomic function. Although DNA sequence-dependent histone octamer binding is important for nucleosome activity, many aspects of this phenomenon have remained elusive. We studied nucleosome structure and stability with diverse DNA sequences, including Widom 601 derivatives with the highest known octamer affinities, to establish a simple model behind the mechanics of sequence dependency. This uncovers the unique but unexpected role of TA dinucleotides and a propensity for G|C-rich sequence elements to conform energetically favourably at most locations around the histone octamer, which rationalizes G|C% as the most predictive factor for nucleosome occupancy in vivo. In addition, our findings reveal dominant constraints on double helix conformation by H3-H4 relative to H2A-H2B binding and DNA sequence context-dependency underlying nucleosome structure, positioning and stability. This provides a basis for improved prediction of nucleosomal properties and the design of tailored DNA constructs for chromatin investigations.     

Nucleosome structural studies

Chromatin plays a fundamental role in eukaryotic genomic regulation, and the increasing awareness of the importance of epigenetic processes in human health and disease emphasizes the need for understanding the structure and function of the nucleosome. Recent advances in chromatin structural studies, including the first structures of nucleosomes containing the Widom 601 sequence and the structure of a chromatin protein-nucleosome assembly, have provided new insight into stretching of nucleosomal DNA, nucleosome positioning, binding of metal ions, drugs and therapeutic candidates to nucleosomes, and nucleosome recognition by nuclear proteins. These discoveries ensure promising future prospects for unravelling structural attributes of chromatin.     

Sequence-dependent Nucleosome Positioning

Eukaryotic DNA is organized into a macromolecular structure called chromatin. The basic repeating unit of chromatin is the nucleosome, which consists of two copies of each of the four core histones and DNA. The nucleosomal organization and the positions of nucleosomes have profound effects on all DNA-dependent processes. Understanding the factors that influence nucleosome positioning is therefore of general interest. Among the many determinants of nucleosome positioning, the DNA sequence has been proposed to have a major role. Here, we analyzed more than 860,000 nucleosomal DNA sequences to identify sequence features that guide the formation of nucleosomes in vivo. We found that both a periodic enrichment of AT base pairs and an out-of-phase oscillating enrichment of GC base pairs as well as the overall preference for GC base pairs are determinants of nucleosome positioning. The preference for GC pairs can be related to a lower energetic cost required for deformation of the DNA to wrap around the histones. In line with this idea, we found that only incorporation of both signal components into a sequence model for nucleosome formation results in maximal predictive performance on a genome-wide scale. In this manner, one achieves greater predictive power than published approaches. Our results confirm the hypothesis that the DNA sequence has a major role in nucleosome positioning in vivo.     

Organization of telomeric nucleosomes: atomic force microscopy imaging and theoretical modeling

Telomeric chromatin has peculiar features with respect to bulk chromatin, which are not fully clarified to date. Nucleosomal arrays, reconstituted on fragments of human telomeric DNA and on tandemly repeated tetramers of 5S rDNA, have been investigated at single-molecule level by atomic force microscopy and Monte Carlo simulations. A satisfactory correlation emerges between experimental and theoretical internucleosomal distance distributions. However, in the case of telomeric nucleosomal arrays containing two nucleosomes, we found significant differences. Our results show that sequence features of DNA are significant in the basic chromatin organization, but are not the only determinant.     

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.     

Comparative analysis of the nucleosomal structure of rye,wheat and their relatives

Analysis of the structure of chromatin in cereal species using micrococcal nuclease (MNase) cleavage showed nucleosomal organization and a ladder with typical nucleosomal spacing of 175–185 bp. Probing with a set of DNA probes localized in the authentic telomeres, subtelomeric regions and bulk chromatin revealed that these chromosomal regions have nucleosomal organization but differ in size of nucleosomes and rate of cleavage between both species and regions. Chromatin from Secale and Dasypyrum cleaved more quickly than that from wheat and barley, perhaps because of their higher content of repetitive sequences with hairpin structures accessible to MNase cleavage. In all species, the telomeric chromatin showed more rapid cleavage kinetics and a shorter nucleosome length (160 bp spacing) than bulk chromatin. Rye telomeric repeat arrays were shortest, ranging from 8 kb to 50 kb while those of wheat ranged from 15 kb up to 175 kb. A gradient of sensitivity to MNase was detected along rye chromosomes. The rye-specific subtelomeric sequences pSc200 and pSc250 have nucleosomes of two lengths, those of the telomeric and of bulk nucleosomes, indicating that the telomeric structure may extended into the chromosomes. More proximal sequences common to rye and wheat, the short tandem-repeat pSc119.2 and rDNA sequence pTa71, showed longer nucleosomal sizes characteristic of bulk chromatin in both species. A strictly defined spacing arrangement (phasing) of nucleosomes was demonstrated along arrays of tandem repeats with different monomer lengths (118, 350 and 550 bp) by combining MNase and restriction enzyme digestion.     

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.     

A critical role for linker DNA in higher-order folding of chromatin fibers

Nucleosome-nucleosome interactions drive the folding of nucleosomal arrays into dense chromatin fibers. A better physical account of the folding of chromatin fibers is necessary to understand the role of chromatin in regulating DNA transactions. Here, we studied the unfolding pathway of regular chromatin fibers as a function of single base pair increments in linker length, using both rigid base-pair Monte Carlo simulations and single-molecule force spectroscopy. Both computational and experimental results reveal a periodic variation of the folding energies due to the limited flexibility of the linker DNA. We show that twist is more restrictive for nucleosome stacking than bend, and find the most stable stacking interactions for linker lengths of multiples of 10 bp. We analyzed nucleosomes stacking in both 1- and 2-start topologies and show that stacking preferences are determined by the length of the linker DNA. Moreover, we present evidence that the sequence of the linker DNA also modulates nucleosome stacking and that the effect of the deletion of the H4 tail depends on the linker length. Importantly, these results imply that nucleosome positioning in vivo not only affects the phasing of nucleosomes relative to DNA but also directs the higher-order structure of chromatin.     

Differential association of linker histones H1 and H5 with telomeric nucleosomes in chicken erythrocytes.

Rat liver telomeric DNA is organised into nucleosomes characterised by a shorter and more homogeneous average nucleosomal repeat than bulk chromatin as shown by Makarov et al. (1). The latter authors were unable to detect the association of any linker histone with the telomeric DNA. We have confirmed these observations but show that in sharp contrast chicken erythrocyte telomeric DNA is organised into nucleosomes whose spacing length and heterogeneity are indistinguishable from those of bulk chromatin. We further show that chicken erythrocyte telomeric chromatin contains chromatosomes which are preferentially associated with histone H1 relative to histone H5. This contrasts with bulk chromatin where histone H5 is the more abundant species. This observation strongly suggests that telomeric DNA condensed into nucleosome core particles has a higher affinity for H1 than H5. We discuss the origin of the discrimination of the lysine rich histones in terms of DNA sequence preferences, telomere nucleosome preferences and particular constraints of the higher order chromatin structure of telomeres.     

No overt nucleosome eviction at deprotected telomeres

Dysfunctional telomeres elicit the canonical DNA damage response, which includes the activation of the ATM or ATR kinase signaling pathways and end processing by nonhomologous end joining (NHEJ) or homologous recombination (HR). The cellular response to DNA double-strand breaks has been proposed to involve chromatin remodeling and nucleosome eviction, but whether dysfunctional telomeres undergo chromatin reorganization is not known. Here, we report on the nucleosomal organization of telomeres that have become deprotected through the deletion of the shelterin components TRF2 or POT1. We found no evidence of changes in the nucleosomal organization of the telomeric chromatin or nucleosome eviction near the telomere terminus. An unaltered chromatin structure was observed at telomeres lacking TRF2, which activate the ATM kinase and are a substrate for NHEJ. Similarly, telomeres lacking POT1a and POT1b, which activate the ATR kinase, showed no overt nucleosome eviction. Finally, telomeres lacking TRF2 and Ku70, which are processed by HR, appeared to maintain their original nucleosomal organization. We conclude that ATM signaling, ATR signaling, NHEJ, and HR at deprotected telomeres can take place in the absence of overt nucleosome eviction.     

Spontaneous access to DNA target sites in folded chromatin fibers

DNA wrapped in nucleosomes is sterically occluded from many protein complexes that must act on it; how such complexes gain access to nucleosomal DNA is not known. In vitro studies on isolated nucleosomes show that they undergo spontaneous partial unwrapping conformational transitions, which make the wrapped nucleosomal DNA transiently accessible. Thus, site exposure might provide a general mechanism allowing access of protein complexes to nucleosomal DNA. However, existing quantitative analyses of site exposure focused on single nucleosomes, while the presence of neighbor nucleosomes and concomitant chromatin folding might significantly influence site exposure. In this work, we carried out quantitative studies on the accessibility of nucleosomal DNA in homogeneous nucleosome arrays. Two striking findings emerged. Organization into chromatin fibers changes the accessibility of nucleosomal DNA only modestly, from ∼ 3-fold decreases to ∼ 8-fold increases in accessibility. This means that nucleosome arrays are intrinsically dynamic and accessible even when they are visibly condensed. In contrast, chromatin folding decreases the accessibility of linker DNA by as much as ∼ 50-fold. Thus, nucleosome positioning dramatically influences the accessibility of target sites located inside nucleosomes, while chromatin folding dramatically regulates access to target sites in linker DNA.     

Organisation of nucleosomal arrays reconstituted with repetitive African green monkey α-satellite DNA as analysed by atomic force microscopy

Alpha-satellite DNA (AS) is part of centromeric DNA and could be relevant for centromeric chromatin structure: its repetitive character may generate a specifically ordered nucleosomal arrangement and thereby facilitate kinetochore protein binding and chromatin condensation. Although nucleosomal positioning on some satellite sequences had been shown, including AS from African green monkey (AGM), the sequence-dependent nucleosomal organisation of repetitive AS of this species has so far not been analysed. We therefore studied the positioning of reconstituted nucleosomes on AGM AS tandemly repeated DNA. Enzymatic analysis of nucleosome arrays formed on an AS heptamer as well as the localisation of mononucleosomes on an AS dimer by atomic force microscopy (AFM) showed one major positioning frame, in agreement with earlier results. The occupancy of this site was in the range of 45–50%, in quite good agreement with published in vivo observations. AFM measurements of internucleosomal distances formed on the heptamer indicated that the nucleosomal arrangement is governed by sequence-specific DNA-histone interactions yielding defined internucleosomal distances, which, nevertheless, are not compatible with a uniform phasing of the nucleosomes with the AGM AS repeats.     

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.     

In vitro assembly of the characteristic chromatin organization at the yeast PHO5 promoter by a replication-independent extract system

An extensive set of analyses of the yeast PHO5 gene, mostly performed in vivo, has made this gene a model for the role of chromatin structure in gene regulation. In the repressed state, the PHO5 promoter shows a characteristic chromatin organization with four positioned nucleosomes and a short hypersensitive site. So far the basis for this nucleosome positioning has remained unresolved. We have therefore decided to complement the in vivo studies by an in vitro approach. As a first step, we have asked whether the characteristic PHO5 promoter chromatin structure depends on the cellular context including replication or higher order nuclear chromatin organization or whether it can be reconstituted in vitro in a cell-free system. To this end we have established an in vitro chromatin assembly system based on yeast extracts. It is capable of generating extensive regular nucleosomal arrays with physiological spacing. Assembly requires supplementation with exogenous histones and is dependent on energy leading to chromatin with dynamic properties due to ATP-dependent activities of the extract. Using the PHO5 promoter sequence as template in this replication independent system, we obtain a nucleosomal pattern over the PHO5 promoter region that is very similar to the in vivo pattern of the repressed state. This shows that the chromatin structure at the PHO5 promoter represents a self-organizing system in cell-free yeast extracts and provides a promising substrate for in vitro studies with a direct in vivo correlate.     

Predicting human nucleosome occupancy from primary sequence

Nucleosomes are the fundamental repeating unit of chromatin and comprise the structural building blocks of the living eukaryotic genome. Micrococcal nuclease (MNase) has long been used to delineate nucleosomal organization. Microarray-based nucleosome mapping experiments in yeast chromatin have revealed regularly-spaced translational phasing of nucleosomes. These data have been used to train computational models of sequence-directed nuclesosome positioning, which have identified ubiquitous strong intrinsic nucleosome positioning signals. Here, we successfully apply this approach to nucleosome positioning experiments from human chromatin. The predictions made by the human-trained and yeast-trained models are strongly correlated, suggesting a shared mechanism for sequence-based determination of nucleosome occupancy. In addition, we observed striking complementarity between classifiers trained on experimental data from weakly versus heavily digested MNase samples. In the former case, the resulting model accurately identifies nucleosome-forming sequences; in the latter, the classifier excels at identifying nucleosome-free regions. Using this model we are able to identify several characteristics of nucleosome-forming and nucleosome-disfavoring sequences. First, by combining results from each classifier applied de novo across the human ENCODE regions, the classifier reveals distinct sequence composition and periodicity features of nucleosome-forming and nucleosome-disfavoring sequences. Short runs of dinucleotide repeat appear as a hallmark of nucleosome-disfavoring sequences, while nucleosome-forming sequences contain short periodic runs of GC base pairs. Second, we show that nucleosome phasing is most frequently predicted flanking nucleosome-free regions. The results suggest that the major mechanism of nucleosome positioning in vivo is boundary-event-driven and affirm the classical statistical positioning theory of nucleosome organization.     

TRF2 controls telomeric nucleosome organization in a cell cycle phase-dependent manner

Mammalian telomeres stabilize chromosome ends as a result of their assembly into a peculiar form of chromatin comprising a complex of non-histone proteins named shelterin. TRF2, one of the shelterin components, binds to the duplex part of telomeric DNA and is essential to fold the telomeric chromatin into a protective cap. Although most of the human telomeric DNA is organized into tightly spaced nucleosomes, their role in telomere protection and how they interplay with telomere-specific factors in telomere organization is still unclear. In this study we investigated whether TRF2 can regulate nucleosome assembly at telomeres.By means of chromatin immunoprecipitation (ChIP) and Micrococcal Nuclease (MNase) mapping assay, we found that the density of telomeric nucleosomes in human cells was inversely proportional to the dosage of TRF2 at telomeres. This effect was not observed in the G1 phase of the cell cycle but appeared coincident of late or post-replicative events. Moreover, we showed that TRF2 overexpression altered nucleosome spacing at telomeres increasing internucleosomal distance. By means of an in vitro nucleosome assembly system containing purified histones and remodeling factors, we reproduced the short nucleosome spacing found in telomeric chromatin. Importantly, when in vitro assembly was performed in the presence of purified TRF2, nucleosome spacing on a telomeric DNA template increased, in agreement with in vivo MNase mapping.Our results demonstrate that TRF2 negatively regulates the number of nucleosomes at human telomeres by a cell cycle-dependent mechanism that alters internucleosomal distance. These findings raise the intriguing possibility that telomere protection is mediated, at least in part, by the TRF2-dependent regulation of nucleosome organization.