Nucleosome plasticity is a critical element of chromatin liquid-liquid phase separation and multivalent nucleosome interactions

Epigenetic histone modifications H3K36me3 and H4K5/8/12/16ac induce open polynucleosome conformations via different mechanisms

Nucleosomes are the basic compaction unit of chromatin and nucleosome structure, and their higher-order assemblies regulate genome accessibility. Many post-translational modifications alter nucleosome dynamics, nucleosome-nucleosome interactions, and ultimately chromatin structure and gene expression. Here, we investigate the role of two post-translational modifications associated with actively transcribed regions, H3K36me3 and H4K5/8/12/16ac, in the contexts of tri-nucleosome arrays that provide a tractable model system for quantitative single-molecule analysis, while enabling us to probe nucleosome-nucleosome interactions. Direct visualization by AFM imaging reveals that H3K36me3 and H4K5/8/12/16ac nucleosomes adopt significantly more open and loose conformations than unmodified nucleosomes. Similarly, magnetic tweezers force spectroscopy shows a reduction in DNA outer turn wrapping and nucleosome-nucleosome interactions for the modified nucleosomes. The results suggest that for H3K36me3 the increased breathing and outer DNA turn unwrapping seen in mononucleosomes propagates to more open conformations in nucleosome arrays. In contrast, the even more open structures of H4K5/8/12/16ac nucleosome arrays do not appear to derive from the dynamics of the constituent mononucleosomes, but are driven by reduced nucleosome-nucleosome interactions, suggesting that stacking interactions can overrule DNA breathing of individual nucleosomes. We anticipate that our methodology will be broadly applicable to reveal the influence of other post-translational modifications and to observe the activity of nucleosome remodelers.

Fundamental Aspects of Phase-Separated Biomolecular Condensates

Biomolecular condensates, formed through phase separation, are upending our understanding in much of molecular, cell, and developmental biology. There is an urgent need to elucidate the physicochemical foundations of the behaviors and properties of biomolecular condensates. Here we aim to fill this need by writing a comprehensive, critical, and accessible review on the fundamental aspects of phase-separated biomolecular condensates. We introduce the relevant theoretical background, present the theoretical basis for the computation and experimental measurement of condensate properties, and give mechanistic interpretations of condensate behaviors and properties in terms of interactions at the molecular and residue levels.

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

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

Multiscale Bayesian simulations reveal functional chromatin condensation of gene loci

Chromatin, the complex assembly of DNA and associated proteins, plays a pivotal role in orchestrating various genomic functions. To aid our understanding of the principles underlying chromatin organization, we introduce Hi-C metainference, a Bayesian approach that integrates Hi-C contact frequencies into multiscale prior models of chromatin. This approach combines both bottom-up (the physics-based prior) and top-down (the data-driven posterior) strategies to characterize the 3D organization of a target genomic locus. We first demonstrate the capability of this method to accurately reconstruct the structural ensemble and the dynamics of a system from contact information. We then apply the approach to investigate the Sox2, Pou5f1, and Nanog loci of mouse embryonic stem cells using a bottom-up chromatin model at 1 kb resolution. We observe that the studied loci are conformationally heterogeneous and organized as crumpled globules, favoring contacts between distant enhancers and promoters. Using nucleosome-resolution simulations, we then reveal how the Nanog gene is functionally organized across the multiple scales of chromatin. At the local level, we identify diverse tetranucleosome folding motifs with a characteristic distribution along the genome, predominantly open at cis-regulatory elements and compact in between. At the larger scale, we find that enhancer-promoter contacts are driven by the transient condensation of chromatin into compact domains stabilized by extensive internucleosome interactions. Overall, this work highlights the condensed, but dynamic nature of chromatin in vivo, contributing to a deeper understanding of gene structure-function relationships.

Chromatin compaction by Polycomb group proteins revisited

The chromatin compaction activity of Polycomb group proteins has traditionally been considered essential for transcriptional repression. However, there is very little information on how Polycomb group proteins compact chromatin at the molecular level and no causal link between the compactness of chromatin and transcriptional repression. Recently, a more complete picture of Polycomb-dependent chromatin architecture has started to emerge, owing to advanced methods for imaging and chromosome conformation capture. Discoveries into Polycomb-driven phase separation add another layer of complexity. Recent observations generally imply that Polycomb group proteins modulate chromatin structure at multiple scales to reduce its dynamics and segregate it from active domains. Hence, it is reasonable to hypothesise that Polycomb group proteins maintain the energetically favourable state of compacted chromatin, rather than actively compact it.

Morphological Transformations of SARS-CoV-2 Nucleocapsid Protein Biocondensates Mediated by Antimicrobial Peptides

Recently, the discovery of antimicrobial peptides (AMPs) as excellent candidates for overcoming antibiotic resistance has attracted significant attention. AMPs are short peptides active against bacteria, cancer cells, and viruses. It has been shown that the SARS-CoV-2 nucleocapsid protein (N-P) undergoes liquid-liquid phase separation in the presence of RNA, resulting in biocondensate formation. These biocondensates are crucial for viral replication as they concentrate the viral RNA with the host cell's protein machinery required for viral protein expression. Thus, N-P biocondensates are promising targets to block or slow down viral RNA transcription and consequently virion assembly. We investigated the ability of three AMPs to interfere with N-P/RNA condensates. Using microscopy techniques, supported by biophysical characterization, we found that the AMP LL-III partitions into the condensate, leading to clustering. Instead, the AMP CrACP1 partitions into the droplets without affecting their morphology but reducing their dynamics. Conversely, GKY20 leads to the formation of fibrillar structures after partitioning. It can be expected that such morphological transformation severely impairs the normal functionality of the N-P droplets and thus virion assembly. These results could pave the way for the development of a new class of AMP-based antiviral agents targeting biocondensates.

