Interphase human chromosome exhibits out of equilibrium glassy dynamics

OpenNucleome for high-resolution nuclear structural and dynamical modeling

The intricate structural organization of the human nucleus is fundamental to cellular function and gene regulation. Recent advancements in experimental techniques, including high-throughput sequencing and microscopy, have provided valuable insights into nuclear organization. Computational modeling has played significant roles in interpreting experimental observations by reconstructing high-resolution structural ensembles and uncovering organization principles. However, the absence of standardized modeling tools poses challenges for furthering nuclear investigations. We present OpenNucleome-an open-source software designed for conducting GPU-accelerated molecular dynamics simulations of the human nucleus. OpenNucleome offers particle-based representations of chromosomes at a resolution of 100 KB, encompassing nuclear lamina, nucleoli, and speckles. This software furnishes highly accurate structural models of nuclear architecture, affording the means for dynamic simulations of condensate formation, fusion, and exploration of non-equilibrium effects. We applied OpenNucleome to uncover the mechanisms driving the emergence of 'fixed points' within the nucleus-signifying genomic loci robustly anchored in proximity to specific nuclear bodies for functional purposes. This anchoring remains resilient even amidst significant fluctuations in chromosome radial positions and nuclear shapes within individual cells. Our findings lend support to a nuclear zoning model that elucidates genome functionality. We anticipate OpenNucleome to serve as a valuable tool for nuclear investigations, streamlining mechanistic explorations and enhancing the interpretation of experimental observations.

Polymer physics models reveal structural folding features of single-molecule gene chromatin conformations

Here, we employ polymer physics models of chromatin to investigate the 3D folding of a 2Mb wide genomic region encompassing the human LTN1 gene, a crucial DNA locus involved in key cellular functions. Through extensive Molecular Dynamics simulations, we reconstruct in-silico the ensemble of single-molecule LTN1 3D structures, which we benchmark against recent in-situ Hi-C 2.0 data. The model-derived single molecules are then used to predict structural folding features at the single-cell level, providing testable predictions for super-resolution microscopy experiments.

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.

Cell type- and transcription-independent spatial proximity between enhancers and promoters

Cell type-specific enhancers are critically important for lineage specification. The mechanisms that determine cell-type specificity of enhancer activity, however, are not fully understood. Most current models for how enhancers function invoke physical proximity between enhancer elements and their target genes. Here, we use an imaging-based approach to examine the spatial relationship of cell type-specific enhancers and their target genes with single-cell resolution. Using high-throughput microscopy, we measure the spatial distance from target promoters to their cell type-specific active and inactive enhancers in individual pancreatic cells derived from distinct lineages. We find increased proximity of all promoter-enhancer pairs relative to non-enhancer pairs separated by similar genomic distances. Strikingly, spatial proximity between enhancers and target genes was unrelated to tissue-specific enhancer activity. Furthermore, promoter-enhancer proximity did not correlate with the expression status of target genes. Our results suggest that promoter-enhancer pairs exist in a distinctive chromatin environment but that genome folding is not a universal driver of cell-type specificity in enhancer function.

Nonequilibrium switching of segmental states can influence compaction of chromatin

Knowledge about the dynamic nature of chromatin organization is essential to understand the regulation of processes like DNA transcription and repair. The existing models of chromatin assume that protein organization and chemical states along chromatin are static and the 3D organization is purely a result of protein-mediated intra-chromatin interactions. Here we present a new hypothesis that certain nonequilibrium processes, such as switching of chemical and physical states due to nucleosome assembly/disassembly or gene repression/activation, can also simultaneously influence chromatin configurations. To understand the implications of this inherent nonequilibrium switching, we present a block copolymer model of chromatin, with switching of its segmental states between two states, mimicking active/repressed or protein unbound/bound states. We show that competition between switching timescale Tt, polymer relaxation timescale τp, and segmental relaxation timescale τs can lead to non-trivial changes in chromatin organization, leading to changes in local compaction and contact probabilities. As a function of the switching timescale, the radius of gyration of chromatin shows a non-monotonic behavior with a prominent minimum when Tt ≈ τp and a maximum when Tt ≈ τs. We find that polymers with a small segment length exhibit a more compact structure than those with larger segment lengths. We also find that the switching can lead to higher contact probability and better mixing of far-away segments. Our study also shows that the nature of the distribution of chromatin clusters varies widely as we change the switching rate.

Modelling Heterogeneous Anomalous Dynamics of Radiation-Induced Double-Strand Breaks in DNA during Non-Homologous End-Joining Pathway

The process of end-joining during nonhomologous repair of DNA double-strand breaks (DSBs) after radiation damage is considered. Experimental evidence has revealed that the dynamics of DSB ends exhibit subdiffusive motion rather than simple diffusion with rare directional movement. Traditional models often overlook the rare long-range directed motion. To address this limitation, we present a heterogeneous anomalous diffusion model consisting of subdiffusive fractional Brownian motion interchanged with short periods of long-range movement. Our model sheds light on the underlying mechanisms of heterogeneous diffusion in DSB repair and could be used to quantify the DSB dynamics on a time scale inaccessible to single particle tracking analysis. The model predicts that the long-range movement of DSB ends is responsible for the misrepair of DSBs in the form of dicentric chromosome lesions.

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.

Correlative single molecule lattice light sheet imaging reveals the dynamic relationship between nucleosomes and the local chromatin environment

In the nucleus, biological processes are driven by proteins that diffuse through and bind to a meshwork of nucleic acid polymers. To better understand this interplay, we present an imaging platform to simultaneously visualize single protein dynamics together with the local chromatin environment in live cells. Together with super-resolution imaging, new fluorescent probes, and biophysical modeling, we demonstrate that nucleosomes display differential diffusion and packing arrangements as chromatin density increases whereas the viscoelastic properties and accessibility of the interchromatin space remain constant. Perturbing nuclear functions impacts nucleosome diffusive properties in a manner that is dependent both on local chromatin density and on relative location within the nucleus. Our results support a model wherein transcription locally stabilizes nucleosomes while simultaneously allowing for the free exchange of nuclear proteins. Additionally, they reveal that nuclear heterogeneity arises from both active and passive processes and highlight the need to account for different organizational principles when modeling different chromatin environments.

