Kinetics of chromosome condensation in the presence of topoisomerases: a phantom chain model

Bridging-mediated compaction of mitotic chromosomes

Within living cells, chromosome shapes undergo a striking morphological transition, from loose and uncondensed fibers during interphase to compacted and cylindrical structures during mitosis. ATP driven loop extrusion performed by a specialized protein complex, condensin, has recently emerged as a key driver of this transition. However, while this mechanism can successfully recapitulate the compaction of chromatids during the early stages of mitosis, it cannot capture structures observed after prophase. Here we hypothesize that a condensin bridging activity plays an additional important role, and review evidence - obtained largely through molecular dynamics simulations - that, in combination with loop extrusion, it can generate compact metaphase cylinders. Additionally, the resulting model qualitatively explains the unusual elastic properties of mitotic chromosomes observed in micromanipulation experiments and provides insights into the role of condensins in the formation of abnormal chromosome structures associated with common fragile sites.

The chromosome folding problem and how cells solve it

Every cell must solve the problem of how to fold its genome. We describe how the folded state of chromosomes is the result of the combined activity of multiple conserved mechanisms. Homotypic affinity-driven interactions lead to spatial partitioning of active and inactive loci. Molecular motors fold chromosomes through loop extrusion. Topological features such as supercoiling and entanglements contribute to chromosome folding and its dynamics, and tethering loci to sub-nuclear structures adds additional constraints. Dramatically diverse chromosome conformations observed throughout the cell cycle and across the tree of life can be explained through differential regulation and implementation of these basic mechanisms. We propose that the first functions of chromosome folding are to mediate genome replication, compaction, and segregation and that mechanisms of folding have subsequently been co-opted for other roles, including long-range gene regulation, in different conditions, cell types, and species.

Mitotic chromosomes are self-entangled and disentangle through a topoisomerase-II-dependent two-stage exit from mitosis

The topological state of chromosomes determines their mechanical properties, dynamics, and function. Recent work indicated that interphase chromosomes are largely free of entanglements. Here, we use Hi-C, polymer simulations, and multi-contact 3C and find that, by contrast, mitotic chromosomes are self-entangled. We explore how a mitotic self-entangled state is converted into an unentangled interphase state during mitotic exit. Most mitotic entanglements are removed during anaphase/telophase, with remaining ones removed during early G1, in a topoisomerase-II-dependent process. Polymer models suggest a two-stage disentanglement pathway: first, decondensation of mitotic chromosomes with remaining condensin loops produces entropic forces that bias topoisomerase II activity toward decatenation. At the second stage, the loops are released, and the formation of new entanglements is prevented by lower topoisomerase II activity, allowing the establishment of unentangled and territorial G1 chromosomes. When mitotic entanglements are not removed in experiments and models, a normal interphase state cannot be acquired.

Homologous chromosome recognition via nonspecific interactions

In many organisms, most notably Drosophila, homologous chromosomes in somatic cells associate with each other, a phenomenon known as somatic homolog pairing. Unlike in meiosis, where homology is read out at the level of DNA sequence complementarity, somatic homolog pairing takes place without double strand breaks or strand invasion, thus requiring some other mechanism for homologs to recognize each other. Several studies have suggested a "specific button" model, in which a series of distinct regions in the genome, known as buttons, can associate with each other, presumably mediated by different proteins that bind to these different regions. Here we consider an alternative model, which we term the "button barcode" model, in which there is only one type of recognition site or adhesion button, present in many copies in the genome, each of which can associate with any of the others with equal affinity. An important component of this model is that the buttons are non-uniformly distributed, such that alignment of a chromosome with its correct homolog, compared with a non-homolog, is energetically favored; since to achieve nonhomologous alignment, chromosomes would be required to mechanically deform in order to bring their buttons into mutual register. We investigated several types of barcodes and examined their effect on pairing fidelity. We found that high fidelity homolog recognition can be achieved by arranging chromosome pairing buttons according to an actual industrial barcode used for warehouse sorting. By simulating randomly generated non-uniform button distributions, many highly effective button barcodes can be easily found, some of which achieve virtually perfect pairing fidelity. This model is consistent with existing literature on the effect of translocations of different sizes on homolog pairing. We conclude that a button barcode model can attain highly specific homolog recognition, comparable to that seen in actual cells undergoing somatic homolog pairing, without the need for specific interactions. This model may have implications for how meiotic pairing is achieved.

Topology-driven spatial organization of ring polymers under confinement

Entropic repulsion between DNA ring polymers under confinement is a key mechanism governing the spatial segregation of bacterial DNA before cell division. Here we establish that "internal" loops within a modified-ring polymer architecture enhance entropic repulsion between two overlapping polymers confined in a cylinder. Interestingly, they also induce entropy-driven spatial organization of polymer segments as seen in vivo. Here we design polymers of different architectures in our simulations by introducing a minimal number of cross-links between specific monomers along the ring-polymer contour. The cross-links are likely induced by various bridging proteins inside living cells. We investigate the segregation of two polymers with modified topologies confined in a cylinder, which initially had spatially overlapping configurations. This helps us to identify the architectures that lead to higher success rates of segregation. We also establish the mechanism that leads to localization of specific polymer segments. We use the blob model to provide a theoretical understanding of why certain architectures lead to enhanced entropic repulsive forces between the polymers. Lastly, we establish a correspondence between the organizational patterns of the chromosome of the C.crescentus bacterium and our results for a specifically designed polymer architecture. However, the principles outlined here pertaining to the organization of polymeric segments are applicable to both synthetic and biological polymers.

Modeling cell biological features of meiotic chromosome pairing to study interlock resolution

