Understanding the dynamics of rings in the melt in terms of the annealed tree model

Unified understanding of the impact of semiflexibility, concentration, and molecular weight on macromolecular-scale ring diffusion

Conformationally fluctuating, globally compact macromolecules such as polymeric rings, single-chain nanoparticles, microgels, and many-arm stars display complex dynamic behaviors due to their rich topological structure and intermolecular organization. Synthetic rings are hybrid objects with conformations that display both ideal random walk and compact globular features, which can serve as models of genomic DNA. To date, emphasis has been placed on the effect of ring molecular weight on their unusual behaviors. Here, we combine simulations and a microscopic force-level theory to build a unified understanding for how key aspects of ring dynamics depend on different tunable molecular properties including backbone rigidity, monomer concentration, degree of traditional entanglement, and molecular weight. Our large-scale molecular dynamics simulations of ring melts with very different backbone stiffnesses reveal unanticipated behaviors which agree well with our generalized theory. This includes a universal master curve for center-of-mass diffusion constants as a function of molecular weight scaled by a chemistry and thermodynamic state-dependent critical molecular weight that generalizes the concept of an entanglement cross-over for linear chains. The key physics is how backbone rigidity and monomer concentration induced changes of the entanglement length, interring packing, degree of interpenetration, and liquid compressibility slow down space-time dynamic-force correlations on macromolecular scales. A power law decay of the center-of-mass diffusion constant with inverse molecular weight squared is the first consequence, followed by an ultraslow activated hopping transport regime. Our results set the stage to address slow dynamics and kinetic arrest in different families of compact synthetic and biological polymeric systems.

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

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

Topology-Based Detection and Tracking of Deadlocks Reveal Aging of Active Ring Melts

Connecting the viscoelastic behavior of stressed ring melts to the various forms of entanglement that can emerge in such systems is still an open challenge. Here, we consider active ring melts, where stress is generated internally, and introduce a topology-based method to detect and track consequential forms of ring entanglements, namely, deadlocks. We demonstrate that, as stress accumulates, more and more rings are co-opted in a growing web of deadlocks that entrap many other rings by threading, bringing the system to a standstill. The method ought to help the study of topological aging in more general polymer contexts.

Mixing Linear Polymers with Rings and Catenanes: Bulk and Interfacial Behavior

We derive and parameterize effective interaction potentials between a multitude of different types of ring polymers and linear chains, varying the bending rigidity and solvent quality for the former species. We further develop and apply a density functional treatment for mixtures of both disconnected (chain-ring) and connected (chain-polycatenane) mixtures of the same, drawing coexistence binodals and exploring the ensuing response functions as well as the interface and wetting behavior of the mixtures. We show that worsening of the solvent quality for the rings brings about a stronger propensity for macroscopic phase separation in the linear-polycatenane mixtures, which is predominantly of the demixing type between phases of similar overall particle density. We formulate a simple criterion based on the effective interactions, allowing us to determine whether any specific linear-ring mixture will undergo a demixing phase separation.

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.

Non-Rouse behavior of short ring polymers in melts by molecular dynamics simulations

Short ring polymers are expected to behave nearly Rouse-like due to the little effect of topological constraints of non-knot and non-concatenation. However, this notion is questioned because of several simulation and experiment findings in recent times, which requires a further more quantitative study. Therefore, we perform a deep investigation of statics and dynamics of flexible short ring polymers (N < 2Ne) in melts via molecular dynamics simulations by further taking linear analogues as well as all-crossing ring and linear polymers with switched off topological constraints for comparisons and demonstrate the noticeable deviations from the Rouse model in terms of local and global scales. Although the overall size is compact, the subchains are swollen, which is traced back to the deeper "segmental correlation hole" effect. The same scaling relationship of the non-Gaussian deviation of the static structure factor holds, but the deviation magnitude of rings is larger than that of linear analogues. By checking the non-Gaussian parameter and autocorrelation function of center-of-mass velocity, the physical origin of anomalous sub-diffusions of short rings is identified as unscreened viscoelastic hydrodynamic interactions and not correlation hole effects, like linear analogues.

