X-ray diffraction characterization of the dense phases formed by nucleosome core particles

Multiscale modeling reveals the ion-mediated phase separation of nucleosome core particles

Due to the vast length scale inside the cell nucleus, multiscale models are required to understand chromatin folding, structure, and dynamics and how they regulate genomic activities such as DNA transcription, replication, and repair. We study the interactions and structure of condensed phases formed by the universal building block of chromatin, the nucleosome core particle (NCP), using bottom-up multiscale coarse-grained (CG) simulations with a model extracted from all-atom MD simulations. In the presence of the multivalent cations Mg(H2O)62+ or CoHex3+, we analyze the internal structures of the NCP aggregates and the contributions of histone tails and ions to the aggregation patterns. We then derive a "super" coarse-grained (SCG) NCP model to study the macroscopic scale phase separation of NCPs. The SCG simulations show the formation of NCP aggregates with Mg(H2O)62+ concentration-dependent densities and sizes. Variation of the CoHex3+ concentrations results in highly ordered lamellocolumnar and hexagonal columnar phases in agreement with experimental data. The results give detailed insights into nucleosome interactions and for understanding chromatin folding in the cell nucleus.

A Bottom-Up Coarse-Grained Model for Nucleosome-Nucleosome Interactions with Explicit Ions

The nucleosome core particle (NCP) is a large complex of 145-147 base pairs of DNA and eight histone proteins and is the basic building block of chromatin that forms the chromosomes. Here, we develop a coarse-grained (CG) model of the NCP derived through a systematic bottom-up approach based on underlying all-atom MD simulations to compute the necessary CG interactions. The model produces excellent agreement with known structural features of the NCP and gives a realistic description of the nucleosome-nucleosome attraction in the presence of multivalent cations (Mg(H₂O)₆²⁺ or Co(NH₃)₆³⁺) for systems comprising 20 NCPs. The results of the simulations reveal structural details of the NCP-NCP interactions unavailable from experimental approaches, and this model opens the prospect for the rigorous modeling of chromatin fibers.

Linker histone defines structure and self-association behaviour of the 177 bp human chromatosome

Linker histones play essential roles in the regulation and maintenance of the dynamic chromatin structure of higher eukaryotes. The influence of human histone H1.0 on the nucleosome structure and biophysical properties of the resulting chromatosome were investigated and compared with the 177-bp nucleosome using Cryo-EM and SAXS. The 4.5 Å Cryo-EM chromatosome structure showed that the linker histone binds at the nucleosome dyad interacting with both linker DNA arms but in a tilted manner leaning towards one of the linker sides. The chromatosome is laterally compacted and rigid in the dyad and linker DNA area, in comparison with the nucleosome where linker DNA region is more flexible and displays structural variability. In solution, the chromatosomes appear slightly larger than the nucleosomes, with the volume increase compared to the bound linker histone, according to solution SAXS measurements. SAXS X-ray diffraction characterisation of Mg-precipitated samples showed that the different shapes of the 177 chromatosome enabled the formation of a highly ordered lamello-columnar phase when precipitated by Mg²⁺, indicating the influence of linker histone on the nucleosome stacking. The biological significance of linker histone, therefore, may be affected by the change in the polyelectrolyte and DNA conformation properties of the chromatosomes, in comparison to nucleosomes.

Characterization of nucleosome sediments for protein interaction studies by solid-state NMR spectroscopy

Regulation of DNA-templated processes such as gene transcription and DNA repair depend on the interaction of a wide range of proteins with the nucleosome, the fundamental building block of chromatin. Both solution and solid-state NMR spectroscopy have become an attractive approach to study the dynamics and interactions of nucleosomes, despite their high molecular weight of ~ 200 kDa. For solid-state NMR (ssNMR) studies, dilute solutions of nucleosomes are converted to a dense phase by sedimentation or precipitation. Since nucleosomes are known to self-associate, these dense phases may induce extensive interactions between nucleosomes, which could interfere with protein-binding studies. Here, we characterized the packing of nucleosomes in the dense phase created by sedimentation using NMR and small-angle X-ray scattering (SAXS) experiments. We found that nucleosome sediments are gels with variable degrees of solidity, have nucleosome concentration close to that found in crystals, and are stable for weeks under high-speed magic angle spinning (MAS). Furthermore, SAXS data recorded on recovered sediments indicate that there is no pronounced long-range ordering of nucleosomes in the sediment. Finally, we show that the sedimentation approach can also be used to study low-affinity protein interactions with the nucleosome. Together, our results give new insights into the sample characteristics of nucleosome sediments for ssNMR studies and illustrate the broad applicability of sedimentation-based NMR studies.