During meiosis, homologous chromosomes become associated side by side in a process known as homologous chromosome pairing. Pairing requires long range chromosome motion through a nucleus that is full of other chromosomes. It remains unclear how the cell manages to align each pair of chromosomes quickly while mitigating and resolving interlocks. Here, we use a coarse-grained molecular dynamics model to investigate how specific features of meiosis, including motor-driven telomere motion, nuclear envelope interactions, and increased nuclear size, affect the rate of pairing and the mitigation/resolution of interlocks. By creating in silico versions of three yeast strains and comparing the results of our model to experimental data, we find that a more distributed placement of pairing sites along the chromosome is necessary to replicate experimental findings. Active motion of the telomeric ends speeds up pairing only if binding sites are spread along the chromosome length. Adding a meiotic bouquet significantly speeds up pairing but does not significantly change the number of interlocks. An increase in nuclear size slows down pairing while greatly reducing the number of interlocks. Interestingly, active forces increase the number of interlocks, which raises the question: How do these interlocks resolve? Our model gives us detailed movies of interlock resolution events which we then analyze to build a step-by-step recipe for interlock resolution. In our model, interlocks must first translocate to the ends, where they are held in a quasi-stable state by a large number of paired sites on one side. To completely resolve an interlock, the telomeres of the involved chromosomes must come in close proximity so that the cooperativity of pairing coupled with random motion causes the telomeres to unwind. Together our results indicate that computational modeling of homolog pairing provides insight into the specific cell biological changes that occur during meiosis.

Simulation of Different Three-Dimensional Models of Whole Interphase Nuclei Compared to Experiments - A Consistent Scale-Bridging Simulation Framework for Genome Organization

The three-dimensional architecture of chromosomes, their arrangement, and dynamics within cell nuclei are still subject of debate. Obviously, the function of genomes-the storage, replication, and transcription of genetic information-has closely coevolved with this architecture and its dynamics, and hence are closely connected. In this work a scale-bridging framework investigates how of the 30 nm chromatin fibre organizes into chromosomes including their arrangement and morphology in the simulation of whole nuclei. Therefore, mainly two different topologies were simulated with corresponding parameter variations and comparing them to experiments: The Multi-Loop-Subcompartment (MLS) model, in which (stable) small loops form (stable) rosettes, connected by chromatin linkers, and the Random-Walk/Giant-Loop (RW/GL) model, in which large loops are attached to a flexible non-protein backbone, were simulated for various loop and linker sizes. The 30 nm chromatin fibre was modelled as a polymer chain with stretching, bending and excluded volume interactions. A spherical boundary potential simulated the confinement to nuclei with different radii. Simulated annealing and Brownian Dynamics methods were applied in a four-step decondensation procedure to generate from metaphase decondensated interphase configurations at thermodynamical equilibrium. Both the MLS and the RW/GL models form chromosome territories, with different morphologies: The MLS rosettes result in distinct subchromosomal domains visible in electron and confocal laser scanning microscopic images. In contrast, the big RW/GL loops lead to a mostly homogeneous chromatin distribution. Even small changes of the model parameters induced significant rearrangements of the chromatin morphology. The low overlap of chromosomes, arms, and subchromosomal domains observed in experiments agrees only with the MLS model. The chromatin density distribution in CLSM image stacks reveals a bimodal behaviour in agreement with recent experiments. Combination of these results with a variety of (spatial distance) measurements favour an MLS like model with loops and linkers of 63 to 126 kbp. The predicted large spaces between the chromatin fibres allow typically sized biological molecules to reach nearly every location in the nucleus by moderately obstructed diffusion and is in disagreement with the much simplified assumption that defined channels between territories for molecular transport as in the Interchromosomal Domain (ICD) hypothesis exist and are necessary for transport. All this is also in agreement with recent selective high-resolution chromosome interaction capture (T2C) experiments, the scaling behaviour of the DNA sequence, the dynamics of the chromatin fibre, the diffusion of molecules, and other measurements. Also all other chromosome topologies can in principle be excluded. In summary, polymer simulations of whole nuclei compared to experimental data not only clearly favour only a stable loop aggregate/rosette like genome architecture whose local topology is tightly connected to the global morphology and dynamics of the cell nucleus and hence can be used for understanding genome organization also in respect to diagnosis and treatment. This is in agreement with and also leads to a general novel framework of genome emergence, function, and evolution.

Computer simulations of melts of ring polymers with nonconserved topology: A dynamic Monte Carlo lattice model

We present computer simulations of a dynamic Monte Carlo algorithm for polymer chains on a fcc lattice which explicitly takes into account the possibility to overcome topological constraints by controlling the rate at which nearby polymer strands may cross through each other. By applying the method to systems of interacting ring polymers at melt conditions, we characterize their structure and dynamics by measuring, in particular, the amounts of knots and links which are formed during the relaxation process. In comparison with standard melts of unknotted and unconcatenated rings, our simulations demonstrate that the mechanism of strand crossing makes polymer dynamics faster provided the characteristic timescale of the process is smaller than the typical timescale for chain relaxation in the unperturbed state, in agreement with recent experiments employing solutions of DNA rings in the presence of the type II topoisomerase enzyme. In the opposite case of slow rates the melt is shown to become slower, and this prediction may be easily validated experimentally.

Chromosome disentanglement driven via optimal compaction of loop-extruded brush structures

Eukaryote cell division features a chromosome compaction-decompaction cycle that is synchronized with their physical and topological segregation. It has been proposed that lengthwise compaction of chromatin into mitotic chromosomes via loop extrusion underlies the compaction-segregation/resolution process. We analyze this disentanglement scheme via considering the chromosome to be a succession of DNA/chromatin loops-a polymer "brush"-where active extrusion of loops controls the brush structure. Given type-II DNA topoisomerase (Topo II)-catalyzed topology fluctuations, we find that interchromosome entanglements are minimized for a certain "optimal" loop that scales with the chromosome size. The optimal loop organization is in accord with experimental data across species, suggesting an important structural role of genomic loops in maintaining a less entangled genome. Application of the model to the interphase genome indicates that active loop extrusion can maintain a level of chromosome compaction with suppressed entanglements; the transition to the metaphase state requires higher lengthwise compaction and drives complete topological segregation. Optimized genomic loops may provide a means for evolutionary propagation of gene-expression patterns while simultaneously maintaining a disentangled genome. We also find that compact metaphase chromosomes have a densely packed core along their cylindrical axes that explains their observed mechanical stiffness. Our model connects chromosome structural reorganization to topological resolution through the cell cycle and highlights a mechanism of directing Topo II-mediated strand passage via loop extrusion-driven lengthwise compaction.

