Beyond Flory theory: Distribution functions for interacting lattice trees

Amoeba Monte Carlo algorithms for random trees with controlled branching activity: Efficient trial move generation and universal dynamics

The reptation Monte Carlo algorithm is a simple, physically motivated and efficient method for equilibrating semidilute solutions of linear polymers. Here, we propose two simple generalizations for the analog Amoeba algorithm for randomly branching chains, which allow us to efficiently deal with random trees with controlled branching activity. We analyze the rich relaxation dynamics of Amoeba algorithms and demonstrate the existence of an unexpected scaling regime for the tree relaxation. Our results suggest that the equilibration time for Amoeba algorithms scales in general like N^{2}〈n_{lin}〉^{Δ}, where N denotes the number of tree nodes, 〈n_{lin}〉 the mean number of linear segments the trees are composed of, and Δ≃0.4.

DNA supercoiling in bacteria: state of play and challenges from a viewpoint of physics based modeling

DNA supercoiling is central to many fundamental processes of living organisms. Its average level along the chromosome and over time reflects the dynamic equilibrium of opposite activities of topoisomerases, which are required to relax mechanical stresses that are inevitably produced during DNA replication and gene transcription. Supercoiling affects all scales of the spatio-temporal organization of bacterial DNA, from the base pair to the large scale chromosome conformation. Highlighted in vitro and in vivo in the 1960s and 1970s, respectively, the first physical models were proposed concomitantly in order to predict the deformation properties of the double helix. About fifteen years later, polymer physics models demonstrated on larger scales the plectonemic nature and the tree-like organization of supercoiled DNA. Since then, many works have tried to establish a better understanding of the multiple structuring and physiological properties of bacterial DNA in thermodynamic equilibrium and far from equilibrium. The purpose of this essay is to address upcoming challenges by thoroughly exploring the relevance, predictive capacity, and limitations of current physical models, with a specific focus on structural properties beyond the scale of the double helix. We discuss more particularly the problem of DNA conformations, the interplay between DNA supercoiling with gene transcription and DNA replication, its role on nucleoid formation and, finally, the problem of scaling up models. Our primary objective is to foster increased collaboration between physicists and biologists. To achieve this, we have reduced the respective jargon to a minimum and we provide some explanatory background material for the two communities.

Crumpled polymer with loops recapitulates key features of chromosome organization

Chromosomes are exceedingly long topologically-constrained polymers compacted in a cell nucleus. We recently suggested that chromosomes are organized into loops by an active process of loop extrusion. Yet loops remain elusive to direct observations in living cells; detection and characterization of myriads of such loops is a major challenge. The lack of a tractable physical model of a polymer folded into loops limits our ability to interpret experimental data and detect loops. Here, we introduce a new physical model - a polymer folded into a sequence of loops, and solve it analytically. Our model and a simple geometrical argument show how loops affect statistics of contacts in a polymer across different scales, explaining universally observed shapes of the contact probability. Moreover, we reveal that folding into loops reduces the density of topological entanglements, a novel phenomenon we refer as "the dilution of entanglements". Supported by simulations this finding suggests that up to ~ 1 - 2Mb chromosomes with loops are not topologically constrained, yet become crumpled at larger scales. Our theoretical framework allows inference of loop characteristics, draws a new picture of chromosome organization, and shows how folding into loops affects topological properties of crumpled polymers.

The Physical Behavior of Interphase Chromosomes: Polymer Theory and Coarse-Grain Computer Simulations

