The main role of the sequence-dependent DNA elasticity in determining the free energy of nucleosome formation on telomeric DNAs

Using a competitive reconstitution assay, we measured the free energy spent in nucleosome formation of eight telomeric DNAs, differing in sequence and/or in length. The obtained values are in satisfactorily good agreement with those derived from a theoretical model that allows the calculation of the free energy of nucleosome formation on the basis of sequence-dependent DNA elasticity, using a statistical thermodynamic approach. Both theoretical and experimental evaluations show that telomeres are characterized by the highest free energies of nucleosome formation among all the DNA sequences so far studied. The free energy of nucleosome formation varies according to the different telomeric sequences and the length of the fragments. Theoretical analysis and experimental mapping by lambda exonuclease show that telomeric nucleosomes occupy multiple positions spaced every telomeric repeat. Sequence-dependent DNA elasticity appears as the main determinant of the stability of telomeric nucleosomes and their multiple translational positioning.     

A theoretical model for the prediction of sequence-dependent nucleosome thermodynamic stability

A theoretical model for predicting nucleosome thermodynamic stability in terms of DNA sequence is advanced. The model is based on a statistical mechanical approach, which allows the calculation of the canonical ensemble free energy involved in the competitive nucleosome reconstitution. It is based on the hypothesis that nucleosome stability mainly depends on the bending and twisting elastic energy to transform the DNA intrinsic superstructure into the nucleosomal structure. The ensemble average free energy is calculated starting from the intrinsic curvature, obtained by integrating the dinucleotide step deviations from the canonical B-DNA and expressed in terms of a Fourier series, in the framework of first-order elasticity. The sequence-dependent DNA flexibility is evaluated from the differential double helix thermodynamic stability. A large number of free-energy experimental data, obtained in different laboratories by competitive nucleosome reconstitution assays, are successfully compared to the theoretical results. They support the hypothesis that the stacking energies are the major factor in DNA rigidity and could be a measure of DNA stiffness. A dual role of DNA intrinsic curvature and flexibility emerges in the determination of nucleosome stability. The difference between the experimental and theoretical (elastic) nucleosome-reconstitution free energy for the whole pool of investigated DNAs suggests a significant role for the curvature-dependent DNA hydration and counterion interactions, which appear to destabilize nucleosomes in highly curved DNAs. This model represents an attempt to clarify the main features of the nucleosome thermodynamic stability in terms of physical-chemical parameters and suggests that in molecular systems with a large degree of complexity, the average molecular properties dominate over the local features, as in a statistical ensemble.     

AFM imaging and theoretical modeling studies of sequence-dependent nucleosome positioning

Telomeric chromatin has different features with respect to bulk chromatin, since nucleosomal repeat along the chain is unusually short. We studied the role of telomeric DNA sequences on nucleosomal spacing in a model system. Nucleosomal arrays, assembled on a 1500-bp-long human telomeric DNA and on a DNA fragment containing 8 copies of the 601 strong nucleosome positioning sequence, have been studied at the single molecule level, by atomic force microscopy imaging. Random nucleosome positioning was found in the case of human telomeric DNA. On the contrary, nucleosome positioning on 601 DNA is characterized by preferential positions of nucleosome dyad axis each 200 bp. The AFM-derived nucleosome organization is in satisfactory agreement with that predicted by theoretical modeling, based on sequence-dependent DNA curvature and flexibility. The reported results show that DNA sequence has a main role, not only in mononucleosome thermodynamic stability, but also in the organization of nucleosomal arrays.     

Dual role of DNA intrinsic curvature and flexibility in determining nucleosome stability