Active transcription and epigenetic reactions synergistically regulate meso-scale genomic organization

In interphase nuclei, chromatin forms dense domains of characteristic sizes, but the influence of transcription and histone modifications on domain size is not understood. We present a theoretical model exploring this relationship, considering chromatin-chromatin interactions, histone modifications, and chromatin extrusion. We predict that the size of heterochromatic domains is governed by a balance among the diffusive flux of methylated histones sustaining them and the acetylation reactions in the domains and the process of loop extrusion via supercoiling by RNAPII at their periphery, which contributes to size reduction. Super-resolution and nano-imaging of five distinct cell lines confirm the predictions indicating that the absence of transcription leads to larger heterochromatin domains. Furthermore, the model accurately reproduces the findings regarding how transcription-mediated supercoiling loss can mitigate the impacts of excessive cohesin loading. Our findings shed light on the role of transcription in genome organization, offering insights into chromatin dynamics and potential therapeutic targets.

100 Years of the Lennard-Jones Potential

It is now 100 years since Lennard-Jones published his first paper introducing the now famous potential that bears his name. It is therefore timely to reflect on the many achievements, as well as the limitations, of this potential in the theory of atomic and molecular interactions, where applications range from descriptions of intermolecular forces to molecules, clusters, and condensed matter.

Advancing biomolecular simulation through exascale HPC, AI and quantum computing

Biomolecular simulation can act as both a digital microscope and a crystal ball; offering the potential for a deeper understanding of experimental observations whilst also presenting a forward-looking avenue for the in silico design and evaluation of hitherto unsynthesized compounds. Indeed, as the intricacy of our scientific inquiries has grown, so too has the computational prowess we seek to deploy in our pursuit of answers. As we enter the Exascale era, this mini-review surveys the computational landscape from both the point of view of the development of new and ever more powerful systems, and the simulations that are run on them. Moreover, as we stand on the cusp of a transformative phase in computational biology, this article offers a contemplative glance into the future, speculating on the profound implications of artificial intelligence and quantum computing for large-scale biomolecular simulations.

Uncovering the molecular interactions underlying MBD2 and MBD3 phase separation

Chromatin organization controls DNA's accessibility to regulatory factors to influence gene expression. Heterochromatin, or transcriptionally silent chromatin enriched in methylated DNA and methylated histone tails, self-assembles through multivalent interactions with its associated proteins into a condensed, but dynamic state. Liquid-liquid phase separation (LLPS) of key heterochromatin regulators, such as heterochromatin protein 1 (HP1), plays an essential role in heterochromatin assembly and function. Methyl-CpG-binding protein 2 (MeCP2), the most studied member of the methyl-CpG-binding domain (MBD) family of proteins, has been recently shown to undergo LLPS in the absence and presence of methylated DNA. These studies provide a new mechanistic framework for understanding the role of methylated DNA and its readers in heterochromatin formation. However, the details of the molecular interactions by which other MBD family members undergo LLPS to mediate genome organization and transcriptional regulation are not fully understood. Here, we focus on two MBD proteins, MBD2 and MBD3, that have distinct but interdependent roles in gene regulation. Using an integrated computational and experimental approach, we uncover the homotypic and heterotypic interactions governing MBD2 and MBD3 phase separation and DNA's influence on this process. We show that despite sharing the highest sequence identity and structural homology among all the MBD protein family members, MBD2 and MBD3 exhibit differing residue patterns resulting in distinct phase separation mechanisms. Understanding the molecular underpinnings of MBD protein condensation offers insights into the higher-order, LLPS-mediated organization of heterochromatin.

Two-way feedback between chromatin compaction and histone modification state explains Saccharomyces cerevisiae heterochromatin bistability

Compact chromatin is closely linked with gene silencing in part by sterically masking access to promoters, inhibiting transcription factor binding and preventing polymerase from efficiently transcribing a gene. However, a broader hypothesis suggests that chromatin compaction can be both a cause and a consequence of the locus histone modification state, with a tight bidirectional interaction underpinning bistable transcriptional states. To rigorously test this hypothesis, we developed a mathematical model for the dynamics of the HMR locus in Saccharomyces cerevisiae, that incorporates activating histone modifications, silencing proteins, and a dynamic, acetylation-dependent, three-dimensional locus size. Chromatin compaction enhances silencer protein binding, which in turn feeds back to remove activating histone modifications, leading to further compaction. The bistable output of the model was in good agreement with prior quantitative data, including switching rates from expressed to silent states (and vice versa), and protein binding/histone modification levels within the locus. We then tested the model by predicting changes in switching rates as the genetic length of the locus was increased, which were then experimentally verified. Such bidirectional feedback between chromatin compaction and the histone modification state may be a widespread and important regulatory mechanism given the hallmarks of many heterochromatic regions: physical chromatin compaction and dimerizing (or multivalent) silencing proteins.

Interplay between charge distribution and DNA in shaping HP1 paralog phase separation and localization