OpenNucleome for high resolution nuclear structural and dynamical modeling

The intricate structural organization of the human nucleus is fundamental to cellular function and gene regulation. Recent advancements in experimental techniques, including high-throughput sequencing and microscopy, have provided valuable insights into nuclear organization. Computational modeling has played significant roles in interpreting experimental observations by reconstructing high-resolution structural ensembles and uncovering organization principles. However, the absence of standardized modeling tools poses challenges for furthering nuclear investigations. We present OpenNucleome-an open-source software designed for conducting GPU-accelerated molecular dynamics simulations of the human nucleus. OpenNucleome offers particle-based representations of chromosomes at a resolution of 100 KB, encompassing nuclear lamina, nucleoli, and speckles. This software furnishes highly accurate structural models of nuclear architecture, affording the means for dynamic simulations of condensate formation, fusion, and exploration of non-equilibrium effects. We applied OpenNucleome to uncover the mechanisms driving the emergence of "fixed points" within the nucleus-signifying genomic loci robustly anchored in proximity to specific nuclear bodies for functional purposes. This anchoring remains resilient even amidst significant fluctuations in chromosome radial positions and nuclear shapes within individual cells. Our findings lend support to a nuclear zoning model that elucidates genome functionality. We anticipate OpenNucleome to serve as a valuable tool for nuclear investigations, streamlining mechanistic explorations and enhancing the interpretation of experimental observations.

Structural basis for the preservation of a subset of topologically associating domains in interphase chromosomes upon cohesin depletion

Compartment formation in interphase chromosomes is a result of spatial segregation between euchromatin and heterochromatin on a few megabase pairs (Mbp) scale. On the sub-Mbp scales, topologically associating domains (TADs) appear as interacting domains along the diagonal in the ensemble averaged Hi-C contact map. Hi-C experiments showed that most of the TADs vanish upon deleting cohesin, while the compartment structure is maintained, and perhaps even enhanced. However, closer inspection of the data reveals that a non-negligible fraction of TADs is preserved (P-TADs) after cohesin loss. Imaging experiments show that, at the single-cell level, TAD-like structures are present even without cohesin. To provide a structural basis for these findings, we first used polymer simulations to show that certain TADs with epigenetic switches across their boundaries survive after depletion of loops. More importantly, the three-dimensional structures show that many of the P-TADs have sharp physical boundaries. Informed by the simulations, we analyzed the Hi-C maps (with and without cohesin) in mouse liver and human colorectal carcinoma cell lines, which affirmed that epigenetic switches and physical boundaries (calculated using the predicted 3D structures using the data-driven HIPPS method that uses Hi-C as the input) explain the origin of the P-TADs. Single-cell structures display TAD-like features in the absence of cohesin that are remarkably similar to the findings in imaging experiments. Some P-TADs, with physical boundaries, are relevant to the retention of enhancer-promoter/promoter-promoter interactions. Overall, our study shows that preservation of a subset of TADs upon removing cohesin is a robust phenomenon that is valid across multiple cell lines.

Transcription-induced active forces suppress chromatin motion

The organization of interphase chromosomes in a number of species is starting to emerge thanks to advances in a variety of experimental techniques. However, much less is known about the dynamics, especially in the functional states of chromatin. Some experiments have shown that the motility of individual loci in human interphase chromosome decreases during transcription and increases upon inhibiting transcription. This is a counterintuitive finding because it is thought that the active mechanical force (F) on the order of ten piconewtons, generated by RNA polymerase II (RNAPII) that is presumably transmitted to the gene-rich region of the chromatin, would render it more open, thus enhancing the mobility. We developed a minimal active copolymer model for interphase chromosomes to investigate how F affects the dynamical properties of chromatin. The movements of the loci in the gene-rich region are suppressed in an intermediate range of F and are enhanced at small F values, which has also been observed in experiments. In the intermediate F, the bond length between consecutive loci increases, becoming commensurate with the distance at the minimum of the attractive interaction between nonbonded loci. This results in a transient disorder-to-order transition, leading to a decreased mobility during transcription. Strikingly, the F-dependent change in the locus dynamics preserves the organization of the chromosome at [Formula: see text]. Transient ordering of the loci, which is not found in the polymers with random epigenetic profiles, in the gene-rich region might be a plausible mechanism for nucleating a dynamic network involving transcription factors, RNAPII, and chromatin.

Steric repulsion introduced by loop constraints modulates the microphase separation of chromatins

Within the confines of a densely populated cell nucleus, chromatin undergoes intricate folding, forming loops, domains, and compartments under the governance of topological constraints and phase separation. This coordinated process inevitably introduces interference between different folding strategies. In this study, we model interphase chromatins as block copolymers with hetero-hierarchical loops within a confined system. Employing dissipative particle dynamics simulations and scaling analysis, we aim to explain how the structure and distribution of loop domains modulate the microphase separation of chromatins. Our results highlight the correlation between the microphase separation of the copolymer and the length, heterogeneity, and hierarchically nested levels of the loop domains. This correlation arises from steric repulsion intrinsic to loop domains. The steric repulsion induces variations in chain stiffness (including local orientation correlations and the persistence length), thereby influencing the degree of phase separation. Through simulations of block copolymers with distinct groups of hetero-hierarchical loop anchors, we successfully reproduce changes in phase separation across diverse cell lines, under fixed interaction parameters. These findings, in qualitative alignment with Hi-C data, suggest that the variations of loop constraints alone possess the capacity to regulate higher-order structures and the gene expressions of interphase chromatins.

Chromatin Remodeling Due to Transient-Link-and-Pass Activity Enhances Subnuclear Dynamics

