Large Scale Chromosome Folding Is Stable against Local Changes in Chromatin Structure

Characterizing the link between small-scale chromatin structure and large-scale chromosome folding during interphase is a prerequisite for understanding transcription. Yet, this link remains poorly investigated. Here, we introduce a simple biophysical model where interphase chromosomes are described in terms of the folding of chromatin sequences composed of alternating blocks of fibers with different thicknesses and flexibilities, and we use it to study the influence of sequence disorder on chromosome behaviors in space and time. By employing extensive computer simulations, we thus demonstrate that chromosomes undergo noticeable conformational changes only on length-scales smaller than 105 basepairs and time-scales shorter than a few seconds, and we suggest there might exist effective upper bounds to the detection of chromosome reorganization in eukaryotes. We prove the relevance of our framework by modeling recent experimental FISH data on murine chromosomes.

A key determining factor in many important cellular processes as DNA transcription, for instance, the specific composition of the chromatin fiber sequence has a major influence on chromosome folding during interphase. Yet, how this is achieved in detail remains largely elusive. In this work, we explore this link by means of a novel quantitative computational polymer model for interphase chromosomes where the associated chromatin filaments are composed of mixtures of fibers with heterogeneous physical properties. Our work suggests a scenario where chromosomes undergo only limited reorganization, namely on length-scales below 10 5 basepairs and time-scales shorter than a few seconds. Our conclusions are supported by recent FISH data on murine chromosomes.

Introduction
Understanding how genomes fold within the crowded environment of the nucleus [1] during interphase represents a necessary step for the comprehension of important cellular processes such as gene expression and regulation [2]. The combined results of high-resolution microscopy [3,4] and mathematical and computer modelling [5,6] seem to suggest that genomes are organized hierarchically [7,8]. Each genome is partitioned into a set of single units, the chromosomes, and each chromosome is made of a single filament of DNA complexed around histone octamers to form a necklace-like fiber % 10 nm thick known as the 10nm chromatin fiber. In in vitro conditions close to the physiological ones, this fiber is observed to fold into a thicker, more compact structure known as the 30nm fiber [1], whose role and existence in vivo are nonetheless still quite debated [9,10]. On larger scales, chromosome conformation capture (3C) techniques [2] have shown that chromosomes appear organized in Topologically Associated Domains (TADs) of sizes ranging from % 0.1 to % 1 megabasepairs (Mbp). Chromosome loci within TADs interact frequently between themselves, but less frequently across different TADs. Finally, chromosomes do not spread inside the whole nucleus, rather they occupy well localized nuclear regions (the so-called "chromosome territories") which play a crucial role in gene expression and regulation [11].
While much progress on the causal relationship between chromosome structure and function has now been made, many fundamental aspects remain still obscure. One of these crucial issues concerns the link between chromosome (re)organization at various length and timescales and the spread of inhomogeneties in the sequence of the chromatin fiber which may arise from (1) selective epigenetic marks induced by chemical modifications of the histone tails [12], (2) nucleosomes arrangements in discrete nanodomains of different sizes [4], and (3) selective stimulation of particular kinds of genes [13] or entire gene families [14]. These events will produce modifications in the local polymer properties of the chromatin fiber (as its persistence length or the local compaction ratio) which might affect in turn the whole hierarchical folding of chromosome organization [9,10].
Motivated by these considerations, we present here the results of Molecular Dynamics computer simulations of a minimal polymer model for interphase chromosomes in order to quantify to what extent the simultaneous presence of chromatin fibers of heterogeneous composition (different thicknesses and flexibilities) is able to generate observable effects on the small-and large-scale structures and motions of the associated chromosome. The proposed model complements previous work by one of the authors [15][16][17][18] concerning the description of chromosome folding in terms of fundamental polymer physics. Similarly to other recent works discussing the explicit role of sequence disorder in chromatin [19] and chromosome behaviors [20,21], here we "deviate" from the description of the chromatin filament as a homopolymer and we discuss sequence effects in space and time through the introduction of controlled amounts of disorder in the chromatin sequence. In this way, we provide a quantitative description for many crucial aspects concerning the structure and dynamics of interphase chromosomes which are "spontaneously" driven by the physical properties of the underlying chromatin sequence with a definite copolymer structure.
In particular, by considering the two "extreme" cases of chromosomes made of: (1) short stretches of a thinner, more flexible fiber randomly interspersed in a "sea" of thicker fiber and (2) chromosomes partitioned into two distinct blocks of thinner and thicker fibers, we show that significative spatial and dynamical rearrangements of chromatin loci appear to be restricted to limited contour lengths (up to % 10 5 basepairs (bp)) and times scales (up to few seconds). Interestingly, there exists a limited range (10 4 −10 5 bp) of chromatin contour lengths where chromosome expansion is accompanied by an increase (rather, than a decrease) of the associated contacts between the fibers. We apply this framework to rationalize the outcome of recent experiments which employ fluorescent microscopy to monitor conformational changes of chromosomes that occur upon transcription activation or chromatin decondensing in mouse embryonic stem cells [13]. Finally, we argue that the effects discussed here are not the consequence of the details of the model, but involve more general aspects of polymer physics.

