Electrostatic origin of salt-induced nucleosome array compaction

Physical modeling of nucleosome clustering in euchromatin resulting from interactions between epigenetic reader proteins

Euchromatin is an accessible phase of genetic material containing genes that encode proteins with increased expression levels. The structure of euchromatin in vitro has been described as a 30-nm fiber formed from ordered nucleosome arrays. However, recent advances in microscopy have revealed an in vivo euchromatin architecture that is much more disordered, characterized by variable-length linker DNA and sporadic nucleosome clusters. In this work, we develop a theoretical model to elucidate factors contributing to the disordered in vivo architecture of euchromatin. We begin by developing a 1D model of nucleosome positioning that captures the interactions between bound epigenetic reader proteins to predict the distribution of DNA linker lengths between adjacent nucleosomes. We then use the predicted linker lengths to construct 3D chromatin configurations consistent with the physical properties of DNA within the nucleosome array, and we evaluate the distribution of nucleosome cluster sizes in those configurations. Our model reproduces experimental cluster-size distributions, which are dramatically influenced by the local pattern of epigenetic marks and the concentration of reader proteins. Based on our model, we attribute the disordered arrangement of euchromatin to the heterogeneous binding of reader proteins and subsequent short-range interactions between bound reader proteins on adjacent nucleosomes. By replicating experimental results with our physics-based model, we propose a mechanism for euchromatin organization in the nucleus that impacts gene regulation and the maintenance of epigenetic marks.

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.

The shelterin component TRF2 mediates columnar stacking of human telomeric chromatin

Telomere repeat binding factor 2 (TRF2) is an essential component of the telomeres and also plays an important role in a number of other non-telomeric processes. Detailed knowledge of the binding and interaction of TRF2 with telomeric nucleosomes is limited. Here, we study the binding of TRF2 to in vitro-reconstituted kilobasepair-long human telomeric chromatin fibres using electron microscopy, single-molecule force spectroscopy and analytical ultracentrifugation sedimentation velocity. Our electron microscopy results revealed that full-length and N-terminally truncated TRF2 promote the formation of a columnar structure of the fibres with an average width and compaction larger than that induced by the addition of Mg2+, in agreement with the in vivo observations. Single-molecule force spectroscopy showed that TRF2 increases the mechanical and thermodynamic stability of the telomeric fibres when stretched with magnetic tweezers. This was in contrast to the result for fibres reconstituted on the 'Widom 601' high-affinity nucleosome positioning sequence, where minor effects on fibre stability were observed. Overall, TRF2 binding induces and stabilises columnar fibres, which may play an important role in telomere maintenance.

Nanopore Sensing of DNA-Histone Complexes on Nucleosome Arrays

The location of nucleosomes in DNA and their structural stability are critical in regulating DNA compaction, site accessibility, and epigenetic gene regulation. Here, we combine the nanopore platform-based fast and label-free single-molecule detection technique with a voltage-dependent force rupture assay to detect distinct structures on nucleosomal arrays and then to induce breakdown of individual nucleosome complexes. Specifically, we demonstrate direct measurement of distinct nucleosome structures present on individual 12-mer arrays. A detailed event analysis showed that nucleosomes are present as a combination of complete and partial structures, during translocation through the pore. By comparing with the voltage-dependent translocation of the mononucleosomes, we find that the partial nucleosomes result from voltage-dependent structural disintegration of nucleosomes. High signal-to-noise detection of heterogeneous levels in translocation of 12-mer array molecules quantifies the heterogeneity and nucleosomal substructure sizes on the arrays. These results facilitate the understanding of electrostatic interactions responsible for the integrity of the nucleosome structure and possible mechanisms of its unraveling by chromatin remodeling enzymes. This study also has potential applications in chromatin profiling.

Chromatin Liquid-Liquid Phase Separation (LLPS) Is Regulated by Ionic Conditions and Fiber Length

The dynamic regulation of the physical states of chromatin in the cell nucleus is crucial for maintaining cellular homeostasis. Chromatin can exist in solid- or liquid-like forms depending on the surrounding ions, binding proteins, post-translational modifications and many other factors. Several recent studies suggested that chromatin undergoes liquid-liquid phase separation (LLPS) in vitro and also in vivo; yet, controversial conclusions about the nature of chromatin LLPS were also observed from the in vitro studies. These inconsistencies are partially due to deviations in the in vitro buffer conditions that induce the condensation/aggregation of chromatin as well as to differences in chromatin (nucleosome array) constructs used in the studies. In this work, we present a detailed characterization of the effects of K⁺, Mg²⁺ and nucleosome fiber length on the physical state and property of reconstituted nucleosome arrays. LLPS was generally observed for shorter nucleosome arrays (15-197-601, reconstituted from 15 repeats of the Widom 601 DNA with 197 bp nucleosome repeat length) at physiological ion concentrations. In contrast, gel- or solid-like condensates were detected for the considerably longer 62-202-601 and lambda DNA (~48.5 kbp) nucleosome arrays under the same conditions. In addition, we demonstrated that the presence of reduced BSA and acetate buffer is not essential for the chromatin LLPS process. Overall, this study provides a comprehensive understanding of several factors regarding chromatin physical states and sheds light on the mechanism and biological relevance of chromatin phase separation in vivo.

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.

Potential of Mean Force for DNA Wrapping Around a Cationic Nanoparticle

