Explicit ion modeling predicts physicochemical interactions for chromatin organization

eLife assessment

The authors have developed a compelling coarse-grained simulation approach for nucleosome-nucleosome interactions within a chromatin array. The data presented are solid and provide new insights that allow for predictions of how chromatin interactions might occur in vivo. The tools presented herein will be valuable for the chromosome biology field.

https://doi.org/10.7554/eLife.90073.3.sa0

Significance of the findings:

Valuable: Findings that have theoretical or practical implications for a subfield

Landmark
Fundamental
Important
Valuable
Useful

Strength of evidence:

Compelling: Evidence that features methods, data and analyses more rigorous than the current state-of-the-art

Solid: Methods, data and analyses broadly support the claims with only minor weaknesses

Exceptional
Compelling
Convincing
Solid
Incomplete
Inadequate

During the peer-review process the editor and reviewers write an eLife Assessment that summarises the significance of the findings reported in the article (on a scale ranging from landmark to useful) and the strength of the evidence (on a scale ranging from exceptional to inadequate). Learn more about eLife Assessments

Abstract
Introduction
Results
Discussion
Materials and methods
Appendix 1
Data availability
References
Article and author information
Metrics

Abstract

Molecular mechanisms that dictate chromatin organization in vivo are under active investigation, and the extent to which intrinsic interactions contribute to this process remains debatable. A central quantity for evaluating their contribution is the strength of nucleosome-nucleosome binding, which previous experiments have estimated to range from 2 to 14 k_BT. We introduce an explicit ion model to dramatically enhance the accuracy of residue-level coarse-grained modeling approaches across a wide range of ionic concentrations. This model allows for de novo predictions of chromatin organization and remains computationally efficient, enabling large-scale conformational sampling for free energy calculations. It reproduces the energetics of protein-DNA binding and unwinding of single nucleosomal DNA, and resolves the differential impact of mono- and divalent ions on chromatin conformations. Moreover, we showed that the model can reconcile various experiments on quantifying nucleosomal interactions, providing an explanation for the large discrepancy between existing estimations. We predict the interaction strength at physiological conditions to be 9 k_BT, a value that is nonetheless sensitive to DNA linker length and the presence of linker histones. Our study strongly supports the contribution of physicochemical interactions to the phase behavior of chromatin aggregates and chromatin organization inside the nucleus.

Introduction

Three-dimensional genome organization plays essential roles in numerous DNA-templated processes (Dekker et al., 2013; Bonev and Cavalli, 2016; Finn and Misteli, 2019; Misteli, 2020; Lin et al., 2021b). Understanding the molecular mechanisms for its establishment could improve our understanding of these processes and facilitate genome engineering. Advancements in high-throughput sequencing and microscopic imaging have enabled genome-wide structural characterization, revealing a striking compartmentalization of chromatin at large scales (Lieberman-Aiden et al., 2009; Quinodoz et al., 2018; Su et al., 2020; Takei et al., 2021). For example, A compartments are enriched with euchromatin and activating post-translational modifications to histone proteins. They are often spatially segregated from B compartments that enclose heterochromatin with silencing histone marks (Gibcus and Dekker, 2013; Finn and Misteli, 2019; Misteli, 2020; Mirny and Dekker, 2022; Xie and Zhang, 2019).

Compartmentalization has been proposed to arise from the microphase separation of different chromatin types as in block copolymer systems (Fujishiro and Sasai, 2022; Jost et al., 2014; Falk et al., 2019; Bajpai et al., 2021; Laghmach et al., 2020; Hu et al., 2013; Lesne et al., 2014; Di Pierro et al., 2016; Xie et al., 2017; Yildirim and Feig, 2018; MacPherson et al., 2018; Shi and Thirumalai, 2021; Brahmachari et al., 2022). However, the molecular mechanisms that drive the microphase separation are not yet fully understood. Protein molecules that recognize specific histone modifications have frequently been found to undergo liquid-liquid phase separation (Larson et al., 2017; Kent et al., 2020; Xie et al., 2022; Leicher et al., 2022; Latham and Zhang, 2021; Lin et al., 2021a; MacPherson et al., 2018), potentially contributing to chromatin demixing. Demixing can also arise from interactions between chromatin and various nuclear landmarks such as nuclear lamina and speckles (Brahmachari et al., 2022; Falk et al., 2019; Mirny and Dekker, 2022; Kamat et al., 2023), as well as active transcriptional processes (Hilbert et al., 2021; Jiang et al., 2022; Brahmachari et al., 2023; Goychuk et al., 2023). Furthermore, recent studies have revealed that nucleosome arrays alone can undergo spontaneous phase separation (Gibson et al., 2019; Strickfaden et al., 2020; Zhang et al., 2022), indicating that compartmentalization may be an intrinsic property of chromatin driven by nucleosome-nucleosome interactions.

The relevance of physicochemical interactions between nucleosomes to chromatin organization in vivo has been constantly debated, partly due to the uncertainty in their strength (Kruithof et al., 2009; Cui and Bustamante, 2000; Kaczmarczyk et al., 2020; Funke et al., 2016). Examining the interactions between native nucleosomes poses challenges due to the intricate chemical modifications that histone proteins undergo within the nucleus and the variations in their underlying DNA sequences (Fenley et al., 2010; Fenley et al., 2018). Many in vitro experiments have opted for reconstituted nucleosomes that lack histone modifications and feature well-positioned 601-sequence DNA (Lowary and Widom, 1998) to simplify the chemical complexity. These experiments aim to establish a fundamental reference point, a baseline for understanding the strength of interactions within native nucleosomes. Nevertheless, even with reconstituted nucleosomes, a consensus regarding the significance of their interactions remains elusive. For example, using force-measuring magnetic tweezers, Kruithof et al. estimated the inter-nucleosome binding energy to be ∼14 k_BT (Kruithof et al., 2009). On the other hand, Funke et al. introduced a DNA-origami-based force spectrometer to directly probe the interaction between a pair of nucleosomes (Funke et al., 2016), circumventing any potential complications from interpretations of single-molecule traces of nucleosome arrays. Their measurement reported a much weaker binding free energy of approximately 2 k_BT. This large discrepancy in the reported reference values complicates a further assessment of the interactions between native nucleosomes and their contribution to chromatin organization in vivo.

Computational modeling is well suited for reconciling the discrepancy across experiments and determining the strength of inter-nucleosome interactions. The high computational cost of atomistic simulations (Winogradoff et al., 2015; Woods et al., 2021; Li et al., 2023) has inspired several groups to calculate the nucleosome binding free energy with coarse-grained models (Moller et al., 2019; Farr et al., 2021). However, the complex distribution of charged amino acids and nucleotides at nucleosome interfaces places a high demand on force field accuracy. In particular, most existing models adopt a mean-field approximation with the Debye-Hückel theory (Phillips, 2012) to describe electrostatic interactions in an implicit-solvent environment (Izadi et al., 2016; Bascom and Schlick, 2018; Moller et al., 2019; Farr et al., 2021), preventing an accurate treatment of the complex salt conditions explored in experiments. Further force field development is needed to improve the accuracy of coarse-grained modeling across different experimental settings (Freeman et al., 2011; Hinckley and de Pablo, 2015; Sun et al., 2022; Hayes et al., 2015).

We introduce a residue-level coarse-grained explicit ion model for simulating chromatin conformations and quantifying inter-nucleosome interactions. We validate our model’s accuracy through extensive simulations, demonstrating that it reproduces the binding affinities of protein-DNA complexes (Privalov et al., 2011) and energetic cost of nucleosomal DNA unwinding (Hall et al., 2009). Further simulations of chromatin at various salt concentrations reproduce experimentally measured sedimentation coefficients (Correll et al., 2012). We also reveal extensive close contacts between histone proteins and DNA across nucleosomes, the perturbation of which explains the discrepancy among various experimental studies. Finally, we determined the binding free energy between a pair of nucleosomes under physiological salt concentrations as ∼9 k_BT. While longer linker DNA would reduce this binding energy, linker histones can more than compensate this reduction to mediate inter-nucleosome interactions with disordered, charged terminal tails. Our study supports the importance of intrinsic physicochemical interactions in chromatin organization in vivo.

Results

Counterion condensation accommodates nucleosomal DNA unwrapping

Various single-molecule studies have been carried out to probe the stability of nucleosomes and the interactions between histone proteins and DNA (Bennink et al., 2001; Cui and Bustamante, 2000; Pope et al., 2005; Bancaud et al., 2007; Hall et al., 2009). The DNA-unzipping experiment performed by Hall et al., 2009, is particularly relevant since the measured forces can be converted into a free energy profile of DNA unwinding at a base-pair resolution, as shown by Forties et al. with a continuous-time Markov model (Forties et al., 2011). The high-resolution quantification of nucleosome energetics is valuable for benchmarking the accuracy of computational models.

We introduce a coarse-grained explicit ion model for chromatin simulations (Figure 1). The model represents each amino acid with one coarse-grained bead and three beads per nucleotide. It resolves the differences among various chemical groups to accurately describe biomolecular interactions with physical chemistry potentials. Our explicit representation of monovalent and divalent ions enables a faithful description of counterion condensation and its impact on electrostatic interactions between protein and DNA molecules. Additional model details are provided in the Materials and methods and Appendix.

Figure 1

Download asset Open asset

Illustration of the residue-level coarse-grained explicit ion model for chromatin simulations.

The left panel presents a snapshot for the simulation box of a 147 bp nucleosome in a solution of 100 mM NaCl and 0.5 mM MgCl₂. The nucleosomal DNA and histone proteins are colored in red and white, respectively. The zoom-in on the right highlights the condensation of ions around the nucleosome, with Na⁺ in cyan and Mg²⁺ in yellow. Negative residues of the histone proteins are colored in pink.

We performed umbrella simulations (Torrie and Valleau, 1977) to determine the free energy profile of nucleosomal DNA unwinding. The experimental buffer condition of 0.10 M NaCl and 0.5 mM MgCl₂ (Hall et al., 2009) was adopted in simulations for direct comparison. As shown in Figure 2B, the simulated values match well with experimental results over a wide range. Furthermore, we computed the binding free energy for a diverse set of protein-DNA complexes and the simulated values again match well with experimental data (Figure 2—figure supplement 1), supporting the model’s accuracy.

Figure 2 with 1 supplement see all

Download asset Open asset

Explicit ion modeling reproduces the energetics of nucleosomal DNA unwrapping.

(A) Illustration of the umbrella simulation setup using the end-to-end distance between two DNA termini as the collective variable. The same color scheme as in Figure 1 is adopted. Only ions close to the nucleosomes are shown for clarity. (B) Comparison between simulated (black) and experimental (red) free energy profile as a function of the unwrapped DNA base pairs. Error bars were computed as the standard deviation of three independent estimates. (C) The average number of Na⁺ ions within 10 Å of the nucleosomal DNA (top) and Cl⁻ions within 10 Å of histone proteins (bottom) are shown as a function of the unwrapped DNA base pairs. Error bars were computed as the standard deviation of three independent estimates.

Counterions are often released upon protein-DNA binding to make room for close contacts at the interface, contributing favorably to the binding free energy in the form of entropic gains (Schiessel, 2003). However, previous studies have shown that the histone-DNA interface in a fully wrapped nucleosome configuration is not tightly sealed but instead permeated with water molecules and mobile ions (Davey et al., 2002; Materese et al., 2009). Given their presence in the bound form, how these counterions contribute to nucleosomal DNA unwrapping remains to be shown. We calculated the number of DNA-bound cations and protein-bound anions as DNA unwraps. Our results, shown in Figure 2C, indicate that only a modest amount of extra Na⁺ and Cl⁻ ions becomes associated with the nucleosome as the outer DNA layer unwraps. However, significantly more ions become bound when the inner layer starts to unwrap (after 73 bp). These findings suggest that counterion release may contribute more significantly to the inner layer wrapping, potentially caused by a tighter protein-DNA interface.

Charge neutralization with Mg²⁺ compacts chromatin

