Critical role of histone tail entropy in nucleosome unwinding

From Nucleosomes to Compartments: Physicochemical Interactions Underlying Chromatin Organization

Chromatin organization plays a critical role in cellular function by regulating access to genetic information. However, understanding chromatin folding is challenging due to its complex, multiscale nature. Significant progress has been made in studying in vitro systems, uncovering the structure of individual nucleosomes and their arrays, and elucidating the role of physicochemical forces in stabilizing these structures. Additionally, remarkable advancements have been achieved in characterizing chromatin organization in vivo, particularly at the whole-chromosome level, revealing important features such as chromatin loops, topologically associating domains, and nuclear compartments. However, bridging the gap between in vitro and in vivo studies remains challenging. The resemblance between in vitro and in vivo chromatin conformations and the relevance of internucleosomal interactions for chromatin folding in vivo are subjects of debate. This article reviews experimental and computational studies conducted at various length scales, highlighting the significance of intrinsic interactions between nucleosomes and their roles in chromatin folding in vivo.

Dynamical phase transition in models that couple chromatin folding with histone modifications

Genomic regions can acquire heritable epigenetic states through unique histone modifications, which lead to stable gene expression patterns without altering the underlying DNA sequence. However, the relationship between chromatin conformational dynamics and epigenetic stability is poorly understood. In this paper, we propose kinetic models to investigate the dynamic fluctuations of histone modifications and the spatial interactions between nucleosomes. Our model explicitly incorporates the influence of chemical modifications on the structural stability of chromatin and the contribution of chromatin contacts to the cooperative nature of chemical reactions. Through stochastic simulations and analytical theory, we have discovered distinct steady-state outcomes in different kinetic regimes, resembling a dynamical phase transition. Importantly, we have validated that the emergence of this transition, which occurs on biologically relevant timescales, is robust against variations in model design and parameters. Our findings suggest that the viscoelastic properties of chromatin and the timescale at which it transitions from a gel-like to a liquidlike state significantly impact dynamic processes that occur along the one-dimensional DNA sequence.

Coarse-Grained Models to Study Protein-DNA Interactions and Liquid-Liquid Phase Separation

Recent advances in coarse-grained (CG) computational models for DNA have enabled molecular-level insights into the behavior of DNA in complex multiscale systems. However, most existing CG DNA models are not compatible with CG protein models, limiting their applications for emerging topics such as protein-nucleic acid assemblies. Here, we present a new computationally efficient CG DNA model. We first use experimental data to establish the model's ability to predict various aspects of DNA behavior, including melting thermodynamics and relevant local structural properties such as the major and minor grooves. We then employ an all-atom hydropathy scale to define nonbonded interactions between protein and DNA sites, to make our DNA model compatible with an existing CG protein model (HPS-Urry), which is extensively used to study protein phase separation, and show that our new model reasonably reproduces the experimental binding affinity for a prototypical protein-DNA system. To further demonstrate the capabilities of this new model, we simulate a full nucleosome with and without histone tails, on a microsecond time scale, generating conformational ensembles and provide molecular insights into the role of histone tails in influencing the liquid-liquid phase separation (LLPS) of HP1α proteins. We find that histone tails interact favorably with DNA, influencing the conformational ensemble of the DNA and antagonizing the contacts between HP1α and DNA, thus affecting the ability of DNA to promote LLPS of HP1α. These findings shed light on the complex molecular framework that fine-tunes the phase transition properties of heterochromatin proteins and contributes to heterochromatin regulation and function. Overall, the CG DNA model presented here is suitable to facilitate micrometer-scale studies with sub-nm resolution in many biological and engineering applications and can be used to investigate protein-DNA complexes, such as nucleosomes, or LLPS of proteins with DNA, enabling a mechanistic understanding of how molecular information may be propagated at the genome level.

Explicit ion modeling predicts physicochemical interactions for chromatin organization

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.

OpenABC enables flexible, simplified, and efficient GPU accelerated simulations of biomolecular condensates