Spatiotemporal coordination of chromatin and subnuclear compartments is crucial for cells. Numerous enzymes act inside nucleus-some of those transiently link and pass two chromatin segments. Here, we study how such an active perturbation affects fluctuating dynamics of an inclusion in the chromatic medium. Using numerical simulations and a versatile effective model, we categorize inclusion dynamics into three distinct modes. The transient-link-and-pass activity speeds up inclusion dynamics by affecting a slow mode related to chromatin remodeling, viz., size and shape of the chromatin meshes.

Computational methods for analysing multiscale 3D genome organization

Recent progress in whole-genome mapping and imaging technologies has enabled the characterization of the spatial organization and folding of the genome in the nucleus. In parallel, advanced computational methods have been developed to leverage these mapping data to reveal multiscale three-dimensional (3D) genome features and to provide a more complete view of genome structure and its connections to genome functions such as transcription. Here, we discuss how recently developed computational tools, including machine-learning-based methods and integrative structure-modelling frameworks, have led to a systematic, multiscale delineation of the connections among different scales of 3D genome organization, genomic and epigenomic features, functional nuclear components and genome function. However, approaches that more comprehensively integrate a wide variety of genomic and imaging datasets are still needed to uncover the functional role of 3D genome structure in defining cellular phenotypes in health and disease.

Modeling the compaction of bacterial chromosomes by biomolecular crowding and the cross-linking protein H-NS

Cells orchestrate the action of various molecules toward organizing their chromosomes. Using a coarse-grained computational model, we study the compaction of bacterial chromosomes by the cross-linking protein H-NS and cellular crowders. In this work, H-NS, modeled as a mobile "binder," can bind to a chromosome-like polymer with a characteristic binding energy. The simulation results reported here clarify the relative role of biomolecular crowding and H-NS in condensing a bacterial chromosome in a quantitative manner. In particular, they shed light on the nature and degree of crowder and H-NS synergetics: while the presence of crowders enhances H-NS binding to a chromosome-like polymer, the presence of H-NS makes crowding effects more efficient, suggesting two-way synergetics in chain compaction. Also, the results show how crowding effects promote clustering of bound H-NS. For a sufficiently large concentration of H-NS, the cluster size increases with the volume fraction of crowders.

Diffusion controls local versus dispersed inheritance of histones during replication and shapes epigenomic architecture

The dynamics of inheritance of histones and their associated modifications across cell divisions can have major consequences on maintenance of the cellular epigenomic state. Recent experiments contradict the long-held notion that histone inheritance during replication is always local, suggesting that active and repressed regions of the genome exhibit fundamentally different histone dynamics independent of transcription-coupled turnover. Here we develop a stochastic model of histone dynamics at the replication fork and demonstrate that differential diffusivity of histones in active versus repressed chromatin is sufficient to quantitatively explain these recent experiments. Further, we use the model to predict patterns in histone mark similarity between pairs of genomic loci that should be developed as a result of diffusion, but cannot originate from either PRC2 mediated mark spreading or transcriptional processes. Interestingly, using a combination of CHIP-seq, replication timing and Hi-C datasets we demonstrate that all the computationally predicted patterns are consistently observed for both active and repressive histone marks in two different cell lines. While direct evidence for histone diffusion remains controversial, our results suggest that dislodged histones in euchromatin and facultative heterochromatin may exhibit some level of diffusion within "Diffusion-Accessible-Domains" (DADs), leading to redistribution of epigenetic marks within and across chromosomes. Preservation of the epigenomic state across cell divisions therefore might be achieved not by passing on strict positional information of histone marks, but by maintaining the marks in somewhat larger DADs of the genome.

Correlative single molecule lattice light sheet imaging reveals the dynamic relationship between nucleosomes and the local chromatin environment

In the nucleus, biological processes are driven by proteins that diffuse through and bind to a meshwork of nucleic acid polymers. To better understand this interplay, we developed an imaging platform to simultaneously visualize single protein dynamics together with the local chromatin environment in live cells. Together with super-resolution imaging, new fluorescent probes, and biophysical modeling, we demonstrated that nucleosomes display differential diffusion and packing arrangements as chromatin density increases whereas the viscoelastic properties and accessibility of the interchromatin space remain constant. Perturbing nuclear functions impacted nucleosome diffusive properties in a manner that was dependent on local chromatin density and supportive of a model wherein transcription locally stabilizes nucleosomes while simultaneously allowing for the free exchange of nuclear proteins. Our results reveal that nuclear heterogeneity arises from both active and passive process and highlights the need to account for different organizational principals when modeling different chromatin environments.

Transcription Factor Dynamics: One Molecule at a Time

Cells must tightly regulate their gene expression programs and yet rapidly respond to acute biochemical and biophysical cues within their environment. This information is transmitted to the nucleus through various signaling cascades, culminating in the activation or repression of target genes. Transcription factors (TFs) are key mediators of these signals, binding to specific regulatory elements within chromatin. While live-cell imaging has conclusively proven that TF-chromatin interactions are highly dynamic, how such transient interactions can have long-term impacts on developmental trajectories and disease progression is still largely unclear. In this review, we summarize our current understanding of the dynamic nature of TF functions, starting with a historical overview of early live-cell experiments. We highlight key factors that govern TF dynamics and how TF dynamics, in turn, affect downstream transcriptional bursting. Finally, we conclude with open challenges and emerging technologies that will further our understanding of transcriptional regulation.

Leveraging polymer modeling to reconstruct chromatin connectivity from live images

Chromosomal dynamics plays a central role in a number of critical biological processes, such as transcriptional regulation, genetic recombination, and DNA replication. However, visualization of chromatin is generally limited to live imaging of a few fluorescently labeled chromosomal loci or high-resolution reconstruction of multiple loci from a single time frame. To aid in mapping the underlying chromosomal structure based on parsimonious experimental measurements, we present an exact analytical expression for the evolution of the polymer configuration based on a flexible-polymer model, and we propose an algorithm that tracks the polymer configuration from live images of chromatin marked with several fluorescent marks. Our theory identifies the resolution of microscopy needed to achieve high-accuracy tracking for a given spacing of markers, establishing the statistical confidence in the assignment of genome identity to the visualized marks. We then leverage experimental data of locus-tracking measurements to demonstrate the validity of our modeling approach and to establish a basis for the design of experiments with a desired resolution. Altogether, this work provides a computational approach founded on polymer physics that vastly improves the interpretation of in vivo measurements of biopolymer dynamics.

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.

