Search

All life forms are miraculous, but some are more inexplicable than others. Trypanosomes are by far one of the most puzzling organisms on Earth: their mitochondrial genome, also called kinetoplast DNA (kDNA) forms an Olympic-ring-like network of interlinked DNA circles, challenging conventional paradigms in both biology and physics. In this review, I will discuss kDNA from the astonished perspective of a polymer physicist and tell a story of how a single sub-celluar structure from a blood-dwelling parasite is inspiring generations of polymer chemists and physicists to create new catenated materials., Comment: review paper on kDNA

In this Letter we study the phase transition between amorphous ices and the nature of the hysteresis cycle separating them. We discover that a topological transition takes place as the system transforms from low-density amorphous ice (LDA) at low pressures, to high-density amorphous ice (HDA) at high pressures. Specifically, we uncover that the hydrogen bond network (HBN) displays qualitatively different topologies in the LDA and HDA phases: the former characterised by disentangled loop motifs, while the latter displaying topologically complex and metastable Hopf-linked and knotted configurations. At the phase transition, the transient opening of the HBN topological motifs yields mechanical fragility on the macroscale. Our results provide a detailed microscopic description of the topological nature of the phase transition and the hysteresis cycle between amorphous ices. We argue that the topological transition discovered in this work may not only improve our understanding of amorphous ices but represent a generic mechanism for the densification of network-forming materials.

Loop extrusion is one of the main processes shaping chromosome organisation across the cell cycle, yet its role in regulating DNA entanglement and nucleoplasm viscoelasticity remains overlooked. We simulate entangled solutions of linear polymers under the action of generic Loop Extruding Factors (LEF) with a model that fully accounts for topological constraints and LEF-DNA uncrossability. We discover that extrusion drives the formation of bottle-brush-like structures which significantly lower the entanglement and effective viscosity of the system through an active fluidification mechanism. Interestingly, this fluidification displays an optimum at one LEF every 300-3000 basepairs. In marked contrast with entangled linear chains, the viscosity of extruded chains scales linearly with polymer length, yielding up to 1000-fold fluidification. Our results illuminate how loop extrusion contributes to actively modulate genome entanglement and viscoelasticity in vivo.

The kinetoplast DNA (kDNA) is the archetype of a two-dimensional Olympic network, composed of thousands of DNA minicircles and found in the mitochondrion of certain parasites. The evolution, replication and self-assembly of this structure are fascinating open questions in biology that can also inform us how to realise synthetic Olympic networks in vitro. To obtain a deeper understanding of the structure and assembly of kDNA networks, we sequenced the Crithidia fasciculata kDNA genome and performed high-resolution Atomic Force Microscopy (AFM) and analysis of kDNA networks that had been partially digested by selected restriction enzymes. We discovered that these topological perturbations lead to networks with significantly different geometrical features and morphologies with respect to the unperturbed kDNA, and that these changes are strongly dependent on the class of DNA circles targeted by the restriction enzymes. Specifically, cleaving maxicircles leads to a dramatic reduction in network size once adsorbed onto the surface, whilst cleaving both maxicircles and a minor class of minicircles yields non-circular and deformed structures. We argue that our results are a consequence of a precise positioning of the maxicircles at the boundary of the network, and we discuss our findings in the context of kDNA biogenesis, design of artificial Olympic networks and detection of in vivo perturbations.

DNA loop extrusion is arguably one of the most important players in genome organization. The precise mechanism by which loop extruding factors (LEFs) work is still unresolved and much debated. One of the major open questions in this field is how do LEFs establish and maintain unidirectional motion along DNA. In this paper, we use High-Speed AFM data to show that condensin hinge domain displays a structural, geometric constraint on the angle within which it can extend with respect to the DNA-bound domains. Using computer simulations, we then show that such a geometrical constraint results in a local symmetry breaking and is enough to rectify the extrusion process, yielding unidirectional loop extrusion along DNA. Our work highlights an overlooked geometric aspect of the loop extrusion process that may have a universal impact on SMC function across organisms.

Understanding the connection between a material's microscopic network topology and its macroscopic properties is a fundamental open challenge across a wide range of physical systems. Here, we investigate the elasticity of DNA nanostar hydrogels - a model system for networks with limited valence - by coupling rheology measurements, confocal imaging and molecular dynamics simulations. We discover that these networks display a large degree of interpenetration and that loops within the network are topologically linked, forming a percolating network-within-network structure. Strikingly, we discover that the onset of topological links between shortest loops fully determines the elasticity of these physical gels. We thus argue that the elasticity of limited valence gels is not dictated by their fractal microstructure, but by the topology of looped motifs. Our findings highlight the emergence of "topological elasticity" as a previously overlooked mechanism in generic network-forming liquids and physical gels and will inspire the design of topologically-controllable material behaviours.

In this paper we investigate the effects of varying cation valency and concentration on the rheology of entangled lambda DNA solutions. We show that monovalent cations moderately increase the viscoelasticty of the solutions mainly by stabilising linear condensation of lambda DNA ``monomers'' via hybridisation of their sticky ends. On the contrary, divalent cations have a far more complex and dramatic effect on the rheology of the solution and we observe evidence of inter-molecular DNA-DNA bridging by Mg2+. We argue that these results may be interesting in the context of dense solutions of single and double stranded DNA, e.g. in vivo or in biotechnology applications such as DNA origami and DNA hydrogels., Comment: 8 pages, 5 figures