Results
In this section, we present the main results of our model by addressing in particular the specific question of how modifications in the small-scale properties of the chromatin fiber turn to affect chromosome behaviour on much larger scales. In order to model these modifications, we represented the linear chromosome filament as a heterogeneous copolymer made of two types of fibers with physical properties matching those of the 10nm and the 30nm chromatin fibers: 10nm fibers are modelled as completely flexible while 30nm fibers have a persistence length = 150 nm as in previous works [15][16][17][18]. Then, we considered several situations: either short (3000 bp) fragments of 10nm fibers are distributed randomly along the whole chromatin filament, or we studied chromosome structures where a single continuous 10nm filament is centered around the chromatin portion closer to or farther from the chromosome center of mass. For more details and technical issues, we invite the reader to look into the "Materials and Methods" section.
With its true existence in vivo appearing more and more debatable [9,10], employment of 30nm fibers might look inappropriate. However, our choice of physical parameters for the model was finally motivated by their previous use and proven consistency [15][16][17][18] with experimental data. We posit though that our results should remain qualitatively valid for more general models of non-homogeneous polymers.
In order to quantify chromosome changes, we have considered either spatial relationships between distal fragments along the chromatin sequence or the dynamic behaviour of chromatin loci. The former aspects were investigated by focusing mainly on the following two observables: (1) the mean-squared internal distance (hR 2 (L)i) and (2) the average contact frequency hp c (L)i between all possible pairs of genomic loci at any given genomic separation L. For hp c (L)i, we have adopted the choice that two monomers are in contact whenever their spatial distance is smaller than a conventional cut-off distance = 60 nm, corresponding to twice the thickness of the 30nm fiber. These two quantities are of particular experimental interest: in fact, hR 2 (L)i can be measured through fluorescence in-situ hybridization (FISH) [11], while hp c (L)i is the result of chromosome conformation capture (3C) techniques [22,23]. Moreover, hR 2 (L)i and hp c (L)i are useful to distinguish between complementary aspects of chromosome conformation, as was recently highlighted by Williamson et al. [24]. It is then interesting to monitor our systems by employing both tools. Dynamical aspects were instead discussed in terms of the mean-square displacement hδr 2 (τ)i of chromatin loci at lag time τ. This is also a quantity of notable experimental interest, as specific chromatin loci can now be followed in vivo by, e.g., fluorescent microscopy on GFP tagged chromosome sequences [25,26].
All results are presented as functions of the basepair content of the 10nm fibers present in the model chromosomes vs. the total chromosome basepair content, expressed in percent (see section "Materials and Methods" for details). For comparison, model chromosomes with no (0%) or entirely made of (100%) 10nm fibers are also discussed.