The heterochromatin protein 1 (HP1) family is a crucial component of heterochromatin with diverse functions in gene regulation, cell cycle control, and cell differentiation. In humans, there are three paralogs, HP1α, HP1β, and HP1γ, which exhibit remarkable similarities in their domain architecture and sequence properties. Nevertheless, these paralogs display distinct behaviors in liquid-liquid phase separation (LLPS), a process linked to heterochromatin formation. Here, we employ a coarse-grained simulation framework to uncover the sequence features responsible for the observed differences in LLPS. We highlight the significance of the net charge and charge patterning along the sequence in governing paralog LLPS propensities. We also show that both highly conserved folded and less-conserved disordered domains contribute to the observed differences. Furthermore, we explore the potential co-localization of different HP1 paralogs in multicomponent assemblies and the impact of DNA on this process. Importantly, our study reveals that DNA can significantly reshape the stability of a minimal condensate formed by HP1 paralogs due to competitive interactions of HP1α with HP1β and HP1γ versus DNA. In conclusion, our work highlights the physicochemical nature of interactions that govern the distinct phase-separation behaviors of HP1 paralogs and provides a molecular framework for understanding their role in chromatin organization.

Resolving the dynamic properties of entangled linear polymers in non-equilibrium coarse grain simulation with a priori scaling factors

The molecular weight of polymers can influence the material properties, but the molecular weight at the experiment level sometimes can be a huge burden for property prediction with full-atomic simulations. The traditional bottom-up coarse grain (CG) simulation can reduce the computation cost. However, the dynamic properties predicted by the CG simulation can deviate from the full-atomic simulation result. Usually, in CG simulations, the diffusion is faster and the viscosity and modulus are much lower. The fast dynamics in CG are usually solved by a posteriori scaling on time, temperature, or potential modifications, which usually have poor transferability to other non-fitted physical properties because of a lack of fundamental physics. In this work, a priori scaling factors were calculated by the loss of degrees of freedom and implemented in the iterative Boltzmann inversion. According to the simulation results on 3 different CG levels at different temperatures and loading rates, such a priori scaling factors can help in reproducing some dynamic properties of polycaprolactone in CG simulation more accurately, such as heat capacity, Young's modulus, and viscosity, while maintaining the accuracy in the structural distribution prediction. The transferability of entropy-enthalpy compensation and a dissipative particle dynamics thermostat is also presented for comparison. The proposed method reveals the huge potential for developing customized CG thermostats and offers a simple way to rebuild multiphysics CG models for polymers with good transferability.

The molecular basis for cellular function of intrinsically disordered protein regions

Intrinsically disordered protein regions exist in a collection of dynamic interconverting conformations that lack a stable 3D structure. These regions are structurally heterogeneous, ubiquitous and found across all kingdoms of life. Despite the absence of a defined 3D structure, disordered regions are essential for cellular processes ranging from transcriptional control and cell signalling to subcellular organization. Through their conformational malleability and adaptability, disordered regions extend the repertoire of macromolecular interactions and are readily tunable by their structural and chemical context, making them ideal responders to regulatory cues. Recent work has led to major advances in understanding the link between protein sequence and conformational behaviour in disordered regions, yet the link between sequence and molecular function is less well defined. Here we consider the biochemical and biophysical foundations that underlie how and why disordered regions can engage in productive cellular functions, provide examples of emerging concepts and discuss how protein disorder contributes to intracellular information processing and regulation of cellular function.

Dynamic PRC1-CBX8 stabilizes a porous structure of chromatin condensates

The compaction of chromatin is a prevalent paradigm in gene repression. Chromatin compaction is commonly thought to repress transcription by restricting chromatin accessibility. However, the spatial organisation and dynamics of chromatin compacted by gene-repressing factors are unknown. Using cryo-electron tomography, we solved the threedimensional structure of chromatin condensed by the Polycomb Repressive Complex 1 (PRC1) in a complex with CBX8. PRC1-condensed chromatin is porous and stabilised through multivalent dynamic interactions of PRC1 with chromatin. Mechanistically, positively charged residues on the internally disordered regions (IDRs) of CBX8 mask negative charges on the DNA to stabilize the condensed state of chromatin. Within condensates, PRC1 remains dynamic while maintaining a static chromatin structure. In differentiated mouse embryonic stem cells, CBX8-bound chromatin remains accessible. These findings challenge the idea of rigidly compacted polycomb domains and instead provides a mechanistic framework for dynamic and accessible PRC1-chromatin condensates.

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.

Probe-free optical chromatin deformation and measurement of differential mechanical properties in the nucleus

The nucleus is highly organized to facilitate coordinated gene transcription. Measuring the rheological properties of the nucleus and its sub-compartments will be crucial to understand the principles underlying nuclear organization. Here, we show that strongly localized temperature gradients (approaching 1°C/µm) can lead to substantial intra-nuclear chromatin displacements (>1 µm), while nuclear area and lamina shape remain unaffected. Using particle image velocimetry (PIV), intra-nuclear displacement fields can be calculated and converted into spatio-temporally resolved maps of various strain components. Using this approach, we show that chromatin displacements are highly reversible, indicating that elastic contributions are dominant in maintaining nuclear organization on the time scale of seconds. In genetically inverted nuclei, centrally compacted heterochromatin displays high resistance to deformation, giving a rigid, solid-like appearance. Correlating spatially resolved strain maps with fluorescent reporters in conventional interphase nuclei reveals that various nuclear compartments possess distinct mechanical identities. Surprisingly, both densely and loosely packed chromatin showed high resistance to deformation, compared to medium dense chromatin. Equally, nucleoli display particularly high resistance and strong local anchoring to heterochromatin. Our results establish how localized temperature gradients can be used to drive nuclear compartments out of mechanical equilibrium to obtain spatial maps of their material responses.

Is euchromatin really open in the cell?