Sharp bending and wrapping of DNA around proteins and nanoparticles (NPs) has been of extensive research interest. Here, we present the potential of mean force (PMF) for wrapping a DNA double helix around a cationic NP using coarse-grained models of a double-stranded DNA and a cationic NP. Starting from a NP wrapped around by DNA, the PMF was calculated along the distance between the center of the NP and one end of the DNA molecule. A relationship between the distance and the extent of DNA wrapping is used to calculate the PMF as a function of DNA wrapping around a NP. In particular, the PMF was compared for two DNA sequences of (AT)₂₅/(AT)₂₅ and (AC)₂₅/(GT)₂₅, for which the persistence lengths are different by ∼10 nm. The simulation results provide solid evidence of the thermodynamic preference for complex formation of a cationic NP with more flexible DNA over the less flexible DNA. Furthermore, we estimated the elastic energy of DNA bending, which was in good order-of-magnitude agreement with the theoretical prediction of elastic rods. This work suggests that the variation of sequence-dependent DNA flexibility can be utilized in DNA nanotechnologies, in which the position and dynamics of NPs are regulated on large-scale DNA structures, or the structural transformation of DNA is triggered by the sequence-dependent binding of NPs.

Effect of DNA Flexibility on Complex Formation of a Cationic Nanoparticle with Double-Stranded DNA

We present extensive molecular dynamics simulations of a cationic nanoparticle and a double-stranded DNA molecule to discuss the effect of DNA flexibility on the complex formation of a cationic nanoparticle with double-stranded DNA. Martini coarse-grained models were employed to describe double-stranded DNA molecules with two different flexibilities and cationic nanoparticles with three different electric charges. As the electric charge of a cationic nanoparticle increases, the degree of DNA bending increases, eventually leading to the wrapping of DNA around the nanoparticle at high electric charges. However, a small increase in the persistence length of DNA by 10 nm requires a cationic nanoparticle with a markedly increased electric charge to bend and wrap DNA around. Thus, a more flexible DNA molecule bends and wraps around a cationic nanoparticle with an intermediate electric charge, whereas a less flexible DNA molecule binds to a nanoparticle with the same electric charge without notable bending. This work provides solid evidence that a small difference in DNA flexibility (as small as 10 nm in persistence length) has a substantial influence on the complex formation of DNA with proteins from a biological perspective and suggests that the variation of sequence-dependent DNA flexibility can be utilized in DNA nanotechnology as a new tool to manipulate the structure of DNA molecules mediated by nanoparticle binding.

Is the H4 histone tail intrinsically disordered or intrinsically multifunctional?

The structural versatility of histone tails is one of the key elements in the organisation of chromatin, which allows for the compact storage of genomic information. However, this structural diversity also complicates experimental and computational studies. Here, the potential and free energy landscape for the isolated and bound H4 histone tail are explored. The landscapes exhibit a set of distinct structural ensembles separated by high energy barriers, with little difference between isolated and bound tails. This consistency is a desirable feature that facilitates the formation of transient interactions, which are required for the liquid-like chromatin organisation. The existence of multiple, distinct structures on a multifunnel energy landscape is likely to be associated with multifunctionality, i.e. a set of evolved, distinct functions. Contrasting it with previously reported results for other disordered peptides, this type of landscape may be associated with a conformational selection based binding mechanism. Given the similarity to other systems exhibiting similar multifunnel energy landscapes, the disorder in histone tails might be better described in context of multifunctionality.

Beyond the Nucleosome: Nucleosome-Protein Interactions and Higher Order Chromatin Structure

The regulation of chromatin biology ultimately depends on the manipulation of its smallest subunit, the nucleosome. The proteins that bind and operate on the nucleosome do so, while their substrate is part of a polymer embedded in the dense nuclear environment. Their molecular interactions must in some way be tuned to deal with this complexity. Due to the rapid increase in the number of high-resolution structures of nucleosome-protein complexes and the increasing understanding of the cellular chromatin structure, it is starting to become clearer how chromatin factors operate in this complex environment. In this review, we analyze the current literature on the interplay between nucleosome-protein interactions and higher-order chromatin structure. We examine in what way nucleosomes-protein interactions can affect and can be affected by chromatin organization at the oligonucleosomal level. In addition, we review the characteristics of nucleosome-protein interactions that can cause phase separation of chromatin. Throughout, we hope to illustrate the exciting challenges in characterizing nucleosome-protein interactions beyond the nucleosome.

Global histone protein surface accessibility in yeast indicates a uniformly loosely packed genome with canonical nucleosomes

Background

The vast majority of methods available to characterize genome-wide chromatin structure exploit differences in DNA accessibility to nucleases or chemical crosslinking. We developed a novel method to gauge genome-wide accessibility of histone protein surfaces within nucleosomes by assessing reactivity of engineered cysteine residues with a thiol-specific reagent, biotin-maleimide (BM).

Results

Yeast nuclei were obtained from cells expressing the histone mutant H2B S116C, in which a cysteine resides near the center of the external flat protein surface of the nucleosome. BM modification revealed that nucleosomes are generally equivalently accessible throughout the S. cerevisiae genome, including heterochromatic regions, suggesting limited, higher-order chromatin structures in which this surface is obstructed by tight nucleosome packing. However, we find that nucleosomes within 500 bp of transcription start sites exhibit the greatest range of accessibility, which correlates with the density of chromatin remodelers. Interestingly, accessibility is not well correlated with RNA polymerase density and thus the level of gene expression. We also investigated the accessibility of cysteine mutations designed to detect exposure of histone surfaces internal to the nucleosome thought to be accessible in actively transcribed genes: H3 102, is at the H2A–H2B dimer/H3–H4 tetramer interface, and H3 A110C, resides at the H3–H3 interface. However, in contrast to the external surface site, we find that neither of these internal sites were found to be appreciably exposed.

