PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of plosonePLoS OneView this ArticleSubmit to PLoSGet E-mail AlertsContact UsPublic Library of Science (PLoS)
 
PLoS One. 2010; 5(8): e12218.
Published online 2010 August 25. doi:  10.1371/journal.pone.0012218
PMCID: PMC2928267

Diffusion-Driven Looping Provides a Consistent Framework for Chromatin Organization

Timothy Ravasi, Editor

Abstract

Chromatin folding inside the interphase nucleus of eukaryotic cells is done on multiple scales of length and time. Despite recent progress in understanding the folding motifs of chromatin, the higher-order structure still remains elusive. Various experimental studies reveal a tight connection between genome folding and function. Chromosomes fold into a confined subspace of the nucleus and form distinct territories. Chromatin looping seems to play a dominant role both in transcriptional regulation as well as in chromatin organization and has been assumed to be mediated by long-range interactions in many polymer models. However, it remains a crucial question which mechanisms are necessary to make two chromatin regions become co-located, i.e. have them in spatial proximity. We demonstrate that the formation of loops can be accomplished solely on the basis of diffusional motion. The probabilistic nature of temporary contacts mimics the effects of proteins, e.g. transcription factors, in the solvent. We establish testable quantitative predictions by deriving scale-independent measures for comparison to experimental data. In this Dynamic Loop (DL) model, the co-localization probability of distant elements is strongly increased compared to linear non-looping chains. The model correctly describes folding into a confined space as well as the experimentally observed cell-to-cell variation. Most importantly, at biological densities, model chromosomes occupy distinct territories showing less inter-chromosomal contacts than linear chains. Thus, dynamic diffusion-based looping, i.e. gene co-localization, provides a consistent framework for chromatin organization in eukaryotic interphase nuclei.

Introduction

The cell nucleus is a main constituent of eukaryotic organisms and yet its complexity prevents detailed knowledge of its function. The genome content is carried by the chromosomes: compactly folded polymers consisting of DNA and histone proteins. While during mitosis chromosomes are found in an extremely condensed state, the chromatin fiber inside the interphase nucleus has a much more decondensed organization. However, at this stage of the cell cycle, highly coordinated processes such as transcription, replication and DNA repair take place, making a random folding of the chromatin fiber very unlikely. A pivotal question is the connection between genome organization and function, which could not be answered in a satisfying way up to now. The organization of the genome in the interphase nucleus of eukaryotic cells is done on multiple scales of length and degrees of compaction. The basic filament is the DNA double helix which is wrapped around histone cores forming the nucleosome. The chromatin fiber is a complex of nucleosomes and linker DNA forming a beads-on-a-string type of filament with a diameter of about 11 nm [1]. In-vitro experiments provide evidence that this structure in turn condenses under certain salt conditions to an even more compact structure of 30 nm, but both its regularity and its existence in living cell nuclei are still under debate [2][5]. Stunningly, even less is known about the structural organization on a scale above 30 nm. Up to now, experimental techniques are limited by the resolution of conventional light microscopes of about 200 nm, requiring indirect assays for investigating chromatin folding. Several experimental techniques have been applied: Fluorescent labeling of large parts of a chromosome yields results on structure, shape and position of chromosomal regions [6] or even of entire individual chromosomes [7]. Labeling two loci of a chromosome with a fluorescent marker was successfully used to establish a relationship between genomic distance An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e001.jpg between these markers and its mean square physical distance in yeast [8], drosophila [9], [10] and human cells [11], [12].

There is now abundant evidence that genome function is tightly related to chromatin folding on several length scales. The one-dimensional distribution of genes along the chromosome is far from being random: the Human Transcriptome Map [13] reveals a clustering of active genes as well as inactive genes into certain domains, which have been named ridges and anti-ridges [6]. Various experiments have shown that the 3D organization of chromatin depends on transcriptional activity: Active genes tend to be located in the nuclear interior while inactive genes are found more often at the nuclear periphery [6], [14], [15], the converse behavior is observed in some experiments [16]. Moreover, a change in the transcriptional state of a gene can have direct influence on its positioning inside the nucleus [17], [18]. Transcriptional active regions (ridges) were observed to have a more open structure than inactive regions (anti-ridges) [6]. Also, the relationship between mean square distance (MSD) of two fluorescent markers and their genomic separation has revealed significant differences in compaction between ridges and anti-ridges [12]. Further, these fluorescence in situ hybridization (FISH) experiments displayed a leveling-off in the MSD for genomic separations above 10 Mb (mega basepairs), the plateau level being in the size range of 2 An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e002.jpgm. On the scale of the nucleus, chromosomes are separated into distinct chromosome territories [7], whose relative positions and ellipsoidal shape varies from cell to cell [19].

Intra-chromosomal as well as inter-chromosomal contacts or loops have been intensively analyzed in the past few years both experimentally and theoretically as a possible mechanism for transcriptional regulation and genome folding. Yet, chromatin loops seem to be an ubiquitous feature of genome organization and genome function. Transcriptional regulation is often controlled by regulatory motifs such as enhancers and silencers. These can be located tens of kb apart from the target gene which they regulate [20], [21]. One possible interaction mechanism is spatial proximity of regulatory element and target gene which requires the looping out of intervening DNA. 3C experiments have demonstrated that this is indeed the case in the An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e003.jpg-globin locus [22]. One idea put forward to explain chromatin loops is the existence of transcription factories or active chromatin hubs, where active polymerases cluster and thereby co-locate genes and regulatory elements [23], [24]. 3C and 4C techniques have since then provided evidence that indeed loops up to several Mb exist in interphase cells [20], [25]. However, the detailed mechanisms and driving forces of looping are still under debate.

We present a polymer model, the Dynamic Loop model, where functional loops are formed solely on the basis of diffusional motion. Importantly, loops are assumed to be dynamic and the sets of loop attachment points change during time. Thus, our loop model is minimal meaning that we do not assume a priori long-range forces and active transport mechanisms. Besides the new motif of dynamic loop formation, this is a major advance with respect to other chromatin models with loops [11], [12], [26][29]. Various other polymer models have been proposed [30], [31] which do not take into account chromatin looping. The assumptions of our model arise from biological evidence: 4C experiments clearly show that loops exist on length scales from several thousand basepairs to tens of Mb [25]. Surely, if looping is related to functional processes like transcriptional regulation and the formation of transcription factories, the cell must be able to control this looping dynamically. Large cell-to-cell variations in FISH distance measurements [12], [32] render such a dynamics a necessary feature of any polymer model.

Our model makes testable predictions on a variety of observable quantities. We predict that chromosomes fold into a confined space and display a different fluctuation regime than non-dynamic looping polymers or linear chains. Importantly, the formation of large loops can be accomplished hierarchically mediated by many loops on the short scale without the assumption of long-range interactions. We demonstrate that the beads of the polymer display sub-diffusive behavior in agreement with experimental data [17] and that chromosome territories are constituted driven by looping; the overlap between different chromosome territories (CTs) depends on the local looping probabilities.

The chromatin model

Our model starts by initially assuming chromatin to consist of a coarse-grained linear polymer chain. Loop formation is achieved on the basis of diffusional motion of the monomers in the following way: Whenever two segments co-localize by diffusional motion, a chromatin loop is formed with a certain probability An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e004.jpg between these two sites. A certain lifetime is assigned to each loop, thus loop attachment points dissolve again during the course of time. Lacking experimental knowledge on the time scales over which chromatin segments remain co-localized, e.g. in transcription factories, different looping lifetimes are considered. Details are described in Materials & Methods.