Model chromatin flows: numerical analysis of linear and nonlinear hydrodynamics inside a sphere

We solve a hydrodynamic model of active chromatin dynamics, within a confined geometry simulating the cell nucleus. Using both analytical and numerical methods, we describe the behavior of the chromatin polymer driven by the activity of motors having polar symmetry, both in the linear response regime as well as in the long-term, fully nonlinear regime of the flows. The introduction of a boundary induces a particular geometry in the flows of chromatin, which we describe using vector spherical harmonics, a tool which greatly simplifies both our analytical and numerical approaches. We find that the long-term behavior of this model in confinement is dominated by steady, transverse flows of chromatin which circulate around the spherical domain. These circulating flows are found to be robust to perturbations, and their characteristic size is set by the size of the domain. This gives us further insight into active chromatin dynamics in the cell nucleus, and provides a foundation for development of further, more complex models of active chromatin dynamics.

Multiscale modelling of chromatin 4D organization in SARS-CoV-2 infected cells

SARS-CoV-2 is able to re-structure chromatin organization and alters the epigenomic landscape of the host genome, though the mechanisms that produce such changes are still poorly understood. Here, we investigate with polymer physics chromatin re-organization of the host genome, in space and time upon SARS-CoV-2 viral infection. We show that re-structuring of A/B compartments is well explained by a re-modulation of intra-compartment homotypic affinities, which leads to the weakening of A-A interactions and enhances A-B mixing. At TAD level, re-arrangements are physically described by a general reduction of the loop extrusion activity coupled with an alteration of chromatin phase-separation properties, resulting in more intermingling between different TADs and spread in space of TADs themselves. In addition, the architecture of loci relevant to the antiviral interferon (IFN) response, such as DDX58 or IFIT, results more variable within the 3D single-molecule population of the infected model, suggesting that viral infection leads to a loss of chromatin structural specificity. Analysis of time trajectories of pairwise gene-enhancer and higher-order contacts reveals that such variability derives from a more fluctuating dynamics in infected case, suggesting that SARS-CoV-2 alters gene regulation by impacting the stability of the contact network in time. Overall, our study provides the first polymer-physics based 4D reconstruction of SARS-CoV-2 infected genome with mechanistic insights on the consequent gene mis-regulation.

Activity-Driven Phase Transition Causes Coherent Flows of Chromatin

We discover a new type of nonequilibrium phase transition in a model of chromatin dynamics, which accounts for the coherent motions that have been observed in experiment. The coherent motion is due to the long-range cooperation of molecular motors tethered to chromatin. Cooperation occurs if each motor acts simultaneously on the polymer and the surrounding solvent, exerting on them equal and opposite forces. This drives the flow of solvent past the polymer, which in turn affects the orientation of nearby motors and, if the drive is strong enough, an active polar ("ferromagnetic") phase of motors can spontaneously form. Depending on boundary conditions, either transverse flows or sustained longitudinal oscillations and waves are possible. Predicted length scales are consistent with experiments. We now have in hand a coarse-grained description of chromatin dynamics which reproduces the directed coherent flows of chromatin seen in experiments. This field-theoretic description can be analytically coupled to other features of the nuclear environment such as fluctuating or porous boundaries, local heterogeneities in the distribution of chromatin or its activity, leading to insights on the effects of activity on the cell nucleus and its contents.

Brownian dynamics simulations of mesoscale chromatin fibers

The relationship between chromatin architecture and function defines a central problem in biology and medicine. Many computational chromatin models with atomic, coarse-grained, mesoscale, and polymer resolution have been used to shed light onto the mechanisms that dictate genome folding and regulation of gene expression. The associated simulation techniques range from Monte Carlo to molecular, Brownian, and Langevin dynamics. Here, we present an efficient Compute Unified Device Architecture (CUDA) implementation of Brownian dynamics (BD) to simulate chromatin fibers at the nucleosome resolution with our chromatin mesoscale model. With the CUDA implementation for computer architectures with graphic processing units (GPUs), we significantly accelerate compute-intensive hydrodynamic tensor calculations in the BD simulations by massive parallelization, boosting the performance a hundred-fold compared with central processing unit calculations. We validate our BD simulation approach by reproducing experimental trends on fiber diffusion and structure as a function of salt, linker histone binding, and histone-tail composition, as well as Monte Carlo equilibrium sampling results. Our approach proves to be physically accurate with performance that makes feasible the study of chromatin fibers in the range of kb or hundreds of nucleosomes (small gene). Such simulations are essential to advance the study of biological processes such as gene regulation and aberrant genome-structure related diseases.

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.

Strong interactions between highly dynamic lamina-associated domains and the nuclear envelope stabilize the 3D architecture of Drosophila interphase chromatin

BACKGROUND

Interactions among topologically associating domains (TADs), and between the nuclear envelope (NE) and lamina-associated domains (LADs) are expected to shape various aspects of three-dimensional (3D) chromatin structure and dynamics; however, relevant genome-wide experiments that may provide statistically significant conclusions remain difficult.

RESULTS

