Predicting chromatin architecture from models of polymer physics

Multiscale molecular modeling of chromatin with MultiMM: From nucleosomes to the whole genome

Motivation: We present a user-friendly 3D chromatin simulation model for the human genome based on OpenMM, addressing the challenges posed by existing models with use-specific implementations. Our approach employs a multi-scale energy minimization strategy, capturing chromatin's hierarchical structure. Initiating with a Hilbert curve-based structure, users can input files specifying nucleosome positioning, loops, compartments, or subcompartments. Results: The model utilizes an energy minimization approach with a large choice of numerical integrators, providing the entire genome's structure within minutes. Output files include the generated structures for each chromosome, offering a versatile and accessible tool for chromatin simulation in bioinformatics studies. Furthermore, MultiMM is capable of producing nucleosome-resolution structures by making simplistic geometric assumptions about the structure and the density of nucleosomes on the DNA. Code availability: Open-source software and the manual are freely available on https://github.com/SFGLab/MultiMM or via pip https://pypi.org/project/MultiMM/.

Chromosome structure modeling tools and their evaluation in bacteria

The three-dimensional (3D) structure of bacterial chromosomes is crucial for understanding chromosome function. With the growing availability of high-throughput chromosome conformation capture (3C/Hi-C) data, the 3D structure reconstruction algorithms have become powerful tools to study bacterial chromosome structure and function. It is highly desired to have a recommendation on the chromosome structure reconstruction tools to facilitate the prokaryotic 3D genomics. In this work, we review existing chromosome 3D structure reconstruction algorithms and classify them based on their underlying computational models into two categories: constraint-based modeling and thermodynamics-based modeling. We briefly compare these algorithms utilizing 3C/Hi-C datasets and fluorescence microscopy data obtained from Escherichia coli and Caulobacter crescentus, as well as simulated datasets. We discuss current challenges in the 3D reconstruction algorithms for bacterial chromosomes, primarily focusing on software usability. Finally, we briefly prospect future research directions for bacterial chromosome structure reconstruction algorithms.

Fundamental insights into the correlation between chromosome configuration and transcription

Eukaryotic chromosomes exhibit a hierarchical organization that spans a spectrum of length scales, ranging from sub-regions known as loops, which typically comprise hundreds of base pairs, to much larger chromosome territories that can encompass a few mega base pairs. Chromosome conformation capture experiments that involve high-throughput sequencing methods combined with microscopy techniques have enabled a new understanding of inter- and intra-chromosomal interactions with unprecedented details. This information also provides mechanistic insights on the relationship between genome architecture and gene expression. In this article, we review the recent findings on three-dimensional interactions among chromosomes at the compartment, topologically associating domain, and loop levels and the impact of these interactions on the transcription process. We also discuss current understanding of various biophysical processes involved in multi-layer structural organization of chromosomes. Then, we discuss the relationships between gene expression and genome structure from perturbative genome-wide association studies. Furthermore, for a better understanding of how chromosome architecture and function are linked, we emphasize the role of epigenetic modifications in the regulation of gene expression. Such an understanding of the relationship between genome architecture and gene expression can provide a new perspective on the range of potential future discoveries and therapeutic research.

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.

A Polymer Physics Model to Dissect Genome Organization in Healthy and Pathological Phenotypes

Novel technologies revealed a nontrivial spatial organization of the chromosomes within the cell nucleus, which includes different levels of compartmentalization and architectural patterns. Notably, such complex three-dimensional structure plays a crucial role in vital biological functions and its alterations can produce extensive rewiring of genomic regulatory regions, thus leading to gene misexpression and disease. Here, we show that theoretical and computational approaches, based on polymer physics, can be employed to dissect chromatin contacts in three-dimensional space and to predict the effects of pathogenic structural variants on the genome architecture. In particular, we discuss the folding of the human EPHA4 and the murine Pitx1 loci as case studies.

The Physical Behavior of Interphase Chromosomes: Polymer Theory and Coarse-Grain Computer Simulations