The stochastic nature of loop formation provides method to effectively incorporate protein-chromatin and chromatin-chromatin interactions. Looping is often thought to be mediated by DNA-binding factors such as CTCF [33], Sat1B [34] or PcG [10] or by regions of increased polymerase concentration, i.e. transcription factories [24]. The probabilistic creation of functional chromatin contacts mimics the effect of protein concentration (there being either proteins binding DNA sites or not) and binding affinity. The detailed nature of the binding affinity thus does not need to be considered explicitely in our model. In the following we denote by “loop” a functional interaction between two parts of a chromatin fiber existing for a certain time as created by the algorithm. In contrast, a “contact” denotes two parts of the chromatin fiber close together by thermal fluctuations without necessarily being an interaction.

A typical human chromosome has a length of about 100 mega basepairs (Mb), rendering a detailed description on the molecular level computationally impossible. Typically, coarse-grained approaches are used, where a long stretch of chromatin is modeled as an effective monomer. Polymer scaling theory [35] tells us that for linear polymers such an approach is well justified above the scale where bending rigidity plays a role. This length scale is established by the persistence length An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e005.jpg, defining the transition from a rod-like to a flexbile polymer. Estimates for the persistence length range from An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e006.jpg nm [36] but are often based on crude approximations by fitting data to linear chain models [8], [30]. Thus, it is reasonable to conduct computer simulations on a coarse-grained scale where it can be securely assumed that the fiber is flexible.

To study the impact of diffusion-based loop formation on the conformational properties isolated from effects of the presence of other chromosomes, we simulate single chromosomes in a dilute solution. In fact, it has been argued that the disentanglement time for the transition from interphase to metaphase chromosomes of size 100 Mb is in the order of 500 years [37], [38], thus requiring the activity of DNA topoisomerase II. Rosa et al. reversed the argument proposing that interphase chromosomes never equilibrate [38]. We ask whether the observed confined folding already arises from the experimentally confirmed loop formation without invoking rather unprecise knowledge of time and length scales. If loop formation turns out to cause confined folding, then the presence of other chromosomes should not alter the conformational properties drastically. That is why we focus first on isolated chains. In a coarse-grained approach we study chain lengths of size An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e007.jpg and An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e008.jpg. We use Monte-Carlo simulations on a cubic lattice employing the well-established bond fluctuation algorithm [39]. The lattice size is chosen to be An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e009.jpg. By using periodic boundary conditions and keeping track of unfolded coordinates we avoid forcing the polymers into a confined space.

While simulations of diluted chromosomes can be used to study the effect of looping isolated from the presence of other chains, simulations of polymers in a dense system are necessary to study the formation of chromosome territories and to answer the question whether density-related effects are observable. Thus, it is a natural next step to perform simulations in a system with many chromosomes. For our simulations we choose a linear simulation box of width An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e010.jpg and a density of An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e011.jpg, which is similar to the conditions in interphase nuclei. A total of 4096 monomers was thus studied.

Most often, simulational studies map coarse-grained monomers to physical length scales, e.g. by assuming a certain persistence length [31], [38]. Thus, a parameter-dependent comparison between the physical distance of two markers with experimental data from FISH measurements can be conducted. To obtain testable quantitative predictions we follow another, more universal, approach. We derive quantities which are independent on the detailed mapping of the model fiber to the biological chromosome, but can be easily evaluated both from simulational data as well as from experimental data. Such quantities comprise dimensionless higher-order moment ratios of the distance distributions as well as scaling exponents. We show that these quantities do not depend on the chosen level of coarse-graining, i.e. the chain length. Thus, without assuming unknown time and length scales, a sensitive comparison between theory and experiment is possible.

Results

Modeling chromosomes with complex interaction patterns results in the need to dramatic simplifications in order to allow sufficient relaxation of the fiber within a feasible computational effort. Therefore, we study the looping dynamics for isolated coarse-grained chromatin fibers first. Although such conditions are not found in in-vivo experiments, the formation of loops and its influence on the chromosome properties can be studied isolated from density-related effects. In the next step, we present the results of simulations of a system of several chromosomes at biological densities. Since the looping results in confined structures, as will be shown below, density-related effects are only minor and, consequently, the formation of chromosome territories is observed.

Mean square distance between chromatin segments

We first show that the Dynamic Loop model is in agreement with experimental data from FISH measurements [12], [40], which provide information about the relative physical distance between two target sites. The mean squared distance value An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e012.jpg between those target sites in relation to genomic distance An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e013.jpg between them can be compared to polymer models. The random walk (RW) and self-avoiding walk (SAW) polymer models predict this mean squared distance to increase monotonically with the distance between two FISH markers,

equation image
(1)

where An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e015.jpg is a model-dependent parameter [41]. In principle, such a scaling is only valid for the end-to-end distances, however, we want to stress that in the absence of other interactions, equation (1) is approximately valid for genomic separations An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e016.jpg of interest much larger than the persistence length An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e017.jpg of the chromatin fiber.

The confined space of the nucleus renders the random walk and self-avoiding walk polymer model inadequate. A 100 Mb chromosome with assumed Kuhn length of approximately 300 nm [8] in a 30 nm fiber packing (300 nm An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e018.jpg 30 kb) would extend on average to 17.3 An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e019.jpgm in the random walk case, whereas the average diameter of a nucleus is of the order of 10 An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e020.jpgm. The globular state model fails for other reasons [32]. Recent experiments [12], however, clearly revealed that while the mean square distance increases monotonically with genomic separation on short distances up to a few Mb, a leveling-off is observed for larger genomic separations. This confined folding is observed on a scale of about 2 An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e021.jpgm, far below the diameter of the nucleus but consistent with the estimated size of chromosome territories [7]. The random loop model [12], [42] explains the behavior by the formation of random loops, without invoking a confined geometry a priori.

We first considered the mean square distance between two beads in the DL model for isolated chains. Given a chain of length An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e022.jpg with monomer positions denoted by An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e023.jpg, the average is calculated over a set An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e024.jpg of independent conformations as well as over different reference points inside the chain

equation image
(2)

Fig. 1 shows the results of the model for a chain of length An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e026.jpg. The looping probability An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e027.jpg is varied such that different values for the average number of loops are obtained. Lacking knowledge of the biological lifetime of the loops, results are shown for three different values of An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e028.jpg depending on the relaxation time of the chromosomes An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e029.jpg (triangles [[filled triangle]] for An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e030.jpg, open diamonds [◊] for An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e031.jpg and filled circles [An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e032.jpg] for An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e033.jpg, see Materials&Methods). The model displays a cross-over from self-avoiding walk behavior (small number of loops) to a leveling-off in the mean square distance. Such a plateau level is recovered if the average number of loops on a coarse-grained chromosome is larger than about An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e034.jpg. This result is independent of the lifetime of the loops as long as the average number of loops remains the same, indicating that the lifetime has no direct influence on the statistical equilibrium properties. These findings clearly show that no long-range interactions are necessary for forcing the polymer to collapse but a purely diffusional motion together with chromatin-chromatin binding affinity suffices to achieve this.

Figure 1
Mean square distance An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e035.jpg in relation to contour length for an isolated fiber.

