Search

The insertion of DNA elements within genomes underpins both genetic diversity and disease when unregulated. Most of DNA insertions are not random and 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 viral DNA or transposons, into naked DNA or chromatin substrates. More specifically, we explore the role of loops within the polymeric substrate and discover that they act as 'geometric catalysts' for DNA integration by reducing the energy barrier for substrate deformation. Additionally, we discover that the 1D pattern and 3D conformation of loops have a marked effect on the distribution of integration sites. Finally, we show that loops may compete with nucleosomes to attract 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 engineering., (© The Author(s) 2024. Published by Oxford University Press on behalf of Nucleic Acids Research.)

In this paper we investigate the effects of varying cation valency and concentration on the rheology of entangled λ DNA solutions. We show that monovalent cations moderately increase the viscoelasticty of the solutions mainly by stabilising linear concatenation of λ 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 Mg 2+ . 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.

Transcription termination pathways mitigate the detrimental consequences of unscheduled promiscuous initiation occurring at hundreds of thousands of genomic cis -regulatory elements. The Restrictor complex, composed of the Pol II-interacting protein WDR82 and the RNA-binding protein ZC3H4, suppresses processive transcription at thousands of extragenic sites in mammalian genomes. Restrictor-driven termination does not involve nascent RNA cleavage, and its interplay with other termination machineries is unclear. Here we show that efficient termination at Restrictor-controlled extragenic transcription units involves the recruitment of the protein phosphatase 1 (PP1) regulatory subunit PNUTS, a negative regulator of the SPT5 elongation factor, and Symplekin, a protein associated with RNA cleavage complexes but also involved in cleavage-independent and phosphatase-dependent termination of noncoding RNAs in yeast. PNUTS and Symplekin act synergistically with, but independently from, Restrictor to dampen processive extragenic transcription. Moreover, the presence of limiting nuclear levels of Symplekin imposes a competition for its recruitment among multiple transcription termination machineries, resulting in mutual regulatory interactions. Hence, by synergizing with Restrictor, Symplekin and PNUTS enable efficient termination of processive, long-range extragenic transcription., (© 2023 Russo et al.; Published by Cold Spring Harbor Laboratory Press.)

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.

In most bacteria, chromosome segregation is driven by the ParABS system where the CTPase protein ParB loads at the parS site to trigger the formation of a large partition complex. Here, we present in vitro studies of the partition complex for Bacillus subtilis ParB, using single-molecule fluorescence microscopy and AFM imaging to show that transient ParB-ParB bridges are essential for forming DNA condensates. Molecular Dynamics simulations confirm that condensation occurs abruptly at a critical concentration of ParB and show that multimerization is a prerequisite for forming the partition complex. Magnetic tweezer force spectroscopy on mutant ParB proteins demonstrates that CTP hydrolysis at the N-terminal domain is essential for DNA condensation. Finally, we show that transcribing RNA polymerases can steadily traverse the ParB-DNA partition complex. These findings uncover how ParB forms a stable yet dynamic partition complex for chromosome segregation that induces DNA condensation and segregation while enabling replication and transcription., (© The Author(s) 2023. Published by Oxford University Press on behalf of Nucleic Acids Research.)

In spite of the nanoscale and single-molecule insights into nucleoid associated proteins (NAPs), their role in modulating the mesoscale viscoelasticity of entangled DNA has been overlooked so far. By combining microrheology and molecular dynamics simulation, we find that the abundant NAP "integration host factor" (IHF) lowers the viscosity of entangled λDNA 20-fold at physiological concentrations and stoichiometries. Our results suggest that IHF may play a previously unappreciated role in resolving DNA entanglements and in turn may be acting as a "genomic fluidizer" for bacterial genomes.

DNA recombination is a ubiquitous process that ensures genetic diversity. Contrary to textbook pictures, DNA recombination, as well as generic DNA translocations, occurs in a confined and highly entangled environment. Inspired by this observation, here, we investigate a solution of semiflexible polymer rings undergoing generic cutting and reconnection operations under spherical confinement. Our setup may be realized using engineered DNA in the presence of recombinase proteins or by considering micelle-like components able to form living (or reversibly breakable) polymer rings. We find that in such systems, there is a topological gelation transition, which can be triggered by increasing either the stiffness or the concentration of the rings. Flexible or dilute polymers break into an ensemble of short, unlinked, and segregated rings, whereas sufficiently stiff or dense polymers self-assemble into a network of long, linked, and mixed loops, many of which are knotted. We predict that the two phases should behave qualitatively differently in elution experiments monitoring the escape dynamics from a permeabilized container. Besides shedding some light on the biophysics and topology of genomes undergoing DNA reconnection in vivo, our findings could be leveraged in vitro to design polymeric complex fluids-e.g., DNA-based complex fluids or living polymer networks-with desired topologies.