Biomolecular condensates are important structures in various cellular processes but are challenging to study using traditional experimental techniques. In silico simulations with residue-level coarse-grained models strike a balance between computational efficiency and chemical accuracy. They could offer valuable insights by connecting the emergent properties of these complex systems with molecular sequences. However, existing coarse-grained models often lack easy-to-follow tutorials and are implemented in software that is not optimal for condensate simulations. To address these issues, we introduce OpenABC, a software package that greatly simplifies the setup and execution of coarse-grained condensate simulations with multiple force fields using Python scripting. OpenABC seamlessly integrates with the OpenMM molecular dynamics engine, enabling efficient simulations with performance on a single GPU that rivals the speed achieved by hundreds of CPUs. We also provide tools that convert coarse-grained configurations to all-atom structures for atomistic simulations. We anticipate that OpenABC will significantly facilitate the adoption of in silico simulations by a broader community to investigate the structural and dynamical properties of condensates.

A coarse-grained DNA model to study protein-DNA interactions and liquid-liquid phase separation

Recent advances in coarse-grained (CG) computational models for DNA have enabled molecular-level insights into the behavior of DNA in complex multiscale systems. However, most existing CG DNA models are not compatible with CG protein models, limiting their applications for emerging topics such as protein-nucleic acid assemblies. Here, we present a new computationally efficient CG DNA model. We first use experimental data to establish the model's ability to predict various aspects of DNA behavior, including melting thermodynamics and relevant local structural properties such as the major and minor grooves. We then employ an all-atom hydropathy scale to define non-bonded interactions between protein and DNA sites, to make our DNA model compatible with an existing CG protein model (HPS-Urry), that is extensively used to study protein phase separation, and show that our new model reasonably reproduces the experimental binding affinity for a prototypical protein-DNA system. To further demonstrate the capabilities of this new model, we simulate a full nucleosome with and without histone tails, on a microsecond timescale, generating conformational ensembles and provide molecular insights into the role of histone tails in influencing the liquid-liquid phase separation (LLPS) of HP1α proteins. We find that histone tails interact favorably with DNA, influencing the conformational ensemble of the DNA and antagonizing the contacts between HP1α and DNA, thus affecting the ability of DNA to promote LLPS of HP1α. These findings shed light on the complex molecular framework that fine-tunes the phase transition properties of heterochromatin proteins and contributes to heterochromatin regulation and function. Overall, the CG DNA model presented here is suitable to facilitate micron-scale studies with sub-nm resolution in many biological and engineering applications and can be used to investigate protein-DNA complexes, such as nucleosomes, or LLPS of proteins with DNA, enabling a mechanistic understanding of how molecular information may be propagated at the genome level.

OpenABC Enables Flexible, Simplified, and Efficient GPU Accelerated Simulations of Biomolecular Condensates

Biomolecular condensates are important structures in various cellular processes but are challenging to study using traditional experimental techniques. In silico simulations with residue-level coarse-grained models strike a balance between computational efficiency and chemical accuracy. They could offer valuable insights by connecting the emergent properties of these complex systems with molecular sequences. However, existing coarse-grained models often lack easy-to-follow tutorials and are implemented in software that is not optimal for condensate simulations. To address these issues, we introduce OpenABC, a software package that greatly simplifies the setup and execution of coarse-grained condensate simulations with multiple force fields using Python scripting. OpenABC seamlessly integrates with the OpenMM molecular dynamics engine, enabling efficient simulations with performances on a single GPU that rival the speed achieved by hundreds of CPUs. We also provide tools that convert coarse-grained configurations to all-atom structures for atomistic simulations. We anticipate that Open-ABC will significantly facilitate the adoption of in silico simulations by a broader community to investigate the structural and dynamical properties of condensates. Open-ABC is available at https://github.com/ZhangGroup-MITChemistry/OpenABC.

Chromatin fiber breaks into clutches under tension and crowding

