Next Article in Journal
Design, Synthesis, and Characterization of Schiff Base Bond-Linked pH-Responsive Doxorubicin Prodrug Based on Functionalized mPEG-PCL for Targeted Cancer Therapy
Previous Article in Journal
A Mild Method for Surface-Grafting PEG Onto Segmented Poly(Ester-Urethane) Film with High Grafting Density for Biomedical Purpose
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Chromatin Loop Extrusion and Chromatin Unknotting

1
Center for Integrative Genomics, University of Lausanne, 1015 Lausanne, Switzerland
2
SIB Swiss Institute of Bioinformatics, 1015 Lausanne, Switzerland
3
Polymer Institute of the Slovak Academy of Sciences, 842 36 Bratislava, Slovakia
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Polymers 2018, 10(10), 1126; https://doi.org/10.3390/polym10101126
Submission received: 14 September 2018 / Revised: 3 October 2018 / Accepted: 8 October 2018 / Published: 11 October 2018

Abstract

:
It has been a puzzle how decondensed interphase chromosomes remain essentially unknotted. The natural expectation is that in the presence of type II DNA topoisomerases that permit passages of double-stranded DNA regions through each other, all chromosomes should reach the state of topological equilibrium. The topological equilibrium in highly crowded interphase chromosomes forming chromosome territories would result in formation of highly knotted chromatin fibres. However, Chromosome Conformation Capture (3C) methods revealed that the decay of contact probabilities with the genomic distance in interphase chromosomes is practically the same as in the crumpled globule state that is formed when long polymers condense without formation of any knots. To remove knots from highly crowded chromatin, one would need an active process that should not only provide the energy to move the system from the state of topological equilibrium but also guide topoisomerase-mediated passages in such a way that knots would be efficiently unknotted instead of making the knots even more complex. We perform coarse-grained molecular dynamics simulations of the process of chromatin loop extrusion involving knotted and catenated chromatin fibres to check whether chromatin loop extrusion may be involved in active unknotting of chromatin fibres. Our simulations show that the process of chromatin loop extrusion is ideally suited to actively unknot, decatenate and demix chromatin fibres in interphase chromosomes.

Graphical Abstract

1. Introduction

Type II DNA topoisomerases that pass DNA segments through each other are required to facilitate various DNA transactions and they are critically important for the segregation of freshly replicated DNA molecules [1,2]. However, DNA–DNA passages occurring in highly crowded DNA molecules in living cells inadvertently lead to formation of knotted DNA [3,4]. Long-lasting knots on DNA are undesirable as tightened DNA knots hinder RNA transcription [5] and lead to DNA breakage [6]. A “proofreading” mechanism is needed to actively remove accidentally introduced knots. There are strong indications that an active knot-removing mechanism operates in highly crowded nuclei of eukaryotic cells. Chromosome conformation capture studies revealed that chromatin fibres in interphase chromosomes have the scaling behaviour characteristic of a fractal globule, which is a condensed polymer state that is free of knots [7]. More recent simulation studies have shown, though, that also knotted chromatin fibres can have the scaling behaviour of the fractal globule [8]. However, to show fractal globule scaling the knotting level in knotted polymers had to be low and, in addition, any interfibre passages were not permitted [8]. Thus, these more recent modelling studies also supported the notion that chromatin fibres in chromosomes are far from reaching topological equilibrium and thus are presumably actively unknotted. Very recent experimental studies based on the determination of the contact maps in individual cells [9] permitted the authors to reconstruct 3D paths of chromatin fibres in individual chromosomes and address the question of whether chromosomes are knotted [10]. These studies confirmed the notion that knots are rare in interphase chromosomes of mammalian cells since, despite their large size, about 20% of chromosomes were unknotted [10]. In addition, only some of the detected knots were observed repeatedly in all independent reconstructions resulting from the same experimental data. Remaining knots were only seen in some of individual reconstructions and therefore could have resulted from errors in the reconstruction procedure of the 3D path taken by chromatin fibres in a given chromosome [10].
Independent, strong argument for the active mechanism of unknotting and decatenation is provided by the fact that individual chromosomes in interphase cells form clearly defined chromosome territories [11]. If chromosomes were permitted to reach topological equilibrium, one would expect that chromatin from all chromosomes would strongly intermingle with each other [12].
The paucity of knots in chromatin, in addition to avoiding problems created by tightened DNA knots, allows unimpeded decondensation of chromosomal regions where genes are expressed or where replication occurs [13,14]. However, active unknotting mechanisms provide significant conceptual challenges [15]. Topoisomerases involved in unknotting would need to be guided to act on those juxtapositions of chromatin fibres where a passage simplifies or removes a knot rather than creates a knot or makes an existing knot more complex [16]. Currently, it is unknown how DNA topoisomerases are guided to unknot chromatin fibres. However, numerous recent experimental studies concluded that interphase chromosomes are shaped by an active process of chromatin loop extrusion [17,18,19]. We therefore decided to check by simulations whether chromatin loop extrusion can be involved in chromatin unknotting. Here we show, using molecular dynamics simulations, that the chromatin loop extrusion process is ideally suited to unknot and decatenate crowded chromatin fibres, as well as to segregate chromosomal territories from each other.