To quantitatively compare the model to experimental data, we assume each bead to represent a 400 kb segment of chromatin with an average extension of 480 nm (in agreement with experimental data [12]). To ensure that the qualitative results are not dependent on chain length, we studied the mean square distance for An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e044.jpg and An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e045.jpg (see Figure S1). In all three cases a leveling-off is observed, indicating that the observed results are independent on scale and the applied coarse-graining is justified.

Self-organized formation of large loops

Loop formation is a central process for the transcriptional regulation in higher eukaryotes. Several studies indicated that co-localization of chromatin segments results in activation or repression of genes [20]. Hypotheses of loop formation range from the attachment to a structure called nuclear matrix [43] to the formation of transcription factories [24], in which transcriptionally active genes come together, forcing the intervening DNA to loop out. It has been proposed that rossette-like loops arise in a self-organized manner due to the heterogeneity of the fiber [44]. Recently it has been shown [28], [29] that loops can promote territory formation with a simple model using fixed loops. However, such a kind of looping does not yield a correct description for the relative positioning of two markers [12]. Rather, it has been shown by 4C experiments [25], that loops exist on scales up to several Mb. 3C/4C/5C and the newly developed Hi-C [45] techniques provide an experimental method to measure loop probabilities and distributions. Therefore, we next investigated how the model alters the distribution and frequency of genes to become co-located. Again, we favor measures that do not depend on the level of coarse-graining and parameters like persistence length. One such measure is the decay of the contact probability and abundance with genomic separation An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e046.jpg. Consider a random walk polymer chain. Clearly, the probability that two beads An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e047.jpg and An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e048.jpg come into contact decreases with the separation An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e049.jpg. More precisely, we obtain a power-law behavior [35]

equation image
(3)

Consider two genes separated by 10 Mb. Assuming a Kuhn segment length of 300 nm [38] consisting of 30 kb chromatin, the probability of co-localization is in the order of An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e051.jpg. How, then, does the cell nucleus manage to co-locate different chromatin segments in a reasonable time? To answer this question, we look at the formation of functional loops in our model and its size distribution An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e052.jpg. Interestingly, the diffusional pathway to loop formation results in a size distribution of functional loops An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e053.jpg which is quite different from the small random contact probabilities of a RW or SAW model (Fig. 2C). Strikingly the probability of having a loop in the size-range of the chain length is enhanced by over two orders of magnitude. The increase in probability for large-scale loops in contrast to small-scale loops can be explained on an intuitive basis: Starting from a linear chain, the diffusional process will bring monomers close together which are not so far away along the contour of the chain. Loop formation will be dominated by small-sized loops as equation (3) still holds. However, as more and more small loops form, even parts of the polymer located further apart come closer together (Fig. 2A), thus enhancing the probability of contact. Figure 2B visualizes for one simulation run the average loop size along the simulation time. We find that this average loop size increases fast and then fluctuates around an equilibrium value. Therefore diffusional looping seems to be a quite fast and effective method of large loop formation.

Figure 2
Size distribution of loops and random contacts.

To allow a comparison to experimental data from 4C and 5C experiments we determine two measures. Firstly, the size-distribution An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e069.jpg of random contacts (as An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e070.jpg C experiments do not only measure functional contacts) between two chromatin segments. Secondly, the specific contact probability An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e071.jpg that two segments at position An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e072.jpg and An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e073.jpg are in contact. From eq. (3) we know that for a random walk the specific contact probability has a power-law behavior depending on the length An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e074.jpg given by An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e075.jpg. A power-law behavior is also found for the self-avoiding walk, where the exponent is determined in Fig. 2D to An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e076.jpg. In fact, scaling theory predicts [35] for self-avoiding or random walk polymers that the contact probability of the end-points of a polymer scales as An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e077.jpg. Our analysis suggests that the contact probability for intra-chain segments decreases more strongly. This is somehow expected, as intra-chain segments have less entropic degrees of freedom and are surrounded by a higher density of adjacent beads than the end points, making contacts with beads further away less likely. Our polymer model as well displays a power-law behavior of the co-localization probability An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e078.jpg (Fig. 2D). However, two different regimes have to be distinguished. For genomic separations An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e079.jpg in the size range of the whole chromosome a different exponent is found as in the size range below about 15% of the fiber length. In the regime of probability-values An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e080.jpg where leveling-off in the mean square distance is observed, we find exponents of about An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e081.jpg for smaller genomic separations in the order of 10 Mb and An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e082.jpg for large genomic separations in the order of 100 Mb (table 1). Intriguingly, the probability of specific contacts between far-apart chromatin segments is increased by over two orders of magnitude compared to the self-avoiding walk. Increasing the looping probability and thus the average number of loops per chain results in smaller exponents An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e083.jpg. Interestingly, this result is in close agreement with recent results from Hi-C data [45] where an exponent of An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e084.jpg has been observed in a region between 500 kb and 7 Mb (see inset of Figure 2D). For genomic separations above 10 Mb we find a scaling exponent of An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e085.jpg, consistent with the smaller scaling exponent found in our model. Similar results are found for other chain lengths (Figure S2 and Figure S3).

Table 1
Decay exponents of the random contact probabilities with genomic separation for direct comparison to 4C and 5C experiments.

5C data provides a detailed map of interactions between chromatin segments without a fixed reference point. Thus, it is more natural to look at the relative abundance of contacts An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e111.jpg of size An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e112.jpg, encompassing all fragments of a certain length An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e113.jpg found in the data independent on their position on the genome. A crude power-law fit An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e114.jpg can be conducted here, too (see Figure S4). We find power-law exponents of An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e115.jpg in the range where leveling-off in the mean square distance is observed (cf. Fig. 1). The exponents both for the specific contact probability An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e116.jpg as well as the size-distribution An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e117.jpg are listed in table 1.

Fig. 3 shows contact maps similar to those obtained by 5C for a An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e118.jpg polymer with different looping probabilities. Contacts between any two beads are marked by a black square. For better visibility, in each map contacts of 4 equilibrated conformations are plotted. Clearly, the self-avoiding walk polymer model (Fig. 3A) only has a few contacts between beads located far apart. Increasing the looping probability (Figs 3B and C) results in a strong increase of both the number of loops as well as the abundance of large loops.

Figure 3
Intra-chromosomal contacts of isolated model polymers.

Theoretically, the change in the scaling exponent can be explained by the topological changes induced to the fiber on introducing loops. For a polymer network, the looping probability generally behaves like An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e120.jpg with the exponent An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e121.jpg depending on the specific topology, i.e. the vertices, of the network [46]. For more complex networks, i.e. loop configurations, the exponent decreases due to the less conformational degrees of freedom available.

Cell-to-cell variation and dynamic fluctuations of the distance distribution

While FISH measurements have been used to establish a connection between the mean square distance of two markers and genomic separation [8], [12], [40], a direct comparison to polymer models requires parameters to map one model bead to physical units like nanometers and base pairs. As these parameters are unknown or based on crude estimates [8], [30], it is desirable to introduce dimensionless quantities not dependent on length scale parameters.

For the the random walk (RW), self-avoiding walk (SAW) and the globular state (GS) model, the following higher-order moments of the distance distribution between two markers turned out to be basically independent of genomic separation [32].

equation image
(4)

An intrinsic advantage of these measures is that they are dimensionless, i.e. both experiments and models yield a numeric value. Even more important, the ratios carry information about the fluctuations, i.e. the cell-to-cell variation of the measurements.