Unexpected Slow Relaxation Dynamics in Pure Ring Polymers Arise from Intermolecular Interactions

Ring polymers have fascinated scientists for decades, but experimental progress has been challenging due to the presence of linear chain contaminants that fundamentally alter dynamics. In this work, we report the unexpected slow stress relaxation behavior of concentrated ring polymers that arises due to ring-ring interactions and ring packing structure. Topologically pure, high molecular weight ring polymers are prepared without linear chain contaminants using cyclic poly(phthalaldehyde) (cPPA), a metastable polymer chemistry that rapidly depolymerizes from free ends at ambient temperatures. Linear viscoelastic measurements of highly concentrated cPPA show slow, non-power-law stress relaxation dynamics despite the lack of linear chain contaminants. Experiments are complemented by molecular dynamics (MD) simulations of unprecedentedly high molecular weight rings, which clearly show non-power-law stress relaxation in good agreement with experiments. MD simulations reveal substantial ring-ring interpenetrations upon increasing ring molecular weight or local backbone stiffness, despite the global collapsed nature of single ring conformation. A recently proposed microscopic theory for unconcatenated rings provides a qualitative physical mechanism associated with the emergence of strong inter-ring caging which slows down center-of-mass diffusion and long wavelength intramolecular relaxation modes originating from ring-ring interpenetrations, governed by the onset variable N/ND, where the crossover degree of polymerization ND is qualitatively predicted by theory. Our work overcomes challenges in achieving ring polymer purity and by characterizing dynamics for high molecular weight ring polymers. Overall, these results provide a new understanding of ring polymer physics.

To thread or not to thread? Effective potentials and threading interactions between asymmetric ring polymers

We use computer simulations to study a system of two unlinked ring polymers, whose length and bending stiffness are systematically varied. We derive the effective potentials between the rings, calculate the areas of minimal surfaces of the same, and characterize the threading between them. When the two rings are of the same kind, threading of a one ring through the surface of the other is immanent for small ring-ring separations. Flexible rings pierce the surface of the other ring several times but only shallowly, as compared to the stiff rings which pierce less frequently but deeply. Typically, the ring that is being threaded swells and flattens up into an oblate-like conformation, while the ring that is threading the other takes a shape of an elongated prolate. The roles of the threader and the threaded ring are being dynamically exchanged. If, on the other hand, the rings are of different kinds, the symmetry is broken and the rings tend to take up roles of the threader and the threaded ring with unequal probabilities. We propose a method how to predict these probabilities based on the parameters of the individual rings. Ultimately, our work captures the interactions between ring polymers in a coarse-grained fashion, opening the way to large-scale modelling of materials such as kinetoplasts, catenanes or topological brushes.

Entanglement Length Scale Separates Threading from Branching of Unknotted and Non-concatenated Ring Polymers in Melts

Current theories on the conformation and dynamics of unknotted and non-concatenated ring polymers in melt conditions describe each ring as a tree-like double-folded object. While evidence from simulations supports this picture on a single ring level, other works show pairs of rings also thread each other, a feature overlooked in the tree theories. Here we reconcile this dichotomy using Monte Carlo simulations of the ring melts with different bending rigidities. We find that rings are double-folded (more strongly for stiffer rings) on and above the entanglement length scale, while the threadings are localized on smaller scales. The different theories disagree on the details of the tree structure, i.e., the fractal dimension of the backbone of the tree. In the stiffer melts we find an indication of a self-avoiding scaling of the backbone, while more flexible chains do not exhibit such a regime. Moreover, the theories commonly neglect threadings and assign different importance to the impact of the progressive constraint release (tube dilation) on single ring relaxation due to the motion of other rings. Despite that each threading creates only a small opening in the double-folded structure, the threading loops can be numerous and their length can exceed substantially the entanglement scale. We link the threading constraints to the divergence of the relaxation time of a ring, if the tube dilation is hindered by pinning a fraction of other rings in space. Current theories do not predict such divergence and predict faster than measured diffusion of rings, pointing at the relevance of the threading constraints in unpinned systems as well. Revision of the theories with explicit threading constraints might elucidate the validity of the conjectured existence of topological glass.

