Genomes of Multicellular Organisms Have Evolved to Attract Nucleosomes to Promoter Regions

The importance of DNA sequence for nucleosome positioning in transcriptional regulation

Nucleosome positioning is a key factor for transcriptional regulation. Nucleosomes regulate the dynamic accessibility of chromatin and interact with the transcription machinery at every stage. Influences to steer nucleosome positioning are diverse, and the according importance of the DNA sequence in contrast to active chromatin remodeling has been the subject of long discussion. In this study, we evaluate the functional role of DNA sequence for all major elements along the process of transcription. We developed a random forest classifier based on local DNA structure that assesses the sequence-intrinsic support for nucleosome positioning. On this basis, we created a simple data resource that we applied genome-wide to the human genome. In our comprehensive analysis, we found a special role of DNA in mediating the competition of nucleosomes with cis-regulatory elements, in enabling steady transcription, for positioning of stable nucleosomes in exons, and for repelling nucleosomes during transcription termination. In contrast, we relate these findings to concurrent processes that generate strongly positioned nucleosomes in vivo that are not mediated by sequence, such as energy-dependent remodeling of chromatin.

Stochastic nucleosome disassembly mediated by remodelers and histone fragmentation

We construct and analyze monomeric and multimeric models of the stochastic disassembly of a single nucleosome. Our monomeric model predicts the time needed for a number of histone-DNA contacts to spontaneously break, leading to dissociation of a non-fragmented histone from DNA. The dissociation process can be facilitated by DNA binding proteins or processing molecular motors that compete with histones for histone-DNA contact sites. Eigenvalue analysis of the corresponding master equation allows us to evaluate histone detachment times under both spontaneous detachment and protein-facilitated processes. We find that competitive DNA binding of remodeling proteins can significantly reduce the typical detachment time but only if these remodelers have DNA-binding affinities comparable to those of histone-DNA contact sites. In the presence of processive motors, the histone detachment rate is shown to be proportional to the product of the histone single-bond dissociation constant and the speed of motor protein procession. Our simple intact-histone model is then extended to allow for multimeric nucleosome kinetics that reveal additional pathways of disassembly. In addition to a dependence of complete disassembly times on subunit-DNA contact energies, we show how histone subunit concentrations in bulk solutions can mediate the disassembly process by rescuing partially disassembled nucleosomes. Moreover, our kinetic model predicts that remodeler binding can also bias certain pathways of nucleosome disassembly, with higher remodeler binding rates favoring intact-histone detachment.

Stochastic nucleosome disassembly mediated by remodelers and histone fragmentation

We construct and analyze monomeric and multimeric models of the stochastic disassembly of a single nucleosome. Our monomeric model predicts the time needed for a number of histone-DNA contacts to spontaneously break, leading to dissociation of a non-fragmented histone from DNA. The dissociation process can be facilitated by DNA binding proteins or processing molecular motors that compete with histones for histone-DNA contact sites. Eigenvalue analysis of the corresponding master equation allows us to evaluate histone detachment times under both spontaneous detachment and protein-facilitated processes. We find that competitive DNA binding of remodeling proteins can significantly reduce the typical detachment time but only if these remodelers have DNA-binding affinities comparable to those of histone-DNA contact sites. In the presence of processive motors, the histone detachment rate is shown to be proportional to the product of the histone single-bond dissociation constant and the speed of motor protein procession. Our simple intact-histone model is then extended to allow for multimeric nucleosome kinetics that reveal additional pathways of disassembly. In addition to a dependence of complete disassembly times on subunit-DNA contact energies, we show how histone subunit concentrations in bulk solution can mediate the disassembly process by rescuing partially disassembled nucleosomes. Moreover, our kinetic model predicts that remodeler binding can also bias certain pathways of nucleosome disassembly, with higher remodeler binding rates favoring intact-histone detachment.

The role of transcript regions and amino acid choice in nucleosome positioning