A statistical mechanistic approach to evaluate the sequence-dependent thermodynamic stability of nucleosomes is proposed. The model is based on the calculation of the DNA intrinsic curvature, obtained by integrating the nucleotide step deviations from the canonical B-DNA structure, and on the evaluation of the first order elastic distortion energy to reach the nucleosomal superstructure. Literature data on the free energy of nucleosome formation as obtained by competitive nucleosome reconstitution of a significant pool of different DNA sequences were compared with the theoretical results, and a satisfactorily good correlation was found. A striking result of the comparison is the emergence of two opposite roles of the DNA intrinsic curvature and flexibility in determining nucleosome stability. Finally, the obtained results suggest that the curvature-dependent DNA hydration should play a relevant role in the sequence-dependent nucleosome stability.     

From the sequence to the superstructural properties of DNAs

A theoretical model for predicting intrinsic and induced DNA superstructures as well as their thermodynamic properties is presented. Intrinsic sequence-dependent superstructures are evaluated by integrating local deviations from the canonical B-DNA of the different dinucleotide steps. Induced superstructures are obtained by adopting the principle of minimum deformation free energy, evaluated in the Fourier space, in the framework of first-order elasticity. Finally dinucleotide stacking energies and melting temperatures are considered to account for local flexibility. In fact the two scales are strongly correlated. The model works very satisfactorily in predicting the sequence-dependent effects on the DNA experimental behavior, such as the gel electrophoresis retardation, the writhe transitions in topologically constrained domains, the thermodynamic constants of circularization reactions as well as the nucleosome thermodynamic stability constants.     

Acetylated nucleosome assembly on telomeric DNAs

The role of histone N-terminal domains on the thermodynamic stability of nucleosomes assembled on several different telomeric DNAs as well as on 'average' sequence DNA and on strong nucleosome positioning sequences, has been studied by competitive reconstitution. We find that histone tails hyperacetylation favors nucleosome formation, in a similar extent for all the examined sequences. On the contrary, removal of histone terminal domains by selective trypsinization causes a decrease of nucleosome stability which is smaller for telomeres compared to the other sequences examined, suggesting that telomeric sequences have only minor interactions with histone tails. Micrococcal nuclease kinetics shows enhanced accessibility of acetylated nucleosomes formed both on telomeric and 'average' sequence DNAs. These results suggest a more complex role for histone acetylation than the decrease of electrostatic interactions between DNA and histones.     

Telomeric nucleosomes are intrinsically mobile

Nucleosomes are no longer considered only static basic units that package eukaryotic DNA but they emerge as dynamic players in all chromosomal processes. Regulatory proteins can gain access to recognition sequences hidden by the histone octamer through the action of ATP-dependent chromatin remodeling complexes that cause nucleosome sliding. In addition, it is known that nucleosomes are able to spontaneously reposition along the DNA due to intrinsic dynamic properties, but it is not clear yet to what extent sequence-dependent dynamic properties contribute to nucleosome repositioning. Here, we study mobility of nucleosomes formed on telomeric sequences as a function of temperature and ionic strength. We find that telomeric nucleosomes are highly intrinsically mobile under physiological conditions, whereas nucleosomes formed on an average DNA sequence mostly remain in the initial position. This indicates that DNA sequence affects not only the thermodynamic stability and the positioning of nucleosomes but also their dynamic properties. Moreover, our findings suggest that the high mobility of telomeric nucleosomes may be relevant to the dynamics of telomeric chromatin.     

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.     

A statistical thermodynamic model applied to experimental AFM population and location data is able to quantify DNA-histone binding strength and internucleosomal interaction differences between acetylated and unacetylated nucleosomal arrays