The insertion of HIV and other DNA elements within genomes underpins both genetic diversity and disease when unregulated. Most of these insertions are not random and occupy specific positions within the genome but the physical mechanisms underlying the integration site selection are poorly understood. Here we perform Molecular Dynamics simulations to study the insertion of DNA elements, such as HIV viral DNA or transposons, into naked DNA or chromatin substrate. More specifically, we explore the role of loops in the DNA substrate and discover that they act as "geometric catalysts" for DNA integration. Additionally, we discover that the 1D and 3D clustering of loops affects the distribution of integration sites. Finally, we show that loops may compete with nucleosomes at attracting DNA integrations. These results may be tested in vitro and they may help to understand patterns of DNA insertions with implications in genome evolution and gene therapy.

Knots are deeply entangled with every branch of science. One of the biggest open challenges in knot theory is to formalise a knot invariant that can unambiguously and efficiently distinguish any two knotted curves. Additionally, the conjecture that the geometrical embedding of a curve encodes information on its underlying topology is, albeit physically intuitive, far from proven. Here we attempt to tackle both these outstanding challenges by proposing a neural network (NN) approach that takes as input a geometric representation of a knotted curve and tries to make predictions of the curve's topology. Intriguingly, we discover that NNs trained with a so-called geometrical "local writhe" representation of a knot can distinguish curves that share one or many topological invariants and knot polynomials, such as mutant and composite knots, and can thus classify knotted curves more precisely than some knot polynomials. Additionally, we also show that our approach can be scaled up to classify all prime knots up to 10-crossings with more than 95\% accuracy. Finally, we show that our NNs can also be trained to solve knot localisation problems on open and closed curves. Our main discovery is that the pattern of "local writhe" is a potentially unique geometric signature of the underlying topology of a curve. We hope that our results will suggest new methods for quantifying generic entanglements in soft matter and even inform new topological invariants., Comment: Accepted in Soft Matter

The process of polymer condensation, i.e. the formation of bonds between reactive end-groups, is ubiquitous in both industry and biology. Here we study generic systems undergoing polymer condensation in competition with cyclisation. Using a generalised Smoluchowski theory, molecular dynamics simulations and experiments using DNA and T4 ligase, we find that this system displays a transition, from a regime with finite-length chains at infinite time and dominated by rings to one dominated by linear polymers that grow in time. Finally, we show that fluids prepared close to the transition may have profoundly different compositions and rheology at large condensation times., Comment: 6 pages, 7 figures (main)

The Kinetoplast DNA (kDNA) is a two-dimensional Olympic-ring-like network of mutually linked 2.5 kb-long DNA minicircles found in certain parasites called Trypanosomes. Understanding the self-assembly and replication of this structure are not only major open questions in biology but can also inform the design of synthetic topological materials. Here we report the first high-resolution, single-molecule study of kDNA network topology using AFM and steered molecular dynamics simulations. We map out the DNA density within the network and the distribution of linking number and valence of the minicircles. We also characterise the DNA hubs that surround the network and show that they cause a buckling transition akin to that of a 2D elastic thermal sheet in the bulk. Intriguingly, we observe a broad distribution of density and valence of the minicircles, indicating heterogeneous network structure and individualism of different kDNA structures. Our findings explain outstanding questions in the field and offer single-molecule insights into the properties of a unique topological material.

Inspired by how certain proteins "sense" knots and entanglements in DNA molecules, here we ask if there exist local geometric features that may be used as a read-out of the underlying topology of generic polymers. We perform molecular simulations of knotted and linked semiflexbile polymers and study four geometric measures to predict topological entanglements: local curvature, local density, local 1D writhe and non-local 3D writhe. We discover that local curvature is a poor predictor of entanglements. In contrast, segments with maximum local density or writhe correlate as much as 90% of the time with the shortest knotted and linked arcs. We find that this accuracy is preserved across different knot types and also under significant spherical confinement, which is known to delocalise essential crossings in knotted polymers. We further discover that non-local 3D writhe is the best geometric read-out of knot location. Finally, we discuss how these geometric features may be used to computationally analyse entanglements in generic polymer melts and gels., Comment: accepted in ACS Polymers Au

In spite of the nanoscale and single-molecule insights into how nucleoid associated proteins (NAPs) interact with DNA, their role in modulating the mesoscale viscoelasticity of the entangled genome in vivo has been overlooked so far. By combining microrheology and molecular dynamics simulation we find that the important NAP called Integration Host Factor (IHF) lowers the viscosity of entangled $\lambda$DNA 20-fold at physiological concentrations and stoichiometries. We argue that IHF may act as a "genomic fluidiser", reducing the effective viscosity of the nucleoid $\sim$200-fold. Our results suggest a previously unappreciated key role of IHF in regulating DNA dynamics and re-organisation in vivo

Phase separation is as familiar as watching vinegar separating from oil in vinegrette. The observation that phase separation of proteins and nucleic acids is widespread in living cells has opened an entire field of research into the biological significance and the biophysical mechanisms of phase separation and protein condensation in biology. Recent evidence indicate that certain proteins and nucleic acids condensates are not simple liquids and instead display both viscous and elastic behaviours, which in turn may have biological significance. The aim of this perspective is to review the state-of-the-art of this quickly emerging field focusing on the material and rheological properties of protein condensates. Finally, we discuss the different techniques that can be employed to quantify the viscoelasticity of condensates and highlight potential future directions and opportunities for interdisciplinary cross-talk between chemists, physicists and biologists., Comment: review paper; accepted in JACS Au