We have developed a coarse-grained dynamical model of D. melanogaster nuclei at TAD resolution that explicitly accounts for four distinct epigenetic classes of TADs and LAD-NE interactions. The model is parameterized to reproduce the experimental Hi-C map of the wild type (WT) nuclei; it describes time evolution of the chromatin over the G1 phase of the interphase. The simulations include an ensemble of nuclei, corresponding to the experimentally observed set of several possible mutual arrangements of chromosomal arms. The model is validated against multiple structural features of chromatin from several different experiments not used in model development. Predicted positioning of all LADs at the NE is highly dynamic-the same LAD can attach, detach and move far away from the NE multiple times during interphase. The probabilities of LADs to be in contact with the NE vary by an order of magnitude, despite all having the same affinity to the NE in the model. These probabilities are mostly determined by a highly variable local linear density of LADs along the genome, which also has the same strong effect on the predicted positioning of individual TADs -- higher probability of a TAD to be near NE is largely determined by a higher linear density of LADs surrounding this TAD. The distribution of LADs along the chromosome chains plays a notable role in maintaining a non-random average global structure of chromatin. Relatively high affinity of LADs to the NE in the WT nuclei substantially reduces sensitivity of the global radial chromatin distribution to variations in the strength of TAD-TAD interactions compared to the lamin depleted nuclei, where a small (0.5 kT) increase of cross-type TAD-TAD interactions doubles the chromatin density in the central nucleus region.

CONCLUSIONS

A dynamical model of the entire fruit fly genome makes multiple genome-wide predictions of biological interest. The distribution of LADs along the chromatin chains affects their probabilities to be in contact with the NE and radial positioning of highly mobile TADs, playing a notable role in creating a non-random average global structure of the chromatin. We conjecture that an important role of attractive LAD-NE interactions is to stabilize global chromatin structure against inevitable cell-to-cell variations in TAD-TAD interactions.

Disruption of polyhomeotic polymerization decreases nucleosome occupancy and alters genome accessibility

Chromatin attains its three-dimensional (3D) conformation by establishing contacts between different noncontiguous regions. Sterile Alpha Motif (SAM)-mediated polymerization of the polyhomeotic (PH) protein regulates subnuclear clustering of Polycomb Repressive Complex 1 (PRC1) and chromatin topology. The mutations that perturb the ability of the PH to polymerize, disrupt long-range chromatin contacts, alter Hox gene expression, and lead to developmental defects. To understand the underlying mechanism, we combined the experiments and theory to investigate the effect of this SAM domain mutation on nucleosome occupancy and accessibility on a genome wide scale. Our data show that disruption of PH polymerization because of SAM domain mutation decreases nucleosome occupancy and alters accessibility. Polymer simulations investigating the interplay between distant chromatin contacts and nucleosome occupancy, both of which are regulated by PH polymerization, suggest that nucleosome density increases when contacts between different regions of chromatin are established. Taken together, it appears that SAM domain-mediated PH polymerization biomechanically regulates the organization of chromatin at multiple scales from nucleosomes to chromosomes and we suggest that higher order organization can have a top-down causation effect on nucleosome occupancy.

Structural changes in chromosomes driven by multiple condensin motors during mitosis

We create a computational framework that utilizes loop extrusion (LE) by multiple condensin I/II motors to predict changes in chromosome organization during mitosis. The theory accurately reproduces the experimental contact probability profiles for the mitotic chromosomes in HeLa and DT40 cells. The LE rate is smaller at the start of mitosis and increases as the cells approach metaphase. Condensin II-mediated mean loop size is about six times larger than loops because of condensin I. The loops, which overlap each other, are stapled to a central dynamically changing helical scaffold formed by the motors during the LE process. A polymer physics-based data-driven method that uses the Hi-C contact map as the only input shows that the helix is characterized as random helix perversions (RHPs) in which the handedness changes randomly along the scaffold. The theoretical predictions, which are testable using imaging experiments, do not contain any parameters.

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.

Phase transitions in heterochromatin organization

Heterochromatin formation has been proposed to involve phase transitions on the level of the three-dimensional folding of heterochromatin regions and the liquid-liquid demixing of heterochromatin proteins. Here, I outline the hallmarks of such transitions and the current challenges to detect them in living cells. I further discuss the abundance and properties of prominent heterochromatin proteins and relate them to their potential role in driving phase transitions. Recent data from mouse fibroblasts indicate that pericentric heterochromatin is organized via a reordering transition on the level of heterochromatin regions that does not necessarily involve liquid-liquid demixing of heterochromatin proteins. Evaluating key hallmarks of the different candidate phase transition mechanisms across cell types and species will be needed to complete the current picture.

Compartmentalization with nuclear landmarks yields random, yet precise, genome organization

The 3D organization of eukaryotic genomes plays an important role in genome function. While significant progress has been made in deciphering the folding mechanisms of individual chromosomes, the principles of the dynamic large-scale spatial arrangement of all chromosomes inside the nucleus are poorly understood. We use polymer simulations to model the diploid human genome compartmentalization relative to nuclear bodies such as nuclear lamina, nucleoli, and speckles. We show that a self-organization process based on a cophase separation between chromosomes and nuclear bodies can capture various features of genome organization, including the formation of chromosome territories, phase separation of A/B compartments, and the liquid property of nuclear bodies. The simulated 3D structures quantitatively reproduce both sequencing-based genomic mapping and imaging assays that probe chromatin interaction with nuclear bodies. Importantly, our model captures the heterogeneous distribution of chromosome positioning across cells while simultaneously producing well-defined distances between active chromatin and nuclear speckles. Such heterogeneity and preciseness of genome organization can coexist due to the nonspecificity of phase separation and the slow chromosome dynamics. Together, our work reveals that the cophase separation provides a robust mechanism for us to produce functionally important 3D contacts without requiring thermodynamic equilibration that can be difficult to achieve.

4D epigenomics: deciphering the coupling between genome folding and epigenomic regulation with biophysical modeling

Recent experimental observations suggest a strong coupling between the 3D nuclear chromosome organization and epigenomics. However, the mechanistic and functional bases of such interplay remain elusive. In this review, we describe how biophysical modeling has been instrumental in characterizing how genome folding may impact the formation of epigenomic domains and, conversely, how epigenomic marks may affect chromosome conformation. Finally, we discuss how this mutual feedback loop between chromatin organization and epigenome regulation, via the formation of physicochemical nanoreactors, may represent a key functional role of 3D compartmentalization in the assembly and maintenance of stable - but yet plastic - epigenomic landscapes.