In addition to contributing to the stability of individual nucleosomes, counterions can also impact higher-order chromatin organization. Numerous groups have characterized the structures of nucleosome arrays (Widom, 1986; Schwarz et al., 1996; Engelhardt, 2004; Correll et al., 2012; Grigoryev et al., 2009; Allahverdi et al., 2015), revealing a strong dependence of chromatin folding on the concentration and valence of cations.

To further understand the role of counterions in chromatin organization, we studied a 12-mer with 20-bp-long linker DNA under different salt conditions. We followed the experiment setup by Correll et al., 2012, that immerses chromatin in solutions with 5 mM NaCl, 150 mM NaCl, 0.6 mM MgCl₂, or 1 mM MgCl₂. To facilitate conformational sampling, we carried out umbrella simulations with a collective variable that quantifies the similarity between a given configuration and a reference two-start helical structure. Simulation details and the precise definition of the collective variable are provided in the Materials and methods and Appendix. Data from different umbrella windows were combined together with proper reweighting (Kumar et al., 1992) for analysis.

As shown in Figure 3A, the average sedimentation coefficients determined from our simulations match well with experimental values. Specifically, the simulations reproduce the strong contrast in chromatin size between the two systems with different NaCl concentrations. Chromatin under 5 mM NaCl features an extended configuration with minimal stacking between one and three nucleosomes (Figure 3B). On the other hand, the compaction is evident at 150 mM NaCl. Notably, in agreement with previous studies (Ding et al., 2021; Liu et al., 2022; Cai et al., 2018; Dombrowski et al., 2022), we observe tri-nucleosome configurations as chromatin extends. Finally, the simulations also support that divalent ions are more effective in packaging chromatin than NaCl. Even in the presence of 0.6 mM MgCl₂, the chromatin sedimentation coefficient is comparable to that obtained at 150 mM of NaCl.

Figure 3 with 1 supplement see all

Download asset Open asset

Explicit ion modeling predicts salt-dependent conformations of a 12-mer nucleosome array.

(A) Top: Comparison of simulated and experimental (Correll et al., 2012) sedimentation coefficients of chromatin at different salt concentrations. Bottom: Number of DNA charges neutralized by bound cations (yellow, left y-axis label) and the fraction of ions bound to DNA (red, right y-axis label) at different salt concentrations. The error bars were estimated from the standard deviation of simulated probability distributions (Figure 3—figure supplement 1). (B) Representative chromatin structures with sedimentation coefficients around the mean values at different salt concentrations.

We further characterized ions that are in close contact with DNA to understand their impact on chromatin organization. Our simulations support the condensation of cations, especially for divalent ions (Figure 3A, bottom) as predicted by the Manning theory (Manning, 1978; Clark and Kimura, 1990). Ion condensation weakens the repulsion among DNA segments that prevents chromatin from collapsing. Notably, the fraction of bound Mg²⁺ is much higher than Na⁺. Correspondingly, the amount of neutralized negative charges is always greater in systems with divalent ions, despite the significantly lower salt concentrations. The difference between the two types of ions arises from the more favorable interactions between Mg²⁺ and phosphate groups that more effectively offset the entropy loss due to ion condensation (Clark and Kimura, 1990). While higher concentrations of NaCl do not dramatically neutralize more charges, the excess ions provide additional screening to weaken the repulsion among DNA segments, stabilizing chromatin compaction.

Close contacts drive nucleosome binding free energy

Encouraged by the explicit ion model’s accuracy in reproducing experimental measurements of single-nucleosomes and nucleosome arrays, we moved to directly quantify the strength of inter-nucleosomes interactions. We once again focus on reconstituted nucleosomes for a direct comparison with in vitro experiments. These experiments have yielded a wide range of values, ranging from 2 to 14 k_BT (Funke et al., 2016; Cui and Bustamante, 2000; Kruithof et al., 2009). Accurate quantification will offer a reference value for conceptualizing the significance of physicochemical interactions among native nucleosomes in chromatin organization in vivo.

To reconcile the discrepancy among various experimental estimations, we directly calculated the binding free energy between a pair of nucleosomes with umbrella simulations. We adopted the same ionic concentrations as in the experiment performed by Funke et al., 2016, with 35 mM NaCl and 11 mM MgCl₂. We focus on this study since the experiment directly measured the inter-nucleosomal interactions, allowing straightforward comparison with simulation results. Furthermore, the reported value for nucleosome binding free energy deviates the most from other studies. In one set of umbrella simulations, we closely mimicked the DNA-origami device employed by Funke et al. to move nucleosomes along a predefined path for disassociation (Figure 4A, A1 to A3). For example, neither nucleosome can freely rotate (Figure 4—figure supplement 1); the first nucleosome is restricted to the initial position, and the second nucleosome can only move within the Y-Z plane along the arc 15 nm away from the origin. For comparison, we performed a second set of independent simulations without imposing any restrictions on nucleosome orientations. Additional simulation details can be found in Materials and methods and Appendix.

Figure 4 with 4 supplements see all

Download asset Open asset

Close contacts give rise to strong inter-nucleosomal interactions.

(A) Illustration of the simulation protocol employed to mimic the nucleosome unbinding pathway dictated by the DNA-origami device (Funke et al., 2016). The three configurations, A1, A2, and A3, corresponding to the three cyan dots in part B at distances 62.7, 80.2, and 96.3 Å. For comparison, a tightly bound configuration uncovered in simulations without any restraints of nucleosome movement is shown as A1’. The number of contacts formed by histone tails and DNA (Htail-DNA) and by histone core and DNA (Hcore-DNA) from different nucleosomes is shown for A1 and A1’. (B) Free energy profile as a function of the distance between the geometric centers of the two nucleosomes, computed from unrestrained (black) and DNA-origami-restrained simulations (red). Error bars were computed as the standard deviation of three independent estimates. (C) Average inter-nucleosomal contacts between DNA and histone tail (orange) and core (blue) residues, computed from unrestrained and DNA-origami-restrained simulations. Error bars were computed as the standard deviation of three independent estimates.

Strikingly, the two sets of simulations produced dramatically different binding free energies. Restricting nucleosome orientations produced a binding free energy of ∼2 k_BT, reproducing the experimental value (Figure 4B, Figure 4—figure supplement 2). On the other hand, the binding free energy increased to 15 k_BT upon removing the constraints.

Further examination of inter-nucleosomal contacts revealed the origin of the dramatic difference in nucleosome binding free energies. As shown in Figure 4C, the average number of contacts formed between histone tails and DNA from different nucleosomes is around 150 and 10 in the two sets of simulations. A similar trend is observed for histone core-DNA contacts across nucleosomes. The differences are most dramatic at small distances (Figure 4B, Figure 4—figure supplement 3) and are clearly visible in the most stable configurations. For example, from the unrestricted simulations, the most stable binding mode corresponds to a configuration in which the two nucleosomes are almost parallel to each other (see configuration A1’ in Figure 4A), with the angle between the two nucleosome planes close to zero (Figure 4B, Figure 4—figure supplement 4). However, the inherent design of the DNA-origami device renders this binding mode inaccessible, and the smallest angle between the two nucleosome planes is around 23° (see configuration A1 in Figure 4A). Therefore, a significant loss of inter-nucleosomal contacts caused the small binding free energy seen experimentally.

Modulation of nucleosome binding free energy by in vivo factors

The predicted strength for unrestricted inter-nucleosome interactions supports their significant contribution to chromatin organization in vivo. However, the salt concentration studied above and in the DNA-origami experiment is much higher than the physiological value (Kaczmarczyk et al., 2020). To further evaluate the in vivo significance of inter-nucleosome interactions, we computed the binding free energy at the physiological salt concentration of 150 mM NaCl and 2 mM of MgCl₂.

We observe a strong dependence of nucleosome orientations on the inter-nucleosome distance. A collective variable, $θ$ , was introduced to quantify the angle between the two nucleosomal planes (Figure 5A). As shown in two-dimensional binding free energy landscape of inter-nucleosome distance, $r$ , and $θ$ (Figure 5B), at small distances (∼60 Å), the two nucleosomes prefer a face-to-face binding mode with small $θ$ values. As the distance increases, the nucleosomes will almost undergo a 90° rotation to adopt perpendicular positions. Such orientations allow the nucleosomes to remain in contact and is more energetically favorable. The orientation preference gradually diminishes at large distances once the two nucleosomes are completely detached. Importantly, we observed a strong inter-nucleosomal interaction with two nucleosomes wrapped by 147 bp 601-sequence DNA (∼9 k_BT).

Figure 5 with 3 supplements see all

Download asset Open asset

Simulations predict significant inter-nucleosome interactions at physiological conditions.

(A) Illustration of the collective variable, $θ$ , defined as the angle between two nucleosomal planes, and $r$ defined as the distance between the nucleosome geometric centers. $\vec{w_{1}}$ and $\vec{w_{2}}$ represent the axes perpendicular to the nucleosomal planes. (B) The 2D binding free energy profile as a function of $θ$ and $r$ at the physiological salt condition (150 mM NaCl and 2 mM MgCl₂) for nucleosomes with the 601 sequence. (C) Dependence of nucleosome binding free energy on nucleosome repeat length (NRL) and linker histone H1.0. Error bars were computed as the standard deviation of three independent estimates. (D) Representative structure showing linker histones (red and blue) mediating inter-nucleosomal contacts.

Furthermore, we found that the nucleosome binding free energy is minimally impacted by the precise DNA sequence. For example, when the 601 sequence is replaced with poly-dA:dT or poly-dG:dC, the free energy only varied by ∼2 k_BT (Figure 5—figure supplement 1). However, the poly-dA:dT sequence produced stronger binding while poly-dG:dC weakened the interactions. The sequence specific effects are potentially due to the increased stiffness of poly-dA:dT DNA (Ortiz and de Pablo, 2011), which causes the DNA to unwrap more frequently, increasing cross-nucleosome contacts at larger distances (Figure 5—figure supplement 2).

In addition to variations in DNA sequences, in vivo nucleosomes also feature different linker lengths. We performed simulations that extend the 601 sequence with 10 extra base pairs of poly-dA:dT sequence at each end, reaching a nucleosome repeat length (NRL) of 167 bp. Consistent with previous studies (Mangenot et al., 2002; Correll et al., 2012; Huang et al., 2018), increasing the NRL weakened inter-nucleosomal interactions (Figure 5C and Figure 5—figure supplement 3), reducing the binding free energy to ∼6 k_BT.

Importantly, we found that the weakened interactions upon extending linker DNA can be more than compensated for by the presence of histone H1 proteins. This is demonstrated in Figure 5C and Figure 5—figure supplement 3, where the free energy cost for tearing apart two nucleosomes with 167 bp DNA in the presence of linker histones (blue) is significantly higher than the curve for bare nucleosomes (red). Notably, at larger inter-nucleosome distances, the values even exceed those for 147 bp nucleosomes (black). A closer examination of the simulation configurations suggests that the disordered C-terminal tail of linker histones can extend and bind the DNA from the second nucleosome, thereby stabilizing the inter-nucleosomal contacts (as shown in Figure 5D). Our results are consistent with prior studies that underscore the importance of linker histones in chromatin compaction (Finch and Klug, 1976; Zhou et al., 2021), particularly in eukaryotic cells with longer linker DNA (Routh et al., 2008; Dombrowski et al., 2022).

Discussion

We introduced a residue-level coarse-grained model with explicit ions to accurately account for electrostatic contributions to chromatin organization. The model achieves quantitative accuracy in reproducing experimental values for the binding affinity of protein-DNA complexes, the energetics of nucleosomal DNA unwinding, nucleosome binding free energy, and the sedimentation coefficients of nucleosome arrays. It captures the counterion atmosphere around the nucleosome core particle as seen in all-atom simulations (Materese et al., 2009) and highlights the contribution of counterions to nucleosome stability. The coarse-grained model also succeeds in resolving the difference between monovalent and divalent ions, supporting the efficacy of divalent ions in neutralizing negative charges and offsetting repulsive interactions among DNA segments.

