Dynamical Entanglement and Cooperative Dynamics in Entangled Solutions of Ring and Linear Polymers

Enzymatic cleaving of entangled DNA rings drives scale-dependent rheological trajectories

DNA, which naturally occurs in linear, ring, and supercoiled topologies, frequently undergoes enzyme-driven topological conversion and fragmentation in vivo, enabling it to perform a variety of functions within the cell. In vitro, highly concentrated DNA polymers form entanglements that yield viscoelastic properties dependent on the topologies and lengths of the DNA. Enzyme-driven alterations of DNA size and shape therefore offer a means of designing active materials with programmable viscoelastic properties. Here, we incorporate multi-site restriction endonucleases into dense DNA solutions to linearize and fragment circular DNA molecules. We pair optical tweezers microrheology with differential dynamic microscopy and single-molecule tracking to measure the linear and nonlinear viscoelastic response and transport properties of entangled DNA solutions over a wide range of spatiotemporal scales throughout the course of enzymatic digestion. We show that, at short timescales, relative to the relaxation timescales of the polymers, digestion of these 'topologically-active' fluids initially causes an increase in elasticity and relaxation times followed by a gradual decrease. Conversely, for long timescales, linear viscoelastic moduli exhibit signatures of increasing elasticity. DNA diffusion, likewise, becomes increasingly slowed, in direct opposition to the short-time behavior. We hypothesize that this scale-dependent rheology arises from the population of small DNA fragments, which increases as digestion proceeds, driving self-association of larger fragments via depletion interactions, giving rise to slow relaxation modes of clusters of entangled chains, interspersed among shorter unentangled fragments. While these slow modes likely dominate at long times, they are presumably frozen out in the short-time limit, which instead probes the faster relaxation modes of the unentangled population.

Unraveling the Glass-like Dynamic Heterogeneity in Ring Polymer Melts: From Semiflexible to Stiff Chain

Ring polymers are an intriguing class of polymers with unique physical properties, and understanding their behavior is important for developing accurate theoretical models. In this study, we investigate the effect of chain stiffness and monomer density on the static and dynamic behaviors of ring polymer melts using molecular dynamics simulations. Our first focus is on the non-Gaussian parameter of center-of-mass displacement as a measure of dynamic heterogeneity, which is commonly observed in glass-forming liquids. We find that the non-Gaussianity in the displacement distribution increases with the monomer density and stiffness of the polymer chains, suggesting that excluded volume interactions between centers of mass have a strong effect on the dynamics of ring polymers. We then analyze the relationship between the radius of gyration and monomer density for semiflexible and stiff ring polymers. Our results indicate that the relationship between the two varies with chain stiffness, which can be attributed to the competition between repulsive forces inside the ring and from adjacent rings. Finally, we study the dynamics of bond-breakage virtually connected between the centers of mass of rings to analyze the exchanges of intermolecular networks of bonds. Our results demonstrate that the dynamic heterogeneity of bond-breakage is coupled with the non-Gaussianity in ring polymer melts, highlighting the importance of the bond-breaking method in determining the intermolecular dynamics of ring polymer melts. Overall, our study sheds light on the factors that govern the dynamic behaviors of ring polymers.

Cooperative Rheological State-Switching of Enzymatically-Driven Composites of Circular DNA And Dextran

Polymer topology, which plays a principal role in the rheology of polymeric fluids, and non-equilibrium materials, which exhibit time-varying rheological properties, are topics of intense investigation. Here, composites of circular DNA and dextran are pushed out-of-equilibrium via enzymatic digestion of DNA rings to linear fragments. These time-resolved rheology measurements reveal discrete state-switching, with composites undergoing abrupt transitions between dissipative and elastic-like states. The gating time and lifetime of the elastic-like states, and the magnitude and sharpness of the transitions, are surprisingly decorrelated from digestion rates and non-monotonically depend on the DNA fraction. These results are modeled using sigmoidal two-state functions to show that bulk state-switching can arise from continuous molecular-level activity due to the necessity for cooperative percolation of entanglements to support macroscopic stresses. This platform, coupling the tunability of topological composites with the power of enzymatic reactions, may be leveraged for diverse material applications from wound-healing to self-repairing infrastructure.

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.

Entropic Mixing of Ring/Linear Polymer Blends