A maximum-entropy model to predict 3D structural ensembles of chromatin from pairwise distances with applications to interphase chromosomes and structural variants

The principles that govern the organization of genomes, which are needed for an understanding of how chromosomes are packaged and function in eukaryotic cells, could be deciphered if the three-dimensional (3D) structures are known. Recently, single-cell imaging techniques have been developed to determine the 3D coordinates of genomic loci in vivo. Here, we introduce a computational method (Distance Matrix to Ensemble of Structures, DIMES), based on the maximum entropy principle, with experimental pairwise distances between loci as constraints, to generate a unique ensemble of 3D chromatin structures. Using the ensemble of structures, we quantitatively account for the distribution of pairwise distances, three-body co-localization, and higher-order interactions. The DIMES method can be applied to both small and chromosome-scale imaging data to quantify the extent of heterogeneity and fluctuations in the shapes across various length scales. We develop a perturbation method in conjunction with DIMES to predict the changes in 3D structures from structural variations. Our method also reveals quantitative differences between the 3D structures inferred from Hi-C and those measured in imaging experiments. Finally, the physical interpretation of the parameters extracted from DIMES provides insights into the origin of phase separation between euchromatin and heterochromatin domains.

Unveiling the Machinery behind Chromosome Folding by Polymer Physics Modeling

Understanding the mechanisms underlying the complex 3D architecture of mammalian genomes poses, at a more fundamental level, the problem of how two or multiple genomic sites can establish physical contacts in the nucleus of the cells. Beyond stochastic and fleeting encounters related to the polymeric nature of chromatin, experiments have revealed specific, privileged patterns of interactions that suggest the existence of basic organizing principles of folding. In this review, we focus on two major and recently proposed physical processes of chromatin organization: loop-extrusion and polymer phase-separation, both supported by increasing experimental evidence. We discuss their implementation into polymer physics models, which we test against available single-cell super-resolution imaging data, showing that both mechanisms can cooperate to shape chromatin structure at the single-molecule level. Next, by exploiting the comprehension of the underlying molecular mechanisms, we illustrate how such polymer models can be used as powerful tools to make predictions in silico that can complement experiments in understanding genome folding. To this aim, we focus on recent key applications, such as the prediction of chromatin structure rearrangements upon disease-associated mutations and the identification of the putative chromatin organizing factors that orchestrate the specificity of DNA regulatory contacts genome-wide.

Interphase chromosomes of the Aedes aegypti mosquito are liquid crystalline and can sense mechanical cues

We use data-driven physical simulations to study the three-dimensional architecture of the Aedes aegypti genome. Hi-C maps exhibit both a broad diagonal and compartmentalization with telomeres and centromeres clustering together. Physical modeling reveals that these observations correspond to an ensemble of 3D chromosomal structures that are folded over and partially condensed. Clustering of the centromeres and telomeres near the nuclear lamina appears to be a necessary condition for the formation of the observed structures. Further analysis of the mechanical properties of the genome reveals that the chromosomes of Aedes aegypti, by virtue of their atypical structural organization, are highly sensitive to the deformation of the nuclei. This last finding provides a possible physical mechanism linking mechanical cues to gene regulation.

The spatial organization of transcriptional control

In animals, the sequences for controlling gene expression do not concentrate just at the transcription start site of genes, but are frequently thousands to millions of base pairs distal to it. The interaction of these sequences with one another and their transcription start sites is regulated by factors that shape the three-dimensional (3D) organization of the genome within the nucleus. Over the past decade, indirect tools exploiting high-throughput DNA sequencing have helped to map this 3D organization, have identified multiple key regulators of its structure and, in the process, have substantially reshaped our view of how 3D genome architecture regulates transcription. Now, new tools for high-throughput super-resolution imaging of chromatin have directly visualized the 3D chromatin organization, settling some debates left unresolved by earlier indirect methods, challenging some earlier models of regulatory specificity and creating hypotheses about the role of chromatin structure in transcriptional regulation.

Nuclear position modulates long-range chromatin interactions

The human genome is non-randomly organized within the cell nucleus. Spatial mapping of genome folding by biochemical methods and imaging has revealed extensive variation in locus interaction frequencies between cells in a population and between homologs within an individual cell. Commonly used mapping approaches typically examine either the relative position of genomic sites to each other or the position of individual loci relative to nuclear landmarks. Whether the frequency of specific chromatin-chromatin interactions is affected by where in the nuclear space a locus is located is unknown. Here, we have simultaneously mapped at the single cell level the interaction frequencies and radial position of more than a hundred locus pairs using high-throughput imaging to ask whether the location within the nucleus affects interaction frequency. We find strong enrichment of many interactions at specific radial positions. Position-dependency of interactions was cell-type specific, correlated with local chromatin type, and cell-type-specific enriched associations were marked by increased variability, sometimes without a significant decrease in mean spatial distance. These observations demonstrate that the folding of the chromatin fiber, which brings genomically distant loci into proximity, and the position of that chromatin fiber relative to nuclear landmarks, are closely linked.

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.

DNA damage reduces heterogeneity and coherence of chromatin motions

Chromatin motions depend on and may regulate genome functions, in particular the DNA damage response. In yeast, DNA double-strand breaks (DSBs) globally increase chromatin diffusion, whereas in higher eukaryotes the impact of DSBs on chromatin dynamics is more nuanced. We mapped the motions of chromatin microdomains in mammalian cells using diffractive optics and photoactivatable chromatin probes and found a high level of spatial heterogeneity. DNA damage reduces heterogeneity and imposes spatially defined shifts in motions: Distal to DNA breaks, chromatin motions are globally reduced, whereas chromatin retains higher mobility at break sites. These effects are driven by context-dependent changes in chromatin compaction. Photoactivated lattices of chromatin microdomains are ideal to quantify microscale coupling of chromatin motion. We measured correlation distances up to 2 µm in the cell nucleus, spanning chromosome territories, and speculate that this correlation distance between chromatin microdomains corresponds to the physical separation of A and B compartments identified in chromosome conformation capture experiments. After DNA damage, chromatin motions become less correlated, a phenomenon driven by phase separation at DSBs. Our data indicate tight spatial control of chromatin motions after genomic insults, which may facilitate repair at the break sites and prevent deleterious contacts of DSBs, thereby reducing the risk of genomic rearrangements.