Genomic DNA is wrapped around a core histone octamer and forms a nucleosome. In higher eukaryotic cells, strings of nucleosomes are irregularly folded as chromatin domains that act as functional genome units. According to a typical textbook model, chromatin can be categorized into two types, euchromatin and heterochromatin, based on its degree of compaction. Euchromatin is open, while heterochromatin is closed and condensed. However, is euchromatin really open in the cell? New evidence from genomics and advanced imaging studies has revealed that euchromatin consists of condensed liquid-like domains. Condensed chromatin seems to be the default chromatin state in higher eukaryotic cells. We discuss this novel view of euchromatin in the cell and how the revealed organization is relevant to genome functions.

Through the lens of phase separation: intrinsically unstructured protein and chromatin looping

The establishment, maintenance and dynamic regulation of three-dimensional (3D) chromatin structures provide an important means for partitioning of genome into functionally distinctive domains, which helps to define specialized gene expression programs associated with developmental stages and cell types. Increasing evidence supports critical roles for intrinsically disordered regions (IDRs) harbored within transcription factors (TFs) and chromatin-modulatory proteins in inducing phase separation, a phenomenon of forming membrane-less condensates through partitioning of biomolecules. Such a process is also critically involved in the establishment of high-order chromatin structures and looping. IDR- and phase separation-driven 3D genome (re)organization often goes wrong in disease such as cancer. This review discusses about recent advances in understanding how phase separation of intrinsically disordered proteins (IDPs) modulates chromatin looping and gene expression.

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.

Catalytic and non-catalytic mechanisms of histone H4 lysine 20 methyltransferase SUV420H1

SUV420H1 di- and tri-methylates histone H4 lysine 20 (H4K20me2/H4K20me3) and plays crucial roles in DNA replication, repair, and heterochromatin formation. It is dysregulated in several cancers. Many of these processes were linked to its catalytic activity. However, deletion and inhibition of SUV420H1 have shown distinct phenotypes, suggesting that the enzyme likely has uncharacterized non-catalytic activities. Our cryoelectron microscopy (cryo-EM), biochemical, biophysical, and cellular analyses reveal how SUV420H1 recognizes its nucleosome substrates, and how histone variant H2A.Z stimulates its catalytic activity. SUV420H1 binding to nucleosomes causes a dramatic detachment of nucleosomal DNA from the histone octamer, which is a non-catalytic activity. We hypothesize that this regulates the accessibility of large macromolecular complexes to chromatin. We show that SUV420H1 can promote chromatin condensation, another non-catalytic activity that we speculate is needed for its heterochromatin functions. Together, our studies uncover and characterize the catalytic and non-catalytic mechanisms of SUV420H1, a key histone methyltransferase that plays an essential role in genomic stability.

Fluorescent protein lifetimes report densities and phases of nuclear condensates during embryonic stem-cell differentiation

Fluorescent proteins (FP) are frequently used for studying proteins inside cells. In advanced fluorescence microscopy, FPs can report on additional intracellular variables. One variable is the local density near FPs, which can be useful in studying densities within cellular bio-condensates. Here, we show that a reduction in fluorescence lifetimes of common monomeric FPs reports increased levels of local densities. We demonstrate the use of this fluorescence-based variable to report the distribution of local densities within heterochromatin protein 1α (HP1α) in mouse embryonic stem cells (ESCs), before and after early differentiation. We find that local densities within HP1α condensates in pluripotent ESCs are heterogeneous and cannot be explained by a single liquid phase. Early differentiation, however, induces a change towards a more homogeneous distribution of local densities, which can be explained as a liquid-like phase. In conclusion, we provide a fluorescence-based method to report increased local densities and apply it to distinguish between homogeneous and heterogeneous local densities within bio-condensates.

Role of Nucleosome Sliding in the Protein Target Search for Covered DNA Sites

Associations of transcription factors (TFs) with specific sites on DNA initiate major cellular processes. But DNA in eukaryotic cells is covered by nucleosomes which prevent TFs from binding. However, nucleosome structures on DNA are not static and exhibit breathing and sliding. We develop a theoretical framework to investigate the effect of nucleosome sliding on a protein target search. By analysis of a discrete-state stochastic model of nucleosome sliding, search dynamics are explicitly evaluated. It is found that for long sliding lengths the target search dynamics are faster for normal TFs that cannot enter the nucleosomal DNA. But for more realistic short sliding lengths, the so-called pioneer TFs, which can invade nucleosomal DNA, locate specific sites faster. It is also suggested that nucleosome breathing, which is a faster process, has a stronger effect on protein search dynamics than that of nucleosome sliding. Theoretical arguments to explain these observations are presented.

Heterogeneous non-canonical nucleosomes predominate in yeast cells in situ

Nuclear processes depend on the organization of chromatin, whose basic units are cylinder-shaped complexes called nucleosomes. A subset of mammalian nucleosomes in situ (inside cells) resembles the canonical structure determined in vitro 25 years ago. Nucleosome structure in situ is otherwise poorly understood. Using cryo-electron tomography (cryo-ET) and 3D classification analysis of budding yeast cells, here we find that canonical nucleosomes account for less than 10% of total nucleosomes expected in situ. In a strain in which H2A-GFP is the sole source of histone H2A, class averages that resemble canonical nucleosomes both with and without GFP densities are found ex vivo (in nuclear lysates), but not in situ. These data suggest that the budding yeast intranuclear environment favors multiple non-canonical nucleosome conformations. Using the structural observations here and the results of previous genomics and biochemical studies, we propose a model in which the average budding yeast nucleosome's DNA is partially detached in situ.

Theoretical and Data-Driven Approaches for Biomolecular Condensates