One prominent feature of FISH measurements in interphase chromatin is that the fluctuations of the distance distributions are larger than expected from a random walk or self-avoiding walk polymer model [42]. Recently it has been shown that this holds true for the case of compact polymers as well [32], where the fluctuations are even smaller. The ratios given in eq. (4) for experimental data sets from Mateos-Langerak et al. [12] as well as Jhunjhunwala et al [27] are presented in Fig. 4A. The figure contains FISH data from human chromosomes 1 and 11 [12], separately measured for ridges (green squares) and anti-ridges (red squares) as well as data from the murine Igh locus [27], which was kindly provided by K. Murre.

Figure 4
Higher-order moments of the distance distributions for experimental data (A) and for the chromatin model (B) according to eq. (4).

The results for the model treated in this paper are shown in Fig. 4B. Model polymers of different length (An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e129.jpg and An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e130.jpg) have been equilibrated and averaged over a huge ensemble of conformations encompassing various configurations of loop attachment points. The data is plotted against the average number of loops per monomer to allow for a comparison between different chain lengths. For small looping probabilities, i.e. small average number of loops, the self-avoiding walk behavior is recovered, whereas increasing the looping probability leads to a strong increase in the fluctuations of the system. The higher-order moment ratios markedly exceed the random walk value in the range of loop numbers between An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e131.jpg and An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e132.jpg. One would expect for a random-walk polymer to have larger fluctuations than a polymer constraint to excluded volume interactions and topological constraints. In our model, the large fluctuations are induced by the dynamic formation of loops, which thus seems to be an important characteristics of chromatin organization. However, it has to be noted that the fluctuations of the model are still too small to explain the moment ratios of most of the experimental data. We will discuss this in more detail in the discussion section.

The dynamics of looping chromosomes

Finally we study the dynamics of the looping chromatin fibers. The center-of-mass motion of a polymer is measured by An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e133.jpg. For a self-avoiding walk polymer it shows normal diffusion behavior, i.e. An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e134.jpg. As can be seen in Fig. 5 the chromatin model shows subdiffusive motion An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e135.jpg on time scales smaller than the relaxation time of the polymer. The actual diffusion exponent An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e136.jpg depends on the looping probability An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e137.jpg. For times larger than the relaxation time one recovers diffusive motion, i.e. An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e138.jpg, however, this motion is slower than for a normal self-avoiding walk (Figure S5). This is consistent with experimental results showing that chromosome territories do not move significantly [47].

Figure 5
Dynamics of the center of mass and the central monomers.

It has to be noted that the regime of large times is not very sensitive for a comparison to experimental data as here the confinement by other polymers comes into play which is not incorporated into the simulations of a single polymer. It is more instructive to look at the motion of the central monomers of a chain on short time scales. The mean square displacement An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e146.jpg displays a distinct behavior for three different time regimes, which are related to the relaxation time An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e147.jpg of a chromosome. For An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e148.jpg there is a pronounced subdiffusive behavior. The anomalous diffusion exponents range from An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e149.jpg in the regime where leveling-off is observed for the mean square displacement (cf. Fig. 1). For An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e150.jpg the predictions of classical polymer dynamics become valid again and we find An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e151.jpg similar to the self-avoiding walk. On large time scales (An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e152.jpg), the monomer motion follows the motion of the center of mass, displaying normal Brownian motion with An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e153.jpg.

While at intermediate and large time scales the motion can be described by classical polymer theory, i.e. Rouse dynamics [41], the scaling exponents on the short time scale An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e154.jpg are unexpected. Following the argument in Refs. [37], [38], this time scale is the prevailing one concerning interphase chromosomes. Clearly, such exponents arise due to the constraints induced by looping, which temporarily slows down the motion of chromatin segments at the loop attachment points. Although experimental data is rare, this is consistent with findings of Cabal et al. [17] in yeast. This study showed that the motion of a labelled spot scales like An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e155.jpg up to An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e156.jpg. Interestingly they found the exponent to depend on the transcriptional state of the GAL genes. This is in support of our conjecture put forward in another publication [12] that the local looping probability may be related to transcriptional activity.

The formation of aspherical chromosome territories

Polymer theory predicts that equilibrated polymers with a large molecular weight in a semi-dilute solution are strongly intermingling [35]. Various studies, however, indicate that chromosomes occupy discrete functional domains [7], [48], [49]. It was shown above that our model polymers adopt a confined structure by virtue of dynamic looping. Amazingly, this result was obtained without subjecting the system to a confined space (in contrast to Refs. [30], [31]) and without introducing long-range interactions (in contrast the polymer models in Refs. [11], [42], [50]).

Surely, simulating isolated chromatin fibers does not yield complete information about the folding in a dense system as in the nucleus, e.g. the formation of chromosome territories (CTs). To investigate whether probabilistic loops are the reason for the formation of chromosome territories, we set up simulations of chromosomes in a box of width An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e157.jpg lattice units and a length of An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e158.jpg. The density of the system An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e159.jpg was chosen close to the estimates of chromatin in cell nuclei. Similar values were used in other publications [28].

An established measure of territory formation is the number of contacts displayed in the contact map [28]. Figure 6 shows such contact maps for different looping probabilities An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e160.jpg. Each map displays contacts between the beads of a subset of 10 model chromosomes out of the system. The beads are numbered consecutively, i.e. bead 0 to 127 belong to chain 1, bead 128 to 256 belong to chain 2 etc. Subsequent chains are alternatingly marked by black and white bars. We find that linear self-avoiding walks (Fig. 6A) without loops display a relatively large number of inter-chromosomal in comparison to intra-chromosomal contacts: 15.4% of the contacts are found to be with other chromosomes. With increasing looping probability An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e161.jpg, the percentage of contacts between different chromosomes decreases. Fig. 6B displays chains with an average number of 45 loops (blue triangle in Fig. 7A). Here we find that only 1.8% of the contacts are inter-chromosomal and in Fig. 6C (92 loops per chain on average, green triangle in Fig. 7A) this value reduces further to An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e162.jpg. Thus, the level of intermingling between CT's strongly depends on the local looping probabilities. As different local looping probabilities seem to play a dominant role in chromatin organization, this finding could explain different levels of intermingling found in several studies [51], [52]. Branco et al. [51] found out that about 20% of the chromosomes are in contact with other chromosomes. To be able to compare these experiments with our model, we determined the fraction of monomers of one chromosome which are close to neighbouring chromosomes. For linear chains without loops (Figure 6A), we find an overlap fraction of about 50%, i.e. the polymers intermingle strongly. For chromosomes with 45 loops per chain on average (Figure 6B), the overlap fraction is about 15%, i.e. slightly smaller than the average observed for chromosomes. Branco et al. indeed observed different overlap fractions dependent on chromosome activity, thus different local levels of looping might indeed mediate such different overlap fractions.

Figure 6
Contact maps and illustrations of chromosomes with different looping probabilities.
Figure 7
Properties of looping polymers in a dense system.

Thus, a disentangling of the fibers, which has been estimated to require a huge amount of time or the action of topoisomerase II [37] is not necessary. Rather, loop formation alone induces a strong repulsive interaction between different chromosomes; a finding which has been quantified for ring polymers [53] and rosette structures [54]. In mitosis, chromatin adopts a compact state, where different chromosomes are unentangled and well-separated. At the onset of interphase, the loop formation forces the chromosomes to a more open, but confined structure, which results in the formation of CTs without requiring the assumption of unequilibrated polymers [38].