2. Materials and Methods

Description of the model. Molecular dynamics simulations presented in Figure 1 and Figure 2 were performed using Extensible Simulation Package for Research on Soft Matter (ESPResSo) [20] and those, results of which are presented in Figure 3 and Figure 4, were performed using a general-purpose particle simulation toolkit HOOMD-blue [21,22,23]. More details about applied simulation methods can be found in our recent publications [24,25]. However, the obtained results, i.e., that chromatin loop extrusion drives unknotting, are general and will not depend on specific settings of such parameters as the persistence length or the effective diameter of modelled chromatin fibres. Chromatin chains were modelled as standard, beaded self-avoiding chains that were composed of 133 beads for the chains presented in Figure 1, and of 400 beads for the construct presented in Figure 2. In these coarse-grained simulations, each bead is assumed to have 10 nm diameter and corresponds to a portion of chromatin fibre containing 400 bp [26]. The large simulated system presented in Figure 4 consisted of 8 chains, each with the length corresponding to 4000 beads placed after each other and thus each chain represented 1.6 Mb large chromatin portions. Nonstandard features built into our models included active movements of cohesin handcuffs (Figure 1 and Figure 2) and progressive shifting of the points of close contacts between contacting chains (Figure 4). Cohesin handcuffs embracing modelled chromatin fibres were composed of 2 small rings (each with 7 beads). These rings were connected with each other using two additional bonds involving two neighbouring beads of each ring. The cohesin rings were advancing along embraced chains by a walking mechanism where bonds with the spring-like potential and the starting length of up to 2 (where 1 is the diameter of beads) were dynamically formed between three nonconsecutive beads in each cohesin ring with progressively shifting beads within embraced chromatin fibres. As the spring-like bonds were very strong (K = 20 kBT) and their rest length was set to 1, they were rapidly shrinking, which resulted in the advancement of cohesin rings along embraced chromatin fibres. Once these spring-like bonds have shrunk to about their rest length, they were replaced with extended bonds connecting cohesin rings with next-neighbour beads positioned further along modelled chromatin fibres. Due to the tightness of cohesin rings and their excluded volume potential, as well as the excluded volume potential of beads forming chromatin fibres, the advancing cohesin rings were effectively pushing forward all encountered entanglements on modelled chromatin fibres. These entanglements were, however, free to untangle when they were pushed towards topologically associating domain (TAD) borders, where short sections of chains without excluded volume mimicked the presence of Top2B at TAD borders [27]. In the simulated system composed of 8 large chains, the process of chromatin loop extrusion was modelled in a simpler way. Initially, a small loop was generated by establishing a bond with a strong spring-like potential (K = 100 kBT) between two nonconsecutive beads. Once the two beads connected with that spring-like bond were brought close together, the spring-like bond was moved to bring together the next pair of beads, thus increasing the size of the formed loop with each move.
Voronoi space tessellation. The extent of mutual intermingling of separate chains was evaluated by using the Voronoi space tessellation approach [28]. Using this approach, the sphere containing 8 chains with 4000 beads each is dissected creating a polyhedron around each of the polymer beads. These polyhedrons represent regions in which every point is closer to the centre of a given bead than to the centre of any other bead. Each of the polyhedrons is created with four or more facets with a given surface area. The total area of the Voronoi facets separating beads belonging to different chains represents the total interfacial area (Figure 4). The interfacial area increases with increasing extent of mutual intermingling between different chains. The tessellation was constructed in a nonperiodic space with a spherical wall. The coordinates of Voronoi vertices, indices, facets with respective areas of the interfaces were obtained by using VORO++ software [28]. The visualizations were made using POVRay.
Knots identification. To determine and present the knot types formed after topological equilibration of 8 large rings, we entered xyz coordinates of each individual ring into Knoto-ID [29] to simplify trajectories of formed knots, while maintaining their original topology.

3. Results and Discussion