Biomolecular condensation processes are increasingly recognized as a fundamental mechanism that living cells use to organize biomolecules in time and space. These processes can lead to the formation of membraneless organelles that enable cells to perform distinct biochemical processes in controlled local environments, thereby supplying them with an additional degree of spatial control relative to that achieved by membrane-bound organelles. This fundamental importance of biomolecular condensation has motivated a quest to discover and understand the molecular mechanisms and determinants that drive and control this process. Within this molecular viewpoint, computational methods can provide a unique angle to studying biomolecular condensation processes by contributing the resolution and scale that are challenging to reach with experimental techniques alone. In this Review, we focus on three types of dry-lab approaches: theoretical methods, physics-driven simulations and data-driven machine learning methods. We review recent progress in using these tools for probing biomolecular condensation across all three fields and outline the key advantages and limitations of each of the approaches. We further discuss some of the key outstanding challenges that we foresee the community addressing next in order to develop a more complete picture of the molecular driving forces behind biomolecular condensation processes and their biological roles in health and disease.

Predicting scale-dependent chromatin polymer properties from systematic coarse-graining

Simulating chromatin is crucial for predicting genome organization and dynamics. Although coarse-grained bead-spring polymer models are commonly used to describe chromatin, the relevant bead dimensions, elastic properties, and the nature of inter-bead potentials are unknown. Using nucleosome-resolution contact probability (Micro-C) data, we systematically coarse-grain chromatin and predict quantities essential for polymer representation of chromatin. We compute size distributions of chromatin beads for different coarse-graining scales, quantify fluctuations and distributions of bond lengths between neighboring regions, and derive effective spring constant values. Unlike the prevalent notion, our findings argue that coarse-grained chromatin beads must be considered as soft particles that can overlap, and we derive an effective inter-bead soft potential and quantify an overlap parameter. We also compute angle distributions giving insights into intrinsic folding and local bendability of chromatin. While the nucleosome-linker DNA bond angle naturally emerges from our work, we show two populations of local structural states. The bead sizes, bond lengths, and bond angles show different mean behavior at Topologically Associating Domain (TAD) boundaries and TAD interiors. We integrate our findings into a coarse-grained polymer model and provide quantitative estimates of all model parameters, which can serve as a foundational basis for all future coarse-grained chromatin simulations.

G-Quadruplexes in Nuclear Biomolecular Condensates

G-quadruplexes (G4s) have long been implicated in the regulation of chromatin packaging and gene expression. These processes require or are accelerated by the separation of related proteins into liquid condensates on DNA/RNA matrices. While cytoplasmic G4s are acknowledged scaffolds of potentially pathogenic condensates, the possible contribution of G4s to phase transitions in the nucleus has only recently come to light. In this review, we summarize the growing evidence for the G4-dependent assembly of biomolecular condensates at telomeres and transcription initiation sites, as well as nucleoli, speckles, and paraspeckles. The limitations of the underlying assays and the remaining open questions are outlined. We also discuss the molecular basis for the apparent permissive role of G4s in the in vitro condensate assembly based on the interactome data. To highlight the prospects and risks of G4-targeting therapies with respect to the phase transitions, we also touch upon the reported effects of G4-stabilizing small molecules on nuclear biomolecular condensates.

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.

Evidence for low nanocompaction of heterochromatin in living embryonic stem cells

Despite advances in the identification of chromatin regulators and genome interactions, the principles of higher-order chromatin structure have remained elusive. Here, we applied FLIM-FRET microscopy to analyse, in living cells, the spatial organisation of nanometre range proximity between nucleosomes, which we called "nanocompaction." Both in naive embryonic stem cells (ESCs) and in ESC-derived epiblast-like cells (EpiLCs), we find that, contrary to expectations, constitutive heterochromatin is much less compacted than bulk chromatin. The opposite was observed in fixed cells. HP1α knockdown increased nanocompaction in living ESCs, but this was overridden by loss of HP1β, indicating the existence of a dynamic HP1-dependent low compaction state in pluripotent cells. Depletion of H4K20me2/3 abrogated nanocompaction, while increased H4K20me3 levels accompanied the nuclear reorganisation during EpiLCs induction. Finally, the knockout of the nuclear cellular-proliferation marker Ki-67 strongly reduced both interphase and mitotic heterochromatin nanocompaction in ESCs. Our data indicate that, contrary to prevailing models, heterochromatin is not highly compacted at the nanoscale but resides in a dynamic low nanocompaction state that depends on H4K20me2/3, the balance between HP1 isoforms, and Ki-67.

Condensed but liquid-like domain organization of active chromatin regions in living human cells

In eukaryotes, higher-order chromatin organization is spatiotemporally regulated as domains, for various cellular functions. However, their physical nature in living cells remains unclear (e.g., condensed domains or extended fiber loops; liquid-like or solid-like). Using novel approaches combining genomics, single-nucleosome imaging, and computational modeling, we investigated the physical organization and behavior of early DNA replicated regions in human cells, which correspond to Hi-C contact domains with active chromatin marks. Motion correlation analysis of two neighbor nucleosomes shows that nucleosomes form physically condensed domains with ~150-nm diameters, even in active chromatin regions. The mean-square displacement analysis between two neighbor nucleosomes demonstrates that nucleosomes behave like a liquid in the condensed domain on the ~150 nm/~0.5 s spatiotemporal scale, which facilitates chromatin accessibility. Beyond the micrometers/minutes scale, chromatin seems solid-like, which may contribute to maintaining genome integrity. Our study reveals the viscoelastic principle of the chromatin polymer; chromatin is locally dynamic and reactive but globally stable.

Studies of the Mechanism of Nucleosome Dynamics: A Review on Multifactorial Regulation from Computational and Experimental Cases