Large scale chromosome behaviour does not depend on the fiber composition
The results of our three case studies are summarized in Fig 1. It is visible that only length-scales smaller than L % 0.1 megabasepairs (Mbp) are affected with hR 2 (L)i expanding sensibly more than in the situation where chromosomes are composed only of 30nm fiber (panels A, C, E), while the behaviour at large scales remains unaffected. Moreover, in the case where the unfolded sequences are grouped into a single cluster (middle and bottom panels), there is no dependence on their spatial positioning with respect to the center of mass of the corresponding chromosome territory. Insensitivity of large scales to changes at small ones is also confirmed (panels B, D, F) by the analysis of contact frequencies, hp c (L)i, whose trend remains, in particular, compatible with the experimentally observed power-law hp c (L)i * L −1 [27]. Interestingly, instead of decreasing as expected from panels A, C and E, corresponding contact frequencies in the limited range [0.01 Mbp − 0.1 Mbp] also increase as a function of the 10nm fiber content (see in particular panel B). In order to understand this result, we use the mathematical relationship [16] between average contact frequencies hp c (L)i and distribution function p L (R) of internal distances R(L) given by: where r c = 60 nm is the cut-off distance for two monomers to form a contact. In particular, eq 1 implies that p L (R) should also increase as a function of the 10nm fiber content at given L. The inset of panel B confirms this behavior for L = 0.015 Mbp. This specific trend is due to the decreasing of the chromatin "effective persistence length" (from % 10 4 to % 10 2 basepairs for model chromosomes entirely made of 30nm and 10nm fibers, respectively) following from increasing amounts of 10nm fibers, which allows genomic loci to contact each other with enhanced probability. For a more quantitative view on chromosome reorganization at small scales, the same data for hR 2 (L)i and hp c (L)i were normalized to corresponding values for chromosome conformations made of 30nm fiber, see S1 Fig The latter, in particular, constitutes a result with non trivial experimental implications. First, in connection to some recently proposed protocols for the reconstruction of chromosome conformations based on 3C (reviewed in [28]) where a monotonous relationship between chromatin distances and contacts is often assumed, our finding demonstrates that some caution is needed in order to avoid systematic bias in the final reconstructed structure. Second, the recent puzzling result [24] where chromatin domains with a high propensity to form 3C contacts seem to undergo rather pronounced decondensation when monitored by using FISH can be understood by considering that the two experimental techniques sample very different intervals of the corresponding p L (R)'s: close to the average value or the median for FISH, around to the lower tail for 3C techniques, as shown in the insets of Fig 1B, 1D and 1F. Our work thus supports the important conclusion of references [24,29], namely that one needs to take into account both kinds of data to reconstruct correctly the shape of chromatin domains.
We complete the discussion by considering the full distributions p L (R) for L = 0.003, 0.015 and 3 Mbp, see Fig 2. For L = 0.003 and L = 0.015 Mbp (panels A-F) there are quantitative differences between the cases where the 10nm fiber is located at sparse random positions along the chromosome and where they form a single chromatin cluster, while for L = 3 Mbp (panels G-I) all distributions show no noticeable difference. In panels A-F, the largest peak corresponds to the most probable value of spatial distances between loci on the 30nm fiber region. Additionally, for random locations of 10nm fiber (panels A and D) corresponding p L (R)'s show a broader population of R values, while in the other cases there exist smaller peaks corresponding to the most probable distance between loci on the 10nm fiber region. Interestingly, similar distribution functions for spatial distances between chromatin loci seem to have been reported in yeast (see  [30]). Although a direct comparison between these experimental results and our data is not possible (our setup applies to large chromosomes, like mammalian ones), we are tempted to speculate that the results reported in reference [30] are a manifestation of the presence of chromatin fibers of different compositions. Alternatively, we may interpret the observed shape of the p L (R)'s in terms of a bias towards large spatial distances: in that respect, we report that a similar feature has been observed in human and mouse chromosomes [14,24] (including the experimental data that will be discussed in this work).
The mechanical behavior of the thinner fibers is influenced by the physical properties of the rest of the chromosome We have then considered those chromosome configurations with single long sequences of 10nm fiber and we have calculated hR 2 (L)i and hp c (L)i on these sequences, see Figs 3 and 4. Our results demonstrate that hR 2 (L)i increases systematically with respect to the analogous quantity for model chromosomes made of 30nm fiber, with no dependence on the specific location of 10nm fiber along the chromosome. This is also shown by the decreased contact probability. This insensitivity to positioning follows from the proposed picture [31] that chromosomes resemble a uniformly dense, semi-dilute solution of branched polymer rings. We expect then chromatin filaments to be similarly constrained or accessible regardless of their position along the genomic sequence.
In order to prove that the reported volume increase is not a fortuitous coincidence of the chosen sequence, we have also repeated the analysis on the same genomic region for the model chromosome with no 10nm fiber (0%) and the model chromosome entirely made of 10nm fiber (100%). Both hR 2 (L)i and hp c (L)i differ quantitatively from the case in which just one part of the chromosome is allowed to swell. Similarly to the previous section, the same results can be recast in terms of ratios to the corresponding quantities calculated for chromosome conformations with no 10nm fiber, see S2 and S3 Figs, which allows to appreciate to what extent chromosomes may effectively reorganize.
Because of the copolymer structure of our model chromosomes, the two types of fibers have different stiffnesses and thicknesses which cause the thinner fiber to move away from the thicker one. Consequently, hp c (L)i calculated for the 10nm fiber sequence tends to decrease and shows a scaling behavior *L −1.18±0.05 for L in the interval 0.01 − 1 Mbp which, being slightly steeper than the one for the entire chromosome, suggests an increase in the volume spanned by the fiber. Qualitatively, volume differences between chromosome regions with different transcription activities have been reported in a recent study on Drosophila [12]. Here the authors have studied three regions of the Drosophila genome, which, according to the histone modifications that they were bearing, were classified as active, inactive and polycombrepressed. They found that the polycomb-repressed region was the most compact, followed by the inactive region, while the active region was the least compact. These results encourage us to speculate that the difference between the active and inactive regions could be due to different polymer physical properties induced by the histone modifications.