Random Knotting in Fractal Ring Polymers

Many ring polymer systems of physical and biological interest exhibit both pronounced topological effects and nontrivial self-similarity, but the relationship between these two phenomena has not yet been clearly established. Here, we use theory and simulation to formulate such a connection by studying a fundamental topological property—the random knotting probability—for ring polymers with varying fractal dimension, d_f. Using straightforward scaling arguments, we generalize a classic mathematical result, showing that the probability of a trivial knot decays exponentially with chain size, N, for all fractal dimensions: P₀(N) ∝ exp(−N/N₀). However, no such simple considerations can account for the dependence of the knotting length, N₀, on d_f, necessitating a more involved analytical calculation. This analysis reveals a complicated double-exponential dependence, which is well supported by numerical data. By contrast, functional forms typical of simple scaling theories fail to adequately describe the observations. These findings are equally valid for two-dimensional ring polymer systems, where “knotting” is defined as the intersection of any two segments.

Branched droplet clusters and the Kramers theorem

Scaling laws inherent for polymer molecules are checked for the linear and branched chains constituting two-dimensional (2D) levitating microdroplet clusters condensed above the locally heated layer of water. We demonstrate that the dimensionless averaged end-to-end distance of the droplet chain r[over ¯] normalized by the averaged distance between centers of the adjacent droplets l[over ¯] scales as r[over ¯]/l[over ¯]∼n^{0.76}, where n is the number of links in the chain, which is close to the power exponent ¾, predicted for 2D polymer chains with excluded volume in the dilution limit. The values of the dimensionless Kuhn length b[over ̃]≅2.12±0.015 and of the averaged absolute value of the bond angle of the droplet chains |θ|[over ¯]=22.0±0.5^{0} are determined. Using these values we demonstrate that the predictions of the Kramers theorem for the gyration radius of branched polymers are valid also for the branched droplets' chains. We discuss physical interactions that explain both the high value of the power exponent and the applicability of the Kramers theorem including the effects of the excluded volume, surrounding droplet monomers, and the prohibition of extreme values of the bond angle.

Formulation of Chromatin Mobility as a Function of Nuclear Size during C. elegans Embryogenesis Using Polymer Physics Theories

During early embryogenesis of the nematode, Caenorhabditis elegans, the chromatin motion markedly decreases. Despite its biological implications, the underlying mechanism for this transition was unclear. By combining theory and experiment, we analyze the mean-square displacement (MSD) of the chromatin loci, and demonstrate that MSD-vs-time relationships in various nuclei collapse into a single master curve by normalizing them with the mesh size and the corresponding time scale. This enables us to identify the onset of the entangled dynamics with the size of tube diameter of chromatin polymer in the C. elegans embryo. Our dynamical scaling analysis predicts the transition between unentangled and entangled dynamics of chromatin polymers, the quantitative formula for MSD as a function of nuclear size and timescale, and provides testable hypotheses on chromatin mobility in other cell types and species.

Theory of Transient Localization, Activated Dynamics, and a Macromolecular Glass Transition in Ring Polymer Liquids

We construct a segmental scale force level theory for the center-of-mass diffusion constant and corresponding relaxation time for globally compact unconcatenated ring polymer solutions and melts (degree of polymerization N). The approach is based on slowly decaying macromolecular scale intermolecular force dynamic correlations as the origin of their unusual dynamics. Unentangled Rouse, weakly caged, and activated regimes are predicted. The barrier of the activated regime scales linearly with N and as a power law of concentration, which drives a kinetic glass transition on the radius-of-gyration scale. The values of N at the two dynamic crossovers (Rouse to weakly caged, weakly caged to activated) are proportional, with nonuniversality entering mainly via macromolecular volume fraction and dimensionless compressibility. Quantitative comparisons with simulation data reveal good agreement. Aspects of intermediate time dynamics are analyzed, and predictions are made for the conditions required to observe a macromolecular glass transition in the laboratory and on the computer.

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.

Topological Constraints with Optimal Length Promote the Formation of Chromosomal Territories at Weakened Degree of Phase Separation