Heterogeneous interactions and polymer entropy decide organization and dynamics of chromatin domains

Chromatin is known to be organized into multiple domains of varying sizes and compaction. While these domains are often imagined as static structures, they are highly dynamic and show cell-to-cell variability. Since processes such as gene regulation and DNA replication occur in the context of these domains, it is important to understand their organization, fluctuation, and dynamics. To simulate chromatin domains, one requires knowledge of interaction strengths among chromatin segments. Here, we derive interaction-strength parameters from experimentally known contact maps and use them to predict chromatin organization and dynamics. Taking two domains on the human chromosome as examples, we investigate its three-dimensional organization, size/shape fluctuations, and dynamics of different segments within a domain, accounting for hydrodynamic effects. Considering different cell types, we quantify changes in interaction strengths and chromatin shape fluctuations in different epigenetic states. Perturbing the interaction strengths systematically, we further investigate how epigenetic-like changes can alter the spatio-temporal nature of the domains. Our results show that heterogeneous weak interactions are crucial in determining the organization of the domains. Computing effective stiffness and relaxation times, we investigate how perturbations in interactions affect the solid- and liquid-like nature of chromatin domains. Quantifying dynamics of chromatin segments within a domain, we show how the competition between polymer entropy and interaction energy influence the timescales of loop formation and maintenance of stable loops.

A simulation model of heterochromatin formation at submolecular detail

Heterochromatin is a physical state of the chromatin fiber that maintains gene repression during cell development. Although evidence exists on molecular mechanisms involved in heterochromatin formation, a detailed structural mechanism of heterochromatin formation needs a better understanding. We made use of a simple Monte Carlo simulation model with explicit representation of key molecular events to observe molecular self-organization leading to heterochromatin formation. Our simulations provide a structural interpretation of several important traits of the heterochromatinization process. In particular, this study provides a depiction of how small amounts of HP1 are able to induce a highly condensed chromatin state through HP1 dimerization and bridging of sequence-remote nucleosomes. It also elucidates structural roots of a yet poorly understood phenomenon of a nondeterministic nature of heterochromatin formation and subsequent gene repression. Experimental chromatin in vivo assay provides an unbiased estimate of time scale of repressive response to a heterochromatin-triggering event.

Loop-extrusion and polymer phase-separation can co-exist at the single-molecule level to shape chromatin folding

Loop-extrusion and phase-separation have been proposed as mechanisms that shape chromosome spatial organization. It is unclear, however, how they perform relative to each other in explaining chromatin architecture data and whether they compete or co-exist at the single-molecule level. Here, we compare models of polymer physics based on loop-extrusion and phase-separation, as well as models where both mechanisms act simultaneously in a single molecule, against multiplexed FISH data available in human loci in IMR90 and HCT116 cells. We find that the different models recapitulate bulk Hi-C and average multiplexed microscopy data. Single-molecule chromatin conformations are also well captured, especially by phase-separation based models that better reflect the experimentally reported segregation in globules of the considered genomic loci and their cell-to-cell structural variability. Such a variability is consistent with two main concurrent causes: single-cell epigenetic heterogeneity and an intrinsic thermodynamic conformational degeneracy of folding. Overall, the model combining loop-extrusion and polymer phase-separation provides a very good description of the data, particularly higher-order contacts, showing that the two mechanisms can co-exist in shaping chromatin architecture in single cells.

Uncovering the statistical physics of 3D chromosomal organization using data-driven modeling

In recent years, much effort has been devoted to understanding the three-dimensional (3D) organization of the genome and how genomic structure mediates nuclear function. The development of experimental techniques that combine DNA proximity ligation with high-throughput sequencing, such as Hi-C, have substantially improved our knowledge about chromatin organization. Numerous experimental advancements, not only utilizing DNA proximity ligation but also high-resolution genome imaging (DNA tracing), have required theoretical modeling to determine the structural ensembles consistent with such data. These 3D polymer models of the genome provide an understanding of the physical mechanisms governing genome architecture. Here, we present an overview of the recent advances in modeling the ensemble of 3D chromosomal structures by employing the maximum entropy approach combined with polymer physics. Particularly, we discuss the minimal chromatin model (MiChroM) along with the "maximum entropy genomic annotations from biomarkers associated with structural ensembles" (MEGABASE) model, which have been remarkably successful in the accurate modeling of chromosomes consistent with both Hi-C and DNA-tracing data.

Chromatin dynamics controls epigenetic domain formation

In multicellular organisms, nucleosomes carry epigenetic information that defines distinct patterns of gene expression, which are inherited over multiple generations. The enhanced capacity for information storage arises by nucleosome modifications, which are triggered by enzymes. Modified nucleosomes can transfer the mark to others that are in proximity by a positive-feedback (modification begets modification) mechanism. We created a generic polymer model, referred to as 3DSpreader, in which each bead, representing a nucleosome, stochastically switches between unmodified (U) and modified (M) states depending on the states of the neighbors. Modification begins at a specific nucleation site (NS) that is permanently in the M state, and could spread to other loci that is dictated by chromatin dynamics. Transfer of marks among the non-nucleation loci occurs stochastically as chromatin evolves in time. If the spreading rate is slower than the chromatin relaxation rate, which is biologically pertinent, then finite-sized domains form, driven by contacts between nucleosomes through a three-dimensional looping mechanism. Surprisingly, simulations based on the 3DSpreader model result in finite bounded domains that arise without the need for any boundary elements. Maintenance of spatially and temporally stable domains requires the presence of the NS, whose removal eliminates finite-sized modified domains. The theoretical predictions are in excellent agreement with experimental data for H3K9me3 spreading in mouse embryonic stem cells.