It is now broadly accepted that topologically associating domains (TADs) in interphase chromosomes of mammalian cells are formed by chromatin loop extrusion mechanism [30,31]. During chromatin loop extrusion, cohesin rings load on chromatin fibres in such a way that they initially form local, small chromatin loops, which are then progressively enlarged by active sliding of cohesin rings along the looped chromatin fibre [30,31]. Until now, chromatin loop extrusion was only seen as a mechanism ensuring that border elements of TADs approach each other, resulting in formation of loop-like chromatin regions with increased frequency of internal contacts [17,18,19,32,33] and it was not considered that chromatin loop extrusion can be involved in active unknotting of chromatin fibres.
To investigate whether chromatin loop extrusion can be involved in chromatin unknotting, we performed coarse-grained molecular dynamics simulations of chromatin loop extrusion involving knotted chromatin fibres (see Section 2 for more detailed description of the simulation procedure). In our proof of principle simulations, we started with knotted chromatin fibres forming closed loops. The closure of loops is justified by the fact that TADs spend a significant fraction of time as closed chromatin loops [17]. In addition, loop closure prevents the introduced knots from escaping. Figure 1 shows snapshots from a simulation of a loop extrusion occurring in knotted chromatin fibre that approximates a loop-forming TAD [17]. The simulated chromatin fibre forms a trefoil knot, which is the most frequently observed knot forming in chromatin rings in vivo [3]. Simulated chromatin loops contain a region mimicking functional elements of the TAD border. There are beads mimicking the presence of bound CTCF proteins, which stop the progression of cohesin rings [17,18,19]. There is also a short region with reduced excluded volume potential that permits passages of other portions of modelled chromatin fibres through that region. Having such a region in modelled chromatin fibres allows us to mimic the action of DNA Top2B that is known to be bound near CTCF and ensures that there is the possibility of chromatin fibres passing through each other [27]. Although type II DNA topoisomerase can permit intersegmental chromatin passages occurring also at other places, the majority of Topo II-mediated passages occur at borders of TADs [34] and that is why we introduced the possibility of Topo II-mediated passages occurring there. After a short thermal equilibration, the process of loop extrusion is started by progressively threading chromatin fibres through the cohesin handcuffs (green-coloured in Figure 1) and thus progressively enlarging the loop spanned by the cohesin handcuff (Figure 1b). Simulations show that as the loop extrusion progresses, cohesin handcuffs push before them the entanglements forming a knot and, therefore, the extruded part (brick-coloured), which is behind cohesin handcuffs, is free from knots, whereas the remaining part (blue) contains the original knot. As the loop extrusion proceeds, the knot gets pushed towards the TAD border, where it is unknotted by a passage involving the region with a reduced excluded volume potential (Figure 1c,d). In our simulations, we initiated chromatin loop extrusion at a position that is equally distant from the left and the right TAD border. In a biological setting, this may be a rare situation. However, when loop extrusion starts from any other site and the cohesin ring that first reaches the TAD border is stopped by bound CTCF, the second cohesin ring still progresses [19]; therefore, the end result is the same, i.e., the knot is pushed towards the TAD border and is unknotted there. Simulations aimed to reproduce experimental contact maps [18,19] and experiments investigating the dynamics of cohesin and CTCF in vivo [35] indicate that chromatin loop extrusion occurs repeatedly, i.e., cohesin handcuffs reach the border of TADs and stay there for some time but eventually dissociate and can start a new round. Therefore, all knots that may spontaneously form in the extruded chromatin loop as a result of Topo II action assisting transcription or replication will be again unknotted by the next round of chromatin loop extrusion.
Next, we wanted to test what happens when a more complex knot forms on a chromatin fibre. Since complex torus knots can be formed as products of site-specific recombination [36], we investigated how chromatin loop extrusion could unknot such complex knots. Figure 1e–h shows loop extrusion occurring in a TAD containing a complex torus knot with 14 crossings. As in the case of the trefoil knot, this more complicated knot gets also unknotted. This indicates that, in principle, chromatin loop extrusion is capable of promoting unknotting of even very complex knots.
Simulations presented in Figure 1 were of simplified systems where individual TADs were modelled as independent chromatin loops. In reality, TADs are positioned one after the other along the chromosome and are separated by short “spacers” [17]. In addition, chromatin knots do not need to be very local and confined within individual TADs, but can be larger, involving chromatin regions larger than one TAD. To model such a situation, we constructed a system where a chromatin portion containing three TADs was closed into a circle forming a large trefoil knot. Closure into a circle is needed to maintain the knot in this relatively short portion of modelled chromosome. In real chromosomes with thousands of TADs, closure would not be needed to keep the knots from slipping over the ends. Figure 2 shows what happens in such simulations where chromatin loop extrusion occurs in each of the constituent TADs. In order to better visualize the process in these proof-of-principle simulations, chromatin fibres were not equilibrated before chromatin loop extrusion has started. This approach allows us to initially maintain the quasi-symmetric configuration of modelled chromatin portion with three TADs, which is a remnant of the starting configuration. This setting of the simulation makes it easier to see that entanglements are excluded from extruded portions of TADs and accumulate in the inter-TAD portions, which are known to be sites where Top2B is located [27]. As already mentioned, the activity of Top2B is accounted for in our simulations by placing short regions with reduced excluded volume potential. As chromatin loop extrusion progresses, the entanglements of knots concentrate in the vicinity of regions with reduced self-avoidance and unknotting occurs at these regions (Figure 2d). In the simulation run presented in Figure 2, the three loop extrusions were proceeding synchronously. However, the same effect would be achieved when the loop extrusions were not synchronized but ended up in all loops being extruded. For unknotting, it is important that knot gets concentrated in a short nonextruded part between the extruded loops, as shown in Figure 2c. For this to happen, one can have the first loop extrusion going to completion, then the second, and then the third.
The simulations shown in Figure 1 and Figure 2 were of relatively simple systems where simulated TADs were not affected by crowding resulting from high concentration of chromatin fibres in chromosome territories [37]. To account for the crowding and the large length/diameter ratio of chromatin fibres forming TADs, we performed simulations of highly concentrated long chromatin rings undergoing loop extrusion. Our simulated system was composed of 8 chromatin rings having each size of about 1.6 Mb. Individual rings thus had sizes of large TADs. The concentration of modelled chromatin fibres was set to 30% by adequately adjusting the radius of the confining sphere where these 8 chromatin rings were placed; 30% volume occupation corresponds to the concentration of chromatin in nuclei of eukaryotic cells [38]. Subsequently, the system with 8 confined rings was topologically equilibrated by switching off the excluded volume potential and thus permitting thermally fluctuating chains to freely pass through themselves and through each other. After that step, we introduced excluded volume potential acting within and between all 8 chains. At the end of this procedure, each modelled chromatin ring was highly knotted (Figure 3) and each ring was also catenated with at least 3 other rings. The step of producing topological equilibrium is only to test what would happen if one could start chromatin loop extrusion from topologically equilibrated situation. As already mentioned, Hi-C studies established that chromatin in interphase cells is essentially unknotted.
Starting from this topologically very complex system composed of 8 knotted and catenated with each other rings that are in addition crowded (see Figure 4a), we initiated simultaneous loop extrusion in each of the 8 rings. In each ring, the point where the loop extrusion started was the most distant from the region with reduced extruded volume that was used to mimic the presence of Top2B at borders of TADs [27]. Panels a and b in Figure 4 show simulated system just before the loop extrusion has started (a) and just after the loop extrusion was completed (b). Each of the simulated rings is presented in a different color. Importantly, although the system started as a mixture of highly knotted and catenated rings, by the end of loop extrusion, each ring was both unknotted and decatenated from all other rings. By comparing panel a and b, one gets an impression that loop extrusion process demixes individual rings from each other and segregates territories occupied by each of the rings. To quantitatively measure the extent of demixing, we have used a Voronoi space tessellation approach [28] to subdivide the volume of the sphere containing the 8 chains into polyhedra enclosing each bead forming the 8 modelled chromatin rings (see methods section for more detail). We subsequently calculated the area of interface between all Voronoi polyhedra-containing beads belonging to different chains. Figure 4c shows the surface of interface between Voronoi polyhedra centred at beads of one individual chain and these centred at beads of other chains. We call these types of surfaces Voronoi envelopes. Figure 4d,e shows the total surface of all Voronoi envelopes for the system before and after chromatin loop extrusion. The graph in Figure 4f shows how the interfacial area between all Voronoi envelopes decreases during chromatin loop extrusion. Therefore, we can conclude that chromatin loop extrusion not only unknots and decatenates individual rings but also segregates modelled chromatin rings and decreases their intermingling.
Experimental studies of chromosomal territories revealed their limited intermingling [11,39]. It has been a question of how this limited intermingling is achieved in the presence of DNA topoisomerases [12]. Uncontrolled action of DNA topoisomerases is expected to result in a highly intermingled state of chromatin fibres from all decondensed chromosomes during the interphase stage of a cell cycle [7]. We showed here that chromatin loop extrusion is likely to be the driving force for unknotting, decatenation, and segregation of chromatin fibres forming different chromosomes. It is important to add here that unknotting, decatenation, and segregation investigated here just depend on the progressive movement of cohesin rings towards TAD borders. That movement might be driven by still undemonstrated translocation activity of cohesin [31] or by transcription-induced supercoiling [40].
In our simulations presented in Figure 4, there are no knots remaining in the modelled topological domain after chromatin loop extrusion is completed. However, this is the consequence of the fact that, in our model, knots could not form behind advancing cohesin rings as there were no possibilities of intersegmental passages there. In vivo, topoisomerases would be expected to act also within TAD and not only at their borders and this would result in inadvertent formation of knots. In the presence of repeated rounds of chromatin loop extrusion, an equilibrium of knotting would form where chromatin crowding and increasing number of sites permitting intersegmental passages would augment knotting, whereas increasing frequency of loop extrusion cycles would decrease knotting.
It needs to be stressed here that the mechanism of chromatin unknotting, decatenation, and demixing by chromatin loop extrusion is very much different from the mechanism of chromatin condensation needed for the formation of mitotic chromosomes. Chromatin condensation mediated by the action of condensin helps to disentangle sister chromatids from each other but this is due to effective shortening of chromatids changing them from long, thin, and floppy filaments that easily entangle with each other to short, thick, and rigid rods that are unable to entangle with each other [26,41]. In contrast to condensin-mediated condensation of chromosomes, cohesin-mediated chromatin loop extrusion does not induce condensation of loops that passed through cohesin rings but just moves chromatin fibres through cohesin rings.
It is a popular saying that one cannot unscramble an egg. Apparently, the chromatin loop extrusion process can very well unscramble scrambled chromatin fibres forming neighbouring chromosome territories. Chromatin loop extrusion can also maintain chromosomes as essentially unknotted and uncatenated with each other even if some knotting and catenation events may occur between sequential rounds of chromatin loop extrusion. The rare knots observed in single cell Hi-C [10] are presumably results of knotting events that occurred after the latest round of chromatin loop extrusion and before cell fixation for Hi-C.
After this work was completed, a very similar idea of chromatin loop extrusion involvement in chromatin unknotting and decatenation was proposed in a preprint by Orlandini et al. [42].