Single clusters of thinner fiber increase the mobility of corresponding genomic loci
Here we discuss the impact of chromatin unfolding on the dynamics of the corresponding genomic loci. Specifically, we have considered the mean-square displacement  at time t and we implicitly assume average over specific monomer positions along the chromatin chain. In fact, this is an observable widely employed in many experiments monitoring the dynamic activity of specific chromatin loci, being especially suitable for comparing genome behavior in response to changes of the environment [32], or when the cell is targeted with drugs which are able to activate selectively certain types of genes [14]. In both cases, we notice a general increase of chromatin mobility as larger and larger portions of 30nm fiber unfold, and, at larger times, a trend which does not substantially depend on the small scale details of chromatin fiber. Not surprisingly, data at short times reflect in part the discussed results for chromatin structure: in particular, we notice that differences in chromatin mobility before and after chromatin unfolding in random locations are only visible below time-scales of about 5 seconds (panel A). Slightly larger discrepancies are observed in the other situation where chromosomes are organized as two separate domains (panel B). In this latter case, unfolded chromatin loci move on average more than folded ones, the latter displaying the same motion than in the case of the homogeneously folded chromosome (compare green vs. blue lines).   [13] (left) and simulation (right) data for the distribution functions, p L (R), of spatial distances R between the ends of the murine Ptn, Nrp1 and Sox6 genes. L * 0.22, 0.21, 0.55 Mbp is the genomic separation between the FISH probes placed at the end of the three considered genes, respectively. Experimental distances are measured before (eGFP, red) and after (DEL, black) stimulated chromatin decondensation. Numerical results come from three independent molecular dynamics runs simulating the unfolding of chromatin filaments of corresponding sizes. The rest of the chromosome remains in the 30nm fiber state. doi:10.1371/journal.pcbi.1004987.g006