Nucleosome conformational variability in solution and in interphase nuclei evidenced by cryo-electron microscopy of vitreous sections

In Eukaryotes, DNA is wound around the histone octamer forming the basic chromatin unit, the nucleosome. Atomic structures have been obtained from crystallography and single particle cryo-electron microscopy (cryoEM) of identical engineered particles. But native nucleosomes are dynamical entities with diverse DNA sequence and histone content, and little is known about their conformational variability, especially in the cellular context. Using cryoEM and tomography of vitreous sections we analyse native nucleosomes, both in vitro, using purified particles solubilized at physiologically relevant concentrations (25–50%), and in situ, within interphase nuclei. We visualize individual nucleosomes at a level of detail that allows us to measure the distance between the DNA gyres wrapped around. In concentrated solutions, we demonstrate a salt-dependent transition, with a high salt compact conformation resembling the canonical nucleosome and an open low salt one, closer to nuclear nucleosomes. Although further particle characterization and cartography are needed to understand the relationship between this conformational variability and chromatin functional states, this work opens a route to chromatin exploration in situ.

Frozen-hydrated chromatin from metaphase chromosomes has an interdigitated multilayer structure

Cryo-electron tomography and small-angle X-ray scattering were used to investigate the chromatin folding in metaphase chromosomes. The tomographic 3D reconstructions show that frozen-hydrated chromatin emanated from chromosomes is planar and forms multilayered plates. The layer thickness was measured accounting for the contrast transfer function fringes at the plate edges, yielding a width of ~ 7.5 nm, which is compatible with the dimensions of a monolayer of nucleosomes slightly tilted with respect to the layer surface. Individual nucleosomes are visible decorating distorted plates, but typical plates are very dense and nucleosomes are not identifiable as individual units, indicating that they are tightly packed. Two layers in contact are ~ 13 nm thick, which is thinner than the sum of two independent layers, suggesting that nucleosomes in the layers interdigitate. X-ray scattering of whole chromosomes shows a main scattering peak at ~ 6 nm, which can be correlated with the distance between layers and between interdigitating nucleosomes interacting through their faces. These observations support a model where compact chromosomes are composed of many chromatin layers stacked along the chromosome axis.

Order and interactions in DNA arrays: Multiscale molecular dynamics simulation

While densely packed DNA arrays are known to exhibit hexagonal and orthorhombic local packings, the detailed mechanism governing the associated phase transition remains rather elusive. Furthermore, at high densities the atomistic resolution is paramount to properly account for fine details, encompassing the DNA molecular order, the contingent ordering of counterions and the induced molecular ordering of the bathing solvent, bringing together electrostatic, steric, thermal and direct hydrogen-bonding interactions, resulting in the observed osmotic equation of state. We perform a multiscale simulation of dense DNA arrays by enclosing a set of 16 atomistically resolved DNA molecules within a semi-permeable membrane, allowing the passage of water and salt ions, and thus mimicking the behavior of DNA arrays subjected to external osmotic stress in a bathing solution of monovalent salt and multivalent counterions. By varying the DNA density, local packing symmetry, and counterion type, we obtain osmotic equation of state together with the hexagonal-orthorhombic phase transition, and full structural characterization of the DNA subphase in terms of its positional and angular orientational fluctuations, counterion distributions, and the solvent local dielectric response profile with its order parameters that allow us to identify the hydration force as the primary interaction mechanism at high DNA densities.

The Influence of Ionic Environment and Histone Tails on Columnar Order of Nucleosome Core Particles