Imaging of nucleosomal arrays by atomic force microscopy allows a determination of the exact statistical distributions for the numbers of nucleosomes per array and the locations of nucleosomes on the arrays. This precision makes such data an excellent reference for testing models of nucleosome occupation on multisite DNA templates. The approach presented here uses a simple statistical thermodynamic model to calculate theoretical population and positional distributions and compares them to experimental distributions previously determined for 5S rDNA nucleosomal arrays (208-12,172-12). The model considers the possible locations of nucleosomes on the template, and takes as principal parameters an average free energy of interaction between histone octamers and DNA, and an average wrapping length of DNA around the octamers. Analysis of positional statistics shows that it is possible to consider interactions between nucleosomes and positioning effects as perturbations on a random positioning noninteracting model. Analysis of the population statistics is used to determine histone-DNA association constants and to test for differences in the free energies of nucleosome formation with different types of histone octamers, namely acetylated or unacetylated, and different DNA templates, namely 172-12 or 208-12 5S rDNA multisite templates. The results show that the two template DNAs bind histones with similar affinities but histone acetylation weakens the association of histones with both templates. Analysis of locational statistics is used to determine the strength of specific nucleosome positioning tendencies by the DNA templates, and the strength of the interactions between neighboring nucleosomes. The results show only weak positioning tendencies and that unacetylated nucleosomes interact much more strongly with one another than acetylated nucleosomes; in fact acetylation appears to induce a small anticooperative occupation effect between neighboring nucleosomes.     

Reconstitution experiments show that sequence-specific histone-DNA interactions are the basis for nucleosome phasing on mouse satellite DNA

The molecular basis underlying the sequence-specific positioning of nucleosomes on DNA was investigated. We previously showed that histone octamers occupy multiple specific positions on mouse satellite DNA in vivo and have now reconstituted the 234 bp mouse satellite repeat unit with pure core histones into mononucleosomes. Histones from mouse liver or chicken erythrocytes bind to the DNA in multiple precisely defined frames in perfect phase with a diverged 9 bp subrepeat of the satellite DNA. This is the first time that nucleosome positions on a DNA in vivo have been compared to those found on the same DNA by in vitro reconstitution. Most of the nucleosomes occupy identical positions in vivo and in vitro. There are, however, some characteristic differences. We conclude that sequence-dependent histone-DNA interactions play a decisive role in the positioning of nucleosomes in vivo, but that the nucleosome locations in native chromatin are subject to additional constraints.     

The core histone tail domains contribute to sequence-dependent nucleosome positioning

The precise positioning of nucleosomes plays a critical role in the regulation of gene expression by modulating the DNA binding activity of trans-acting factors. However, molecular determinants responsible for positioning are not well understood. We examined whether the removal of the core histone tail domains from nucleosomes reconstituted with specific DNA fragments led to alteration of translational positions. Remarkably, we find that removal of tail domains from a nucleosome assembled on a DNA fragment containing a Xenopus borealis somatic-type 5S RNA gene results in repositioning of nucleosomes along the DNA, including two related major translational positions that move about 20 bp further upstream with respect to the 5S gene. In a nucleosome reconstituted with a DNA fragment containing the promoter of a Drosophila alcohol dehydrogenase gene, several translational positions shifted by about 10 bp along the DNA upon tail removal. However, the positions of nucleosomes assembled with a DNA fragment known to have one of the highest binding affinities for core histone proteins in the mouse genome were not altered by removal of core histone tail domains. Our data support the notion that the basic tail domains bind to nucleosomal DNA and influence the selection of the translational position of nucleosomes and that once tails are removed movement between translational positions occurs in a facile manner on some sequences. However, the effect of the N-terminal tails on the positioning and movement of a nucleosome appears to be dependent on the DNA sequence such that the contribution of the tails can be masked by very high affinity DNA sequences. Our results suggest a mechanism whereby sequence-dependent nucleosome positioning can be specifically altered by regulated changes in histone tail-DNA interactions in chromatin.     

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.     

Biophysical characterization of the centromere-specific nucleosome from budding yeast

The centromeric DNA of all eukaryotes is assembled upon a specialized nucleosome containing a histone H3 variant known as CenH3. Despite the importance and conserved nature of this protein, the characteristics of the centromeric nucleosome are still poorly understood. In particular, the stoichiometry and DNA-binding properties of the CenH3 nucleosome have been the subject of some debate. We have characterized the budding yeast centromeric nucleosome by biochemical and biophysical methods and show that it forms a stable octamer containing two copies of the Cse4 protein and wraps DNA in a left-handed supercoil, similar to the canonical H3 nucleosome. The DNA-binding properties of the recombinant nucleosome are identical to those observed in vivo demonstrating that the octameric structure is physiologically relevant.     