The arrangement of nucleosomes inside chromatin is of extensive interest. While in vitro experiments have revealed the formation of 30 nm fibers, most in vivo studies have failed to confirm their presence in cell nuclei. To reconcile the diverging experimental findings, we characterized chromatin organization using a residue-level coarse-grained model. The computed force–extension curve matches well with measurements from single-molecule experiments. Notably, we found that a dodeca-nucleosome in the two-helix zigzag conformation breaks into structures with nucleosome clutches and a mix of trimers and tetramers under tension. Such unfolded configurations can also be stabilized through trans interactions with other chromatin chains. Our study suggests that unfolding from chromatin fibers could contribute to the irregularity of in vivo chromatin configurations. We further revealed that chromatin segments with fibril or clutch structures engaged in distinct binding modes and discussed the implications of these inter-chain interactions for a potential sol–gel phase transition.

Phase Separation and Correlated Motions in Motorized Genome

The human genome is arranged in the cell nucleus nonrandomly, and phase separation has been proposed as an important driving force for genome organization. However, the cell nucleus is an active system, and the contribution of nonequilibrium activities to phase separation and genome structure and dynamics remains to be explored. We simulated the genome using an energy function parametrized with chromosome conformation capture (Hi-C) data with the presence of active, nondirectional forces that break the detailed balance. We found that active forces that may arise from transcription and chromatin remodeling can dramatically impact the spatial localization of heterochromatin. When applied to euchromatin, active forces can drive heterochromatin to the nuclear envelope and compete with passive interactions among heterochromatin that tend to pull them in opposite directions. Furthermore, active forces induce long-range spatial correlations among genomic loci beyond single chromosome territories. We further showed that the impact of active forces could be understood from the effective temperature defined as the fluctuation-dissipation ratio. Our study suggests that nonequilibrium activities can significantly impact genome structure and dynamics, producing unexpected collective phenomena.

On the stability and layered organization of protein-DNA condensates

Multi-component phase separation is emerging as a key mechanism for the formation of biological condensates that play essential roles in signal sensing and transcriptional regulation. The molecular factors that dictate these condensates' stability and spatial organization are not fully understood, and it remains challenging to predict their microstructures. Using a near-atomistic, chemically accurate force field, we studied the phase behavior of chromatin regulators that are crucial for heterochromatin organization and their interactions with DNA. Our computed phase diagrams recapitulated previous experimental findings on different proteins. They revealed a strong dependence of condensate stability on the protein-DNA mixing ratio as a result of balancing protein-protein interactions and charge neutralization. Notably, a layered organization was observed in condensates formed by mixing HP1, histone H1, and DNA. This layered organization may be of biological relevance, as it enables cooperative DNA packaging between the two chromatin regulators: histone H1 softens the DNA to facilitate the compaction induced by HP1 droplets. Our study supports near-atomistic models as a valuable tool for characterizing the structure and stability of biological condensates.

DNA methylation cues in nucleosome geometry, stability and unwrapping

Cytosine methylation at the 5-carbon position is an essential DNA epigenetic mark in many eukaryotic organisms. Although countless structural and functional studies of cytosine methylation have been reported, our understanding of how it influences the nucleosome assembly, structure, and dynamics remains obscure. Here, we investigate the effects of cytosine methylation at CpG sites on nucleosome dynamics and stability. By applying long molecular dynamics simulations on several microsecond time scale, we generate extensive atomistic conformational ensembles of full nucleosomes. Our results reveal that methylation induces pronounced changes in geometry for both linker and nucleosomal DNA, leading to a more curved, under-twisted DNA, narrowing the adjacent minor grooves, and shifting the population equilibrium of sugar-phosphate backbone geometry. These DNA conformational changes are associated with a considerable enhancement of interactions between methylated DNA and the histone octamer, doubling the number of contacts at some key arginines. H2A and H3 tails play important roles in these interactions, especially for DNA methylated nucleosomes. This, in turn, prevents a spontaneous DNA unwrapping of 3-4 helical turns for the methylated nucleosome with truncated histone tails, otherwise observed in the unmethylated system on several microseconds time scale.

