Depletion Effects Massively Change Chromatin Properties and Influence Genome Folding

We present a Monte Carlo model for genome folding at the 30-nm scale with focus on linker-histone and nucleosome depletion effects. We find that parameter distributions from experimental data do not lead to one specific chromatin fiber structure, but instead to a distribution of structures in the chromatin phase diagram. Depletion of linker histones and nucleosomes affects, massively, the flexibility and the extension of chromatin fibers. Increasing the amount of nucleosome skips (i.e., nucleosome depletion) can lead either to a collapse or to a swelling of chromatin fibers. These opposing effects are discussed and we show that depletion effects may even contribute to chromatin compaction. Furthermore, we find that predictions from experimental data for the average nucleosome skip rate lie exactly in the regime of maximum chromatin compaction. Finally, we determine the pair distribution function of chromatin. This function reflects the structure of the fiber, and its Fourier-transform can be measured experimentally. Our calculations show that even in the case of fibers with depletion effects, the main dominant peaks (characterizing the structure and the length scales) can still be identified.     

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.     

Nucleosome geometry and internucleosomal interactions control the chromatin fiber conformation

Based on model structures with atomic resolution, a coarse-grained model for the nucleosome geometry was implemented. The dependence of the chromatin fiber conformation on the spatial orientation of nucleosomes and the path and length of the linker DNA was systematically explored by Monte Carlo simulations. Two fiber types were analyzed in detail that represent nucleosome chains without and with linker histones, respectively: two-start helices with crossed-linker DNA (CL conformation) and interdigitated one-start helices (ID conformation) with different nucleosome tilt angles. The CL conformation was derived from a tetranucleosome crystal structure that was extended into a fiber. At thermal equilibrium, the fiber shape persisted but relaxed into a structure with a somewhat lower linear mass density of 3.1 ± 0.1 nucleosomes/11 nm fiber. Stable ID fibers required local nucleosome tilt angles between 40° and 60°. For these configurations, much higher mass densities of up to 7.9 ± 0.2 nucleosomes/11 nm fiber were obtained. A model is proposed, in which the transition between a CL and ID fiber is mediated by relatively small changes of the local nucleosome geometry. These were found to be in very good agreement with changes induced by linker histone H1 binding as predicted from the high resolution model structures.     

Changing Chromatin Fiber Conformation by Nucleosome Repositioning

Chromatin conformation is dynamic and heterogeneous with respect to nucleosome positions, which can be changed by chromatin remodeling complexes in the cell. These molecular machines hydrolyze ATP to translocate or evict nucleosomes, and establish loci with regularly and more irregularly spaced nucleosomes as well as nucleosome-depleted regions. The impact of nucleosome repositioning on the three-dimensional chromatin structure is only poorly understood. Here, we address this issue by using a coarse-grained computer model of arrays of 101 nucleosomes considering several chromatin fiber models with and without linker histones, respectively. We investigated the folding of the chain in dependence of the position of the central nucleosome by changing the length of the adjacent linker DNA in basepair steps. We found in our simulations that these translocations had a strong effect on the shape and properties of chromatin fibers: i), Fiber curvature and flexibility at the center were largely increased and long-range contacts between distant nucleosomes on the chain were promoted. ii), The highest destabilization of the fiber conformation occurred for a nucleosome shifted by two basepairs from regular spacing, whereas effects of linker DNA changes of ∼10 bp in phase with the helical twist of DNA were minimal. iii), A fiber conformation can stabilize a regular spacing of nucleosomes inasmuch as favorable stacking interactions between nucleosomes are facilitated. This can oppose nucleosome translocations and increase the energetic costs for chromatin remodeling. Our computational modeling framework makes it possible to describe the conformational heterogeneity of chromatin in terms of nucleosome positions, and thus advances theoretical models toward a better understanding of how genome compaction and access are regulated within the cell.     

Nucleosome positioning in relation to nucleosome spacing and DNA sequence-specific binding of a protein