Measurement of histone-DNA interaction free energy in nucleosomes

Nucleosome positioning DNA sequences are of increasing interest because of their proposed roles in gene regulation and other chromosome functions in vivo, and because they have revealed new insights into the sequence-dependent structures and mechanics of DNA itself. Here, we describe methods to quantify the relative affinities of histone-DNA interactions in nucleosomes, i.e., the nucleosome positioning power of differing DNA sequences. We review methods developed by others and then discuss in detail our own approach to measurement of histone-DNA interaction free energies. Compared to earlier methods, our dialysis-based approach reduces the possibility that non-equilibrium or irreproducible results could be obtained. It facilitates a direct comparison of free energies for many sequences at the same time and it allows analysis of DNAs having a wide range of relative affinities.     

Effect of DNA length and H4 acetylation on the thermal stability of reconstituted nucleosome particles

To examine the factors involved with nucleosome stability, we reconstituted nonacetylated particles containing various lengths (192, 162, and 152 base pairs) of DNA onto the Lytechinus variegatus nucleosome positioning sequence in the absence of linker histone. We characterized the particles and examined their thermal stability. DNA of less than chromatosome length (168 base pairs) produces particles with altered denaturation profiles, possibly caused by histone rearrangement in those core-like particles. We also examined the effects of tetra-acetylation of histone H4 on the thermal stability of reconstituted nucleosome particles. Tetra-acetylation of H4 reduces the nucleosome thermal stability by 0.8 degrees C as compared with nonacetylated particles. This difference is close to values published comparing bulk nonacetylated nucleosomes and core particles to ones enriched for core histone acetylation, suggesting that H4 acetylation has a dominant effect on nucleosome particle energetics.     

DNA structure, nucleosome placement and chromatin remodelling: a perspective

A major question in chromatin biology is to what extent the sequence of DNA directly determines the genetic and chromatin organization of a eukaryotic genome? We consider two aspects to this question: the DNA sequence-specified positioning of nucleosomes and the determination of NDRs (nucleosome-depleted regions) or barriers. We argue that, in budding yeast, while DNA sequence-specified nucleosome positioning may contribute to positions flanking the regions lacking nucleosomes, DNA thermodynamic stability is a major component determinant of the genetic organization of this organism.     

A Comparison of In Vitro Nucleosome Positioning Mapped with Chicken,Frog and a Variety of Yeast Core Histones

Using high-throughput sequencing, we have mapped sequence-directed nucleosome positioning in vitro on four plasmid DNAs containing DNA fragments derived from the genomes of sheep, drosophila, human and yeast. Chromatins were prepared by reconstitution using chicken, frog and yeast core histones. We also assembled yeast chromatin in which histone H3 was replaced by the centromere-specific histone variant, Cse4. The positions occupied by recombinant frog and native chicken histones were found to be very similar. In contrast, nucleosomes containing the canonical yeast octamer or, in particular, the Cse4 octamer were assembled at distinct populations of locations, a property that was more apparent on particular genomic DNA fragments. The factors that may contribute to this variation in nucleosome positioning and the implications of the behavior are discussed.     

Sequence-dependent nucleosome structural and dynamic polymorphism. Potential involvement of histone H2B N-terminal tail proximal domain