Model verification: Experimental observation of local chromatin unfolding
In a recent work [13], Therizols et al. used FISH to show that, in embryonic stem cells, chromosome decondensation is sufficient to alter nuclear organization. Two sets of experiments were done: first, a viral transactivator was used to activate transcription of three genes (Ptn1, Nrp1 and Sox6). Alternatively, repositioning of the same genes was also observed after the treatment with an artificial peptide (DELQPASIDP) which decondenses chromatin without inducing transcription. By comparing the measured spatial distances between the two ends of each selected sequence, the authors concluded that nuclear organization is driven mainly by chromatin remodeling rather than transcription.
We have tested the predictions of our model by using data from these experiments. We have simulated the unfolding of a chromatin region of genomic length corresponding to the size of the specific gene we wanted to mimic. This situation corresponds to the second case studied in this work, namely the unfolding of a chromatin cluster. For a fair comparison, we have processed the experimental data as follows: we have reconstructed first the probability density distribution function for spatial distances between the ends of each gene, for the control condition (denoted as eGFP) and cells in which DELQPASIDP is recruited to the chromosome loci (denoted as DEL). Since FISH distances are recorded as two-dimensional vectors projected on the confocal plane, for the purpose of comparison a (large) set of three-dimensional distances with equivalent 2d projections was generated numerically by assuming random orientations of the 3d vectors in relation to the axis orthogonal to the confocal plane. It can be observed in Fig 6 panels A, C and E that, upon chromatin decondensation, the peak and the shape of the distribution change dramatically, in particular the distributions become wider. The same effect is displayed by our simulations (panels B, D and F). This comparison thus validates our result that larger genes tend to expand when unfolded.
While our model describes reasonably well the trend of the experimental data, the median values and the part of the distribution covering the larger distances, it seems to perform worse for the part of the distribution describing the small distances. This could be due to several factors which, to maintain the model simple, have been neglected. These include either protein linkers that physically bridge two sites along the genes forming a loop [20,33,34], either chemical modifications which could change the charge on some histones, thus inducing an effective electrostatic attraction [21] between regions of the gene that would make it more compact. Moreover, for the experiments using DEL, not all cells receive the exact same quantity of plasmid. That is, the experimental distributions of distances could be biased towards small R values because, in some cells, there was not enough plasmid to induce decondensation.

Polymer physics aspects
The choice for the physical parameters employed in this work as fibers flexibilities, thicknesses and excluded volume interactions (see section "Materials and Methods") was mainly motivated by reasons of simplicity. Understanding to what extent our results depend on these specific choices requires a deeper analysis. In order to isolate and quantify the effects of one single parameter, we have considered a simpler model chromosome where monomers have fixed nearest neighbor distance along the sequence equal to 30 nm. The polymer was then split in two complementary domains: the fiber in one domain has the same features of the 30nm fiber, while the physical properties of the fiber in the other domain are changed, in turns, to: (1) completely flexible fiber; (2) thickness equal to 10 nm; (3) larger strength for the corresponding excluded volume interaction. Each of these changes was applied, in turns, to a domain occupying 20%, 50% and respectively 100% of the polymer. S4 Fig (top panels) shows that, not surprisingly, the persistence length influences hR 2 (L)i and hp c (L)i only at small scales up to 0.1 Mbp, suggesting that one of the main effects of having a heterogeneous chromatin composition consists in a global change of the overall persistence length of the chromatin fiber. By varying instead the range of the excluded volume interactions without modifying the nominal persistence length we produce variations in hR 2 (L)i and hp c (L)i similar to the ones reported in Figs 3 and 4 (middle panels). Viceversa, when varying the strength of the interaction no change is observed (bottom panels). These observations suggest a leading role for the effective range of the excluded volume interaction: in our model, the thinner 10nm fiber tends to occupy a larger region because its own self-repulsion becomes less important when compared to the repulsion from the thicker 30nm fiber. Quite interestingly, a similar effect concerning the change in the excluded volume was observed in equilibrium simulations of generic block copolymers [35]. Finally, we stress that this effect is less visible in Fig 1 because of the base pair content assigned to the monomers: since the base pair content of the 10nm fiber region is in general low compared to 30nm fiber region, its effect is not visible when averaging over the entire chromosome.