Chromatin network retards nucleoli coalescence

Nuclear bodies are membraneless condensates that may form via liquid-liquid phase separation. The viscoelastic chromatin network could impact their stability and may hold the key for understanding experimental observations that defy predictions of classical theories. However, quantitative studies on the role of the chromatin network in phase separation have remained challenging. Using a diploid human genome model parameterized with chromosome conformation capture (Hi-C) data, we study the thermodynamics and kinetics of nucleoli formation. Dynamical simulations predict the formation of multiple droplets for nucleolar particles that experience specific interactions with nucleolus-associated domains (NADs). Coarsening dynamics, surface tension, and coalescence kinetics of the simulated droplets are all in quantitative agreement with experimental measurements for nucleoli. Free energy calculations further support that a two-droplet state, often observed for nucleoli in somatic cells, is metastable and separated from the single-droplet state with an entropic barrier. Our study suggests that nucleoli-chromatin interactions facilitate droplets' nucleation but hinder their coarsening due to the coupled motion between droplets and the chromatin network: as droplets coalesce, the chromatin network becomes increasingly constrained. Therefore, the chromatin network supports a nucleation and arrest mechanism to stabilize the multi-droplet state for nucleoli and possibly for other nuclear bodies.

Binding of regulatory proteins to nucleosomes is modulated by dynamic histone tails

Little is known about the roles of histone tails in modulating nucleosomal DNA accessibility and its recognition by other macromolecules. Here we generate extensive atomic level conformational ensembles of histone tails in the context of the full nucleosome, totaling 65 microseconds of molecular dynamics simulations. We observe rapid conformational transitions between tail bound and unbound states, and characterize kinetic and thermodynamic properties of histone tail-DNA interactions. Different histone types exhibit distinct binding modes to specific DNA regions. Using a comprehensive set of experimental nucleosome complexes, we find that the majority of them target mutually exclusive regions with histone tails on nucleosomal/linker DNA around the super-helical locations ± 1, ± 2, and ± 7, and histone tails H3 and H4 contribute most to this process. These findings are explained within competitive binding and tail displacement models. Finally, we demonstrate the crosstalk between different histone tail post-translational modifications and mutations; those which change charge, suppress tail-DNA interactions and enhance histone tail dynamics and DNA accessibility.

The intrinsic disorder of histone tails poses challenges in their characterization. Here the authors apply extensive molecular dynamics simulations of the full nucleosome to show reversible binding to DNA with specific binding modes of different types of histone tails, where charge-altering modifications suppress tail-DNA interactions and may boost interactions between nucleosomes and nucleosome-binding proteins.

Cooperative DNA looping by PRC2 complexes

Polycomb repressive complex 2 (PRC2) is an essential protein complex that silences gene expression via post-translational modifications of chromatin. This paper combined homology modeling, atomistic and coarse-grained molecular dynamics simulations, and single-molecule force spectroscopy experiments to characterize both its full-length structure and PRC2-DNA interactions. Using free energy calculations with a newly parameterized protein-DNA force field, we studied a total of three potential PRC2 conformations and their impact on DNA binding and bending. Consistent with cryo-EM studies, we found that EZH2, a core subunit of PRC2, provides the primary interface for DNA binding, and its curved surface can induce DNA bending. Our simulations also predicted the C2 domain of the SUZ12 subunit to contact DNA. Multiple PRC2 complexes bind with DNA cooperatively via allosteric communication through the DNA, leading to a hairpin-like looped configuration. Single-molecule experiments support PRC2-mediated DNA looping and the role of AEBP2 in regulating such loop formation. The impact of AEBP2 can be partly understood from its association with the C2 domain, blocking C2 from DNA binding. Our study suggests that accessory proteins may regulate the genomic location of PRC2 by interfering with its DNA interactions.

