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 kBT . 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 kBT, 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.
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, but some of the claims should be tempered. The tools will be valuable for the chromosome biology field.
Three-dimensional genome organization plays essential roles in numerous DNA-templated processes.1–5 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.6–9 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. 3, 4, 10–12
Compartmentalization has been proposed to arise from the microphase separation of different chromatin types as in block copolymer systems.13–16, 16–26 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,23, 27–32 potentially contributing to chromatin demixing. Demixing can also arise from interactions between chromatin and various nuclear landmarks such as nuclear lamina and speckles,11, 15, 26, 33 as well as active transcriptional processes.34–37 Furthermore, recent studies have revealed that nucleosome arrays alone can undergo spontaneous phase separation,38–40 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.41–44 For example, using reconstituted nucleosome arrays with well-positioning 601-sequence DNA and force-measuring magnetic tweezers, Kruithof et al. estimated the internucleosome binding energy to be ∼ 14 kBT .41 On the other hand, Funke et al. introduced a DNA origami-based force spectrometer to directly probe the interaction between a pair of nucleosomes,44 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 kBT . This large discrepancy in reported values complicates the assessment of the contribution of inter-nucleosomal interactions to chromatin organization in vivo.
Computational modeling is well suited for reconciling the discrepancy across experiments and determining the strength of internucleosome interactions. The high computational cost of atomistic simulations45–47 have inspired several groups to calculate the nucleosome binding free energy with coarse-grained models.48, 49 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 theory50 to describe electrostatic interactions in an implicit-solvent environment,48, 49, 51, 52 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.53–56
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 complexes57 and energetic cost of nucleosomal DNA unwinding.58 Further simulations of chromatin at various salt concentrations reproduce experimentally measured sedimentation coefficients.59 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 kBT . 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.
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.42, 58, 60–62 The DNA-unzipping experiment performed by Hall et al.58 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.63 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 (Fig. 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 counter ion condensation and its impact on electrostatic interactions between protein and DNA molecules. Additional model details are provided in the Materials and Methods and the Supporting Information.
We performed umbrella simulations64 to determine the free energy profile of nucleosomal DNA unwinding. The experimental buffer condition of 0.10M NaCl and 0.5mM MgCl2 58 was adopted in simulations for direct comparison. As shown in Fig. 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 (Fig. S1), supporting the model’s accuracy.
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.65 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. 66, 67 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 Fig. 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 Mg2+ 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,59, 68–72 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.59 that immerses chromatin in solutions with 5mM NaCl, 150mM NaCl, 0.6mM MgCl2, or 1mM MgCl2. 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 Supporting Information. Data from different umbrella windows were combined together with proper reweighting73 for analysis. As shown in Fig. 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 1-3 nucleosomes (Fig. 3B). On the other hand, the compaction is evident at 150 mM NaCl.
Notably, in agreement with previous studies,74–77 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 MgCl2, the chromatin sedimentation coefficient is comparable to that obtained at 150 mM of NaCl.
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 (Fig. 3A bottom) as predicted by the Manning theory.78, 79 Ion condensation weakens the repulsion among DNA segments that prevents chromatin from collapsing. Notably, the fraction of bound Mg2+ 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 Mg2+ and phosphate groups that more effectively offset the entropy loss due to ion condensation.79 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. Accurate quantification is essential for conceptualizing the importance of physicochemical interactions for chromatin organization in vivo, since prior studies have proved inconclusive, producing a wide range of values from 2 to 14 kBT .41, 42, 44 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.44 with 35 mM NaCl and 11 mM MgCl2. 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 (Fig. 4A, A1 to A3). For example, neither nucleosome can freely rotate (Fig. S3); 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 the Supporting Information.
Strikingly, the two sets of simulations produced dramatically different binding free energies. Restricting nucleosome orientations produced a binding free energy of ∼ 2kBT, reproducing the experimental value (Fig. 4B). On the other hand, the binding free energy increased to 15 kBT 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 Fig. 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 (Fig. S4) 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 Fig. 4A), with the angle between the two nucleosome planes close to zero (Fig. S5C). 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 Fig. 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.43 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 MgCl2.
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 (Fig. 5A). As shown in two-dimensional binding free energy landscape of internucleosome distance, r, and θ (Fig. 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 kBT).
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 kBT (Fig. S6). 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,80 which causes the DNA to unwrap more frequently, increasing cross nucleosome contacts at larger distances (Fig. S7).
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,59, 81, 82 increasing the NRL weakened inter-nucleosomal interactions (Fig. 5C and Fig. S8), reducing the binding free energy to ∼ 6 kBT . However, the weakened interactions can be more than compensated with the presence of linker histones (Fig. 5C and Fig. S8). As shown in the example configuration (Fig. 5D), the disordered C-terminal tail of linker histones can extend and bind the DNA from the second nucleosome to stabilize the internucleosomal contacts.
Conclusions and Discussion
We introduced a residual-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 simulations67 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 kBT . 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. It also agrees with the high forces found in single-molecule pulling experiments that are needed for chromatin unfolding.41, 43, 83 We further demonstrate that this strong inter-nucleosomal interaction is mostly preserved at longer NRL with the presence of linker histone proteins. 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,40, 84–86 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.75 The inter-chain contacts can indeed be driven by favorable inter-nucleosome interactions.
Materials and Methods
Coarse-grained modeling of chromatin
The large system size of chromatin and the slow timescale for its conformational relaxation necessitates coarse-grained modeling. Following previous studies,5, 32, 74, 75, 87 we adopted a residue-level coarse-grained model for efficient simulations of chromatin. The structure-based model88, 89 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,90 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 Ref. 91, 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 Ref. 74 and the Supporting Information.
Coarse-grained modeling of counter ions
To accurately account for electrostatic interactions, we included explicit particle-based representations for monovalent and divalent ions (Fig. 1). The ions interact with nucleotide particles following the potential introduced in Ref. 53, with detailed expressions provided in the Supporting Information. Parameters in these potentials were tuned 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 S1.
Details of molecular dynamics simulations
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: 1KX566 and 3LZ1,92 or a tetranucleosome with PDB ID: 1ZBB.93 We used the 3DNA software94 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 Supporting Information. 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 to neutralize the system. The number of ions and the size of simulation boxes are provided in Table S2.
All simulations were performed at constant temperature and constant volume (NVT) using the software package LAMMPS.95 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.96 A Nośe-Hoover style algorithm97 was used to maintain the system temperature at 300K 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.74, 75 Umbrella simulations were used to enhance the sampling of the conformational space,64 and details of the collective variables employed in these simulations are provided in the Supporting Information. All the results presented in the main text are reweighted from the biased simulations by the weighted histogram algorithm.73
This work was supported by the National Institutes of Health (Grant R35GM133580) and the National Science Foundation (Grant MCB-2042362).
- (1)Exploring the Three-Dimensional Organization of Genomes: Interpreting Chromatin Interaction DataNat Rev Genet 14:390–403
- (2)Organization and Function of the 3D GenomeNat Rev Genet 17:661–678
- (3)Molecular Basis and Biological Function of Variability in Spatial Genome OrganizationScience 365
- (4)The Self-Organizing Genome: Principles of Genome Architecture and FunctionCell 183:28–45
- (5)Multiscale Modeling of Genome Organization with Maximum Entropy OptimizationJ. Chem. Phys 155
- (6)Comprehensive Mapping of Long-Range Interactions Reveals Folding Principles of the Human GenomeScience 326:289–293
- (7)Higher-Order Inter-chromosomal Hubs Shape 3D Genome Organization in the NucleusCell 174:744–757
- (8)Genome-Scale Imaging of the 3D Organization and Transcriptional Activity of ChromatinCell 182:1641–1659
- (9)Integrated Spatial Genomics Reveals Global Architecture of Single NucleiNature 590:344–350
- (10)The Hierarchy of the 3D GenomeMolecular Cell 49:773–782
- (11)Mechanisms of Chromosome Folding and Nuclear Organization: Their Interplay and Open QuestionsCold Spring Harb Perspect Biol 14
- (12)Learning the Formation Mechanism of Domain-Level Chromatin States with Epigenomics DataBiophysical Journal 116:2047–2056
- (13)Generation of Dynamic Three-Dimensional Genome Structure through Phase Separation of ChromatinProc. Natl. Acad. Sci. U.S.A 119
- (14)Modeling Epigenome Folding: Formation and Dynamics of Topologically Associated Chromatin DomainsNucleic Acids Res 42:9553–9561
- (15)Heterochromatin Drives Compartmentalization of Inverted and Conventional NucleiNature 570:395–399
- (16)Mesoscale Phase Separation of Chromatin in the NucleuseLife 10
- (17)Mesoscale Liquid Model of Chromatin Recapitulates Nuclear Order of EukaryotesBiophysical Journal 118:2130–2140
- (18)Bayesian Inference of Spatial Organizations of ChromosomesPLoS Comput Biol 9
- (19)3D Genome Reconstruction from Chromosomal ContactsNat Methods 11:1141–1143
- (20)Transferable Model for Chromosome ArchitectureProc. Natl. Acad. Sci. U.S.A 113:12168–12173
- (21)Structural Modeling of Chromatin Integrates Genome Features and Reveals Chromosome Folding PrincipleSci Rep 7
- (22)High-Resolution 3D Models of Caulobacter Crescentus Chromosome Reveal Genome Structural Variability and OrganizationNucleic Acids Res 46:3937–3952
- (23)Bottom–up Modeling of Chromatin Segregation Due to Epigenetic ModificationsProc. Natl. Acad. Sci. U.S.A 115:12739–12744
- (24)From Hi-C Contact Map to Three-Dimensional Organization of Interphase Human ChromosomesPhys. Rev. X 11
- (25)Liquid–Liquid Phase Separation Driven Compartmentalization of Reactive NucleoplasmPhys. Biol 18
- (26)Shaping the Genome via Lengthwise Compaction, Phase Separation, and Lamina AdhesionNucleic Acids Res 50:4258–4271
- (27)Liquid Droplet Formation by HP1α Suggests a Role for Phase Separation in HeterochromatinNature 547:236–240
- (28)Phase-Separated Transcriptional Condensates Accelerate Target-Search Process Revealed by Live-Cell Single-Molecule ImagingCell Reports 33
- (29)BRD2 Compartmentalizes the Accessible GenomeNat Genet 54:481–491
- (30)Single-Stranded Nucleic Acid Binding and Coacervation by Linker Histone H1Nat Struct Mol Biol 29:463–471
- (31)Consistent Force Field Captures Homologue-Resolved HP1 Phase SeparationJ. Chem. Theory Comput 17:3134–3144
- (32)Cooperative DNA Looping by PRC2 ComplexesNucleic Acids Res 49:6238–6248
- (33)Compartmentalization with Nuclear Landmarks Yields Random, yet Precise, Genome OrganizationBiophysical Journal 122:1376–1389
- (34)Transcription Organizes Euchromatin via Microphase SeparationNat Commun 12
- (35)Phase Separation and Correlated Motions in Motorized GenomeJ. Phys. Chem. B 126:5619–5628
- (36)Temporally Correlated Active Forces Drive Chromosome Structure and Dynamics
- (37)Polymer Folding through Active Processes Recreates Features of Genome OrganizationProc. Natl. Acad. Sci. U.S.A 120
- (38)Organization of Chromatin by Intrinsic and Regulated Phase SeparationCell 179:470–484
- (39)Condensed Chromatin Behaves like a Solid on the Mesoscale In Vitro and in Living CellsCell 183:1772–1784
- (40)Molecular Organization of the Early Stages of Nucleosome Phase Separation Visualized by Cryo-Electron TomographyMolecular Cell 82:3000–3014
- (41)Single-Molecule Force Spectroscopy Reveals a Highly Compliant Helical Folding for the 30-Nm Chromatin FiberNat Struct Mol Biol 16:534–540
- (42)Pulling a Single Chromatin Fiber Reveals the Forces That Maintain Its Higher-Order StructureProc. Natl. Acad. Sci. U.S.A 97:127–132
- (43)Chromatin Fibers Stabilize Nucleosomes under Torsional StressNat Commun 11
- (44)Uncovering the Forces between Nucleosomes Using DNA OrigamiSci. Adv 2
- (45)Shearing of the CENP-A Dimerization Interface Mediates Plasticity in the Octameric Centromeric NucleosomeSci Rep 5
- (46)The Dynamic Influence of Linker Histone Saturation within the Poly-Nucleosome ArrayJournal of Molecular Biology 433
- (47)Histone Variant H2A.Z Modulates Nucleosome Dynamics to Promote DNA AccessibilityNat Commun 14
- (48)The Free Energy Landscape of Internucleosome Interactions and Its Relation to Chromatin Fiber StructureACS Cent. Sci 5:341–348
- (49)Nucleosome Plasticity Is a Critical Element of Chromatin Liquid–Liquid Phase Separation and Multivalent Nucleosome InteractionsNat Commun 12
- (50)Physical Biology of the CellGarland Science
- (51)Implicit Solvent Model for Million-Atom Atomistic Simulations: Insights into the Organization of 30-Nm Chromatin FiberJ. Chem. Theory Comput 12:5946–5959
- (52)Nuclear Architecture and Dynamics:123–147
- (53)A Coarse-Grain Three-Site-per-Nucleotide Model for DNA with Explicit IonsThe Journal of Chemical Physics 135
- (54)Coarse-Grained Ions for Nucleic Acid ModelingJ. Chem. Theory Comput 11:5436–5446
- (55)A Bottom-Up Coarse-Grained Model for Nucleosome–Nucleosome Interactions with Explicit IonsJ. Chem. Theory Comput 18:3948–3960
- (56)Generalized Manning Condensation Model Captures the RNA Ion AtmospherePhys. Rev. Lett 114
- (57)Interpreting Protein/DNA Interactions: Distinguishing Specific from Non-Specific and Electrostatic from Non-Electrostatic ComponentsNucleic Acids Research 39:2483–2491
- (58)High-Resolution Dynamic Mapping of Histone-DNA Interactions in a NucleosomeNat Struct Mol Biol 16:124–129
- (59)Short Nucleosome Repeats Impose Rotational Modulations on Chromatin Fibre FoldingEMBO J 31:2416–2426
- (60)Unfolding Individual Nucleosomes by Stretching Single Chromatin Fibers with Optical TweezersNat Struct Biol 8:606–610
- (61)Single Chromatin Fiber Stretching Reveals Physically Distinct Populations of Disassembly EventsBiophysical Journal 88:3572–3583
- (62); Conde e Silva, N.; Lavelle, C.; Wong, H.; Mozziconacci, J.; Barbi, M.; Sivolob, A.; Le Cam, E.; Mouawad, L.; Viovy, J.-L.; Victor, J.-M.; PrunellA. Nucleosome Chiral Transition under Positive Torsional Stress in Single Chromatin Fibers. Molecular Cell 27:135–147
- (63)A Quantitative Model of Nucleosome DynamicsNucleic Acids Res 39:8306–8313
- (64)Nonphysical Sampling Distributions in Monte Carlo Free-Energy Estimation: Umbrella SamplingJournal of Computational Physics 23:187–199
- (65)The Physics of ChromatinJ. Phys.: Condens. Matter 15:R699–R774
- (66)Solvent Mediated Interactions in the Structure of the Nucleosome Core Particle at 1.9Å ResolutionJournal of Molecular Biology 319:1097–1113
- (67)Counterion Atmosphere and Hydration Patterns near a Nucleosome Core ParticleJ. Am. Chem. Soc 131:15005–15013
- (68)Physicochemical Studies of the Folding of the 100 A Nucleosome Filament into the 300 A Filament. Cation DependenceJ Mol Biol 190:411–424
- (69)Reversible Oligonucleosome Self-Association: Dependence on Divalent Cations and Core Histone Tail DomainsBiochemistry 35:4009–4015
- (70)Condensation of Chromatin in Situ by Cation-Dependent Charge Shielding and AggregationBiochem Biophys Res Commun 324:1210–1214
- (71)Evidence for Heteromorphic Chromatin Fibers from Analysis of Nucleosome InteractionsProc. Natl. Acad. Sci. U.S.A 106:13317–13322
- (72)Chromatin Compaction under Mixed Salt Conditions: Opposite Effects of Sodium and Potassium Ions on Nucleosome Array FoldingSci Rep 5
- (73)THE Weighted Histogram Analysis Method for Free-Energy Calculations on BiomoleculesI. The Method. J. Comput. Chem 13:1011–1021
- (74)Stability and Folding Pathways of Tetra-Nucleosome from Six-Dimensional Free Energy SurfaceNat Commun 12
- (75)Chromatin Fiber Breaks into Clutches under Tension and CrowdingNucleic Acids Res 50:9738–9747
- (76)The in Situ Structures of Mono-Di-, and Trinucleosomes in Human Heterochromatin. MBoC 29:2450–2457
- (77)Histone H1 Binding to Nucleosome Arrays Depends on Linker DNA Length and TrajectoryNat Struct Mol Biol 29:493–501
- (78)The Molecular Theory of Polyelectrolyte Solutions with Applications to the Electrostatic Properties of PolynucleotidesQ Rev Biophys 11:179–246
- (79)Electrostatic Mechanism of Chromatin FoldingJ Mol Biol 211:883–896
- (80)Molecular Origins of DNA Flexibility: Sequence Effects on Conformational and Mechanical PropertiesPhys. Rev. Lett 106
- (81)Salt-Induced Conformation and Interaction Changes of Nucleosome Core ParticlesBiophys J 82:345–356
- (82)The Effect of Linker DNA on the Structure and Interaction of Nucleosome Core ParticlesSoft Matter 14:9096–9106
- (83)Quantitative Analysis of Single-Molecule Force Spectroscopy on Folded Chromatin FibersNucleic Acids Res 43:3578–3590
- (84)Nucleosomal Arrays Self-Assemble into Supramolecular Globular Structures Lacking 30-Nm FibersEMBO J 35:1115–1132
- (85)Chromatin Fibers Are Formed by Heterogeneous Groups of Nucleosomes in VivoCell 160:1145–1158
- (86)ChromEMT: Visualizing 3D Chromatin Structure and Compaction in Interphase and Mitotic CellsScience 357
- (87)Single-Molecule and in Silico Dissection of the Interaction between Polycomb Repressive Complex 2 and ChromatinProc. Natl. Acad. Sci. U.S.A 117:30465–30475
- (88)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 ProteinsJournal of Molecular Biology 298:937–953
- (89)SMOG 2: A Versatile Software Package for Generating Structure-Based ModelsPLoS Comput Biol 12
- (90)Estimation of Effective Interresidue Contact Energies from Protein Crystal Structures: Quasi-Chemical ApproximationMacromolecules 18:534–552
- (91)Coarse-Grained Modeling of DNA CurvatureThe Journal of Chemical Physics 141
- (92)Crystal Structures of Nucleosome Core Particles Containing the ’601’ Strong Positioning SequenceJ Mol Biol 403:1–10
- (93)X-Ray Structure of a Tetranu-cleosome and Its Implications for the Chromatin FibreNature 436:138–141
- (94)3DNA: A Software Package for the Analysis, Rebuilding and Visualization of Three-Dimensional Nucleic Acid StructuresNucleic Acids Res 31:5108–5121
- (95)Fast Parallel Algorithms for Short-Range Molecular DynamicsJ. Comput. Phys 117:1–19
- (96)Computer Simulation Using Particles
- (97)Rapid Estimation of Elastic Constants by Molecular Dynamics Simulation under Constant StressPhys. Rev. B 69