Eukaryotic DNA is organized and compacted in a string of nucleosomes, DNA-wrapped protein cylinders. The positions of nucleosomes along DNA are not random but show well-known base pair sequence preferences that result from the sequence-dependent elastic and geometric properties of the DNA double helix. Here, we focus on DNA around transcription start sites, which are known to typically attract nucleosomes in multicellular life forms through their high GC content. We aim to understand how these GC signals, as observed in genome-wide averages, are produced and encoded through different genomic regions (mainly 5' UTRs, coding exons, and introns). Our study uses a bioinformatics approach to decompose the genome-wide GC signal into between-region and within-region signals. We find large differences in GC signal contributions between vertebrates and plants and, remarkably, even between closely related species. Introns contribute most to the GC signal in vertebrates, while in plants the exons dominate. Further, we find signal strengths stronger on DNA than on mRNA, suggesting a biological function of GC signals along the DNA itself, as is the case for nucleosome positioning. Finally, we make the surprising discovery that both the choice of synonymous codons and amino acids contribute to the nucleosome positioning signal.

Multiplexing mechanical and translational cues on genes

The genetic code gives precise instructions on how to translate codons into amino acids. Due to the degeneracy of the genetic code-18 out of 20 amino acids are encoded for by more than one codon-more information can be stored in a basepair sequence. Indeed, various types of additional information have been discussed in the literature, e.g., the positioning of nucleosomes along eukaryotic genomes and the modulation of the translating efficiency in ribosomes to influence cotranslational protein folding. The purpose of this study is to show that it is indeed possible to carry more than one additional layer of information on top of a gene. In particular, we show how much translation efficiency and nucleosome positioning can be adjusted simultaneously without changing the encoded protein. We achieve this by mapping genes on weighted graphs that contain all synonymous genes, and then finding shortest paths through these graphs. This enables us, for example, to readjust the disrupted translational efficiency profile after a gene has been introduced from one organism (e.g., human) into another (e.g., yeast) without greatly changing the nucleosome landscape intrinsically encoded by the DNA molecule.

Galaxy Dnpatterntools for Computational Analysis of Nucleosome Positioning Sequence Patterns

Nucleosomes are basic units of DNA packing in eukaryotes. Their structure is well conserved from yeast to human and consists of the histone octamer core and 147 bp DNA wrapped around it. Nucleosomes are bound to a majority of the eukaryotic genomic DNA, including its regulatory regions. Hence, they also play a major role in gene regulation. For the latter, their precise positioning on DNA is essential. In the present paper, we describe Galaxy dnpatterntools—software package for nucleosome DNA sequence analysis and mapping. This software will be useful for computational biologists practitioners to conduct more profound studies of gene regulatory mechanisms.

The Active Mechanism of Nucleosome Depletion by Poly(dA:dT) Tracts In Vivo

Poly(dA:dT) tracts cause nucleosome depletion in many species, e.g., at promoters and replication origins. Their intrinsic biophysical sequence properties make them stiff and unfavorable for nucleosome assembly, as probed by in vitro nucleosome reconstitution. The mere correlation between nucleosome depletion over poly(dA:dT) tracts in in vitro reconstituted and in in vivo chromatin inspired an intrinsic nucleosome exclusion mechanism in vivo that is based only on DNA and histone properties. However, we compile here published and new evidence that this correlation does not reflect mechanistic causation. (1) Nucleosome depletion over poly(dA:dT) in vivo is not universal, e.g., very weak in S. pombe. (2) The energy penalty for incorporating poly(dA:dT) tracts into nucleosomes is modest (<10%) relative to ATP hydrolysis energy abundantly invested by chromatin remodelers. (3) Nucleosome depletion over poly(dA:dT) is much stronger in vivo than in vitro if monitored without MNase and (4) actively maintained in vivo. (5) S. cerevisiae promoters evolved a strand-biased poly(dA) versus poly(dT) distribution. (6) Nucleosome depletion over poly(dA) is directional in vivo. (7) The ATP dependent chromatin remodeler RSC preferentially and directionally displaces nucleosomes towards 5′ of poly(dA). Especially distribution strand bias and displacement directionality would not be expected for an intrinsic mechanism. Together, this argues for an in vivo mechanism where active and species-specific read out of intrinsic sequence properties, e.g., by remodelers, shapes nucleosome organization.

Coupling between Sequence-Mediated Nucleosome Organization and Genome Evolution

The nucleosome is a major modulator of DNA accessibility to other cellular factors. Nucleosome positioning has a critical importance in regulating cell processes such as transcription, replication, recombination or DNA repair. The DNA sequence has an influence on the position of nucleosomes on genomes, although other factors are also implicated, such as ATP-dependent remodelers or competition of the nucleosome with DNA binding proteins. Different sequence motifs can promote or inhibit the nucleosome formation, thus influencing the accessibility to the DNA. Sequence-encoded nucleosome positioning having functional consequences on cell processes can then be selected or counter-selected during evolution. We review the interplay between sequence evolution and nucleosome positioning evolution. We first focus on the different ways to encode nucleosome positions in the DNA sequence, and to which extent these mechanisms are responsible of genome-wide nucleosome positioning in vivo. Then, we discuss the findings about selection of sequences for their nucleosomal properties. Finally, we illustrate how the nucleosome can directly influence sequence evolution through its interactions with DNA damage and repair mechanisms. This review aims to provide an overview of the mutual influence of sequence evolution and nucleosome positioning evolution, possibly leading to complex evolutionary dynamics.

Multiple roles of H2A.Z in regulating promoter chromatin architecture in human cells

Chromatin accessibility of a promoter is fundamental in regulating transcriptional activity. The histone variant H2A.Z has been shown to contribute to this regulation, but its role has remained poorly understood. Here, we prepare high-depth maps of the position and accessibility of H2A.Z-containing nucleosomes for all human Pol II promoters in epithelial, mesenchymal and isogenic cancer cell lines. We find that, in contrast to the prevailing model, many different types of active and inactive promoter structures are observed that differ in their nucleosome organization and sensitivity to MNase digestion. Key aspects of an active chromatin structure include positioned H2A.Z MNase resistant nucleosomes upstream or downstream of the TSS, and a MNase sensitive nucleosome at the TSS. Furthermore, the loss of H2A.Z leads to a dramatic increase in the accessibility of transcription factor binding sites. Collectively, these results suggest that H2A.Z has multiple and distinct roles in regulating gene expression dependent upon its location in a promoter.

On the role of transcription in positioning nucleosomes

Nucleosome positioning is crucial for the genome’s function. Though the role of DNA sequence in positioning nucleosomes is well understood, a detailed mechanistic understanding on the impact of transcription remains lacking. Using numerical simulations, we investigated the dependence of nucleosome density profiles on transcription level across multiple species. We found that the low nucleosome affinity of yeast, but not mouse, promoters contributes to the formation of phased nucleosomes arrays for inactive genes. For the active genes, a heterogeneous distribution of +1 nucleosomes, caused by a tug-of-war between two types of remodeling enzymes, is essential for reproducing their density profiles. In particular, while positioning enzymes are known to remodel the +1 nucleosome and align it toward the transcription start site (TSS), spacer enzymes that use a pair of nucleosomes as their substrate can shift the nucleosome array away from the TSS. Competition between these enzymes results in two types of nucleosome density profiles with well- and ill-positioned +1 nucleosome. Finally, we showed that Pol II assisted histone exchange, if occurring at a fast speed, can abolish the impact of remodeling enzymes. By elucidating the role of individual factors, our study reconciles the seemingly conflicting results on the overall impact of transcription in positioning nucleosomes across species.

Nucleosome positioning plays a key role in the genome’s function by regulating the accessibility of protein binding sites as well as higher-order chromatin organization. Though significant progress has been made towards studying the role of DNA sequence in positioning the nucleosomes, our understanding on the impact of transcription lags behind. Our study uses kinetic simulations to explore the role of DNA sequence specificity, transcription factor binding, enzyme remodeling, and Pol II elongation in positioning nucleosomes. It suggests that the differences in nucleosome density profiles observed at various transcription levels in yeast and mouse embryonic stem cells can be understood from a tug-of-war between two types of remodeling enzymes.

Translational nucleosome positioning: A computational study

About three-quarters of eukaryotic DNA is wrapped into nucleosomes; DNA spools with a protein core. The affinity of a given DNA stretch to be incorporated into a nucleosome is known to depend on the base-pair sequence-dependent geometry and elasticity of the DNA double helix. This causes the rotational and translational positioning of nucleosomes. In this study we ask the question whether the latter can be predicted by a simple coarse-grained DNA model with sequence-dependent elasticity, the rigid base-pair model. Whereas this model is known to be rather robust in predicting rotational nucleosome positioning, we show that the translational positioning is a rather subtle effect that is dominated by the guanine-cytosine content dependence of entropy rather than energy. A correct qualitative prediction within the rigid base-pair framework can only be achieved by assuming that DNA elasticity effectively changes on complexation into the nucleosome complex. With that extra assumption we arrive at a model which gives an excellent quantitative agreement to experimental in vitro nucleosome maps, under the additional assumption that nucleosomes equilibrate their positions only locally.

Role of cell-type specific nucleosome positioning in inducible activation of mammalian promoters

The organization of nucleosomes across functional genomic elements represents a critical layer of control. Here, we present a strategy for high-resolution nucleosome profiling at selected genomic features, and use this to analyse dynamic nucleosome positioning at inducible and cell-type-specific mammalian promoters. We find that nucleosome patterning at inducible promoters frequently resembles that at active promoters, even before stimulus-driven activation. Accordingly, the nucleosome profile at many inactive inducible promoters is sufficient to predict cell-type-specific responsiveness. Induction of gene expression is generally not associated with major changes to nucleosome patterning, and a subset of inducible promoters can be activated without stable nucleosome depletion from their transcription start sites. These promoters are generally dependent on remodelling enzymes for their inducible activation, and exhibit transient nucleosome depletion only at alleles undergoing transcription initiation. Together, these data reveal how the responsiveness of inducible promoters to activating stimuli is linked to cell-type-specific nucleosome patterning.

Nucleosome organisation plays important roles in regulating functional genomic elements. Here, the authors use high-resolution profiling to analyse dynamic nucleosome positioning at inducible and cell-type-specific promoters, providing a global view of chromatin architecture at inducible promoters.

Nucleosome positioning sequence patterns as packing or regulatory

Nucleosome positioning DNA sequence patterns (NPS)-usually distributions of particular dinucleotides or other sequence elements in nucleosomal DNA-at least partially determine chromatin structure and arrangements of nucleosomes that in turn affect gene expression. Statistically, NPS are defined as oscillations of the dinucleotide periodicity of about 10 base pairs (bp) which reflects the double helix period. We compared the nucleosomal DNA patterns in mouse, human and yeast organisms and observed few distinctive patterns that can be termed as packing and regulatory referring to distinctive modes of chromatin function. For the first time the NPS patterns in nucleus accumbens cells (NAC) in mouse brain were characterized and compared to the patterns in human CD4+ and apoptotic lymphocyte cells and well studied patterns in yeast. The NPS patterns in human CD4+ cells and mouse brain cells had very high positive correlation. However, there was no correlation between them and patterns in human apoptotic lymphocyte cells and yeast, but the latter two were highly correlated with each other. By their dinucleotide arrangements the analyzed NPS patterns classified into stable canonical WW/SS (W = A or T and S = C or G dinucleotide) and less stable RR/YY (R = A or G and Y = C or T dinucleotide) patterns and anti-patterns. In the anti-patterns positioning of the dinucleotides is flipped compared to those in the regular patterns. Stable canonical WW/SS patterns and anti-patterns are ubiquitously observed in many organisms and they had high resemblance between yeast and human apoptotic cells. Less stable RR/YY patterns had higher positive correlation between mouse and normal human cells. Our analysis and evidence from scientific literature lead to idea that various distinct patterns in nucleosomal DNA can be related to the two roles of the chromatin: packing (WW/SS) and regulatory (RR/YY and "anti").

Ensembles of Breathing Nucleosomes: A Computational Study

About three-fourths of the human DNA molecules are wrapped into nucleosomes, protein spools with DNA. Nucleosomes are highly dynamic, transiently exposing their DNA through spontaneous unspooling. Recent experiments allowed to observe the DNA of an ensemble of such breathing nucleosomes through x-ray diffraction with contrast matching between the solvent and the protein core. In this study, we calculate such an ensemble through a Monte Carlo simulation of a coarse-grained nucleosome model with sequence-dependent DNA mechanics. Our analysis gives detailed insights into the sequence dependence of nucleosome breathing observed in the experiment and allows us to determine the adsorption energy of the DNA bound to the protein core as a function of the ionic strength. Moreover, we predict the breathing behavior of other potentially interesting sequences and compare the findings to earlier related experiments.

In vitro versus in vivo compositional landscapes of histone sequence preferences in eucaryotic genomes

Motivation

Although the nucleosome occupancy along a genome can be in part predicted by in vitro experiments, it has been recently observed that the chromatin organization presents important differences in vitro with respect to in vivo. Such differences mainly regard the hierarchical and regular structures of the nucleosome fiber, whose existence has long been assumed, and in part also observed in vitro, but that does not apparently occur in vivo. It is also well known that the DNA sequence has a role in determining the nucleosome occupancy. Therefore, an important issue is to understand if, and to what extent, the structural differences in the chromatin organization between in vitro and in vivo have a counterpart in terms of the underlying genomic sequences.

Results

We present the first quantitative comparison between the in vitro and in vivo nucleosome maps of two model organisms (S. cerevisiae and C. elegans). The comparison is based on the construction of weighted k-mer dictionaries. Our findings show that there is a good level of sequence conservation between in vitro and in vivo in both the two organisms, in contrast to the abovementioned important differences in chromatin structural organization. Moreover, our results provide evidence that the two organisms predispose themselves differently, in terms of sequence composition and both in vitro and in vivo, for the nucleosome occupancy. This leads to the conclusion that, although the notion of a genome encoding for its own nucleosome occupancy is general, the intrinsic histone k-mer sequence preferences tend to be species-specific.

Availability and implementation

The files containing the dictionaries and the main results of the analysis are available at http://math.unipa.it/rombo/material.

Supplementary information

Supplementary data are available at Bioinformatics online.

The nucleosome position-encoding WW/SS sequence pattern is depleted in mammalian genes relative to other eukaryotes

Nucleosomal DNA sequences generally follow a well-known pattern with ∼10-bp periodic WW (where W is A or T) dinucleotides that oscillate in phase with each other and out of phase with SS (where S is G or C) dinucleotides. However, nucleosomes with other DNA patterns have not been systematically analyzed. Here, we focus on an opposite pattern, namely anti-WW/SS pattern, in which WW dinucleotides preferentially occur at DNA sites that bend into major grooves and SS (where S is G or C) dinucleotides are often found at sites that bend into minor grooves. Nucleosomes with the anti-WW/SS pattern are widespread and exhibit a species- and context-specific distribution in eukaryotic genomes. Unlike non-mammals (yeast, nematode and fly), there is a positive correlation between the enrichment of anti-WW/SS nucleosomes and RNA Pol II transcriptional levels in mammals (mouse and human). Interestingly, such enrichment is not due to underlying DNA sequence. In addition, chromatin remodeling complexes have an impact on the abundance but not on the distribution of anti-WW/SS nucleosomes in yeast. Our data reveal distinct roles of cis- and trans-acting factors in the rotational positioning of nucleosomes between non-mammals and mammals. Implications of the anti-WW/SS sequence pattern for RNA Pol II transcription are discussed.

Shortest paths through synonymous genomes

The elasticity of the DNA double helix varies with the underlying base pair sequence. This allows one to put mechanical cues into sequences that in turn influence the packaging of DNA into nucleosomes, DNA-wrapped protein cylinders. Nucleosomes dictate a broad range of biological processes, ranging from gene regulation, recombination, and replication to chromosome condensation. Here we map base pair sequences onto graphs and use shortest paths algorithms to determine which DNA stretches are easiest or hardest to bend inside a nucleosome. We further demonstrate how genetic and mechanical information can be multiplexed by studying paths through graphs of synonymous codons. Using this method we find that nucleosomes can be placed by mechanical cues nearly everywhere on the genome of baker's yeast (Saccharomyces cerevisiae).

DNA sequence encodes the position of DNA supercoils

The three-dimensional organization of DNA is increasingly understood to play a decisive role in vital cellular processes. Many studies focus on the role of DNA-packaging proteins, crowding, and confinement in arranging chromatin, but structural information might also be directly encoded in bare DNA itself. Here, we visualize plectonemes (extended intertwined DNA structures formed upon supercoiling) on individual DNA molecules. Remarkably, our experiments show that the DNA sequence directly encodes the structure of supercoiled DNA by pinning plectonemes at specific sequences. We develop a physical model that predicts that sequence-dependent intrinsic curvature is the key determinant of pinning strength and demonstrate this simple model provides very good agreement with the data. Analysis of several prokaryotic genomes indicates that plectonemes localize directly upstream of promoters, which we experimentally confirm for selected promotor sequences. Our findings reveal a hidden code in the genome that helps to spatially organize the chromosomal DNA.

Structure and function of archaeal histones

The genomes of all organisms throughout the tree of life are compacted and organized in chromatin by association of chromatin proteins. Eukaryotic genomes encode histones, which are assembled on the genome into octamers, yielding nucleosomes. Post-translational modifications of the histones, which occur mostly on their N-terminal tails, define the functional state of chromatin. Like eukaryotes, most archaeal genomes encode histones, which are believed to be involved in the compaction and organization of their genomes. Instead of discrete multimers, in vivo data suggest assembly of “nucleosomes” of variable size, consisting of multiples of dimers, which are able to induce repression of transcription. Based on these data and a model derived from X-ray crystallography, it was recently proposed that archaeal histones assemble on DNA into “endless” hypernucleosomes. In this review, we discuss the amino acid determinants of hypernucleosome formation and highlight differences with the canonical eukaryotic octamer. We identify archaeal histones differing from the consensus, which are expected to be unable to assemble into hypernucleosomes. Finally, we identify atypical archaeal histones with short N- or C-terminal extensions and C-terminal tails similar to the tails of eukaryotic histones, which are subject to post-translational modification. Based on the expected characteristics of these archaeal histones, we discuss possibilities of involvement of histones in archaeal transcription regulation.

Both Archaea and eukaryotes express histones, but whereas the tertiary structure of histones is conserved, the quaternary structure of histone–DNA complexes is very different. In a recent study, the crystal structure of the archaeal hypernucleosome was revealed to be an “endless” core of interacting histones that wraps the DNA around it in a left-handed manner. The ability to form a hypernucleosome is likely determined by dimer–dimer interactions as well as stacking interactions between individual layers of the hypernucleosome. We analyzed a wide variety of archaeal histones and found that most but not all histones possess residues able to facilitate hypernucleosome formation. Among these are histones with truncated termini or extended histone tails. Based on our analysis, we propose several possibilities of archaeal histone involvement in transcription regulation.

Evidence for DNA Sequence Encoding of an Accessible Nucleosomal Array across Vertebrates

Nucleosome-depleted regions around which nucleosomes order following the "statistical" positioning scenario were recently shown to be encoded in the DNA sequence in human. This intrinsic nucleosomal ordering strongly correlates with oscillations in the local GC content as well as with the interspecies and intraspecies mutation profiles, revealing the existence of both positive and negative selection. In this letter, we show that these predicted nucleosome inhibitory energy barriers (NIEBs) with compacted neighboring nucleosomes are indeed ubiquitous to all vertebrates tested. These 1 kb-sized chromatin patterns are widely distributed along vertebrate chromosomes, overall covering more than a third of the genome. We have previously observed in human deviations from neutral evolution at these genome-wide distributed regions, which we interpreted as a possible indication of the selection of an open, accessible, and dynamic nucleosomal array to constitutively facilitate the epigenetic regulation of nuclear functions in a cell-type-specific manner. As a first, very appealing observation supporting this hypothesis, we report evidence of a strong association between NIEB borders and the poly(A) tails of Alu sequences in human. These results suggest that NIEBs provide adequate chromatin patterns favorable to the integration of Alu retrotransposons and, more generally to various transposable elements in the genomes of primates and other vertebrates.

The Latest Twists in Chromatin Remodeling

In its most restrictive interpretation, the notion of chromatin remodeling refers to the action of chromatin-remodeling enzymes on nucleosomes with the aim of displacing and removing them from the chromatin fiber (the effective polymer formed by a DNA molecule and proteins). This local modification of the fiber structure can have consequences for the initiation and repression of the transcription process, and when the remodeling process spreads along the fiber, it also results in long-range effects essential for fiber condensation. There are three regulatory levels of relevance that can be distinguished for this process: the intrinsic sequence preference of the histone octamer, which rules the positioning of the nucleosome along the DNA, notably in relation to the genetic information coded in DNA; the recognition or selection of nucleosomal substrates by remodeling complexes; and, finally, the motor action on the nucleosome exerted by the chromatin remodeler. Recent work has been able to provide crucial insights at each of these three levels that add new twists to this exciting and unfinished story, which we highlight in this perspective.

The role of DNA sequence in nucleosome breathing

Roughly 3/4 of human genomes are sequestered by nucleosomes, DNA spools with a protein core, dictating a broad range of biological processes, ranging from gene regulation, recombination, and replication, to chromosome condensation. Nucleosomes are dynamical structures and temporarily expose wrapped DNA through spontaneous unspooling from either end, a process called site exposure or nucleosome breathing. Here we ask how this process is influenced by the mechanical properties of the wrapped DNA, which is known to depend on the underlying base pair sequence. Using a coarse-grained nucleosome model we calculate the accessibility profiles for site exposure. We find that the process is very sensitive to sequence effects, so that evolution could potentially tune the accessibility of nucleosomal DNA and would only need a small number of mutations to do so.

Physics behind the mechanical nucleosome positioning code

The positions along DNA molecules of nucleosomes, the most abundant DNA-protein complexes in cells, are influenced by the sequence-dependent DNA mechanics and geometry. This leads to the "nucleosome positioning code", a preference of nucleosomes for certain sequence motives. Here we introduce a simplified model of the nucleosome where a coarse-grained DNA molecule is frozen into an idealized superhelical shape. We calculate the exact sequence preferences of our nucleosome model and find it to reproduce qualitatively all the main features known to influence nucleosome positions. Moreover, using well-controlled approximations to this model allows us to come to a detailed understanding of the physics behind the sequence preferences of nucleosomes.

Single-molecule force spectroscopy on histone H4 tail-cross-linked chromatin reveals fiber folding

The eukaryotic genome is highly compacted into a protein-DNA complex called chromatin. The cell controls access of transcriptional regulators to chromosomal DNA via several mechanisms that act on chromatin-associated proteins and provide a rich spectrum of epigenetic regulation. Elucidating the mechanisms that fold chromatin fibers into higher-order structures is therefore key to understanding the epigenetic regulation of DNA accessibility. Here, using histone H4-V21C and histone H2A-E64C mutations, we employed single-molecule force spectroscopy to measure the unfolding of individual chromatin fibers that are reversibly cross-linked through the histone H4 tail. Fibers with covalently linked nucleosomes featured the same folding characteristics as fibers containing wild-type histones but exhibited increased stability against stretching forces. By stabilizing the secondary structure of chromatin, we confirmed a nucleosome repeat length (NRL)-dependent folding. Consistent with previous crystallographic and cryo-EM studies, the obtained force-extension curves on arrays with 167-bp NRLs best supported an underlying structure consisting of zig-zag, two-start fibers. For arrays with 197-bp NRLs, we previously inferred solenoidal folding, which was further corroborated by force-extension curves of the cross-linked fibers. The different unfolding pathways exhibited by these two types of arrays and reported here extend our understanding of chromatin structure and its potential roles in gene regulation. Importantly, these findings imply that chromatin compaction by nucleosome stacking protects nucleosomal DNA from external forces up to 4 piconewtons.

Designing nucleosomal force sensors

About three quarters of our DNA is wrapped into nucleosomes: DNA spools with a protein core. It is well known that the affinity of a given DNA stretch to be incorporated into a nucleosome depends on the geometry and elasticity of the basepair sequence involved, causing the positioning of nucleosomes. Here we show that DNA elasticity can have a much deeper effect on nucleosomes than just their positioning: it affects their "identities". Employing a recently developed computational algorithm, the mutation Monte Carlo method, we design nucleosomes with surprising physical characteristics. Unlike any other nucleosomes studied so far, these nucleosomes are short-lived when put under mechanical tension whereas other physical properties are largely unaffected. This suggests that the nucleosome, the most abundant DNA-protein complex in our cells, might more properly be considered a class of complexes with a wide array of physical properties, and raises the possibility that evolution has shaped various nucleosome species according to their genomic context.

Benchmarking and refining probability-based models for nucleosome-DNA interaction

Background

In investigations of nucleosome positioning preferences, a model that assigns an affinity to a given sequence is necessary to make predictions. One important class of models, which treats a nucleosome sequence as a Markov chain, has been applied with success when informed with experimentally measured nucleosomal sequence preferences.

Results

We find that we can also use such models as a fast approximative scheme for computationally expensive biophysical models, vastly increasing their reach. Employing these models in this way also allows us to benchmark them for the first time. Doing so for the approximative in silico models indirectly tells us about the accuracy we can expect of them when applied to real data.

Conclusion

We find that models presented in the literature should perform well, but this performance depends on factors such as the order of the Markov model, the preprocessing of the probability distributions on which the model is based, and the size and quality of the sequence ensemble from which those distributions are calculated.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-017-1569-0) contains supplementary material, which is available to authorized users.