Statistics of randomly cross-linked polymer models to interpret chromatin conformation capture data

Chromatin phase separated nanoregions explored by polymer cross-linker models and reconstructed from single particle trajectories

Phase separated domains (PSDs) are ubiquitous in cell biology, representing nanoregions of high molecular concentration. PSDs appear at diverse cellular domains, such as neuronal synapses but also in eukaryotic cell nucleus, limiting the access of transcription factors and thus preventing gene expression. We develop a generalized cross-linker polymer model, to study PSDs: we show that increasing the number of cross-linkers induces a polymer condensation, preventing access of diffusing molecules. To investigate how the PSDs restrict the motion of diffusing molecules, we compute the mean residence and first escaping times. Finally, we develop a method based on mean-square-displacement of single particle trajectories to reconstruct the properties of PSDs from the continuum range of anomalous exponents. We also show here that PSD generated by polymers do not induces a long-range attracting field (potential well), in contrast with nanodomains at neuronal synapses. To conclude, PSDs can result from condensed chromatin organization, where the number of cross-linkers controls molecular access.

Crumpled polymer with loops recapitulates key features of chromosome organization

Chromosomes are exceedingly long topologically-constrained polymers compacted in a cell nucleus. We recently suggested that chromosomes are organized into loops by an active process of loop extrusion. Yet loops remain elusive to direct observations in living cells; detection and characterization of myriads of such loops is a major challenge. The lack of a tractable physical model of a polymer folded into loops limits our ability to interpret experimental data and detect loops. Here, we introduce a new physical model - a polymer folded into a sequence of loops, and solve it analytically. Our model and a simple geometrical argument show how loops affect statistics of contacts in a polymer across different scales, explaining universally observed shapes of the contact probability. Moreover, we reveal that folding into loops reduces the density of topological entanglements, a novel phenomenon we refer as "the dilution of entanglements". Supported by simulations this finding suggests that up to ~ 1 - 2Mb chromosomes with loops are not topologically constrained, yet become crumpled at larger scales. Our theoretical framework allows inference of loop characteristics, draws a new picture of chromosome organization, and shows how folding into loops affects topological properties of crumpled polymers.

Multi-feature clustering of CTCF binding creates robustness for loop extrusion blocking and Topologically Associating Domain boundaries

Topologically Associating Domains (TADs) separate vertebrate genomes into insulated regulatory neighborhoods that focus genome-associated processes. TADs are formed by Cohesin-mediated loop extrusion, with many TAD boundaries consisting of clustered binding sites of the CTCF insulator protein. Here we determine how this clustering of CTCF binding contributes to the blocking of loop extrusion and the insulation between TADs. We identify enrichment of three features of CTCF binding at strong TAD boundaries, consisting of strongly bound and closely spaced CTCF binding peaks, with a further enrichment of DNA-binding motifs within these peaks. Using multi-contact Nano-C analysis in cells with normal and perturbed CTCF binding, we establish that individual CTCF binding sites contribute to the blocking of loop extrusion, but in an incomplete manner. When clustered, individual CTCF binding sites thus create a stepwise insulation between neighboring TADs. Based on these results, we propose a model whereby multiple instances of temporal loop extrusion blocking create strong insulation between TADs.

A guide for single-particle chromatin tracking in live cell nuclei

The emergence of labeling strategies and live cell imaging methods enables the imaging of chromatin in living cells at single digit nanometer resolution as well as milliseconds temporal resolution. These technical breakthroughs revolutionize our understanding of chromatin structure, dynamics and functions. Single molecule tracking algorithms are usually preferred to quantify the movement of these intranucleus elements to interpret the spatiotemporal evolution of the chromatin. In this review, we will first summarize the fluorescent labeling strategy of chromatin in live cells which will be followed by a sys-tematic comparison of live cell imaging instrumentation. With the proper microscope, we will discuss the image analysis pipelines to extract the biophysical properties of the chromatin. Finally, we expect to give practical suggestions to broad biologists on how to select methods and link to the model properly according to different investigation pur-poses. This article is protected by copyright. All rights reserved.