Loop extrusion convincingly describes how certain Structural Maintenance of Chromosome (SMC) proteins mediate the formation of large DNA loops. Yet, most of the existing computational models cannot reconcile recent in vitro observations showing that condensins can traverse each other, bypass large roadblocks and perform steps longer than its own size. To fill this gap, we propose a three-dimensional (3D) "trans-grabbing" model for loop extrusion which not only reproduces the experimental features of loop extrusion by one SMC complex, but also predicts the formation of so-called "Z-loops" via the interaction of two or more SMCs extruding along the same DNA substrate. By performing Molecular Dynamics simulations of this model we discover that the experimentally observed asymmetry in the different types of Z-loops is a natural consequence of the DNA tethering in vitro. Intriguingly, our model predicts this bias to disappear in absence of tethering and a third type of Z-loop, which has not yet been identified in experiments, to appear. Our model naturally explains road-block bypassing and the appearance of steps larger than the SMC size as a consequence of non-contiguous DNA grabbing. Finally, it is the first to our knowledge to address how Z-loops and bypassing might occur in a way that is broadly consistent with existing cis-only 1D loop extrusion models., Comment: to appear in biophysical journal

How type 2 Topoisomerase (TopoII) proteins relax and simplify the topology of DNA molecules is one of the most intriguing open questions in biophysics. Most of the existing models neglect the dynamics of TopoII which is characteristics for proteins searching their targets via facilitated diffusion. Here, we show that dynamic binding of TopoII speeds up the topological relaxation of knotted substrates by enhancing the search of the knotted arc. Intriguingly, this in turn implies that the timescale of topological relaxation is virtually independent of the substrate length. We then discover that considering binding biases due to facilitated diffusion on looped substrates steers the sampling of the topological space closer to the boundaries between different topoisomers yielding an optimally fast topological relaxation. We discuss our findings in the context of topological simplification in vitro and in vivo.

Understanding and controlling the rheology of polymeric complex fluids that are pushed out-of-equilibrium is a fundamental problem in both industry and biology. For example, to package, repair, and replicate DNA, cells use enzymes to constantly manipulate DNA topology, length, and structure. Inspired by this, here we engineer and study DNA-based complex fluids that undergo enzymatically-driven topological and architectural alterations via restriction endonuclease (RE) reactions. We show that these systems display time-dependent rheological properties that depend on the concentrations and properties of the comprising DNA and REs. Through time-resolved microrheology experiments and Brownian Dynamics simulations, we show that conversion of supercoiled to linear DNA topology leads to a monotonic increase in viscosity. On the other hand, the viscosity of entangled linear DNA undergoing fragmentation displays a universal decrease that we rationalize using living polymer theory. Finally, to showcase the tunability of these behaviours, we design a DNA fluid that exhibits a time-dependent increase, followed by a temporally-gated decrease, of its viscosity. Our results present a class of polymeric fluids that leverage naturally occurring enzymes to drive diverse time-varying rheology by performing architectural alterations to the constituents.

Polymer composites are ubiquitous in biology and industry alike, owing to their emergent desirable mechanical properties not attainable in single-species systems. At the same time, polymer topology has been shown to play a key role in tuning the rheology of polymeric fluids. However, how topology impacts the rheology of composites remains poorly understood. Here, we create composites of rigid rods (microtubules) polymerized within entangled solutions of flexible linear and ring polymers (DNA). We couple linear and nonlinear optical tweezers microrheology with confocal microscopy and scaled particle theory to show that composites of linear DNA and microtubules exhibit a strongly non-monotonic dependence of elasticity and stiffness on microtubule concentration due to depletion-driven polymerization and flocculation of microtubules. In contrast, composites of ring DNA and microtubules show a much more modest monotonic increase in elastic strength with microtubule concentration, which we demonstrate arises from the increased ability of rings to mix with microtubules., Comment: 23 pages

Understanding how entanglements affect the behaviour of polymeric complex fluids is an open challenge in many fields. To elucidate the nature and consequence of entanglements in dense polymer solutions, we propose a novel method: a "dynamical entanglement analysis" (DEA) to extract spatio-temporal entanglement structures from the pair-wise displacement correlation of entangled chains. By applying this method to large-scale Molecular Dynamics simulations of linear and unknotted, nonconcatenated ring polymers, we find a strong and unexpected cooperative dynamics: the footprint of mutual entrainment between entangled chains. We show that DEA is a powerful and sensitive probe that reveals previously unnoticed, and architecture-dependent, spatio-temporal structures of dynamical entanglement in polymeric solutions. We also propose a mean-field approximation of our analysis which provides previously under-appreciated physical insights into the dynamics of generic entangled polymers. We envisage DEA will be useful to analyse the dynamical evolution of entanglements in generic polymeric systems such as blends and composites., Comment: To appear in ACS MacroLetters

Topological entanglements are abundant, and often detrimental, in polymeric systems in biology and materials science. Here we theoretically investigate the topological simplification of knots by diffusing slip-links (SLs), which may represent biological or synthetic molecules, such as proteins on the genome or cyclodextrines in slide-ring gels. We find that SLs entropically compete with knots and can localise them, greatly facilitating their downstream simplification by transient strand-crossing. We further show that the efficiency of knot localisation strongly depends on the topology of the SL network and, informed by our findings, discuss potential strategies to control the topology of biological and synthetic materials.