Inspired by how certain proteins "sense" knots and entanglements in DNA molecules, here, we ask if local geometric features that may be used as a readout of the underlying topology of generic polymers exist. We perform molecular simulations of knotted and linked semiflexible polymers and study four geometric measures to predict topological entanglements: local curvature, local density, local 1D writhe, and nonlocal 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 delocalize essential crossings in knotted polymers. We further discover that nonlocal 3D writhe is the best geometric readout of the knot location. Finally, we discuss how these geometric features may be used to computationally analyze entanglements in generic polymer melts and gels., Competing Interests: The authors declare no competing financial interest., (© 2022 The Authors. Published by American Chemical Society.)

The ParAB S system is essential for prokaryotic chromosome segregation. After loading at parS on the genome, ParB (partition protein B) proteins rapidly redistribute to distances of ~15 kilobases from the loading site. It has remained puzzling how this large-distance spreading can occur along DNA loaded with hundreds of proteins. Using in vitro single-molecule fluorescence imaging, we show that ParB from Bacillus subtilis can load onto DNA distantly of parS , as loaded ParB molecules themselves are found to be able to recruit additional ParB proteins from bulk. Notably, this recruitment can occur in cis but also in trans, where, at low tensions within the DNA, newly recruited ParB can bypass roadblocks as it gets loaded to spatially proximal but genomically distant DNA regions. The data are supported by molecular dynamics simulations, which show that cooperative ParB-ParB recruitment can enhance spreading. ParS -independent recruitment explains how ParB can cover substantial genomic distance during chromosome segregation, which is vital for the bacterial cell cycle.

In a recent issue of Science, Gabriele et al. have, for the first time, quantified the dynamics of a topologically associating domain (TAD) in live cells by coupling super-resolution imaging and computational modelling, concluding that a TAD spends most of its life in a "partially extruded state" and that CTCF-CTCF loops are rare., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2022 Elsevier Inc. All rights reserved.)

Phase separation is as familiar as watching vinegar separating from oil in vinaigrette. 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 indicates that certain proteins and nucleic acids condensates are not simple liquids and instead display both viscous and elastic behaviors, 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., Competing Interests: The authors declare no competing financial interest., (© 2022 The Authors. Published by American Chemical Society.)

How type 2 Topoisomerase (TopoII) proteins relax and simplify the topology of DNA molecules is one of the most intriguing open questions in genome and DNA biophysics. Most of the existing models neglect the dynamics of TopoII which is expected of 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., (© The Author(s) 2022. Published by Oxford University Press on behalf of Nucleic Acids Research.)