Relaxation of nucleosomes on an homologous series (pBR) of ca 350-370 bp DNA minicircles originating from plasmid pBR322 was recently used as a tool to study their structure and dynamics. These nucleosomes thermally fluctuated between three distinct DNA conformations within a histone N-terminal tail-modulated equilibrium: one conformation was canonical, with 1.75 turn wrapping and negatively crossed entering and exiting DNAs; another was also "closed", but with these DNAs positively crossed; and the third was "open", with a lower than 1.5 turn wrapping and uncrossed DNAs. In this work, a new minicircle series (5S) of similar size was used, which contained the 5S nucleosome positioning sequence. Results showed that DNA in pBR nucleosomes was untwisted by approximately 0.2 turn relative to 5S nucleosomes, which DNase I footprinting confirmed in revealing a approximately 1 bp untwisting at each of the two dyad-distal sites where H2B N-terminal tails pass between the two gyres. In contrast, both nucleosomes showed untwistings at the dyad-proximal sites, i.e. on the other gyre, which were also observed in the high-resolution crystal structure. 5S nucleosomes also differ with respect to their dynamics: they hardly accessed the positively crossed conformation, but had an easier access to the negatively crossed conformation. Simulation showed that such reverse effects on the conformational free energies could be simply achieved by slightly altering the trajectories of entering and exiting DNAs. We propose that this is accomplished by H2B tail untwisting at the distal sites through action at a distance ( approximately 20 bp) on H3-tail interactions with the small groove at the nucleosome entry-exit. These results may help to gain a first glimpse into the two perhaps most intriguing features of the high-resolution structure: the alignment of the grooves on the two gyres and the passage of H2B and H3 N-terminal tails between them.     

Nucleosomal locations of dominant DNA sequence motifs for histone-DNA interactions and nucleosome positioning

DNA sequence is an important determinant of the positioning, stability, and activity of nucleosomes, yet the molecular basis of these effects remains elusive. A "consensus DNA sequence" for nucleosome positioning has not been reported and, while certain DNA sequence preferences or motifs for nucleosome positioning have been discovered, how they function is not known. Here, we report that an unexpected observation concerning the reassembly of nucleosomes during salt gradient dialysis has allowed a breakthrough in our efforts to identify the nucleosomal locations of the DNA sequence motifs that dominate histone-DNA interactions and nucleosome positioning. We conclude that a previous selection experiment for high-affinity, nucleosome-forming DNA sequences exerted selective pressure chiefly on the central stretch of the nucleosomal DNA. This observation implies that algorithms for aligning the selected DNA sequences should seek to optimize the alignment over much less than the full 147 bp of nucleosomal DNA. A new alignment calculation implemented these ideas and successfully aligned 19 of the 41 sequences in a non-redundant database of selected high-affinity, nucleosome-positioning sequences. The resulting alignment reveals strong conservation of several stretches within a central 71 bp of the nucleosomal DNA. The alignment further reveals an inherent palindromic symmetry in the selected DNAs; it makes testable predictions of nucleosome positioning on the aligned sequences and for the creation of new positioning sequences, both of which are upheld experimentally; and it suggests new signals that may be important in translational nucleosome positioning.     

Unique translational positioning of nucleosomes on synthetic DNAs.

A computational study was previously carried out to analyze DNA sequences that are known to position histone octamers at single translational sites. A conserved pattern of intrinsic DNA curvature was uncovered that was proposed to direct the formation of nucleosomes to unique positions. The pattern consists of two regions of curved DNA separated by preferred lengths of non-curved DNA. In the present study, 11 synthetic DNAs were constructed which contain two regions of curved DNA of the form [(A5.T5)(G/C)5]4 separated by non-curved regions of variable length. Translational mapping experiments of in vitro reconstituted mononucleosomes using exonuclease III, micrococcal nuclease and restriction enzymes demonstrated that two of the fragments positioned nucleosomes at a single site while the remaining fragments positioned octamers at multiple sites spaced at 10 base intervals. The synthetic molecules that positioned nucleosomes at a single site contain non-curved central regions of the same lengths that were seen in natural nucleosome positioning sequences. Hydroxyl radical and DNase I digests of the synthetic DNAs in reconstituted nucleosomes showed that the synthetic curved element on one side of the nucleosomal dyad assumed a rotational orientation where narrow minor grooves of the A-tracts faced the histone surface with all molecules. In contrast, the curved element on the other side of the nucleosome displayed variable rotational orientations between molecules which appeared to be related to the positioning effect. These results suggest that asymmetry between the two halves of nucleosomal DNA may facilitate translational positioning.