Multiscale modeling of genome organization with maximum entropy optimization

Three-dimensional (3D) organization of the human genome plays an essential role in all DNA-templated processes, including gene transcription, gene regulation, and DNA replication. Computational modeling can be an effective way of building high-resolution genome structures and improving our understanding of these molecular processes. However, it faces significant challenges as the human genome consists of over 6 × 109 base pairs, a system size that exceeds the capacity of traditional modeling approaches. In this perspective, we review the progress that has been made in modeling the human genome. Coarse-grained models parameterized to reproduce experimental data via the maximum entropy optimization algorithm serve as effective means to study genome organization at various length scales. They have provided insight into the principles of whole-genome organization and enabled de novo predictions of chromosome structures from epigenetic modifications. Applications of these models at a near-atomistic resolution further revealed physicochemical interactions that drive the phase separation of disordered proteins and dictate chromatin stability in situ. We conclude with an outlook on the opportunities and challenges in studying chromosome dynamics.

Asymmetry in histone rotation in forced unwrapping and force quench rewrapping in a nucleosome

Single molecule pulling experiments have shown that DNA in the nucleosomes unwraps in two stages from the histone protein core (HPC). The first stage, attributed to the rupture of the outer DNA turn, occurs between 3 and 5 pNs, and is reversible. The inner DNA turn ruptures irreversibly at forces between 9 and 15 pNs (or higher) in the second stage. Molecular simulations using the Self-Organized Polymer model capture the experimental findings. The unwrapping of the outer DNA turn is independent of the pulling direction. The rupture of the DNA inner turn depends on the pulling direction and involves overcoming substantial energetic (most likely electrostatic in origin) and kinetic barriers. They arise because the mechanical force has to generate sufficient torque to rotate the HPC by 180°. On the other hand, during the rewrapping process, HPC rotation is stochastic, with force playing no role. The assembly of the outer DNA wrap upon force quench nearly coincides with the unwrapping process, confirming the reversibility of the outer turn rupture. The asymmetry in HPC rotation during unwrapping and rewrapping explains the observed hysteresis in the stretch-release cycles in experiments. We propose experiments to test the prediction that HPC rotation produces kinetic barriers in the unwrapping process.

Interplay among transacting factors around promoter in the initial phases of transcription

The initiation signals are raised around the promoter by one of the general transcription factors, triggering a sequence of events that lead to mRNA transcript formation from target genes. Both specific noncoding DNA regions and transacting, macromolecular assemblies are intricately involved and indispensable. The transition between the subsequent transcriptional stages is accompanied by stage-specific signals and structural changes in the macromolecular assemblies and facilitated by the recruitment/removal of other chromatin and transcription-associated elements. Here, we discuss the choreography of transacting factors around promoter in the establishment and effectuation of the initial phases of transcription such as NDR formation, +1 nucleosome positioning, promoter DNA opening, and RNAPII promoter escape from a structural viewpoint.

Histone tails as signaling antennas of chromatin

Histone tails, representing the N-terminal or C-terminal regions flanking the histone core, play essential roles in chromatin signaling networks. Intrinsic disorder of histone tails and their propensity for post-translational modifications allow them to serve as hubs in coordination of epigenetic processes within the nucleosomal context. Deposition of histone variants with distinct histone tail properties further enriches histone tails' repertoire in epigenetic signaling. Given the advances in experimental techniques and in silico modelling, we review the most recent data on histone tails' effects on nucleosome stability and dynamics, their function in regulating chromatin accessibility and folding. Finally, we discuss different molecular mechanisms to understand how histone tails are involved in nucleosome recognition by binding partners and formation of higher-order chromatin structures.

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

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

Inhibition of CRISPR-Cas12a DNA targeting by nucleosomes and chromatin