Lamina-associated domains (LADs) cover a large part of the human genome and are thought to play a major role in shaping the nuclear architectural landscape. Here, we perform polymer simulations, microscopy and mass spectrometry to dissect the roles played by heterochromatin- and lamina-mediated interactions in nuclear organisation. Our model explains the conventional organisation of heterochromatin and euchromatin in growing cells and the pathological organisation found in oncogene-induced senescence and progeria. We show that the experimentally observed changes in the locality of contacts in senescent and progeroid cells can be explained as arising due to phase transitions in the system. Within our simulations LADs are highly stochastic, as in experiments. Our model suggests that, once established, the senescent phenotype should be metastable even if lamina-mediated interactions were reinstated. Overall, our simulations uncover a generic physical mechanism that can regulate heterochromatin segregation and LAD formation in a wide range of mammalian nuclei., Comment: 20 pages, 8 figures

DNA is increasingly employed for bio and nanotechnology thanks to its exquisite versatility and designability. Most of its use is limited to linearised and torsionally relaxed DNA but non-trivial architectures and torsionally constrained -- or supercoiled -- DNA plasmids are largely neglected; this is partly due to the limited understanding of how supercoiling affects the rheology of entangled DNA. To address this open question we perform large scale Molecular Dynamics (MD) simulations of entangled solutions of DNA plasmids modelled as twistable chains. We discover that, contrarily to what generally assumed in the literature, larger supercoiling increases the average size of plasmids in the entangled regime. At the same time, we discover that this is accompanied by an unexpected increase of diffusivity. We explain our findings as due to a decrease in inter-plasmids threadings and entanglement

Systems of "living" polymers are ubiquitous in industry and are traditionally realised using surfactants. Here I review the state-of-the-art of living polymers and discuss non-equilibrium extensions that may be realised with advanced synthetic chemistry or DNA functionalised by proteins. These systems are not only interesting in order to realise novel "living" soft matter but can also shed insight into how genomes are (topologically) regulated in vivo., Comment: 6 pages, 4 figures

Topological constraints (TCs) between polymers determine the behaviour of complex fluids such as creams, oils and plastics. Most of the polymer solutions used every day life employ linear chains; their behaviour is accurately captured by the reptation and tube theories which connect microscopic TCs to macroscopic viscoelasticity. On the other hand, polymers with non-trivial topology, such as rings, hold great promise for new technology but pose a challenging problem as they do not obey standard theories; additionally, topological invariance -- i.e. the fact that rings must remain unknotted and unlinked if prepared so -- precludes any serious analytical treatment. Here we propose an unambiguous, parameter-free algorithm to characterise TCs in polymeric solutions and show its power in characterising TCs of entnagled rings. We analyse large-scale molecular dynamics (MD) simulations via persistent homology, a key mathematical tool to extract robust topological information from large datasets. This method allows us to identify ring-specific TCs which we call "homological threadings" (H-threadings) and to connect them to the polymers' behaviour. It also allows us to identify, in a physically appealing and unambiguous way, scale-dependent loops which have eluded precise quantification so far. We discover that while threaded neighbours slowly grow with the rings' length, the ensuing TCs are extensive also in the asymptotic limit. Our proposed method is not restricted to ring polymers and can find broader applications for the study of TCs in generic polymeric materials., Comment: Accepted for publication in Physical Review Research. 10 pages, 4 figures in main text + 2 figures in appendix

We use Brownian dynamics simulations and advanced topological profiling methods to characterize the out-of-equilibrium evolution of self-entanglement in linear polymers confined into nano-channels and under periodic compression. We distinguish two main forms of entanglement, geometrical and topological. The latter is measured by the number of (essential) crossings of the physical knot detected after a suitable bridging of the chain termini. The former is instead measured as the average number of times a linear chain appears to cross itself when viewed under all projections, and is irrespective of the physical knotted state. The key discovery of our work is that these two forms of entanglement are uncoupled and evolve with distinct dynamics. While geometrical entanglement is typically in phase with the compression-elongation cycles and it is primarily sensitive to its force f, the topological measure is mildly sensitive to cyclic modulation but strongly depends on both compression force f and duration k. The findings could assist the interpretation of experiments using fluorescence molecular tracers to track physical knots in polymers. Furthermore, we identify optimal regions in the experimentally-controllable parameter space in which to obtain more/less topological and geometrical entanglement; this may help designing polymers with targeted topology.

The mechanism through which cells determine their fate is intimately related to the spreading of certain biochemical (so-called epigenetic) marks along their genome. The mechanisms behind mark spreading and maintenance are not yet fully understood, and current models often assume a long-range infection-like process for the dynamics of marks, due to the polymeric nature of the chromatin fibre which allows looping between distant sites. While these existing models typically consider antagonising marks, here we propose a qualitatively different scenario which analyses the spreading of a single mark. We define a 1D stochastic model in which mark spreading/infection occurs as a long-range process whereas mark erasure/recovery is a local process, with an enhanced rate at boundaries of infected domains. In the limiting case where our model exhibits absorbing states, we find a first-order-like transition separating the marked/infected phase from the unmarked/recovered phase. This suggests that our model, in this limit, belongs to the long-range compact directed percolation universality class. The abrupt nature of the transition is retained in a more biophysically realistic situation when a basal infection/recovery rate is introduced (thereby removing absorbing states). Close to the transition there is a range of bistability where both the marked/infected and unmarked/recovered states are metastable and long lived, which provides a possible avenue for controlling fate decisions in cells. Increasing the basal infection/recovery rate, we find a second transition between a coherent (marked or unmarked) phase, and a mixed, or random, one., Comment: 11 pages, 7 figures, 2 appendices

