<title>Introduction</title> <p>Eukaryotic genomes are organized in sets of chromosomes which are made up by a single continuous piece of DNA and associated proteins ##UREF##0##[1]##. During cell division (mitosis) chromosomes adopt a compact form which is suitable for transport and which can be discerned in a light microscope. During periods of normal cell activity (interphase), chromosomes decondense. More than 100 years ago, Rabl discovered that interphase chromosomes in newt and <italic>Drosophila</italic> remain organized in distinct territories ##UREF##1##[2]##. During the last twenty years similar territories of various shapes have been observed in many organisms ##REF##16920639##[3]##, a notable exception being budding yeast whose chromosomes appear to mix freely ##REF##17251970##[4]##,##REF##8943041##[5]##.</p> <p>The function of these territories, the mechanism responsible for their formation, and the reasons for the differences between species are still unclear ##REF##17251970##[4]##,##REF##17213034##[6]##. In this paper we investigate, if the observed interphase structure and dynamics are the consequence of a generic polymer effect, the preservation of the local topological state in solutions of entangled chain molecules undergoing Brownian motion. This effect plays an important role for the viscoelastic properties of polymeric systems ##UREF##2##[7]##,##REF##14764875##[8]##. In the present context, Sikorav and Jannink ##REF##8011915##[9]## assumed that interphase nuclei behave as equilibrated polymer solutions and estimated the disentanglement time <italic>τ</italic>
<sub>d</sub> of condensing metaphase chromosomes as <italic>τ</italic>
<sub>d</sub> = 1.5×10<sup>−7</sup> (#nucleosomes)<sup>3</sup> s, where “#nucleosomes” is the total number of nucleosomes in a chromosome. A human chromosome of typical size ≈100 mega-basepairs (Mbp) has ≈500,000 nucleosomes ##UREF##0##[1]##, i.e., <italic>τ</italic>
<sub>d</sub>≈2×10<sup>10</sup> s (≈500 years). From this prohibitively high estimate Sikorav and Jannink concluded that the process requires substantial topoisomerase-II (topo-II) activity.</p> <p>Here we reverse the argument. We suggest that interphase nuclei never equilibrate and behave like semi-dilute solutions of <italic>un</italic>entangled ring polymers which are known to segregate due to topological constraints ##UREF##3##[10]##. Within these territories, individual genomic sites are highly mobile and accessible. However, the structure of interphase and metaphase chromosomes remains largely identical from a topological point of view. Thus, instead of being a problem to be overcome by evolution, slow equilibration of long chromosomes <italic>accelerates</italic> the reverse process of chromosome condensation.</p> <title>Experimental Evidence and Polymer Theory</title> <p>Nowadays, the large-length scale structure of decondensed chromosomes can be experimentally studied using Fluorescence in situ Hybridization (FISH): nucleic acids are chemically modified to incorporate fluorescent probes and specific sequences on single chromosomes can be detected ##REF##11283701##[11]##. In particular, it is possible to mark different portions of the genome (chromosome painting) and to determine locations of and spatial distances between targeted sites ##REF##11283701##[11]##. Chromosome painting indicates that chromosome territories in human nuclei have an ellipsoidal shape with radii of the order of 1 µm ##REF##17251970##[4]##. In contrast and as already discovered by Rabl, the interphase nuclei of organisms like newt or <italic>Drosophila</italic> are organized in elongated territories oriented between two poles of the nucleus ##UREF##1##[2]##,##REF##16920639##[3]##. Furthemore, there are also organisms such as budding yeast whose chromosomes appear to mix freely or, at least, considerably less organized ##REF##17251970##[4]##,##REF##8943041##[5]##. The localization of territories inside the nucleus exhibits regular patterns: gene-rich chromosomes in human lymphocytes preferably locate in the nuclear interior while gene-poor chromosomes are typically found closer to the periphery ##REF##15111398##[12]##,##REF##11159939##[13]##; in contrast, in human fibroblasts positioning of territories was shown to correlate with chromosome size and not with its gene content ##REF##15839726##[14]##. In general, interactions between specific chromosome regions and structural elements within the nuclear envelope, such as nuclear pores or nuclear lamina, are believed to shape chromatin organization ##REF##11882311##[15]##.</p> <p>Data on the (relative) position and motion of target sites provide further insight into the organization of interphase chromosomes. In ##FIG##0##Figure 1A## we show average spatial distances between targeted sites as a function of their genomic separation. The figure contains FISH data for yeast chromosomes 6 and 14 (Chr6 and Chr14, brown ○) ##REF##15545610##[16]##, human chromosome 4 (Chr4, blue ○ and ◊) ##REF##7708711##[17]## and <italic>Drosophila</italic> chromosome 2L (Chr2L, orange and green ○) ##REF##15371546##[18]##. In the latter case, orange symbols refer to embryos in DS5 phase and green symbols to the DS1 phase which appears later in the cell cycle ##REF##9531544##[19]##. Two-dimensional spatial distances between sites on Chr4 measured in fibroblasts cells fixed on microscope slides ##REF##7708711##[17]## were here rescaled by 3/2 to obtain the corresponding 3 d distances.</p> <p>Observations for the various organisms agree on short length scales and coincide with the known properties of the (30 nm) chromatin fiber ##REF##15545610##[16]##. Given further the rather structureless appearance of interphase nuclei in the light microscope, a useful starting point for a theoretical description is a confined, equilibrated semi-dilute solution of “worm-like” chromatin fibers ##REF##15111398##[12]##. Within the worm-like chain (WLC) model, the fiber statistics can be calculated analytically ##UREF##2##[7]##. It is characterized by a crossover from rigid rod to random coil behavior at a characteristic length scale, the Kuhn length <italic>l</italic>
<sub>K</sub>≈300 nm of the 30 nm chromatin fiber ##REF##15545610##[16]##,##UREF##4##[20]##. Consider two points located at <italic>N</italic>
<sub>1</sub> and <italic>N</italic>
<sub>2</sub> Mbp from one chosen end of the fiber. They are separated by <italic>L</italic> = |<italic>N</italic>
<sub>2</sub>−<italic>N</italic>
<sub>1</sub>|×10 µm Mbp<sup>−1</sup> along the contour of the chromatin fiber ##REF##15545610##[16]##. The fiber is essentially stiff with a mean square spatial distance <italic>R</italic>
<sup>2</sup>(<italic>L</italic>≫<italic>l</italic>
<sub>K</sub>) = <italic>L</italic>
<sup>2</sup> on small scales and bent by thermal fluctuations on large scales with <italic>R</italic>
<sup>2</sup>(<italic>L</italic>) = <italic>l</italic>
<sub>K</sub>
<italic>L</italic>. The full crossover is described by ##UREF##5##[21]##
(black continuous line, ##FIG##0##Figure 1##). In particular, Equation 1 holds in the bulk of semi-dilute solutions where chains strongly overlap. Given the typical contour length of <italic>L</italic>
<sub>c</sub> = 1 mm of the chromatin fiber of a human chromosome with ≈100 Mbp, the expected chain extension of  largely exceeds the nuclear radius of 5 µm. In an equilibrated solution, the fibers should fill the nucleus homogeneously with mean-square internal distances saturating at a limiting value (black dashed line, ##FIG##0##Figure 1A##). (The exact probability distribution function of the square internal distances <italic>R</italic>
<sup>2</sup>(|<italic>N</italic>
<sub>2</sub>−<italic>N</italic>
<sub>1</sub>|) of a polymer without self-interactions obeys diffusion equation ##UREF##2##[7]## with null boundary conditions (in our case the boundary is the sphere which models the nucleus).) In contrast, the smaller yeast (<italic>S. cerevisiae</italic>) chromosome (≈1 Mbp, <italic>L</italic>
<sub>c</sub>≈0.01 mm, ) should be only weakly affected by a confinement to its nucleus of ≈1 µm radius ##REF##11739961##[22]## while <italic>Drosophila</italic> chromosomes (≈20 Mbp, <italic>L</italic>
<sub>c</sub>≈0.2 mm, ) in embryonic cells (for which FISH data are avalaible ##REF##15371546##[18]##) are confined inside nuclei whose radius grows from ≈2 µm to ≈5 µm in ≈30 minutes. Not surprisingly, the large-scale statistics of human and <italic>Drosophila</italic> chromosomes does not agree at all with the predictions of a WLC model assuming confinement at the scale of the nucleus (##FIG##0##Figure 1A##). Rather, the data reflect the different territory shapes observed for the two species. Note, however, that confinement on large scales alone cannot explain the unexpectedly small distances on intermediate scales |<italic>N</italic>
<sub>2</sub>−<italic>N</italic>
<sub>1</sub>|&gt;4 Mbp for Chr4 (blue ◊).</p> <p>There is less data available for the dynamics of interphase chromosomes. In mammalian cells chromatin domains of ∼1 µm diameter display little or no motion in a period of several hours ##REF##11356361##[23]##. Cabal et al. ##REF##16760982##[24]## followed the motion of a marked active gene (<italic>GAL</italic>) in in-vivo yeast nuclei. They observed a mean-square displacement (msd) <italic>g</italic>
<sub>1</sub>(<italic>t</italic> = 100 s)≈0.1 µm<sup>2</sup> for their largest observation interval, i.e., much less than the typical territory size in organisms with larger chromosomes. In particular, the authors reported anomalous diffusion with <italic>g</italic>
<sub>1</sub>(<italic>t</italic>)∼<italic>t</italic>
<sup>0.4</sup>.</p> <p>To rationalize this result, it is again useful to consider “worm-like” chromatin fibers in equilibrated semi-dilute solutions at typical nuclear densities. Neglecting entanglement effects, <italic>g</italic>
<sub>1</sub>(<italic>t</italic>) displays crossovers between different regimes: (1) <italic>g</italic>
<sub>1</sub>(<italic>t</italic>)∼<italic>t</italic>
<sup>0.75</sup> up to length scales of ≈1 Kuhn length ##UREF##6##[25]##; (2) <italic>g</italic>
<sub>1</sub>(<italic>t</italic>)∼<italic>t</italic>
<sup>0.5</sup> (Rouse behavior) up to length scales of the chain radius of gyration 
##UREF##2##[7]##; and (3) <italic>g</italic>
<sub>1</sub>(<italic>t</italic>)∼<italic>t</italic> at larger times, when the monomer motion is dominated by the center-of-mass diffusion (cyan line, ##FIG##1##Figure 2##).</p> <p>In semidilute solutions, linear chains with a contour length exceeding a characteristic value, <italic>L</italic>≫<italic>L</italic>
<sub>e</sub>, become mutually entangled, leading to confinement to a tube-like region following the coarse-grained chain contour and a drastically altered, “reptation” dynamics ##UREF##2##[7]##. Estimating <italic>L</italic>
<sub>e</sub> is not a trivial task. How strongly linear polymers entangle with each other depends on their stiffness and on the contour length density of the polymer melt or solution ##REF##18247995##[26]##. The latter is most suitably expressed in terms of the density of Kuhn segments, <italic>ρ</italic>
<sub>K</sub>. In loosely entangled systems with  the mean-free chain length between collisions is larger than the Kuhn length, leading to random coil behavior between entanglement points. In contrast, for  filaments are tightly entangled and exhibit only small bending fluctuations between entanglement points. For a solution of chromatin fibers at a typical nuclear density of ≈0.012 bp/nm<sup>3</sup> and a Kuhn length of 300 nm (##TAB##0##Table 1##) , i.e. the system is loosely entangled, but close to the crossover between the limiting cases. The entanglement contour length, <italic>L</italic>
<sub>e</sub>, can be estimated via ##REF##18247995##[26]##
, yielding <italic>L</italic>
<sub>e</sub>≈1.2 µm or four times the Kuhn length. To a first approximation, chains can thus be considered to be flexible on the tube scale, i.e., we expect around a msd of  a crossover from Rouse behavior to a <italic>g</italic>
<sub>1</sub>(<italic>t</italic>)∼<italic>t</italic>
<sup>0.25</sup> regime characteristic of reptation ##UREF##2##[7]##. Interestingly, this estimate coincides with the observations of Cabal et al. ##REF##16760982##[24]##, who reported an intermediate effective power law <italic>g</italic>
<sub>1</sub>(<italic>t</italic>)∼<italic>t</italic>
<sup>0.4</sup> for msd 0.05 µm<sup>2</sup>≤<italic>g</italic>
<sub>1</sub>(<italic>t</italic>)≤0.17 µm<sup>2</sup>. Using their data, we can obtain estimates for the entanglement time, <italic>τ</italic>
<sub>e</sub>≃32 s, as well as for the disentanglement times, <italic>τ</italic>
<sub>d</sub>≈<italic>τ</italic>
<sub>e</sub>(<italic>L</italic>
<sub>c</sub>/<italic>L</italic>
<sub>e</sub>)<sup>3</sup>
##UREF##2##[7]##, of the order of <italic>τ</italic>
<sub>d</sub>≃2×10<sup>4</sup> s, 2×10<sup>8</sup> s (≈5 years) and 2×10<sup>10</sup> s (≈500 years) for yeast, <italic>Drosophila</italic> and human chromosomes, respectively. Since this exceeds the life time of the entire organism (not to mention the much shorter cell cycle of most animal cells ##UREF##0##[1]##), <italic>Drosophila</italic> and human chromosomes do not have the time to equilibrate during interphase. (This conclusion does not change, if we take into account entanglement relaxation via topo-II discussed in ##REF##8011915##[9]##. At best, this mechanism could completely remove the barrier for chain crossing, thus converting the system to a solution of phantom chains whose relaxation time is given by the Rouse time <italic>τ</italic>
<sub>R</sub>≈<italic>τ</italic>
<sub>e</sub>(<italic>L</italic>
<sub>c</sub>/<italic>L</italic>
<sub>e</sub>)<sup>2</sup>
##UREF##2##[7]##. Yeast, <italic>Drosophila</italic> and human chromosomes would relax in, respectively, 2×10<sup>3</sup> s, 10<sup>6</sup> s (≈10 days), and 2×10<sup>7</sup> s (≈250 days).)</p> <p>While the discussion up to this point has shed some light on various aspects of the structure and dynamics of interphase chromosomes, we have so far evaded the central question, the origin of the observed chromosome territories. A priori, segregation or (micro) phase separation due to small chemical differences between chains is a common phenomenon in polymeric systems ##UREF##5##[21]##. Organisms could, in principle, render different chromosomes immiscible by a labeling technique akin to chromosome painting. In practice, it is difficult to conceive a corresponding, self-organizing molecular mechanism. Here we propose that the formation of chromosome territories could be related to a different, less well-known effect, the segregation of <italic>unentangled</italic> ring polymers in concentrated solutions due to topological barriers ##UREF##3##[10]##,##REF##16885211##[27]##. Well-separated metaphase chromosomes are clearly unentangled at the onset of interphase. Initially, decondensing chains can only rearrange locally and spread out uniformly without changing the global topological state. Up to the extremely long relaxation times for large chromosomes, interphase nuclei should therefore show a behavior similar to concentrated solutions of unentangled ring polymers. In particular, the chromosomes should remain segregated!</p> <p>It is instructive to compare this explanation to previously published models describing interphase chromosomes as <italic>equilibrium</italic> structures. The unexpectedly small distances on intermediate scales |<italic>N</italic>
<sub>2</sub>−<italic>N</italic>
<sub>1</sub>|&gt;4 Mbp for Chr4 (blue ◊) were rationalized in terms of giant loops of fibers departing from an underlying (protein) backbone ##REF##7708711##[17]## or alternatively, in terms of random loops on all length scales resulting from specific chromatin-chromatin interactions ##UREF##7##[28]##. Simulations of a multi-loop subcompartment polymer model reproduced the experimental observations on human Chr4, by <italic>imposing</italic> (and hence not explaining) confinement to a spherical territory ##UREF##4##[20]##,##REF##9887267##[29]##. We do not exclude the possibility of such contacts. However, we claim that territories should also form, if the involved proteins are disabled. For the inverse test—to keep the linking proteins, but to equilibrate a nucleus with disabled local topology conservation—it would be instructive to investigate the structure of nuclei in long-living cells arrested in interphase and to devise ways to maximize the efficiency of topo-II. (This conclusion does not change, if we take into account entanglement relaxation via topo-II discussed in ##REF##8011915##[9]##. At best, this mechanism could completely remove the barrier for chain crossing, thus converting the system to a solution of phantom chains whose relaxation time is given by the Rouse time <italic>τ</italic>
<sub>R</sub>≈<italic>τ</italic>
<sub>e</sub>(<italic>L</italic>
<sub>c</sub>/<italic>L</italic>
<sub>e</sub>)<sup>2</sup>
##UREF##2##[7]##. Yeast, <italic>Drosophila</italic> and human chromosomes would relax in, respectively, 2×10<sup>3</sup> s, 10<sup>6</sup> s (≈10 days), and 2×10<sup>7</sup> s (≈250 days).) We note that a few cross-links or attachment points to a residual skeleton would be sufficient to suppress chromosome equilibration via reptation ##UREF##8##[30]##. Long-lived contacts could thus stabilize the observed structures without being at their origin.</p> <p>How much of the experimental observations can be explained by this topology-conserving, parameter-<italic>free</italic>, minimal model of decondensing chromosomes? Unfortunately, it is difficult to derive quantitative predictions from an analytical theory due to the non-trivial initial conformation, the simultaneous presence of various crossovers (stiff/flexible, loosely/tightly entangled), and the lack of a theory describing the conformational statistics and dynamics of the unentangled ring polymer melts. We have therefore resorted to Molecular Dynamics (MD) computer simulations as a tool which allows us to study the model without further approximations.</p> <title>The Model</title> <p>With a spatial discretization of 30 nm (corresponding to the bead diameter), the employed generic bead-spring polymer model ##UREF##9##[31]## accounts for the linear connectivity, self-avoidance and bending stiffness of the chromatin fiber (<xref ref-type="sec" rid="s4">Materials and Methods</xref>). In particular, there is an energy barrier of 70<italic>K</italic>
<sub>B</sub>
<italic>T</italic> to prevent chain crossing ##UREF##10##[32]##. We emphasize that our description does not invoke any protein-like machinery as the nuclear matrix ##REF##10712500##[33]##. Furthermore, we neglect local changes of the chromatin fiber as they occur, e.g., as a result of chromatin remodeling during transcription ##REF##10364545##[34]##, because these processes do not alter the local topological state of the fiber and are therefore irrelevant for the phenomenon we discuss. This argument does not hold for the action of topo-II whose role is precisely to (dis)entangle DNA allowing strands to cross ##REF##8011915##[9]##,##UREF##4##[20]##. Non-directed topology changes with a particular rate could be included by suitable modifications of the energy barrier for chain crossing ##REF##10044759##[35]##. Similarly, it is straightforward to include (protein-mediated) interactions between specific DNA sites or effects such as confinement by or anchoring to the nuclear envelope ##REF##11283701##[11]##,##REF##10712500##[33]##,##REF##15568979##[36]##. However, here we concentrate on the generic case of decondensing long, internally and mutually unentangled polymers at total concentrations far above the overlap concentration.</p> <p>As initial states of our simulations we chose linear or ring-shaped helical structures remnant of metaphase chromosomes (<xref ref-type="sec" rid="s4">Materials and Methods</xref>). Given the anisotropic shape of our “metaphase” chromosomes, we were interested to see how the decondensation is affected by the presence of other chains. The l.h.s. panel in ##FIG##2##Figure 3## shows the initial chromosome conformations in our simulations on a common scale, indicated by a typical human nuclear radius of 5 µm. For <italic>Drosophila</italic> (marked “B”, only one chromosome is shown for clarity) we assumed that 8 Chr2L model chromosomes are initially aligned along the common axis of a rectangular simulation box (<italic>nematic</italic> orientation). In the case of yeast (marked “C”) and of the human (marked “A”), we followed the decondensation of 6 respectively 4 chromosomes of equal size which were oriented randomly in the simulation box ##REF##15839726##[14]##. For comparison we have also studied ring shaped chromosomes (see inset of ##FIG##0##Fig. 1F##) of different length under conditions corresponding to those of the human cell nucleus. 27 small rings (<italic>L</italic>
<sub>c</sub> = 2×2.7 Mbp) were randomly oriented inside the simulation box, while for larger rings (<italic>L</italic>
<sub>c</sub> = 2×48.6 Mbp and <italic>L</italic>
<sub>c</sub> = 2×97.2 Mbp) we limited ourselves to simulations of single chains in contact with their periodic copies in adjacent simulation cells. The setup as a ring allows us to eliminate chain end effects which otherwise play an important role.</p> <p>All simulations were performed in a constant isotropic pressure ensemble using rectangular simulation boxes with three independently fluctuating linear dimensions. The imposed pressure leads to the final density corresponding to the experimental value of ≈0.012 bp/nm<sup>3</sup> for human nuclei or 10% of volume fraction of chromatin (##TAB##0##Table 1##). This appears a reasonable choice because the experimental density in yeast nuclei is only two times lower (≈0.006 bp/nm<sup>3</sup>, ##TAB##0##Table 1##), while the rapid growing size of <italic>Drosophila</italic> embryos nuclei ##REF##9531544##[19]## does not allow for a univocal choice. We emphasize that the employed periodic boundary conditions do <italic>not</italic> introduce confinement to the finite volume of the simulation box. Using properly unfolded coordinates, chains can extend over arbitrarily large distances (see ##FIG##3##Figure 4## for the example of a MD simulation using a similar model but with a strongly reduced barrier for chain crossing). To give an idea of the computational effort, we consider the example of Chr4, where we simulated four model chromosomes of half of the actual length of Chr4. Each chromosome is modeled as a chain of 32,400 beads with a total contour length of 10<sup>−3</sup> m or 97.2 Mbp. Using ≈7×10<sup>4</sup> single-processor CPU-hours on a CRAY XD1 parallel computer, we followed the dynamics over 12×10<sup>6</sup> MD time steps. The comparison to the measured single-site mobility for yeast ##REF##16760982##[24]## in ##FIG##1##Figure 2## suggests the value of <italic>τ</italic> ≈2×10<sup>−2</sup> s used throughout the paper. According to this estimate, we followed the chain dynamics over 240,000 s (≈3 days) of real time. However, it is clear that more experimental data are needed to reliably fix the absolute time scale of our simulations.</p>
<title>Materials and Methods</title> <title>The Bead-Spring Polymer Model</title> <p>To model chromatin fiber we used the generic bead-spring polymer model of Kremer and Grest ##UREF##9##[31]##. Chains are composed of interacting beads of diameter <italic>σ</italic>. There are three types of interactions: <italic>U</italic>
<sub>LJ</sub>, <italic>U</italic>
<sub>FENE</sub>, and <italic>U</italic>
<sub>stiff</sub>. <italic>U</italic>
<sub>LJ</sub> is a shifted, purely repulsive Lennard-Jones potentialbetween any two monomers. The potentialgives the additional interaction between nearest neighbours along the chain. Finally, the stiffness of the fiber is modeled bywhere <italic>θ</italic> is the angle formed by the oriented unit vectors of two consecutive bonds. The bead diameter <italic>σ</italic> equals 30 nm, thus each bead corresponds to 3,000 bp ##UREF##1##[2]##. The other parameters are given by <italic>R</italic>
<sub>0</sub> = 1.5<italic>σ</italic>, <italic>k</italic> = 30.0<italic>ε</italic>/<italic>σ</italic>
<sup>2</sup>, and temperature <italic>K</italic>
<sub>B</sub>
<italic>T</italic> = 1.0<italic>ε</italic>
##UREF##9##[31]##. Since the Kuhn's length of the 30-nm fiber is <italic>l</italic>
<sub>K</sub> = 300 nm = 10<italic>σ</italic>
##REF##15545610##[16]##,##UREF##4##[20]##, the stiffness constant <italic>βk<sub>θ</sub></italic> is taken = 5 ##UREF##12##[42]##.</p> <title>Design of the Initial Configuration</title> <p>Experimental evidence suggests that metaphase chromosomes are folded into loops 30–100 kbp long (rosettes), arranged radially along the axis of the chromatid (see ##REF##8011915##[9]## and references therein). Metaphase chromosome are ≈700 nm thick and the length of each chromosome is related to its size ##REF##17213034##[6]##. On average, a typical human chromosome has 10<sup>8</sup> bp, i.e., a contour length <italic>L</italic>
<sub>c</sub> = 10<sup>6</sup> nm and a length <italic>h</italic>
<sub>chr</sub>≈5,000 nm ##REF##4894149##[43]##.</p> <p>As a starting configuration, we have placed chain beads along the generalized helix described by the equation:where <italic>r</italic>
<sub>chr</sub> = 12<italic>σ</italic>, <italic>p</italic> = <italic>σ</italic>, and <italic>h</italic>
<sub>chr</sub> = 170<italic>σ</italic>. With this choice of parameters, the length of each turn is approximately = 200<italic>σ</italic>. Given an average loop length of 50 kbp ≈17<italic>σ</italic>, we have ≈12 loops/turn. That fixes the remaining parameters <italic>k</italic> = 6 and <italic>x</italic> = 0.38.</p> <p>The contour lengths of the <italic>simulated</italic> human Chr4 and <italic>Drosophila</italic> Chr2L are, respectively, <italic>L</italic>
<sub>c</sub> = 97.2 Mbp and <italic>L</italic>
<sub>c</sub> = 21.6 Mbp, which corresponds to chains composed of 32,400 and 7,200 beads.</p> <p>The ring setup is described by the following equation:where the period , , and <italic>r<sub>t</sub></italic> = 42<italic>σ</italic>.</p> <title>Details of the Simulations</title> <p>The simulations have been performed in a <italic>constant isotropic pressure</italic> ensemble. Since the value of the pressure which must be imposed to the system is not known <italic>a priori</italic>, we have designed the following procedure: the decondensation of a ring chain of contour length <italic>L</italic>
<sub>c</sub> = 5.4 Mbp (1, 800 beads) has been simulated in a <italic>constant volume</italic> ensemble and the average diagonal components of the pressure tensor <italic>P<sub>αβ</sub></italic> (<italic>α</italic>,<italic>β</italic> = <italic>x</italic>,<italic>y</italic>,<italic>z</italic>) ##UREF##13##[44]## have been calculated. We have found <italic>P<sub>xx</sub></italic> = <italic>P<sub>yy</sub></italic> = <italic>P<sub>zz</sub></italic> = 0.01 and this value has been used throughout the paper. However, simulations of yeast chromosomes <italic>dynamics</italic> have been performed in a <italic>constant volume</italic> ensemble, because in the constant pressure ensemble the small system size leads to large <italic>unphysical</italic> fluctuations of the simulation box. In this constant volume ensemble, simulated yeast chromosomes have been initially arranged in an <italic>equilibrated</italic> configuration.</p> <p>We have chosen the integration time <italic>t</italic>
<sub>int</sub> = 0.012<italic>τ</italic>, where <italic>τ</italic> = <italic>σ</italic>(<italic>m</italic>/<italic>ε</italic>)<sup>1/2</sup> is the Lennard-Jones time and <italic>m</italic> is the <italic>mass</italic> of each bead ##UREF##9##[31]##. Each simulation runs up to time 10<sup>9</sup>
<italic>t</italic>
<sub>int</sub> = 12×10<sup>6</sup>
<italic>τ</italic>. Since we have sampled each 10<sup>6</sup>
<italic>t</italic>
<sub>int</sub>, each running produces 10<sup>3</sup> configurations.</p> <p>Notice that the time behavior of the msd of the 6 inner beads (<italic>g</italic>
<sub>1</sub>(<italic>t</italic>)) (##FIG##1##Figure 2##) has been calculated after shifting to the frame where the center of mass of the <italic>whole</italic> system is at rest. For human and <italic>Drosophila</italic> chromosomes, <italic>g</italic>
<sub>1</sub>(<italic>t</italic>) and  have been calculated neglecting the first 6×10<sup>5</sup>
<italic>τ</italic>≈12,000 s of the simulated trajectory.</p> <title>Gyration Tensor</title> <p>The gyration tensor <italic>T</italic> of an object composed of <italic>N</italic> beads is the 3×3 symmetric matrix whose elements are , where <italic>r<sub>l</sub></italic> is the vector pointing at the <italic>l</italic>th bead,  is the center of mass of the beads and <italic>i</italic>,<italic>j</italic> = 1,2,3 are the three indices for cartesian components. The trace of <italic>T</italic>,  where  is the square of the gyration radius of the object ##UREF##2##[7]##. It also equals Λ<sub>1</sub>+Λ<sub>2</sub>+Λ<sub>3</sub>, where Λ<italic><sub>i</sub></italic> is the <italic>i</italic>th eigenvalue of <italic>T</italic>. For objects with spherical symmetry, Λ<sub>1</sub> = Λ<sub>2</sub> = Λ<sub>3</sub>. Then, differences between the eigenvalues measure the anisotropy of the object ##UREF##11##[38]##.</p>
<title>Results</title> <p>Since there are no attractive interactions in our model of the chromatin fiber, the bent and kinked initial state is unstable and unfolds rapidly. The initial rapid expansion stops when chromosomes come into contact with others, including their periodic replicas in adjacent simulation cells. Our simulation time is sufficient to mix and equilibrate the short (1 Mbp) yeast chromosomes (##FIG##2##Figure 3##). Fast equilibration of yeast chromosomes explains why apparently there is no territorial organization in yeast nuclei ##REF##8943041##[5]##. (Most chromosomes in yeast have a size smaller than 1 Mbp, corresponding to a disentanglement time comparable to the time duration of the relative interphase (∼1 hour ##REF##6392855##[37]##).) In this case memory of the initial condition is rapidly lost: a simulation where the chains are initially prepared as rods oriented along the same direction produces similar results (data not shown).</p> <p>The much longer <italic>Drosophila</italic> and human chromosomes remain confined to distinct territories (##FIG##2##Figure 3##). For the nematically ordered initial state we assumed for <italic>Drosophila</italic>, we observed that decondensation leads to the formation of Rabl-like elongated territories. In contrast, in isotropically arranged copies of the human Chr4, the preferred axial expansion is suppressed and the resulting territory shapes resemble elongated ellipsoids. Our ring chromosomes essentially reproduce the latter behavior. More quantitatively, the shape of the human Chr4 territory is described by the average of the 3 eigenvalues Λ<sub>1</sub>, Λ<sub>2</sub>, and Λ<sub>3</sub> of the corresponding gyration tensor (##UREF##11##[38]## and <xref ref-type="sec" rid="s4">Materials and Methods</xref>) and the ratios of the two largest eigenvalues Λ<sub>1</sub> and Λ<sub>2</sub> over Λ<sub>3</sub> are two quantities which could experimentally be tested. We have found that averaging over the configurations of all the possible sections of half the total ring size gives Λ<sub>1</sub>∶Λ<sub>2</sub>∶Λ<sub>3</sub> = 6.4(±1.4)∶1.9(±0.4)∶1.0 (2×97.2 Mbp), Λ<sub>1</sub>∶Λ<sub>2</sub>∶Λ<sub>3</sub> = 5.5(±1.2)∶2.1(±0.5)∶1.0 (2×48.6 Mbp) and Λ<sub>1</sub>∶Λ<sub>2</sub>∶Λ<sub>3</sub> = 6.9(±0.7)∶2.2(±0.2)∶1.0 (2×2.7 Mbp), while averaging over all the 4-Chr4 configurations gives Λ<sub>1</sub>∶Λ<sub>2</sub>∶Λ<sub>3</sub> = 11.0(±1.2)∶1.5(±0.3)∶1.0. The ∼2 times larger value found in the latter case is probably an artifact of the setup (see also below).</p> <p>In ##FIG##0##Figure 1## (panels B to E) we compare the simulation results for mean-square spatial distances between marked sites on the chromosomes to the experimental findings shown in ##FIG##0##Figure 1A## and discussed in the introduction. Gray lines represent spatial distances between sites in the initial, compact “metaphase” configuration. To give an impression of the time dependence of the results, we have averaged the <italic>R</italic>
<sup>2</sup>(<italic>N</italic>
<sub>2</sub>−<italic>N</italic>
<sub>1</sub>|) curves over three exponentially spaced time windows: 240 s&lt;<italic>t</italic>&lt;2,400 s, 2,400 s&lt;<italic>t</italic>&lt;24,000 s, 24,000 s&lt;<italic>t</italic>&lt;240,000 s (dark red, magenta and cyan lines respectively). In panels B–D we show results averaged over the entire length of the simulated <italic>Drosophila</italic> chromosome Chr2L, yeast Chr6 and Chr14 and human Chr4. While the former two are in excellent agreement with the experimental data, this is not the case for our first set of results for the human Chr4. Here simulation and experimental data agree quantitatively only on short length scales. It turns out, that there are different explanations for the deviations on intermediate and on large length scales.</p> <p>
##FIG##0##Figure 1F## shows the corresponding comparison to our data for ring chromosomes. In this case, the experimentally observed conformational statistics of human Chr4 on large scales is perfectly reproduced. In fact, when we reanalyzed data for the linear chromosome assuming the existence of a “centromere-hinge,” we found nearly perfect agreement with the ring data (not shown). This suggests to interpret the (nearly linear) large scale behavior of our simulation results in ##FIG##0##Figure 1D## as an artifact of the straight initial configuration.</p> <p>Interestingly, the simulation data follow the experimentally observed effective power law <italic>R</italic>
<sup>2</sup>∼<italic>L</italic>
<sup>2<italic>ν</italic></sup> with <italic>ν</italic>≈0.32 ##REF##9887267##[29]## already on intermediate scales (<italic>L</italic>&gt;1 Mbp). (We note that the relation between the square of the gyration radius  and the mean square internal distances of a polymer <italic>R</italic>
<sup>2</sup>(|<italic>N</italic>
<sub>2</sub>−<italic>N</italic>
<sub>1</sub>|), 
##UREF##2##[7]##, is compatible with chromosome territories being compact objects with . However, the reverse conclusion ##UREF##4##[20]##,##REF##9887267##[29]## is incorrect: globular polymer conformations also follow , but do not have a fractal structure where the same exponent characterises the entire chain conformation (see, for example the dashed line in ##FIG##0##Figure 1##). Simple confinement alone cannot explain the chain structure.) This behavior seems to be robust, since all our simulation data for linear chains and rings of different size beautifully collapse onto each other. Similar, quasi-fractal structures were reported in ##UREF##3##[10]##. Taken together this suggests that our ring samples are relatively well-equilibrated and that (in agreement with our working hypothesis) long, unentangled linear chains initially relax to a very similar structure. However, we still require an explanation for the deviations between this apparently quite robust prediction and the experimental data in ##FIG##0##Figure 1D##.</p> <p>Reptation theory ##UREF##8##[30]## would suggest that the further equilibration of linear chromosomes proceeds by a very slow escape of the terminal parts of the chain from their initial environment. Qualitatively, this effect is directly observable in ##FIG##2##Figure 3## where we have marked the terminal parts of our model chromosomes in red. Interestingly, the experimental data for the spatial distances between sites with genomic distances in the Mbp range on human Chr4 were determined in the ≈4.5 Mbp 4p16.3 region which is located at the <italic>end</italic> of <italic>p</italic>-arm ##REF##7708711##[17]##. A good way to quantify the consequences is to measure <italic>R</italic>
<sup>2</sup>(|<italic>N</italic>
<sub>2</sub>−<italic>N</italic>
<sub>1</sub>|,<italic>N</italic>
<sub>1</sub> = 0), i.e., mean-square spatial distances between the chain <italic>ends</italic> and points along the fiber (##FIG##0##Figure 1E## and ##FIG##4##Figure 5##). These distances show a stronger time dependence than results averaged over the entire chain. In particular, they follow the WLC prediction up to much larger contour length distances before crossing over to the bulk averages. The point of departure from the WLC prediction can be used to estimate up to which distance from the end the chains are equilibrated after a certain time. (The temporal behavior of the ratio between the escaped portion of the chain and the whole contour length <italic>L</italic>
<sub>c</sub> at short times <italic>t</italic> is compatible with the power-law ∼<italic>t</italic>
<sup>1/2</sup> predicted by reptation theory (##UREF##8##[30]##, data not shown).) The comparison to the experimental data in ##FIG##0##Figure 1E## suggests that the 4p16.3 region of the human Chr4 was nearly equilibrated in the experimental situation. We emphasize that we expect spatial distances between marked sites in the interior of long chromosomes to fall onto the corresponding simulation data in ##FIG##0##Figure 1D and 1F##. This is at least qualitatively supported by a remark in ##REF##1388286##[39]## where van den Engh et al. report the more centrally located 6p21 region on human Chr6 to be more compact than the 4p16.3 region near the end of Chr4.</p> <p>As a final point, we turn to the dynamics of chromosomes during interphase. ##FIG##1##Figure 2## shows the msd of the 6 inner beads (<italic>g</italic>
<sub>1</sub>(<italic>t</italic>)) after the complete (yeast) and initial (human, <italic>Drosophila</italic>) relaxation in comparison to the respective gyration radii . By adjusting the time scale of the simulations, the simulation data can be mapped on the experimental results from ##REF##16760982##[24]##. The good agreement indicates that our model provides a simple, quantitative explanation for the reported anomalous diffusion. In particular, the model reproduces the dynamic (entanglement) length scale with no adjustable parameters. Moreover, the dynamic range of the simulation data (0.1 s&lt;<italic>t</italic>&lt;3 days) significantly exceeds the observation window in ##REF##16760982##[24]##, allowing us to extrapolate to longer times. For comparison, in ##FIG##1##Figure 2## we have included data for equilibrated yeast chromosome solutions from simulation of a model without excluded volume interactions and topological barriers (cyan line). All simulations exhibit identical short time behavior in agreement with theoretical expectations ##UREF##6##[25]##. A Rouse regime for  is only observable in simulations without topological barriers. The yeast chromosomes in <italic>equilibrated</italic> entangled solutions exhibit instead <italic>g</italic>
<sub>1</sub>(<italic>t</italic>)∼<italic>t</italic>
<sup>0.25</sup> reptation dynamics. Interestingly, our data for human and <italic>Drosophila</italic> chromosomes show the same behavior in spite of the very different microscopic topological state and the (on these scales) weakly perturbed chain statistics. (The small deviations from the yeast data are artifacts of the constant-pressure simulations used for human and <italic>Drosophila</italic> chromosomes.) In our simulations the asymptotic free diffusion regime—where the center of mass has moved farther than the chain size ##UREF##2##[7]##—is reached only for yeast chromosomes. (Note that the corresponding simulation data cannot be compared directly to experiments, since we have neglected the nuclear confinement in the present study.) While human and <italic>Drosophila</italic> chromosomes remain confined to their territories and do not equilibrate, individual sites are extremely dynamic. Cabal et al. ##REF##16760982##[24]## reported that invidual loci on yeast chromosomes explore regions of linear size ∼0.4 µm. The simulations indicate that msd's of ∼1 µm<sup>2</sup> are reached on the time scale of ∼5 hours.</p>
<title>Discussion</title> <p>We have studied the decondensation, structure and dynamics of interphase chromosomes using Molecular Dynamics simulations of a bead-spring model of the 30 nm chromatin fiber. Our results suggest that <italic>for sufficiently long</italic> chromosomes territories form as a consequence of a generic polymer effect, the preservation of the local topological state in solutions of long chain molecules undergoing Brownian motion. In fact, we argue that such interphase nuclei never equilibrate and behave like concentrated solutions of unentangled ring polymers, which segregate due to topological constraints ##UREF##3##[10]##. Such cases are also know from material science where they result in unusual material properties ##REF##16041376##[40]##. The slow kinetics leads to memory effects. For example, different chromosome arrangements in the nucleus at the end of metaphase provide a rationale for the different territory shapes observed in humans and flies. Similarly, the negligible relative motion of territories provides a natural explanation for the tendency of chromosomes to “reappear” at the end of interphase at similar relative positions as those occupied at the end of the preceeding anaphase ##REF##12654243##[41]##. Our simulations confirm this tendency: the centers of mass of the large human chromosomes remain confined to small regions of linear size ≈0.1 µm and retain their relative positions (##FIG##5##Figure 6##). In contrast, individual sites are extremely dynamic inside the territories and explore much larger regions up to a linear size of ≈1 µm (##FIG##1##Figure 2##). We emphasize that conservation of the local topology during <italic>decondensation</italic> discussed in the present work considerably simplifies the reverse process of chromosome <italic>condensation</italic> at the end of interphase, a process which takes only about 1 hour in most animal cells ##UREF##0##[1]## and which is difficult to conceive for fully equilibrated nuclei ##REF##8011915##[9]##.</p> <p>Obviously, there is more to the structure and dynamics of eukaryotic nuclei than can be captured by the present model in its basic form. However, our results suggest that effects such as active transport ##REF##11739961##[22]##, chromosome anchoring to the nuclear envelope ##REF##15568979##[36]##, replication ##REF##10364545##[34]## and homologous pairing ##REF##9531544##[19]## should be investigated in the framework of the polymer description presented here. As we have shown, computer simulations along the present lines can now reach the relevant time and length scales.</p>
<p>Conceived and designed the experiments: AR RE. Performed the experiments: AR RE. Analyzed the data: AR RE. Contributed reagents/materials/analysis tools: AR RE. Wrote the paper: AR RE.</p> <p>During interphase chromosomes decondense, but fluorescent in situ hybridization experiments reveal the existence of distinct territories occupied by individual chromosomes inside the nuclei of most eukaryotic cells. We use computer simulations to show that the existence and stability of territories is a kinetic effect that can be explained without invoking an underlying nuclear scaffold or protein-mediated interactions between DNA sequences. In particular, we show that the experimentally observed territory shapes and spatial distances between marked chromosome sites for human, <italic>Drosophila</italic>, and budding yeast chromosomes can be reproduced by a parameter-free minimal model of decondensing chromosomes. Our results suggest that the observed interphase structure and dynamics are due to generic polymer effects: confined Brownian motion conserving the local topological state of long chain molecules and segregation of mutually unentangled chains due to topological constraints.</p> <title>Author Summary</title> <p>Eukaryotic genomes are organized in sets of chromosomes. Each chromosome consists of a single continuous DNA double-helix and associated proteins that organize locally in the form of a <italic>chromatin</italic> fiber. During cell division (mitosis) chromosomes adopt a compact form that is suitable for transport. During periods of normal cell activity (interphase), they decondense inside the cell nucleus. Being long-chain molecules (in the case of human chromosomes the contour length of the chromatin fiber is on the order of 1 mm), the random thermal motion of interphase chromatin fibers is hindered by entanglements, similar to those restricting the manipulation of a knotted ball of wool. We have studied the consequences of this effect using computer simulations. Most importantly, we find that entanglement effects cause sufficiently long chromosomes to remain segregated during interphase and to form “territories.” Our model (1) reproduces currently avaliable experimental results for the existence and shape of territories as well as for the internal chromosome structure and dynamics in interphase nuclei and (2) explains why entanglement effects do <italic>not</italic> interfere with the reverse process of chromosome condensation at the end of interphase.</p>
<p>AR thanks C. Svaneborg for stimulating discussions.</p>
<fig id="pcbi-1000153-g001" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1000153.g001</object-id><label>Figure 1</label><caption><title>Experimental FISH data for spatial distances <italic>R</italic>
<sup>2</sup>(|<italic>N</italic>
<sub>2</sub>−<italic>N</italic>
<sub>1</sub>|) between targeted chromosome sites compared to the estimates based on the WLC model (A) and results from our simulations (B–F).</title><p>
<italic>Brown ○</italic>: <italic>Saccharomyces cerevisiae</italic> Chr6 and Chr14 ##REF##15545610##[16]##. <italic>Blue ○ and ◊</italic>: <italic>Homo Sapiens</italic> Chr4, |<italic>N</italic>
<sub>2</sub>−<italic>N</italic>
<sub>1</sub>|&lt;4.5 Mbp and |<italic>N</italic>
<sub>2</sub>−<italic>N</italic>
<sub>1</sub>|&gt;4.5 Mbp, respectively ##REF##7708711##[17]##. <italic>Orange and green ○</italic>: <italic>Drosophila melanogaster</italic> Chr2L, DS5 and DS1 embryos respectively ##REF##15371546##[18]##. DS5 and DS1 are two phases of cell cycle. DS1 appears later. <italic>Black continuous line</italic>: Mean-square internal distances predicted by the WLC model, Equation 1. (A) <italic>Black dashed line</italic>: Mean-square internal distances of an ideal polymer chain inside a spherical nucleus of 5 µm radius ##UREF##2##[7]##. (The exact probability distribution function of the square internal distances <italic>R</italic>
<sup>2</sup>(|<italic>N</italic>
<sub>2</sub>−<italic>N</italic>
<sub>1</sub>|) of a polymer without self-interactions obeys diffusion equation ##UREF##2##[7]## with null boundary conditions (in our case the boundary is the sphere which models the nucleus).)While data for Chr4 and Chr2L show a reasonable agreement at short-length scales, the apparent large-length scale Chr4 behavior <italic>L</italic>
<sup>2/3</sup>
##REF##9887267##[29]## contrasts with the observed <italic>L</italic>
<sup>2</sup> for Chr2L. The insets show two schematic drawings of the Chr4 territory in a human nucleus (blue) and of Chr2L in Rabl phase in a <italic>Drosophila</italic> nucleus (orange). (B–E) Gray lines represent internal distances in the initial, “metaphase-like” chromosome configuration (<xref ref-type="sec" rid="s4">Materials and Methods</xref>). Internal distances in simulated chromosomes have been averaged over 3 time windows of exponentially growing size: 240 s&lt;<italic>t</italic>&lt;2,400 s (dark red line), 2,400 s&lt;<italic>t</italic>&lt;24,000 s (magenta line) and 24,000 s&lt;<italic>t</italic>&lt;240,000 s (cyan line). Since yeast chromosomes rapidly equilibrate only averages over the first 24,000 s are here reported (panel C). In panel E, <italic>N</italic>
<sub>1</sub> = 0, i.e., has been fixed at the origin of the chain to make equilibration of the chain ends evident. (F) Data from simulations of three ring polymers of decreasing half-size <italic>L</italic>
<sub>c</sub> = 97.2, 48.6, and 2.7 Mbp (green, magenta and red lines respectively). Mean distances seem to extrapolate to an effective power law ∼<italic>L</italic>
<sup>2/3</sup>. <italic>Inset</italic>: Initial (left) and final (right) conformation of a (randomly chosen) half of the largest (2×97.2 Mbp) simulated ring chromosome.</p></caption></fig> <fig id="pcbi-1000153-g002" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1000153.g002</object-id><label>Figure 2</label><caption><title>Time behavior of the msd of the six inner beads (<italic>g</italic>
<sub>1</sub>(<italic>t</italic>), continuous lines), compared to the average square gyration radius  (horizontal dashed lines) of the whole chromosome and measurements of the msd of the active <italic>GAL</italic> gene inside in vivo yeast nuclei ##REF##16760982##[24]## (purple dots).</title><p>For comparison, <italic>g</italic>
<sub>1</sub>(<italic>t</italic>) for yeast chromosomes <italic>without</italic> topological constraints has been shown (cyan line). On short time scales, our model reproduces the typical dynamics of semi-flexible polymers with <italic>g</italic>
<sub>1</sub>(<italic>t</italic>)∼<italic>t</italic>
<sup>0.75</sup>
##UREF##6##[25]##. For the model with constraints, there is no extended Rouse regime due to the insufficient separation of Kuhn and entanglement length. Nevertheless, we observe the characteristic <italic>g</italic>
<sub>1</sub>(<italic>t</italic>)∼<italic>t</italic>
<sup>0.25</sup> regime for entangled, flexible polymers ##UREF##2##[7]##.</p></caption></fig> <fig id="pcbi-1000153-g003" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1000153.g003</object-id><label>Figure 3</label><caption><title>Initial (“metaphase-like”, left) and final (right) configurations of human Chr4 (A), of <italic>Drosophila</italic> Chr2L (B) and of yeast Chr6 and Chr14 (C) shown together with the spherical nucleus (black circle) of 10 µm in diameter and the corresponding simulation boxes (in black).</title><p>For the blue configuration in A and for the configuration B, we have highlighted in red the two terminal parts up to 4.5 Mbp. In Chr4, this corresponds to the terminal 4p16.3 region ##REF##7708711##[17]##. (A) Simultaneous decondensation of 4 model chromosomes half the size the human Chr4. (B) Decondensation of 1 model chromosome the size the <italic>Drosophila</italic> Chr2L. The final elongated shape qualitatively resembles a Rabl-like territory. (C) Simultaneous decondensation of 6 model chromosomes the size the yeast Chr6 and Chr14. Arrows points at magnified versions of the same configurations. Lack of chromosome territoriality is evident.</p></caption></fig> <fig id="pcbi-1000153-g004" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1000153.g004</object-id><label>Figure 4</label><caption><title>Human Chr4 territories are less stable if the energy barrier against chain crossing is switched off.</title><p>The swelling from the initial “metaphase” configuration is monitored through the time behavior of the gyration radius 
##UREF##2##[7]##, where <italic>r<sub>l</sub></italic>(<italic>t</italic>) is the position vector of the <italic>l</italic>th bead and <italic>r</italic>
<sub>cm</sub>(<italic>t</italic>) is the center of mass of the configuration at time <italic>t</italic>. Without barrier, chromosomes swell easier and have larger size (green and red lines, (A)). Comparison amongst internal distances between two sites located at <italic>N</italic>
<sub>1</sub> and <italic>N</italic>
<sub>2</sub> Mbp from one chosen end of the fiber and avalaible experimental data reflects this behavior (B). We have averaged over 3 time windows of exponentially growing size: 240 s&lt;<italic>t</italic>&lt;2,400 s (dark red line), 2,400 s&lt;<italic>t</italic>&lt;24,000 s (magenta line) and 24,000 s&lt;<italic>t</italic>&lt;240,000 s (cyan line). In particular, we notice that the fortuitous agreement of the magenta line with the data is lost due to the fast relaxation to equilibrium. The gray line corresponds to internal distances in the initial configuration. As expected (C), the final configuration of human Chr4 without energy barrier occupies a larger volume and is more random-walk-like than the ones where the energy barrier has been included.</p></caption></fig> <fig id="pcbi-1000153-g005" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1000153.g005</object-id><label>Figure 5</label><caption><title>Mean square spatial distances <italic>R</italic>
<sup>2</sup>(|<italic>N</italic>
<sub>2</sub>−<italic>N</italic>
<sub>1</sub>|) between a site of the fiber located at <italic>N</italic>
<sub>2</sub> Mbp from one chosen end of the chain and the end (here located at <italic>N</italic>
<sub>1</sub> = 0): comparison between simulated and the avalaible experimental data on <italic>Drosophila</italic> Chr2L (left) and yeast Chr6 and Chr14 (right).</title><p>Gray lines represent internal distances in the initial, “metaphase-like” chromosome configuration (<xref ref-type="sec" rid="s4">Materials and Methods</xref>). Internal distances have been averaged over 3 time windows of exponentially growing size: 240 s&lt;<italic>t</italic>&lt;2,400 s (dark red line), 2,400 s&lt;<italic>t</italic>&lt;24,000 s (magenta line) and 24,000 s&lt;<italic>t</italic>&lt;240,000 s (cyan line). Since yeast chromosomes rapidly equilibrate only averages over the first 24,000 s are here reported. The black continuous line is the plot of the average internal distances predicted by the WLC model, Equation 1.</p></caption></fig> <fig id="pcbi-1000153-g006" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1000153.g006</object-id><label>Figure 6</label><caption><title>Three dimensional spatial trajectories of the centers of mass of the 4 simulated human Chr4.</title><p>The color code used corresponds to the snapshots A shown in ##FIG##1##Figure 2##. Motion resembles confined diffusion.</p></caption></fig>
<table-wrap id="pcbi-1000153-t001" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1000153.t001</object-id><label>Table 1</label><caption><title>Summary of the relevant physical parameters for the polymer model of interphase chromosomes.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Parameter</td><td align="left" rowspan="1" colspan="1">Value</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">Typical nuclear radius of a human cell ##UREF##4##[20]##
</td><td align="left" rowspan="1" colspan="1">5 µm</td></tr><tr><td align="left" rowspan="1" colspan="1">Radius of the yeast (<italic>S. Cerevisiae</italic>) nucleus ##REF##11739961##[22]##
</td><td align="left" rowspan="1" colspan="1">1 µm</td></tr><tr><td align="left" rowspan="1" colspan="1">Length of the diploid human genome<xref ref-type="table-fn" rid="nt101">a</xref>
</td><td align="left" rowspan="1" colspan="1">6×10<sup>9</sup> bp</td></tr><tr><td align="left" rowspan="1" colspan="1">Length of the diploid <italic>Drosophila</italic> genome<xref ref-type="table-fn" rid="nt101">a</xref>
</td><td align="left" rowspan="1" colspan="1">3×10<sup>8</sup> bp</td></tr><tr><td align="left" rowspan="1" colspan="1">Length of the diploid yeast (<italic>S. Cerevisiae</italic>) genome<xref ref-type="table-fn" rid="nt101">a</xref>
</td><td align="left" rowspan="1" colspan="1">2×10<sup>7</sup> bp</td></tr><tr><td align="left" rowspan="1" colspan="1">Compaction ratio of chromatin ##REF##15545610##[16]##
</td><td align="left" rowspan="1" colspan="1">10<sup>2</sup> bp/nm</td></tr><tr><td align="left" rowspan="1" colspan="1">Kuhn length of chromatin ##REF##15545610##[16]##, <italic>l</italic>
<sub>K</sub>
</td><td align="left" rowspan="1" colspan="1">300 nm</td></tr><tr><td align="left" rowspan="1" colspan="1">Volume fraction of chromatin</td><td align="left" rowspan="1" colspan="1">10%</td></tr></tbody></table></alternatives></table-wrap>
<disp-formula><label>(1)</label></disp-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <disp-formula></disp-formula> <disp-formula></disp-formula> <disp-formula></disp-formula> <disp-formula></disp-formula> <disp-formula></disp-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula>
<table-wrap-foot><fn id="nt101"><label>a</label><p>See, e.g., the website <ext-link ext-link-type="uri" xlink:href="http://www.ensembl.org/index.html">http://www.ensembl.org/index.html</ext-link>.</p></fn></table-wrap-foot> <fn-group><fn fn-type="COI-statement"><p>The authors have declared that no competing interests exist.</p></fn><fn fn-type="financial-disclosure"><p>RE is supported by a chair of excellence grant from the Agence Nationale de Recherche (France).</p></fn></fn-group>
<graphic xlink:href="pcbi.1000153.g001"/> <graphic xlink:href="pcbi.1000153.e001.jpg" mimetype="image" position="float"/> <inline-graphic xlink:href="pcbi.1000153.e002.jpg" mimetype="image"/> <inline-graphic xlink:href="pcbi.1000153.e003.jpg" mimetype="image"/> <inline-graphic xlink:href="pcbi.1000153.e004.jpg" mimetype="image"/> <inline-graphic xlink:href="pcbi.1000153.e005.jpg" mimetype="image"/> <inline-graphic xlink:href="pcbi.1000153.e006.jpg" mimetype="image"/> <graphic xlink:href="pcbi.1000153.g002"/> <inline-graphic xlink:href="pcbi.1000153.e007.jpg" mimetype="image"/> <inline-graphic xlink:href="pcbi.1000153.e008.jpg" mimetype="image"/> <inline-graphic xlink:href="pcbi.1000153.e009.jpg" mimetype="image"/> <inline-graphic xlink:href="pcbi.1000153.e010.jpg" mimetype="image"/> <inline-graphic xlink:href="pcbi.1000153.e011.jpg" mimetype="image"/> <graphic id="pcbi-1000153-t001-1" xlink:href="pcbi.1000153.t001"/> <graphic xlink:href="pcbi.1000153.g003"/> <inline-graphic xlink:href="pcbi.1000153.e012.jpg" mimetype="image"/> <graphic xlink:href="pcbi.1000153.g004"/> <inline-graphic xlink:href="pcbi.1000153.e013.jpg" mimetype="image"/> <inline-graphic xlink:href="pcbi.1000153.e014.jpg" mimetype="image"/> <inline-graphic xlink:href="pcbi.1000153.e015.jpg" mimetype="image"/> <inline-graphic xlink:href="pcbi.1000153.e016.jpg" mimetype="image"/> <graphic xlink:href="pcbi.1000153.g005"/> <inline-graphic xlink:href="pcbi.1000153.e017.jpg" mimetype="image"/> <inline-graphic xlink:href="pcbi.1000153.e018.jpg" mimetype="image"/> <graphic xlink:href="pcbi.1000153.g006"/> <graphic xlink:href="pcbi.1000153.e019.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000153.e020.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000153.e021.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000153.e022.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000153.e023.jpg" mimetype="image" position="float"/> <inline-graphic xlink:href="pcbi.1000153.e024.jpg" mimetype="image"/> <inline-graphic xlink:href="pcbi.1000153.e025.jpg" mimetype="image"/> <inline-graphic xlink:href="pcbi.1000153.e026.jpg" mimetype="image"/> <inline-graphic xlink:href="pcbi.1000153.e027.jpg" mimetype="image"/> <inline-graphic xlink:href="pcbi.1000153.e028.jpg" mimetype="image"/> <inline-graphic xlink:href="pcbi.1000153.e029.jpg" mimetype="image"/> <inline-graphic xlink:href="pcbi.1000153.e030.jpg" mimetype="image"/>
[{"label": ["1"], "element-citation": ["\n"], "surname": ["Alberts"], "given-names": ["B"], "year": ["2002"], "source": ["Molecular Biology of the Cell. 4th edition"], "publisher-loc": ["New York"], "publisher-name": ["Garland Science"]}, {"label": ["2"], "element-citation": ["\n"], "surname": ["Rabl"], "given-names": ["C"], "year": ["1885"], "article-title": ["\u00dcber Zellteilung."], "source": ["Morph Jb"], "volume": ["10"], "fpage": ["214"], "lpage": ["330"]}, {"label": ["7"], "element-citation": ["\n"], "surname": ["Doi", "Edwards"], "given-names": ["M", "SF"], "year": ["1986"], "source": ["The Theory of Polymer Dynamics"], "publisher-loc": ["New York"], "publisher-name": ["Oxford University Press"]}, {"label": ["10"], "element-citation": ["\n"], "surname": ["M\u00fcller", "Wittmer", "Cates"], "given-names": ["M", "JP", "ME"], "year": ["2000"], "article-title": ["Topological effects in ring polymers. II. Influence of persistence length."], "source": ["Phys Rev E"], "volume": ["61"], "fpage": ["4078"], "lpage": ["4089"]}, {"label": ["20"], "element-citation": ["\n"], "surname": ["M\u00fcnkel", "Langowski"], "given-names": ["C", "J"], "year": ["1998"], "article-title": ["Chromosome structure predicted by a polymer model."], "source": ["Phys Rev E"], "volume": ["57"], "fpage": ["5888"], "lpage": ["5896"]}, {"label": ["21"], "element-citation": ["\n"], "surname": ["Rubinstein", "Colby"], "given-names": ["M", "RH"], "year": ["2003"], "source": ["Polymer Physics"], "publisher-loc": ["New York"], "publisher-name": ["Oxford University Press"]}, {"label": ["25"], "element-citation": ["\n"], "surname": ["Farge", "Maggs"], "given-names": ["E", "AC"], "year": ["1993"], "article-title": ["Dynamic scattering from semiflexible polymers."], "source": ["Macromolecules"], "volume": ["26"], "fpage": ["5041"], "lpage": ["5044"]}, {"label": ["28"], "element-citation": ["\n"], "surname": ["Bohn", "Heermann", "van Driel"], "given-names": ["M", "DW", "R"], "year": ["2007"], "article-title": ["Random loop model for long polymers."], "source": ["Phys Rev E"], "volume": ["76"], "fpage": ["051805"]}, {"label": ["30"], "element-citation": ["\n"], "surname": ["de Gennes"], "given-names": ["PG"], "year": ["1971"], "article-title": ["Reptation of a polymer chain in presence of fixed obstacles."], "source": ["J Chem Phys"], "volume": ["55"], "fpage": ["572"], "lpage": ["579"]}, {"label": ["31"], "element-citation": ["\n"], "surname": ["Kremer", "Grest"], "given-names": ["K", "GS"], "year": ["1990"], "article-title": ["Dynamics of entangled linear polymer melts: a molecular-dynamics simulation."], "source": ["J Chem Phys"], "volume": ["92"], "fpage": ["5057"], "lpage": ["5086"]}, {"label": ["32"], "element-citation": ["\n"], "surname": ["Sukumaran", "Grest", "Kremer", "Everaers"], "given-names": ["SK", "GS", "K", "R"], "year": ["2005"], "article-title": ["Identifying the primitive path mesh in entangled polymer liquids."], "source": ["J Polym Sci B Polym Phys"], "volume": ["43"], "fpage": ["917"], "lpage": ["933"]}, {"label": ["38"], "element-citation": ["\n"], "surname": ["Bishop", "Michels"], "given-names": ["M", "JPJ"], "year": ["1986"], "article-title": ["Scaling in three-dimensional linear and ring polymers."], "source": ["J Chem Phys"], "volume": ["84"], "fpage": ["444"], "lpage": ["446"]}, {"label": ["42"], "element-citation": ["\n"], "surname": ["Auhl", "Everaers", "Grest", "Kremer", "Plimpton"], "given-names": ["R", "R", "GS", "K", "SJ"], "year": ["2003"], "article-title": ["Equilibration of long chain polymer melts in computer simulations."], "source": ["J Chem Phys"], "volume": ["119"], "fpage": ["12718"], "lpage": ["12728"]}, {"label": ["44"], "element-citation": ["\n"], "surname": ["Frenkel", "Smith"], "given-names": ["D", "B"], "year": ["2002"], "source": ["Understanding molecular simulation: from Algorithms to Applications"], "publisher-loc": ["New York"], "publisher-name": ["Academic Press"]}]
PLoS Comput Biol. 2008 Aug 22; 4(8):e1000153