Fluorescence in situ hybridization and chromosome conformation capture methods point to the same conclusion: that chromosomes appear to the external observer as compact structures with a highly nonrandom three-dimensional organization. In this work, we recapitulate the efforts made by us and other groups to rationalize this behavior in terms of the mathematical language and tools of polymer physics. After a brief introduction dedicated to some crucial experiments dissecting the structure of interphase chromosomes, we discuss at a nonspecialistic level some fundamental aspects of theoretical and numerical polymer physics. Then, we inglobe biological and polymer aspects into a polymer model for interphase chromosomes which moves from the observation that mutual topological constraints, such as those typically present between polymer chains in ordinary melts, induce slow chain dynamics and "constraint" chromosomes to resemble double-folded randomly branched polymer conformations. By explicitly turning these ideas into a multi-scale numerical algorithm which is described here in full details, we can design accurate model polymer conformations for interphase chromosomes and offer them for systematic comparison to experiments. The review is concluded by discussing the limitations of our approach and pointing to promising perspectives for future work.

Understanding 3D genome organization by multidisciplinary methods

Understanding how chromatin is folded in the nucleus is fundamental to understanding its function. Although 3D genome organization has been historically difficult to study owing to a lack of relevant methodologies, major technological breakthroughs in genome-wide mapping of chromatin contacts and advances in imaging technologies in the twenty-first century considerably improved our understanding of chromosome conformation and nuclear architecture. In this Review, we discuss methods of 3D genome organization analysis, including sequencing-based techniques, such as Hi-C and its derivatives, Micro-C, DamID and others; microscopy-based techniques, such as super-resolution imaging coupled with fluorescence in situ hybridization (FISH), multiplex FISH, in situ genome sequencing and live microscopy methods; and computational and modelling approaches. We describe the most commonly used techniques and their contribution to our current knowledge of nuclear architecture and, finally, we provide a perspective on up-and-coming methods that open possibilities for future major discoveries.

Divide and Rule: Phase Separation in Eukaryotic Genome Functioning

The functioning of a cell at various organizational levels is determined by the interactions between macromolecules that promote cellular organelle formation and orchestrate metabolic pathways via the control of enzymatic activities. Although highly specific and relatively stable protein-protein, protein-DNA, and protein-RNA interactions are traditionally suggested as the drivers for cellular function realization, recent advances in the discovery of weak multivalent interactions have uncovered the role of so-called macromolecule condensates. These structures, which are highly divergent in size, composition, function, and cellular localization are predominantly formed by liquid-liquid phase separation (LLPS): a physical-chemical process where an initially homogenous solution turns into two distinct phases, one of which contains the major portion of the dissolved macromolecules and the other one containing the solvent. In a living cell, LLPS drives the formation of membrane-less organelles such as the nucleolus, nuclear bodies, and viral replication factories and facilitates the assembly of complex macromolecule aggregates possessing regulatory, structural, and enzymatic functions. Here, we discuss the role of LLPS in the spatial organization of eukaryotic chromatin and regulation of gene expression in normal and pathological conditions.

Computational approaches for inferring 3D conformations of chromatin from chromosome conformation capture data

Chromosome conformation capture (3C) and its variants are powerful experimental techniques for probing intra- and inter-chromosomal interactions within cell nuclei at high resolution and in a high-throughput, quantitative manner. The contact maps derived from such experiments provide an avenue for inferring the 3D spatial organization of the genome. This review provides an overview of the various computational methods developed in the past decade for addressing the very important but challenging problem of deducing the detailed 3D structure or structure population of chromosomal domains, chromosomes, and even entire genomes from 3C contact maps.

Chromatin Compaction Multiscale Modeling: A Complex Synergy Between Theory, Simulation, and Experiment