Epigenetics is a driving force of important and ubiquitous phenomena in nature such as cell differentiation or even metamorphosis. Oppositely to its widespread role, understanding the biophysical principles that allow epigenetics to control and rewire gene regulatory networks remains an open challenge. In this work we study the effects of epigenetic modifications on the spatial folding of chromosomes -- and hence on the expression of the underlying genes -- by mapping the problem to a class of models known as magnetic polymers. In this work we show that a first-order phase transition underlies the simultaneous spreading of certain epigenetic marks and the conformational collapse of a chromosome. Further, we describe Brownian Dynamics simulations of the model in which the topology of the polymer and thermal fluctuations are fully taken into account and that confirms our mean-field predictions. Extending our models to allow for non-equilibrium terms yields new stable phases which qualitatively agrees with observations in vivo. Our results show that statistical mechanics techniques applied to models of magnetic polymers can be successfully exploited to rationalize the outcomes of experiments designed to probe the interplay between a dynamic epigenetic landscape and chromatin organization., Comment: to appear in PRE

Chromosome organisation is increasingly recognised as an essential component of genome regulation, cell fate and cell health. Within the realm of transposable elements (TEs) however, the spatial information of how genomes are folded is still only rarely integrated in experimental studies or accounted for in modelling. Here, we propose a new predictive modelling framework for the study of the integration patterns of TEs based on extensions of widely employed polymer models for genome organisation. Whilst polymer physics is recognised as an important tool to understand the mechanisms of genome folding, we now show that it can also offer orthogonal and generic insights into the integration and distribution profiles (or "topography") of TEs across organisms. Here, we present polymer physics arguments and molecular dynamics simulations on TEs inserting into heterogeneously flexible polymers and show with a simple model that polymer folding and local flexibility affects TE integration patterns. The preliminary discussion presented herein lay the foundations for a large-scale analysis of TE integration dynamics and topography as a function of the three-dimensional host genome.

Understanding the spatial organisation of the genome in the cell nucleus is one of the current grand challenges in biophysics. Certain biochemical -- or epigenetic -- marks that are deposited along the genome are thought to play an important, yet poorly understood, role in determining genome organisation and cell identity. The physical principles underlying the interplay between epigenetic dynamics and genome folding remain elusive. Here we propose and study a theory that assumes a coupling between epigenetic mark and genome densities, and which can be applied at the scale of the whole nucleus. We show that equilibrium models are not compatible with experiments and a qualitative agreement is recovered by accounting for non-equilibrium processes which can stabilise microphase separated epigenomic domains. We finally discuss the potential biophysical origin of these terms., Comment: To appear in PRL; The quality of figures in this version is downgraded to comply with the arXiv storage limits

We use Brownian dynamics simulations to study the formation of chromatin loops through diffusive sliding of slip-link-like proteins, mimicking the behaviour of cohesin molecules. We recently proposed that diffusive sliding is sufficient to explain the extrusion of chromatin loops of hundreds of kilo-base-pairs (kbp), which may then be stabilised by interactions between cohesin and CTCF proteins. Here we show that the flexibility of the chromatin fibre strongly affects this dynamical process, and find that diffusive loop extrusion is more efficient on stiffer chromatin regions. We also show that the dynamics of loop formation are faster in confined and collapsed chromatin conformations but that this enhancement is counteracted by the increased crowding. We provide a simple theoretical argument explaining why stiffness and collapsed conformations favour diffusive extrusion. In light of the heterogeneous physical and conformational properties of eukaryotic chromatin, we suggest that our results are relevant to understand the looping and organisation of interphase chromosomes in vivo.

The relationship between polymer topology and bulk rheology remains a key question in soft matter physics. Architecture-specific constraints (or threadings) are thought to control the dynamics of ring polymers in ring-linear blends, which thus affects the viscosity to range between that of the pure rings and a value larger, but still comparable to, that of the pure linear melt. Here we consider qualitatively different systems of linear and ring polymers, fused together in "chimeric" architectures. The simplest example of this family is a "tadpole"-shaped polymer - a single ring fused to the end of a single linear chain. We show that polymers with this architecture display a threading-induced dynamical transition that substantially slows chain relaxation. Our findings shed light on how threadings control dynamics and may inform design principles for chimeric polymers with topologically-tunable bulk rheological properties., Comment: Accepted in ACS Macroletters

Polymer composites are ideal candidates for next generation biomimetic soft materials because of their exquisite bottom-up designability. However, the richness of behaviours comes at a price: the need for precise and extensive characterisation of material properties over a highly-dimensional parameter space, as well as a quantitative understanding of the physical principles underlying desirable features. Here we couple large-scale Molecular Dynamics simulations with optical tweezers microrheology to characterise the viscoelastic response of DNA-actin composites. We discover that the previously observed non-monotonic stress-stiffening of these composites is robust, yet tunable, in a broad range of the parameter space that spans two orders of magnitude in DNA length. Importantly, we discover that the most pronounced stiffening is achieved when the species are maximally coupled, i.e. have similar number of entanglements, and not when the number of entanglements per DNA chain is largest. We further report novel dynamical oscillations of the microstructure of the composites, alternating between mixed and bundled phases, opening the door to future investigations. The generic nature of our system renders our results applicable to the behaviour of a broad class of polymer composites., Comment: Accepted in Soft Matter

