Unexpected entanglement dynamics in semidilute blends of supercoiled and ring DNA

From Bioreactor to Bulk Rheology: Achieving Scalable Production of Highly Concentrated Circular DNA

DNA serves as a model system in polymer physics due to its ability to be obtained as a uniform polymer with controllable topology and nonequilibrium behavior. Currently, a major obstacle in the widespread adoption of DNA is obtaining it on a scale and cost basis that accommodates bulk rheology and high-throughput screening. To address this, recent advancements in bioreactor-based plasmid DNA production is coupled with anion exchange chromatography producing a unified approach to generating gram-scale quantities of monodisperse DNA. With this method, 1.1 grams of DNA is obtained per batch to generate solutions with concentrations up to 116 mg mL-1. This solution of uniform supercoiled and relaxed circular plasmid DNA, is roughly 69 times greater than the overlap concentration. The utility of this method is demonstrated by performing bulk rheology measurements at sample volumes up to 1 mL on DNA of different lengths, topologies, and concentrations. The measured elastic moduli are orders of magnitude larger than those previously reported for DNA and allowed for the construction of a time-concentration superposition curve that spans 12 decades of frequency. Ultimately, these results can provide important insights into the dynamics of ring polymers and the nature of highly condensed DNA dynamics.

Switchable microscale stress response of actin-vimentin composites emerges from scale-dependent interactions

The mechanical properties of the mammalian cell regulate many cellular functions and are largely dictated by the cytoskeleton, a composite network of protein filaments, including actin, microtubules, and intermediate filaments. Interactions between these distinct filaments give rise to emergent mechanical properties that are difficult to generate synthetically, and recent studies have made great strides in advancing our understanding of the mechanical interplay between actin and microtubule filaments. While intermediate filaments play critical roles in the stress response of cells, their effect on the rheological properties of the composite cytoskeleton remains poorly understood. Here, we use optical tweezers microrheology to measure the linear viscoelastic properties and nonlinear stress response of composites of actin and vimentin with varying molar ratios of actin to vimentin. We reveal a surprising, nearly opposite effect of actin-vimentin network mechanics compared to single-component networks in the linear versus nonlinear regimes. Namely, the linear elastic plateau modulus and zero-shear viscosity are markedly reduced in composites compared to single-component networks of actin or vimentin, whereas the initial response force and stiffness are maximized in composites versus single-component networks in the nonlinear regime. While these emergent trends are indicative of distinct interactions between actin and vimentin, nonlinear stiffening and longtime stress response appear to both be dictated primarily by actin, at odds with previous bulk rheology studies. We demonstrate that these complex, scale-dependent effects arise from the varied contributions of network density, filament stiffness, non-specific interactions, and poroelasticity to the mechanical response at different spatiotemporal scales. Cells may harness this complex behavior to facilitate distinct stress responses at different scales and in response to different stimuli to allow for their hallmark multifunctionality.

Supercoiled ring polymers under shear flow

We apply monomer-resolved computer simulations of supercoiled ring polymers under shear, taking full account of the hydrodynamic interactions, accompanied, in parallel, by simulations in which these are switched off. The combination of bending and torsional rigidities inherent in these polymers, in conjunction with hydrodynamics, has a profound impact on their flow properties. In contrast to their flexible counterparts, which dramatically deform and inflate under shear [Liebetreu et al., Commun. Mater. 2020, 1, 4], supercoiled rings undergo only weak changes in their overall shape and they display both a reduced propensity to tumbling (at fixed Weissenberg number) and a much stronger orientational resistance with respect to their flexible counterparts. In the presence of hydrodynamic interactions, the coupling of the polymer to solvent flow is capable of bringing about a topological transformation of writhe to twist at strong shear upon conservation of the overall linking number.

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.

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.

Optical-Tweezers-integrating-Differential-Dynamic-Microscopy maps the spatiotemporal propagation of nonlinear strains in polymer blends and composites