Discussion
In this article, we have presented results of extensive Molecular Dynamics computer simulations of a coarse-grained polymer model for interphase chromosomes that extends previous numerical work [15][16][17][18] by introducing a crucial ingredient which was neglected before: namely, the presence of two kinds of chromatin fibers of different thickness and stiffness mutually interacting inside their own chromosome territory. Starting from a chromosome configuration made of a single and homogeneous 30nm fiber like filament, we have monitored chromosome spatial and temporal behaviors when this conformation is altered by the introduction of increasing amounts of 10nm-like fiber through the controlled unfolding of the thicker 30nm fiber.
The work shows that there exists detectable chromosome (re)organization for spatial scales smaller than 0.1 Mbp (Fig 1) and time-scales shorter than just a few seconds (Fig 5). Quite interestingly, these findings appear systematic and do not depend on the size of the chromosome region affected by the phenomenon of local unfolding or by its location along the chromosome. Interestingly, these results tend to suggest that experimental methodologies like FISH or HiC might be of little or no help in distinguishing between fibers of different compaction, unless they investigate genomic distances smaller than 0.1 Mbp. This prediction can be tested, for instance, by employing the recently developed oligonucleotide based FISH probes which seem to provide the necessary fine resolution [36,37]. An important conclusion of our work is that, although our model uses parameter values that can be associated to the traditional "10nm/ 30nm" chromatin fiber paradigm [9,10], our results reflect a generic physical effect [35] which ought to be observable in more general systems of crumpled polymers constituted of fibers with different thickness and/or stiffness (see section "Polymer physics aspects" and S4 Fig).
Our simulation protocol compares qualitatively well (see Fig 6) with experimental results on chromosome reorganization in mouse embryos treated with a synthetic transcription factor [13] which produces selective activation of a specific gene. In particular, we predict the observed shifting for distribution functions of spatial distances. Of course, it is quite possible that other mechanisms may explain the experimental results equally well. In fact, intentionally our model tends to neglect other important aspects of chromosome folding which have been highlighted recently by other authors: sequence-specific attractive interactions [21], protein linkers between chromatin fibers [20], mechanisms of active regulation [38,39], "loop extrusion" [33,34] involved in the reorganization of small chromosome domains, or the anchoring to the nuclear envelope or other nuclear organelles [40][41][42]. In the lack of a more quantitative analysis, we can only speculate on the fact that the inclusion of these mechanisms into our model might alter significantly the conclusions sketched here. On the other hand, this should represent an important stimulus for investigating further and more quantitatively the delicate relationship between chromosome structure and function.

Simulation protocol I. Force field
In this work, the chromatin fiber is modeled as a coarse-grained polymer chain, with monomer-monomer interactions described by analytical expressions similar to the ones used in our previous works [15][16][17] and suitably adapted to take into account the different sizes of 30nm and 10nm monomers.
The full Hamiltonian governing the system, H, consists of three terms: where N (see Table 1) is the total number of monomers constituting the ring polymer modeling the chromosome (see Section "Construction of model chromosome conformation" for details on this point) and i and j run over the indexes of the monomers. The latter are assumed to be numbered consecutively along the ring from one chosen reference monomer. The modulo-N indexing is implicitly assumed because of the ring periodicity. By taking the nominal monomer diameter of the 30nm chromatin fiber, σ = 30 nm = 3000 bp [16], as our unit of length, the vector position of the ith monomer,r i , the pairwise vector distance between monomers i and j,d i;j ¼r j Àr i , and its norm, d i, j , the energy terms in eq 2 are given by the following expressions: 1. The chain connectivity term, U FENE (i, i + 1) is expressed as: where R 0 = 1.5σ, k = 30.0/σ 2 and the thermal energy k B T equals 1.0 [43]. The shift distance Δ is used to tune the length of the bond between nearest neighbors monomers: (a) Δ = 0 if   both beads model the 30nm fiber, (b) Δ = σ/3 = 10 nm if one bead models the 30nm fiber and the other the 10nm fiber, (c) Δ = 2σ/3 = 20 nm if both beads model the 10nm fiber.
2. The bending energy term has the standard Kratky-Porod form (discretized worm-like chain): The persistence length between triplets of 30nm beads is the same as in our previous works [15][16][17], ξ p = 5.0σ = 150nm. In the other cases, we model the fiber as completely flexible, ξ p = 0. This is in agreement with measurements showing that the persistence length for unfolded chromosomes is close to % 6 nm [30].
3. The excluded volume interaction between distinct monomers (including consecutive ones) corresponds to a purely repulsive Lennard-Jones potential: where the shift distances Δ take the same values as the ones adopted in the U FENE term described above.