We present a novel mechanism for resolving the mechanical rigidity of nanoscopic circular polymers that flow in a complex environment. The emergence of a regime of negative differential mobility induced by topological interactions between the rings and the substrate is the key mechanism for selective sieving of circular polymers with distinct flexibilities. A simple model accurately describes the sieving process observed in molecular dynamics simulations and yields experimentally verifiable analytical predictions, which can be used as a reference guide for improving filtration procedures of circular filaments. The topological sieving mechanism we propose ought to be relevant also in probing the microscopic details of complex substrates., Comment: Main text

Topological entanglements severely interfere with important biological processes. For this reason, genomes must be kept unknotted and unlinked during most of a cell cycle. Type II Topoisomerase (TopoII) enzymes play an important role in this process but the precise mechanisms yielding systematic disentanglement of DNA in vivo are not clear. Here we report computational evidence that Structural Maintenance of Chromosomes (SMC) proteins -- such as cohesins and condensins -- can cooperate with TopoII to establish a synergistic mechanism to resolve topological entanglements. SMC-driven loop extrusion (or diffusion) induces the spatial localisation of essential crossings in turn catalysing the simplification of knots and links by TopoII enzymes even in crowded and confined conditions. The mechanism we uncover is universal in that it does not qualitatively depend on the specific substrate, whether DNA or chromatin, or on SMC processivity; we thus argue that this synergy may be at work across organisms and throughout the cell cycle., Comment: main text + SI; accepted in PNAS; supplementary movies can be found at https://www2.ph.ed.ac.uk/~dmichiel/Gallery.html

The genetic instructions stored in the genome require an additional layer of information to robustly determine cell fate. This additional regulation is provided by the interplay between chromosome-patterning biochemical ("epigenetic") marks and three-dimensional genome folding. Yet, the physical principles underlying the dynamical coupling between three-dimensional genomic organisation and one-dimensional epigenetic patterns remain elusive. To shed light on this issue, here we study by mean field theory and Brownian dynamics simulations a magnetic polymer model for chromosomes, where each monomer carries a dynamic epigenetic mark. At the single chromosome level, we show that a first order transition describes the unlimited spreading of epigenetic marks, a phenomenon that is often observed in vivo. At the level of the whole nucleus, experiments suggest chromosomes form micro-phase separated compartments with distinct epigenetic marks. We here discover that for a melt of magnetic polymers such a morphology is thermodynamically unstable, but can be stabilised by a non- equilibrium and ATP-mediated epigenetic switch between different monomer states., Comment: main text + SI; Supplementary Movies can be found at https://www2.ph.ed.ac.uk/~dmichiel/MP.html

Composites of flexible and rigid polymers are ubiquitous in biology and industry alike, yet the physical principles determining their mechanical properties are far from understood. Here, we couple force spectroscopy with large-scale Brownian Dynamics simulations to elucidate the unique viscoelastic properties of custom-engineered blends of entangled flexible DNA molecules and semiflexible actin filaments. We show that composites exhibit enhanced stress-stiffening and prolonged mechano-memory compared to systems of actin or DNA alone, and that these nonlinear features display a surprising non-monotonic dependence on the fraction of actin in the composite. Simulations reveal that these counterintuitive results arise from synergistic microscale interactions between the two biopolymers. Namely, DNA entropically drives actin filaments to form bundles that stiffen the network but reduce the entanglement density, while a percolating actin network is required to reinforce the DNA network against yielding and flow. The competition between bundling and percolation triggers an unexpected stress response that leads equal mass actin-DNA composites to exhibit the most pronounced stress-stiffening and the most long-lived entanglements.

Certain retroviruses, including HIV, insert their DNA in a non-random fraction of the host genome via poorly understood selection mechanisms. Here, we develop a biophysical model for retroviral integrations as stochastic and quasi-equilibrium topological reconnections between polymers. We discover that physical effects, such as DNA accessibility and elasticity, play important and universal roles in this process. Our simulations predict that integration is favoured within nucleosomal and flexible DNA, in line with experiments, and that these biases arise due to competing energy barriers associated with DNA deformations. By considering a long chromosomal region in human T-cells during interphase, we discover that at these larger scales integration sites are predominantly determined by chromatin accessibility. Finally, we propose and solve a reaction-diffusion problem that recapitulates the distribution of HIV hot-spots within T-cells. With few generic assumptions, our model can rationalise experimental observations and identifies previously unappreciated physical contributions to retroviral integration site selection., Comment: Accepted in Nat Comm. SI and Movies can be found at https://www2.ph.ed.ac.uk/~dmichiel/

Reconciling the stability of epigenetic patterns with the rapid turnover of histone modifications and their adaptability to external stimuli is an outstanding challenge. Here, we propose a new biophysical mechanism that can establish and maintain robust yet plastic epigenetic domains via genomic bookmarking (GBM). We model chromatin as a recolourable polymer whose segments bear non-permanent histone marks (or colours) which can be modified by "writer" proteins. The three-dimensional chromatin organisation is mediated by protein bridges, or "readers", such as Polycomb Repressive Complexes and Transcription Factors. The coupling between readers and writers drives spreading of biochemical marks and sustains the memory of local chromatin states across replication and mitosis. In contrast, GBM-targeted perturbations destabilise the epigenetic patterns. Strikingly, we demonstrate that GBM alone can explain the full distribution of Polycomb marks in a whole Drosophila chromosome. We finally suggest that our model provides a starting point for an understanding of the biophysics of cellular differentiation and reprogramming., Comment: Published in Nucleic Acids Research; Supplementary Movies can be found at this url: https://www2.ph.ed.ac.uk/~dmichiel/ShapingEpigeneticMemory.html or https://www.youtube.com/watch?v=xY-DNP58yBU