The nucleosome, which organizes the long coil of genomic DNA in a highly condensed, polymeric way, is thought to be the basic unit of chromosomal structure. As the most important protein–DNA complex, its structural and dynamic features have been successively revealed in recent years. However, its regulatory mechanism, which is modulated by multiple factors, still requires systemic discussion. This study summarizes the regulatory factors of the nucleosome’s dynamic features from the perspective of histone modification, DNA methylation, and the nucleosome-interacting factors (transcription factors and nucleosome-remodeling proteins and cations) and focuses on the research exploring the molecular mechanism through both computational and experimental approaches. The regulatory factors that affect the dynamic features of nucleosomes are also discussed in detail, such as unwrapping, wrapping, sliding, and stacking. Due to the complexity of the high-order topological structures of nucleosomes and the comprehensive effects of regulatory factors, the research on the functional modulation mechanism of nucleosomes has encountered great challenges. The integration of computational and experimental approaches, the construction of physical modes for nucleosomes, and the application of deep learning techniques will provide promising opportunities for further exploration.

An associative memory Hamiltonian model for DNA and nucleosomes

A model for DNA and nucleosomes is introduced with the goal of studying chromosomes from a single base level all the way to higher-order chromatin structures. This model, dubbed the Widely Editable Chromatin Model (WEChroM), reproduces the complex mechanics of the double helix including its bending persistence length and twisting persistence length, and their respective temperature dependence. The WEChroM Hamiltonian is composed of chain connectivity, steric interactions, and associative memory terms representing all remaining interactions leading to the structure, dynamics, and mechanical characteristics of the B-DNA. Several applications of this model are discussed to demonstrate its applicability. WEChroM is used to investigate the behavior of circular DNA in the presence of positive and negative supercoiling. We show that it recapitulates the formation of plectonemes and of structural defects that relax mechanical stress. The model spontaneously manifests an asymmetric behavior with respect to positive or negative supercoiling, similar to what was previously observed in experiments. Additionally, we show that the associative memory Hamiltonian is also capable of reproducing the free energy of partial DNA unwrapping from nucleosomes. WEChroM is designed to emulate the continuously variable mechanical properties of the 10nm fiber and, by virtue of its simplicity, is ready to be scaled up to molecular systems large enough to investigate the structural ensembles of genes. WEChroM is implemented in the OpenMM simulation toolkits and is freely available for public use.

Catalytic and non-catalytic mechanisms of histone H4 lysine 20 methyltransferase SUV420H1

The intricate regulation of chromatin plays a key role in controlling genome architecture and accessibility. Histone lysine methyltransferases regulate chromatin by catalyzing the methylation of specific histone residues but are also hypothesized to have equally important non-catalytic roles. SUV420H1 di- and tri-methylates histone H4 lysine 20 (H4K20me2/me3) and plays crucial roles in DNA replication, repair, and heterochromatin formation, and is dysregulated in several cancers. Many of these processes were linked to its catalytic activity. However, deletion and inhibition of SUV420H1 have shown distinct phenotypes suggesting the enzyme likely has uncharacterized non-catalytic activities. To characterize the catalytic and non-catalytic mechanisms SUV420H1 uses to modify chromatin, we determined cryo- EM structures of SUV420H1 complexes with nucleosomes containing histone H2A or its variant H2A.Z. Our structural, biochemical, biophysical, and cellular analyses reveal how both SUV420H1 recognizes its substrate and H2A.Z stimulates its activity, and show that SUV420H1 binding to nucleosomes causes a dramatic detachment of nucleosomal DNA from histone octamer. We hypothesize that this detachment increases DNA accessibility to large macromolecular complexes, a prerequisite for DNA replication and repair. We also show that SUV420H1 can promote chromatin condensates, another non-catalytic role that we speculate is needed for its heterochromatin functions. Together, our studies uncover and characterize the catalytic and non-catalytic mechanisms of SUV420H1, a key histone methyltransferase that plays an essential role in genomic stability.

Genome modeling: From chromatin fibers to genes

The intricacies of the 3D hierarchical organization of the genome have been approached by many creative modeling studies. The specific model/simulation technique combination defines and restricts the system and phenomena that can be investigated. We present the latest modeling developments and studies of the genome, involving models ranging from nucleosome systems and small polynucleosome arrays to chromatin fibers in the kb-range, chromosomes, and whole genomes, while emphasizing gene folding from first principles. Clever combinations allow the exploration of many interesting phenomena involved in gene regulation, such as nucleosome structure and dynamics, nucleosome-nucleosome stacking, polynucleosome array folding, protein regulation of chromatin architecture, mechanisms of gene folding, loop formation, compartmentalization, and structural transitions at the chromosome and genome levels. Gene-level modeling with full details on nucleosome positions, epigenetic factors, and protein binding, in particular, can in principle be scaled up to model chromosomes and cells to study fundamental biological regulation.

Conformational Change of Nucleosome Arrays prior to Phase Separation

Chromatin phase transition serves as a regulatory mechanism for eukaryotic transcription. Understanding this process requires the characterization of the nucleosome array structure in response to external stimuli prior to phase separation. However, the intrinsic flexibility and heterogeneity hinders the arrays' structure determination. Here we exploit advances in cryogenic electron tomography (cryo-ET) to determine the three-dimensional (3D) structure of each individual particle of mono-, di-, tri-, and tetranucleosome arrays. Statistical analysis reveals the ionic strength changes the angle between the DNA linker and nucleosome core particle (NCP), which regulate the overall morphology of nucleosome arrays. The finding that one-third of the arrays in the presence of H1 contain an NCP invaded by foreign DNA suggests an alternative function of H1 in constructing nucleosomal networks. The new insights into the nucleosome conformational changes prior to the intermolecular interaction stage extends our understanding of chromatin phase separation regulation.