It is generally agreed that the nuclei of eukaryotic cells at interphase are partitioned into disjointed territories, with distinct regions occupied by certain chromosomes. However, the underlying mechanism for such territorialization is still under debate. Here we model chromosomes as coarse-grained block copolymers and to investigate the effect of loop domains (LDs) on the formation of compartments and territories based on dissipative particle dynamics. A critical length of LDs, which depends sensitively on the length of polymeric blocks, is obtained to minimize the degree of phase separation. This also applies to the two-polymer system: The critical length not only maximizes the degree of territorialization but also minimizes the degree of phase separation. Interestingly, by comparing with experimental data, we find the critical length for LDs and the corresponding length of blocks to be respectively very close to the mean length of topologically associating domains (TADs) and chromosomal segments with different densities of CpG islands for human chromosomes. The results indicate that topological constraints with optimal length can contribute to the formation of territories by weakening the degree of phase separation, which likely promotes the chromosomal flexibility in response to genetic regulations.

Multiscale Approaches for Confined Ring Polymer Solutions

We apply a hierarchy of multiscale modeling approaches to investigate the structure of ring polymer solutions under planar confinement. In particular, we employ both monomer-resolved (MR-DFT) and a coarse-grained (CG-DFT) density functional theories for fully flexible ring polymers, with the former based on a flexible tangent hard-sphere model and the latter based on an effective soft-colloid representation, to elucidate the ring polymer organization within slits of variable width in different concentration regimes. The predicted monomer and polymer center-of-mass densities in confinement, as well as the surface tension at the solution-wall interface, are compared to explicit molecular dynamics (MD) simulations. The approaches yield quantitative (MR-DFT) or semiquantitative (CG-DFT) agreement with MD. In addition, we provide a systematic comparison between confined linear and ring polymer solutions. When compared to their linear counterparts, the rings are found to feature a higher propensity to structure in confinement that translates into a distinct shape of the depletion potentials between two walls immersed into a polymer solution. The depletion potentials that we extract from CG-DFT and MR-DFT are in semiquantitative agreement with each other. Overall, we find consistency among all approaches as regards the shapes, trends, and qualitative characteristics of density profiles and depletion potentials induced on hard walls by linear and cyclic polymers.

Self-Similar Dynamics of Large Polymer Rings: A Neutron Spin Echo Study

This work clarifies the self-similar dynamics of large polymer rings using pulsed-field gradient nuclear magnetic resonance and neutron spin echo spectroscopy. We find center of mass diffusion taking place in three dynamic regimes starting (i) with a strongly subdiffusive domain ⟨r^{2}(t)⟩_{com}∼t^{α} (0.4≤α≤0.65); (ii) a second subdiffusive region ⟨r^{2}(t)⟩_{com}∼t^{0.75} that (iii) finally crosses over to Fickian diffusion. While the t^{0.75} range previously has been found in simulations and was predicted by theory, we attribute the first to the effect of cooperative dynamics resulting from the correlation hole potential. The internal dynamics at scales below the elementary loop size is well described by ring Rouse motion. At larger scales the dynamics is self-similar and follows very well the predictions of the scaling models with preference for the self-consistent fractal loopy globule model.

Effects of Topological Constraints on Penetration Structures of Semi-Flexible Ring Polymers

The effects of topological constraints on penetration structures of semi-flexible ring polymers in a melt are investigated using molecular dynamics simulations, considering simultaneously the effects of the chain stiffness. Three topology types of rings are considered: 01-knot (the unknotted), 31-knot and 61-knot ring polymers, respectively. With the improved algorithm to detect and quantify the inter-ring penetration (or inter-ring threading), the degree of ring threading does not increase monotonously with the chain stiffness, existing a peak value at the intermediate stiffness. It indicates that rings interpenetrate most at intermediate stiffness where there is a balance between coil expansion (favoring penetrations) and stiffness (inhibiting penetrations). Meanwhile, the inter-ring penetration would be suppressed with the knot complexity of the rings. The analysis of effective potential between the rings provides a better understanding for this non-monotonous behavior in inter-ring penetration.