One significant finding from our study is the predicted strong inter-nucleosome interactions under the physiological salt environment, reaching approximately 9 $k_{B} T$ . We showed that the much lower value reported in a previous DNA-origami experiment is due to the restricted nucleosomal orientation inherent to the device design. Unrestricted nucleosomes allow more close contacts to stabilize binding. A significant nucleosome binding free energy also agrees with the high forces found in single-molecule pulling experiments that are needed for chromatin unfolding (Kruithof et al., 2009; Meng et al., 2015; Kaczmarczyk et al., 2020). We also demonstrate that this strong inter-nucleosomal interaction is largely preserved at longer NRL in the presence of linker histone proteins. While post-translational modifications of histone proteins may influence inter-nucleosomal interactions, their effects are limited, as indicated by Ding et al. (Ding et al., 2021), and are unlikely to completely abolish the significant interactions reported here. Therefore, we anticipate that, in addition to molecular motors, chromatin regulators, and other molecules inside the nucleus, intrinsic inter-nucleosome interactions are important players in chromatin organization in vivo.

We focused our study on single chromatin chains. Strong inter-nucleosome interactions support the compaction and stacking of chromatin, promoting the formation of fibril-like structures. However, as shown in many studies (Maeshima et al., 2016; Ricci et al., 2015; Ou et al., 2017; Zhang et al., 2022), such fibril configurations can hardly be detected in vivo. It is worth emphasizing that this lack of fibril configurations does not contradict our conclusion on the significance of inter-nucleosome interactions. In a prior paper, we found that many in vivo factors, most notably crowding, could disrupt fibril configurations in favor of inter-chain contacts (Liu et al., 2022). The inter-chain contacts can indeed be driven by favorable inter-nucleosome interactions.

Several aspects of the coarse-grained model presented here can be further improved. For instance, the introduction of specific protein-DNA interactions could help address the differences in non-bonded interactions between amino acids and nucleotides beyond electrostatics (Lin et al., 2021a). Such a modification would enhance the model’s accuracy in predicting interactions between chromatin and chromatin proteins. Additionally, the single-bead-per-amino-acid representation used in this study encounters challenges when attempting to capture the influence of histone modifications, which are known to be prevalent in native nucleosomes. Multiscale simulation approaches may be necessary (Collepardo-Guevara et al., 2015). One could first assess the impact of these modifications on the conformation of disordered histone tails using atomistic simulations. By incorporating these conformational changes into the coarse-grained model, systematic investigations of histone modifications on nucleosome interactions and chromatin organization can be conducted. Such a strategy may eventually enable the direct quantification of interactions among native nucleosomes and even the prediction of chromatin organization in vivo.

Materials and methods

Coarse-grained modeling of chromatin

Request a detailed protocol

The large system size of chromatin and the slow timescale for its conformational relaxation necessitates coarse-grained modeling. Following previous studies (Leicher et al., 2020; Ding et al., 2021; Lin et al., 2021a; Lin et al., 2021b; Liu et al., 2022), we adopted a residue-level coarse-grained model for efficient simulations of chromatin. The structure-based model (Clementi et al., 2000; Noel et al., 2016) was applied to represent the histone proteins with one bead per amino acid and to preserve the tertiary structure of the folded regions. The disordered histone tails were kept flexible without tertiary structure biases. A sequence-specific potential, in the form of the Lennard-Jones (LJ) potential and with the strength determined from the Miyazwa-Jernigan (MJ) potential (Miyazawa and Jernigan, 1985), was added to describe the interactions between amino acids. The 3SPN.2C model was adopted to represent each nucleotide with three beads and interactions among DNA beads follow the potential outlined in Freeman et al., 2014, except that the charge of each phosphate site was switched from –0.6 to –1.0 to account for the presence of explicit ions. The Coulombic potential was applied between charged protein and DNA particles. In addition, a weak, non-specific LJ potential was used to account for the excluded volume effect among all protein-DNA beads. Detail expressions for protein-protein and protein-DNA interaction potentials can be found in Ding et al., 2021, and the Appendix section ‘Coarse-grained protein-DNA model’.

We observe that residue-level coarse-grained models have been extensively utilized in prior studies to examine the free energy penalty associated with nucleosomal DNA unwinding (Lequieu et al., 2016; Parsons and Zhang, 2019; Zhang et al., 2016), sequence-dependent nucleosome sliding (Lequieu et al., 2017; Brandani et al., 2018), binding free energy between two nucleosomes (Moller et al., 2019), chromatin folding (Ding et al., 2021; Liu et al., 2022), the impact of histone modifications on tri-nucleosome structures (Chang and Takada, 2016), and protein-chromatin interactions (Watanabe et al., 2018; Leicher et al., 2020). The frequent quantitative agreement between simulation and experimental results supports the utility of such models in chromatin studies. Our introduction of explicit ions, as detailed in Appendix section ‘Coarse-grained explicit ion model’, further extends the applicability of these models to explore the dependence of chromatin conformations on salt concentrations.

Coarse-grained modeling of counterions

Request a detailed protocol

Explicit particle-based representations for monovalent and divalent ions are needed to accurately account for electrostatic interactions (Freeman et al., 2011; Hinckley and de Pablo, 2015; Hayes et al., 2015; Denesyuk and Thirumalai, 2015; Denesyuk et al., 2018; Wang et al., 2022; Sun et al., 2022). We followed Freeman et al., 2011, to introduce explicit ions (see Figure 1) and adopted their potentials to describe the interactions between ions and nucleotide particles, with detailed expressions provided in the Appendix section ‘Coarse-grained explicit ion model’. Parameters in these potentials were tuned by Freeman et al., 2011, to reproduce the radial distribution functions and the potential of mean force between ion pairs determined from all-atom simulations.

This explicit ion model was originally introduced for nucleic acid simulations. We generalized the model for protein simulations by approximating the interactions between charged amino acids and ions with parameters tuned for phosphate sites. Parameter values for ion-amino acid interactions are provided in Table 1 and Table 2.

Table 1

Summary of parameters used to describe interactions between ions and charged particles.

See text section ‘Coarse-grained explicit ion model’ for definitions of various parameters.

Coarse-grained pair	$ϵ$ (kcal/mol)	$σ$ (Å)	$r_{m ϵ}$ (Å)	$σ_{ϵ}$ (Å)	$H_{1}$ (kcal/mol)	$r_{m h 1}$ (Å)	$σ_{h 1}$ (Å)	$H_{2}$ (kcal/mol)	$r_{m h 2}$ (Å)	$σ_{h 2}$ (Å)
P-P	0.18379	6.86	6.86	0.5	–	–	–	–	–	–
Na⁺-P	0.02510	4.14	3.44	1.25	3.15488	4.1	0.57	0.47801	6.5	0.4
Na⁺-AA^+*	0.239	4.065	3.44	1.25	3.15488	4.1	0.57	–	–	–
Na⁺-AA^−†	0.239	4.065	3.44	1.25	3.15488	4.1	0.57	0.47801	6.5	0.4
Mg²⁺-P	0.1195	4.87	3.75	1.0	1.29063	6.1	0.5	0.97992	8.3	1.2
Mg²⁺-AA⁺	0.239	3.556	3.75	1.0	1.29063	6.1	0.5	–	–	–
Mg²⁺-AA⁻	0.239	3.556	3.75	1.0	1.29063	6.1	0.5	0.97992	8.3	1.2
Cl⁻-P	0.08121	5.5425	4.2	0.5	0.83652	6.7	1.5	–	–	–
Cl⁻-AA⁺	0.239	4.8725	4.2	0.5	0.83652	6.7	1.5	0.47801	5.6	0.4
Cl⁻-AA⁻	0.239	4.8725	4.2	0.5	0.83652	6.7	1.5	–	–	–
Na⁺-Na⁺	0.01121	2.43	2.7	0.57	0.17925	5.8	0.57	–	–	–
Na⁺-Mg²⁺	0.04971	2.37	2.37	0.5	–	–	–	–	–	–
Na⁺-Cl⁻	0.08387	3.1352	3.9	2.06	5.49713	3.3	0.57	0.47801	5.6	0.4
Mg²⁺-Mg²⁺	0.89460	1.412	1.412	0.5	–	–	–	–	–	–
Mg²⁺-Cl⁻	0.49737	4.74	4.48	0.57	1.09943	5.48	0.44	0.05975	8.16	0.35
Cl⁻-Cl⁻	0.03585	4.045	4.2	0.56	0.23901	6.2	0.5	–	–	–

*

Positive amino acids.
†

Negative amino acids.

Table 2

Summary of parameters used to describe the WCA interactions between ions and neutral particles.

See text section ‘Coarse-grained explicit ion model’ for definitions of various parameters.

Coarse-grained pair	$ϵ$ (kcal/mol)	$σ$ (Å)
Na⁺-S*	0.239	4.315
Na⁺-A^†	0.239	3.915
Na⁺-T^‡	0.239	4.765
Na⁺-G^§	0.239	3.665
Na⁺-C^¶	0.239	4.415
Na⁺-AA**	0.239	4.065
Mg²⁺-S	0.239	3.806
Mg²⁺-A	0.239	3.406
Mg²⁺-T	0.239	4.256
Mg²⁺-G	0.239	3.156
Mg²⁺-C	0.239	3.906
Mg²⁺-AA**	0.239	3.556
Cl⁻-S	0.239	5.1225
Cl⁻-A	0.239	4.7225
Cl⁻-T	0.239	5.5725
Cl⁻-G	0.239	4.4725
Cl⁻-C	0.239	5.2225
Cl⁻-AA**	0.239	4.8725

*

Sugar.
†

Adenine base.
‡

Thymine base.
§

Guanine base.
¶

Cytosine base.
**

Non-charged amino acids.

Details of molecular dynamics simulations

Request a detailed protocol

We simulated various chromatin systems, including a single-nucleosome, two-nucleosomes, and a 12-mer nucleosome array. The initial configurations for the molecular dynamics simulations were constructed based on the crystal structure of a single nucleosome with PDB ID: 1KX5 (Davey et al., 2002) and 3LZ1 (Vasudevan et al., 2010), or a tetranucleosome with PDB ID: 1ZBB (Schalch et al., 2005). We used the 3DNA software (Lu and Olson, 2003) to add additional DNA, connect and align nucleosomes, and extend the chain length as necessary. Further details on constructing the initial configurations are provided in the Appendix section ‘Ionic dependence of the conformation for a 12-mer nucleosomal array’. Chromatin was positioned at the center of a cubic box with a length selected to avoid interactions between nucleosomes and their periodic images. Counterions were added on a uniformly spaced grid to achieve the desired salt concentrations and neutralize the system. The number of ions and the size of simulation boxes are provided in Table 3.

Table 3

Summary of simulation setups used in this study.

Additional simulation details can be found in text section ‘Molecular dynamics simulation details’.

Studies	Box size (nm³)	Number of Na⁺	Number of Mg²⁺	Number of Cl⁻
Single nucleosome 100 mM NaCl+0.5 mM MgCl₂	216,000	13,017	65	13,003
Twelve nucleosomes 5 mM NaCl	1,331,000	6196	0	4006
Twelve nucleosomes 150 mM NaCl	216,000	21,695	0	19,505
Twelve nucleosomes 0.6 mM MgCl₂	3,375,000	0	2314	2438
Twelve nucleosomes 1 mM MgCl₂	3,375,000	0	3127	4064
Two 147 bp 601-seq nucleosomes 35 mM NaCl+11 mM MgCl₂	125,000	2922	828	4290
Two 147 bp 601-seq nucleosomes 150 mM NaCl+2 mM MgCl₂	216,000	19,505	260	19,737
Two 147 bp poly-dA:dT nucleosomes 150 mM NaCl+2 mM MgCl₂	216,000	19,505	260	19,737
Two 147 bp poly-dG:dC nucleosomes 150 mM NaCl+2 mM MgCl₂	216,000	19,505	260	19,737
Two 167 bp 601-seq nucleosomes 150 mM NaCl+2 mM MgCl₂	216,000	19,505	260	19,657
Two 167 bp 601-seq nucleosomes with H1.0 150 mM NaCl+2 mM MgCl₂	216,000	19,505	260	19,763

All simulations were performed at constant temperature and constant volume (NVT) using the software package LAMMPS (Plimpton, 1995). The electrostatic interactions were implemented with the particle-particle particle-mesh solver, with the relative root-mean-square error in per-atom force set to 0.0001 (Hockney and Eastwood, 2021). A Nosé-Hoover style algorithm (Shinoda et al., 2004) was used to maintain the system temperature at 300 K with a damping parameter of 1 ps. We further modeled the histone core and the inner layer of the nucleosomal DNA together as a rigid body to improve computational efficiency. This approximation does not affect the thermodynamic properties of chromatin (Ding et al., 2021; Liu et al., 2022). Umbrella simulations were used to enhance the sampling of the conformational space (Torrie and Valleau, 1977), and details of the collective variables employed in these simulations are provided in the Appendix section ‘Molecular dynamics simulation details’. All the results presented in the main text are reweighted from the biased simulations by the weighted histogram algorithm (Kumar et al., 1992).

Appendix 1

Coarse-grained protein-DNA model

The force fields describing protein-protein, protein-DNA, and DNA-DNA interactions followed previous studies (Ding et al., 2021; Liu et al., 2022; Freeman et al., 2014; Zhang et al., 2016). Detailed expressions for the potential energies can be found in these references. Specifically, for DNA-DNA interactions, we followed the approach described in Freeman et al., 2011, and with the parameters updated to the latest version of the DNA model 3SPN.2C (Freeman et al., 2014).

Protein-protein interactions include structure-based terms extracted from the initial configuration and generic terms for specific amino acid interactions. We first generated the bonded and non-bonded structure-based interactions within histone proteins using the Shadow algorithm (Noel et al., 2012) implemented by the SMOG software package (Noel et al., 2016). We further scaled the non-bonded interaction strength (Noel et al., 2016) by 2.5 to prevent proteins from unfolding at 300 K. Interactions between histone proteins from different nucleosomes were described using the MJ potential (Miyazawa and Jernigan, 1985) scaled by a factor of 0.4. We have shown in our previous studies that the scaled MJ potential gives a balanced modeling of the radius of gyration for both ordered and disordered proteins (Ding et al., 2021).

Protein-DNA interactions include Coulombic interactions and the excluded volume effect. Unlike the previous Debye-Hückel treatment of electrostatic interactions in an implicit-solvent environment, we modeled the electrostatics between proteins and DNA using the Coulombic potential

U_{e l e c} = \frac{1}{4 π ϵ_{0}} \frac{q_{i} q_{j}}{ϵ_{0} r},

where $ϵ_{0} = 78.0$ is the dielectric constant of the bulk solvent. $q_{i}$ and $q_{j}$ correspond to the charges of the two particles. The excluded volume effect was modeled using the WCA potential with the following form