We find that the predictions from the study of isolated model chromosomes are still valid for a dense system of chromatin. Amongst others, this is a direct consequence of loop-based segregation observed in Fig. 6. Fig. 7A shows the mean square distance between two model beads in relation to contour length An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e180.jpg (in biological terms: genomic distance An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e181.jpg). Similar to the results of Fig. 1, the mean square distance displays a leveling-off for average loop numbers larger than about 45 loops per coarse-grained monomer. Obviously, for small looping probabilities An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e182.jpg (red curve) or self-avoiding walks (An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e183.jpg, not shown), polymers do not level off, thus they do not form separate territories. The behavior of territory formation and segregation is a distinct result of loop formation.

While the mean square distance An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e184.jpg displays a leveling-off for several polymer models (e.g. globular state [31], [32], random walk in a confined space [30], etc.), a more sensitive measure are again the dimensionless ratios of higher order moments given by eq. (4). As the fluctuation regime could possibly change under the transition from isolated polymers to a dense system, we investigate the ratio An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e185.jpg. Fig. 7B shows that fluctuations are larger than predicted by the random walk, self-avoiding walk or globular state model. In the regime where a leveling-off is obtained in Fig. 7A, i.e. the average loop number is larger than about 45, the moment ratios are approximately in the range An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e186.jpg.

The relative abundance of contacts An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e187.jpg is displayed in Fig. 7C for polymers in a dense system. Again, the co-localization frequency is greatly enhanced by the formation of functional loops. A crude power-law fit An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e188.jpg results in exponents of An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e189.jpg and smaller in the region where a leveling-off is observed in the mean square distance. Similar results are found for the specific contact probability An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e190.jpg (Fig. 7C), which display a biphasic behaviour already observed in the case of isolated chromosomes (Fig. 2). In the size range of large genomic separations in the order of the entire chromosome, the contact probability decreases with a power-law An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e191.jpg with exponents starting from An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e192.jpg in the self-avoiding walk model to An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e193.jpg in the parameter range where leveling-off is observed. On intermediate scales (An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e194.jpg Mb), for biologically relevant looping probabilities, an exponent of An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e195.jpg is found. Amazingly, a similar value of An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e196.jpg has been recently found by Hi-C experiments [45] on a scale between 500 kb to 7 Mb.

Isolated model chromosomes displayed a pronounced conformational asphericity (Figure S6). A similar behavior is also found for chromosomes in a dense system. In fact, deviations from a sphere-like shape are expected for the self-avoiding walk as well as the random walk model [55], however, not for a compact globular state polymer [32]. Whereas looping polymers can adopt a highly compacted state, their properties differ clearly from a globular state. Indeed, the shape of simulated chromosomes territories is not spherical as one would expect for compact polymers, rather we find that the gyration ellipsoid has a prolate shape. The ratios of the gyration tensor's eigenvalues are listed in table 2. In the parameter range where the mean square distance displays a leveling off, we find ratios of the eigenvalues An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e197.jpg in the regime between An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e198.jpg and An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e199.jpg. These values are smaller than those of Rosa et al. [38] for ring polymers and consistent with those of looping polymers [28]. A non-spherical shape of CTs has also been found in experimental studies [56], [57]. Mouse chromosomes exhibit an aspherical shape approximated by ellipsoids with axis ratios An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e200.jpg. A one-to-one correspondence of these numbers with results from the shape of the gyration ellipsoid, however, can not be established.

Table 2
Shape parameters of simulated chromosomes.

Discussion

In this study, a polymer model was presented where loops form dynamically on the basis of diffusional collisions. We use Monte Carlo simulations to demonstrate the effect of such a kind of loop formation. While loops have been recognized as an ubiquitous feature in transcriptional regulation, the pathways of its formation remained unclear and most polymer models proposing loops so far did not explain the transport mechanisms by which two parts of chromatin become co-located. Our results suggest that even large loops can arise without active transport mechanisms. Our model neither assumes a confined geometry nor any long-range interactions. Loop formation is based on the diffusional motion of the fiber. Collisions lead to a probabilistic chromatin-chromatin interaction which forces the participating regions to be co-located for a certain time. The probabilistic nature of the interactions is meant to mimic the effect of chromatin binding factors on chromatin-chromatin interactions. Although this Dynamic Loop model is kept minimal, it reproduces many experimental results quantitatively, highlighting the possibility that chromatin folding is tightly related to function through the loop formation process.

One of our major results is that dynamic loop formation drives chromosomes into an entropically segregated state. Indeed, linear polymers intermingle freely (Fig. 6A) in agreement with polymer theory [35]. Looping polymers, in contrast, fold into a confined space (Fig. 1). Such a confinement is also observed for the globular state polymer model, which, however, displays a markedly different fluctuation regime than the experimental data [32]. The importance of looping on the formation of chromosome territories has been investigated recently by Cook et al. [28]. In their qualitative study, rosette-structures with fixed loop attachment points are used. While this model can be used as a simple model for studying entropic effects of looping, it does not explain fluctuations in FISH data [12], [42].

The second important result of this study concerns the pathway of the formation of large loops. Results from 3C/4C/5C experiments reveal that loops are abundant on the short scale [58]. Nevertheless, functional loops on the scale of several mega basepairs have been detected in 4C experiments [25]. While the probability of specific random contacts An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e211.jpg decreases strongly with site separation An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e212.jpg for linear polymers (given by a power-law behavior An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e213.jpg), the contact probability is increased by over two orders of magnitude when introducing loops (Fig. 2). Obviously, small functional loops which can easily co-localize by diffusional motion strongly support the formation of long-distance contacts. The contact distribution An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e214.jpg displays power-law exponents of An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e215.jpg, the contact probability An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e216.jpg exponents in the size range between An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e217.jpg. This is in agreement with recent experimental data by Lieberman-Aiden et al. [45], however, their interpretation in terms of a fractional globule differs from ours. Clearly, a fractional globule, where the physical distance between two loci scales with An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e218.jpg [59], is in contrast to experimental findings from FISH data [12].

The impossibility to perform Monte-Carlo simulations on a very detailed scale requires a coarse-graining procedure. Looking at large-scale features above the persistence length An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e219.jpg, such an approach is well-justified [35]. For linear polymers, scaling laws provide a simple way of rescaling a polymer. For a model with loops, the connection between chain length, bond length and looping probability An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e220.jpg is non-trivial. In fact, even for chromatin models using linear chains (see Refs. [30], [31], [38]), the establishment of a correspondence between simulational units and biological units requires the knowledge of the persistence length of chromatin. The latter has been estimated by fitting a random walk model [30] or a worm-like chain model [8] to FISH data. Estimates on the persistence length range from 40–220 nm [36]. For a quantitative comparison to experimental data we derived measures independent on both the level of coarsening and unknown biological parameters. These can be easily evaluated both for experimental data as well as polymer models. Amongst others, these measures comprise the power-law exponent of the contact distribution (Fig. 2 and table 1), the dimensionless higher-order moment ratios of the distance distribution between two FISH markers (Fig. 4), the asphericity of chromosomes and finally the diffusion exponents (Fig. 5).