Genome engineering nucleases must access chromatinized DNA. Here, we investigate how AsCas12a cleaves DNA within human nucleosomes and phase-condensed nucleosome arrays. Using quantitative kinetics approaches, we show that dynamic nucleosome unwrapping regulates target accessibility to Cas12a and determines the extent to which both steps of binding-PAM recognition and R-loop formation-are inhibited by the nucleosome. Relaxing DNA wrapping within the nucleosome by reducing DNA bendability, adding histone modifications, or introducing target-proximal dCas9 enhances DNA cleavage rates over 10-fold. Unexpectedly, Cas12a readily cleaves internucleosomal linker DNA within chromatin-like, phase-separated nucleosome arrays. DNA targeting is reduced only ~5-fold due to neighboring nucleosomes and chromatin compaction. This work explains the observation that on-target cleavage within nucleosomes occurs less often than off-target cleavage within nucleosome-depleted genomic regions in cells. We conclude that nucleosome unwrapping regulates accessibility to CRISPR-Cas nucleases and propose that increasing nucleosome breathing dynamics will improve DNA targeting in eukaryotic cells.

Stability and folding pathways of tetra-nucleosome from six-dimensional free energy surface

The three-dimensional organization of chromatin is expected to play critical roles in regulating genome functions. High-resolution characterization of its structure and dynamics could improve our understanding of gene regulation mechanisms but has remained challenging. Using a near-atomistic model that preserves the chemical specificity of protein-DNA interactions at residue and base-pair resolution, we studied the stability and folding pathways of a tetra-nucleosome. Dynamical simulations performed with an advanced sampling technique uncovered multiple pathways that connect open chromatin configurations with the zigzag crystal structure. Intermediate states along the simulated folding pathways resemble chromatin configurations reported from in situ experiments. We further determined a six-dimensional free energy surface as a function of the inter-nucleosome distances via a deep learning approach. The zigzag structure can indeed be seen as the global minimum of the surface. However, it is not favored by a significant amount relative to the partially unfolded, in situ configurations. Chemical perturbations such as histone H4 tail acetylation and thermal fluctuations can further tilt the energetic balance to stabilize intermediate states. Our study provides insight into the connection between various reported chromatin configurations and has implications on the in situ relevance of the 30 nm fiber.

The three-dimensional organization of chromatin plays critical roles in regulating genome function. Here the authors apply a near atomistic model to study the structure and dynamics of the chromatin folding unit - the tetra-nucleosome - to provide insight into how chromatin folds.

Quantifying the Stability of Coupled Genetic and Epigenetic Switches With Variational Methods

The Waddington landscape provides an intuitive metaphor to view development as a ball rolling down the hill, with distinct phenotypes as basins and differentiation pathways as valleys. Since, at a molecular level, cell differentiation arises from interactions among the genes, a mathematical definition for the Waddington landscape can, in principle, be obtained by studying the gene regulatory networks. For eukaryotes, gene regulation is inextricably and intimately linked to histone modifications. However, the impact of such modifications on both landscape topography and stability of attractor states is not fully understood. In this work, we introduced a minimal kinetic model for gene regulation that combines the impact of both histone modifications and transcription factors. We further developed an approximation scheme based on variational principles to solve the corresponding master equation in a second quantized framework. By analyzing the steady-state solutions at various parameter regimes, we found that histone modification kinetics can significantly alter the behavior of a genetic network, resulting in qualitative changes in gene expression profiles. The emerging epigenetic landscape captures the delicate interplay between transcription factors and histone modifications in driving cell-fate decisions.

On the role of transcription in positioning nucleosomes