Understanding the mechanisms that trigger chromatin compaction, its patterns, and the factors they depend on, is a fundamental and still open question in Biology. Chromatin compacts and reinforces DNA and is a stable but dynamic structure, to make DNA accessible to proteins. In recent years, computational advances have provided larger amounts of data and have made large-scale simulations more viable. Experimental techniques for the extraction and reconstitution of chromatin fibers have improved, reinvigorating theoretical and experimental interest in the topic and stimulating debate on points previously considered as certainties regarding chromatin. A great assortment of approaches has emerged, from all-atom single-nucleosome or oligonucleosome simulations to various degrees of coarse graining, to polymer models, to fractal-like structures and purely topological models. Different fiber-start patterns have been studied in theory and experiment, as well as different linker DNA lengths. DNA is a highly charged macromolecule, making ionic and electrostatic interactions extremely important for chromatin topology and dynamics. Indeed, the repercussions of varying ionic concentration have been extensively examined at the computational level, using all-atom, coarse-grained, and continuum techniques. The presence of high-curvature AT-rich segments in DNA can cause conformational variations, attesting to the fact that the role of DNA is both structural and electrostatic. There have been some tentative attempts to describe the force fields governing chromatin conformational changes and the energy landscapes of these transitions, but the intricacy of the system has hampered reaching a consensus. The study of chromatin conformations is an intrinsically multiscale topic, influenced by a wide range of biological and physical interactions, spanning from the atomic to the chromosome level. Therefore, powerful modeling techniques and carefully planned experiments are required for an overview of the most relevant phenomena and interactions. The topic provides fertile ground for interdisciplinary studies featuring a synergy between theoretical and experimental scientists from different fields and the cross-validation of respective results, with a multi-scale perspective. Here, we summarize some of the most representative approaches, and focus on the importance of electrostatics and solvation, often overlooked aspects of chromatin modeling.

Physical and data structure of 3D genome

With the textbook view of chromatin folding based on the 30-nm fiber being challenged, it has been proposed that interphase DNA has an irregular 10-nm nucleosome polymer structure whose folding philosophy is unknown. Nevertheless, experimental advances suggest that this irregular packing is associated with many nontrivial physical properties that are puzzling from a polymer physics point of view. Here, we show that the reconciliation of these exotic properties necessitates modularizing three-dimensional genome into tree data structures on top of, and in striking contrast to, the linear topology of DNA double helix. These functional modules need to be connected and isolated by an open backbone that results in porous and heterogeneous packing in a quasi-self-similar manner, as revealed by our electron and optical imaging. Our multiscale theoretical and experimental results suggest the existence of higher-order universal folding principles for a disordered chromatin fiber to avoid entanglement and fulfill its biological functions.

Higher-order Chromosome Structures Investigated by Polymer Physics in Cellular Morphogenesis and Differentiation

Experimental advances in Molecular Biology demonstrated that chromatin architecture and gene regulation are deeply related. Hi-C data, for instance, returned a scenario where chromosomes form a complex pattern of interactions, including TADs, metaTADs, and compartments, correlated with genomic and epigenomic features. Here, we discuss the emerging hierarchical organization of chromatin and show how it remains partially conserved during mouse neuronal differentiation with changes highly related to modifications in gene expression. In this scenario, models of polymer physics, such as the Strings & Binders (SBS) model, can be a crucial instrument to understand the molecular mechanisms underlying the formation of such a higher order 3D structure. In particular, we focus on the case study of the murine Pitx1 genomic region. At this locus, two alternative spatial conformations take place in the hindlimb and forelimb tissues, corresponding to two different transcriptional states of Pitx1. We finally show how the structural variants can affect the locus 3D organization leading to ectopic gene expression and limb malformations.

Merging 1D and 3D genomic information: Challenges in modelling and validation

Genome organization in eukaryotes during interphase stems from the delicate balance between non-random correlations present in the DNA polynucleotide linear sequence and the physico/chemical reactions which shape continuously the form and structure of DNA and chromatin inside the nucleus of the cell. It is now clear that these mechanisms have a key role in important processes like gene regulation, yet the detailed ways they act simultaneously and, eventually, come to influence each other even across very different length-scales remain largely unexplored. In this paper, we recapitulate some of the main results concerning gene regulatory and physical mechanisms, in relation to the information encoded in the 1D sequence and the 3D folding structure of DNA. In particular, we stress how reciprocal crossfeeding between 1D and 3D models may provide original insight into how these complex processes work and influence each other. This article is part of a Special Issue entitled: Transcriptional Profiles and Regulatory Gene Networks edited by Dr. Dr. Federico Manuel Giorgi and Dr. Shaun Mahony.

Whole-Cell Models and Simulations in Molecular Detail

Comprehensive data about the composition and structure of cellular components have enabled the construction of quantitative whole-cell models. While kinetic network-type models have been established, it is also becoming possible to build physical, molecular-level models of cellular environments. This review outlines challenges in constructing and simulating such models and discusses near- and long-term opportunities for developing physical whole-cell models that can connect molecular structure with biological function.

Microrheology of interphase chromosomes with spatial constraints: a computational study

