LoopSage: An energy-based Monte Carlo approach for the loop extrusion modeling of chromatin

The connection between the patterns observed in 3C-type experiments and the modeling of polymers remains unresolved. This paper presents a simulation pipeline that generates thermodynamic ensembles of 3D structures for topologically associated domain (TAD) regions by loop extrusion model (LEM). The simulations consist of two main components: a stochastic simulation phase, employing a Monte Carlo approach to simulate the binding positions of cohesins, and a dynamical simulation phase, utilizing these cohesins' positions to create 3D structures. In this approach, the system's total energy is the combined result of the Monte Carlo energy and the molecular simulation energy, which are iteratively updated. The structural maintenance of chromosomes (SMC) protein complexes are represented as loop extruders, while the CCCTC-binding factor (CTCF) locations on DNA sequence are modeled as energy minima on the Monte Carlo energy landscape. Finally, the spatial distances between DNA segments from ChIA-PET experiments are compared with the computer simulations, and we observe significant Pearson correlations between predictions and the real data. LoopSage model offers a fresh perspective on chromatin loop dynamics, allowing us to observe phase transition between sparse and condensed states in 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.

Large-scale data-driven and physics-based models offer insights into the relationships among the structures, dynamics, and functions of chromosomes

The organized three-dimensional chromosome architecture in the cell nucleus provides scaffolding for precise regulation of gene expression. When the cell changes its identity in the cell-fate decision-making process, extensive rearrangements of chromosome structures occur accompanied by large-scale adaptations of gene expression, underscoring the importance of chromosome dynamics in shaping genome function. Over the last two decades, rapid development of experimental methods has provided unprecedented data to characterize the hierarchical structures and dynamic properties of chromosomes. In parallel, these enormous data offer valuable opportunities for developing quantitative computational models. Here, we review a variety of large-scale polymer models developed to investigate the structures and dynamics of chromosomes. Different from the underlying modeling strategies, these approaches can be classified into data-driven ('top-down') and physics-based ('bottom-up') categories. We discuss their contributions to offering valuable insights into the relationships among the structures, dynamics, and functions of chromosomes and propose the perspective of developing data integration approaches from different experimental technologies and multidisciplinary theoretical/simulation methods combined with different modeling strategies.

The dynamic role of cohesin in maintaining human genome architecture

Recent advances in genomic and imaging techniques have revealed the complex manner of organizing billions of base pairs of DNA necessary for maintaining their functionality and ensuring the proper expression of genetic information. The SMC proteins and cohesin complex primarily contribute to forming higher-order chromatin structures, such as chromosomal territories, compartments, topologically associating domains (TADs) and chromatin loops anchored by CCCTC-binding factor (CTCF) protein or other genome organizers. Cohesin plays a fundamental role in chromatin organization, gene expression and regulation. This review aims to describe the current understanding of the dynamic nature of the cohesin-DNA complex and its dependence on cohesin for genome maintenance. We discuss the current 3C technique and numerous bioinformatics pipelines used to comprehend structural genomics and epigenetics focusing on the analysis of Cohesin-centred interactions. We also incorporate our present comprehension of Loop Extrusion (LE) and insights from stochastic modelling.

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.

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.

Topological Constraints with Optimal Length Promote the Formation of Chromosomal Territories at Weakened Degree of Phase Separation

It is generally agreed that the nuclei of eukaryotic cells at interphase are partitioned into disjointed territories, with distinct regions occupied by certain chromosomes. However, the underlying mechanism for such territorialization is still under debate. Here we model chromosomes as coarse-grained block copolymers and to investigate the effect of loop domains (LDs) on the formation of compartments and territories based on dissipative particle dynamics. A critical length of LDs, which depends sensitively on the length of polymeric blocks, is obtained to minimize the degree of phase separation. This also applies to the two-polymer system: The critical length not only maximizes the degree of territorialization but also minimizes the degree of phase separation. Interestingly, by comparing with experimental data, we find the critical length for LDs and the corresponding length of blocks to be respectively very close to the mean length of topologically associating domains (TADs) and chromosomal segments with different densities of CpG islands for human chromosomes. The results indicate that topological constraints with optimal length can contribute to the formation of territories by weakening the degree of phase separation, which likely promotes the chromosomal flexibility in response to genetic regulations.

Multiscale modeling of genome organization with maximum entropy optimization

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

Coarse-grained force-field for large scale molecular dynamics simulations of polyacrylamide and polyacrylamide-gels based on quantum mechanics

Developing a coarse-grained force field for polyacrylamide based on quantum mechanics equation of state.