Nucleosome positioning is an important mechanism for the regulation of eukaryotic gene expression. Folding of the chromatin fiber can influence nucleosome positioning, whereas similar electrostatic mechanisms govern the nucleosome repeat length and chromatin fiber folding in vitro. The position of the nucleosomes is directed either by the DNA sequence or by the boundaries created due to the binding of certain trans-acting factors to their target sites in the DNA. Increasing ionic strength results in an increase in nucleosome spacing on the chromatin assembled by the S-190 extract of Drosophila embryos. In this study, a mutant lac repressor protein R3 was used to find the mechanisms of nucleosome positioning on a plasmid with three R3-binding sites. With increasing ionic strength in the presence of R3, the number of positioned nucleosomes in the chromatin decreased, whereas the internucleosomal spacings of the positioned nucleosomes in a single register did not change. The number of the positioned nucleosomes in the chromatin assembled in vitro over different plasmid DNAs with 1-3 lac operators changed with the relative position and number of the R3-binding sites. We found that in the presence of R3, nucleosomes were positioned in the salt gradient method of the chromatin assembly, even in the absence of a nucleosome-positioning sequence. Our results show that nucleosome-positioning mechanisms are dominant, as the nucleosomes can be positioned even in the absence of regular spacing mechanisms. The protein-generated boundaries are more effective when more than one binding site is present with a minimum distance of approximately 165 bp, greater than the nucleosome core DNA length, between them.     

Exploring the Conformational Space of Chromatin Fibers and Their Stability by Numerical Dynamic Phase Diagrams

The three-dimensional structure of chromatin affects DNA accessibility and is therefore a key regulator of gene expression. However, the path of the DNA between consecutive nucleosomes, and the resulting chromatin fiber organization remain controversial. The conformational space available for the folding of the nucleosome chain has been analytically described by phase diagrams with a two-angle model, which describes the chain trajectory by a DNA entry-exit angle at the nucleosome and a torsion angle between consecutive nucleosomes. Here, a novel type of numerical phase diagrams is introduced that relates the geometric phase space to the energy associated with a given chromatin conformation. The resulting phase diagrams revealed differences in the energy landscape that reflect the probability of a given conformation to form in thermal equilibrium. Furthermore, we investigated the effects of entropy and additional degrees of freedom in the dynamic phase diagrams by performing Monte Carlo simulations of the initial chain trajectories. Using our approach, we were able to demonstrate that conformations that initially were geometrically impossible could evolve into energetically favorable states in thermal equilibrium due to DNA bending and torsion. In addition, dynamic phase diagrams were applied to identify chromatin fibers that reflect certain experimentally determined features.     

Changing Chromatin Fiber Conformation by Nucleosome Repositioning

Chromatin conformation is dynamic and heterogeneous with respect to nucleosome positions, which can be changed by chromatin remodeling complexes in the cell. These molecular machines hydrolyze ATP to translocate or evict nucleosomes, and establish loci with regularly and more irregularly spaced nucleosomes as well as nucleosome-depleted regions. The impact of nucleosome repositioning on the three-dimensional chromatin structure is only poorly understood. Here, we address this issue by using a coarse-grained computer model of arrays of 101 nucleosomes considering several chromatin fiber models with and without linker histones, respectively. We investigated the folding of the chain in dependence of the position of the central nucleosome by changing the length of the adjacent linker DNA in basepair steps. We found in our simulations that these translocations had a strong effect on the shape and properties of chromatin fibers: i), Fiber curvature and flexibility at the center were largely increased and long-range contacts between distant nucleosomes on the chain were promoted. ii), The highest destabilization of the fiber conformation occurred for a nucleosome shifted by two basepairs from regular spacing, whereas effects of linker DNA changes of ∼10 bp in phase with the helical twist of DNA were minimal. iii), A fiber conformation can stabilize a regular spacing of nucleosomes inasmuch as favorable stacking interactions between nucleosomes are facilitated. This can oppose nucleosome translocations and increase the energetic costs for chromatin remodeling. Our computational modeling framework makes it possible to describe the conformational heterogeneity of chromatin in terms of nucleosome positions, and thus advances theoretical models toward a better understanding of how genome compaction and access are regulated within the cell.     