Statistical Properties of Lasso-Shape Polymers and Their Implications for Complex Lasso Proteins Function

The shape and properties of closed loops depend on various topological factors. One of them is loop-threading, which is present in complex lasso proteins. In this work, we analyze the probability of loop-threading by the tail and its influence on the shape of the loop measured by the radius of gyration, distention, asphericity, and prolateness. In particular, we show that the probability of a trivial lasso for phantom polymer is non-zero even for an infinite structure, as well as that the threading flattens the loop by restricting its motion in one dimension. These results are further used to show that there are fewer non-trivial protein lassos than expected and select potentially functional complex lasso proteins.

Slow Dynamics of Ring Polymer Melts by Asymmetric Interaction of Threading Configuration: Monte Carlo Study of a Dynamically Constrained Lattice Model

Abnormally slower diffusional processes than its internal structure relaxation have been observed in ring polymeric melt systems recently. A key structural feature in ring polymer melts is topological constraints which allow rings to assume a threading configuration in the melt phase. In this work, we constructed a lattice model under the assumption of asymmetric diffusivity between two threading rings, and investigated a link between the structural correlation and its dynamic behavior via Monte Carlo simulations. We discovered that the hierarchical threading configurations render the whole system to exhibit abnormally slow dynamics. By analyzing statistical distributions of timescales of threading configurations, we found that the decoupling between internal structure relaxation and diffusion is crucial to understand the threading effects on the dynamics of a ring melt. In particular, in the limit of small but threaded rings, scaling exponents of the diffusion coefficient D and timescale τ diff with respect to the degree of polymerization N agree well with that of the annealed tree model as well as our mean-field analysis. As N increases, however, the ring diffusion abruptly slows down to the glassy behavior, which is supported by a breakdown of the Stokes⁻Einstein relation.

Threading of Unconcatenated Ring Polymers at High Concentrations: Double-Folded vs Time-Equilibrated Structures

Unconcatenated ring polymers in concentrated solutions and melt are remarkably well described as double-folded conformations on randomly branched primitive trees. This picture though contrasts recent evidence for extensive intermingling between close-by rings in the form of long-lived topological constraints or threadings. Here, we employ the concept of ring minimal surface to quantify the extent of threadings in polymer solutions of the double-folded rings vs rings in equilibrated molecular dynamics computer simulations. Our results show that the double-folded ring polymers are significantly less threaded compared to their counterparts at equilibrium. Second, threadings form through a slow process whose characteristic time-scale is of the same order of magnitude as that of the diffusion of the rings in solution. These findings are robust, being based on universal (model-independent) observables as the average fraction of threaded length or the total penetrations between close-by rings and the corresponding distribution functions.

Glassiness and Heterogeneous Dynamics in Dense Solutions of Ring Polymers

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.

Ring Polymers: Threadings, Knot Electrophoresis and Topological Glasses

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.

Nanoparticle Motion in Entangled Melts of Linear and Nonconcatenated Ring Polymers

The motion of nanoparticles (NPs) in entangled melts of linear polymers and nonconcatenated ring polymers are compared by large-scale molecular dynamics simulations. The comparison provides a paradigm for the effects of polymer architecture on the dynamical coupling between NPs and polymers in nanocomposites. Strongly suppressed motion of NPs with diameter d larger than the entanglement spacing a is observed in a melt of linear polymers before the onset of Fickian NP diffusion. This strong suppression of NP motion occurs progressively as d exceeds a and is related to the hopping diffusion of NPs in the entanglement network. In contrast to the NP motion in linear polymers, the motion of NPs with d > a in ring polymers is not as strongly suppressed prior to Fickian diffusion. The diffusion coefficient D decreases with increasing d much slower in entangled rings than in entangled linear chains. NP motion in entangled nonconcatenated ring polymers is understood through a scaling analysis of the coupling between NP motion and the self-similar entangled dynamics of ring polymers.

Beyond Flory theory: Distribution functions for interacting lattice trees