The condensin SMC protein complex organizes chromosomal structure by extruding loops of DNA. Its ATP-dependent motor mechanism remains unclear but likely involves steps associated with large conformational changes within the ∼50 nm protein complex. Here, using high-resolution magnetic tweezers, we resolve single steps in the loop extrusion process by individual yeast condensins. The measured median step sizes range between 20-40 nm at forces of 1.0-0.2 pN, respectively, comparable with the holocomplex size. These large steps show that, strikingly, condensin typically reels in DNA in very sizeable amounts with ∼200 bp on average per single extrusion step at low force, and occasionally even much larger, exceeding 500 bp per step. Using Molecular Dynamics simulations, we demonstrate that this is due to the structural flexibility of the DNA polymer at these low forces. Using ATP-binding-impaired and ATP-hydrolysis-deficient mutants, we find that ATP binding is the primary step-generating stage underlying DNA loop extrusion. We discuss our findings in terms of a scrunching model where a stepwise DNA loop extrusion is generated by an ATP-binding-induced engagement of the hinge and the globular domain of the SMC complex., (© The Author(s) 2021. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Polymer topology has been shown to play a key role in tuning the dynamics of complex fluids and gels. At the same time, polymer composites, ubiquitous in everyday life, have been shown to exhibit emergent desirable mechanical properties not attainable in single-species systems. Yet, how topology impacts the dynamics and structure of polymer composites remains poorly understood. Here, we create composites of rigid rods (microtubules) polymerized within entangled solutions of flexible linear and ring polymers (DNA) of equal length. We couple optical tweezers microrheology with confocal microscopy and scaled particle theory to show that composites with linear DNA exhibit a strongly nonmonotonic dependence of elasticity and stiffness on microtubule concentration due to depletion-driven polymerization and flocculation of microtubules. In contrast, composites containing ring DNA show a much more modest monotonic increase in elastic strength with microtubule concentration, which we demonstrate arises from the decreased conformational size and increased miscibility of rings.

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 their 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 the absence of tethering and a third type of Z-loop, which has not yet been identified in experiments, to appear. Our model naturally explains roadblock bypassing and the appearance of steps larger than the SMC size as a consequence of non-contiguous DNA grabbing. Finally, this study 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., (Copyright © 2021 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

Organization of the genome into transcriptionally active euchromatin and silenced heterochromatin is essential for eukaryotic cell function. Phase-separation has been implicated in heterochromatin formation, but it is unclear how phase-separated condensates can contribute to stable repression, particularly for heritable epigenetic changes. Polycomb complex PRC1 is key for heterochromatin formation, but the multitude of Polycomb proteins has hindered our understanding of their collective contribution to chromatin repression. Here, we show that PRC1 forms multicomponent condensates through hetero-oligomerization. They preferentially seed at H3K27me3 marks, and subsequently write H2AK119Ub marks. We show that inducing Polycomb phase-separation can cause chromatin compaction, but polycomb condensates are dispensable for maintenance of the compacted state. Our data and simulations are consistent with a model in which the time integral of Polycomb phase-separation is progressively recorded in repressive histone marks, which subsequently drive compaction. These findings link the equilibrium thermodynamics of phase-separation with the fundamentally non-equilibrium concept of epigenetic memory., (© 2021. The Author(s).)

Ring polymers in dense solutions are among the most intriguing problems in polymer physics. Because of its natural occurrence in circular form, DNA has been extensively used as a proxy to study the fundamental physics of ring polymers in different topological states. Yet, torsionally constrained-such as supercoiled-topologies have been largely neglected so far. The applicability of existing theoretical models to dense supercoiled DNA is thus unknown. Here, we address this gap by coupling large-scale molecular dynamics simulations with differential dynamic microscopy of entangled supercoiled DNA plasmids. We find that, unexpectedly, larger supercoiling increases the size of entangled plasmids and concomitantly induces an enhancement in DNA mobility. These findings are reconciled as due to supercoiling-driven asymmetric and double-folded plasmid conformations that reduce interplasmid entanglements and threadings. Our results suggest a way to topologically tune DNA mobility via supercoiling, thus enabling topological control over the (micro)rheology of DNA-based complex fluids., (Copyright © 2021 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works. Distributed under a Creative Commons Attribution NonCommercial License 4.0 (CC BY-NC).)

Understanding biological function requires the identification and characterization of complex patterns of molecules. Single-molecule localization microscopy (SMLM) can quantitatively measure molecular components and interactions at resolutions far beyond the diffraction limit, but this information is only useful if these patterns can be quantified and interpreted. We provide a new approach for the analysis of SMLM data that develops the concept of structures and super-structures formed by interconnected elements, such as smaller protein clusters. Using a formal framework and a parameter-free algorithm, (super-)structures formed from smaller components are found to be abundant in classes of nuclear proteins, such as heterogeneous nuclear ribonucleoprotein particles (hnRNPs), but are absent from ceramides located in the plasma membrane. We suggest that mesoscopic structures formed by interconnected protein clusters are common within the nucleus and have an important role in the organization and function of the genome. Our algorithm, SuperStructure, can be used to analyze and explore complex SMLM data and extract functionally relevant information., (© 2021 Marenda et al.)

We study the effect of transcription on the kinetics of DNA supercoiling in three dimensions by means of Brownian dynamics simulations of a single-nucleotide-resolution coarse-grained model for double-stranded DNA. By explicitly accounting for the action of a transcribing RNA polymerase (RNAP), we characterize the geometry and nonequilibrium dynamics of the ensuing twin supercoiling domains. Contrary to the typical textbook picture, we find that the generation of twist by RNAP results in the formation of plectonemes (writhed DNA) some distance away. We further demonstrate that this translates into an "action at a distance" on DNA-binding proteins; for instance, positive supercoils downstream of an elongating RNAP destabilize nucleosomes long before the transcriptional machinery reaches the histone octamer. We also analyze the relaxation dynamics of supercoiled double-stranded DNA, and characterize the widely different timescales of twist diffusion, which is a simple and fast process, and writhe relaxation, which is much slower and entails multiple steps., Competing Interests: The authors declare no competing interest.

Structural maintenance of chromosome (SMC) protein complexes are able to extrude DNA loops. While loop extrusion constitutes a fundamental building block of chromosomes, other factors may be equally important. Here, we show that yeast cohesin exhibits pronounced clustering on DNA, with all the hallmarks of biomolecular condensation. DNA-cohesin clusters exhibit liquid-like behavior, showing fusion of clusters, rapid fluorescence recovery after photobleaching and exchange of cohesin with the environment. Strikingly, the in vitro clustering is DNA length dependent, as cohesin forms clusters only on DNA exceeding 3 kilo-base pairs. We discuss how bridging-induced phase separation, a previously unobserved type of biological condensation, can explain the DNA-cohesin clustering through DNA-cohesin-DNA bridges. We confirm that, in yeast cells in vivo, a fraction of cohesin associates with chromatin in a manner consistent with bridging-induced phase separation. Biomolecular condensation by SMC proteins constitutes a new basic principle by which SMC complexes direct genome organization., (Copyright © 2021 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works. Distributed under a Creative Commons Attribution License 4.0 (CC BY).)

Understanding how entanglements affect the behavior 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 spatiotemporal entanglement structures from the pairwise 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 spatiotemporal structures of dynamical entanglement in polymeric solutions. We also propose a mean-field approximation of our analysis that provides previously under-appreciated physical insights into the dynamics of generic entangled polymers. We envisage DEA will be useful to analyze the dynamical evolution of entanglements in generic polymeric systems such as blends and composites.

The study and the investigation of structural and dynamical properties of complex systems have attracted considerable interest among scientists in general and physicists and biologists in particular. The present review paper represents a broad overview of the research performed over the nonlinear dynamics of DNA, devoted to some different aspects of DNA physics and including analytical, quantum and computational tools to understand nonlinear DNA physics. We review in detail the semi-discrete approximation within helicoidal Peyrard-Bishop model and show that localized modulated solitary waves, usually called breathers, can emerge and move along the DNA. Since living processes occur at submolecular level, we then discuss a quantum treatment to address the problem of how charge and energy are transported on DNA and how they may play an important role for the functioning of living cells. While this problem has attracted the attention of researchers for a long time, it is still poorly understood how charge and energy transport can occur at distances comparable to the size of macromolecules. Here, we review a theory based on the mechanism of 'self-trapping' of electrons due to their interaction with mechanical (thermal) oscillation of the DNA structure. We also describe recent computational models that have been developed to capture nonlinear mechanics of DNA in vitro and in vivo , possibly under topological constraints. Finally, we provide some conjectures on potential future directions for this field., Competing Interests: We declare we have no competing interests., (© 2020 The Authors.)

Systems of "living" polymers are ubiquitous in industry and are traditionally realised using surfactants. Here I first review the theoretical state-of-the-art of living polymers and then 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.

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 nanochannels and under periodic compression. By introducing suitable observables, we can distinguish two main forms of entanglement that we term 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 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., Competing Interests: The authors declare no competing financial interest.

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 infectionlike process for the dynamics of marks, due to the polymeric nature of the chromatin fiber which allows looping between distant sites. While these existing models typically consider antagonizing marks, here we propose a qualitatively different scenario which analyses the spreading of a single mark. We define a one-dimensional 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.

Understanding the spatial organization 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 organization 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 nonequilibrium processes that can stabilize microphase separated epigenomic domains. We finally discuss the potential biophysical origin of these terms.

Epigenetics is a driving force of important and ubiquitous phenomena in nature such as cell differentiation or even metamorphosis. Opposite 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 confirm our mean field predictions. Extending our models to allow for nonequilibrium 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.

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 organization. Our model explains the conventional organization of heterochromatin and euchromatin in growing cells and the pathological organization 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., (Copyright © 2019 The Authors. Published by Elsevier Inc. All rights reserved.)

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.

The importance of three-dimensional chromatin organisation in genome regulation has never been clearer. But in spite of the enormous technological advances to probe chromatin organisation in vivo, there is still a lack of mechanistic understanding of how such an arrangement is achieved. Here we review emerging evidence pointing to an intriguing role of nuclear RNA in shaping large-scale chromatin structure and regulating genome function. We suggest this role may be achieved through the formation of a dynamic nuclear mesh that can exploit ATP-driven processes and phase separation of RNA-binding proteins to tune its assembly and material properties., (Copyright © 2019 The Authors. Published by Elsevier Ltd.. All rights reserved.)

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 localization of essential crossings, in turn catalyzing 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., Competing Interests: The authors declare no conflict of interest., (Copyright © 2019 the Author(s). Published by PNAS.)

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 mechanomemory compared to systems of actin or DNA alone, and that these nonlinear features display a surprising nonmonotonic 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 uniform well-connected actin network is required to reinforce the DNA network against yielding and flow. The competition between bundling and connectivity triggers an unexpected stress response that leads equal mass DNA-actin composites to exhibit the most pronounced stress stiffening and the most long-lived entanglements.

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.

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., (© The Author(s) 2017. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Understanding how topological constraints affect the dynamics of polymers in solution is at the basis of any polymer theory and it is particularly needed for melts of rings. These polymers fold as crumpled and space-filling objects and, yet, they display a large number of topological constraints. To understand their role, here we systematically probe the response of solutions of rings at various densities to "random pinning" perturbations. We show that these perturbations trigger non-Gaussian and heterogeneous dynamics, eventually leading to nonergodic and glassy behavior. We then derive universal scaling relations for the values of solution density and polymer length marking the onset of vitrification in unperturbed solutions. Finally, we directly connect the heterogeneous dynamics of the rings with their spatial organization and mutual interpenetration. Our results suggest that deviations from the typical behavior observed in systems of linear polymers may originate from architecture-specific (threading) topological constraints.

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., Competing Interests: The authors declare no conflict of interest.

Fluorescence microscopy reveals that the contents of many (membrane-free) nuclear bodies exchange rapidly with the soluble pool while the underlying structure persists; such observations await a satisfactory biophysical explanation. To shed light on this, we perform large-scale Brownian dynamics simulations of a chromatin fiber interacting with an ensemble of (multivalent) DNA-binding proteins able to switch between an "on" (binding) and an "off" (nonbinding) state. This system provides a model for any DNA-binding protein that can be posttranslationally modified to change its affinity for DNA (e.g., through phosphorylation). Protein switching is a nonequilibrium process, and it leads to the formation of clusters of self-limiting size, where individual proteins in a cluster exchange with the soluble pool with kinetics similar to those seen in photobleaching experiments. This behavior contrasts sharply with that exhibited by nonswitching proteins, which are permanently in the on-state; when these bind to DNA nonspecifically, they form clusters that grow indefinitely in size. To explain these findings, we propose a mean-field theory from which we obtain a scaling relation between the typical cluster size and the protein switching rate. Protein switching also reshapes intrachromatin contacts to give networks resembling those seen in topologically associating domains, as switching markedly favors local (short-range) contacts over distant ones. Our results point to posttranslational modification of chromatin-bridging proteins as a generic mechanism driving the self-assembly of highly dynamic, nonequilibrium, protein clusters with the properties of nuclear bodies., (Copyright © 2017 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

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 interpenetrate 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 both demonstrates the existence of long-lived inter-ring penetrations and realizes a novel, topologically induced, glass transition.

Gel electrophoresis is a powerful experimental method to probe the topology of DNA and other biopolymers. Although there is a large body of experimental work that 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 nonmonotonic mobility at moderate fields. Consequently, in 2D 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.

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 three, as found experimentally for the kinetoplast DNA (kDNA). 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 kDNA 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.

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.

We study 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 a 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.

We perform large scale three-dimensional molecular dynamics simulations of unlinked and unknotted ring polymers diffusing through a background gel, here a three-dimensional cubic lattice. Taking advantage of this architecture, we propose a new method to unambiguously identify and quantify inter-ring threadings (penetrations) and to relate these to the dynamics of the ring polymers. We find that both the number and the persistence time of the threadings increase with the length of the chains, ultimately leading to a percolating network of inter-ring penetrations. We discuss the implications of these findings for the possible emergence of a topological jammed state of very long rings.

Gap junctions, also referred to as electrical synapses, are expressed along the entire central nervous system and are important in mediating various brain rhythms in both normal and pathological states. These connections can form between the dendritic trees of individual cells. Many dendrites express membrane channels that confer on them a form of sub-threshold resonant dynamics. To obtain insight into the modulatory role of gap junctions in tuning networks of resonant dendritic trees, we generalise the "sum-over-trips" formalism for calculating the response function of a single branching dendrite to a gap junctionally coupled network. Each cell in the network is modelled by a soma connected to an arbitrary structure of dendrites with resonant membrane. The network is treated as a single extended tree structure with dendro-dendritic gap junction coupling. We present the generalised "sum-over-trips" rules for constructing the network response function in terms of a set of coefficients defined at special branching, somatic and gap-junctional nodes. Applying this framework to a two-cell network, we construct compact closed form solutions for the network response function in the Laplace (frequency) domain and study how a preferred frequency in each soma depends on the location and strength of the gap junction.