Modeling meiotic chromosome pairing: a tug of war between telomere forces and a pairing-based Brownian ratchet leads to increased pairing fidelity

Meiotic homolog pairing involves associations between homologous DNA regions scattered along the length of a chromosome. When homologs associate, they tend to do so by a processive zippering process, which apparently results from avidity effects. Using a computational model, we show that this avidity-driven processive zippering reduces the selectivity of pairing. When active random forces are applied to telomeres, this drop in selectivity is eliminated in a force-dependent manner. Further simulations suggest that active telomere forces are engaged in a tug-of-war against zippering, which can be interpreted as a Brownian ratchet with a stall force that depends on the dissociation constant of pairing. When perfectly homologous regions of high affinity compete with homeologous regions of lower affinity, the affinity difference can be amplified through this tug of war effect provided the telomere force acts in a range that is strong enough to oppose zippering of homeologs while still permitting zippering of correct homologs. The degree of unzippering depends on the radius of the nucleus, such that complete unzippering of homeologous regions can only take place if the nucleus is large enough to pull the two chromosomes completely apart. A picture of meiotic pairing thus emerges that is fundamentally mechanical in nature, possibly explaining the purpose of active telomere forces, increased nuclear diameter, and the presence of 'Maverick' chromosomes in meiosis.

The Rabl configuration limits topological entanglement of chromosomes in budding yeast

The three dimensional organization of genomes remains mostly unknown due to their high degree of condensation. Biophysical studies predict that condensation promotes the topological entanglement of chromatin fibers and the inhibition of function. How organisms balance between functionally active genomes and a high degree of condensation remains to be determined. Here we hypothesize that the Rabl configuration, characterized by the attachment of centromeres and telomeres to the nuclear envelope, helps to reduce the topological entanglement of chromosomes. To test this hypothesis we developed a novel method to quantify chromosome entanglement complexity in 3D reconstructions obtained from Chromosome Conformation Capture (CCC) data. Applying this method to published data of the yeast genome, we show that computational models implementing the attachment of telomeres or centromeres alone are not sufficient to obtain the reduced entanglement complexity observed in 3D reconstructions. It is only when the centromeres and telomeres are attached to the nuclear envelope (i.e. the Rabl configuration) that the complexity of entanglement of the genome is comparable to that of the 3D reconstructions. We therefore suggest that the Rabl configuration is an essential player in the simplification of the entanglement of chromatin fibers.

Local loop opening in untangled ring polymer melts: a detailed "Feynman test" of models for the large scale structure

The conformational statistics of ring polymers in melts or dense solutions is strongly affected by their quenched microscopic topological state. The effect is particularly strong for untangled (i.e. non-concatenated and unknotted) rings, which are known to crumple and segregate. Here we study these systems using a computationally efficient multi-scale approach, where we combine massive simulations on the fiber level with the explicit construction of untangled ring melt configurations based on theoretical ideas for their large scale structure. We find (i) that topological constraints may be neglected on scales below the standard entanglement length, Le, (ii) that rings with a size 1 ≤ Lr/Le ≤ 30 exhibit nearly ideal lattice tree behavior characterized by primitive paths which are randomly branched on the entanglement scale, and (iii) that larger rings are compact with gyration radii Rg2(Lr) ∝ Lr2/3. The detailed comparison between equilibrated and constructed ensembles allows us to perform a "Feynman test" of our understanding of untangled rings: can we convert ideas for the large scale ring structure into algorithms for constructing (nearly) equilibrated ring melt samples? We show that most structural observables are quantitatively reproduced by two different construction schemes: hierarchical crumpling and ring melts derived from the analogy to interacting branched polymers. However, the latter fail the "Feynman test" with respect to the magnetic radius, Rm, which we have defined based on an analogy to magnetostatics. While Rm is expected to vanish for double-folded structures, the observed values of Rm2(Lr) ∝ Rg2(Lr) provide a simple and computationally convenient measure of the presence of a non-negligible amount of local loop opening in crumpled rings.

Quantitative Proteomics of the Mitotic Chromosome Scaffold Reveals the Association of BAZ1B with Chromosomal Axes

In mitosis, chromosomes achieve their characteristic shape through condensation, an essential process for proper segregation of the genome during cell division. A classical model for mitotic chromosome condensation proposes that non-histone proteins act as a structural framework called the chromosome scaffold. The components of the chromosome scaffold, such as DNA topoisomerase IIα (TOP2A) and structural maintenance of chromosomes protein 2 (SMC2), are necessary to generate stable mitotic chromosomes; however, the existence of this scaffold remains controversial. The aim of this study was to determine the protein composition of the chromosome scaffold. We used the DT40 chicken cell line to isolate mitotic chromosomes and extract the associated protein fraction, which could contain the chromosome scaffold. MS revealed a novel component of the chromosome scaffold, bromodomain adjacent to zinc finger 1B (BAZ1B), which was localized to the mitotic chromosome axis. Knocking out BAZ1B caused prophase delay because of altered chromosome condensation timing and mitosis progression errors, and the effect was aggravated if BAZ1A, a BAZ1B homolog, was simultaneously knocked out; however, protein composition of prometaphase chromosomes was normal. Our results suggest that BAZ1 proteins are essential for timely chromosome condensation at mitosis entry. Further characterization of the functional role of BAZ1 proteins would provide new insights into the timing of chromosome condensation.

Chain organization of human interphase chromosome determines the spatiotemporal dynamics of chromatin loci