Conclusions

Overall, our finding that nucleosomes surfaces within S. cerevisiae chromatin are equivalently accessible genome-wide is consistent with a globally uncompacted chromatin structure lacking substantial higher-order organization. However, we find modest differences in accessibility that correlate with chromatin remodelers but not transcription, suggesting chromatin poised for transcription is more accessible than actively transcribed or intergenic regions. In contrast, we find that two internal sites remain inaccessible, suggesting that such non-canonical nucleosome species generated during transcription are rapidly and efficiently converted to canonical nucleosome structure and thus not widely present in native chromatin.

Abiotic stress-mediated modulation of the chromatin landscape in Arabidopsis thaliana

Limited information is available on abiotic stress-mediated alterations of chromatin conformation influencing gene expression in plants. In order to characterize the effect of abiotic stresses on changes in chromatin conformation, we employed FAIRE-seq (formaldehyde-assisted isolation of regulatory element sequencing) and DNase-seq to isolate accessible regions of chromatin from Arabidopsis thaliana seedlings exposed to either heat, cold, salt, or drought stress. Approximately 25% of regions in the Arabidopsis genome were captured as open chromatin, the majority of which included promoters and exons. A large proportion of chromatin regions apparently did not change their conformation in response to any of the four stresses. Digital footprints present within these regions had differential enrichment of motifs for binding of 43 different transcription factors. Further, in contrast to drought and salt stress, both high and low temperature treatments resulted in increased accessibility of the chromatin. Also, pseudogenes attained increased chromatin accessibility in response to cold and drought stresses. The highly accessible and inaccessible chromatin regions of seedlings exposed to drought stress correlated with the Ser/Thr protein kinases (MLK1 and MLK2)-mediated reduction and increase in H3 phosphorylation (H3T3Ph), respectively. The presented results provide a deeper understanding of abiotic stress-mediated chromatin modulation in plants.

The effect of linker DNA on the structure and interaction of nucleosome core particles

In eukaryotes, the compaction of chromatin fibers composed of nucleosome core particles (NCPs) connected by a linker DNA into chromosomes is highly efficient; however, the underlying folding mechanisms remain elusive. We used small angle X-ray scattering (SAXS) to investigate the influence of linker DNA length on the local structure and the interparticle interactions of the NCPs. In the presence of the linker DNA of 30 bp or less in length, the results suggest partial unwrapping of nucleosomal DNA on the NCP irrespective of the linker DNA length. Moreover, the presence of 15 bp linker DNA alleviated the electrostatic repulsion between the NCPs and prevented the formation of an ordered columnar hexagonal phase, demonstrating that the linker DNA plays an active role in chromatin folding.

Structure of an H1-Bound 6-Nucleosome Array Reveals an Untwisted Two-Start Chromatin Fiber Conformation

Chromatin adopts a diversity of regular and irregular fiber structures in vitro and in vivo. However, how an array of nucleosomes folds into and switches between different fiber conformations is poorly understood. We report the 9.7 Å resolution crystal structure of a 6-nucleosome array bound to linker histone H1 determined under ionic conditions that favor incomplete chromatin condensation. The structure reveals a flat two-start helix with uniform nucleosomal stacking interfaces and a nucleosome packing density that is only half that of a twisted 30-nm fiber. Hydroxyl radical footprinting indicates that H1 binds the array in an on-dyad configuration resembling that observed for mononucleosomes. Biophysical, cryo-EM, and crosslinking data validate the crystal structure and reveal that a minor change in ionic environment shifts the conformational landscape to a more compact, twisted form. These findings provide insights into the structural plasticity of chromatin and suggest a possible assembly pathway for a 30-nm fiber.

Compaction of Single-Molecule Megabase-Long Chromatin under the Influence of Macromolecular Crowding

The megabase-sized length of chromatin is highly relevant to the state of chromatin in vivo, where it is subject to a highly crowded environment and is organized in topologically associating domains of similar dimension. We developed an in vitro experimental chromatin model system reconstituted from T4 DNA (approximately 166 kbp) and histone octamers and studied the monomolecular compaction of this megabase-sized chromatin fiber under the influence of macromolecular crowding. We used single-molecule fluorescence microscopy and observed compaction in aqueous solutions containing poly(ethylene glycol) in the presence of monovalent (Na+ and K+) and divalent (Mg2+) cations. Both DNA and chromatin demonstrated compaction under comparable conditions in the presence of poly(ethylene glycol) and Na+ or Mg2+ salt. However, the mechanism of the compaction changed from a first-order phase transition for DNA to a continuous folding for megabase-sized chromatin fibers. A more efficient and pronounced chromatin compaction was observed in the presence of Na+ compared to K+. A flow-stretching technique to unfold DNA and chromatin coils was used to gain further insight into the morphology of partially folded chromatin fibers. The results revealed a distribution of partially folded chromatin fibers. This variability is likely the result of the heterogeneous distribution of nucleosomes on the DNA chain. The packaging of DNA in the form of chromatin in the crowded nuclear environment appears essential to ensure gradual conformational changes of DNA.

Site-Specific Studies of Nucleosome Interactions by Solid-State NMR Spectroscopy