The nucleosome core particle (NCP) is the basic building block of chromatin. Nucleosome-nucleosome interactions are instrumental in chromatin compaction, and understanding NCP self-assembly is important for understanding chromatin structure and dynamics. Recombinant NCPs aggregated by multivalent cations form various ordered phases that can be studied by x-ray diffraction (small-angle x-ray scattering). In this work, the effects on the supramolecular structure of aggregated NCPs due to lysine histone H4 tail acetylations, histone H2A mutations (neutralizing the acidic patch of the histone octamer), and the removal of histone tails were investigated. The formation of ordered mainly hexagonal columnar NCP phases is in agreement with earlier studies; however, the highly homogeneous recombinant NCP systems used in this work display a more compact packing. The long-range order of the NCP columnar phase was found to be abolished or reduced by acetylation of the H4 tails, acidic patch neutralization, and removal of the H3 and H2B tails. Loss of nucleosome stacking upon removal of the H3 tails in combination with other tails was observed. In the absence of the H2A tails, the formation of an unknown highly ordered phase was observed.

Multiscale coarse-grained modelling of chromatin components: DNA and the nucleosome

To model large biomolecular systems, such as cell and organelles an atomistic description is not currently achievable and is not generally practical. Therefore, simplified coarse-grained (CG) modelling becomes a necessity. One of the most important cellular components is chromatin, a large DNA-protein complex where DNA is highly compacted. Recent progress in coarse graining modelling of the major chromatin components, double helical DNA and the nucleosome core particle (NCP) is presented. First, general principles and approaches allowing rigorous bottom-to-top generation of interaction potentials in the CG models are presented. Then, recent CG models of DNA are reviewed and their adequacy is benchmarked against experimental data on the salt dependence of DNA flexibility (persistence length). Furthermore, a few recent CG models of the NCP are described and their application for studying salt-dependent NCP-NCP interaction is discussed. An example of a multiscale approach to CG modelling of chromatin is presented where interactions and self-assembly of thousands of NCPs in solution are observed.

Alignment and Graphene-Assisted Decoration of Lyotropic Chromonic Liquid Crystals Containing DNA Origami Nanostructures

Composites of DNA origami nanostructures dispersed in a lyotropic chromonic liquid crystal are studied by polarizing optical microscopy. The homogeneous aqueous dispersions can be uniformly aligned by confinement between two glass substrates, either parallel to the substrates owing to uniaxial rubbing or perpendicular to the substrates using ozonized graphene layers. These opportunities of uniform alignment may pave the way for tailored anisometric plasmonic DNA nanostructures to photonic materials. In addition, a decorated texture with nonuniform orientation is observed on substrates coated with pristine graphene. When the water is allowed to evaporate slowly, microscopic crystal needles appear, which are aligned along the local orientation of the director. This decoration method can be used for studying the local orientational order and the defects in chromonic liquid crystals.

Multiscale modelling of nucleosome core particle aggregation

The nucleosome core particle (NCP) is the basic building block of chromatin. Under the influence of multivalent cations, isolated mononucleosomes exhibit a rich phase behaviour forming various columnar phases with characteristic NCP-NCP stacking. NCP stacking is also a regular element of chromatin structure in vivo. Understanding the mechanism of nucleosome stacking and the conditions leading to self-assembly of NCPs is still incomplete. Due to the complexity of the system and the need to describe electrostatics properly by including the explicit mobile ions, novel modelling approaches based on coarse-grained (CG) methods at the multiscale level becomes a necessity. In this work we present a multiscale CG computer simulation approach to modelling interactions and self-assembly of solutions of NCPs induced by the presence of multivalent cations. Starting from continuum simulations including explicit three-valent cobalt(III)hexammine (CoHex(3+)) counterions and 20 NCPs, based on a previously developed advanced CG NCP model with one bead per amino acid and five beads per two DNA base pair unit (Fan et al 2013 PLoS One 8 e54228), we use the inverse Monte Carlo method to calculate effective interaction potentials for a 'super-CG' NCP model consisting of seven beads for each NCP. These interaction potentials are used in large-scale simulations of up to 5000 NCPs, modelling self-assembly induced by CoHex(3+). The systems of 'super-CG' NCPs form a single large cluster of stacked NCPs without long-range order in agreement with experimental data for NCPs precipitated by the three-valent polyamine, spermidine(3+).