We investigate spatiotemporal dynamics of human interphase chromosomes by employing a heteropolymer model that incorporates the information of human chromosomes inferred from Hi-C data. Despite considerable heterogeneities in the chromosome structures generated from our model, chromatins are organized into crumpled globules with space-filling (SF) statistics characterized by a single universal scaling exponent (ν = 1/3), and this exponent alone can offer a quantitative account of experimentally observed, many different features of chromosome dynamics. The local chromosome structures, whose scale corresponds to that of topologically associated domains (∼ 0.1 − 1 Mb), display dynamics with a fast relaxation time (≲ 1 − 10 sec); in contrast, the long-range spatial reorganization of the entire chromatin (≳O(102) Mb) occurs on a much slower time scale (≳ hour), providing the dynamic basis of cell-to-cell variability and glass-like behavior of chromosomes. Biological activities, modeled using stronger isotropic white noises added to active loci, accelerate the relaxation dynamics of chromatin domains associated with the low frequency modes and induce phase segregation between the active and inactive loci. Surprisingly, however, they do not significantly change the dynamics at local scales from those obtained under passive conditions. Our study underscores the role of chain organization of chromosome in determining the spatiotemporal dynamics of chromatin loci.

Chromosomes are giant chain molecules made of hundreds of megabase-long DNA intercalated with proteins. Structure and dynamics of interphase chromatin in space and time hold the key to understanding the cell type-dependent gene regulation. In this study, we establish that the crumpled and space-filling (SF) organization of chromatin fiber in the chromosome territory, characterized by a single scaling exponent, is sufficient to explain the complex spatiotemporal hierarchy in chromatin dynamics as well as the subdiffusive motion of the chromatin loci. While seemingly a daunting problem at a first glance, our study shows that relatively simple principles, rooted in polymer physics, can be used to grasp the essence of dynamical properties of the interphase chromatin.

Simulation of different three-dimensional polymer models of interphase chromosomes compared to experiments-an evaluation and review framework of the 3D genome organization

Despite all the efforts the three-dimensional higher-order architecture and dynamics in the cell nucleus are still debated. The regulation of genes, their transcription, replication, as well as differentiation in Eukarya is, however, closely connected to this architecture and dynamics. Here, an evaluation and review framework is setup to investigate the folding of a 30 nm chromatin fibre into chromosome territories by comparing computer simulations of two different chromatin topologies to experiments: The Multi-Loop-Subcompartment (MLS) model, in which small loops form rosettes connected by chromatin linkers, and the Random-Walk/Giant-Loop (RW/GL) model, in which large loops are attached to a flexible non-protein backbone, were simulated for various loop, rosette, and linker sizes. The 30 nm chromatin fibre was modelled as a polymer chain with stretching, bending, and excluded volume interactions. A spherical boundary potential simulated the confinement by other chromosomes and the nuclear envelope. Monte Carlo and Brownian Dynamics methods were applied to generate chain configurations at thermodynamic equilibrium. Both the MLS and the RW/GL models form chromosome territories, with different morphologies: The MLS rosettes form distinct subchromosomal domains, compatible in size as those from light microscopic observations. In contrast, the big RW/GL loops lead to a more homogeneous chromatin distribution. Only the MLS model agrees with the low overlap of chromosomes, their arms, and subchromosomal domains found experimentally. A review of experimental spatial distance measurements between genomic markers labelled by FISH as a function of their genomic separation from different publications and comparison to simulated spatial distances also favours an MLS-like model with loops and linkers of 63 to 126 kbp. The chromatin folding topology also reduces the apparent persistence length of the chromatin fibre to a value significantly lower than the free solution persistence length, explaining the low persistence lengths found various experiments. The predicted large spaces between the chromatin fibres allow typically sized biological molecules to reach nearly every location in the nucleus by moderately obstructed diffusion and disagrees with the much simplified assumption that defined channels between territories for molecular transport as in the Interchromosomal Domain (ICD) hypothesis exist. All this is also in agreement with recent selective high-resolution chromosome interaction capture (T2C) experiments, the scaling behaviour of the DNA sequence, the dynamics of the chromatin fibre, the nuclear diffusion of molecules, as well as other experiments. In summary, this polymer simulation framework compared to experimental data clearly favours only a quasi-chromatin fibre forming a stable multi-loop aggregate/rosette like genome organization and dynamics whose local topology is tightly connected to the global morphology and dynamics of the cell nucleus.

The elusive quest for RNA knots

Physical entanglement, and particularly knots arise spontaneously in equilibrated polymers that are sufficiently long and densely packed. Biopolymers are no exceptions: knots have long been known to occur in proteins as well as in encapsidated viral DNA. The rapidly growing number of RNA structures has recently made it possible to investigate the incidence of physical knots in this type of biomolecule, too. Strikingly, no knots have been found to date in the known RNA structures. In this Point of View Article we discuss the absence of knots in currently available RNAs and consider the reasons why knots in RNA have not yet been found, despite the expectation that they should exist in Nature. We conclude by singling out a number of RNA sequences that, based on the properties of their predicted secondary structures, are good candidates for knotted RNAs.

Hi-C-constrained physical models of human chromosomes recover functionally-related properties of genome organization

Combining genome-wide structural models with phenomenological data is at the forefront of efforts to understand the organizational principles regulating the human genome. Here, we use chromosome-chromosome contact data as knowledge-based constraints for large-scale three-dimensional models of the human diploid genome. The resulting models remain minimally entangled and acquire several functional features that are observed in vivo and that were never used as input for the model. We find, for instance, that gene-rich, active regions are drawn towards the nuclear center, while gene poor and lamina associated domains are pushed to the periphery. These and other properties persist upon adding local contact constraints, suggesting their compatibility with non-local constraints for the genome organization. The results show that suitable combinations of data analysis and physical modelling can expose the unexpectedly rich functionally-related properties implicit in chromosome-chromosome contact data. Specific directions are suggested for further developments based on combining experimental data analysis and genomic structural modelling.