Nucleosome positioning is crucial for the genome’s function. Though the role of DNA sequence in positioning nucleosomes is well understood, a detailed mechanistic understanding on the impact of transcription remains lacking. Using numerical simulations, we investigated the dependence of nucleosome density profiles on transcription level across multiple species. We found that the low nucleosome affinity of yeast, but not mouse, promoters contributes to the formation of phased nucleosomes arrays for inactive genes. For the active genes, a heterogeneous distribution of +1 nucleosomes, caused by a tug-of-war between two types of remodeling enzymes, is essential for reproducing their density profiles. In particular, while positioning enzymes are known to remodel the +1 nucleosome and align it toward the transcription start site (TSS), spacer enzymes that use a pair of nucleosomes as their substrate can shift the nucleosome array away from the TSS. Competition between these enzymes results in two types of nucleosome density profiles with well- and ill-positioned +1 nucleosome. Finally, we showed that Pol II assisted histone exchange, if occurring at a fast speed, can abolish the impact of remodeling enzymes. By elucidating the role of individual factors, our study reconciles the seemingly conflicting results on the overall impact of transcription in positioning nucleosomes across species.

Nucleosome positioning plays a key role in the genome’s function by regulating the accessibility of protein binding sites as well as higher-order chromatin organization. Though significant progress has been made towards studying the role of DNA sequence in positioning the nucleosomes, our understanding on the impact of transcription lags behind. Our study uses kinetic simulations to explore the role of DNA sequence specificity, transcription factor binding, enzyme remodeling, and Pol II elongation in positioning nucleosomes. It suggests that the differences in nucleosome density profiles observed at various transcription levels in yeast and mouse embryonic stem cells can be understood from a tug-of-war between two types of remodeling enzymes.

Single-molecule and in silico dissection of the interaction between Polycomb repressive complex 2 and chromatin

Polycomb repressive complex 2 (PRC2) installs and spreads repressive histone methylation marks on eukaryotic chromosomes. Because of the key roles that PRC2 plays in development and disease, how this epigenetic machinery interacts with DNA and nucleosomes is of major interest. Nonetheless, the mechanism by which PRC2 engages with native-like chromatin remains incompletely understood. In this work, we employ single-molecule force spectroscopy and molecular dynamics simulations to dissect the behavior of PRC2 on polynucleosome arrays. Our results reveal an unexpectedly diverse repertoire of PRC2 binding configurations on chromatin. Besides reproducing known binding modes in which PRC2 interacts with bare DNA, mononucleosomes, and adjacent nucleosome pairs, our data also provide direct evidence that PRC2 can bridge pairs of distal nucleosomes. In particular, the "1-3" bridging mode, in which PRC2 engages two nucleosomes separated by one spacer nucleosome, is a preferred low-energy configuration. Moreover, we show that the distribution and stability of different PRC2-chromatin interaction modes are modulated by accessory subunits, oncogenic histone mutations, and the methylation state of chromatin. Overall, these findings have implications for the mechanism by which PRC2 spreads histone modifications and compacts chromatin. The experimental and computational platforms developed here provide a framework for understanding the molecular basis of epigenetic maintenance mediated by Polycomb-group proteins.

Quantifying epigenetic stability with minimum action paths

Chromatin can adopt multiple stable, heritable states with distinct histone modifications and varying levels of gene expression. Insight on the stability and maintenance of such epigenetic states can be gained by mathematical modeling of stochastic reaction networks for histone modifications. Analytical results for the kinetic networks are particularly valuable. Compared to computationally demanding numerical simulations, they often are more convenient at evaluating the robustness of conclusions with respect to model parameters. In this communication, we developed a second-quantization-based approach that can be used to analyze discrete stochastic models with a fixed, finite number of particles using a representation of the SU(2) algebra. We applied the approach to a kinetic model of chromatin states that captures the feedback between nucleosomes and the enzymes conferring histone modifications. Using a path-integral expression for the transition probability, we computed the epigenetic landscape that helps to identify the emergence of bistability and the most probable path connecting the two steady states. We anticipate the generalizability of the approach will make it useful for studying more complicated models that couple epigenetic modifications with transcription factors and chromatin structure.

Mesoscale Modeling of Nucleosome-Binding Antibody PL2-6: Mono- versus Bivalent Chromatin Complexes

