Topological gelation of reconnecting polymers.

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.

Topological Spectra and Entropy of Chromatin Loop Networks.

The 3D folding of a mammalian gene can be studied by a polymer model, where the chromatin fiber is represented by a semiflexible polymer which interacts with multivalent proteins, representing complexes of DNA-binding transcription factors and RNA polymerases. This physical model leads to the natural emergence of clusters of proteins and binding sites, accompanied by the folding of chromatin into a set of topologies, each associated with a different network of loops. Here, we combine numerics and analytics to first classify these networks and then find their relative importance or statistical weight, when the properties of the underlying polymer are those relevant to chromatin. Unlike polymer networks previously studied, our chromatin networks have finite average distances between successive binding sites, and this leads to giant differences between the weights of topologies with the same number of edges and nodes but different wiring. These weights strongly favor rosettelike structures with a local cloud of loops with respect to more complicated nonlocal topologies. Our results suggest that genes should overwhelmingly fold into a small fraction of all possible 3D topologies, which can be robustly characterized by the framework we propose here.

Being Heterogeneous Is Advantageous: Extreme Brownian Non-Gaussian Searches.

Redundancy in biology may be explained by the need to optimize extreme searching processes, where one or few among many particles are requested to reach the target like in human fertilization. We show that non-Gaussian rare fluctuations in Brownian diffusion dominates such searches, introducing drastic corrections to the known Gaussian behavior. Our demonstration entails different physical systems and pinpoints the relevance of diversity within redundancy to boost fast targeting. We sketch an experimental context to test our results: polydisperse systems.

Folding kinetics of an entangled protein.

The possibility of the protein backbone adopting lasso-like entangled motifs has attracted increasing attention. After discovering the surprising abundance of natively entangled protein domain structures, it was shown that misfolded entangled subpopulations might become thermosensitive or escape the homeostasis network just after translation. To investigate the role of entanglement in shaping folding kinetics, we introduce a novel indicator and analyze simulations of a coarse-grained, structure-based model for two small single-domain proteins. The model recapitulates the well-known two-state folding mechanism of a non-entangled SH3 domain. However, despite its small size, a natively entangled antifreeze RD1 protein displays a rich refolding behavior, populating two distinct kinetic intermediates: a short-lived, entangled, near-unfolded state and a longer-lived, non-entangled, near-native state. The former directs refolding along a fast pathway, whereas the latter is a kinetic trap, consistently with known experimental evidence of two different characteristic times. Upon trapping, the natively entangled loop folds without being threaded by the N-terminal residues. After trapping, the native entangled structure emerges by either backtracking to the unfolded state or threading through the already formed but not yet entangled loop. Along the fast pathway, trapping does not occur because the native contacts at the closure of the lasso-like loop fold after those involved in the N-terminal thread, confirming previous predictions. Despite this, entanglement may appear already in unfolded configurations. Remarkably, a longer-lived, near-native intermediate, with non-native entanglement properties, recalls what was observed in cotranslational folding.

Dynamic and facilitated binding of topoisomerase accelerates topological relaxation.

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.

Synergy of topoisomerase and structural-maintenance-of-chromosomes proteins creates a universal pathway to simplify genome topology.

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.

The rise and fall of branching: A slowing down mechanism in relaxing wormlike micellar networks.

A mean-field kinetic model suggests that the relaxation dynamics of wormlike micellar networks is a long and complex process due to the problem of reducing the number of free end-caps (or dangling ends) while also reaching an equilibrium level of branching after an earlier overgrowth. The model is validated against mesoscopic molecular dynamics simulations and is based on kinetic equations accounting for scission and synthesis processes of blobs of surfactants. A long relaxation time scale is reached with both thermal quenches and small perturbations of the system. The scaling of this relaxation time is exponential with the free energy of an end cap and with the branching free energy. We argue that the subtle end-recombination dynamics might yield effects that are difficult to detect in rheology experiments, with possible underestimates of the typical time scales of viscoelastic fluids.

Aging of living polymer networks: a model with patchy particles.

Microrheology experiments show that viscoelastic media composed by wormlike micellar networks display complex relaxations lasting seconds even at the scale of micrometers. By mapping a model of patchy colloids with suitable mesoscopic elementary motifs to a system of worm-like micelles, we are able to simulate its relaxation dynamics, upon a thermal quench, spanning many decades, from microseconds up to tens of seconds. After mapping the model to real units and to experimental scission energies, we show that the relaxation process develops through a sequence of non-local and energetically challenging arrangements. These adjustments remove undesired structures formed as a temporary energetic solution for stabilizing the thermodynamically unstable free caps of the network. We claim that the observed scale-free nature of this stagnant process may complicate the correct quantification of experimentally relevant time scales as the Weissenberg number.

Shaping epigenetic memory via genomic bookmarking.

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.

Nonequilibrium Theory of Epigenomic Microphase Separation in the Cell Nucleus.

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.

Spatial confinement induces hairpins in nicked circular DNA.

In living cells, DNA is highly confined in space with the help of condensing agents, DNA binding proteins and high levels of supercoiling. Due to challenges associated with experimentally studying DNA under confinement, little is known about the impact of spatial confinement on the local structure of the DNA. Here, we have used well characterized slits of different sizes to collect high resolution atomic force microscopy images of confined circular DNA with the aim of assessing the impact of the spatial confinement on global and local conformational properties of DNA. Our findings, supported by numerical simulations, indicate that confinement imposes a large mechanical stress on the DNA as evidenced by a pronounced anisotropy and tangent-tangent correlation function with respect to non-constrained DNA. For the strongest confinement we observed nanometer sized hairpins and interwound structures associated with the nicked sites in the DNA sequence. Based on these findings, we propose that spatial DNA confinement in vivo can promote the formation of localized defects at mechanically weak sites that could be co-opted for biological regulatory functions.