Compaction and segregation of sister chromatids via active loop extrusion

The mechanism by which chromatids and chromosomes are segregated during mitosis and meiosis is a major puzzle of biology and biophysics. Using polymer simulations of chromosome dynamics, we show that a single mechanism of loop extrusion by condensins can robustly compact, segregate and disentangle chromosomes, arriving at individualized chromatids with morphology observed in vivo. Our model resolves the paradox of topological simplification concomitant with chromosome 'condensation', and explains how enzymes a few nanometers in size are able to control chromosome geometry and topology at micron length scales. We suggest that loop extrusion is a universal mechanism of genome folding that mediates functional interactions during interphase and compacts chromosomes during mitosis.

Modeling meiotic chromosome pairing: nuclear envelope attachment, telomere-led active random motion, and anomalous diffusion

The recognition and pairing of homologous chromosomes during meiosis is a complex physical and molecular process involving a combination of polymer dynamics and molecular recognition events. Two highly conserved features of meiotic chromosome behavior are the attachment of telomeres to the nuclear envelope and the active random motion of telomeres driven by their interaction with cytoskeletal motor proteins. Both of these features have been proposed to facilitate the process of homolog pairing, but exactly what role these features play in meiosis remains poorly understood. Here we investigate the roles of active motion and nuclear envelope tethering using a Brownian dynamics simulation in which meiotic chromosomes are represented by a Rouse polymer model subjected to tethering and active forces at the telomeres. We find that tethering telomeres to the nuclear envelope slows down pairing relative to the rates achieved by unattached chromosomes, but that randomly directed active forces applied to the telomeres speed up pairing dramatically in a manner that depends on the statistical properties of the telomere force fluctuations. The increased rate of initial pairing cannot be explained by stretching out of the chromosome conformation but instead seems to correlate with anomalous diffusion of sub-telomeric regions.

Confinement-Induced Glassy Dynamics in a Model for Chromosome Organization

Recent experiments showing scaling of the intrachromosomal contact probability, P(s)∼s(-1) with the genomic distance s, are interpreted to mean a self-similar fractal-like chromosome organization. However, scaling of P(s) varies across organisms, requiring an explanation. We illustrate dynamical arrest in a highly confined space as a discriminating marker for genome organization, by modeling chromosomes inside a nucleus as a homopolymer confined to a sphere of varying sizes. Brownian dynamics simulations show that the chain dynamics slows down as the polymer volume fraction (ϕ) inside the confinement approaches a critical value ϕ(c). The universal value of ϕ(c)(∞)≈0.44 for a sufficiently long polymer (N≫1) allows us to discuss genome dynamics using ϕ as the sole parameter. Our study shows that the onset of glassy dynamics is the reason for the segregated chromosome organization in humans (N≈3×10(9), ϕ≳ϕ(c)(∞)), whereas chromosomes of budding yeast (N≈10(8), ϕ<ϕ(c)(∞)) are equilibrated with no clear signature of such organization.

Current theoretical models fail to predict the topological complexity of the human genome

Understanding the folding of the human genome is a key challenge of modern structural biology. The emergence of chromatin conformation capture assays (e.g., Hi-C) has revolutionized chromosome biology and provided new insights into the three dimensional structure of the genome. The experimental data are highly complex and need to be analyzed with quantitative tools. It has been argued that the data obtained from Hi-C assays are consistent with a fractal organization of the genome. A key characteristic of the fractal globule is the lack of topological complexity (knotting or inter-linking). However, the absence of topological complexity contradicts results from polymer physics showing that the entanglement of long linear polymers in a confined volume increases rapidly with the length and with decreasing volume. In vivo and in vitro assays support this claim in some biological systems. We simulate knotted lattice polygons confined inside a sphere and demonstrate that their contact frequencies agree with the human Hi-C data. We conclude that the topological complexity of the human genome cannot be inferred from current Hi-C data.

Topology, structures, and energy landscapes of human chromosomes

Chromosome conformation capture experiments provide a rich set of data concerning the spatial organization of the genome. We use these data along with a maximum entropy approach to derive a least-biased effective energy landscape for the chromosome. Simulations of the ensemble of chromosome conformations based on the resulting information theoretic landscape not only accurately reproduce experimental contact probabilities, but also provide a picture of chromosome dynamics and topology. The topology of the simulated chromosomes is probed by computing the distribution of their knot invariants. The simulated chromosome structures are largely free of knots. Topologically associating domains are shown to be crucial for establishing these knotless structures. The simulated chromosome conformations exhibit a tendency to form fibril-like structures like those observed via light microscopy. The topologically associating domains of the interphase chromosome exhibit multistability with varying liquid crystalline ordering that may allow discrete unfolding events and the landscape is locally funneled toward "ideal" chromosome structures that represent hierarchical fibrils of fibrils.

Ring polymers in the melt state: the physics of crumpling

The conformational statistics of ring polymers in melts or dense solutions is strongly affected by their quenched microscopic topological state. The effect is particularly strong for nonconcatenated unknotted rings, which are known to crumple and segregate and which have been implicated as models for the generic behavior of interphase chromosomes. Here we use a computationally efficient multiscale approach to show that melts of rings of total contour length Lr can be quantitatively mapped onto melts of interacting lattice trees with gyration radii ⟨R(g)(2)(Lr)⟩ ∝ L(r)(2ν) and ν = 0.32 ± 0.01.

From a melt of rings to chromosome territories: the role of topological constraints in genome folding

We review pro and contra of the hypothesis that generic polymer properties of topological constraints are behind many aspects of chromatin folding in eukaryotic cells. For that purpose, we review, first, recent theoretical and computational findings in polymer physics related to concentrated, topologically simple (unknotted and unlinked) chains or a system of chains. Second, we review recent experimental discoveries related to genome folding. Understanding in these fields is far from complete, but we show how looking at them in parallel sheds new light on both.

Organization of the mitotic chromosome