Chromatin fibers within the interior of the nucleus of the cell make stable interactions with the nucleoskeleton, an ensemble of 'extra-chromatin' structures which help ensuring genome stability. Although the role of these interactions appears crucial to the correct behavior of the cell, their impact on chromatin structure and dynamics remains to be elucidated. In order to tackle this important issue, in this work we introduce a simple polymer model for chromatin fibers in interphase which takes into account the two generic properties of chain-versus-chain mutual uncrossability and the presence of stable binding interactions to an extra-chromatin nuclear matrix. To study how these constraints affect chromatin structure from small to large scales, we employ extensive molecular dynamics computer simulations and we monitor the motion of nanoprobes of different sizes embedded within the polymer medium. Our results demonstrate that nanoprobes show hampered motion whenever their linear size becomes larger than chromatin stiffness. This transition is also displaying features which usually belong to the realm of glassy systems, namely long-tail correlations in the distribution functions of nanoprobe spatial displacements and heterogeneous behavior accompanied by ergodicity breaking.

Transcription-driven genome organization: a model for chromosome structure and the regulation of gene expression tested through simulations

Current models for the folding of the human genome see a hierarchy stretching down from chromosome territories, through A/B compartments and topologically-associating domains (TADs), to contact domains stabilized by cohesin and CTCF. However, molecular mechanisms underlying this folding, and the way folding affects transcriptional activity, remain obscure. Here we review physical principles driving proteins bound to long polymers into clusters surrounded by loops, and present a parsimonious yet comprehensive model for the way the organization determines function. We argue that clusters of active RNA polymerases and their transcription factors are major architectural features; then, contact domains, TADs and compartments just reflect one or more loops and clusters. We suggest tethering a gene close to a cluster containing appropriate factors—a transcription factory—increases the firing frequency, and offer solutions to many current puzzles concerning the actions of enhancers, super-enhancers, boundaries and eQTLs (expression quantitative trait loci). As a result, the activity of any gene is directly influenced by the activity of other transcription units around it in 3D space, and this is supported by Brownian-dynamics simulations of transcription factors binding to cognate sites on long polymers.

How the Genome Folds: The Biophysics of Four-Dimensional Chromatin Organization

The genetic information that instructs transcription and other cellular functions is carried by the chromosomes, polymers of DNA in complex with histones and other proteins. These polymers are folded inside nuclei five orders of magnitude smaller than their linear length, and many facets of this folding correlate with or are causally related to transcription and other cellular functions. Recent advances in sequencing and imaging-based techniques have enabled new views into several layers of chromatin organization. These experimental findings are accompanied by computational modeling efforts based on polymer physics that can provide mechanistic insights and quantitative predictions. Here, we review current knowledge of the main levels of chromatin organization, from the scale of nucleosomes to the entire nucleus, our current understanding of their underlying biophysical and molecular mechanisms, and some of their functional implications.

Mapping the spectrum of 3D communities in human chromosome conformation capture data

Several experiments show that the three dimensional (3D) organization of chromosomes affects genetic processes such as transcription and gene regulation. To better understand this connection, researchers developed the Hi-C method that is able to detect the pairwise physical contacts of all chromosomal loci. The Hi-C data show that chromosomes are composed of 3D compartments that range over a variety of scales. However, it is challenging to systematically detect these cross-scale structures. Most studies have therefore designed methods for specific scales to study foremost topologically associated domains (TADs) and A/B compartments. To go beyond this limitation, we tailor a network community detection method that finds communities in compact fractal globule polymer systems. Our method allows us to continuously scan through all scales with a single resolution parameter. We found: (i) polymer segments belonging to the same 3D community do not have to be in consecutive order along the polymer chain. In other words, several TADs may belong to the same 3D community. (ii) CTCF proteins-a loop-stabilizing protein that is ascribed a big role in TAD formation-are well correlated with community borders only at one level of organization. (iii) TADs and A/B compartments are traditionally treated as two weakly related 3D structures and detected with different algorithms. With our method, we detect both by simply adjusting the resolution parameter. We therefore argue that they represent two specific levels of a continuous spectrum 3D communities, rather than seeing them as different structural entities.

Spontaneous domain formation in disordered copolymers as a mechanism for chromosome structuring