Author Contributions

A.S. conceived the study. D.R. and F.B. performed simulations and implemented new features into molecular dynamics simulation packages. D.G. analyzed the topology of detected knots. A.S. wrote the manuscript with inputs from all authors.

Funding

This research was funded by the Leverhulme Trust, grant number RP2013-K-017 and by the Swiss National Science Foundation, grant number 31003A_166684. The APC was funded by Swiss National Science Foundation.

Acknowledgments

We thank Dorothy Buck and Sean D. Colloms for their comments and suggestions concerning this work.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Bates, A.D.; Maxwell, A. DNA Topology; Oxford University Press: Oxford, UK, 2005; p. 198. [Google Scholar]
  2. Bjorkegren, C.; Baranello, L. DNA Supercoiling, Topoisomerases, and Cohesin: Partners in Regulating Chromatin Architecture? Int. J. Mol. Sci. 2018, 19. [Google Scholar] [CrossRef] [PubMed]
  3. Valdes, A.; Segura, J.; Dyson, S.; Martinez-Garcia, B.; Roca, J. DNA knots occur in intracellular chromatin. Nucleic Acids Res. 2018, 46, 650–660. [Google Scholar] [CrossRef] [PubMed]
  4. O’Donnol, D.; Stasiak, A.; Buck, D. Two convergent pathways of DNA knotting in replicating DNA molecules as revealed by theta-curve analysis. Nucleic Acids Res. 2018. [Google Scholar] [CrossRef] [PubMed]
  5. Portugal, J.; Rodriguez-Campos, A. T7 RNA polymerase cannot transcribe through a highly knotted DNA template. Nucleic Acids Res. 1996, 24, 4890–4894. [Google Scholar] [CrossRef] [PubMed]
  6. Deibler, R.W.; Mann, J.K.; Sumners, D.W.; Zechiedrich, L. Hin-mediated DNA knotting and recombining promote replicon dysfunction and mutation. BMC Mol. Biol. 2007, 8, 44. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Lieberman-Aiden, E.; van Berkum, N.L.; Williams, L.; Imakaev, M.; Ragoczy, T.; Telling, A.; Amit, I.; Lajoie, B.R.; Sabo, P.J.; Dorschner, M.O.; et al. Comprehensive mapping of long-range interactions reveals folding principles of the human genome. Science 2009, 326, 289–293. [Google Scholar] [CrossRef] [PubMed]
  8. Arsuaga, J.; Jayasinghe, R.G.; Scharein, R.G.; Segal, M.R.; Stolz, R.H.; Vazquez, M. Current theoretical models fail to predict the topological complexity of the human genome. Front. Mol. Biosci. 2015, 2, 48. [Google Scholar] [CrossRef] [PubMed]
  9. Stevens, T.J.; Lando, D.; Basu, S.; Atkinson, L.P.; Cao, Y.; Lee, S.F.; Leeb, M.; Wohlfahrt, K.J.; Boucher, W.; O’Shaughnessy-Kirwan, A.; et al. 3D structures of individual mammalian genomes studied by single-cell Hi-C. Nature 2017, 544, 59–64. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Siebert, J.; Kivel, A.; Atkinson, L.; Stevens, T.; Laue, E.; Virnau, P. Are there knots in chromosomes? Polymers 2017, 9, 317. [Google Scholar] [CrossRef]
  11. Bolzer, A.; Kreth, G.; Solovei, I.; Koehler, D.; Saracoglu, K.; Fauth, C.; Muller, S.; Eils, R.; Cremer, C.; Speicher, M.R.; et al. Three-dimensional maps of all chromosomes in human male fibroblast nuclei and prometaphase rosettes. PLoS Biol. 2005, 3, e157. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Dorier, J.; Stasiak, A. Topological origins of chromosomal territories. Nucleic Acids Res. 2009, 37, 6316–6322. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Grosberg, A.; Rabin, Y.; Havlin, S.; Neer, A. Crumpled globule model of the three-dimensional structure of DNA. Eur. Lett. 1993, 23, 373. [Google Scholar] [CrossRef]
  14. Mirny, L.A. The fractal globule as a model of chromatin architecture in the cell. Chromosome Res. 2011, 19, 37–51. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Yan, J.; Magnasco, M.O.; Marko, J.F. A kinetic proofreading mechanism for disentanglement of DNA by topoisomerases. Nature 1999, 401, 932–935. [Google Scholar] [CrossRef] [PubMed]
  16. Buck, G.R.; Zechiedrich, E.L. DNA disentangling by type-2 topoisomerases. J. Mol. Biol. 2004, 340, 933–939. [Google Scholar] [CrossRef] [PubMed]
  17. Rao, S.S.; Huntley, M.H.; Durand, N.C.; Stamenova, E.K.; Bochkov, I.D.; Robinson, J.T.; Sanborn, A.L.; Machol, I.; Omer, A.D.; Lander, E.S.; et al. A 3D map of the human genome at kilobase resolution reveals principles of chromatin looping. Cell 2014, 159, 1665–1680. [Google Scholar] [CrossRef] [PubMed]
  18. Sanborn, A.L.; Rao, S.S.; Huang, S.C.; Durand, N.C.; Huntley, M.H.; Jewett, A.I.; Bochkov, I.D.; Chinnappan, D.; Cutkosky, A.; Li, J.; et al. Chromatin extrusion explains key features of loop and domain formation in wild-type and engineered genomes. Proc. Natl. Acad. Sci. USA 2015, 112, E6456–E6465. [Google Scholar] [CrossRef] [PubMed]
  19. Fudenberg, G.; Imakaev, M.; Lu, C.; Goloborodko, A.; Abdennur, N.; Mirny, L.A. Formation of chromosomal domains by loop extrusion. Cell Rep. 2016, 15, 2038–2049. [Google Scholar] [CrossRef] [PubMed]
  20. Limbach, H.J.; Arnold, A.; Mann, B.A.; Holm, C. ESPResSo—An extensible simulation package for research on soft matter systems. Comp. Phys. Comm. 2006, 174, 704–727. [Google Scholar] [CrossRef]
  21. Anderson, J.A.; Lorenz, C.D.; Travesset, A. General purpose molecular dynamics simulations fully implemented on graphics processing units. J. Comput. Phys. 2008, 227, 5342–5359. [Google Scholar] [CrossRef] [Green Version]
  22. Glaser, J.; Nguyen, T.D.; Anderson, J.A.; Liu, P.; Spiga, F.; Millan, J.A.; Morse, D.C.; Glotzer, S.C. Strong scaling of general-purpose molecular dynamics simulations on GPUs. Comput. Phys. Commun. 2015, 192, 97–107. [Google Scholar] [CrossRef] [Green Version]
  23. Reith, D.; Mirny, L.; Virnau, P. GPU based molecular dynamics simulations of polymer rings in concentrated solution: Structure and scaling. Prog. Theor. Phys. Suppl. 2011, 191, 135–145. [Google Scholar] [CrossRef]
  24. Racko, D.; Benedetti, F.; Dorier, J.; Burnier, Y.; Stasiak, A. Generation of supercoils in nicked and gapped DNA drives DNA unknotting and postreplicative decatenation. Nucleic Acids Res. 2015, 43, 7229–7236. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Racko, D.; Benedetti, F.; Dorier, J.; Burnier, Y.; Stasiak, A. Molecular Dynamics Simulation of Supercoiled, Knotted, and Catenated DNA Molecules, Including Modeling of Action of DNA Gyrase. Methods Mol. Biol. 2017, 1624, 339–372. [Google Scholar] [CrossRef] [PubMed]
  26. Goloborodko, A.; Imakaev, M.V.; Marko, J.F.; Mirny, L. Compaction and segregation of sister chromatids via active loop extrusion. eLife 2016, 5. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Uuskula-Reimand, L.; Hou, H.; Samavarchi-Tehrani, P.; Rudan, M.V.; Liang, M.; Medina-Rivera, A.; Mohammed, H.; Schmidt, D.; Schwalie, P.; Young, E.J.; et al. Topoisomerase II beta interacts with cohesin and CTCF at topological domain borders. Genome Biol. 2016, 17, 182. [Google Scholar] [CrossRef] [PubMed]
  28. Rycroft, C.H.; Grest, G.S.; Landry, J.W.; Bazant, M.Z. Analysis of granular flow in a pebble-bed nuclear reactor. Phys. Rev. E Stat. Nonlinear Soft Matter Phys. 2006, 74, 021306. [Google Scholar] [CrossRef] [PubMed]
  29. Dorier, J.; Goundaroulis, D.; Benedetti, F.; Stasiak, A. Knoto-ID: A tool to study the entanglement of open protein chains using the concept of knotoids. Bioinformatics 2018, 34, 3402–3404. [Google Scholar] [CrossRef] [PubMed]
  30. Schwarzer, W.; Abdennur, N.; Goloborodko, A.; Pekowska, A.; Fudenberg, G.; Loe-Mie, Y.; Fonseca, N.A.; Huber, W.; Haering, C.H.; Mirny, L.; et al. Two independent modes of chromatin organization revealed by cohesin removal. Nature 2017, 551, 51–56. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Fudenberg, G.; Abdennur, N.; Imakaev, M.; Goloborodko, A.; Mirny, L.A. Emerging evidence of chromosome folding by loop extrusion. Cold Spring Harb. Symp. Quant. Biol. 2017, 82, 45–55. [Google Scholar] [CrossRef] [PubMed]
  32. Dixon, J.R.; Selvaraj, S.; Yue, F.; Kim, A.; Li, Y.; Shen, Y.; Hu, M.; Liu, J.S.; Ren, B. Topological domains in mammalian genomes identified by analysis of chromatin interactions. Nature 2012, 485, 376–380. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Nora, E.P.; Lajoie, B.R.; Schulz, E.G.; Giorgetti, L.; Okamoto, I.; Servant, N.; Piolot, T.; van Berkum, N.L.; Meisig, J.; Sedat, J.; et al. Spatial partitioning of the regulatory landscape of the X-inactivation centre. Nature 2012, 485, 381–385. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Canela, A.; Maman, Y.; Jung, S.; Wong, N.; Callen, E.; Day, A.; Kieffer-Kwon, K.R.; Pekowska, A.; Zhang, H.; Rao, S.S.P.; et al. Genome Organization Drives Chromosome Fragility. Cell 2017, 170, 507–521.e518. [Google Scholar] [CrossRef] [PubMed]
  35. Hansen, A.S.; Pustova, I.; Cattoglio, C.; Tjian, R.; Darzacq, X. CTCF and cohesin regulate chromatin loop stability with distinct dynamics. eLife 2017, 6. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Spengler, S.J.; Stasiak, A.; Cozzarelli, N.R. The Stereostructure of Knots and Catenanes Produced by Phage-Lambda Integrative Recombination—Implications for Mechanism and DNA-Structure. Cell 1985, 42, 325–334. [Google Scholar] [CrossRef]
  37. Meaburn, K.J.; Misteli, T. Cell biology: Chromosome territories. Nature 2007, 445, 379–781. [Google Scholar] [CrossRef] [PubMed]
  38. Sugawara, T.; Kimura, A. Physical properties of the chromosomes and implications for development. Dev. Growth Differ. 2017, 59, 405–414. [Google Scholar] [CrossRef] [PubMed]
  39. Branco, M.R.; Pombo, A. Intermingling of chromosome territories in interphase suggests role in translocations and transcription-dependent associations. PLoS Biol. 2006, 4, e138. [Google Scholar] [CrossRef] [PubMed]
  40. Racko, D.; Benedetti, F.; Dorier, J.; Stasiak, A. Transcription-induced supercoiling as the driving force of chromatin loop extrusion during formation of TADs in interphase chromosomes. Nucleic Acids Res. 2018, 46, 1648–1660. [Google Scholar] [CrossRef] [PubMed]
  41. Alipour, E.; Marko, J.F. Self-organization of domain structures by DNA-loop-extruding enzymes. Nucleic Acids Res. 2012, 40, 11202–11212. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Orlandini, E.; Marenduzzo, D.; Michieletto, D. Synergy of SMC and topoisomerase creates a universal pathway to simplify genome topology. arXiv, 2018; arXiv:1809.01267. [Google Scholar]