KymoKnot: A web server and software package to identify and locate knots in trajectories of linear or circular polymers.

The KymoKnot software package and web server identifies and locates physical knots or proper knots in a series of polymer conformations. It is mainly intended as an analysis tool for trajectories of linear or circular polymers, but it can be used on single instances too, e.g. protein structures in PDB format. A key element of the software package is the so-called minimally interfering chain closure algorithm that is used to detect physical knots in open chains and to locate the knotted region in both open and closed chains. The web server offers a user-friendly graphical interface that identifies the knot type and highlights the knotted region on each frame of the trajectory, which the user can visualize interactively from various viewpoints. The dynamical evolution of the knotted region along the chain contour is presented as a kymograph. All data can be downloaded in text format. The KymoKnot package is licensed under the BSD 3-Clause licence. The server is publicly available at http://kymoknot.sissa.it/kymoknot/interactive.php .

Topological patterns in two-dimensional gel electrophoresis of DNA knots.

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.

Non-monotonic knotting probability and knot length of semiflexible rings: the competing roles of entropy and bending energy.

We consider self-avoiding rings of up to 1000 beads and study, by Monte Carlo techniques, how their equilibrium knotting properties depend on the bending rigidity. When the rings are taken from the rigid to fully-flexible limit, their average compactness increases, as expected. However, this progressive compactification is not parallelled by a steady increase of the abundance of knots. In fact the knotting probability, Pk, has a prominent maximum when the persistence length is a few times larger than the bead size. At similar bending rigidities, the knot length has, instead, a minimum. We show that the observed non-monotonicity of Pk arises from the competition between two effects. The first one is the entropic cost of introducing a knot. The second one is the gain in bending energy due to the presence of essential crossings. These, in fact, constrain the knotted region and keep it less bent than average. The two competing effects make knots maximally abundant when the persistence length is 5-10 times larger than the bead size. At such intermediate bending rigidities, knots in the chains of 500 and 1000 beads are 40 times more likely than in the fully-flexible limit.

Sorting ring polymers by knot type with modulated nanochannels.

In this theoretical study we discuss a novel method for sorting ring polymers according to their topological, knotted state. The proposed approach harnesses the rich dynamical behaviour of polymers confined inside spatially-modulated nanochannels. The longitudinal mobility of the rings is shown to have two key properties that are ideally suited for knot sorting. First, at fixed topology, the mobility has an intriguing oscillatory dependence on chain length. Second, the mobility ranking of different knot types is inverted upon increasing the chain length. We show that this complex interplay of channel geometry, chain length and topology can be rationalised within a simple theoretical framework based on Fick-Jacobs's diffusive theory. The results and the interpretative scheme ought to be useful for designing microfluidic devices with optimal topological sorting capabilities.

Single-File Escape of Colloidal Particles from Microfluidic Channels.

Single-file diffusion is a ubiquitous physical process exploited by living and synthetic systems to exchange molecules with their environment. It is paramount to quantify the escape time needed for single files of particles to exit from constraining synthetic channels and biological pores. This quantity depends on complex cooperative effects, whose predominance can only be established through a strict comparison between theory and experiments. By using colloidal particles, optical manipulation, microfluidics, digital microscopy, and theoretical analysis we uncover the self-similar character of the escape process and provide closed-formula evaluations of the escape time. We find that the escape time scales inversely with the diffusion coefficient of the last particle to leave the channel. Importantly, we find that at the investigated microscale, bias forces as tiny as 10^{-15} N determine the magnitude of the escape time by drastically reducing interparticle collisions. Our findings provide crucial guidelines to optimize the design of micro- and nanodevices for a variety of applications including drug delivery, particle filtering, and transport in geometrical constrictions.

Curvature-driven positioning of Turing patterns in phase-separating curved membranes.

We introduce a new finite difference scheme to study the dynamics of Turing patterns of a two-species activator-inhibitor system embedded on a phase-separating curved membrane, modelling for instance a lipid bilayer. We show that the underlying binary fluid can strongly affect both the dynamical and the steady state properties of the ensuing Turing patterns. Furthermore, geometry plays a key role, as a large enough local membrane curvature can both arrest the coarsening of the lipid domains and position the patterns selectively at areas of high or small local curvature. The physical phenomena we observe are due to a minimal coupling, between the diffusivity of the Turing components and the local membrane composition. While our study is theoretical in nature, it can provide a framework within which to address intracellular pattern formation in systems of interacting membrane proteins.

Correction: Curvature-driven positioning of Turing patterns in phase-separating curved membranes.

Correction for 'Curvature-driven positioning of Turing patterns in phase-separating curved membranes' by Giulio Vandin et al., Soft Matter, 2016, DOI: .

Topological friction strongly affects viral DNA ejection.

Bacteriophages initiate infection by releasing their double-stranded DNA into the cytosol of their bacterial host. However, what controls and sets the timescales of DNA ejection? Here we provide evidence from stochastic simulations which shows that the topology and organization of DNA packed inside the capsid plays a key role in determining these properties. Even with similar osmotic pressure pushing out the DNA, we find that spatially ordered DNA spools have a much lower effective friction than disordered entangled states. Such spools are only found when the tendency of nearby DNA strands to align locally is accounted for. This topological or conformational friction also depends on DNA knot type in the packing geometry and slows down or arrests the ejection of twist knots and very complex knots. We also find that the family of (2, 2k+1) torus knots unravel gradually by simplifying their topology in a stepwise fashion. Finally, an analysis of DNA trajectories inside the capsid shows that the knots formed throughout the ejection process mirror those found in gel electrophoresis experiments for viral DNA molecules extracted from the capsids.

Is the kinetoplast DNA a percolating network of linked rings at its critical point?

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.