Chromatin function depends on a dense network of interactions between nucleosomes and a wide range of proteins. A detailed description of these protein-nucleosome interactions is required to reach a full molecular understanding of chromatin function in both genetics and epigenetics. Herein, we show that the structure, dynamics, and interactions of nucleosomes can be interrogated in a residue-specific manner by using state-of-the-art solid-state NMR spectroscopy. Using sedimented nucleosomes, high-resolution spectra were obtained for both flexible histone tails and the non-mobile histone core. Through co-sedimentation of a nucleosome-binding peptide, we demonstrate that protein-binding sites on the nucleosome surface can be determined. We believe that this approach holds great promise as it is generally applicable, extendable to include the structure and dynamics of the bound proteins, and scalable to interactions of proteins with higher-order chromatin structures, including isolated and cellular chromatin.

Nucleosome-level 3D organization of the genome

Nucleosomes are the unitary structures of chromosome folding, and their arrangements are intimately coupled to the regulation of genome activities. Conventionally, structural analyses using electron microscopy and X-ray crystallography have been used to study such spatial nucleosome arrangements. In contrast, recent improvements in the resolution of sequencing-based methods allowed investigation of nucleosome arrangements separately at each genomic locus, enabling exploration of gene-dependent regulation mechanisms. Here, we review recent studies on nucleosome folding in chromosomes from these two methodological perspectives: conventional structural analyses and DNA sequencing, and discuss their implications for future research.

Modulation of nucleosomal DNA accessibility via charge-altering post-translational modifications in histone core

BACKGROUND

Controlled modulation of nucleosomal DNA accessibility via post-translational modifications (PTM) is a critical component to many cellular functions. Charge-altering PTMs in the globular histone core-including acetylation, phosphorylation, crotonylation, propionylation, butyrylation, formylation, and citrullination-can alter the strong electrostatic interactions between the oppositely charged nucleosomal DNA and the histone proteins and thus modulate accessibility of the nucleosomal DNA, affecting processes that depend on access to the genetic information, such as transcription. However, direct experimental investigation of the effects of these PTMs is very difficult. Theoretical models can rationalize existing observations, suggest working hypotheses for future experiments, and provide a unifying framework for connecting PTMs with the observed effects.

RESULTS

A physics-based framework is proposed that predicts the effect of charge-altering PTMs in the histone core, quantitatively for several types of lysine charge-neutralizing PTMs including acetylation, and qualitatively for all phosphorylations, on the nucleosome stability and subsequent changes in DNA accessibility, making a connection to resulting biological phenotypes. The framework takes into account multiple partially assembled states of the nucleosome at the atomic resolution. The framework is validated against experimentally known nucleosome stability changes due to the acetylation of specific lysines, and their effect on transcription. The predicted effect of charge-altering PTMs on DNA accessibility can vary dramatically, from virtually none to a strong, region-dependent increase in accessibility of the nucleosomal DNA; in some cases, e.g., H4K44, H2AK75, and H2BK57, the effect is significantly stronger than that of the extensively studied acetylation sites such H3K56, H3K115 or H3K122. Proximity to the DNA is suggestive of the strength of the PTM effect, but there are many exceptions. For the vast majority of charge-altering PTMs, the predicted increase in the DNA accessibility should be large enough to result in a measurable modulation of transcription. However, a few possible PTMs, such as acetylation of H4K77, counterintuitively decrease the DNA accessibility, suggestive of the repressed chromatin. A structural explanation for the phenomenon is provided. For the majority of charge-altering PTMs, the effect on DNA accessibility is simply additive (noncooperative), but there are exceptions, e.g., simultaneous acetylation of H4K79 and H3K122, where the combined effect is amplified. The amplification is a direct consequence of the nucleosome-DNA complex having more than two structural states. The effect of individual PTMs is classified based on changes in the accessibility of various regions throughout the nucleosomal DNA. The PTM's resulting imprint on the DNA accessibility, "PTMprint," is used to predict effects of many yet unexplored PTMs. For example, acetylation of H4K44 yields a PTMprint similar to the PTMprint of H3K56, and thus acetylation of H4K44 is predicted to lead to a wide range of strong biological effects.

CONCLUSION

Charge-altering post-translational modifications in the relatively unexplored globular histone core may provide a precision mechanism for controlling accessibility to the nucleosomal DNA.

Single-molecule compaction of megabase-long chromatin molecules by multivalent cations

To gain insight into the conformational properties and compaction of megabase-long chromatin molecules, we reconstituted chromatin from T4 phage DNA (165 kb) and recombinant human histone octamers (HO). The unimolecular compaction, induced by divalent Mg2+ or tetravalent spermine4+ cations, studied by single-molecule fluorescence microscopy (FM) and dynamic light scattering (DLS) techniques, resulted in the formation of 250-400 nm chromatin condensates. The compaction on this scale of DNA size is comparable to that of chromatin topologically associated domains (TAD) in vivo. Variation of HO loading revealed a number of unique features related to the efficiency of chromatin compaction by multivalent cations, the mechanism of compaction, and the character of partly compact chromatin structures. The observations may be relevant for how DNA accessibility in chromatin is maintained. Compaction of saturated chromatin, in turn, is accompanied by an intra-chain segregation at the level of single chromatin molecules, suggesting an intriguing scenario of selective activation/deactivation of DNA as a result of chromatin fiber heterogeneity due to the nucleosome positioning. We suggest that this chromatin, reconstituted on megabase-long DNA because of its large size, is a useful model of eukaryotic chromatin.

Effects of charge-modifying mutations in histone H2A α3-domain on nucleosome stability assessed by single-pair FRET and MD simulations