While Flory theories [J. Isaacson and T. C. Lubensky, J. Physique Lett. 41, 469 (1980)JPSLBO0302-072X10.1051/jphyslet:019800041019046900; M. Daoud and J. F. Joanny, J. Physique 42, 1359 (1981)JOPQAG0302-073810.1051/jphys:0198100420100135900; A. M. Gutin et al., Macromolecules 26, 1293 (1993)MAMOBX0024-929710.1021/ma00058a016] provide an extremely useful framework for understanding the behavior of interacting, randomly branching polymers, the approach is inherently limited. Here we use a combination of scaling arguments and computer simulations to go beyond a Gaussian description. We analyze distribution functions for a wide variety of quantities characterizing the tree connectivities and conformations for the four different statistical ensembles, which we have studied numerically in [A. Rosa and R. Everaers, J. Phys. A: Math. Theor. 49, 345001 (2016)1751-811310.1088/1751-8113/49/34/345001 and J. Chem. Phys. 145, 164906 (2016)JCPSA60021-960610.1063/1.4965827]: (a) ideal randomly branching polymers, (b) 2d and 3d melts of interacting randomly branching polymers, (c) 3d self-avoiding trees with annealed connectivity, and (d) 3d self-avoiding trees with quenched ideal connectivity. In particular, we investigate the distributions (i) p_{N}(n) of the weight, n, of branches cut from trees of mass N by severing randomly chosen bonds; (ii) p_{N}(l) of the contour distances, l, between monomers; (iii) p_{N}(r[over ⃗]) of spatial distances, r[over ⃗], between monomers, and (iv) p_{N}(r[over ⃗]|l) of the end-to-end distance of paths of length l. Data for different tree sizes superimpose, when expressed as functions of suitably rescaled observables x[over ⃗]=r[over ⃗]/sqrt[〈r^{2}(N)〉] or x=l/〈l(N)〉. In particular, we observe a generalized Kramers relation for the branch weight distributions (i) and find that all the other distributions (ii-iv) are of Redner-des Cloizeaux type, q(x[over ⃗])=C|x|^{θ}exp(-(K|x|)^{t}). We propose a coherent framework, including generalized Fisher-Pincus relations, relating most of the RdC exponents to each other and to the contact and Flory exponents for interacting trees.

LassoProt: server to analyze biopolymers with lassos

The LassoProt server, http://lassoprot.cent.uw.edu.pl/, enables analysis of biopolymers with entangled configurations called lassos. The server offers various ways of visualizing lasso configurations, as well as their time trajectories, with all the results and plots downloadable. Broad spectrum of applications makes LassoProt a useful tool for biologists, biophysicists, chemists, polymer physicists and mathematicians. The server and our methods have been validated on the whole PDB, and the results constitute the database of proteins with complex lassos, supported with basic biological data. This database can serve as a source of information about protein geometry and entanglement-function correlations, as a reference set in protein modeling, and for many other purposes.

Minimal Surfaces on Unconcatenated Polymer Rings in Melt

In order to quantify the effect of mutual threading on conformations and dynamics of unconcatenated and unknotted rings in the melt we computationally examine their minimal surfaces. We found a linear scaling of the surface area with the ring length. Minimal surfaces allow for an unambiguous algorithmic definition of mutual threading between rings. Based on it, we found that, although ring threading is frequent, majority of cases correspond to short loops. These findings explain why approximate theories that neglect threading are so unexpectedly successful despite having no small parameter justification. We also examine threading dynamics and identify the threading order parameter that reflects the ring diffusivity.

Sequence Dependence of Viral RNA Encapsidation