Fluorescence in situ hybridization and chromosome conformation capture methods point to the same conclusion: that chromosomes appear to the external observer as compact structures with a highly nonrandom three-dimensional organization. In this work, we recapitulate the efforts made by us and other groups to rationalize this behavior in terms of the mathematical language and tools of polymer physics. After a brief introduction dedicated to some crucial experiments dissecting the structure of interphase chromosomes, we discuss at a nonspecialistic level some fundamental aspects of theoretical and numerical polymer physics. Then, we inglobe biological and polymer aspects into a polymer model for interphase chromosomes which moves from the observation that mutual topological constraints, such as those typically present between polymer chains in ordinary melts, induce slow chain dynamics and "constraint" chromosomes to resemble double-folded randomly branched polymer conformations. By explicitly turning these ideas into a multi-scale numerical algorithm which is described here in full details, we can design accurate model polymer conformations for interphase chromosomes and offer them for systematic comparison to experiments. The review is concluded by discussing the limitations of our approach and pointing to promising perspectives for future work.

Monte Carlo simulation of a lattice model for the dynamics of randomly branching double-folded ring polymers

Supercoiled DNA, crumpled interphase chromosomes, and topologically constrained ring polymers often adopt treelike, double-folded, randomly branching configurations. Here we study an elastic lattice model for tightly double-folded ring polymers, which allows for the spontaneous creation and deletion of side branches coupled to a diffusive mass transport, which is local both in space and on the connectivity graph of the tree. We use Monte Carlo simulations to study systems falling into three different universality classes: ideal double-folded rings without excluded volume interactions, self-avoiding double-folded rings, and double-folded rings in the melt state. The observed static properties are in good agreement with exact results, simulations, and predictions of Flory theory for randomly branching polymers. For example, in the melt state rings adopt compact configurations and exhibit territorial behavior. In particular, we show that the emergent dynamics is in excellent agreement with a recent scaling theory and illustrate the qualitative differences with the familiar reptation dynamics of linear chains.

Influence of Branching on the Configurational and Dynamical Properties of Entangled Polymer Melts

We probe the influence of branching on the configurational, packing, and density correlation function properties of polymer melts of linear and star polymers, with emphasis on molecular masses larger than the entanglement molecular mass of linear chains. In particular, we calculate the conformational properties of these polymers, such as the hydrodynamic radius R h , packing length p, pair correlation function g ( r ) , and polymer center of mass self-diffusion coefficient, D, with the use of coarse-grained molecular dynamics simulations. Our simulation results reproduce the phenomenology of simulated linear and branched polymers, and we attempt to understand our observations based on a combination of hydrodynamic and thermodynamic modeling. We introduce a model of "entanglement" phenomenon in high molecular mass polymers that assumes polymers can viewed in a coarse-grained sense as "soft" particles and, correspondingly, we model the emergence of heterogeneous dynamics in polymeric glass-forming liquids to occur in a fashion similar to glass-forming liquids in which the molecules have soft repulsive interactions. Based on this novel perspective of polymer melt dynamics, we propose a functional form for D that can describe our simulation results for both star and linear polymers, covering both the unentangled to entangled polymer melt regimes.

Influence of knot complexity on glass-formation in low molecular mass ring polymer melts

We perform molecular dynamics simulations on a coarse-grained polymer melt to study the dynamics of glass-formation in ring polymer melts of variable knot complexity. After generating melts of non-concatenated polymeric rings having a range of minimum crossing number values, mc, we compute the coherent intermediate scattering function, the segmental α-relaxation time, fragility, and the glass transition temperature as a function of mc. Variation of knot complexity is found to have a pronounced effect on the dynamics of polymer melts since both molecular rigidity and packing are altered, primary physical factors governing glass-formation in polymeric materials.

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

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

Conformational statistics of randomly branching double-folded ring polymers

The conformations of topologically constrained double-folded ring polymers can be described as wrappings of randomly branched primitive trees. We extend previous work on the tree statistics under different (solvent) conditions to explore the conformational statistics of double-folded rings in the limit of tight wrapping. In particular, we relate the exponents characterizing the ring statistics to those describing the primitive trees and discuss the distribution functions [Formula: see text] and [Formula: see text] for the spatial distance, [Formula: see text], and tree contour distance, L, between monomers as a function of their ring contour distance, [Formula: see text].