An advanced coarse-grained nucleosome core particle model for computer simulations of nucleosome-nucleosome interactions under varying ionic conditions

In the eukaryotic cell nucleus, DNA exists as chromatin, a compact but dynamic complex with histone proteins. The first level of DNA organization is the linear array of nucleosome core particles (NCPs). The NCP is a well-defined complex of 147 bp DNA with an octamer of histones. Interactions between NCPs are of paramount importance for higher levels of chromatin compaction. The polyelectrolyte nature of the NCP implies that nucleosome-nucleosome interactions must exhibit a great influence from both the ionic environment as well as the positively charged and highly flexible N-terminal histone tails, protruding out from the NCP. The large size of the system precludes a modelling analysis of chromatin at an all-atom level and calls for coarse-grained approximations. Here, a model of the NCP that include the globular histone core and the flexible histone tails described by one particle per each amino acid and taking into account their net charge is proposed. DNA wrapped around the histone core was approximated at the level of two base pairs represented by one bead (bases and sugar) plus four beads of charged phosphate groups. Computer simulations, using a Langevin thermostat, in a dielectric continuum with explicit monovalent (K(+)), divalent (Mg(2+)) or trivalent (Co(NH(3))(6) (3+)) cations were performed for systems with one or ten NCPs. Increase of the counterion charge results in a switch from repulsive NCP-NCP interaction in the presence of K(+), to partial aggregation with Mg(2+) and to strong mutual attraction of all 10 NCPs in the presence of CoHex(3+). The new model reproduced experimental results and the structure of the NCP-NCP contacts is in agreement with available data. Cation screening, ion-ion correlations and tail bridging contribute to the NCP-NCP attraction and the new NCP model accounts for these interactions.

Influence of histone tails and H4 tail acetylations on nucleosome-nucleosome interactions

Nucleosome-nucleosome interaction plays a fundamental role in chromatin folding and self-association. The cation-induced condensation of nucleosome core particles (NCPs) displays properties similar to those of chromatin fibers, with important contributions from the N-terminal histone tails. We study the self-association induced by addition of cations [Mg(2+), Ca(2+), cobalt(III)hexammine(3+), spermidine(3+) and spermine(4)(+)] for NCPs reconstituted with wild-type unmodified histones and with globular tailless histones and for NCPs with the H4 histone tail having lysine (K) acetylations or lysine-to-glutamine mutations at positions K5, K8, K12 and K16. In addition, the histone construct with the single H4K16 acetylation was investigated. Acetylated histones were prepared by a semisynthetic native chemical ligation method. The aggregation behavior of NCPs shows a general cation-dependent behavior similar to that of the self-association of nucleosome arrays. Unlike nucleosome array self-association, NCP aggregation is sensitive to position and nature of the H4 tail modification. The tetra-acetylation in the H4 tail significantly weakens the nucleosome-nucleosome interaction, while the H4 K→Q tetra-mutation displays a more modest effect. The single H4K16 acetylation also weakens the self-association of NCPs, which reflects the specific role of H4K16 in the nucleosome-nucleosome stacking. Tailless NCPs can aggregate in the presence of oligocations, which indicates that attraction also occurs by tail-independent nucleosome-nucleosome stacking and DNA-DNA attraction in the presence of cations. The experimental data were compared with the results of coarse-grained computer modeling for NCP solutions with explicit presence of mobile ions.

Cation-induced polyelectrolyte-polyelectrolyte attraction in solutions of DNA and nucleosome core particles

The paper reviews our current studies on the experimentally induced cation compaction and aggregation in solutions of DNA and nucleosome core particles and the theoretical modelling of these processes using coarse-grained continuum models with explicit mobile ions and with all-atom molecular dynamics (MD) simulations. Recent experimental results on DNA condensation by cationic oligopeptides and the effects of added salt are presented. The results of MD simulations modelling DNA-DNA attraction due to the presence of multivalent ions including the polyamine spermidine and fragments of histone tails, which exhibit bridging between adjacent DNA molecules, are discussed. Experimental data on NCP aggregation, using recombinantly prepared systems are summarized. Literature data and our results of studying of the NCP solutions are compared with predictions of coarse-grained MD simulations, including the important ion correlation as well as bridging mechanisms. The importance of the results to chromatin folding and aggregation is discussed.