We develop a Flory mean-field theory for viral RNA (vRNA) molecules that extends the current RNA folding algorithms to include interactions between different sections of the secondary structure. The theory is applied to sequence-selective vRNA encapsidation. The dependence on sequence enters through a single parameter: the largest eigenvalue of the Kramers matrix of the branched polymer obtained by coarse graining the secondary structure. Differences between the work of encapsidation of vRNA molecules and of randomized isomers are found to be in the range of 20 kBT, more than sufficient to provide a strong bias in favor of vRNA encapsidation. The method is applied to a packaging competition experiment where large vRNA molecules compete for encapsidation with two smaller RNA species that together have the same nucleotide sequence as the large molecule. We encounter a substantial, generic free energy bias, that also is of the order of 20 kBT, in favor of encapsidating the single large RNA molecule. The bias is mainly the consequence of the fact that dividing up a large vRNA molecule involves the release of stored elastic energy. This provides an important, nonspecific mechanism for preferential encapsidation of single larger vRNA molecules over multiple smaller mRNA molecules with the same total number of nucleotides. The result is also consistent with recent RNA packaging competition experiments by Comas-Garcia et al.1 Finally, the Flory method leads to the result that when two RNA molecules are copackaged, they are expected to remain segregated inside the capsid.

A topologically driven glass in ring polymers

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.

Torque-winding interdependence for a flexible polymer chain wound around a cylinder in the presence of obstacles

Polymer chains winding around each other or around other objects occur in many natural systems; the physical consequences of this winding are therefore of significant interest. A polymer chain could be surrounded by various bulky objects (referred as obstacles), such as other macromolecules or macromolecular aggregates. Here we show that for a long flexible polymer chain wound around a cylinder, the presence of obstacles could modify the winding-torque interdependence, in some cases leading to phase-transition-like behavior in which the winding occurs only when the torque exceeds some critical value. Possible implications of this effect are discussed in relation to the biophysics of nucleic acids.

Self-Similar Conformations and Dynamics in Entangled Melts and Solutions of Nonconcatenated Ring Polymers

A scaling model of self-similar conformations and dynamics of nonconcatenated entangled ring polymers is developed. Topological constraints force these ring polymers into compact conformations with fractal dimension d_f = 3 that we call fractal loopy globules (FLGs). This result is based on the conjecture that the overlap parameter of subsections of rings on all length scales is the same and equal to the Kavassalis-Noolandi number O_KN ≈ 10-20. The dynamics of entangled rings is self-similar and proceeds as loops of increasing sizes are rearranged progressively at their respective diffusion times. The topological constraints associated with smaller rearranged loops affect the dynamics of larger loops through increasing the effective friction coefficient but have no influence on the entanglement tubes confining larger loops. As a result, the tube diameter defined as the average spacing between relevant topological constraints increases with time t, leading to "tube dilation". Analysis of the primitive paths in molecular dynamics simulations suggests a complete tube dilation with the tube diameter on the order of the time-dependent characteristic loop size. A characteristic loop at time t is defined as a ring section that has diffused a distance equal to its size during time t. We derive dynamic scaling exponents in terms of fractal dimensions of an entangled ring and the underlying primitive path and a parameter characterizing the extent of tube dilation. The results reproduce the predictions of different dynamic models of a single nonconcatenated entangled ring. We demonstrate that traditional generalization of single-ring models to multi-ring dynamics is not self-consistent and develop a FLG model with self-consistent multi-ring dynamics and complete tube dilation. This selfconsistent FLG model predicts that the longest relaxation time of nonconcatenated entangled ring polymers scales with their degree of polymerization N as τ_relax ~ N^7/3, while the diffusion coefficient of these rings scales as D_3d ~ N^-5/3. For the entangled solutions and melts of rings, we predict power law stress relaxation function G(t) ~ t^-3/7 at t < τ_relax without a rubbery plateau and the corresponding viscosity scaling with the degree of polymerization N as η ~ N^4/3. These theoretical predictions are in good agreement with recent computer simulations and are consistent with experiments of melts of nonconcatenated entangled rings.

Anomalous diffusion in fractal globules

The fractal globule state is a popular model for describing chromatin packing in eukaryotic nuclei. Here we provide a scaling theory and dissipative particle dynamics computer simulation for the thermal motion of monomers in the fractal globule state. Simulations starting from different entanglement-free initial states show good convergence which provides evidence supporting the existence of a unique metastable fractal globule state. We show monomer motion in this state to be subdiffusive described by ⟨X(2)(t)⟩∼t(αF) with αF close to 0.4. This result is in good agreement with existing experimental data on the chromatin dynamics, which makes an additional argument in support of the fractal globule model of chromatin packing.