Single-Molecule Human Nucleosome Spontaneously Ruptures under the Stress of Compressive Force: A New Perspective on Gene Stability and Epigenetic Pathways

Force manipulation on the biological entities from living cells to protein molecules has revealed many mechanical details of cell biology from resolving folding and unfolding pathways to finding molecular interaction forces. A nucleosome is the basic repeating unit of chromatin where the histone octamer is wrapped by DNA, important for gene stability and regulation. How the inner side of the DNA gets accessed by other DNA binding molecules has been a puzzle that has been intensively studied and debated, important to epigenetics, gene stability, and regulations. Here we report our observation of spontaneous ruptures of human nucleosomes under pico-Newton (pN) compressive force. The amplitude of the compressive force, a squeezing rather than pulling force, involved in our experiment is tens of pN, which can be thermally available by biological force fluctuation at room temperature and under physiological conditions. This kind of structural rupture can loosen up the DNA around the histone, which in turn makes the DNA accessible to transcription and epigenetic modifications.

The Chromatin Regulator HMGA1a Undergoes Phase Separation in the Nucleus

The protein high mobility group A1 (HMGA1) is an important regulator of chromatin organization and function. However, the mechanisms by which it exerts its biological function are not fully understood. Here, we report that the HMGA isoform, HMGA1a, nucleates into foci that display liquid-like properties in the nucleus, and that the protein readily undergoes phase separation to form liquid condensates in vitro. By bringing together machine-leaning modelling, cellular and biophysical experiments and multiscale simulations, we demonstrate that phase separation of HMGA1a is promoted by protein-DNA interactions, and has the potential to be modulated by post-transcriptional effects such as phosphorylation. We further show that the intrinsically disordered C-terminal tail of HMGA1a significantly contributes to its phase separation through electrostatic interactions via AT hooks 2 and 3. Our work sheds light on HMGA1 phase separation as an emergent biophysical factor in regulating chromatin structure.

Multiscale Modeling of Protein-RNA Condensation in and Out of Equilibrium

A vast number of intracellular membraneless bodies also known as biomolecular condensates form through a liquid-liquid phase separation (LLPS) of biomolecules. To date, phase separation has been identified as the main driving force for a membraneless organelles such as nucleoli, Cajal bodies, stress granules, and chromatin compartments. Recently, the protein-RNA condensation is receiving increased attention, because it is closely related to the biological function of cells such as transcription, translation, and RNA metabolism. Despite the multidisciplinary efforts put forth to study the biophysical properties of protein-RNA condensates, there are many fundamental unanswered questions regarding the mechanism of formation and regulation of protein-RNA condensates in eukaryotic cells. Major challenges in studying protein-RNA condensation stem from (i) the molecular heterogeneity and conformational flexibility of RNA and protein chains and (ii) the nonequilibrium nature of transcription and cellular environment. Computer simulations, bioinformatics, and mathematical models are uniquely positioned for shedding light on the microscopic nature of protein-RNA phase separation. To this end, there is an urgent need for innovative models with the right spatiotemporal resolution for confronting the experimental observables in a comprehensive and physics-based manner. In this chapter, we will summarize the currently emerging research efforts, which employ atomistic and coarse-grained molecular models and field theoretical models to understand equilibrium and nonequilibrium aspects of protein-RNA condensation.

Nucleosome breathing facilitates cooperative binding of pluripotency factors Sox2 and Oct4 to DNA

Critical lineage commitment events are staged by multiple transcription factors (TFs) binding to their cognate motifs, often positioned at nucleosome-enriched regions of chromatin. The underlying mechanism remains elusive due to difficulty in disentangling the heterogeneity in chromatin states. Using a novel coarse-grained model and molecular dynamics simulations, here we probe the association of Sox2 and Oct4 proteins that show clustered binding at the entry-exit region of a nucleosome. The model captures the conformational heterogeneity of nucleosome breathing dynamics that features repeated wrap-unwrap transitions of a DNA segment from one end of the nucleosome. During the dynamics, DNA forms bulges that diffuse stochastically and may regulate the target search dynamics of a protein by nonspecifically interacting with it. The overall search kinetics of the TF pair follows a "dissociation-compensated-association" mechanism, where Oct4 binding is facilitated by the association of Sox2. The cooperativity stems from a change in entropy caused by an alteration in the nucleosome dynamics upon TF binding. The binding pattern is consistent with a live-cell single-particle tracking experiment, suggesting the mechanism observed for clustered binding of a TF pair, which is a hallmark of cis-regulatory elements, has broader implications in understanding gene regulation in a complex chromatin environment.

Molecular dynamics simulations for the study of chromatin biology

The organization of Eukaryotic DNA into chromatin has profound implications for the processing of genetic information. In the past years, molecular dynamics (MD) simulations proved to be a powerful tool to investigate the mechanistic basis of chromatin biology. We review recent all-atom and coarse-grained MD studies revealing how the structure and dynamics of chromatin underlie its biological functions. We describe the latest method developments; the structural fluctuations of nucleosomes and the various factors affecting them; the organization of chromatin fibers, with particular emphasis on its liquid-like character; the interactions and dynamics of transcription factors on chromatin; and how chromatin organization is modulated by molecular motors acting on DNA.

The dawn of interoperating spatial models in cell biology