Elucidating the physics of a concentrated suspension of ring polymers, or of an ensemble of ring polymers in a complex environment, is an important outstanding question in polymer physics. Many of the characteristic features of these systems arise due to topological interactions between polymers, or between the polymers and the environment, and it is often challenging to describe this quantitatively. Here we review recent research which suggests that a key role is played by inter-ring threadings (or penetrations), which become more abundant as the ring size increases. As we discuss, the physical consequences of such threadings are far-reaching: for instance, they lead to a topologically-driven glassy behaviour of ring polymer melts under pinning perturbations, while they can also account for the shape of experimentally observed patterns in two-dimensional gel electrophoresis of DNA knots., Comment: Topical Review; Published in "Polymers" (http://www.mdpi.com/2073-4360/9/8/349). Special issue on "Knotted and Catenated Polymers"

The spreading and regulation of epigenetic marks on chromosomes is crucial to establish and maintain cellular identity. Nonetheless, the dynamical mechanism leading to the establishment and maintenance of a given, cell-line specific, epigenetic pattern is still poorly understood. In this work we propose, and investigate in silico, a possible experimental strategy to illuminate the interplay between 3D chromatin structure and epigenetic dynamics. We consider a set-up where a reconstituted chromatin fibre is stretched at its two ends (e.g., by laser tweezers), while epigenetic enzymes (writers) and chromatin-binding proteins (readers) are flooded into the system. We show that, by tuning the stretching force and the binding affinity of the readers for chromatin, the fibre undergoes a sharp transition between a stretched, epigenetically disordered, state and a crumpled, epigenetically coherent, one. We further investigate the case in which a knot is tied along the chromatin fibre, and find that the knotted segment enhances local epigenetic order, giving rise to "epigenetic solitons" which travel and diffuse along chromatin. Our results point to an intriguing coupling between 3D chromatin topology and epigenetic dynamics, which may be investigated via single molecule experiments., Comment: Accepted version; Supplementary movies can be found at http://www2.ph.ed.ac.uk/~dmichiel/KnottedSolitons.html and https://www.youtube.com/watch?v=Osghh9nEhek

One of the most challenging problems in polymer physics is providing a theoretical description for the behaviour of rings in dense solutions and melts. Although it is nowadays well established that the overall size of a ring in these conditions scales like that of a collapsed globule, there is compelling evidence that rings may exhibit ramified and tree-like conformations. In this work I show how to characterise these local tree-like structures by measuring the local writhing of the rings' segments and by identifying the patterns of intra-chain contacts. These quantities reveal two major topological structures: loops and terminal branches which strongly suggest that the strictly double-folded "lattice animal" picture for rings in the melt may be replaced by a more relaxed tree-like structure accommodating loops. In particular, I show that one can identify hierarchically looped structures whose degree increases linearly with the size of a ring, and that terminal branches are found to store about 30% of the whole ring mass, irrespectively of its length. Finally, I draw an analogy between rings in the melt and slip-linked chains, where contact points are enforced by mobile slip-links and for which a field-theoretic treatment can be employed to get some insight into their typical conformations. These findings are ultimately discussed in the light of recent works on the static structure of rings and on the existence of inter-ring threadings., Comment: Accepted version. To appear in RSC Soft Matter. Suppl. Movie can be found at http://www2.ph.ed.ac.uk/~dmichiel/

One of the most important problems in development is how epigenetic domains can be first established, and then maintained, within cells. To address this question, we propose a framework which couples 3D chromatin folding dynamics, to a "recolouring" process modelling the writing of epigenetic marks. Because many intra-chromatin interactions are mediated by bridging proteins, we consider a "two-state" model with self-attractive interactions between two epigenetic marks which are alike (either active or inactive). This model displays a first-order-like transition between a swollen, epigenetically disordered, phase, and a compact, epigenetically coherent, chromatin globule. If the self-attraction strength exceeds a threshold, the chromatin dynamics becomes glassy, and the corresponding interaction network freezes. By modifying the epigenetic read-write process according to more biologically-inspired assumptions, our polymer model with recolouring recapitulates the ultrasensitive response of epigenetic switches to perturbations, and accounts for multi-domain conformations, strikingly similar to the topologically-associating-domains observed in eukaryotic chromosomes., Comment: Accepted version. To appear in Physical Review X. Combined main text + SI; Suppl. Movies at http://www2.ph.ed.ac.uk/~dmichiel/

Three-dimensional interphase organization of metazoan genomes has been linked to cellular identity. However, the principles governing 3D interphase genome architecture and its faithful transmission through disruptive events of cell-cycle, like mitosis, are not fully understood. By using Brownian dynamics simulations of Drosophila chromosome 3R up to time-scales of minutes, we show that chromatin binding profile of Polycomb-repressive-complex-1 robustly predicts a sub-set of topologically associated domains (TADs), and inclusion of other factors recapitulates the profile of all TADs, as observed experimentally. Our simulations show that chromosome 3R attains interphase organization from mitotic state by a two-step process in which formation of local TADs is followed by long-range interactions. Our model also explains statistical features and tracks the assembly kinetics of polycomb subnuclear clusters. In conclusion, our approach can be used to predict structural and kinetic features of 3D chromosome folding and its associated proteins in biological relevant genomic and time scales., Comment: 20 pages; text + SI; submitted version; supplementary movies can be found at http://www2.ph.ed.ac.uk/~dmichiel/

The static and dynamic properties of ring polymers in concentrated solutions remains one of the last deep unsolved questions in Polymer Physics. At the same time, the nature of the glass transition in polymeric systems is also not well understood. In this work we study a novel glass transition in systems made of circular polymers by exploiting the topological constraints that are conjectured to populate concentrated solutions of rings. We show that such rings strongly inter-penetrate through one another, generating an extensive network of topological interactions that dramatically affects their dynamics. We show that a kinetically arrested state can be induced by randomly pinning a small fraction of the rings. This occurs well above the classical glass transition temperature at which microscopic mobility is lost. Our work demonstrates both the existence of long-lived inter-ring penetrations and also realises a novel, topologically-induced, glass transition., Comment: Accepted version; To appear in PNAS. 13 pages, main text+SI. Supplementary movies can be found at the URL: http://www2.ph.ed.ac.uk/~dmichiel/

Gel electrophoresis is a powerful experimental method to probe the topology of DNA and other biopolymers. While there is a large body of experimental work which allows us to accurately separate different topoisomers of a molecule, a full theoretical understanding of these experiments has not yet been achieved. Here we show that the mobility of DNA knots depends crucially and subtly on the physical properties of the gel, and in particular on the presence of dangling ends. The topological interactions between these and DNA molecules can be described in terms of an "entanglement number", and yield a non-monotonic mobility at moderate fields. Consequently, in two-dimensional electrophoresis, gel bands display a characteristic arc pattern; this turns into a straight line when the density of dangling ends vanishes. We also provide a novel framework to accurately predict the shape of such arcs as a function of molecule length and topological complexity, which may be used to inform future experiments., Comment: To appear in PNAS; 12 pages: text+SI

In this paper we study the role of topology in DNA gel electrophoresis experiments via molecular dynamics simulations. The gel is modelled as a 3D array of obstacles from which half edges are removed at random with probability p, thereby generating a disordered environment. Changes in the microscopic structure of the gel are captured by measuring the electrophoretic mobility of ring polymers moving through the medium, while their linear counterparts provide a control system as we show they are insensitive to these changes. We show that ring polymers provide a novel non-invasive way of exploiting topology to sense microscopic disorder. Finally, we compare the results from the simulations with an analytical model for the non-equilibrium differential mobility, and find a striking agreement between simulation and theory, Comment: 8 pages,6 figures Accepted for publication in Soft Matter

We study theoretically the replication of Kinetoplast DNA consisting of several thousands separate mini-circles found in organisms of the class Kinetoplastida. When the cell is not actively dividing these are topologically connected in a marginally linked network of rings with only one connected component. During cell division each mini-circle is removed from the network, duplicated and then re-attached, along with its progeny. We study this process under the hypothesis that there is a coupling between the topological state of the mini-circles and the expression of genetic information encoded on them, leading to the production of Topoisomerase. This model describes a self-regulating system capable of full replication that reproduces several previous experimental findings. We find that the fixed point of the system depends on a primary free parameter of the model: the ratio between the rate of removal of mini-circles from the network (R) and their (re)attachment rate (A). The final topological state is found to be that of a marginally linked network structure in which the fraction of mini-circles linked to the largest connected component approaches unity as R/A decreases. Finally we discuss how this may suggest an evolutionary trade-off between the speed of replication and the accuracy with which a fully topologically linked state is produced., Comment: 14 pages, 6 figures

In this work we present a computational study of the Kinetoplast genome, modelled as a large number of semiflexible unknotted loops, which are allowed to link with each other. As the DNA density increases, the systems shows a percolation transition between a gas of unlinked rings and a network of linked loops which spans the whole system. Close to the percolation transition, we find that the mean valency of the network, i.e. the average number of loops which are linked to any one loop, is around 3 as found experimentally for the Kinetoplast DNA. Even more importantly, by simulating the digestion of the network by a restriction enzyme, we show that the distribution of oligomers, i.e. structures formed by a few loops which remain linked after digestion, quantitatively matches experimental data obtained from gel electrophoresis, provided that the density is, once again, close to the percolation transition. With respect to previous work, our analysis builds on a reduced number of assumptions, yet can still fully explain the experimental data. Our findings suggest that the Kinetoplast DNA can be viewed as a network of linked loops positioned very close to the percolation transition, and we discuss the possible biological implications of this remarkable fact., Comment: 13p, 5f, merged SI

We study of the dynamics of ring polymers confined to diffuse in a background gel at low concentrations. We do this in order to probe the inter-play between topology and dynamics in ring polymers. We develop an algorithm that takes into account the possibility that the rings hinder their own motion by passing through themselves, i.e. "self-threading". Our results suggest that the number of self-threadings scales extensively with the length of the rings and that this is substantially independent of the details of the model. The slowing down of the rings' dynamics is found to be related to the fraction of segments that can contribute to the motion. Our results give a novel perspective on the motion of ring polymers in gel, for which a complete theory is still lacking, and may help us to understand the irreversible trapping of ring polymers in gel electrophoresis experiments., Comment: 11 pages,14 figures, Soft Matter 2014