Uncovering the Principles of Genome Folding by 3D Chromatin Modeling

Our understanding of how genomic DNA is tightly packed inside the nucleus, yet is still accessible for vital cellular processes, has grown dramatically over recent years with advances in microscopy and genomics technologies. Computational methods have played a pivotal role in the structural interpretation of experimental data, which helped unravel some organizational principles of genome folding. Here, we give an overview of current computational efforts in mechanistic and data-driven 3D chromatin structure modeling. We discuss strengths and limitations of different methods and evaluate the added value and benefits of computational approaches to infer the 3D structural and dynamic properties of the genome and its underlying mechanisms at different scales and resolution, ranging from the dynamic formation of chromatin loops and topological associated domains to nuclear compartmentalization of chromatin and nuclear bodies.

Dynamic and facilitated binding of topoisomerase accelerates topological relaxation

How type 2 Topoisomerase (TopoII) proteins relax and simplify the topology of DNA molecules is one of the most intriguing open questions in genome and DNA biophysics. Most of the existing models neglect the dynamics of TopoII which is expected of proteins searching their targets via facilitated diffusion. Here, we show that dynamic binding of TopoII speeds up the topological relaxation of knotted substrates by enhancing the search of the knotted arc. Intriguingly, this in turn implies that the timescale of topological relaxation is virtually independent of the substrate length. We then discover that considering binding biases due to facilitated diffusion on looped substrates steers the sampling of the topological space closer to the boundaries between different topoisomers yielding an optimally fast topological relaxation. We discuss our findings in the context of topological simplification in vitro and in vivo.

Effects of alignment activity on the collapse kinetics of a flexible polymer

The dynamics of various biological filaments can be understood within the framework of active polymer models. Here we consider a bead-spring model for a flexible polymer chain in which the active interaction among the beads is introduced via an alignment rule adapted from the Vicsek model. Following quenching from the high-temperature coil phase to a low-temperature state point, we study the coarsening kinetics via molecular dynamics (MD) simulations using the Langevin thermostat. For the passive polymer case the low-temperature equilibrium state is a compact globule. The results from our MD simulations reveal that though the globular state is also the typical final state in the active case, the nonequilibrium pathways to arrive at such a state differ from the picture for the passive case due to the alignment interaction among the beads. We notice that deviations from the intermediate "pearl-necklace"-like arrangement, which is observed in the passive case, and the formation of more elongated dumbbell-like structures increase with increasing activity. Furthermore, it appears that while a small active force on the beads certainly makes the coarsening process much faster, there exists a nonmonotonic dependence of the collapse time on the strength of active interaction. We quantify these observations by comparing the scaling laws for the collapse time and growth of pearls with the passive case.

Active Topological Glass Confined within a Spherical Cavity

We study active topological glass under spherical confinement, allowing us to exceed the chain lengths simulated previously and determine the critical exponents of the arrested conformations. We find a previously unresolved "tank-treading" dynamic mode of active segments along the ring contour. This mode can enhance active-passive phase separation in the state of active topological glass when both diffusional and conformational relaxation of the rings are significantly suppressed. Within the observational time, we see no systematic trends in the positioning of the separated active domains within the confining sphere. The arrested state exhibits coherent stochastic rotations. We discuss possible connections of the conformational and dynamic features of the system to chromosomes enclosed in the nucleus of a living cell.

A Liquid State Perspective on Dynamics of Chromatin Compartments

The interior of the eukaryotic cell nucleus has a crowded and heterogeneous environment packed with chromatin polymers, regulatory proteins, and RNA molecules. Chromatin polymer, assisted by epigenetic modifications, protein and RNA binders, forms multi-scale compartments which help regulate genes in response to cellular signals. Furthermore, chromatin compartments are dynamic and tend to evolve in size and composition in ways that are not fully understood. The latest super-resolution imaging experiments have revealed a much more dynamic and stochastic nature of chromatin compartments than was appreciated before. An emerging mechanism explaining chromatin compartmentalization dynamics is the phase separation of protein and nucleic acids into membraneless liquid condensates. Consequently, concepts and ideas from soft matter and polymer systems have been rapidly entering the lexicon of cell biology. In this respect, the role of computational models is crucial for establishing a rigorous and quantitative foundation for the new concepts and disentangling the complex interplay of forces that contribute to the emergent patterns of chromatin dynamics and organization. Several multi-scale models have emerged to address various aspects of chromatin dynamics, ranging from equilibrium polymer simulations, hybrid non-equilibrium simulations coupling protein binding and chromatin folding, and mesoscopic field-theoretic models. Here, we review these emerging theoretical paradigms and computational models with a particular focus on chromatin’s phase separation and liquid-like properties as a basis for nuclear organization and dynamics.

Spatial organization of chromosomes leads to heterogeneous chromatin motion and drives the liquid- or gel-like dynamical behavior of chromatin

Chromosome organization and dynamics are involved in regulating many fundamental processes such as gene transcription and DNA repair. Experiments unveiled that chromatin motion is highly heterogeneous inside cell nuclei, ranging from a liquid-like, mobile state to a gel-like, rigid regime. Using polymer modeling, we investigate how these different physical states and dynamical heterogeneities may emerge from the same structural mechanisms. We found that the formation of topologically associating domains (TADs) is a key driver of chromatin motion heterogeneity. In particular, we showed that the local degree of compaction of the TAD regulates the transition from a weakly compact, fluid state of chromatin to a more compact, gel state exhibiting anomalous diffusion and coherent motion. Our work provides a comprehensive study of chromosome dynamics and a unified view of chromatin motion enabling interpretation of the wide variety of dynamical behaviors observed experimentally across different biological conditions, suggesting that the "liquid" or "solid" state of chromatin are in fact two sides of the same coin.