How local stresses propagate through polymeric fluids, and, more generally, how macromolecular dynamics give rise to viscoelasticity are open questions vital to wide-ranging scientific and industrial fields. Here, to unambiguously connect polymer dynamics to force response, and map the deformation fields that arise in macromolecular materials, we present Optical-Tweezers-integrating-Differential -Dynamic-Microscopy (OpTiDMM) that simultaneously imposes local strains, measures resistive forces, and analyzes the motion of the surrounding polymers. Our measurements with blends of ring and linear polymers (DNA) and their composites with stiff polymers (microtubules) uncover an unexpected resonant response, in which strain alignment, superdiffusivity, and elasticity are maximized when the strain rate is comparable to the entanglement rate. Microtubules suppress this resonance, while substantially increasing elastic storage, due to varying degrees to which the polymers buildup, stretch and flow along the strain path, and configurationally relax induced stress. More broadly, the rich multi-scale coupling of mechanics and dynamics afforded by OpTiDDM, empowers its interdisciplinary use to elucidate non-trivial phenomena that sculpt stress propagation dynamics-critical to commercial applications and cell mechanics alike.

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.

DNA Conformation Dictates Strength and Flocculation in DNA-Microtubule Composites

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

i-RheoFT: Fourier transforming sampled functions without artefacts

In this article we present a new open-access code named “i-RheoFT” that implements the analytical method first introduced in [PRE, 80, 012501 (2009)] and then enhanced in [New J Phys 14, 115032 (2012)], which allows to evaluate the Fourier transform of any generic time-dependent function that vanishes for negative times, sampled at a finite set of data points that extend over a finite range, and need not be equally spaced. I-RheoFT has been employed here to investigate three important experimental factors: (i) the ‘density of initial experimental points’ describing the sampled function, (ii) the interpolation function used to perform the “virtual oversampling” procedure introduced in [New J Phys 14, 115032 (2012)], and (iii) the detrimental effect of noises on the expected outcomes. We demonstrate that, at relatively high signal-to-noise ratios and density of initial experimental points, all three built-in MATLAB interpolation functions employed in this work (i.e., Spline, Makima and PCHIP) perform well in recovering the information embedded within the original sampled function; with the Spline function performing best. Whereas, by reducing either the number of initial data points or the signal-to-noise ratio, there exists a threshold below which all three functions perform poorly; with the worst performance given by the Spline function in both the cases and the least worst by the PCHIP function at low density of initial data points and by the Makima function at relatively low signal-to-noise ratios. We envisage that i-RheoFT will be of particular interest and use to all those studies where sampled or time-averaged functions, often defined by a discrete set of data points within a finite time-window, are exploited to gain new insights on the systems’ dynamics.

Chain-End Effects on Supramolecular Poly(ethylene glycol) Polymers

In this work we present a fundamental analysis based on small-angle scattering, linear rheology and differential scanning calorimetry (DSC) experiments of the role of different hydrogen bonding (H-bonding) types on the structure and dynamics of chain-end modified poly(ethylene glycol) (PEG) in bulk. As such bifunctional PEG with a molar mass below the entanglement mass Me is symmetrically end-functionalized with three different hydrogen bonding (H-bonding) groups: thymine-1-acetic acid (thy), diamino-triazine (dat) and 2-ureido-4[1H]-pyrimidinone (upy). A linear block copolymer structure and a Newtonian-like dynamics is observed for PEG-thy/dat while results for PEG-upy structure and dynamics reveal a sphere and a network-like behavior, respectively. These observations are concomitant with an increase of the Flory-Huggins interaction parameter from PEG-thy/dat to PEG-upy that is used to quantify the difference between the H-bonding types. The upy association into spherical clusters is established by the Percus-Yevick approximation that models the inter-particle structure factor for PEG-upy. Moreover, the viscosity study reveals for PEG-upy a shear thickening behavior interpreted in terms of the free path model and related to the time for PEG-upy to dissociate from the upy clusters, seen as virtual crosslinks of the formed network. Moreover, a second relaxation time of different nature is also obtained from the complex shear modulus measurements of PEG-upy by the inverse of the angular frequency where G' and G'' crosses from the network-like to glass-like transition relaxation time, which is related to the segmental friction of PEG-upy polymeric network strands. In fact, not only do PEG-thy/dat and PEG-upy have different viscoelastic properties, but the relaxation times found for PEG-upy are much slower than the ones for PEG-thy/dat. However, the activation energy related to the association dynamics is very similar for both PEG-thy/dat and PEG-upy. Concerning the segmental dynamics, the glass transition temperature obtained from both rheological and calorimetric analysis is similar and increases for PEG-upy while for PEG-thy/dat is almost independent of association behavior. Our results show how supramolecular PEG properties vary by modifying the H-bonding association type and changing the molecular Flory-Huggins interaction parameter, which can be further explored for possible applications.