Figure 1. Chromatin loop extrusion unknots individual loop-forming TADs modelled as chromatin rings. (ad) Unknotting of a simple trefoil knot. (eh) Unknotting of a complex knot with 14 crossings. (c,d) show highly compressed knots being pushed towards a region with reduced excluded volume potential, which in our simulations mimic the action of Top2B, known to be present near TADs borders. Inset in (c) shows a close-up of a strongly confined knot, which is about to get unknotted by an intersegmental passage involving a region with reduced excluded volume potential. (d,h) show that at the end of simulations the knots are already unknotted. In our simulations, cohesin rings (green) are initially arbitrarily placed opposite to borders of modelled TADs. Once simulations are started, stacked rings of cohesin take the form of a quasi-planar handcuff and actively translocate in that form. Borders of TADs are recognizable here by yellow coloured beads that in our simulations stop the progression of cohesin rings and thus mimic the biological action of CTCF proteins located at borders of TADs. TAD: topologically associating domain.
Figure 1. Chromatin loop extrusion unknots individual loop-forming TADs modelled as chromatin rings. (ad) Unknotting of a simple trefoil knot. (eh) Unknotting of a complex knot with 14 crossings. (c,d) show highly compressed knots being pushed towards a region with reduced excluded volume potential, which in our simulations mimic the action of Top2B, known to be present near TADs borders. Inset in (c) shows a close-up of a strongly confined knot, which is about to get unknotted by an intersegmental passage involving a region with reduced excluded volume potential. (d,h) show that at the end of simulations the knots are already unknotted. In our simulations, cohesin rings (green) are initially arbitrarily placed opposite to borders of modelled TADs. Once simulations are started, stacked rings of cohesin take the form of a quasi-planar handcuff and actively translocate in that form. Borders of TADs are recognizable here by yellow coloured beads that in our simulations stop the progression of cohesin rings and thus mimic the biological action of CTCF proteins located at borders of TADs. TAD: topologically associating domain.
Polymers 10 01126 g001
Figure 2. Chromatin loop extrusion unknots a delocalized trefoil knot spanning a modelled chromosome portion containing three TADs. Three TADs, each with two cohesin rings forming handcuffs (green), beads mimicking action of CTCF (yellow) and TopII sites (green semi-transparent) are modelled as presented in Figure 1. (a) The starting configuration shows that the formed trefoil knot is not localized within an individual TAD, but diffuses over three TADs. (bd) Chromatin loop extrusion proceeds in all three TADs. The knot becomes excluded from the extruded chromatin portions and is concentrated in the chromatin portions with TopII sites of action, which are modelled as regions with reduced excluded volume potential. As the knot becomes more confined, it is more likely to interact with topoisomerases localized outside of TADs and is eventually unknotted there. Insets show close-ups of regions near TADs borders. For better visibility, the diameter of modelled chromatin fibres and cohesin handcuffs is decreased. In the inset shown in (c), the arrow indicates an unknotting passage involving a region with reduced excluded volume potential (semitransparent) and a chromatin portion that did not pass yet through cohesin rings (blue). In the inset shown in (d), the extrusion is nearly finished as cohesin rings are approaching sites with bound CTCF (yellow). The portion of chromatin, which will not be extruded further, is already free from knots.
Figure 2. Chromatin loop extrusion unknots a delocalized trefoil knot spanning a modelled chromosome portion containing three TADs. Three TADs, each with two cohesin rings forming handcuffs (green), beads mimicking action of CTCF (yellow) and TopII sites (green semi-transparent) are modelled as presented in Figure 1. (a) The starting configuration shows that the formed trefoil knot is not localized within an individual TAD, but diffuses over three TADs. (bd) Chromatin loop extrusion proceeds in all three TADs. The knot becomes excluded from the extruded chromatin portions and is concentrated in the chromatin portions with TopII sites of action, which are modelled as regions with reduced excluded volume potential. As the knot becomes more confined, it is more likely to interact with topoisomerases localized outside of TADs and is eventually unknotted there. Insets show close-ups of regions near TADs borders. For better visibility, the diameter of modelled chromatin fibres and cohesin handcuffs is decreased. In the inset shown in (c), the arrow indicates an unknotting passage involving a region with reduced excluded volume potential (semitransparent) and a chromatin portion that did not pass yet through cohesin rings (blue). In the inset shown in (d), the extrusion is nearly finished as cohesin rings are approaching sites with bound CTCF (yellow). The portion of chromatin, which will not be extruded further, is already free from knots.
Polymers 10 01126 g002
Figure 3. Topology of knotted rings formed upon topological equilibration of 8 modelled chromatin rings. Knot types of all 8 rings forming a system presented in Figure 4. Alexander–Briggs notation of knots uses two numbers to indicate the knot type. The first number written in normal font indicates the minimal number of crossings of a given knot and the second number written as a subscript indicates the tabular position of this knot among the knots with the same minimal crossing number. So for example, 911 indicates a knot that in tables of knots is placed at 11th position among knots which minimal crossing number is 9. The # sign indicates that a knot is formed by a composition or merging together of simpler (prime) knots that are listed as its components. For prime knots with a large number of minimal crossings, e.g., 12, as is the case of ring 4, the knots are additionally divided into these with alternating pattern of crossings (indicated as a) and with nonalternating pattern of crossings (indicated as n). Standard tables of knots do not include knots with more than 12 crossings and we indicate these knots as unclassified.
Figure 3. Topology of knotted rings formed upon topological equilibration of 8 modelled chromatin rings. Knot types of all 8 rings forming a system presented in Figure 4. Alexander–Briggs notation of knots uses two numbers to indicate the knot type. The first number written in normal font indicates the minimal number of crossings of a given knot and the second number written as a subscript indicates the tabular position of this knot among the knots with the same minimal crossing number. So for example, 911 indicates a knot that in tables of knots is placed at 11th position among knots which minimal crossing number is 9. The # sign indicates that a knot is formed by a composition or merging together of simpler (prime) knots that are listed as its components. For prime knots with a large number of minimal crossings, e.g., 12, as is the case of ring 4, the knots are additionally divided into these with alternating pattern of crossings (indicated as a) and with nonalternating pattern of crossings (indicated as n). Standard tables of knots do not include knots with more than 12 crossings and we indicate these knots as unclassified.
Polymers 10 01126 g003
Figure 4. Chromatin loop extrusion demixes topologically equilibrated chromatin loops. (a) A snapshot of topologically equilibrated system composed of 8 large chromatin circles (each shown in a different colour) that are knotted and catenated with each other. (b) A snapshot showing how the simulated system changes after each of the rings has undergone chromatin loop extrusion. (c) One of the 8 chromatin circles is shown together with its Voronoi envelope marking the interface of that chromatin circle with other circles. (d,e) Voronoi envelopes marking interfacial area between different chromatin circles before and after chromatin loop, respectively. (f) As chromatin loop extrusion progresses, the interfacial area between all modelled loops decreases. The snapshots present the confined chains together with their Voronoi envelopes at the moments when chromatin loop extrusion has started and finished, respectively. Red, encircled points indicate the interfacial area measures for the shown snapshots.
Figure 4. Chromatin loop extrusion demixes topologically equilibrated chromatin loops. (a) A snapshot of topologically equilibrated system composed of 8 large chromatin circles (each shown in a different colour) that are knotted and catenated with each other. (b) A snapshot showing how the simulated system changes after each of the rings has undergone chromatin loop extrusion. (c) One of the 8 chromatin circles is shown together with its Voronoi envelope marking the interface of that chromatin circle with other circles. (d,e) Voronoi envelopes marking interfacial area between different chromatin circles before and after chromatin loop, respectively. (f) As chromatin loop extrusion progresses, the interfacial area between all modelled loops decreases. The snapshots present the confined chains together with their Voronoi envelopes at the moments when chromatin loop extrusion has started and finished, respectively. Red, encircled points indicate the interfacial area measures for the shown snapshots.
Polymers 10 01126 g004

Share and Cite

MDPI and ACS Style

Racko, D.; Benedetti, F.; Goundaroulis, D.; Stasiak, A. Chromatin Loop Extrusion and Chromatin Unknotting. Polymers 2018, 10, 1126. https://doi.org/10.3390/polym10101126

AMA Style

Racko D, Benedetti F, Goundaroulis D, Stasiak A. Chromatin Loop Extrusion and Chromatin Unknotting. Polymers. 2018; 10(10):1126. https://doi.org/10.3390/polym10101126

Chicago/Turabian Style

Racko, Dusan, Fabrizio Benedetti, Dimos Goundaroulis, and Andrzej Stasiak. 2018. "Chromatin Loop Extrusion and Chromatin Unknotting" Polymers 10, no. 10: 1126. https://doi.org/10.3390/polym10101126

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop