An all-atom model of the chromatin fiber containing linker histones reveals a versatile structure tuned by the nucleosomal repeat length

In the nucleus of eukaryotic cells, histone proteins organize the linear genome into a functional and hierarchical architecture. In this paper, we use the crystal structures of the nucleosome core particle, B-DNA and the globular domain of H5 linker histone to build the first all-atom model of compact chromatin fibers. In this 3D jigsaw puzzle, DNA bending is achieved by solving an inverse kinematics problem. Our model is based on recent electron microscopy measurements of reconstituted fiber dimensions. Strikingly, we find that the chromatin fiber containing linker histones is a polymorphic structure. We show that different fiber conformations are obtained by tuning the linker histone orientation at the nucleosomes entry/exit according to the nucleosomal repeat length. We propose that the observed in vivo quantization of nucleosomal repeat length could reflect nature's ability to use the DNA molecule's helical geometry in order to give chromatin versatile topological and mechanical properties.     

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.     

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.     

Local geometry and elasticity in compact chromatin structure

The hierarchical packaging of DNA into chromatin within a eukaryotic nucleus plays a pivotal role in both the accessibility of genomic information and the dynamics of replication. Our work addresses the role of nanoscale physical and geometric properties in determining the structure of chromatin at the mesoscale level. We study the packaging of DNA in chromatin fibers by optimization of regular helical morphologies, considering the elasticity of the linker DNA as well as steric packing of the nucleosomes and linkers. Our model predicts a broad range of preferred helix structures for a fixed linker length of DNA; changing the linker length alters the predicted ensemble. Specifically, we find that the twist registry of the nucleosomes, as set by the internucleosome repeat length, determines the preferred angle between the nucleosomes and the fiber axis. For moderate to long linker lengths, we find a number of energetically comparable configurations with different nucleosome-nucleosome interaction patterns, indicating a potential role for kinetic trapping in chromatin fiber formation. Our results highlight the key role played by DNA elasticity and local geometry in regulating the hierarchical packaging of the genome.     

Effect of the length of internucleosome regions of DNA on the conformation of chromatid fibrils

Models of chromatin fibers structures with linear regions of linker DNA were analysed. Limitations put by end dimensions of linker DNA and nucleosomes are considered. Good agreement between the structural properties of model and real chromatin fibers was obtained. It has been shown that the models with three and more configurations of closely located nucleosomes have linker DNA of 19-53 base pairs length, which is characteristic of real chromatin of the majority of somatic cells.     

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.     

Quantitative analysis of single-molecule force spectroscopy on folded chromatin fibers

Single-molecule techniques allow for picoNewton manipulation and nanometer accuracy measurements of single chromatin fibers. However, the complexity of the data, the heterogeneity of the composition of individual fibers and the relatively large fluctuations in extension of the fibers complicate a structural interpretation of such force-extension curves. Here we introduce a statistical mechanics model that quantitatively describes the extension of individual fibers in response to force on a per nucleosome basis. Four nucleosome conformations can be distinguished when pulling a chromatin fiber apart. A novel, transient conformation is introduced that coexists with single wrapped nucleosomes between 3 and 7 pN. Comparison of force-extension curves between single nucleosomes and chromatin fibers shows that embedding nucleosomes in a fiber stabilizes the nucleosome by 10 k_BT. Chromatin fibers with 20- and 50-bp linker DNA follow a different unfolding pathway. These results have implications for accessibility of DNA in fully folded and partially unwrapped chromatin fibers and are vital for understanding force unfolding experiments on nucleosome arrays.     

Preferentially quantized linker DNA lengths in Saccharomyces cerevisiae

The exact lengths of linker DNAs connecting adjacent nucleosomes specify the intrinsic three-dimensional structures of eukaryotic chromatin fibers. Some studies suggest that linker DNA lengths preferentially occur at certain quantized values, differing one from another by integral multiples of the DNA helical repeat, approximately 10 bp; however, studies in the literature are inconsistent. Here, we investigate linker DNA length distributions in the yeast Saccharomyces cerevisiae genome, using two novel methods: a Fourier analysis of genomic dinucleotide periodicities adjacent to experimentally mapped nucleosomes and a duration hidden Markov model applied to experimentally defined dinucleosomes. Both methods reveal that linker DNA lengths in yeast are preferentially periodic at the DNA helical repeat ( approximately 10 bp), obeying the forms 10n+5 bp (integer n). This 10 bp periodicity implies an ordered superhelical intrinsic structure for the average chromatin fiber in yeast.     

Chromatin fibers are left-handed double helices with diameter and mass per unit length that depend on linker length

Four classes of models have been proposed for the internal structure of eukaryotic chromosome fibers--the solenoid, twisted-ribbon, crossed-linker, and superbead models. We have collected electron image and x-ray scattering data from nuclei, and isolated chromatin fibers of seven different tissues to distinguish between these models. The fiber diameters are related to the linker lengths by the equation: D(N) = 19.3 + 0.23 N, where D(N) is the external diameter (nm) and N is the linker length (base pairs). The number of nucleosomes per unit length of the fibers is also related to linker length. Detailed studies were done on the highly regular chromatin from erythrocytes of Necturus (mud puppy) and sperm of Thyone (sea cucumber). Necturus chromatin fibers (N = 48 bp) have diameters of 31 nm and have 7.5 +/- 1 nucleosomes per 10 nm along the axis. Thyone chromatin fibers (N = 87 bp) have diameters of 39 nm and have 12 +/- 2 nucleosomes per 10 nm along the axis. Fourier transforms of electron micrographs of Necturus fibers showed left-handed helical symmetry with a pitch of 25.8 +/- 0.8 nm and pitch angle of 32 +/- 3 degrees, consistent with a double helix. Comparable conclusions were drawn from the Thyone data. The data do not support the solenoid, twisted-ribbon, or supranucleosomal particle models. The data do support two crossed-linker models having left-handed double-helical symmetry and conserved nucleosome interactions.     

The Dynamic Influence of Linker Histone Saturation within the Poly-Nucleosome Array

Linker histones bind to nucleosomes and modify chromatin structure and dynamics as a means of epigenetic regulation. Biophysical studies have shown that chromatin fibers can adopt a plethora of conformations with varying levels of compaction. Linker histone condensation, and its specific binding disposition, has been associated with directly tuning this ensemble of states. However, the atomistic dynamics and quantification of this mechanism remains poorly understood. Here, we present molecular dynamics simulations of octa-nucleosome arrays, based on a cryo-EM structure of the 30-nm chromatin fiber, with and without the globular domains of the H1 linker histone to determine how they influence fiber structures and dynamics. Results show that when bound, linker histones inhibit DNA flexibility and stabilize repeating tetra-nucleosomal units, giving rise to increased chromatin compaction. Furthermore, upon the removal of H1, there is a significant destabilization of this compact structure as the fiber adopts less strained and untwisted states. Interestingly, linker DNA sampling in the octa-nucleosome is exaggerated compared to its mono-nucleosome counterparts, suggesting that chromatin architecture plays a significant role in DNA strain even in the absence of linker histones. Moreover, H1-bound states are shown to have increased stiffness within tetra-nucleosomes, but not between them. This increased stiffness leads to stronger long-range correlations within the fiber, which may result in the propagation of epigenetic signals over longer spatial ranges. These simulations highlight the effects of linker histone binding on the internal dynamics and global structure of poly-nucleosome arrays, while providing physical insight into a mechanism of chromatin compaction.     

Linker histone tails and N-tails of histone H3 are redundant: scanning force microscopy studies of reconstituted fibers.

The mechanisms responsible for organizing linear arrays of nucleosomes into the three-dimensional structure of chromatin are still largely unknown. In a companion paper (Leuba, S. H., et al. 1998. Biophys. J. 74:2823-2829), we study the contributions of linker histone domains and the N-terminal tail of core histone H3 to extended chromatin fiber structure by scanning force microscopy imaging of mildly trypsinized fibers. Here we complement and extend these studies by scanning force microscopy imaging of selectively reconstituted chromatin fibers, which differ in subtle but distinctive ways in their histone composition. We demonstrate an absolute requirement for the globular domain of the linker histones and a structural redundancy of the tails of linker histones and of histone H3 in determining conformational stability.     

Small angle x-ray scattering of chromatin. Radius and mass per unit length depend on linker length.

Analyses of low angle x-ray scattering from chromatin, isolated by identical procedures but from different species, indicate that fiber diameter and number of nucleosomes per unit length increase with the amount of nucleosome linker DNA. Experiments were conducted at physiological ionic strength to obtain parameters reflecting the structure most likely present in living cells. Guinier analyses were performed on scattering from solutions of soluble chromatin from Necturus maculosus erythrocytes (linker length 48 bp), chicken erythrocytes (linker length 64 bp), and Thyone briareus sperm (linker length 87 bp). The results were extrapolated to infinite dilution to eliminate interparticle contributions to the scattering. Cross-sectional radii of gyration were found to be 10.9 +/- 0.5, 12.1 +/- 0.4, and 15.9 +/- 0.5 nm for Necturus, chicken, and Thyone chromatin, respectively, which are consistent with fiber diameters of 30.8, 34.2, and 45.0 nm. Mass per unit lengths were found to be 6.9 +/- 0.5, 8.3 +/- 0.6, and 11.8 +/- 1.4 nucleosomes per 10 nm for Necturus, chicken, and Thyone chromatin, respectively. The geometrical consequences of the experimental mass per unit lengths and radii of gyration are consistent with a conserved interaction among nucleosomes. Cross-linking agents were found to have little effect on fiber external geometry, but significant effect on internal structure. The absolute values of fiber diameter and mass per unit length, and their dependencies upon linker length agree with the predictions of the double-helical crossed-linker model. A compilation of all published x-ray scattering data from the last decade indicates that the relationship between chromatin structure and linker length is consistent with data obtained by other investigators.     

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.     

Computer simulation of the 30-nanometer chromatin fiber

A new Monte Carlo model for the structure of chromatin is presented here. Based on our previous work on superhelical DNA and polynucleosomes, it reintegrates aspects of the "solenoid" and the "zig-zag" models. The DNA is modeled as a flexible elastic polymer chain, consisting of segments connected by elastic bending, torsional, and stretching springs. The electrostatic interaction between the DNA segments is described by the Debye-Hückel approximation. Nucleosome core particles are represented by oblate ellipsoids; their interaction potential has been parameterized by a comparison with data from liquid crystals of nucleosome solutions. DNA and chromatosomes are linked either at the surface of the chromatosome or through a rigid nucleosome stem. Equilibrium ensembles of 100-nucleosome chains at physiological ionic strength were generated by a Metropolis-Monte Carlo algorithm. For a DNA linked at the nucleosome stem and a nucleosome repeat of 200 bp, the simulated fiber diameter of 32 nm and the mass density of 6.1 nucleosomes per 11 nm fiber length are in excellent agreement with experimental values from the literature. The experimental value of the inclination of DNA and nucleosomes to the fiber axis could also be reproduced. Whereas the linker DNA connects chromatosomes on opposite sides of the fiber, the overall packing of the nucleosomes leads to a helical aspect of the structure. The persistence length of the simulated fibers is 265 nm. For more random fibers where the tilt angles between two nucleosomes are chosen according to a Gaussian distribution along the fiber, the persistence length decreases to 30 nm with increasing width of the distribution, whereas the other observable parameters such as the mass density remain unchanged. Polynucleosomes with repeat lengths of 212 bp also form fibers with the expected experimental properties. Systems with larger repeat length form fibers, but the mass density is significantly lower than the measured value. The theoretical characteristics of a fiber with a repeat length of 192 bp where DNA and nucleosomes are connected at the core particle are in agreement with the experimental values. Systems without a stem and a repeat length of 217 bp do not form fibers.     

Pulling chromatin fibers: computer simulations of direct physical micromanipulations

A low-resolution molecular model, which combines the known mechanical properties of protein-free DNA with the accumulating picture of chromatosome structure, has been developed to account for the stretching of single chromatin fibers by an imposed external force. Force-extension characteristics of sets of chains accumulated by Monte Carlo sampling are consistent with recently observed findings in the non-destructive regime (<20 pN imposed force), where the structure of the chromatosome remains intact. The correspondence between simulation and the relaxation phase of the experiment limits the equilibrium entry-exit angle of linker DNA on the chromatosome to W=50(+/-10) degrees and the effective DNA linker length to L(eff)=40(+/-5) bp. The computed force-extension characteristics are relatively insensitive to other parameters of the model, precluding their accurate estimation. The introduction of an attractive potential between closely spaced nucleosomes reproduces the added initial resistance of single fibers to extension at high salt conditions. The consideration of elastic linkers also improves the fitting of assorted classical measurements of unstressed chromatin structure in solution. The overall picture of chromatin that emerges is an irregular, fluctuating, three-dimensional, zig-zag structure with intact, mechanically stable chromatosome units and deformable linkers. The modeled fiber undergoes large-scale configurational rearrangements without significant perturbation of the constituent chromatosome beads, collapsing into a highly condensed form in response to small (<2kT) inter-nucleosomal attractions.