Computer modeling reveals that modifications of the histone tail charges define salt-dependent interaction of the nucleosome core particles

Coarse-grained Langevin molecular dynamics computer simulations were conducted for systems that mimic solutions of nucleosome core particles (NCPs). The NCP was modeled as a negatively charged spherical particle representing the complex of DNA and the globular part of the histones combined with attached strings of connected charged beads modeling the histone tails. The size, charge, and distribution of the tails relative to the core were built to match real NCPs. Three models of NCPs were constructed to represent different extents of covalent modification on the histone tails: (nonmodified) recombinant (rNCP), acetylated (aNCP), and acetylated and phosphorylated (paNCP). The simulation cell contained 10 NCPs in a dielectric continuum with explicit mobile counterions and added salt. The NCP-NCP interaction is decisively dependent on the modification state of the histone tails and on salt conditions. Increasing the monovalent salt concentration (KCl) from salt-free to physiological concentration leads to NCP aggregation in solution for rNCP, whereas NCP associates are observed only occasionally in the system of aNCPs. In the presence of divalent salt (Mg(2+)), rNCPs form dense stable aggregates, whereas aNCPs form aggregates less frequently. Aggregates are formed via histone-tail bridging and accumulation of counterions in the regions of NCP-NCP contacts. The paNCPs do not show NCP-NCP interaction upon addition of KCl or in the presence of Mg(2+). Simulations for systems with a gradual substitution of K(+) for Mg(2+), to mimic the Mg(2+) titration of an NCP solution, were performed. The rNCP system showed stronger aggregation that occurred at lower concentrations of added Mg(2+), compared to the aNCP system. Additional molecular dynamics simulations performed with a single NCP in the simulation cell showed that detachment of the tails from the NCP core was modest under a wide range of salt concentrations. This implies that salt-induced tail dissociation of the histone tails from the globular NCP is not in itself a major factor in NCP-NCP aggregation. The approximation of coarse-graining, with respect to the description of the NCP as a sphere with uniform charge distribution, was tested in control simulations. A more detailed description of the NCP did not change the main features of the results. Overall, the results of this work are in agreement with experimental data reported for NCP solutions and for chromatin arrays.

The folding and unfolding of eukaryotic chromatin

In vivo, chromatin exists as fibres with differing degrees of compaction. We argue here that the packing density of the chromatin fibre is an important parameter, such that fibres with six nucleosomes/11 nm are enriched in 'euchromatin' while more highly compacted forms with higher packing densities correspond to some heterochromatic regions. The fibre forms differ in the extent of nucleosome stacking-in the '30 nm' fibre stacking is suboptimal while in 'heterochromatic' fibres optimal stacking allows a greater compaction. One factor affecting the choice of different endpoints in fibre formation depends on the homogeneity and optimisation of linker length within a nucleosomal array. The '30 nm' fibre can accommodate some variation in linker length while formation of the more compact forms requires that linker lengths be homogeneous and optimal. In vivo, chromatin remodelling machines and histone tail modifications would mediate and regulate this optimisation.

Computational modeling of the chromatin fiber

The packing of the genomic DNA in the living cell is essential for its biological function. While individual aspects of the genome architecture, such as DNA and nucleosome structure or the arrangement of chromosome territories are well studied, much information is missing for a unified description of cellular DNA at all its structural levels. Computer modeling can contribute to such a description. We present here some typical approaches to models of the chromatin fiber, including different amounts of detail in the description of the local nucleosome structure. The main results from our simulations are that the physical properties of the chromatin fiber can be well described by a simplified model consisting of cylinder-like nucleosomes connected by flexible DNA segments, with a geometry determined by the bending and twisting angles between nucleosomes. Randomness in the local geometry - such as random absence of linker histone H1 - leads to a dramatic increase in the chromatin fiber flexibility. Furthermore, we show that chromatin is much more flexible to bending than to stretching, and that the structure of the chromatin fiber favors the formation of sharp bends.