Post-replicative pairing of sister ter regions in Escherichia coli involves multiple activities of MatP

The ter region of the bacterial chromosome, where replication terminates, is the last to be segregated before cell division in Escherichia coli. Delayed segregation is controlled by the MatP protein, which binds to specific sites (matS) within ter, and interacts with other proteins such as ZapB. Here, we investigate the role of MatP by combining short-time mobility analyses of the ter locus with biochemical approaches. We find that ter mobility is similar to that of a non ter locus, except when sister ter loci are paired after replication. This effect depends on MatP, the persistence of catenanes, and ZapB. We characterise MatP/DNA complexes and conclude that MatP binds DNA as a tetramer, but bridging matS sites in a DNA-rich environment remains infrequent. We propose that tetramerisation of MatP links matS sites with ZapB and/or with non-specific DNA to promote optimal pairing of sister ter regions until cell division.

Protein, MatP, binds to and delays segregation of the ter region of the bacterial chromosome before cell division. Here, the authors show that MatP displays multiple activities to promote optimal pairing of sister ter regions until cell division.

A Lamin-Associated Chromatin Model for Chromosome Organization

We propose a simple model for chromatin organization based on the interaction of the chromatin fibers with lamin proteins along the nuclear membrane. Lamin proteins are known to be a major factor that influences chromatin organization and hence gene expression in the cells. We provide a quantitative understanding of lamin-associated chromatin organization in a crowded macromolecular environment by systematically varying the heteropolymer segment distribution and the strength of the lamin-chromatin attractive interaction. Our minimal polymer model reproduces the formation of lamin-associated-domains and provides an in silico tool for quantifying domain length distributions for different distributions of heteropolymer segments. We show that a Gaussian distribution of heteropolymer segments, coupled with strong lamin-chromatin interactions, can qualitatively reproduce observed length distributions of lamin-associated-domains. Further, lamin-mediated interaction can enhance the formation of chromosome territories as well as the organization of chromatin into tightly packed heterochromatin and the loosely packed gene-rich euchromatin regions.

Conformational heterogeneity in human interphase chromosome organization reconciles the FISH and Hi-C paradox

Hi-C experiments are used to infer the contact probabilities between loci separated by varying genome lengths. Contact probability should decrease as the spatial distance between two loci increases. However, studies comparing Hi-C and FISH data show that in some cases the distance between one pair of loci, with larger Hi-C readout, is paradoxically larger compared to another pair with a smaller value of the contact probability. Here, we show that the FISH-Hi-C paradox can be resolved using a theory based on a Generalized Rouse Model for Chromosomes (GRMC). The FISH-Hi-C paradox arises because the cell population is highly heterogeneous, which means that a given contact is present in only a fraction of cells. Insights from the GRMC is used to construct a theory, without any adjustable parameters, to extract the distribution of subpopulations from the FISH data, which quantitatively reproduces the Hi-C data. Our results show that heterogeneity is pervasive in genome organization at all length scales, reflecting large cell-to-cell variations.

Studies comparing Hi-C and FISH data show that in some cases the distance between one pair of loci is paradoxically larger compared to another pair with a smaller value of the contact probability. Here the authors use a theory based on a Generalized Rouse Model for Chromosomes to resolve this paradox.

Advances Using Single-Particle Trajectories to Reconstruct Chromatin Organization and Dynamics

Chromatin organization remains complex and far from understood. In this article, we review recent statistical methods of extracting biophysical parameters from in vivo single-particle trajectories of loci to reconstruct chromatin reorganization in response to cellular stress such as DNA damage. We look at methods for analyzing both single locus and multiple loci tracked simultaneously and explain how to quantify and describe chromatin motion using a combination of extractable parameters. These parameters can be converted into information about chromatin dynamics and function. Furthermore, we discuss how the timescale of recurrent encounter between loci can be extracted and interpreted. We also discuss the effect of sampling rate on the estimated parameters. Finally, we review a polymer method to reconstruct chromatin structure using crosslinkers between chromatin sites. We list and refer to some software packages that are now publicly available to simulate polymer motion. To conclude, chromatin organization and dynamics can be reconstructed from locus trajectories and predicted based on polymer models.