Mitotic chromosomes are among the most recognizable structures in the cell, yet for over a century their internal organization remains largely unsolved. We applied chromosome conformation capture methods, 5C and Hi-C, across the cell cycle and revealed two distinct three-dimensional folding states of the human genome. We show that the highly compartmentalized and cell type-specific organization described previously for nonsynchronous cells is restricted to interphase. In metaphase, we identified a homogenous folding state that is locus-independent, common to all chromosomes, and consistent among cell types, suggesting a general principle of metaphase chromosome organization. Using polymer simulations, we found that metaphase Hi-C data are inconsistent with classic hierarchical models and are instead best described by a linearly organized longitudinally compressed array of consecutive chromatin loops.

Exploring the three-dimensional organization of genomes: interpreting chromatin interaction data

How DNA is organized in three dimensions inside the cell nucleus and how this affects the ways in which cells access, read and interpret genetic information are among the longest standing questions in cell biology. Using newly developed molecular, genomic and computational approaches based on the chromosome conformation capture technology (such as 3C, 4C, 5C and Hi-C), the spatial organization of genomes is being explored at unprecedented resolution. Interpreting the increasingly large chromatin interaction data sets is now posing novel challenges. Here we describe several types of statistical and computational approaches that have recently been developed to analyse chromatin interaction data.

Structure of metaphase chromosomes: a role for effects of macromolecular crowding

In metaphase chromosomes, chromatin is compacted to a concentration of several hundred mg/ml by mechanisms which remain elusive. Effects mediated by the ionic environment are considered most frequently because mono- and di-valent cations cause polynucleosome chains to form compact ~30-nm diameter fibres in vitro, but this conformation is not detected in chromosomes in situ. A further unconsidered factor is predicted to influence the compaction of chromosomes, namely the forces which arise from crowding by macromolecules in the surrounding cytoplasm whose measured concentration is 100-200 mg/ml. To mimic these conditions, chromosomes were released from mitotic CHO cells in solutions containing an inert volume-occupying macromolecule (8 kDa polyethylene glycol, 10.5 kDa dextran, or 70 kDa Ficoll) in 100 µM K-Hepes buffer, with contaminating cations at only low micromolar concentrations. Optical and electron microscopy showed that these chromosomes conserved their characteristic structure and compaction, and their volume varied inversely with the concentration of a crowding macromolecule. They showed a canonical nucleosomal structure and contained the characteristic proteins topoisomerase IIα and the condensin subunit SMC2. These observations, together with evidence that the cytoplasm is crowded in vivo, suggest that macromolecular crowding effects should be considered a significant and perhaps major factor in compacting chromosomes. This model may explain why ~30-nm fibres characteristic of cation-mediated compaction are not seen in chromosomes in situ. Considering that crowding by cytoplasmic macromolecules maintains the compaction of bacterial chromosomes and has been proposed to form the liquid crystalline chromosomes of dinoflagellates, a crowded environment may be an essential characteristic of all genomes.

The ancestral role of ATP hydrolysis in type II topoisomerases: prevention of DNA double-strand breaks

Type II DNA topoisomerases (topos) catalyse changes in DNA topology by passing one double-stranded DNA segment through another. This reaction is essential to processes such as replication and transcription, but carries with it the inherent danger of permanent double-strand break (DSB) formation. All type II topos hydrolyse ATP during their reactions; however, only DNA gyrase is able to harness the free energy of hydrolysis to drive DNA supercoiling, an energetically unfavourable process. A long-standing puzzle has been to understand why the majority of type II enzymes consume ATP to support reactions that do not require a net energy input. While certain type II topos are known to 'simplify' distributions of DNA topoisomers below thermodynamic equilibrium levels, the energy required for this process is very low, suggesting that this behaviour is not the principal reason for ATP hydrolysis. Instead, we propose that the energy of ATP hydrolysis is needed to control the separation of protein-protein interfaces and prevent the accidental formation of potentially mutagenic or cytotoxic DSBs. This interpretation has parallels with the actions of a variety of molecular machines that catalyse the conformational rearrangement of biological macromolecules.

Concentrated DNA Rheology and Microrheology

We present mechanical measurements of the frequency-dependent linear viscoelastic storage and loss moduli, G′(ω) and G″(ω), and the yield stress, τy, and yield strain, γy, for calf thymus DNA (13 kbp) over a range of mitotically relevant concentrations from CDNA = 1 to 10 mg/ml. For large CDNA, we find a dominant plateau elasticity, G′p, at high ω. As ω decreases, G′ falls until it is equal to G′ at the crossover frequency, ωc, below which G″ dominates. We measure G′p ∼ CDNA2.25 and ωc ∼ CDNA−2.4, consistent with scaling exponents for classical polymer solutions. The mechanical |G*(ω)| agree well with those measured using a new microrheological technique based on video tracking microscopy of thermally-driven fluorescent colloidal spheres and a frequency-dependent Stokes-Einstein equation. We have developed this technique to probe how enzymes, typically available in small quantities, can affect the rheology of the DNA. Using it, we report preliminary measurements of a higher ωc for a DNA network in which the ATP-powered enzyme Topoisomerase II transiently cuts and rebinds the DNA, thereby relaxing entanglements.

Repulsive forces between looping chromosomes induce entropy-driven segregation

One striking feature of chromatin organization is that chromosomes are compartmentalized into distinct territories during interphase, the degree of intermingling being much smaller than expected for linear chains. A growing body of evidence indicates that the formation of loops plays a dominant role in transcriptional regulation as well as the entropic organization of interphase chromosomes. Using a recently proposed model, we quantitatively determine the entropic forces between chromosomes. This Dynamic Loop Model assumes that loops form solely on the basis of diffusional motion without invoking other long-range interactions. We find that introducing loops into the structure of chromatin results in a multi-fold higher repulsion between chromosomes compared to linear chains. Strong effects are observed for the tendency of a non-random alignment; the overlap volume between chromosomes decays fast with increasing loop number. Our results suggest that the formation of chromatin loops imposes both compartmentalization as well as order on the system without requiring additional energy-consuming processes.