Myosin-driven actin-microtubule networks exhibit self-organized contractile dynamics

The cytoskeleton is a dynamic network of proteins, including actin, microtubules, and their associated motor proteins, that enables essential cellular processes such as motility, division, and growth. While actomyosin networks are extensively studied, how interactions between actin and microtubules, ubiquitous in the cytoskeleton, influence actomyosin activity remains an open question. Here, we create a network of co-entangled actin and microtubules driven by myosin II. We combine dynamic differential microscopy, particle image velocimetry, and particle tracking to show that both actin and microtubules undergo ballistic contraction with unexpectedly indistinguishable characteristics. This contractility is distinct from faster disordered motion and rupturing that active actin networks exhibit. Our results suggest that microtubules enable self-organized myosin-driven contraction by providing flexural rigidity and enhanced connectivity to actin networks. Beyond the immediate relevance to cytoskeletal dynamics, our results shed light on the design of active materials that can be precisely tuned by the network composition.

Hydrodynamic interactions in topologically linked ring polymers

Despite decades of interdisciplinary research on topologically linked ring polymers, their dynamics remain largely unstudied. These systems represent a major scientific challenge as they are often subject to both topological and hydrodynamic interactions (HI), which render dynamical solutions either mathematically intractable or computationally prohibitive. Here we circumvent these limitations by preaveraging the HI of linked rings. We show that the symmetry of ring polymers leads to a hydrodynamic decoupling of ring dynamics. This decoupling is valid even for nonideal polymers and nonequilibrium conditions. Physically, our findings suggest that the effects of topology and HI are nearly independent and do not act cooperatively to influence polymer dynamics. We use this result to develop highly efficient Brownian dynamics algorithms that offer enormous performance improvements over conventional methods and apply these algorithms to simulate catenated ring polymers at equilibrium, confirming the independence of topological effects and HI. The methods developed here can be used to study and simulate large systems of linked rings with HI, including kinetoplast DNA, Olympic gels, and poly[n]catenanes.

Stress relaxation in symmetric ring-linear polymer blends at low ring fractions

We combine linear viscoelastic measurements and modelling in order to explore the dynamics of blends of the same-molecular-weight ring and linear polymers in the regime of the low volume fraction (0.3 or lower) of the ring component. The stress relaxation modulus is affected by the constraint release (CR) of both rings and linear components due to the motion of linear chains. We develop a CR-based model of ring-linear blends that predicts the stress relaxation function in the low fraction regime of ring component in excellent agreement with experiments. Rings trapped by their entanglements with linear chains can only relax by linear-chain-induced constraint release, resulting in much slower relaxation of rings than of linear chains. The relative viscosity η ( ϕ R * ) / η L of the blend with respect to the linear melt viscosity η _L at ring overlap volume fraction ϕ R * increases proportionally to the square root of ring molecular weight M w , R . Our experimental results clearly demonstrate that it is possible to enhance the viscosity and simultaneously the structural relaxation time of linear polymer melts by adding a small fraction of ring polymers. These results not only provide fundamental insights into the physics of the CR process but also suggest ways to fine-tune the flow properties of linear polymers by means of adding rings.

On the diffusivity of ring polymers

A free volume theory of linear polyethylene (PE) melts [Soft Matter, 2019, 15, 9300-9309] was used to describe the crossover in the size dependence of the center-of-mass diffusion coefficients of ring PE melts from Dcm ∼ N-1.03±0.15 to Dcm ∼ N-1.88±0.14 at Nc ≈ 100. Parameters in the theory were determined theoretically by either the polymer reference interaction site model (PRISM) theory or molecular dynamics (MD) simulation data. The former method is computationally less expensive, whereas the latter gives more accurate results in calculating the temperature dependence of Dcm as different intermolecular and intramolecular potentials were explicitly included in the equation of motion. Both approaches were able to describe the dynamics below and above the crossover in the size dependence of Dcm.