Simulation protocol II. Molecular Dynamics simulations
As in our previous work [15][16][17][18], chromosome dynamics was studied by performing fixed-volume Molecular Dynamics (MD) simulations with periodic boundary conditions at near-physiological fixed chromatin density ρ = 0.012bp/nm 3 . Note that periodic boundary conditions do not introduce confinement to the simulation box: using properly unfolded coordinates, the model chromatin fibers can extend over arbitrarily large distances [15].
The system dynamics was integrated by using LAMMPS [44] with Langevin thermostat in order to keep the temperature of the system fixed to 1.0k B T. Given the unit mass m 30 = 1 of the 30nm-bead, we fixed the mass of the 10nm-bead to m 10 ¼ 1

27
. The integration time step was fixed to t int = 0.001τ MD , where t MD ¼ s m 30 À Á 1=2 is the elementary Lennard-Jones time. γ = 0.5/τ MD is the friction coefficient [43] which takes into account the corresponding interaction with the background implicit solvent. The total length of each MD simulation run is = 3 Á 10 5 τ MD , with an overall computational effort ranging from a minimum of % 10 3 to a maximum of % 9 Á 10 4 hours of single CPU for "0%" and "100%" model chromosome conformations, respectively. Single chromosome conformations were sampled every 10 3 τ MD , implying 300 configurations per each run.
We have verified that our rings are well equilibrated by considering the mean-square distances, hR 2 (L)i, and mean contact probabilities, hp c (L)i, calculated on different sets of configurations log-spaced in time. As shown in S5 Fig for a 20% amount of 10nm fiber, all curves give the same results. Similar results are obtained for all simulated systems.