The topological constraints of nonconcatenated ring polymers force them to form compact loopy globular conformations with much lower entropy than unconstrained ideal rings. The closed-loop structure of ring polymers also enables them to be threaded by linear polymers in ring/linear blends, resulting in less compact ring conformations with higher entropy. This conformational entropy increase promotes mixing rings with linear polymers. Here, using molecular dynamics simulations for bead-spring chains, ring/linear blends are shown to be significantly more miscible than linear/linear blends and that there is an entropic mixing, negative χ, for ring/linear blends compared to linear/linear and ring/ring blends. In analogy with small angle neutron scattering, the static structure function S(q) is measured, and the resulting data are fit to the random phase approximation model to determine χ. In the limit that the two components are the same, χ = 0 for the linear/linear and ring/ring blends as expected, while χ < 0 for the ring/linear blends. With increasing chain stiffness, χ for the ring/linear blends becomes more negative, varying reciprocally with the number of monomers between entanglements. Ring/linear blends are also shown to be more miscible than either ring/ring or linear/linear blends and stay in single phase for a wider range of increasing repulsion between the two components.

Topological digestion drives time-varying rheology of entangled DNA fluids

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

Understanding and controlling the rheology of polymeric complex fluids is of fundamental importance in both industry and biology. Here, Michieletto et al. show how to achieve time-dependent rheology of DNA solutions via enzymatically-driven architectural alterations by restriction endonucleases.

Tuning Protein Hydrogel Mechanics through Modulation of Nanoscale Unfolding and Entanglement in Postgelation Relaxation

Globular folded proteins are versatile nanoscale building blocks to create biomaterials with mechanical robustness and inherent biological functionality due to their specific and well-defined folded structures. Modulating the nanoscale unfolding of protein building blocks during network formation (in situ protein unfolding) provides potent opportunities to control the protein network structure and mechanics. Here, we control protein unfolding during the formation of hydrogels constructed from chemically cross-linked maltose binding protein using ligand binding and the addition of cosolutes to modulate protein kinetic and thermodynamic stability. Bulk shear rheology characterizes the storage moduli of the bound and unbound protein hydrogels and reveals a correlation between network rigidity, characterized as an increase in the storage modulus, and protein thermodynamic stability. Furthermore, analysis of the network relaxation behavior identifies a crossover from an unfolding dominated regime to an entanglement dominated regime. Control of in situ protein unfolding and entanglement provides an important route to finely tune the architecture, mechanics, and dynamic relaxation of protein hydrogels. Such predictive control will be advantageous for future smart biomaterials for applications which require responsive and dynamic modulation of mechanical properties and biological function.

Geometric Predictors of Knotted and Linked Arcs

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.

Decoding polymer self-dynamics using a two-step approach

The self-correlation function and corresponding self-intermediate scattering function in Fourier space are important quantities for describing the molecular motions of liquids. This work draws attention to a largely overlooked issue concerning the analysis of these space-time density-density correlation functions of polymers. We show that the interpretation of non-Gaussian behavior of polymers is generally complicated by intrachain averaging of distinct self-dynamics of different segments. By the very nature of the mathematics involved, the averaging process not only conceals critical dynamical information, but also contributes to the observed non-Gaussian dynamics. To fully expose this issue and provide a thorough benchmark of polymer self-dynamics, we perform analyses of coarse-grained molecular dynamics simulations of linear and ring polymer melts as well as several theoretical models using a "two-step" approach, where interchain and intrachain averagings of segmental self-dynamics are separated. While past investigations primarily focused on the average behavior, our results indicate that a more nuanced approach to polymer self-dynamics is clearly required.

Superstretchable Elastomer from Cross-linked Ring Polymers

The stretchability of polymeric materials is critical to many applications such as flexible electronics and soft robotics, yet the stretchability of conventional cross-linked linear polymers is limited by the entanglements between polymer chains. We show using molecular dynamics simulations that cross-linked ring polymers are significantly more stretchable than cross-linked linear polymers. Compared to linear polymers, the entanglements between ring polymers do not act as effective cross-links. As a result, the stretchability of cross-linked ring polymers is determined by the maximum extension of polymer strands between cross-links, rather than between trapped entanglements as in cross-linked linear polymers. The more compact conformation of ring polymers before deformation also contributes to the increase in stretchability.

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.

Effects of chain length on Rouse modes and non-Gaussianity in linear and ring polymer melts

The dynamics of ring polymer melts are studied via molecular dynamics simulations of the Kremer-Grest bead-spring model. Rouse mode analysis is performed in comparison with linear polymers by changing the chain length. Rouse-like behavior is observed in ring polymers by quantifying the chain length dependence of the Rouse relaxation time, whereas a crossover from Rouse to reptation behavior is observed in linear polymers. Furthermore, the non-Gaussian parameters of the monomer bead displacement and chain center-of-mass displacement are analyzed. It is found that the non-Gaussianity of ring polymers is remarkably suppressed with slight growth for the center-of-mass dynamics at long chain length, which is in contrast to the growth in linear polymers for both the monomer bead and center-of-mass dynamics.

Topological tuning of DNA mobility in entangled solutions of supercoiled plasmids

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.