The Ising model in swollen vs. compact polymers: Mean-field approach and computer simulations

We study the properties of the classical Ising model with nearest-neighbor interaction for spins located at the monomers of long polymer chains in 2 and 3 dimensions. We compare results for two ensembles of polymers with very different single chain properties: 1) swollen, self-avoiding linear polymer chains in good solvent conditions and 2) compact, space-filling randomly branching polymers in melt. By employing a mean-field approach and Monte Carlo computer simulations, we show that swollen polymers cannot sustain an ordered phase. On the contrary, compact polymers may indeed produce an observable phase transition. Finally, we briefly consider the statistical properties of the ordered phase by comparing polymer chains within the same universality class but characterized by very different shapes.

A comparative study of thermodynamic, conformational, and structural properties of bottlebrush with star and ring polymer melts

Thermodynamic, conformational, and structural properties of bottlebrush polymer melts are investigated with molecular dynamics simulations and compared to linear, regular star, and unknotted ring polymer melts to gauge the influence of molecular topology on polymer melt properties. We focus on the variation of the backbone chain length, the grafting density along the backbone, and the length of the side chains at different temperatures above the melt glass transition temperature. Based on these comparisons, we find that the segmental density, isothermal compressibility, and isobaric thermal expansion of bottlebrush melts are quantitatively similar to unknotted ring polymer melts and star polymer melts having a moderate number ( f = 5 to 6) of arms. These similarities extend to the mass scaling of the chain radius of gyration. Our results together indicate that the configurational properties of bottlebrush polymers in their melt state are more similar to randomly branched polymers than linear polymer chains. We also find that the average shape of bottlebrush polymers having short backbone chains with respect to the side chain length is also rather similar to the unknotted ring and moderately branched star polymers in their melt state. As a general trend, the molecular shape of bottlebrush polymers becomes more spherically symmetric when the length of the side chains has a commensurate length as the backbone chain. Finally, we calculate the partial static structure factor of the backbone segments and we find the emergence of a peak at the length scales that characterizes the average distance between the backbone chains. This peak is absent when we calculate the full static structure factor. We characterize the scaling of this peak with parameters characterizing the bottlebrush molecular architecture to aid in the experimental characterization of these molecules by neutron scattering.

Marginally compact hyperbranched polymer trees

Assuming Gaussian chain statistics along the chain contour, we generate by means of a proper fractal generator hyperbranched polymer trees which are marginally compact. Static and dynamical properties, such as the radial intrachain pair density distribution ρ_pair(r) or the shear-stress relaxation modulus G(t), are investigated theoretically and by means of computer simulations. We emphasize that albeit the self-contact density diverges logarithmically with the total mass N, this effect becomes rapidly irrelevant with increasing spacer length S. In addition to this it is seen that the standard Rouse analysis must necessarily become inappropriate for compact objects for which the relaxation time τ_p of mode p must scale as τ_p ∼ (N/p)^5/3 rather than the usual square power law for linear chains.

Flory theory of randomly branched polymers

Randomly branched polymer chains (or trees) are a classical subject of polymer physics with connections to the theory of magnetic systems, percolation and critical phenomena. More recently, the model has been reconsidered for RNA, supercoiled DNA and the crumpling of topologically-constrained polymers. While solvable in the ideal case, little is known exactly about randomly branched polymers with volume interactions. Flory theory provides a simple, unifying description for a wide range of branched systems, including isolated trees in good and θ-solvent, and tree melts. In particular, the approach provides a common framework for the description of randomly branched polymers with quenched connectivity and for randomly branching polymers with annealed connectivity. Here we review the Flory theory for interacting trees in the asymptotic limit of very large polymerization degree for good solvent, θ-solutions and melts, and report its predictions for annealed connectivity in θ-solvents. We compare the predictions of Flory theory for randomly branched polymers to a wide range of available analytical and numerical results and conclude that they are qualitatively excellent and quantitatively good in most cases.