Nucleosomes are important for chromatin compaction and gene regulation; their integrity depends crucially on the structural properties of the histone tails. Recent all-atom molecular dynamics simulations revealed that removal of the N-terminal tails of histone H3, known to destabilize nucleosomes, causes a rearrangement of two arginines of histone H2A, namely R81 and R88 by altering the electrostatic environment of the H2A α3 domain. Whether this rearrangement is the cause or the effect of decreased stability, is unclear. Here, we emulate the altered electrostatic environment that was found after H3 tail clipping through charge-modifying mutations to decouple its impact on intranucleosomal interactions from that of the histone tails. Förster resonance energy transfer experiments on recombinant nucleosomes and all-atom molecular dynamics simulations reveal a compensatory role of those amino acids in nucleosome stability. The simulations indicate a weakened interface between H2A-H2B dimers and the (H3-H4)₂ tetramer, as well as between dimers and DNA. These findings agree with the experimental observations of position and charge dependent decreased nucleosome stability induced by the introduced mutations. This work highlights the importance of the H2A α3 domain and suggests allosteric effects between this domain and the outer DNA gyre as well as the H3 N-terminal tail.

Basic surface features of nuclear FKBPs facilitate chromatin binding

The nucleoplasmin family of histone chaperones is identified by a pentamer-forming domain and multiple acidic tracts that mediate histone binding and chaperone activity. Within this family, a novel domain organization was recently discovered that consists of an N-terminal nucleoplasmin-like (NPL) domain and a C-terminal FKBP peptidyl-proline isomerase domain. Saccharomyces cerevisiae Fpr4 is one such protein. Here we report that in addition to its known histone prolyl isomerase activities, the Fpr4 FKBP domain binds to nucleosomes and nucleosome arrays in vitro. This ability is mediated by a collection of basic patches that enable the enzyme to stably associate with linker DNA. The interaction of the Fpr4 FKBP with recombinant chromatin complexes condenses nucleosome arrays independently of its catalytic activity. Based on phylogenetic comparisons we propose that the chromatin binding ability of ‘basic’ FKBPs is shared amongst related orthologues present in fungi, plants, and insects. Thus, a subclass of FKBP prolyl isomerase enzymes is recruited to linker regions of chromatin.

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.

Confinement-driven organization of a histone-complexed DNA molecule in a dense array of nanoposts

The first step in the controlled storage of lengthy DNA molecules is to keep DNA molecules separated while integrated in micrometer-sized space. Herein, we present hybrid Monte Carlo simulations of a histone-complexed DNA (hcDNA) molecule confined in a dense array of nanoposts. Depending on the nanopost dimension, a single, 8.7 kilobase pair hcDNA molecule was either localized and elongated in a single inter-post space surrounded by four nanoposts or spread over several inter-post spaces through passages between two neighboring nanoposts. The conformational change of a hcDNA molecule is interpreted in terms of competitive effects of confinements in the inter-post and passage spaces. We propose that, by elaborately designing nanopost arrays, the competitive confinement effects can be adjusted such that each hcDNA molecule is localized in a single inter-post space, and thereby multiple hcDNA molecules can be physically separated from each other while stored together in the nanopost array.

Biophysical studies of cholesterol effects on chromatin

Changes in chromatin structure regulate gene expression and genome maintenance. Molecules that bind to the nucleosome, the complex of DNA and histone proteins, are key modulators of chromatin structure. Previous work indicated that cholesterol, a ubiquitous cellular lipid, may bind to chromatin in vivo, suggesting a potential function for lipids in modulating chromatin architecture. However, the molecular mechanisms of cholesterol’s action on chromatin structure have remained unclear. Here, we explored the biophysical impact of cholesterol on nucleosome and chromatin fibers reconstituted in vitro and characterized in silico the cholesterol binding to the nucleosome. Our findings support that cholesterol assists 10 and 30 nm chromatin formation and induces folding of long chromatin fibers as a result of direct interaction of the cholesterol to six nucleosomal binding sites.

Regulation of Nucleosome Stacking and Chromatin Compaction by the Histone H4 N-Terminal Tail-H2A Acidic Patch Interaction

Chromatin folding and dynamics are critically dependent on nucleosome-nucleosome interactions with important contributions from internucleosome binding of the histone H4 N-terminal tail K16-R23 domain to the surface of the H2A/H2B dimer. The H4 Lys16 plays a pivotal role in this regard. Using in vitro reconstituted 12-mer nucleosome arrays, we have investigated the mechanism of the H4 N-terminal tail in maintaining nucleosome-nucleosome stacking and mediating intra- and inter-array chromatin compaction, with emphasis on the role of K16 and the positive charge region, R17-R23. Analytical ultracentrifugation sedimentation velocity experiments and precipitation assays were employed to analyze effects on chromatin folding and self-association, respectively. Effects on chromatin folding caused by various mutations and modifications at position K16 in the H4 histone were studied. Additionally, using charge-quenching mutations, we characterized the importance of the interaction of the residues within the H4 positive charge region R17-R23 with the H2A acidic patch of the adjacent nucleosome. Furthermore, crosslinking experiments were conducted to establish the proximity of the basic tail region to the acidic patch. Our data indicate that the positive charge and length of the side chain of H4 K16 are important for its access to the adjacent nucleosome in the process of nucleosome-nucleosome stacking and array folding. The location and orientation of the H4 R17-R23 domain on the H2A/H2B dimer surface of the neighboring nucleosome core particle (NCP) in the compacted chromatin fiber were established. The dominance of electrostatic interactions in maintaining intra-array interaction was demonstrated.