Motivated by the problem of domain formation in chromosomes, we studied a co-polymer model where only a subset of the monomers feel attractive interactions. These monomers are displaced randomly from a regularly-spaced pattern, thus introducing some quenched disorder in the system. Previous work has shown that in the case of regularly-spaced interacting monomers this chain can fold into structures characterized by multiple distinct domains of consecutive segments. In each domain, attractive interactions are balanced by the entropy cost of forming loops. We show by advanced replica-exchange simulations that adding disorder in the position of the interacting monomers further stabilizes these domains. The model suggests that the partitioning of the chain into well-defined domains of consecutive monomers is a spontaneous property of heteropolymers. In the case of chromosomes, evolution could have acted on the spacing of interacting monomers to modulate in a simple way the underlying domains for functional reasons.

Enhancer functions in three dimensions: beyond the flat world perspective

Transcriptional enhancers constitute a subclass of regulatory elements that facilitate transcription. Such regions are generally organized by short stretches of DNA enriched in transcription factor-binding sites but also can include very large regions containing clusters of enhancers, termed super-enhancers. These regions increase the probability or the rate (or both) of transcription generally in cis and sometimes over very long distances by altering chromatin states and the activity of Pol II machinery at promoters. Although enhancers were discovered almost four decades ago, their inner workings remain enigmatic. One important opening into the underlying principle has been provided by observations that enhancers make physical contacts with their target promoters to facilitate the loading of the RNA polymerase complex. However, very little is known about how such chromatin loops are regulated and how they govern transcription in the three-dimensional context of the nuclear architecture. Here, we present current themes of how enhancers may boost gene expression in three dimensions and we identify currently unresolved key questions.

Molecular Dynamics simulations of the Strings and Binders Switch model of chromatin

In recent years interest has grown on the applications of polymer physics to model chromatin folding in order to try to make sense of the complexity of experimental data emerging from new technologies such as Hi-C or GAM, in a principled way. Here we review the methods employed to efficiently implement Molecular Dynamics computer simulations of polymer models, focusing in particular on the String&Binders Switch (SBS) model. The constant improvement of such methods and computer power is returning increasingly more accurate insights on the structure and molecular mechanisms underlying the spatial organization of chromosomes in the cell nucleus. We aim to provide an account of the state of the art of computational techniques employed in this type of investigations and to review recent applications of such methods to the description of real genomic loci, such as the Sox9 locus in mESC.

Structure of the human chromosome interaction network

New Hi-C technologies have revealed that chromosomes have a complex network of spatial contacts in the cell nucleus of higher organisms, whose organisation is only partially understood. Here, we investigate the structure of such a network in human GM12878 cells, to derive a large scale picture of nuclear architecture. We find that the intensity of intra-chromosomal interactions is power-law distributed. Inter-chromosomal interactions are two orders of magnitude weaker and exponentially distributed, yet they are not randomly arranged along the genomic sequence. Intra-chromosomal contacts broadly occur between epigenomically homologous regions, whereas inter-chromosomal contacts are especially associated with regions rich in highly expressed genes. Overall, genomic contacts in the nucleus appear to be structured as a network of networks where a set of strongly individual chromosomal units, as envisaged in the ‘chromosomal territory’ scenario derived from microscopy, interact with each other via on average weaker, yet far from random and functionally important interactions.

A Polymer Physics Investigation of the Architecture of the Murine Orthologue of the 7q11.23 Human Locus

In the last decade, the developments of novel technologies, such as Hi-C or GAM methods, allowed to discover that chromosomes in the nucleus of mammalian cells have a complex spatial organization, encompassing the functional contacts between genes and regulators. In this work, we review recent progresses in chromosome modeling based on polymer physics to understand chromatin structure and folding mechanisms. As an example, we derive in mouse embryonic stem cells the full 3D structure of the Bmp7 locus, a genomic region that plays a key role in osteoblastic differentiation. Next, as an application to Neuroscience, we present the first 3D model for the mouse orthologoue of the Williams-Beuren syndrome 7q11.23 human locus. Deletions and duplications of the 7q11.23 region generate neurodevelopmental disorders with multi-system involvement and variable expressivity, and with autism. Understanding the impact of such mutations on the rewiring of the interactions of genes and regulators could be a new key to make sense of their related diseases, with potential applications in biomedicine.