Spatial simulations are becoming an increasingly ubiquitous component in the cycle of discovery, experimentation, and communication across the sciences. In cell biology, many researchers share a vision of developing multiscale models that recapitulate observable behaviors spanning from atoms to cells to tissues. For this dream to become a reality, however, simulation technologies must provide a means for integration and interoperability as they advance. Already, the field has developed numerous methods that span scales of length, time, and complexity to create an extensive body of effective simulation approaches, and although these approaches rarely interoperate, they collectively cover a large spectrum of knowledge that future models may handle in a more unified manner. Here, we discuss the importance of making the data, workflows, and outputs of spatial simulations shareable and interoperable; and how democratization could encourage diverse biologists to participate more easily in developing models to advance our understanding of biological systems.

Chromatin accessibility: methods, mechanisms, and biological insights

Access to DNA is a prerequisite to the execution of essential cellular processes that include transcription, replication, chromosomal segregation, and DNA repair. How the proteins that regulate these processes function in the context of chromatin and its dynamic architectures is an intensive field of study. Over the past decade, genome-wide assays and new imaging approaches have enabled a greater understanding of how access to the genome is regulated by nucleosomes and associated proteins. Additional mechanisms that may control DNA accessibility in vivo include chromatin compaction and phase separation - processes that are beginning to be understood. Here, we review the ongoing development of accessibility measurements, we summarize the different molecular and structural mechanisms that shape the accessibility landscape, and we detail the many important biological functions that are linked to chromatin accessibility.

Contrastive Learning of Coarse-Grained Force Fields

Coarse-grained models have proven helpful for simulating complex systems over long time scales to provide molecular insights into various processes. Methodologies for systematic parametrization of the underlying energy function or force field that describes the interactions among different components of the system are of great interest for ensuring simulation accuracy. We present a new method, potential contrasting, to enable efficient learning of force fields that can accurately reproduce the conformational distribution produced with all-atom simulations. Potential contrasting generalizes the noise contrastive estimation method with umbrella sampling to better learn the complex energy landscape of molecular systems. When applied to the Trp-cage protein, we found that the technique produces force fields that thoroughly capture the thermodynamics of the folding process despite the use of only α-carbons in the coarse-grained model. We further showed that potential contrasting could be applied over large data sets that combine the conformational ensembles of many proteins to improve force field transferability. We anticipate potential contrasting as a powerful tool for building general-purpose coarse-grained force fields.

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

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

Reconstituted TAD-size chromatin fibers feature heterogeneous nucleosome clusters

Large topologically associated domains (TADs) contain irregularly spaced nucleosome clutches, and interactions between such clutches are thought to aid the compaction of these domains. Here, we reconstituted TAD-sized chromatin fibers containing hundreds of nucleosomes on native source human and lambda-phage DNA and compared their mechanical properties at the single-molecule level with shorter ‘601’ arrays with various nucleosome repeat lengths. Fluorescent imaging showed increased compaction upon saturation of the DNA with histones and increasing magnesium concentration. Nucleosome clusters and their structural fluctuations were visualized in confined nanochannels. Force spectroscopy revealed not only similar mechanical properties of the TAD-sized fibers as shorter fibers but also large rupture events, consistent with breaking the interactions between distant clutches of nucleosomes. Though the arrays of native human DNA, lambda-phage and ‘601’ DNA featured minor differences in reconstitution yield and nucleosome stability, the fibers’ global structural and mechanical properties were similar, including the interactions between nucleosome clutches. These single-molecule experiments quantify the mechanical forces that stabilize large TAD-sized chromatin domains consisting of disordered, dynamically interacting nucleosome clutches and their effect on the condensation of large chromatin domains.

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.

Molecular organization of the early stages of nucleosome phase separation visualized by cryo-electron tomography

It has been proposed that the intrinsic property of nucleosome arrays to undergo liquid-liquid phase separation (LLPS) in vitro is responsible for chromatin domain organization in vivo. However, understanding nucleosomal LLPS has been hindered by the challenge to characterize the structure of the resulting heterogeneous condensates. We used cryo-electron tomography and deep-learning-based 3D reconstruction/segmentation to determine the molecular organization of condensates at various stages of LLPS. We show that nucleosomal LLPS involves a two-step process: a spinodal decomposition process yielding irregular condensates, followed by their unfavorable conversion into more compact, spherical nuclei that grow into larger spherical aggregates through accretion of spinodal materials or by fusion with other spherical condensates. Histone H1 catalyzes more than 10-fold the spinodal-to-spherical conversion. We propose that this transition involves exposure of nucleosome hydrophobic surfaces causing modified inter-nucleosome interactions. These results suggest a physical mechanism by which chromatin may transition from interphase to metaphase structures.

Identification of Two Non-Peptidergic Small Molecule Inhibitors of CBX2 Binding to K27 Trimethylated Oligonucleosomes

The dysregulation of the PRC1/2 complex plays a key role in lineage plasticity in prostate cancer and may be required to maintain neuroendocrine phenotype. [1] CBX2, a key component of the canonical PRC1 complex, is an epigenetic reader, recognizing trimethylated lysine on histone 3 (H3K27me3) [2] and is overexpressed in metastatic neuroendocrine prostate cancer. [3,4] We implemented a screening strategy using nucleosome substrates to identify inhibitors of CBX2 binding to chromatin. Construct design and phosphorylation state of CBX2 were critical for successful implementation and execution of an HTS library screen. A rigorous screening funnel including counter and selectivity assays allowed us to quickly focus on true positive hit matter. Two distinct non-peptide-like chemotypes were identified and confirmed in orthogonal biochemical and biophysical assays demonstrating disruption of CBX2 binding to nucleosomes and direct binding to purified CBX2, respectively.