The DL model studied here displays a pronounced aspherical elongated shape (table 2) which is also found in experiments [6]. Consistent with experimental data in yeast, the motion of single monomers is subdiffusive (Fig. 5); the actual subdiffusion exponent depends on the looping probability, which was suggested to be closely related to transcriptional activity [12]. A good agreement with experimental data is obtained for the higher-order moments of the distance distribution. Surprisingly, the moment ratios (given in eq. (4)), which display fluctuations of the distance distributions, exceed the random walk value for looping probabilities An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e221.jpg in the range where a leveling off in the mean square distance is observed. This is not necessarily expected for a model with excluded volume which restricts the degrees of freedom and therefore shows less fluctuations. The increase of the fluctuations with respect to the self-avoiding walk is due to the dynamic formation of loops. However, several independent experiments [12], [27] consistently show even larger fluctuations. We suspect two major reasons for this: First, the chromatin fiber is not a homogeneous polymer and there is evidence that looping probabilities vary depending on the transcriptional state [12]. Secondly, inside the nucleus, topoisomerase-II might effectively counteract excluded volume interactions, resulting in an underestimate of the fluctuations in our model.

While our model suggests that chromosome segregation might be driven by the diffusional formation of loops, Rosa and Everaers suggested [38] that segregation is a consequence of large entanglement times. The entanglement times, however, might be strongly reduced by the effect of topoisomerase-II [60]. Notwithstanding that time-scales play an important role, this study reveals that loop formation provides a complementary and fully sufficient mechanism for CT formation.

Importantly, the Dynamic Loop model displays a distinctly different behaviour than the equilibrium or fractal globular state (GS) models, which have been proposed for chromatin organization [31], [45]. First, the GS model does not explain the large fluctuations observed in several experiments (Figure 4), rather it predicts fluctuations being even smaller than for a random walk or self-avoiding walk model. Second, the equilibrium globular state predicts contact probabilities to behave as An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e222.jpg [61], while Hi-C data shows An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e223.jpg, consistent with our model. On the other hand, the fractal globular state predicts the correct contact probabilities but not the leveling-off in the mean square displacement. Third, the globular state predicts rather spherical territories, while our model as well as experimental data shows aspherical chromosome territories. The differences indeed arise by the introduction of probabilistic loops which among others leads to the large fluctuations [42] observed in experiments.

Clearly the chromatin model proposed here does not capture all details of the complex nuclear organization. For reasons of simplicity we neglect the heterogeneity of the chromatin fiber here and assume the same looping probability and chromatin affinity along the complete chromosome. Using such a simple model allowed us to derive the basic properties of chromatin fibers without introducing unnecessarily many parameters. However, experiments clearly show that loop formation is strongly dependent on the differentation state of the cell [25] as well as gene activity [12]. In future work such differences might be incorporated into our model either by adjusting the looping probabilities locally or by distributing specific binding sites along the chromosomes.

Materials and Methods

The biological model is implemented using Monte Carlo simulations [62]. These Monte Carlo simulations are performed on a lattice in order to simplify the handling of excluded volume. Calculation of excluded volume interactions thereby is reduced to checking whether one lattice site is already occupied or not. Instead of using a simple local-move algorithm on a cubic lattice we employ the bond-fluctuation method introduced by Carmesin [39]. This method has the advantage over other lattice models of allowing 108 different bond vectors; the length of a bond can take the values An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e224.jpg [63]. The bond-fluctuation model is especially suited for dense and compact systems where a lattice algorithm would no longer be feasible due to high rejection rates during the Monte Carlo process. It has successfully been applied to several studies on the static and dynamical properties of polymer systems [63][66]. The simulation method fulfills the following important features: (i) it produces unbiased results, i.e. each possible conformation out of the ensemble is sampled with equal probability, (ii) it takes into account excluded volume interactions, i.e. two monomers are not allowed to occupy the same region in space and (iii) using some restrictions on the moves and bond vectors it ensures that no bond crossings can occur during a Monte Carlo step, i.e. it preserves the topological state of the conformation. The algorithm conducts only local moves in order to resemble the dynamics of real polymers [39]. Using a coarse-grained lattice approach is reasonable as we are only interested in features of looping chromatin independent on local structure. Coarse-graining allows us to abstract from the complex environment and highlight the main driving forces and effects of chromatin folding.

Simulations for single polymers are performed on a lattice of size An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e225.jpg. Periodic boundary conditions are used, but the algorithm always keeps track of unfolded coordinates, such that the polymer does not feel any confined volume. The lattice size An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e226.jpg is chosen larger than the radius of gyration An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e227.jpg of the chains studied such that effects of the backfolding are negligible.

A dense system of model chromosomes is simulated in a system of size An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e228.jpg and chain length An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e229.jpg. The total number of monomers is 4096, the density An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e230.jpg.

In order to obtain thermodynamical equilibrated conformations we perform the Metropolis Monte Carlo method. Chromosomes are initially equilibrated as self-avoiding walks using local moves of a monomer to one of the nearest neighbors on the lattice. After the initial equilibration steps, the Monte Carlo algorithm allows for the formation of loops. After each Monte Carlo trial move, one monomer is selected at random. It is then checked whether another monomer on the same chain is in the neighborhood, i.e. co-localized. The co-localization condition is fulfilled whenever the distance between the monomers is less than 3 lattice units. If the two monomers are co-localized, then a loop is formed with a certain probability An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e231.jpg. If the loop An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e232.jpg is created, it is assigned a certain lifetime An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e233.jpg which is drawn from a Poissonian distribution

equation image
(5)

where the parameter An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e235.jpg determines the average lifetime of the loops. In the simulations we use three different values of An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e236.jpg:

equation image
(6)

where An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e238.jpg is the integrated autocorrelation time (see below) of the squared radius of gyration for the corresponding self-avoiding walk system. Loop lifetimes are chosen relative to the relaxation time to make results for different values of the other parameters comparable. As an example, we want to give a few numbers here. For a single chain of length An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e239.jpg, the equilibration time is An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e240.jpg MCS; thus, one corresponding loop lifetime would be An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e241.jpg MCS. A looping probability of An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e242.jpg yields an average number of 19 loops per chain, a looping probability of An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e243.jpg yields an average number of 131 loops per chain.

Since subsequent conformations in the Markov chain created by the Monte Carlo algorithm are highly correlated, one has to perform a certain number of Monte Carlo steps to obtain two independent conformations. For each set of parameters (chain length An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e244.jpg, looping probability An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e245.jpg and lifetime of loops An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e246.jpg) we determine the autocorrelation function An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e247.jpg (see e.g. Ref. [67]) of the squared radius of gyration An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e248.jpg

equation image
equation image

To obtain a reasonable result for the integrated autocorrelation time we have to sample about An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e251.jpg Monte Carlo steps. Then An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e252.jpg itself is approximated by the following algorithm, which is often called the windowing procedure and was introduced by Sokal [68].

  1. Determine An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e253.jpg.
  2. Choose An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e254.jpg to be the smallest integer such that An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e255.jpg. Then An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e256.jpg is supposed to be the best approximation of the autocorrelation time.

Another criterion to ensure the uncorrelatedness of subsequent conformations is given by the motion of the center of mass. This method has been used for example by Mueller et al. in his study on ring polymers [65]. Here we determine the function

equation image

The time of interest, An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e258.jpg, after which the center of mass has moved at about one radius of gyration, is defined by

equation image
(7)

We consider two subsequent conformations as uncorrelated after An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e260.jpg Monte Carlo steps. Actually, for each set of parameters considered here, we found that after this time the center of mass has moved on average at least by one radius of gyration, i.e. An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e261.jpg.