Atomic force microscopy of chromatin arrays reveal non-monotonic salt dependence of array compaction in solution

Compaction of DNA in chromatin is a hallmark of the eukaryotic cell and unravelling its structure is required for an understanding of DNA involving processes. Despite strong experimental efforts, many questions concerning the DNA packing are open. In particular, it is heavily debated whether an ordered structure referred to as the “30 nm fibre” exist in vivo. Scanning probe microscopy has become a cutting edge technology for the high-resolution imaging of DNA- protein complexes. Here, we perform high-resolution atomic force microscopy of non-cross-linked chromatin arrays in liquid, under different salt conditions. A statistical analysis of the data reveals that array compaction is salt dependent in a non-monotonic fashion. A simple physical model can qualitatively explain the observed findings due to the opposing effects of salt dependent stiffening of DNA, nucleosome stability and histone-histone interactions. While for different salt concentrations different compaction states are observed, our data do not provide support for the existence of regular chromatin fibres. Our studies add new insights into chromatin structure, and with that contribute to a further understanding of the DNA condensation.

Calcium ions function as a booster of chromosome condensation

Chromosome condensation is essential for the faithful transmission of genetic information to daughter cells during cell division. The depletion of chromosome scaffold proteins does not prevent chromosome condensation despite structural defects. This suggests that other factors contribute to condensation. Here we investigated the contribution of divalent cations, particularly Ca2+, to chromosome condensation in vitro and in vivo. Ca2+ depletion caused defects in proper mitotic progression, particularly in chromosome condensation after the breakdown of the nuclear envelope. Fluorescence lifetime imaging microscopy-Förster resonance energy transfer and electron microscopy demonstrated that chromosome condensation is influenced by Ca2+. Chromosomes had compact globular structures when exposed to Ca2+ and expanded fibrous structures without Ca2+. Therefore, we have clearly demonstrated a role for Ca2+ in the compaction of chromatin fibres.

Protein/DNA interactions in complex DNA topologies: expect the unexpected

DNA supercoiling results in compacted DNA structures that can bring distal sites into close proximity. It also changes the local structure of the DNA, which can in turn influence the way it is recognised by drugs, other nucleic acids and proteins. Here, we discuss how DNA supercoiling and the formation of complex DNA topologies can affect the thermodynamics of DNA recognition. We then speculate on the implications for transcriptional control and the three-dimensional organisation of the genetic material, using examples from our own simulations and from the literature. We introduce and discuss the concept of coupling between the multiple length-scales associated with hierarchical nuclear structural organisation through DNA supercoiling and topology.

Protein/DNA interactions in complex DNA topologies: expect the unexpected

DNA supercoiling results in compacted DNA structures that can bring distal sites into close proximity. It also changes the local structure of the DNA, which can in turn influence the way it is recognised by drugs, other nucleic acids and proteins. Here, we discuss how DNA supercoiling and the formation of complex DNA topologies can affect the thermodynamics of DNA recognition. We then speculate on the implications for transcriptional control and the three-dimensional organisation of the genetic material, using examples from our own simulations and from the literature. We introduce and discuss the concept of coupling between the multiple length-scales associated with hierarchical nuclear structural organisation through DNA supercoiling and topology.

Protein-DNA and ion-DNA interactions revealed through contrast variation SAXS

Understanding how DNA carries out its biological roles requires knowledge of its interactions with biological partners. Since DNA is a polyanionic polymer, electrostatic interactions contribute significantly. These interactions are mediated by positively charged protein residues or charge compensating cations. Direct detection of these partners and/or their effect on DNA conformation poses challenges, especially for monitoring conformational dynamics in real time. Small-angle x-ray scattering (SAXS) is uniquely sensitive to both the conformation and local environment (i.e. protein partner and associated ions) of the DNA. The primary challenge of studying multi-component systems with SAXS lies in resolving how each component contributes to the measured scattering. Here, we review two contrast variation (CV) strategies that enable targeted studies of the structures of DNA or its associated partners. First, solution contrast variation enables measurement of DNA conformation within a protein-DNA complex by masking out the protein contribution to the scattering profile. We review a specific example, in which the real-time unwrapping of DNA from a nucleosome core particle is measured during salt-induced disassembly. The second method, heavy atom isomorphous replacement, reports the spatial distribution of the cation cloud around duplex DNA by exploiting changes in the scattering strength of cations with varying atomic numbers. We demonstrate the application of this approach to provide the spatial distribution of monovalent cations (Na+, K+, Rb+, Cs+) around a standard 25-base pair DNA. The CV strategies presented here are valuable tools for understanding DNA interactions with its biological partners.

Molecular dynamics simulations demonstrate the regulation of DNA-DNA attraction by H4 histone tail acetylations and mutations