Control of the flow properties of DNA by topoisomerase II and its targeting inhibitor

The flow properties of DNA are important for understanding cell division and, indirectly, cancer therapy. DNA topology controlling enzymes such as topoisomerase II are thought to play an essential role. We report experiments showing how double-strand passage facilitated by topoisomerase II controls DNA rheology. For this purpose, we have measured the elastic storage and viscous loss moduli of a model system comprising bacteriophage λ-DNA and human topoisomerase IIα using video tracking of the Brownian motion of colloidal probe particles. We found that the rheology is critically dependent on the formation of temporal entanglements among the DNA molecules with a relaxation time of ∼1 s. We observed that topoisomerase II effectively removes these entanglements and transforms the solution from an elastic physical gel to a viscous fluid depending on the consumption of ATP. A second aspect of this study is the effect of the generic topoisomerase II inhibitor adenylyl-imidodiphosphate (AMP-PNP). In mixtures of AMP-PNP and ATP, the double-strand passage reaction gets blocked and progressively fewer entanglements are relaxed. A total replacement of ATP by AMP-PNP results in a temporal increase in elasticity at higher frequencies, but no transition to an elastic gel with fixed cross-links.

Diffusion-driven looping provides a consistent framework for chromatin organization

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

Biopolymer organization upon confinement

Biopolymers in vivo are typically subject to spatial restraints, either as a result of molecular crowding in the cellular medium or of direct spatial confinement. DNA in living organisms provides a prototypical example of a confined biopolymer. Confinement prompts a number of biophysics questions. For instance, how can the high level of packing be compatible with the necessity to access and process the genomic material? What mechanisms can be adopted in vivo to avoid the excessive geometrical and topological entanglement of dense phases of biopolymers? These and other fundamental questions have been addressed in recent years by both experimental and theoretical means. A review of the results, particularly of those obtained by numerical studies, is presented here. The review is mostly devoted to DNA packaging inside bacteriophages, which is the best studied example both experimentally and theoretically. Recent selected biophysical studies of the bacterial genome organization and of chromosome segregation in eukaryotes are also covered.

Topological origins of chromosomal territories

Using freely jointed polymer model we compare equilibrium properties of crowded polymer chains whose segments are either permeable or not permeable for other segments to pass through. In particular, we addressed the question whether non-permeability of long chain molecules, in the absence of excluded volume effect, is sufficient to compartmentalize highly crowded polymer chains, similarly to what happens during formation of chromosomal territories in interphase nuclei. Our results indicate that even polymers without excluded volume compartmentalize and show strongly reduced intermingling when they are mutually non-permeable. Judging from the known fact that chromatin fibres originating from different chromosomes show very limited intermingling in interphase nuclei, we propose that regular chromatin fibres during chromosome decondensation can hardly serve as a substrate of cellular type II DNA topoisomerases.

Entropic organization of interphase chromosomes

Chromosomes are not distributed randomly in nuclei. Appropriate positioning can activate (or repress) genes by bringing them closer to active (or inactive) compartments like euchromatin (or heterochromatin), and this is usually assumed to be driven by specific local forces (e.g., involving H bonds between nucleosomes or between nucleosomes and the lamina). Using Monte Carlo simulations, we demonstrate that nonspecific (entropic) forces acting alone are sufficient to position and shape self-avoiding polymers within a confining sphere in the ways seen in nuclei. We suggest that they can drive long flexible polymers (representing gene-rich chromosomes) to the interior, compact/thick ones (and heterochromatin) to the periphery, looped (but not linear) ones into appropriately shaped (ellipsoidal) territories, and polymers with large terminal beads (representing centromeric heterochromatin) into peripheral chromocenters. Flexible polymers tend to intermingle less than others, which is in accord with observations that gene-dense (and so flexible) chromosomes make poor translocation partners. Thus, entropic forces probably participate in the self-organization of chromosomes within nuclei.

The why and how of DNA unlinking

The nucleotide sequence of DNA is the repository of hereditary information. Yet, it is now clear that the DNA itself plays an active role in regulating the ability of the cell to extract its information. Basic biological processes, including control of gene transcription, faithful DNA replication and segregation, maintenance of the genome and cellular differentiation are subject to the conformational and topological properties of DNA in addition to the regulation imparted by the sequence itself. How do these DNA features manifest such striking effects and how does the cell regulate them? In this review, we describe how misregulation of DNA topology can lead to cellular dysfunction. We then address how cells prevent these topological problems. We close with a discussion on recent theoretical advances indicating that the topological problems, themselves, can provide the cues necessary for their resolution by type-2 topoisomerases.

Structure and dynamics of interphase chromosomes

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, Drosophila, 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.

Hin-mediated DNA knotting and recombining promote replicon dysfunction and mutation

Background

The genetic code imposes a dilemma for cells. The DNA must be long enough to encode for the complexity of an organism, yet thin and flexible enough to fit within the cell. The combination of these properties greatly favors DNA collisions, which can knot and drive recombination of the DNA. Despite the well-accepted propensity of cellular DNA to collide and react with itself, it has not been established what the physiological consequences are.

Results

Here we analyze the effects of recombined and knotted plasmids in E. coli using the Hin site-specific recombination system. We show that Hin-mediated DNA knotting and recombination (i) promote replicon loss by blocking DNA replication; (ii) block gene transcription; and (iii) cause genetic rearrangements at a rate three to four orders of magnitude higher than the rate for an unknotted, unrecombined plasmid.

Conclusion

These results show that DNA reactivity leading to recombined and knotted DNA is potentially toxic and may help drive genetic evolution.

Topological information embodied in local juxtaposition geometry provides a statistical mechanical basis for unknotting by type-2 DNA topoisomerases