Statistics of chromatin organization during cell differentiation revealed by heterogeneous cross-linked polymers

Chromatin of mammalian nucleus folds into discrete contact enriched regions such as Topologically Associating Domains (TADs). Folding hierarchy and internal organization of TADs is highly dynamic throughout cellular differentiation, and are correlated with gene activation and silencing. To account for multiple interacting TADs, we developed a parsimonious randomly cross-linked (RCL) polymer model that maps high frequency Hi-C encounters within and between TADs into direct loci interactions using cross-links at a given base-pair resolution. We reconstruct three TADs of the mammalian X chromosome for three stages of differentiation. We compute the radius of gyration of TADs and the encounter probability between genomic segments. We found 1) a synchronous compaction and decompaction of TADs throughout differentiation and 2) high order organization into meta-TADs resulting from weak inter-TAD interactions. Finally, the present framework allows to infer transient properties of the chromatin from steady-state statistics embedded in the Hi-C/5C data.

Chromatin is folded into Topologically Associating domains (TADs), with the organization and folding hierarchy of the TADs being highly dynamic. Here the authors develop a parsimonious randomly cross-linked (RCL) polymer model that maps high frequency encounters present in Hi-C data within and between TADs and reconstruct TADs across cell differentiation, revealing local chromatin re-organization.

Redundancy principle and the role of extreme statistics in molecular and cellular biology

The paradigm of chemical activation rates in cellular biology has been shifted from the mean arrival time of a single particle to the mean of the first among many particles to arrive at a small activation site. The activation rate is set by extremely rare events, which have drastically different time scales from the mean times between activations, and depends on different structural parameters. This shift calls for reconsideration of physical processes used in deterministic and stochastic modeling of chemical reactions that are based on the traditional forward rate, especially for fast activation processes in living cells. Consequently, the biological activation time is not necessarily exponentially distributed. We review here the physical models, the mathematical analysis and the new paradigm of setting the scale to be the shortest time for activation that clarifies the role of population redundancy in selecting and accelerating transient cellular search processes. We provide examples in cellular transduction, gene activation, cell senescence activation or spermatozoa selection during fertilization, where the rate depends on numbers. We conclude that the statistics of the minimal time to activation set kinetic laws in biology, which can be very different from the ones associated to average times.

Interphase human chromosome exhibits out of equilibrium glassy dynamics

Fingerprints of the three-dimensional organization of genomes have emerged using advances in Hi-C and imaging techniques. However, genome dynamics is poorly understood. Here, we create the chromosome copolymer model (CCM) by representing chromosomes as a copolymer with two epigenetic loci types corresponding to euchromatin and heterochromatin. Using novel clustering techniques, we establish quantitatively that the simulated contact maps and topologically associating domains (TADs) for chromosomes 5 and 10 and those inferred from Hi-C experiments are in good agreement. Chromatin exhibits glassy dynamics with coherent motion on micron scale. The broad distribution of the diffusion exponents of the individual loci, which quantitatively agrees with experiments, is suggestive of highly heterogeneous dynamics. This is reflected in the cell-to-cell variations in the contact maps. Chromosome organization is hierarchical, involving the formation of chromosome droplets (CDs) on genomic scale, coinciding with the TAD size, followed by coalescence of the CDs, reminiscent of Ostwald ripening.

Two loci single particle trajectories analysis: constructing a first passage time statistics of local chromatin exploration

Stochastic single particle trajectories are used to explore the local chromatin organization. We present here a statistical analysis of the first contact time distributions between two tagged loci recorded experimentally. First, we extract the association and dissociation times from data for various genomic distances between loci, and we show that the looping time occurs in confined nanometer regions. Second, we characterize the looping time distribution for two loci in the presence of multiple DNA damages. Finally, we construct a polymer model, that accounts for the local chromatin organization before and after a double-stranded DNA break (DSB), to estimate the level of chromatin decompaction. This novel passage time statistics method allows extracting transient dynamic at scales varying from one to few hundreds of nanometers, it predicts the local changes in the number of binding molecules following DSB and can be used to characterize the local dynamic of the chromatin.