Simulations of even small polymers are very time-consuming due to the looping interaction and the resulting compactness of the polymers. Furthermore, simulation runs have to be quite long to capture the dynamics of loop formation. We have used the Helics2- and bwGrid parallel computing facilities at the Interdisciplinary Center for Scientific Research (IWR) at the University of Heidelberg. For each set of parameters An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e262.jpg we created 10 000–100 000 independent conformations. We study polymers of lengths An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e263.jpg and An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e264.jpg. The looping probabilities are chosen such that the average number of loops in the resulting conformational ensemble is between zero and An external file that holds a picture, illustration, etc.
Object name is pone.0012218.e265.jpg. The lifetimes of the loops are chosen from the set given in eq. (6).

Supporting Information

Figure S1

Mean square distance for various N.

(0.19 MB PDF)

Figure S2

Loop size distribution and specific contact probabilities for N  = 128.

(0.11 MB PDF)

Figure S3

Loop size distribution and specific contact probabilities for N  = 512.

(0.12 MB PDF)

Figure S4

Relative abundance of contacts h(l) in relation to genomic separation l.

(0.15 MB PDF)

Figure S5

Dynamics of the center of mass and motion of the central monomers.

(0.12 MB PDF)

Figure S6

The shape of chromosomes.

(1.41 MB PDF)

Footnotes

Competing Interests: The authors have declared that no competing interests exist.

Funding: M.B. has received funding from the Landesgraduiertenforderung Baden-Wurttemberg. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

References