Topoisomerases may unknot by recognizing specific DNA juxtapositions. The physical basis of this hypothesis is investigated by considering single-loop conformations in a coarse-grained polymer model. We determine the statistical relationship between the local geometry of a juxtaposition of two chain segments and whether the loop is knotted globally, and ascertain how the knot/unknot topology is altered by a topoisomerase-like segment passage at the juxtaposition. Segment passages at a "free" juxtaposition tend to increase knot probability. In contrast, segment passages at a "hooked" juxtaposition cause more transitions from knot to unknot than vice versa, resulting in a steady-state knot probability far lower than that at topological equilibrium. The reduction in knot population by passing chain segments through a hooked juxtaposition is more prominent for loops of smaller sizes, n, but remains significant even for larger loops: steady-state knot probability is only approximately 2%, and approximately 5% of equilibrium, respectively, for n=100 and 500 in the model. An exhaustive analysis of approximately 6000 different juxtaposition geometries indicates that the ability of a segment passage to unknot correlates strongly with the juxtaposition's "hookedness". Remarkably, and consistent with experiments on type-2 topoisomerases from different organisms, the unknotting potential of a juxtaposition geometry in our polymer model correlates almost perfectly with its corresponding decatenation potential. These quantitative findings suggest that it is possible for topoisomerases to disentangle by acting selectively on juxtapositions with "hooked" geometries.

A polymer model for the structural organization of chromatin loops and minibands in interphase chromosomes

A quantitative model of interphase chromosome higher-order structure is presented based on the isochore model of the genome and results obtained in the field of copolymer research. G1 chromosomes are approximated in the model as multiblock copolymers of the 30-nm chromatin fiber, which alternately contain two types of 0.5- to 1-Mbp blocks (R and G minibands) differing in GC content and DNA-bound proteins. A G1 chromosome forms a single-chain string of loop clusters (micelles), with each loop approximately 1-2 Mbp in size. The number of approximately 20 loops per micelle was estimated from the dependence of geometrical versus genomic distances between two points on a G1 chromosome. The greater degree of chromatin extension in R versus G minibands and a difference in the replication time for these minibands (early S phase for R versus late S phase for G) are explained in this model as a result of the location of R minibands at micelle cores and G minibands at loop apices. The estimated number of micelles per nucleus is close to the observed number of replication clusters at the onset of S phase. A relationship between chromosomal and nuclear sizes for several types of higher eukaryotic cells (insects, plants, and mammals) is well described through the micelle structure of interphase chromosomes. For yeast cells, this relationship is described by a linear coil configuration of chromosomes.

Topological complexity of different populations of pBR322 as visualized by two-dimensional agarose gel electrophoresis

Neutral/neutral two-dimensional (2D) agarose gelelectrophoresis was used to investigate populations of the different topological conformations that pBR322 can adopt in vivo in bacterial cells as well as in Xenopus egg extracts. To help in interpretation and identification of all the different signals, undigested as well as DNA samples pretreated with DNase I, topoisomerase I and topoisomerase II were analyzed. The second dimension of the 2D gel system was run with or without ethidium bromide to account for any possible changes in the migration behavior of DNA molecules caused by intercalation of this planar agent. Finally, DNA samples were isolated from a recA-strain of Escherichia coli , as well as after direct labeling of the replication intermediates in extracts of Xenopus laevis eggs. Altogether, the results obtained demonstrated that 2D gels can be readily used to identify most of the complex topological populations that circular molecules can adopt in vivo in both bacteria and eukaryotic cells.

Polymer models of meiotic and mitotic chromosomes

Polymers tied together by constraints exhibit an internal pressure; this idea is used to analyze physical properties of the bottle-brush-like chromosomes of meiotic prophase that consist of polymer-like flexible chromatin loops, attached to a central axis. Using a minimal number of experimental parameters, semiquantitative predictions are made for the bending rigidity, radius, and axial tension of such brushes, and the repulsion acting between brushes whose bristles are forced to overlap. The retraction of lampbrush loops when the nascent transcripts are stripped away, the oval shape of diplotene bivalents between chiasmata, and the rigidity of pachytene chromosomes are all manifestations of chromatin pressure. This two-phase (chromatin plus buffer) picture that suffices for meiotic chromosomes has to be supplemented by a third constituent, a chromatin glue to understand mitotic chromosomes, and explain how condensation can drive the resolution of entanglements. This process resembles a thermal annealing in that a parameter (the affinity of the glue for chromatin and/or the affinity of the chromatin for buffer) has to be tuned to achieve optimal results. Mechanical measurements to characterize this protein-chromatin matrix are proposed. Finally, the propensity for even slightly chemically dissimilar polymers to phase separate (cluster like with like) can explain the apparent segregation of the chromatin into A + T- and G + C-rich regions revealed by chromosome banding.

Electron tomography of metaphase nucleolar organizer regions: evidence for a twisted-loop organization

Metaphase nucleolar organizer regions (NORs), one of four types of chromosome bands, are located on human acrocentric chromosomes. They contain r-chromatin, i.e., ribosomal genes complexed with proteins such as upstream binding factor and RNA polymerase I, which are argyrophilic NOR proteins. Immunocytochemical and cytochemical labelings of these proteins were used to reveal r-chromatin in situ and to investigate its spatial organization within NORs by confocal microscopy and by electron tomography. For each labeling, confocal microscopy revealed small and large double-spotted NORs and crescent-shaped NORs. Their internal three-dimensional (3D) organization was studied by using electron tomography on specifically silver-stained NORs. The 3D reconstructions allow us to conclude that the argyrophilic NOR proteins are grouped as a fiber of 60-80 nm in diameter that constitutes either one part of a turn or two or three turns of a helix within small and large double-spotted NORs, respectively. Within crescent-shaped NORs, virtual slices reveal that the fiber constitutes several longitudinally twisted loops, grouped as two helical 250- to 300-nm coils, each centered on a nonargyrophilic axis of condensed chromatin. We propose a model of the 3D organization of r-chromatin within elongated NORs, in which loops are twisted and bent to constitute one basic chromatid coil.