Are liquid crystalline properties of nucleosomes involved in chromosome structure and dynamics?

Nucleosome core particles correspond to the structural units of eukaryotic chromatin. They are charged colloids, 101 Angstrom in diameter and 55 Angstrom in length, formed by the coiling of a 146/147 bp DNA fragment (50 nm) around the histone protein octamer. Solutions of these particles can be concentrated, under osmotic pressure, up to the concentrations found in the nuclei of living cells. In the presence of monovalent cations (Na(+)), nucleosomes self-assemble into crystalline or liquid crystalline phases. A lamello-columnar phase is observed at 'low salt' concentrations, while a two-dimensional hexagonal phase and a three-dimensional quasi-hexagonal phase form at 'high salt' concentrations. We followed the formation of these phases from the dilute isotropic solutions to the ordered phases by combining cryoelectron microscopy and X-ray diffraction analyses. The phase diagram is presented as a function of the monovalent salt concentration and applied osmotic pressure. An alternative method to condense nucleosomes is to induce their aggregation upon addition of divalent or multivalent cations (Mg(2+), spermidine(3+) and spermine(4+)). Ordered phases are also found in the aggregates. We also discuss whether these condensed phases of nucleosomes may be relevant from a biological point of view.

Structural plasticity of single chromatin fibers revealed by torsional manipulation

Magnetic tweezers were used to study the mechanical response under torsion of single nucleosome arrays reconstituted on tandem repeats of 5S positioning sequences. Regular arrays are extremely resilient and can reversibly accommodate a large amount of supercoiling without much change in length. This behavior is quantitatively described by a molecular model of the chromatin three-dimensional architecture. In this model, we assume the existence of a dynamic equilibrium between three conformations of the nucleosome, corresponding to different crossing statuses of the entry/exit DNAs (positive, null or negative, respectively). Torsional strain displaces that equilibrium, leading to an extensive reorganization of the fiber's architecture. The model explains a number of long-standing topological questions regarding DNA in chromatin and may provide the basis to better understand the dynamic binding of chromatin-associated proteins.Note: In the supplementary information initially published online to accompany this article, Supplementary Figure 2 was mistakenly replaced by Supplementary Equation 2. The error has been corrected online.

Polymer chain models of DNA and chromatin

Many properties of the genome in the cell nucleus can be understood by modeling DNA and chromatin as a flexible polymer chain. This article introduces into current models for such a coarse-grained description and reviews some recent results from our own group. Examples given are the unrolling of DNA from the histone core and the response of the 30 nm chromatin fiber to mechanical stretching.

Systematic search for compact structures of telomeric nucleosomes

Telomeres are structures functionally and structurally distinct from bulk chromatin. They are constituted of highly conserved 5-7 bp tandemly repeated units, organized into nucleosomes with short linkers, whereas the knowledge of the linker histone role in telomeric chromatin is still fragmentary. Experimental evidence suggests the structural organization of telomeric nucleosomes is different from that of the bulk chromatin. This work presents a systematic search of the telomeric nucleosome arrangements. A low-resolution molecular model was used to evaluate the relative nucleosome packing energy. Structures with favorable energy were found, reducing the possible telomeric chromatin conformations to two different three-dimensional folds.

Phase diagram of nucleosome core particles

We present a phase diagram of the nucleosome core particle (NCP) as a function of the monovalent salt concentration and applied osmotic pressure. Above a critical pressure, NCPs stack on top of each other to form columns that further organize into multiple columnar phases. An isotropic (and in some cases a nematic) phase of columns is observed in the moderate pressure range. Under higher pressure conditions, a lamello-columnar phase and an inverse hexagonal phase form under low salt conditions, whereas a 2D hexagonal phase or a 3D orthorhombic phase is found at higher salt concentration. For intermediate salt concentrations, microphase separation occurs. The richness of the phase diagram originates from the heterogeneous distribution of charges at the surface of the NCP, which makes the particles extremely sensitive to small ionic variations of their environment, with consequences on their interactions and supramolecular organization. We discuss how the polymorphism of NCP supramolecular organization may be involved in chromatin changes in the cellular context.