1. Schiessel H, Gelbart WM, Bruinsma R. DNA folding: structural and mechanical properties of the two-angle model for chromatin. Biophys J. 2001;80:1940–1956. [PubMed]
2. Woodcock C, Frado L, Rattner J. The higher-order structure of chromatin: evidence for a helical ribbon arrangement. J Cell Biol. 1984;99:42–52. [PMC free article] [PubMed]
3. Horowitz RA, Agard DA, Sedat JW, Woodcock CL. The three-dimensional architecture of chromatin in situ: electron tomography reveals fibers composed of a continuously variable zig-zag nucleosomal ribbon. J Cell Biol. 1994;125:1–10. [PMC free article] [PubMed]
4. van Holde K, Zlatanova J. What determines the folding of the chromatin fiber? Proc Natl Acad Sci U S A. 1996;93:10548–10555. [PubMed]
5. Diesinger PM, Heermann DW. Depletion effects massively change chromatin properties and influence genome folding. Biophys J. 2009;97:2146–2153. [PubMed]
6. Goetze S, Mateos-Langerak J, Gierman HJ, de Leeuw W, Giromus O, et al. The three-dimensional structure of human interphase chromosomes is related to the transcriptome map. Mol Cell Biol. 2007;27:4475–4487. [PMC free article] [PubMed]
7. Cremer T, Cremer C. Chromosome territories, nuclear architecture and gene regulation in mammalian cells. Nat Rev Genet. 2001;2:292–301. [PubMed]
8. Bystricky K, Heun P, Gehlen L, Langowski J, Gasser SM. Long-range compaction and flexibility of interphase chromatin in budding yeast analyzed by high-resolution imaging techniques. Proc Natl Acad Sci U S A. 2004;101:16495–16500. [PubMed]
9. Lowenstein MG, Goddard TD, Sedat JW. Long-range interphase chromosome organization in drosophila: A study using color barcoded fluorescence in situ hybridization and structural clustering analysis. Mol Biol Cell. 2004;15:5678–5692. [PMC free article] [PubMed]
10. Bantignies F, Grimaud C, Lavrov S, Gabut M, Cavalli G. Inheritance of Polycomb-dependent chromosomal interactions in Drosophila. Genes & Development. 2003;17:2406–2420. [PubMed]
11. Sachs R, Engh G, Trask B, Yokota H, Hearst J. A Random-Walk/Giant-Loop Model for Interphase Chromosomes. Proc Natl Acad Sci U S A. 1995;92:2710–2714. [PubMed]
12. Mateos-Langerak J, Bohn M, de Leeuw W, Giromus O, Manders EMM, et al. Spatially confined folding of chromatin in the interphase nucleus. Proceedings of the National Academy of Sciences. 2009;106:3812–3817. [PubMed]
13. Versteeg R, van Schaik BDC, van Batenburg MF, Roos M, Monajemi R, et al. The human transcriptome map reveals extremes in gene density, intron length, GC content, and repeat pattern for domains of highly and weakly expressed genes. Genome Res. 2003;13:1998–2004. [PubMed]
14. Scheuermann MO, Tajbakhsh J, Kurz A, Saracoglu K, Eils R, et al. Topology of genes and nontranscribed sequences in human interphase nuclei. Exp Cell Res. 2004;301:266–279. [PubMed]
15. Dietzel S, Zolghadr K, Hepperger C, Belmont AS. Differential large-scale chromatin compaction and intranuclear positioning of transcribed versus non-transcribed transgene arrays containing beta-globin regulatory sequences. J Cell Sci. 2004;117:4603–4614. [PubMed]
16. Janicki SM, Tsukamoto T, Salghetti SE, Tansey WP, Sachidanandam R, et al. From silencing to gene expression: Real-time analysis in single cells. Cell. 2004;116:683–698. [PubMed]
17. Cabal GG, Genovesio A, Rodriguez-Navarro S, Zimmer C, Gadal O, et al. SAGA interacting factors confine sub-diffusion of transcribed genes to the nuclear envelope. Nature. 2006;441:770–773. [PubMed]
18. Zink D, Amaral MD, Englmann A, Lang S, Clarke LA, et al. Transcription-dependent spatial arrangements of CFTR and adjacent genes in human cell nuclei. J Cell Biol. 2004;166:815–825. [PMC free article] [PubMed]
19. Meaburn KJ, Misteli T. Cell biology: Chromosome territories. Nature. 2007;445:379–381. [PubMed]
20. Fraser P. Transcriptional control thrown for a loop. Current Opinion in Genetics & Development. 2006;16:490–495. [PubMed]
21. Tolhuis B, Palstra RJ, Splinter E, Grosveld F, de Laat W. Looping and interaction between hypersensitive sites in the active β-globin locus. Mol Cell. 2002;10:1453–1465. [PubMed]
22. Dekker J. A closer look at long-range chromosomal interactions. Trends Biochem Sci. 2003;28:277–280. [PubMed]
23. Cook PR. Predicting three-dimensional genome structure from transcriptional activity. Nat Genet. 2002;32:347–352. [PubMed]
24. Marenduzzo D, Faro-Trindade I, Cook PR. What are the molecular ties that maintain genomic loops? Trends Genet. 2007;23:126–133. [PubMed]
25. Simonis M, Klous P, Splinter E, Moshkin Y, Willemsen R, et al. Nuclear organization of active and inactive chromatin domains uncovered by chromosome conformation capture-on-chip (4C). Nat Genet. 2006;38:1348–1354. [PubMed]
26. Münkel C, Langowski J. Chromosome structure predicted by a polymer model. Phys Rev E. 1998;57:5888–5896.
27. Jhunjhunwala S, van Zelm MC, Peak MM, Cutchin S, Riblet R, et al. The 3D structure of the immunoglobulin heavy-chain locus: implications for long-range genomic interactions. Cell. 2008;133:265–279. [PMC free article] [PubMed]
28. Cook PR, Marenduzzo D. Entropic organization of interphase chromosomes. J Cell Biol. 2009;186:825–834. [PMC free article] [PubMed]
29. de Nooijer S, Wellink J, Mulder B, Bisseling T. Non-specific interactions are sufficient to explain the position of heterochromatic chromocenters and nucleoli in interphase nuclei. Nucl Acids Res. 2009;37:3558–3568. [PMC free article] [PubMed]
30. Hahnfeldt P, Hearst JE, Brenner DJ, Sachs RK, Hlatky LR. Polymer models for interphase chromosomes. Proc Natl Acad Sci U S A. 1993;90:7854–7858. [PubMed]
31. Emanuel M, Radja NH, Henriksson A, Schiessel H. The physics behind the larger scale organization of dna in eukaryotes. Physical Biology. 2009;6:025008 (11pp). [PubMed]
32. Bohn M, Heermann DW. Conformational properties of compact polymers. J Chem Phys. 2009;130:174901. [PubMed]
33. Phillips JE, Corces VG. CTCF: Master Weaver of the Genome. Cell. 2009;137:1194–1211. [PMC free article] [PubMed]
34. Galande S, Purbey PK, Notani D, Kumar PP. The third dimension of gene regulation: organization of dynamic chromatin loopscape by SATB1. Current Opinion in Genetics & Development. 2007;17:408–414. [PubMed]
35. de Gennes PG. Scaling concepts in polymer physics. Ithaca, N.Y.: Cornell University Press; 1979.
36. Langowski J. Polymer chain models of DNA and chromatin. Eur Phys J E Soft Matter. 2006;19:241–249. [PubMed]
37. Sikorav JL, Jannink G. Kinetics of chromosome condensation in the presence of topoisomerases - a phantom chain model. Biophys J. 1994;66:827–837. [PubMed]
38. Rosa A, Everaers R. Structure and dynamics of interphase chromosomes. PLoS Comput Biol. 2008;4:e1000153. [PMC free article] [PubMed]
39. Carmesin I, Kremer K. The bond fluctuation method: a new effective algorithm for the dynamics of polymers in all spatial dimensions. Macromolecules. 1988;21:2819–2823.
40. Yokota H, van den Engh G, Hearst J, Sachs R, Trask B. Evidence for the organization of chromatin in megabase pair-sized loops arranged along a random walk path in the human G0/G1 interphase nucleus. J Cell Biol. 1995;130:1239–1249. [PMC free article] [PubMed]
41. Grosberg AY, Khokhlov AR. Statistical Physics of Macromolecules. AIP Press; 1994.
42. Bohn M, Heermann DW, van Driel R. Random loop model for long polymers. Phys Rev E. 2007;76:051805. [PubMed]
43. Bode J, Goetze S, Heng H, Krawetz S, Benham C. From DNA structure to gene expression: mediators of nuclear compartmentalization and dynamics. Chromosome Res. 2003;11:435–445. [PubMed]
44. St-Jean P, Vaillant C, Audit B, Arneodo A. Spontaneous emergence of sequence-dependent rosettelike folding of chromatin fiber. Phys Rev E. 2008;77:061923. [PubMed]
45. Lieberman-Aiden E, van Berkum NL, Williams L, Imakaev M, Ragoczy T, et al. Comprehensive mapping of long-range interactions reveals folding principles of the human genome. Science. 2009;326:289–293. [PMC free article] [PubMed]
46. Duplantier B. Statistical mechanics of polymer networks of any topology. Journal of Statistical Physics. 1989;54:581–680.
47. Belmont AS. Visualizing chromosome dynamics with GFP. Trends Cell Biol. 2001;11:250–257. [PubMed]
48. Dundr M, Misteli T. Functional architecture in the cell nucleus. Biochem J. 2001;356:297–310. [PubMed]
49. Gilbert N, Gilchrist S, Bickmore WA. Chromatin organization in the mammalian nucleus. Int Rev Cytol. 2005;242:283–336. [PubMed]
50. Münkel C, Eils R, Dietzel S, Zink D, Mehring C, et al. Compartmentalization of interphase chromosomes observed in simulation and experiment. J Mol Biol. 1999;285:1053–1065. [PubMed]
51. Branco M, Pombo A. Intermingling of chromosome territories in interphase suggests role in translocations and transcription-dependent associations. PLoS Biol. 2006;4:97–104. [PMC free article] [PubMed]
52. Bolzer A, Kreth G, Solovei I, Koehler D, Saracoglu K, et al. Three-dimensional maps of all chromosomes in human male fibroblast nuclei and prometaphase rosettes. PLoS Biol. 2005;5:e157. [PMC free article] [PubMed]
53. Bohn M, Heermann DW. Topological interactions between ring polymers: Implications for chromatin loops. Journal of Chemical Physics. 2010 in press. [PubMed]
54. Marenduzzo D, Orlandini E. Topological and entropic repulsion in biopolymers. Journal of Statistical Mechanics: Theory and Experiment. 2009;2009:L09002.
55. Bruns W. Angular-correlations within coiled polymer-molecules. Makromolekulare Chemie-Theory And Simulations. 1992;1:287–293.
56. Khalil A, Grant J, Caddle L, Atzema E, Mills K, et al. Chromosome territories have a highly nonspherical morphology and nonrandom positioning. Chromosome Res. 2007;15:899–916. [PubMed]
57. Eils R, Dietzel S, Bertin E, Schröck E, Speicher MR, et al. Three-dimensional reconstruction of painted human interphase chromosomes: active and inactive X chromosome territories have similar volumes but differ in shape and surface structure. J Cell Biol. 1996;135:1427–1440. [PMC free article] [PubMed]
58. Dekker J, Rippe K, Dekker M, Kleckner N. Capturing Chromosome Conformation. Science. 2002;295:1306–1311. [PubMed]
59. Grosberg A, Rabin Y, Havlin S, Neer A. Crumpled globule model of the three-dimensional structure of DNA. EPL (Europhysics Letters) 1993;23:373–378.
60. Sikorski A. Monte Carlo study of catenated ring polymers. Polymer. 1994;35:3792–3794.
61. Lua R, Borovinskiy AL, Grosberg AY. Fractal and statistical properties of large compact polymers: a computational study. Polymer. 2004;45:717–731.
62. Binder K, Heermann DW. Monte Carlo Simulations in Statistical Physics. An Introduction. Springer, 4 edition; 2002.
63. Deutsch HP, Binder K. Interdiffusion and self-diffusion in polymer mixtures: A Monte Carlo study. J Chem Phys. 1991;94:2294–2304.
64. Paul W, Binder K, Heermann DW, Kremer K. Crossover scaling in semidilute polymer-solutions - a Monte-Carlo test. Journal De Physique Ii. 1991;1:37–60.
65. Müller M, Wittmer JP, Cates ME. Topological effects in ring polymers: A computer simulation study. Phys Rev E. 1996;53:5063–5074. [PubMed]
66. Müller M, Wittmer JP, Cates ME. Topological effects in ring polymers. II. Influence of persistence length. Phys Rev E. 2000;61:4078–4089. [PubMed]
67. Binder K, editor. Monte Carlo and Molecular Dynamics Simulations in Polymer Science. Oxford University Press; 1995.
68. Sokal A. Monte carlo methods in statistical mechanics: foundations and new algorithms. Functional integration 1996

Articles from PLoS ONE are provided here courtesy of Public Library of Science