The positively charged N-terminal histone tails play a crucial role in chromatin compaction and are important modulators of DNA transcription, recombination, and repair. The detailed mechanism of the interaction of histone tails with DNA remains elusive. To model the unspecific interaction of histone tails with DNA, all-atom molecular dynamics (MD) simulations were carried out for systems of four DNA 22-mers in the presence of 20 or 16 short fragments of the H4 histone tail (variations of the 16-23 a. a. KRHRKVLR sequence, as well as the unmodified fragment a. a.13-20, GGAKRHRK). This setup with high DNA concentration, explicit presence of DNA-DNA contacts, presence of unstructured cationic peptides (histone tails) and K(+) mimics the conditions of eukaryotic chromatin. A detailed account of the DNA interactions with the histone tail fragments, K(+) and water is presented. Furthermore, DNA structure and dynamics and its interplay with the histone tail fragments binding are analysed. The charged side chains of the lysines and arginines play major roles in the tail-mediated DNA-DNA attraction by forming bridges and by coordinating to the phosphate groups and to the electronegative sites in the minor groove. Binding of all species to DNA is dynamic. The structure of the unmodified fully-charged H4 16-23 a.a. fragment KRHRKVLR is dominated by a stretched conformation. The H4 tail a. a. fragment GGAKRHRK as well as the H4 Lys16 acetylated fragment are highly flexible. The present work allows capturing typical features of the histone tail-counterion-DNA structure, interaction and dynamics.

Sedimentation Velocity Analysis of Large Oligomeric Chromatin Complexes Using Interference Detection

Sedimentation velocity experiments measure the transport of molecules in solution under centrifugal force. Here, we describe a method for monitoring the sedimentation of very large biological molecular assemblies using the interference optical systems of the analytical ultracentrifuge. The mass, partial-specific volume, and shape of macromolecules in solution affect their sedimentation rates as reflected in the sedimentation coefficient. The sedimentation coefficient is obtained by measuring the solute concentration as a function of radial distance during centrifugation. Monitoring the concentration can be accomplished using interference optics, absorbance optics, or the fluorescence detection system, each with inherent advantages. The interference optical system captures data much faster than these other optical systems, allowing for sedimentation velocity analysis of extremely large macromolecular complexes that sediment rapidly at very low rotor speeds. Supramolecular oligomeric complexes produced by self-association of 12-mer chromatin fibers are used to illustrate the advantages of the interference optics. Using interference optics, we show that chromatin fibers self-associate at physiological divalent salt concentrations to form structures that sediment between 10,000 and 350,000S. The method for characterizing chromatin oligomers described in this chapter will be generally useful for characterization of any biological structures that are too large to be studied by the absorbance optical system.

The chromatin fiber: multiscale problems and approaches

The structure of chromatin, affected by many factors from DNA linker lengths to posttranslational modifications, is crucial to the regulation of eukaryotic cells. Combined experimental and computational methods have led to new insights into its structural and dynamical features, from interactions due to the flexible core histone tails or linker histones to the physical mechanism driving the formation of chromosomal domains. Here we present a perspective of recent advances in chromatin modeling techniques at the atomic, mesoscopic, and chromosomal scales with a view toward developing multiscale computational strategies to integrate such findings. Innovative modeling methods that connect molecular to chromosomal scales are crucial for interpreting experiments and eventually deciphering the complex dynamic organization and function of chromatin in 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 kBT. 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.

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+).

Chromatin compaction under mixed salt conditions: opposite effects of sodium and potassium ions on nucleosome array folding

It is well known that chromatin structure is highly sensitive to the ionic environment. However, the combined effects of a physiologically relevant mixed ionic environment of K⁺, Mg²⁺ and Na⁺, which are the main cations of the cell cytoplasm, has not been systematically investigated. We studied folding and self-association (aggregation) of recombinant 12-mer nucleosome arrays with 177 bp DNA repeat length in solutions of mixtures of K⁺ and Mg²⁺ or Na⁺ and Mg²⁺. In the presence of Mg²⁺, the addition of sodium ions promotes folding of array into 30-nm fibres, whereas in mixtures of K⁺ and Mg²⁺, potassium ions abrogate folding. We found that self-association of nucleosome arrays in mixed salt solutions is synergistically promoted by Mg²⁺ and monovalent ions, with sodium being slightly more efficient than potassium in amplifying the self-association. The results highlight the importance of a mixed ionic environment for the compaction of chromatin under physiological conditions and demonstrate the complicated nature of the various factors that determine and regulate chromatin compaction in vivo.

Uncovering protein polyamination by the spermine-specific antiserum and mass spectrometric analysis

The polyamines spermidine and spermine, and their precursor putrescine, have been shown to play an important role in cell migration, proliferation, and differentiation. Because of their polycationic property, polyamines are traditionally thought to be involved in DNA replication, gene expression, and protein translation. However, polyamines can also be covalently conjugated to proteins by transglutaminase 2 (TG2). This modification leads to an increase in positive charge in the polyamine-incorporated region which significantly alters the structure of proteins. It is anticipated that protein polyamine conjugation may affect the protein-protein interaction, protein localization, and protein function of the TG2 substrates. In order to investigate the roles of polyamine modification, we synthesized a spermine-conjugated antigen and generated an antiserum against spermine. In vitro TG2-catalyzed spermine incorporation assays were carried out to show that actin, tubulins, heat shock protein 70 and five types of histone proteins were modified with spermine, and modification sites were also identified by liquid chromatography and linear ion trap-orbitrap hybrid mass spectrometry. Subsequent mass spectrometry-based shotgun proteomic analysis also identified 254 polyaminated sites in 233 proteins from the HeLa cell lysate catalyzed by human TG2 with spermine, thus allowing, for the first time, a global appraisal of site-specific protein polyamination. Global analysis of mouse tissues showed that this modification really exists in vivo. Importantly, we have demonstrated that there is a new histone modification, polyamination, in cells. However, the functional significance of histone polyamination demands further investigations.

Elucidating internucleosome interactions and the roles of histone tails