U_{W C A} = {\begin{cases} 4 ϵ^{P D} [(\frac{σ}{r})^{12} - (\frac{σ}{r})^{6}] + ϵ^{P D} & r < r_{c u t} \\ 0 & r > r_{c u t} . \end{cases}

The cutoff distance $r_{c u t}$ was set to $2^{\frac{1}{6}} σ$ , with $σ = 5.7$ Å. The interactive strength $ϵ^{P D}$ was set as 0.02987572 kcal/mol. More details of this potential can be found in Zhang et al., 2016.

Coarse-grained explicit ion model

Following Freeman et al., 2011, we adopted three terms to describe interactions between charged particles and ions: the Coulombic potential for electrostatic interactions, the Gaussian potential for the hydration effect, and the LJ potential for the excluded volume effect. Thus,

U = U_{e l e c} + U_{h y d r} + U_{L J} .

The electrostatic potential

U_{e l e c} = \frac{1}{4 π ϵ_{0}} \frac{q_{i} q_{j}}{ϵ_{D} (r) r},

where $ϵ_{D} (r)$ is a distance-dependent dielectric constant given in the form

ϵ_{D} (r) = (\frac{5.2 + ϵ_{s}}{2}) + (\frac{5.2 + ϵ_{s}}{2}) \tanh [\frac{r - r_{m ϵ}}{σ_{ϵ}}] .

$ϵ_{s} = 78.0$ is the dielectric constant of the bulk solvent, and values for $σ_{ϵ}$ and $r_{m ϵ}$ are provided in Table 1. The distance cutoff of this potential is set at 20.0 Å. Electrostatic interactions outside this cutoff are computed in reciprocal space.

The hydration potential

U_{h y d r} = \frac{H}{σ_{h} \sqrt{2 π}} \exp [- \frac{(r - r_{m h})^{2}}{2 σ_{h}^{2}}] .

$r_{m h}$ , $σ_{h}$ , and $H$ represent the midpoint, the width, and the height of the hydration shell, respectively, and their ion-specific values are provided in Table 1. Any pair of ions experiences one hydration potential defined above. For pairs formed with ions of distinct types, a second hydration potential with a different set of parameters is applied, with parameters provided in Table 1. The distance cutoff of this potential is set at 12.0 Å.

The LJ potential is given by

U_{L J} = 4 ϵ [(\frac{σ}{r})^{12} - (\frac{σ}{r})^{6}] .

The ion-specific values of $ϵ$ and $σ$ are given in Table 1, and the distance cutoff of this potential is set at 12.0 Å.

The interactions between neutral particles and ions are described by the WCA potential

U_{e x c l} = {\begin{cases} 4 ϵ_{e x c l} [(\frac{σ}{r})^{12} - (\frac{σ}{r})^{6}] + ϵ_{e x c l} & r < r_{c u t} \\ 0 & r > r_{c u t} . \end{cases}

$r_{c u t} = 2^{\frac{1}{6}} σ$ is located at the minimum of the corresponding LJ potential. Values for $σ$ and $ϵ_{e x c l}$ follow the parameters given by Freeman et al., 2011, and are presented in Table 2.

Molecular dynamics simulation details

All simulations were carried out using the software Lammps (Plimpton, 1995) with the force fields defined in the previous two sections. Umbrella sampling simulations (Torrie and Valleau, 1977) were performed using the Plumed software package (The PLUMED consortium, 2019). We used the weighted histogram analysis method (Kumar et al., 1992) implemented by the SMOG software package (Noel et al., 2016) to process the simulation data and compute the free energy profiles.

Binding free energy of protein-DNA complexes

We carried out a series of umbrella-sampling simulations to compute the binding free energies of a set of nine protein-DNA complexes with experimentally documented binding dissociation constants (Dragan et al., 2003a; Dragan et al., 2004; Dragan et al., 2003b; Dragan et al., 2006; Privalov et al., 2011). Initial configurations of these simulations were prepared using the crystal structures with the corresponding PDB IDs listed in Figure 2—figure supplement 1.

The simulations were performed under the same experimental conditions of 100 mM monovalent ions. We used a spring constant of 0.01 kcal/mol/Å² to restrain the distance between the geometric centers of protein and DNA. The centers of the umbrella windows were placed on a uniform grid of [0.0:140.0:10.0] Å, and each umbrella trajectory lasts for 7.15 million steps, with a time step of 2.0 fs. We excluded the first 3 million steps when constructing the free energy profile.

Single-nucleosome simulations for DNA unwinding energetics

To study DNA unwinding from a 601-sequence nucleosome, we built the system by combining histone proteins with explicit coordinates for the disordered tails from PDB ID: 1KX5 (Davey et al., 2002) with the DNA structure from PDB ID: 3LZ1 (Vasudevan et al., 2010).

Umbrella simulations with the DNA end-to-end distance as the collective variable was performed to determine the free energy profile. The end-to-end distance was defined as the geometric center distance between the first and last five base pairs. We used a harmonic umbrella potential with a spring constant of 0.001 kcal/(mol·Å²), and the umbrella centers were placed on a uniform grid of [30.0:510.0:30.0] Å. To increase computational efficiency, the histone core proteins and the two nucleotides located on the dyad axis of the nucleosome were rigidified during the simulations. Each umbrella trajectory lasts for 13.65 million steps, with a time step of 10 fs, and we excluded the first 3 million steps when constructing the free energy profile.

The simulation used the same ionic concentration as the experiment, which includes 0.10 M NaCl and 0.5 mM MgCl₂ (Hall et al., 2009). The cubic simulation box size was set to 600 Å. Extra ions were added to neutralize the system. In total, the system includes a total of 13,017 Na⁺, 65 Mg²⁺, and 13,003 Cl⁻ ions.

Ionic dependence of the conformation for a 12-mer nucleosomal array

We constructed a nucleosomal array of 12 nucleosomes with 20 bp linker DNA to study the impact of different ions on the higher-order chromatin organization. Using the protocol outlined in a previous study (Liu et al., 2022), we started with a nucleosome unit extracted from the tetranucleosome crystal structure (PDB ID: 1ZBB) (Schalch et al., 2005). To connect multiple nucleosome units, we left an extra 20 bp linker DNA at the exit site of the nucleosome. We connected 12 nucleosome units to build the 12-mer nucleosomal array, with an additional 20 bp linker DNA at the end. This 20 bp extra linker of the last nucleosome unit was removed to complete the system setup. To ensure complete histone proteins with disordered tails, we replaced the histone proteins of the nucleosome units with those from the crystal structure with PDB ID: 1KX5 (Davey et al., 2002).

To enhance conformational sampling, we performed umbrella simulations with the collective variable $Q$ defined as

Q = \frac{1}{N} \sum_{i}^{N} \exp (- \frac{(r_{i} - d_{0}^{i})^{2}}{2 r_{0}^{2}})

$i$ enumerates all the nucleosome pairs in the system and $r_{i}$ is the distance between the ith pair. N=66 is the number of nucleosome pairs, and $r_{0}$ = 20.0 Å. $d_{0}^{i}$ corresponds to the distance between the ith pair of nucleosomes determined from the reference two-start structure. $Q$ measures the similarity of a given 12-mer configuration to the reference two-start structure, with larger values representing higher similarity. The reference two-start fibril structure was built in our previous study (Liu et al., 2022) by aligning the structure with a template generated by the software fiberModel (Koslover et al., 2010). We placed the umbrella centers at [0.40:0.90,0.1] and used a spring constant of 50.0 kcal/mol in the harmonic potentials. Each umbrella trajectory lasted 7.8 million steps, with a time step of 10 fs. We used the last 4.8 million steps to construct free energy profiles and compute ensemble averages.

We simulated the nucleosome array under four ionic conditions for comparison with experimental measurements (Correll et al., 2012). The 12-mer was placed in the center of a cubic box, and counterions were added to reach the desired concentration. Excess ions were also introduced to ensure the net neutrality of the system. Specific values of the simulation box size and counterion numbers are provided in Table 3 for reference.

Binding free energy between nucleosomes

Simulations at high salt concentrations

To determine the inter-nucleosome interactions and compare them with the DNA-origami-based experiment (Funke et al., 2016), we simulated a pair of 601-sequence nucleosomes. The nucleosomes were built in the same way as the single-nucleosome DNA unwrapping simulation. The simulation box size was 500 Å, and extra ions were added to neutralize the system. The system comprised a total of 2922 Na⁺, 828 Mg²⁺, and 4290 Cl⁻ ions.

We defined the internal coordinate system for each nucleosome to control their relative orientations as follows. The center of each nucleosome was determined using the geometric center of the list of residues: 63–120, 165–217, 263–324, 398–462, 550–607, 652–704, 750–811, and 885–949. The IDs continuously index residues from chain A to chain H of the crystal structure with PDB ID: 1KX5. To define the nucleosome plane, we chose another point based on the geometric center of the nucleosome dyad site (residue ID: 81–131, 568–618). The unit vector pointing from the nucleosome center to the dyad site center is denoted as $\vec{v}$ . We further introduced another unit vector $\vec{u}$ in the nucleosome plane that points from the nucleosome center to a point defined as the geometric center of residues 63–120, 165–217, 750–811, and 885–949. Finally, the unit vector perpendicular to the nucleosome plane, $\vec{w}$ , is determined as the cross product $\vec{u} \times \vec{v}$ . We use ${\vec{w}}_{1}$ and ${\vec{w}}_{2}$ to differentiate the two nucleosomes.

We utilized two collective variables for the system without constraints to perform the umbrella simulations. The first collective variable measures the distance $r$ between the two nucleosome centers. The second collective variable corresponds to the angle $θ$ between the two unit vectors, ${\vec{w}}_{1}$ and ${\vec{w}}_{2}$ . The umbrella centers were placed on a uniform grid of [60.0, 130.0, 10.0] Å × [0.0, 180.0, 30.0] degrees. The spring constants of the harmonic potentials are 0.01 kcal/(mol · Å²) and 0.001 kcal/(mol · ○²).

For the system that mimics the DNA-origami experiments, we imposed several spatial restraints such that the nucleosomes unbind along a specific pathway (Figure 4—figure supplement 1). As detailed below, the restraints ensure that the first nucleosome is fixed on the X-Y plane while the second nucleosome moves along an arc 15 nm away from the origin.

We introduced three virtual sites, denoted as $O, A$ , and $B$ , with Cartesian coordinates as [0, 0, 0], [150, 0, 0], and [0, 150, 0] Å, respectively. The vectors $\vec{O A}$ and $\vec{O B}$ define the X-Y plane. We further denote the centers of the two nucleosomes as $C_{1}$ and $C_{2}$ .
The first nucleosome was restrained at site B using a harmonic potential with a spring constant of 100 kcal/(mol · Å²). In addition, to mimic its attachment to the bottom arm of the DNA origami, we forced this nucleosome to be parallel to the X-Y plane. Specifically, we restrained the angles between ${\vec{w}}_{1}$ and $\vec{O A}$ or $\vec{O B}$ to be 90°. The spring constant of these harmonic restraints was set to 100 kcal/(mol · ○²).
To mimic the attachment of the second nucleosome to the upper arm of the origami, we restrained the distance between $C_{2}$ and site O as 150 Å with a harmonic potential. The spring constant of this potential was set to 100.0 kcal/(mol · Å²). In addition, we ensured that the second nucleosome is parallel to the plane formed by the vector $\vec{O A}$ and the vector connecting site O to $C_{2}$ . Two harmonic potentials were applied on the angles between ${\vec{w}}_{2}$ and $\vec{O A}$ or $\vec{O C_{2}}$ to restrict them to 90°. The spring constant of these restraints was set to $100 k c a l / (m o l \cdot \circ^{2})$ . We further restricted the second nucleosome to only move in the Y-Z plane by biasing the angle between $\vec{O C_{2}}$ and $\vec{O A}$ to 90° with a spring constant of 100 kcal/(mol · ○²).
Finally, we ensured that the dyad axes of the two nucleosomes in our system are at an angle of 78°, as done experimentally (Funke et al., 2016), by applying a harmonic potential on the angle between the $\vec{v}$ of the two nucleosomes. The spring constant of this potential was set to 100 kcal/(mol · ○²).

We used the angle $θ$ between the two vectors ${\vec{w}}_{1}$ and ${\vec{w}}_{2}$ as the collective variable for umbrella simulations. The umbrella centers were placed on a uniform grid of [0.0, 110.0, 5.0] degrees. The spring constant of the harmonic potential was set as 0.01 kcal/(mol · ○²). Each umbrella simulation lasted 13 million steps and a time step of 10 fs. The first 3 million steps of the simulation were discarded as equilibration.

Simulations at the physiological salt concentration

We performed a series of simulations under the physiological salt concentration, i.e., 150 mM NaCl and 2 mM MgCl₂ for nucleosomes with different repeat lengths and DNA sequences. The 601-sequence nucleosomes were built in the same way as the single-nucleosome simulations. To explore the effect of DNA sequences on inter-nucleosomal interactions, we replaced the original nucleosomal DNA with poly-dA:dT and poly-dG:dC sequences.

To investigate the effect of linker DNA, we simulated nucleosomes with a repeat length of 167 bp. We added 10 base pairs of poly-dA:dT sequences on each side of the existing 147 bp 601-nucleosomal DNA using the software X3DNA (Lu and Olson, 2003). Specifically, we generated an 11-base-pair linker DNA of poly-dA:dT sequence. The additional DNA base pair was created to align the linker DNA with the existing nucleosomal DNA. This alignment was performed such that this additional base pair overlapped with the nucleosomal DNA’s first or last base pair, fixing the linker DNA’s orientation. Finally, we deleted the additional base pair of DNA after the alignment.

Additionally, we built nucleosomes with linker histones using a recently resolved chromatosome structure through cryoelectron microscopy (Zhou et al., 2021). The experimentally determined structure (PDB ID: 7K5X) includes a 197 bp 601-sequence nucleosome with the globular domain of H1.0. As the disordered regions of the linker histone were not resolved in the structure, we modeled them based on the protein sequence using the software Modeller (Eswar et al., 2006) and connected the modeled structure to the globular domain. We then replaced the histone proteins with those from the PDB ID: 1KX5 to provide explicit coordinates for the histone tails. Only the central 167 bp of DNA was kept to build a system with 10 bp linker DNA. The globular domain of H1.0 was bound to the nucleosome dyad and simulated with the histone core protein as a rigid body for computational efficiency.

The numbers of ions and box sizes in each simulation are provided in Table 3. We employed the same two collective variables as the unrestrained simulation at high salt concentrations to conduct umbrella simulations. The umbrella centers were placed on a uniform grid of [60.0, 130.0, 10.0] Å × [0.0, 180.0, 30.0] degrees. The spring constants of the harmonic potentials are 0.01 kcal/(mol · Å²) and 0.001 kcal/(mol · ○²). Each simulation lasted 13 million steps with a time step of 10 fs. We excluded the first 3 million steps when constructing the free energy profiles.

Details of simulation analysis

Number of ions bound to DNA and histone proteins

To calculate the number of ions bound to the nucleosomal DNA and histone proteins, we used the COORDINATIONNUMBER command available in the Plumed (The PLUMED consortium, 2019) software package.

For example, for every Na⁺, we computed the coordinate number as $CN = \sum_{i} s (r_{i})$ , where $i$ loops over all coarse-grained DNA sites and $r_{i}$ is the distance between the ion and the ith DNA bead. $s (r)$ is a switching function defined as

s (r) = \frac{1 - (\frac{r - d_{0}}{r_{0}})^{n}}{1 - (\frac{r - d_{0}}{r_{0}})^{m}},

where $d_{0} = 0.0$ , $r_{0} = 10.0$ , $n = 15$ , and $m = 30$ . An ion with a coordination number greater than 1 was considered bound to DNA. We followed the same procedure to calculate the number of ions bound to histone proteins.

The calculations were performed using the open-source, community-developed PLUMED library (The PLUMED consortium, 2019), version 2.4 (Tribello et al., 2014; The PLUMED consortium, 2019).

Number of unwrapped DNA base pairs

We computed the number of unwrapped DNA base pairs using a similar procedure to the one used for calculating the number of bound ions.

First, we computed a coordination number for each DNA base pair to determine whether it was bound to the histone core. The coordinate number was defined as $CN = \sum_{i} \sum_{j} s (r_{i, j})$ , where $i$ loops over all coarse-grained sites of the corresponding DNA base pair and $j$ loops over all coarse-grained sites of the histone core. $s (r)$ is defined in Equation 10 with $d_{0} = 0.0$ , $r_{0} = 8.0$ , $n = 15$ , and $m = 30$ . A DNA base pair with CN greater than 1 was considered bound to histone proteins. As the histone core is not a perfect cylinder, there were several continuous regions of bound DNA interspersed by regions of unbound DNA. To avoid ambiguity, we defined the wrapped base pairs, $N_{w r a p p e d}$ , as those between the first and last bound base pairs. Correspondingly, the number of unwrapped base pairs was $N_{u n w r a p p e d} = 147 - N_{w r a p p e d}$ .

We set $r_{0} = 8.0$ Å when computing the switching function. At larger values for $r_{0}$ , we found that the calculated numbers overestimate the unwrapped base pairs, as seen from visual inspection of the structures (Figure 2, Appendix 1—figure 1).

Sedimentation coefficients of nucleosome arrays

We calculated the sedimentation coefficients for the 12-mer nucleosome array using the HullRad method (Fleming and Fleming, 2018) with the following equation

s = 10^{8} (\frac{M - M \bar{v} ρ_{20, w}}{N_{A} 6 π η_{0} R_{T}}) .

$M$ is the molar mass of the molecule, $N_{A}$ is Avogadro’s number, and $\bar{v}$ is the partial specific volume. $ρ_{20, w}$ is the density of water at 20°C, and $η_{0}$ is the water viscosity at 20°C. R_T is the translational hydrodynamic radius calculated based on the convex hull of the target biomolecule.

Estimation of error bars

We estimated the error bars of the 12-mer simulations based on the standard deviation calculated from the probability distribution of the variables (Figure 3—figure supplement 1), i.e.,

σ (X) = \sqrt{E [X^{2}] - (E [X])^{2}}

where $E [X]$ is the expected value of $X$ .

We divided the trajectories into three equal-length partitions for all the other simulations and computed the free energy profiles independently. The error bars were estimated as the standard deviation of the three independent estimates.

Appendix 1—figure 1

Download asset Open asset

A cutoff value of 8.0 Å produces more accurate values for the number of unwrapped DNA base pairs as determined from visual inspection of representative configurations, related to Figure 2 of the main text.

See text section ‘Number of unwrapped DNA base pairs’ for additional discussions. A typical nucleosome structure with most of the outer layer DNA unwrapped was used to examine the impact of different cutoff values. The histone core is colored in gold, with histone tails in white, the wrapped DNA in blue, and the unwrapped DNA in red. The discrepancy among various cutoff values is evident in the highlight regions enclosed by dotted circles. As shown in the zoom-ins in the middle panel, a cutoff of 10 Å results in three additional base pairs of DNA detected as wrapped in I (highlighted in orange square). However, these extra base pairs not detected with a cutoff of 8 Å are visibly detached from histone proteins. Similarly, 9 Å and 10 Å cutoff values result in five extra base pairs of DNA detected as wrapped in II (highlighted in orange squares in the bottom panel).

Data availability

The current manuscript is a computational study, so no data have been generated for this manuscript. Modelling code is uploaded to GitHub: https://github.com/ZhangGroup-MITChemistry/Explicit_Ion_Chromatin (copy archived at Lin, 2023).

References

(2015) Chromatin compaction under mixed salt conditions: opposite effects of sodium and potassium ions on nucleosome array folding
Scientific Reports 5:8512.

https://doi.org/10.1038/srep08512
- PubMed
- Google Scholar
1. Bajpai G
2. Amiad Pavlov D
3. Lorber D
4. Volk T
5. Safran S
(2021) Mesoscale phase separation of chromatin in the nucleus
eLife 10:e63976.

https://doi.org/10.7554/eLife.63976
- PubMed
- Google Scholar
1. Bancaud A
2. Wagner G
3. Conde E Silva N
4. Lavelle C
5. Wong H
6. Mozziconacci J
7. Barbi M
8. Sivolob A
9. Le Cam E
10. Mouawad L
11. Viovy J-L
12. Victor J-M
13. Prunell A
(2007) Nucleosome chiral transition under positive torsional stress in single chromatin fibers
Molecular Cell 27:135–147.

https://doi.org/10.1016/j.molcel.2007.05.037
- PubMed
- Google Scholar
Book
1. Bascom G
2. Schlick T
(2018) Nuclear Architecture and Dynamics
Elsevier.

https://doi.org/10.1016/C2014-0-03248-0
- Google Scholar
1. Bennink ML
2. Leuba SH
3. Leno GH
4. Zlatanova J
5. de Grooth BG
6. Greve J
(2001) Unfolding individual nucleosomes by stretching single chromatin fibers with optical tweezers
Nature Structural Biology 8:606–610.

https://doi.org/10.1038/89646
- PubMed
- Google Scholar
1. Bonev B
2. Cavalli G
(2016) Organization and function of the 3D genome
Nature Reviews. Genetics 17:661–678.

https://doi.org/10.1038/nrg.2016.112
- PubMed
- Google Scholar
(2022) Shaping the genome via lengthwise compaction, phase separation, and lamina adhesion
Nucleic Acids Research 50:4258–4271.

https://doi.org/10.1093/nar/gkac231
- PubMed
- Google Scholar
Preprint
(2023) Temporally Correlated Active Forces Drive Chromosome Structure and Dynamics
bioRxiv.

https://doi.org/10.1101/2023.04.23.528410
- Google Scholar
1. Brandani GB
2. Niina T
3. Tan C
4. Takada S
(2018) DNA sliding in nucleosomes via twist defect propagation revealed by molecular simulations
Nucleic Acids Research 46:2788–2801.

https://doi.org/10.1093/nar/gky158
- PubMed
- Google Scholar
1. Cai S
2. Böck D
3. Pilhofer M
4. Gan L
(2018) The in situ structures of mono-, di-, and trinucleosomes in human heterochromatin
Molecular Biology of the Cell 29:2450–2457.

https://doi.org/10.1091/mbc.E18-05-0331
- PubMed
- Google Scholar
1. Chang L
2. Takada S
(2016) Histone acetylation dependent energy landscapes in tri-nucleosome revealed by residue-resolved molecular simulations
Scientific Reports 6:34441.

https://doi.org/10.1038/srep34441
- PubMed
- Google Scholar
1. Clark DJ
2. Kimura T
(1990) Electrostatic mechanism of chromatin folding
Journal of Molecular Biology 211:883–896.

https://doi.org/10.1016/0022-2836(90)90081-V
- PubMed
- Google Scholar
(2000) Topological and energetic factors: what determines the structural details of the transition state ensemble and “en-route” intermediates for protein folding? An investigation for small globular proteins
Journal of Molecular Biology 298:937–953.

https://doi.org/10.1006/jmbi.2000.3693
- PubMed
- Google Scholar
(2015) Chromatin unfolding by epigenetic modifications explained by dramatic impairment of internucleosome interactions: A multiscale computational study
Journal of the American Chemical Society 137:10205–10215.

https://doi.org/10.1021/jacs.5b04086
- PubMed
- Google Scholar
(2012) Short nucleosome repeats impose rotational modulations on chromatin fibre folding
The EMBO Journal 31:2416–2426.

https://doi.org/10.1038/emboj.2012.80
- PubMed
- Google Scholar
1. Cui Y
2. Bustamante C
(2000) Pulling a single chromatin fiber reveals the forces that maintain its higher-order structure
PNAS 97:127–132.

https://doi.org/10.1073/pnas.97.1.127
- PubMed
- Google Scholar
1. Davey CA
2. Sargent DF
3. Luger K
4. Maeder AW
5. Richmond TJ
(2002) Solvent mediated interactions in the structure of the nucleosome core particle at 1.9 a resolution
Journal of Molecular Biology 319:1097–1113.

https://doi.org/10.1016/S0022-2836(02)00386-8
- PubMed
- Google Scholar
(2013) Exploring the three-dimensional organization of genomes: interpreting chromatin interaction data
Nature Reviews. Genetics 14:390–403.

https://doi.org/10.1038/nrg3454
- PubMed
- Google Scholar
1. Denesyuk NA
2. Thirumalai D
(2015) How do metal ions direct ribozyme folding?
Nature Chemistry 7:793–801.

https://doi.org/10.1038/nchem.2330
- PubMed
- Google Scholar
(2018) Molecular Simulations of Ion Effects on the Thermodynamics of RNA Folding
The Journal of Physical Chemistry. B 122:11860–11867.

https://doi.org/10.1021/acs.jpcb.8b08142
- PubMed
- Google Scholar
1. Ding X
2. Lin X
3. Zhang B
(2021) Stability and folding pathways of tetra-nucleosome from six-dimensional free energy surface
Nature Communications 12:1091.

https://doi.org/10.1038/s41467-021-21377-z
- PubMed
- Google Scholar
(2016) Transferable model for chromosome architecture
PNAS 113:12168–12173.

https://doi.org/10.1073/pnas.1613607113
- PubMed
- Google Scholar
(2022) Histone H1 binding to nucleosome arrays depends on linker DNA length and trajectory
Nature Structural & Molecular Biology 29:493–501.

https://doi.org/10.1038/s41594-022-00768-w
- PubMed
- Google Scholar
(2003a) DNA binding of a non-sequence-specific HMG-D protein is entropy driven with a substantial non-electrostatic contribution
Journal of Molecular Biology 331:795–813.

https://doi.org/10.1016/s0022-2836(03)00785-x
- PubMed
- Google Scholar
(2003b) The energetics of specific binding of AT-hooks from HMGA1 to target DNA
Journal of Molecular Biology 327:393–411.

https://doi.org/10.1016/s0022-2836(03)00050-0
- PubMed
- Google Scholar
(2004) DNA binding and bending by HMG boxes: energetic determinants of specificity
Journal of Molecular Biology 343:371–393.

https://doi.org/10.1016/j.jmb.2004.08.035
- PubMed
- Google Scholar
1. Dragan AI
2. Li Z
3. Makeyeva EN
4. Milgotina EI
5. Liu Y
6. Crane-Robinson C
7. Privalov PL
(2006) Forces driving the binding of homeodomains to DNA
Biochemistry 45:141–151.

https://doi.org/10.1021/bi051705m
- PubMed
- Google Scholar
1. Engelhardt M
(2004) Condensation of chromatin in situ by cation-dependent charge shielding and aggregation
Biochemical and Biophysical Research Communications 324:1210–1214.

https://doi.org/10.1016/j.bbrc.2004.09.175
- PubMed
- Google Scholar
1. Eswar N
2. Webb B
3. Marti-Renom MA
4. Madhusudhan MS
5. Eramian D
6. Shen MY
7. Pieper U
8. Sali A
(2006) Comparative protein structure modeling using Modeller
Current Protocols in Bioinformatics Chapter 5:bi0506s15.

https://doi.org/10.1002/0471250953.bi0506s15
- PubMed
- Google Scholar
1. Falk M
2. Feodorova Y
3. Naumova N
4. Imakaev M
5. Lajoie BR
6. Leonhardt H
7. Joffe B
8. Dekker J
9. Fudenberg G
10. Solovei I
11. Mirny LA
(2019) Heterochromatin drives compartmentalization of inverted and conventional nuclei
Nature 570:395–399.

https://doi.org/10.1038/s41586-019-1275-3
- PubMed
- Google Scholar
(2021) Nucleosome plasticity is a critical element of chromatin liquid-liquid phase separation and multivalent nucleosome interactions
Nature Communications 12:2883.

https://doi.org/10.1038/s41467-021-23090-3
- PubMed
- Google Scholar
(2010) Charge state of the globular histone core controls stability of the nucleosome
Biophysical Journal 99:1577–1585.

https://doi.org/10.1016/j.bpj.2010.06.046
- PubMed
- Google Scholar
(2018) Modulation of nucleosomal DNA accessibility via charge-altering post-translational modifications in histone core
Epigenetics & Chromatin 11:11.

https://doi.org/10.1186/s13072-018-0181-5
- PubMed
- Google Scholar
1. Finch JT
2. Klug A
(1976) Solenoidal model for superstructure in chromatin
PNAS 73:1897–1901.

https://doi.org/10.1073/pnas.73.6.1897
- PubMed
- Google Scholar
1. Finn EH
2. Misteli T
(2019) Molecular basis and biological function of variability in spatial genome organization
Science 365:eaaw9498.

https://doi.org/10.1126/science.aaw9498
- PubMed
- Google Scholar
1. Fleming PJ
2. Fleming KG
(2018) Hullrad: Fast calculations of folded and disordered protein and nucleic acid hydrodynamic properties
Biophysical Journal 114:856–869.

https://doi.org/10.1016/j.bpj.2018.01.002
- PubMed
- Google Scholar
1. Forties RA
2. North JA
3. Javaid S
4. Tabbaa OP
5. Fishel R
6. Poirier MG
7. Bundschuh R
(2011) A quantitative model of nucleosome dynamics
Nucleic Acids Research 39:8306–8313.

https://doi.org/10.1093/nar/gkr422
- PubMed
- Google Scholar
(2011) A coarse-grain three-site-per-nucleotide model for DNA with explicit ions
The Journal of Chemical Physics 135:165104.

https://doi.org/10.1063/1.3652956
- PubMed
- Google Scholar
(2014) Coarse-grained modeling of DNA curvature
The Journal of Chemical Physics 141:165103.

https://doi.org/10.1063/1.4897649
- PubMed
- Google Scholar
1. Fujishiro S
2. Sasai M
(2022) Generation of dynamic three-dimensional genome structure through phase separation of chromatin
PNAS 119:e2109838119.

https://doi.org/10.1073/pnas.2109838119
- PubMed
- Google Scholar
1. Funke JJ
2. Ketterer P
3. Lieleg C
4. Schunter S
5. Korber P
6. Dietz H
(2016) Uncovering the forces between nucleosomes using DNA origami
Science Advances 2:e1600974.

https://doi.org/10.1126/sciadv.1600974
- PubMed
- Google Scholar
1. Gibcus JH
2. Dekker J
(2013) The hierarchy of the 3D genome
Molecular Cell 49:773–782.

https://doi.org/10.1016/j.molcel.2013.02.011
- PubMed
- Google Scholar
1. Gibson BA
2. Doolittle LK
3. Schneider MWG
4. Jensen LE
5. Gamarra N
6. Henry L
7. Gerlich DW
8. Redding S
9. Rosen MK
(2019) Organization of chromatin by intrinsic and regulated phase separation
Cell 179:470–484.

https://doi.org/10.1016/j.cell.2019.08.037
- PubMed
- Google Scholar
(2023) Polymer folding through active processes recreates features of genome organization
PNAS 120:e2221726120.

https://doi.org/10.1073/pnas.2221726120
- PubMed
- Google Scholar
(2009) Evidence for heteromorphic chromatin fibers from analysis of nucleosome interactions
PNAS 106:13317–13322.

https://doi.org/10.1073/pnas.0903280106
- PubMed
- Google Scholar
1. Hall MA
2. Shundrovsky A
3. Bai L
4. Fulbright RM
5. Lis JT
6. Wang MD
(2009) High-resolution dynamic mapping of histone-DNA interactions in a nucleosome
Nature Structural & Molecular Biology 16:124–129.

https://doi.org/10.1038/nsmb.1526
- PubMed
- Google Scholar
1. Hayes RL
2. Noel JK
3. Mandic A
4. Whitford PC
5. Sanbonmatsu KY
6. Mohanty U
7. Onuchic JN
(2015) Generalized Manning Condensation Model Captures the RNA Ion Atmosphere
Physical Review Letters 114:258105.

https://doi.org/10.1103/PhysRevLett.114.258105
- PubMed
- Google Scholar
(2021) Transcription organizes euchromatin via microphase separation
Nature Communications 12:1360.

https://doi.org/10.1038/s41467-021-21589-3
- PubMed
- Google Scholar
1. Hinckley DM
2. de Pablo JJ
(2015) Coarse-grained ions for nucleic acid modeling
Journal of Chemical Theory and Computation 11:5436–5446.

https://doi.org/10.1021/acs.jctc.5b00341
- PubMed
- Google Scholar
Book
1. Hockney RW
2. Eastwood JW
(2021) Computer Simulation Using Particles
CRC Press.

https://doi.org/10.1201/9780367806934
- Google Scholar
1. Hu M
2. Deng K
3. Qin Z
4. Dixon J
5. Selvaraj S
6. Fang J
7. Ren B
8. Liu JS
9. Tanay A
(2013) Bayesian inference of spatial organizations of chromosomes
PLOS Computational Biology 9:e1002893.

https://doi.org/10.1371/journal.pcbi.1002893
- PubMed
- Google Scholar
1. Huang Y-C
2. Su C-J
3. Korolev N
4. Berezhnoy NV
5. Wang S
6. Soman A
7. Chen C-Y
8. Chen H-L
9. Jeng U-S
10. Nordenskiöld L
(2018) The effect of linker DNA on the structure and interaction of nucleosome core particles
Soft Matter 14:9096–9106.

https://doi.org/10.1039/c8sm00998h
- PubMed
- Google Scholar
(2016) Implicit solvent model for million-atom atomistic simulations: Insights into the organization of 30-nm chromatin fiber
Journal of Chemical Theory and Computation 12:5946–5959.

https://doi.org/10.1021/acs.jctc.6b00712
- PubMed
- Google Scholar
1. Jiang Z
2. Qi Y
3. Kamat K
4. Zhang B
(2022) Phase separation and correlated motions in motorized genome
The Journal of Physical Chemistry. B 126:5619–5628.

https://doi.org/10.1021/acs.jpcb.2c03238
- PubMed
- Google Scholar
(2014) Modeling epigenome folding: formation and dynamics of topologically associated chromatin domains
Nucleic Acids Research 42:9553–9561.

https://doi.org/10.1093/nar/gku698
- PubMed
- Google Scholar
1. Kaczmarczyk A
2. Meng H
3. Ordu O
4. van Noort J
5. Dekker NH
(2020) Chromatin fibers stabilize nucleosomes under torsional stress
Nature Communications 11:126.

https://doi.org/10.1038/s41467-019-13891-y
- PubMed
- Google Scholar
1. Kamat K
2. Lao Z
3. Qi Y
4. Wang Y
5. Ma J
6. Zhang B
(2023) Compartmentalization with nuclear landmarks yields random, yet precise, genome organization
Biophysical Journal 122:1376–1389.

https://doi.org/10.1016/j.bpj.2023.03.003
- PubMed
- Google Scholar
1. Kent S
2. Brown K
3. Yang CH
4. Alsaihati N
5. Tian C
6. Wang H
7. Ren X
(2020) Phase-separated transcriptional condensates accelerate target-search process revealed by live-cell single-molecule imaging
Cell Reports 33:108248.

https://doi.org/10.1016/j.celrep.2020.108248
- PubMed
- Google Scholar
(2010) Local geometry and elasticity in compact chromatin structure
Biophysical Journal 99:3941–3950.

https://doi.org/10.1016/j.bpj.2010.10.024
- PubMed
- Google Scholar
1. Kruithof M
2. Chien F-T
3. Routh A
4. Logie C
5. Rhodes D
6. van Noort J
(2009) Single-molecule force spectroscopy reveals a highly compliant helical folding for the 30-nm chromatin fiber
Nature Structural & Molecular Biology 16:534–540.

https://doi.org/10.1038/nsmb.1590
- PubMed
- Google Scholar
(1992) The weighted histogram analysis method for free‐energy calculations on biomolecules. I. The method
Journal of Computational Chemistry 13:1011–1021.

https://doi.org/10.1002/jcc.540130812
- Google Scholar
(2020) Mesoscale liquid model of chromatin recapitulates nuclear order of eukaryotes
Biophysical Journal 118:2130–2140.

https://doi.org/10.1016/j.bpj.2019.09.013
- Google Scholar
1. Larson AG
2. Elnatan D
3. Keenen MM
4. Trnka MJ
5. Johnston JB
6. Burlingame AL
7. Agard DA
8. Redding S
9. Narlikar GJ
(2017) Liquid droplet formation by HP1α suggests a role for phase separation in heterochromatin
Nature 547:236–240.

https://doi.org/10.1038/nature22822
- PubMed
- Google Scholar
1. Latham AP
2. Zhang B
(2021) Consistent Force Field Captures Homologue-Resolved HP1 Phase Separation
Journal of Chemical Theory and Computation 17:3134–3144.

https://doi.org/10.1021/acs.jctc.0c01220
- PubMed
- Google Scholar
1. Leicher R
2. Ge EJ
3. Lin X
4. Reynolds MJ
5. Xie W
6. Walz T
7. Zhang B
8. Muir TW
9. Liu S
(2020) Single-molecule and in silico dissection of the interaction between Polycomb repressive complex 2 and chromatin
PNAS 117:30465–30475.

https://doi.org/10.1073/pnas.2003395117
- PubMed
- Google Scholar
1. Leicher R
2. Osunsade A
3. Chua GNL
4. Faulkner SC
5. Latham AP
6. Watters JW
7. Nguyen T
8. Beckwitt EC
9. Christodoulou-Rubalcava S
10. Young PG
11. Zhang B
12. David Y
13. Liu S
(2022) Single-stranded nucleic acid binding and coacervation by linker histone H1
Nature Structural & Molecular Biology 29:463–471.

https://doi.org/10.1038/s41594-022-00760-4
- PubMed
- Google Scholar
(2016) Tension-dependent free energies of nucleosome unwrapping
ACS Central Science 2:660–666.

https://doi.org/10.1021/acscentsci.6b00201
- PubMed
- Google Scholar
(2017) In silico evidence for sequence-dependent nucleosome sliding
PNAS 114:E9197–E9205.

https://doi.org/10.1073/pnas.1705685114
- PubMed
- Google Scholar
1. Lesne A
2. Riposo J
3. Roger P
4. Cournac A
5. Mozziconacci J
(2014) 3D genome reconstruction from chromosomal contacts
Nature Methods 11:1141–1143.

https://doi.org/10.1038/nmeth.3104
- Google Scholar
1. Li S
2. Wei T
3. Panchenko AR
(2023) Histone variant H2A.Z modulates nucleosome dynamics to promote DNA accessibility
Nature Communications 14:769.

https://doi.org/10.1038/s41467-023-36465-5
- PubMed
- Google Scholar
1. Lieberman-Aiden E
2. van Berkum NL
3. Williams L
4. Imakaev M
5. Ragoczy T
6. Telling A
7. Amit I
8. Lajoie BR
9. Sabo PJ
10. Dorschner MO
11. Sandstrom R
12. Bernstein B
13. Bender MA
14. Groudine M
15. Gnirke A
16. Stamatoyannopoulos J
17. Mirny LA
18. Lander ES
19. Dekker J
(2009) Comprehensive mapping of long-range interactions reveals folding principles of the human genome
Science 326:289–293.

https://doi.org/10.1126/science.1181369
- PubMed
- Google Scholar
1. Lin X
2. Leicher R
3. Liu S
4. Zhang B
(2021a) Cooperative DNA looping by PRC2 complexes
Nucleic Acids Research 49:6238–6248.

https://doi.org/10.1093/nar/gkab441
- PubMed
- Google Scholar
1. Lin X
2. Qi Y
3. Latham AP
4. Zhang B
(2021b) Multiscale modeling of genome organization with maximum entropy optimization
The Journal of Chemical Physics 155:010901.

https://doi.org/10.1063/5.0044150
- PubMed
- Google Scholar
Software
1. Lin XC
(2023) Explicit_Ion_Chromatin, version swh:1:rev:20839f7ae8ddb3a990132983c100bbf495c8a920
Software Heritage.

https://archive.softwareheritage.org/swh:1:dir:edda538171b11d7f444c2ed8049f2d19db2f2c2e;origin=https://github.com/ZhangGroup-MITChemistry/Explicit_Ion_Chromatin;visit=swh:1:snp:9d80d45d093e17cff8e8f32aad91f7a7d3757802;anchor=swh:1:rev:20839f7ae8ddb3a990132983c100bbf495c8a920
1. Liu S
2. Lin X
3. Zhang B
(2022) Chromatin fiber breaks into clutches under tension and crowding
Nucleic Acids Research 50:9738–9747.

https://doi.org/10.1093/nar/gkac725
- PubMed
- Google Scholar
1. Lowary PT
2. Widom J
(1998) New DNA sequence rules for high affinity binding to histone octamer and sequence-directed nucleosome positioning
Journal of Molecular Biology 276:19–42.

https://doi.org/10.1006/jmbi.1997.1494
- PubMed
- Google Scholar
1. Lu X-J
2. Olson WK
(2003) 3DNA: a software package for the analysis, rebuilding and visualization of three-dimensional nucleic acid structures
Nucleic Acids Research 31:5108–5121.

https://doi.org/10.1093/nar/gkg680
- PubMed
- Google Scholar
(2018) Bottom-up modeling of chromatin segregation due to epigenetic modifications
PNAS 115:12739–12744.

https://doi.org/10.1073/pnas.1812268115
- PubMed
- Google Scholar
1. Maeshima K
2. Rogge R
3. Tamura S
4. Joti Y
5. Hikima T
6. Szerlong H
7. Krause C
8. Herman J
9. Seidel E
10. DeLuca J
11. Ishikawa T
12. Hansen JC
(2016) Nucleosomal arrays self-assemble into supramolecular globular structures lacking 30-nm fibers
The EMBO Journal 35:1115–1132.

https://doi.org/10.15252/embj.201592660
- PubMed
- Google Scholar
(2002) Salt-induced conformation and interaction changes of nucleosome core particles
Biophysical Journal 82:345–356.

https://doi.org/10.1016/S0006-3495(02)75399-X
- PubMed
- Google Scholar
1. Manning GS
(1978) The molecular theory of polyelectrolyte solutions with applications to the electrostatic properties of polynucleotides
Quarterly Reviews of Biophysics 11:179–246.

https://doi.org/10.1017/s0033583500002031
- PubMed
- Google Scholar
(2009) Counterion atmosphere and hydration patterns near a nucleosome core particle
Journal of the American Chemical Society 131:15005–15013.

https://doi.org/10.1021/ja905376q
- PubMed
- Google Scholar
(2015) Quantitative analysis of single-molecule force spectroscopy on folded chromatin fibers
Nucleic Acids Research 43:3578–3590.

https://doi.org/10.1093/nar/gkv215
- PubMed
- Google Scholar
1. Mirny L
2. Dekker J
(2022) Mechanisms of chromosome folding and nuclear organization: Their interplay and open questions
Cold Spring Harbor Perspectives in Biology 14:a040147.

https://doi.org/10.1101/cshperspect.a040147
- PubMed
- Google Scholar
1. Misteli T
(2020) The self-organizing genome: Principles of genome architecture and function
Cell 183:28–45.

https://doi.org/10.1016/j.cell.2020.09.014
- PubMed
- Google Scholar
1. Miyazawa S
2. Jernigan RL
(1985) Estimation of effective interresidue contact energies from protein crystal structures: quasi-chemical approximation
Macromolecules 18:534–552.

https://doi.org/10.1021/ma00145a039
- Google Scholar
(2019) The free energy landscape of internucleosome interactions and its relation to chromatin fiber structure
ACS Central Science 5:341–348.

https://doi.org/10.1021/acscentsci.8b00836
- PubMed
- Google Scholar
(2012) The shadow map: A general contact definition for capturing the dynamics of biomolecular folding and function
The Journal of Physical Chemistry. B 116:8692–8702.

https://doi.org/10.1021/jp300852d
- PubMed
- Google Scholar
1. Noel JK
2. Levi M
3. Raghunathan M
4. Lammert H
5. Hayes RL
6. Onuchic JN
7. Whitford PC
(2016) SMOG 2: A Versatile Software Package for Generating Structure-Based Models
PLOS Computational Biology 12:e1004794.

https://doi.org/10.1371/journal.pcbi.1004794
- PubMed
- Google Scholar
1. Ortiz V
2. de Pablo JJ
(2011) Molecular origins of DNA flexibility: sequence effects on conformational and mechanical properties
Physical Review Letters 106:238107.

https://doi.org/10.1103/PhysRevLett.106.238107
- PubMed
- Google Scholar
1. Ou HD
2. Phan S
3. Deerinck TJ
4. Thor A
5. Ellisman MH
6. O’Shea CC
(2017) ChromEMT: Visualizing 3D chromatin structure and compaction in interphase and mitotic cells
Science 357:eaag0025.

https://doi.org/10.1126/science.aag0025
- PubMed
- Google Scholar
1. Parsons T
2. Zhang B
(2019) Critical role of histone tail entropy in nucleosome unwinding
The Journal of Chemical Physics 150:185103.

https://doi.org/10.1063/1.5085663
- Google Scholar
Book
1. Phillips R
(2012)
Physical Biology of the Cell

Garland Science.
- Google Scholar
1. Plimpton S
(1995) Fast parallel algorithms for short-range molecular dynamics
Journal of Computational Physics 117:1–19.

https://doi.org/10.1006/jcph.1995.1039
- Google Scholar
(2005) Single chromatin fiber stretching reveals physically distinct populations of disassembly events
Biophysical Journal 88:3572–3583.

https://doi.org/10.1529/biophysj.104.053074
- PubMed
- Google Scholar
(2011) Interpreting protein/DNA interactions: distinguishing specific from non-specific and electrostatic from non-electrostatic components
Nucleic Acids Research 39:2483–2491.

https://doi.org/10.1093/nar/gkq984
- PubMed
- Google Scholar
1. Quinodoz SA
2. Ollikainen N
3. Tabak B
4. Palla A
5. Schmidt JM
6. Detmar E
7. Lai MM
8. Shishkin AA
9. Bhat P
10. Takei Y
11. Trinh V
12. Aznauryan E
13. Russell P
14. Cheng C
15. Jovanovic M
16. Chow A
17. Cai L
18. McDonel P
19. Garber M
20. Guttman M
(2018) Higher-Order Inter-chromosomal Hubs Shape 3D Genome Organization in the Nucleus
Cell 174:744–757.

https://doi.org/10.1016/j.cell.2018.05.024
- PubMed
- Google Scholar
(2015) Chromatin fibers are formed by heterogeneous groups of nucleosomes in vivo
Cell 160:1145–1158.

https://doi.org/10.1016/j.cell.2015.01.054
- PubMed
- Google Scholar
(2008) Nucleosome repeat length and linker histone stoichiometry determine chromatin fiber structure
PNAS 105:8872–8877.

https://doi.org/10.1073/pnas.0802336105
- PubMed
- Google Scholar
(2005) X-ray structure of a tetranucleosome and its implications for the chromatin fibre
Nature 436:138–141.

https://doi.org/10.1038/nature03686
- PubMed
- Google Scholar
1. Schiessel H
(2003) The physics of chromatin
Journal of Physics 15:R699–R774.

https://doi.org/10.1088/0953-8984/15/19/203
- Google Scholar
(1996) Reversible oligonucleosome self-association: dependence on divalent cations and core histone tail domains
Biochemistry 35:4009–4015.

https://doi.org/10.1021/bi9525684
- PubMed
- Google Scholar
1. Shi G
2. Thirumalai D
(2021) From Hi-C Contact Map to Three-Dimensional Organization of Interphase Human Chromosomes
Physical Review X 11:011051.

https://doi.org/10.1103/PhysRevX.11.011051
- Google Scholar
(2004) Rapid estimation of elastic constants by molecular dynamics simulation under constant stress
Physical Review B 69:134103.

https://doi.org/10.1103/PhysRevB.69.134103
- Google Scholar
(2020) Condensed chromatin behaves like a solid on the mesoscale in vitro and in living cells
Cell 183:1772–1784.

https://doi.org/10.1016/j.cell.2020.11.027
- PubMed
- Google Scholar
1. Su J-H
2. Zheng P
3. Kinrot SS
4. Bintu B
5. Zhuang X
(2020) Genome-Scale Imaging of the 3D Organization and Transcriptional Activity of Chromatin
Cell 182:1641–1659.

https://doi.org/10.1016/j.cell.2020.07.032
- PubMed
- Google Scholar
(2022) A bottom-up coarse-grained model for nucleosome-nucleosome interactions with explicit ions
Journal of Chemical Theory and Computation 18:3948–3960.

https://doi.org/10.1021/acs.jctc.2c00083
- PubMed
- Google Scholar
1. Takei Y
2. Yun J
3. Zheng S
4. Ollikainen N
5. Pierson N
6. White J
7. Shah S
8. Thomassie J
9. Suo S
10. Eng C-HL
11. Guttman M
12. Yuan G-C
13. Cai L
(2021) Integrated spatial genomics reveals global architecture of single nuclei
Nature 590:344–350.

https://doi.org/10.1038/s41586-020-03126-2
- PubMed
- Google Scholar
1. The PLUMED consortium
(2019) Promoting transparency and reproducibility in enhanced molecular simulations
Nature Methods 16:670–673.

https://doi.org/10.1038/s41592-019-0506-8
- Google Scholar
1. Torrie GM
2. Valleau JP
(1977) Nonphysical sampling distributions in Monte Carlo free-energy estimation: Umbrella sampling
Journal of Computational Physics 23:187–199.

https://doi.org/10.1016/0021-9991(77)90121-8
- Google Scholar
(2014) PLUMED 2: New feathers for an old bird
Computer Physics Communications 185:604–613.

https://doi.org/10.1016/j.cpc.2013.09.018
- Google Scholar
(2010) Crystal structures of nucleosome core particles containing the “601” strong positioning sequence
Journal of Molecular Biology 403:1–10.

https://doi.org/10.1016/j.jmb.2010.08.039
- PubMed
- Google Scholar
1. Wang A
2. Levi M
3. Mohanty U
4. Whitford PC
(2022) Diffuse ions coordinate dynamics in a ribonucleoprotein assembly
Journal of the American Chemical Society 144:9510–9522.

https://doi.org/10.1021/jacs.2c04082
- PubMed
- Google Scholar
1. Watanabe S
2. Mishima Y
3. Shimizu M
4. Suetake I
5. Takada S
(2018) Interactions of HP1 Bound to H3K9me3 Dinucleosome by Molecular Simulations and Biochemical Assays
Biophysical Journal 114:2336–2351.

https://doi.org/10.1016/j.bpj.2018.03.025
- PubMed
- Google Scholar
1. Widom J
(1986) Physicochemical studies of the folding of the 100 A nucleosome filament into the 300 A filament:cation dependence
Journal of Molecular Biology 190:411–424.

https://doi.org/10.1016/0022-2836(86)90012-4
- PubMed
- Google Scholar
(2015) Shearing of the CENP-A dimerization interface mediates plasticity in the octameric centromeric nucleosome
Scientific Reports 5:17038.

https://doi.org/10.1038/srep17038
- PubMed
- Google Scholar
(2021) The dynamic influence of linker histone saturation within the poly-nucleosome array
Journal of Molecular Biology 433:166902.

https://doi.org/10.1016/j.jmb.2021.166902
- PubMed
- Google Scholar
1. Xie WJ
2. Meng L
3. Liu S
4. Zhang L
5. Cai X
6. Gao YQ
(2017) Structural modeling of chromatin integrates genome features and reveals chromosome folding principle
Scientific Reports 7:2818.

https://doi.org/10.1038/s41598-017-02923-6
- PubMed
- Google Scholar
1. Xie WJ
2. Zhang B
(2019) Learning the formation mechanism of domain-level chromatin states with epigenomics data
Biophysical Journal 116:2047–2056.

https://doi.org/10.1016/j.bpj.2019.04.006
- PubMed
- Google Scholar
1. Xie L
2. Dong P
3. Qi Y
4. Hsieh T-HS
5. English BP
6. Jung S
7. Chen X
8. De Marzio M
9. Casellas R
10. Chang HY
11. Zhang B
12. Tjian R
13. Liu Z
(2022) BRD2 compartmentalizes the accessible genome
Nature Genetics 54:481–491.

https://doi.org/10.1038/s41588-022-01044-9
- PubMed
- Google Scholar
1. Yildirim A
2. Feig M
(2018) High-resolution 3D models of Caulobacter crescentus chromosome reveal genome structural variability and organization
Nucleic Acids Research 46:3937–3952.

https://doi.org/10.1093/nar/gky141
- PubMed
- Google Scholar
1. Zhang B
2. Zheng W
3. Papoian GA
4. Wolynes PG
(2016) Exploring the free energy landscape of nucleosomes
Journal of the American Chemical Society 138:8126–8133.

https://doi.org/10.1021/jacs.6b02893
- PubMed
- Google Scholar
1. Zhang M
2. Díaz-Celis C
3. Onoa B
4. Cañari-Chumpitaz C
5. Requejo KI
6. Liu J
7. Vien M
8. Nogales E
9. Ren G
10. Bustamante C
(2022) Molecular organization of the early stages of nucleosome phase separation visualized by cryo-electron tomography
Molecular Cell 82:3000–3014.

https://doi.org/10.1016/j.molcel.2022.06.032
- PubMed
- Google Scholar
1. Zhou B-R
2. Feng H
3. Kale S
4. Fox T
5. Khant H
6. de Val N
7. Ghirlando R
8. Panchenko AR
9. Bai Y
(2021) Distinct structures and dynamics of chromatosomes with different human linker histone isoforms
Molecular Cell 81:166–182.

https://doi.org/10.1016/j.molcel.2020.10.038
- PubMed
- Google Scholar

Article and author information

Author details

Xingcheng Lin

Department of Chemistry, Massachusetts Institute of Technology, Cambridge, United States
Present address
1. Department of Physics, North Carolina State University, Raleigh, United States
2. Bioinformatics Research Center, North Carolina State University, Raleigh, United States
Contribution
Conceptualization, Software, Formal analysis, Validation, Investigation, Visualization, Methodology, Writing – original draft, Writing – review and editing

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0002-9378-6174
Bin Zhang

Department of Chemistry, Massachusetts Institute of Technology, Cambridge, United States

Contribution
Conceptualization, Software, Formal analysis, Supervision, Funding acquisition, Validation, Investigation, Visualization, Methodology, Writing – original draft, Project administration, Writing – review and editing

For correspondence
binz@mit.edu

Competing interests
Reviewing editor, eLife

"This ORCID iD identifies the author of this article:" 0000-0002-3685-7503

Funding

National Institute of General Medical Sciences (R35GM133580)

Bin Zhang

National Science Foundation (MCB-2042362)

Bin Zhang

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

This work was supported by the National Institutes of Health (Grant R35GM133580) and the National Science Foundation (Grant MCB-2042362).

Version history

Preprint posted: May 18, 2023
Sent for peer review: June 26, 2023
Reviewed Preprint version 1: August 21, 2023
Reviewed Preprint version 2: January 10, 2024
Version of Record published: January 30, 2024
Version of Record updated: July 5, 2024

Cite all versions

You can cite all versions using the DOI https://doi.org/10.7554/eLife.90073. This DOI represents all versions, and will always resolve to the latest one.

Copyright

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.

Metrics

889

views
95

downloads
9

citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Mendeley

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

Xingcheng Lin
Bin Zhang

(2024)

Explicit ion modeling predicts physicochemical interactions for chromatin organization

eLife 12:RP90073.

https://doi.org/10.7554/eLife.90073.3

Share this article

Cite this article

Illustration of the residue-level coarse-grained explicit ion model for chromatin simulations.

Explicit ion modeling reproduces the energetics of nucleosomal DNA unwrapping.

Explicit ion modeling predicts salt-dependent conformations of a 12-mer nucleosome array.

Close contacts give rise to strong inter-nucleosomal interactions.

Simulations predict significant inter-nucleosome interactions at physiological conditions.

Summary of parameters used to describe interactions between ions and charged particles.

Summary of parameters used to describe the WCA interactions between ions and neutral particles.

Summary of simulation setups used in this study.

A cutoff value of 8.0 Å produces more accurate values for the number of unwrapped DNA base pairs as determined from visual inspection of representative configurations, related to Figure 2 of the main text.

Author details

Xingcheng Lin

Present address

Contribution

Competing interests

Bin Zhang

Contribution

For correspondence

Competing interests

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

Categories and tags

Further reading