Interactions of chromatin with bivalent immunoglobin nucleosome-binding antibodies and their monovalent (papain-derived) antigen-binding fragment analogs are useful probes for examining chromatin conformational states. To help interpret antibody-chromatin interactions and explore how antibodies might compete for interactions with chromatin components, we incorporate coarse-grained PL2-6 antibody modeling into our mesoscale chromatin model. We analyze interactions and fiber structures for the antibody-chromatin complexes in open and condensed chromatin, with and without H1 linker histone (LH). Despite minimal and transient interactions at physiological salt, we capture significant differences in antibody-chromatin complex configurations in open fibers, with more intense interactions between the bivalent antibody and chromatin compared to monovalent antigen-binding fragments. For these open chromatin fiber morphologies, antibody binding to histone tails is increased and compaction is greater for bivalent compared to monovalent and antibody-free systems. Differences between monovalent and bivalent binding result from antibody competition with internal chromatin fiber components (nucleosome core and linker DNA) for histone tail (H3, H4, H2A, H2B) interactions. This antibody competition for tail contacts reduces tail-core and tail-linker interactions and increases tail-antibody interactions. Such internal structural changes in open fibers resemble mechanisms of LH condensation, driven by charge screening and entropy changes. For condensed fibers at physiological salt, the three systems are much more similar overall, but some subtle tail interaction differences can be noted. Adding LH results in less-dramatic changes for all systems, except that the bivalent complex at physiological salt shows cooperative effects between LH and the antibodies in condensing chromatin fibers. Such dynamic interactions that depend on the internal structure and complex-stabilizing interactions within the chromatin fiber have implications for gene regulation and other chromatin complexes such as with LH, remodeling proteins, and small molecular chaperones that bind and modulate chromatin structure.

Theory of Active Chromatin Remodeling

Nucleosome positioning controls the accessible regions of chromatin and plays essential roles in DNA-templated processes. ATP driven remodeling enzymes are known to be crucial for its establishment in vivo, but their nonequilibrium nature has hindered the development of a unified theoretical framework for nucleosome positioning. Using a perturbation theory, we show that the effect of these enzymes can be well approximated by effective equilibrium models with rescaled temperatures and interactions. Numerical simulations support the accuracy of the theory in predicting both kinetic and steady-state quantities, including the effective temperature and the radial distribution function, in biologically relevant regimes. The energy landscape view emerging from our study provides an intuitive understanding for the impact of remodeling enzymes in either reinforcing or overwriting intrinsic signals for nucleosome positioning, and may help improve the accuracy of computational models for its prediction in silico.

The Oligomerization Landscape of Histones

In eukaryotes, DNA is packaged within nucleosomes. The DNA of each nucleosome is typically centered around an octameric histone protein core: one central tetramer plus two separate dimers. Studying the assembly mechanisms of histones is essential for understanding the dynamics of entire nucleosomes and higher-order DNA packaging. Here, we investigate canonical histone assembly and that of the centromere-specific histone variant, centromere protein A (CENP-A), using molecular dynamics simulations. We quantitatively characterize their thermodynamical and dynamical features, showing that two H3/H4 dimers form a structurally floppy, weakly bound complex, the latter exhibiting large instability around the central interface manifested via a swiveling motion of two halves. This finding is consistent with the recently observed DNA handedness flipping of the tetrasome. In contrast, the variant CENP-A encodes distinctive stability to its tetramer with a rigid but twisted interface compared to the crystal structure, implying diverse structural possibilities of the histone variant. Interestingly, the observed tetramer dynamics alter significantly and appear to reach a new balance when H2A/H2B dimers are present. Furthermore, we found that the preferred structure for the (CENP-A/H4)₂ tetramer is incongruent with the octameric structure, explaining many of the unusual dynamical behaviors of the CENP-A nucleosome. In all, these data reveal key mechanistic insights and structural details for the assembly of canonical and variant histone tetramers and octamers, providing theoretical quantifications and physical interpretations for longstanding and recent experimental observations. Based on these findings, we propose different chaperone-assisted binding and nucleosome assembly mechanisms for the canonical and CENP-A histone oligomers.