Histone tails cooperate to control the breathing of genomic nucleosomes

Genomic DNA is packaged in chromatin, a dynamic fiber variable in size and compaction. In chromatin, repeating nucleosome units wrap 145–147 DNA basepairs around histone proteins. Genetic and epigenetic regulation of genes relies on structural transitions in chromatin which are driven by intra- and inter-nucleosome dynamics and modulated by chemical modifications of the unstructured terminal tails of histones. Here we demonstrate how the interplay between histone H3 and H2A tails control ample nucleosome breathing motions. We monitored large openings of two genomic nucleosomes, and only moderate breathing of an engineered nucleosome in atomistic molecular simulations amounting to 24 μs. Transitions between open and closed nucleosome conformations were mediated by the displacement and changes in compaction of the two histone tails. These motions involved changes in the DNA interaction profiles of clusters of epigenetic regulatory aminoacids in the tails. Removing the histone tails resulted in a large increase of the amplitude of nucleosome breathing but did not change the sequence dependent pattern of the motions. Histone tail modulated nucleosome breathing is a key mechanism of chromatin dynamics with important implications for epigenetic regulation.     

Histone Octamer Helical Tubes Suggest that an Internucleosomal Four-Helix Bundle Stabilizes the Chromatin Fiber

A major question in chromatin involves the exact organization of nucleosomes within the 30-nm chromatin fiber and its structural determinants of assembly. Here we investigate the structure of histone octamer helical tubes via the method of iterative helical real-space reconstruction. Accurate placement of the x-ray structure of the histone octamer within the reconstructed density yields a pseudoatomic model for the entire helix, and allows precise identification of molecular interactions between neighboring octamers. One such interaction that would not be obscured by DNA in the nucleosome consists of a twofold symmetric four-helix bundle formed between pairs of H2B-α3 and H2B-αC helices of neighboring octamers. We believe that this interface can act as an internucleosomal four-helix bundle within the context of the chromatin fiber. The potential relevance of this interface in the folding of the 30-nm chromatin fiber is discussed.     

Salt-dependent compaction of di- and trinucleosomes studied by small-angle neutron scattering

Using small-angle neutron scattering (SANS), we have measured the salt-dependent static structure factor of di- and trinucleosomes from chicken erythrocytes and from COS-7 cells. We also determined the sedimentation coefficients of these dinucleosomes and dinucleosomes reconstituted on a 416-bp DNA containing two nucleosome positioning sequences of the 5S rDNA of Lytechinus variegatus at low and high salt concentrations. The internucleosomal distance d was calculated by simulation as well as Fourier back-transformation of the SANS curves and by hydrodynamic simulation of sedimentation coefficients. Nucleosome dimers from chicken erythrocyte chromatin show a decrease in d from approximately 220 A at 5 mM NaCl to 150 A at 100 mM NaCl. For dinucleosomes from COS-7 chromatin, d decreases from 180 A at 5 mM to 140 A at 100 mM NaCl concentration. Our measurements on trinucleosomes are compatible with a compaction through two different mechanisms, depending on the salt concentration. Between 0 and 20 mM NaCl, the internucleosomal distance between adjacent nucleosomes remains constant, whereas the angle of the DNA strands entering and leaving the central nucleosome decreases. Above 20 mM NaCl, the adjacent nucleosomes approach each other, similar to the compaction of dinucleosomes. The internucleosomal distance of 140-150 A at 100 mM NaCl is in agreement with distances measured by scanning force microscopy and electron microscopy on long chromatin filaments.     

Differences in nanoscale organization of regulatory active and inactive human chromatin