Simulation protocol III. Initial configuration
Construction of the model chromosome conformation. In spite of the complexity of the chromatin fiber and the nuclear medium, three-dimensional chromosome conformations are remarkably well described by generic polymer models [6,45,46]. In particular, it was suggested [15,47] that the experimentally observed [27] crumpled chromosome structure can be understood as the consequence of slow equilibration of chromatin fibers due to mutual chain uncrossability during thermal motion. As a consequence, chromosomes do not behave like equilibrated linear polymers in solution [48,49]. Instead, they appear rather similar to unlinked and unknotted circular (ring) polymers in entangled solution. In fact, under these conditions ring polymers are known to spontaneously segregate and form compact conformations [31,46,47], strikingly similar to images of chromosomes in live cells obtained by fluorescence techniques [11].
Due to the typically large size of mammalian chromosomes (each chromosome contains, on average, *10 8 basepairs of DNA), even minimalistic computational models would require the simulation of large polymer chains, made of tens of thousands of beads [15][16][17]. For these reasons, in this work we resort to our recent mixed Monte-Carlo/Molecular Dynamics multi-scale algorithm [31] in order to design a single, equilibrated ring polymer conformation at the nominal chromatin density of ρ = 0.012 bp/nm 3 (see previous section). The ring is constituted by N = 39154 monomer particles of linear size σ = 30 nm = 3000 basepairs, which corresponds to the average linear size of a typical mammalian chromosome of % 117.46 × 10 6 basepairs. By construction, the adopted protocol guarantees that the polymer has the nominal local features of the 30nm-chromatin fiber (stiffness, density and topology conservation) that have been already employed elsewhere [15][16][17][18]. For the details of the multi-scale protocol, we refer the reader to reference [31].
Local unfolding of chromatin fiber. To model the 30nm fiber, we employ the typical density of % 100 bp/nm [30], which we have already adopted successfully elsewhere [15][16][17][18]. For the 10nm fiber, we assume that the typical repeat length of a nucleosome (% 200 bp) is shared between the nucleosome core particle (147 bp of DNA for 10 nm of diameter) and linker DNA (% 60 bp) assumed completely stretched (% 20 nm). The total compaction of the 10nm fiber results being then % 200 bp/(10 nm + 20 nm)% 7 bp/nm. As shown next, and for practical purposes only, we round this quantity to % 11 bp/nm implying that each 30nm-fiber monomer (=3000 bp) swells into 27 10nm-fiber monomers. In fact, in order to model the effect of having alternate filaments of 10nm and 30nm fibers within the same chromosome, we started from the homogeneous model chromosome constructed as outlined in the previous section and we gradually replaced each single 30nm fiber monomer by 3 × 3 × 3 = 27 beads of smaller diameter = 10 nm arranged on a regular cubic lattice of linear step = 10 nm, as shown schematically in Fig 7 (panels A and B). This step was followed by a short (120 τ MD ) MD run employing the LAMMPS option NVE/LIMIT (http://lammps.sandia. gov) which, by limiting the maximum distance a particle can move in a single time-step, allows the accommodation of the new beads without introducing mechanical stress inside the structure. Next, we performed one more short run of 120 τ MD with the normal thermostat to accommodate the monomers further. The resulting conformation was then employed as the input configuration for the full MD production run.
It is important to stress that the number of 10nm fiber monomers and their content in base pairs are not equivalent (see Table 1). In fact, the unfolding of a given quantity of 30nm fiber monomers produces a comparable much longer (in μm) 10nm fiber: upon unfolding and equilibration, the 3000 basepair content lumped in a 30nm monomer will now extend over a total contour length of 270 nm. The total length (in μm) of the chromosome changes upon decondensation, and a total amount of 2% of 10nm fiber only will now extend over *18% of the total contour length. For the largest content of 30% of 10nm fiber studied in this work, the corresponding ratio to the total chromosome contour length is *80%. Given the fact that, from the physical point of view, is the relative contour length of the blocks rather than the relative base pair content that influences the behaviour, we assume having covered a sufficient range of case studies.
In this work, we have studied the effects of the local unfolding of chromatin fiber in two different scenarios: 1. In the first one, the 30nm fiber beads selected for being decondensed into 10nm fiber beads were chosen at random positions along the chromatin fiber, with corresponding amounts summarized in Table 1. Fig 8A shows a snapshot of a model chromosome conformation with 8% of 10nm fiber content.
2. In the second scenario, they were chosen clustered along a single linear chromatin sequence.
Here, two prototypical situations were considered: either the sequence is centered around the 30nm monomer at the closest distance to the chromosome centre of mass, or the sequence is centered around the monomer at the farthest distance from the center of mass.
In this way, we were able to monitor the two very distinct situations where chromatin  unfolding is located inside the corresponding chromosome territory or at its periphery. A summary of the simulated systems is given in Table 1. Fig 8B shows a snapshot of a model chromosome conformation with 8% of 10nm fiber content localized close to the chromosome centre of mass.
For reference, we have also considered the two limiting cases where all monomers are either 30nm fiber or 10nm fiber like (Table 1).

Data analysis
Quantities discussed in this work as the mean-square distances, hR 2 (L)i, and the average contact frequencies, hp c (L)i, between chromosome loci are plotted as a function of the genomic distance, L. In order to avoid unphysical behavior arising from the ring closure condition, we considered contour lengths L 1/4 of the total contour length of the ring, or L 30 Mbp.
Possible numerical artifacts due to the presence of monomers with different degrees of resolution ("10nm fiber" monomers vs. "30nm fiber" monomers) were removed by averaging over all possible pairs of monomers at fixed genomic separations L and spatial resolution of the 10nm fiber: this was achieved by replacing all 30nm fiber monomers not already decondensed by the 27 equivalent 10nm monomers, as shown in Fig 7 (panels C and D). In this way, each chromosome conformation always contributes with the same number of monomers and genomic distances smaller than 3 kbp can be effectively sampled. It is important to stress that what is described here concerns only the final analysis of the data, and it is not implicated in the motion of the monomers during the MD runs which is always performed as described in previous sections.
Final values for hR 2 (L)i and hp c (L)i at large L's were obtained by averaging further over logspaced intervals centered at the corresponding L's. This procedure improves the accuracy by reducing considerably statistical fluctuations. Error bars were calculated accordingly.