The nucleosome is the first level of genome organization and regulation in eukaryotes where negatively charged DNA is wrapped around largely positively charged histone proteins. Interaction between nucleosomes is dominated by electrostatics at long range and guided by specific contacts at short range, particularly involving their flexible histone tails. We have thus quantified how internucleosome interactions are modulated by salts (KCl, MgCl2) and histone tail deletions (H3, H4 N-terminal), using small-angle x-ray scattering and theoretical modeling. We found that measured effective charges at low salts are ∼1/5th of the theoretically predicted renormalized charges and that H4 tail deletion suppresses the attraction at high salts to a larger extent than H3 tail deletion.

Counting the number of magnesium ions bound to the surface-immobilized thymine oligonucleotides that comprise spherical nucleic acids

Label-free studies carried out under aqueous phase conditions quantify the number of Mg(2+) ions binding to surface-immobilized T40 sequences, the subsequent reordering of DNA on the surface, and the consequences of Mg(2+) binding for DNA-DNA interactions. Second harmonic generation measurements indicate that, within error, 18-20 Mg(2+) ions are bound to the T40 strand at saturation and that the metal-DNA interaction is associated with a near 30% length contraction of the strand. Structural reordering, evaluated using vibrational sum frequency generation, atomic force microscopy, and dynamic light scattering, is attributed to increased charge screening as the Mg(2+) ions bind to the negatively charged DNA, reducing repulsive Coulomb forces between nucleotides and allowing the DNA single strands to collapse or coil upon themselves. The impact of Mg(2+) binding on DNA hybridization and duplex stability is assessed with spherical nucleic acid (SNA) gold nanoparticle conjugates in order to determine an optimal working range of Mg(2+) concentrations for DNA-DNA interactions in the absence of NaCl. The findings are consistent with a charge titration effect in which, in the absence of NaCl, (1) hybridization does not occur at room temperature if an average of 17.5 or less Mg(2+) ions are bound per T40 strand, which is not reached until the bulk Mg(2+) concentration approaches 0.5 mM; (2) hybridization proceeds, albeit with low duplex stability having an average Tm of 31(3)°C, if an average of 17.5-18.0 Mg(2+) ions are bound; and (3) highly stable duplexes having a Tm of 64(2)°C form if 18.5-19.0 Mg(2+) ions are bound, corresponding to saturation of the T40 strand.

Modulation of Higher Order Chromatin Conformation in Mammalian Cell Nuclei Can Be Mediated by Polyamines and Divalent Cations

The organisation of the large volume of mammalian genomic DNA within cell nuclei requires mechanisms to regulate chromatin compaction involving the reversible formation of higher order structures. The compaction state of chromatin varies between interphase and mitosis and is also subject to rapid and reversible change upon ATP depletion/repletion. In this study we have investigated mechanisms that may be involved in promoting the hyper-condensation of chromatin when ATP levels are depleted by treating cells with sodium azide and 2-deoxyglucose. Chromatin conformation was analysed in both live and permeabilised HeLa cells using FLIM-FRET, high resolution fluorescence microscopy and by electron spectroscopic imaging microscopy. We show that chromatin compaction following ATP depletion is not caused by loss of transcription activity and that it can occur at a similar level in both interphase and mitotic cells. Analysis of both live and permeabilised HeLa cells shows that chromatin conformation within nuclei is strongly influenced by the levels of divalent cations, including calcium and magnesium. While ATP depletion results in an increase in the level of unbound calcium, chromatin condensation still occurs even in the presence of a calcium chelator. Chromatin compaction is shown to be strongly affected by small changes in the levels of polyamines, including spermine and spermidine. The data are consistent with a model in which the increased intracellular pool of polyamines and divalent cations, resulting from depletion of ATP, bind to DNA and contribute to the large scale hyper-compaction of chromatin by a charge neutralisation mechanism.

Structure-driven homology pairing of chromatin fibers: the role of electrostatics and protein-induced bridging

Chromatin domains formed in vivo are characterized by different types of 3D organization of interconnected nucleosomes and architectural proteins. Here, we quantitatively test a hypothesis that the similarities in the structure of chromatin fibers (which we call "structural homology") can affect their mutual electrostatic and protein-mediated bridging interactions. For example, highly repetitive DNA sequences in heterochromatic regions can position nucleosomes so that preferred inter-nucleosomal distances are preserved on the surfaces of neighboring fibers. On the contrary, the segments of chromatin fiber formed on unrelated DNA sequences have different geometrical parameters and lack structural complementarity pivotal for stable association and cohesion. Furthermore, specific functional elements such as insulator regions, transcription start and termination sites, and replication origins are characterized by strong nucleosome ordering that might induce structure-driven iterations of chromatin fibers. We propose that shape-specific protein-bridging interactions facilitate long-range pairing of chromatin fragments, while for closely-juxtaposed fibers electrostatic forces can in addition yield fine-tuned structure-specific recognition and pairing. These pairing effects can account for some features observed for mitotic and inter-phase chromatins.

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.

Modelling chromatin structure and dynamics: status and prospects

The packaging of genomic DNA into chromatin in the eukaryotic cell nucleus demands extensive compaction. This requires attractive nucleosome-nucleosome interactions to overcome repulsion between the negatively charged DNA segments as well as other constraints. At the same time, DNA must be dynamically accessible to the cellular machinery that operates on it. Recent progress in the experimental characterisation of the higher order structure and dynamics of well-defined chromatin fibres has stimulated the attempts at theoretical description of chromatin and the nucleosome. Here we review the present status of chromatin modelling, with particular emphasis on coarse-grained computer simulation models, the role of electrostatic interactions, and discuss future perspectives in the field.

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.