Methodological advances in conformation capture techniques have fundamentally changed our understanding of chromatin architecture. However, the nanoscale organization of chromatin and its cell-to-cell variance are less studied. Analyzing genome-wide data from 733 human cell and tissue samples, we identified 2 prototypical regions that exhibit high or absent hypersensitivity to deoxyribonuclease I, respectively. These regulatory active or inactive regions were examined in the lymphoblast cell line K562 by using high-throughput super-resolution microscopy. In both regions, we systematically measured the physical distance of 2 fluorescence in situ hybridization spots spaced by only 5 kb of DNA. Unexpectedly, the resulting distance distributions range from very compact to almost elongated configurations of more than 200-nm length for both the active and inactive regions. Monte Carlo simulations of a coarse-grained model of these chromatin regions based on published data of nucleosome occupancy in K562 cells were performed to understand the underlying mechanisms. There was no parameter set for the simulation model that can explain the microscopically measured distance distributions. Obviously, the chromatin state given by the strength of internucleosomal interaction, nucleosome occupancy, or amount of histone H1 differs from cell to cell, which results in the observed broad distance distributions. This large variability was not expected, especially in inactive regions. The results for the mechanisms for different distance distributions on this scale are important for understanding the contacts that mediate gene regulation. Microscopic measurements show that the inactive region investigated here is expected to be embedded in a more compact chromatin environment. The simulation results of this region require an increase in the strength of internucleosomal interactions. It may be speculated that the higher density of chromatin is caused by the increased internucleosomal interaction strength.     

Organization of internucleosomal DNA in rat liver chromatin

A detailed analysis of the length distribution of DNA in nucleosome dimers trimmed with exonuclease III and S1 nuclease suggests that the previously described variation of internucleosomal distance in rat liver occurs, at least for a subset of the nucleosomes, by integral multiples of the helical repeat of the DNA. Results obtained upon digestion of chromatin with DNase II further suggest that lengths of internucleosomal DNA are integral multiples of the helical repeat of the DNA plus approximately 5 bp. Restraints imposed by these features on the arrangement of nucleosomes along the fiber are discussed.     

Dynamic condensation of linker histone C-terminal domain regulates chromatin structure

The basic and intrinsically disordered C-terminal domain (CTD) of the linker histone (LH) is essential for chromatin compaction. However, its conformation upon nucleosome binding and its impact on chromatin organization remain unknown. Our mesoscale chromatin model with a flexible LH CTD captures a dynamic, salt-dependent condensation mechanism driven by charge neutralization between the LH and linker DNA. Namely, at low salt concentration, CTD condenses, but LH only interacts with the nucleosome and one linker DNA, resulting in a semi-open nucleosome configuration; at higher salt, LH interacts with the nucleosome and two linker DNAs, promoting stem formation and chromatin compaction. CTD charge reduction unfolds the domain and decondenses chromatin, a mechanism in consonance with reduced counterion screening in vitro and phosphorylated LH in vivo. Divalent ions counteract this decondensation effect by maintaining nucleosome stems and expelling the CTDs to the fiber exterior. Additionally, we explain that the CTD folding depends on the chromatin fiber size, and we show that the asymmetric structure of the LH globular head is responsible for the uneven interaction observed between the LH and the linker DNAs. All these mechanisms may impact epigenetic regulation and higher levels of chromatin folding.     

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.     

Relationship Between Chromatin High-order Folding and Nucleosomal Linker Twist in Nuclei of Human HeLa S3 Cells

Abstract

We have used the intercalative agent ethidium bromide to examine the association between chromatin high-order folding and the twist of internucleosomal DNA regions. The analysis was carried out on intact nuclei isolated from human HeLa S3 cells. Our data shows that alterations in the nucleosomal linker twist significantly influence the way in which a chain of nucleosomes folds to form different higher-order structures. The assay used allowed us to identify the existence of two chromatin fractions differing in their extent of high-order folding. We have also found that active gene sequences are preferentially associated with the chromatin fraction corresponding to the more extended conformation. A model is proposed to account for the effect of variations in the nucleosome linker twist on the state of chromatin folding.