Effect of Bending Rigidity on the Knotting of a Polymer under Tension

Compression of a confined semiflexible polymer under direct and oscillating fields

The folding transition of biopolymers from the coil to compact structures has attracted wide research interest in the past and is well studied in polymer physics. Recent seminal works on DNA in confined devices have shown that these long biopolymers tend to collapse under an external field, which is contrary to the previously reported stretching of the chain. In this work, we capture the compression of a confined semiflexible polymer under direct and oscillating fields using a coarse-grained computer simulation model in the presence of long-range hydrodynamics. In the case of a semiflexible polymer chain, the inhomogeneous hydrodynamic drag from the center to the periphery of the coil couples with the chain bending to cause a swirling movement of the chain segments, leading to structural intertwining and compaction. Contrarily, a flexible chain of the same length lacks such structural deformation and forms a well-established tadpole structure. While bending rigidity profoundly influences the chain's folding favorability, we also found that subject to the direct field, chains in stronger confinements exhibit substantial compaction, contrary to the one in moderate confinements or bulk where such compaction is absent. However, an alternating field within an optimum frequency can effectuate this compression even in moderate or no confinement. This field-induced collapse is a quintessential hydrodynamic phenomenon, resulting in intertwined knotted structures even for shorter chains, unlike other spontaneous knotting experiments where it happens exclusively for longer chains.

Phase transition of DNA knotting in spherical space

Knots have been discovered in various biological systems, such as DNA. The knotting probability of DNA in free space depends non-monotonically on its bending rigidity and has a prominent peak. The current work aims to understand the underlying mechanism of the non-monotonic dependence of DNA knotting probability on bending rigidity. Monte Carlo simulations are performed on a closed DNA molecule confined in spherical space described by a worm-like chain model and a flexible kink model, respectively. The closed DNA's contour length and the spherical space radius both increase knotting probability, but also alter the unimodal dependence of knotting probability on bending rigidity. This is generalized using universal phase diagrams based on the two models. Under the flexible kink model, the total knotting probability of closed DNA is obviously increased at a relatively high excited energy. This supports the expectation that the entropy effect of knot size favours knot formation at a relatively low bending rigidity. In a given spherical space, the increasing contour length of closed DNA described by the worm-like chain model results in a visible shift in the knotting probability distribution. At the same time, the gyration radius of non-trivial closed DNA becomes comparable to that of trivial closed DNA, so that their ratio is not anti-correlated with average knot length. For closed DNA of various contour lengths, the relationship between average knot length and bending rigidity has a universal behaviour: the average knot length decreases to a local minimum at a bending rigidity of ∼5 and then gradually increases to a constant value. The existence of the local minimum is determined by the cut-off distance in repulsive Lennard-Jones potential. The bending rigidity corresponding to the beginning of the constant average knot length is consistent with that at the peak in the knotting distribution. At this point, the knot-size effect balances with the fragment free-energy effect and, at an even greater bending rigidity, knot length breathes around the average knot length value.

Free energy of a knotted polymer confined to narrow cylindrical and conical channels

Monte Carlo simulations are used to study the conformational behavior of a semiflexible polymer confined to cylindrical and conical channels. The channels are sufficiently narrow that the conditions for the Odijk regime are marginally satisfied. For cylindrical confinement, we examine polymers with a single knot of topology 3_{1}, 4_{1}, or 5_{1}, as well as unknotted polymers that are capable of forming S loops. We measure the variation of the free energy F with the end-to-end polymer extension length X and examine the effect of varying the polymer topology, persistence length P, and cylinder diameter D on the free-energy functions. Similarly, we characterize the behavior of the knot span along the channel. We find that increasing the knot complexity increases the typical size of the knot. In the regime of low X, where the knot/S-loop size is large, the conformational behavior is independent of polymer topology. In addition, the scaling properties of the free energy and knot span are in agreement with predictions from a theoretical model constructed using known properties of interacting polymers in the Odijk regime. We also examine the variation of F with the position of a knot in conical channels for various values of the cone angle α. The free energy decreases as the knot moves in a direction where the cone widens, and it also decreases with increasing α and with increasing knot complexity. The behavior is in agreement with predictions from a theoretical model in which the dominant contribution to the change in F is the change in the size of the hairpins as the knot moves to the wider region of the channel.

Non-equilibrium effects of molecular motors on polymers

We present a generic coarse-grained model to describe molecular motors acting on polymer substrates, mimicking, for example, RNA polymerase on DNA or kinesin on microtubules. The polymer is modeled as a connected chain of beads; motors are represented as freely diffusing beads which, upon encountering the substrate, bind to it through a short-ranged attractive potential. When bound, motors and polymer beads experience an equal and opposite active force, directed tangential to the polymer; this leads to motion of the motors along the polymer contour. The inclusion of explicit motors differentiates our model from other recent active polymer models. We study, by means of Langevin dynamics simulations, the effect of the motor activity on both the conformational and dynamical properties of the substrate. We find that activity leads, in addition to the expected enhancement of polymer diffusion, to an effective reduction of its persistence length. We discover that this effective "softening" is a consequence of the emergence of double-folded branches, or hairpins, and that it can be tuned by changing the number of motors or the force they generate. Finally, we investigate the effect of the motors on the probability of knot formation. Counter-intuitively our simulations reveal that, even though at equilibrium a more flexible substrate would show an increased knotting probability, motor activity leads to a marked decrease in the occurrence of knotted conformations with respect to equilibrium.

Ring polymer dynamics and tumbling-stretch transitions in planar mixed flows

The properties of dilute polymer solutions are governed by the conformational dynamics of individual polymers which can be perturbed in the presence of an applied flow. Much of our understanding of dilute solutions comes from studying how flows manipulate the molecular features of polymer chains out of equilibrium, primarily focusing on linear polymer chains. Recently there has been an emerging interest in the dynamics of nonlinear architectures, particularly ring polymers, which exhibit surprising out-of-equilibrium dynamics in dilute solutions. In particular, it has been observed that hydrodynamics can couple to topology in planar elongational and shear flows, driving molecular expansion in the nonflow direction that is not observed for linear chains. In this paper, we extend our understanding of dilute ring polymer dynamics to mixed flows, which represent flow profiles intermediate between simple shear or planar elongation. We map the conformational behaviors at a number of flow geometries and strengths, demonstrating transitions between coiled, tumbling, and stretched regimes. Indeed, these observations are consistent with how linear chains respond to mixed flows. For both linear and ring polymers, we observe a marked first-order-like transition between tumbling and stretched polymers that we attribute to a dynamic energy barrier between the two states. This manifests as bimodal extension distributions in a narrow range of flow strengths and geometries, with the primary difference between rings and linear chains being the presence of molecular expansion in the vorticity direction.

Trefoil Knot Hydrodynamic Delocalization on Sheared Ring Polymers

The behavior of unknotted and trefoil-knotted ring polymers under shear flow is here examined by means of mesoscopic simulations. In contrast to most polymers, ring polymers in a hydrodynamic solvent at high shear rates do not get shortened in the vorticity direction. This is a consequence of the backflow produced by the interaction of the sheared solvent with the end-free polymer topology. The extended structures of the ring in the vorticity-flow plane, when they are aligned in a constant velocity plane, favor ring contour fluctuations. This variety of conformations largely suppresses the tank-treading type of rotation with extended conformations in favor of the tumbling type of rotations, where stretched and collapsed conformations alternate. The extension of trefoil knots is also enhanced, so that the knots become delocalized. We anticipate that these effects, which disappear in the absence of hydrodynamic interactions, will have a crucial impact on the rheological properties of concentrated ring solutions, and will also influence the behavior of more complicated systems such as mixtures of polymers with different topologies.

Molecular Knots

The first synthetic molecular trefoil knot was prepared in the late 1980s. However, it is only in the last few years that more complex small-molecule knot topologies have been realized through chemical synthesis. The steric restrictions imposed on molecular strands by knotting can impart significant physical and chemical properties, including chirality, strong and selective ion binding, and catalytic activity. As the number and complexity of accessible molecular knot topologies increases, it will become increasingly useful for chemists to adopt the knot terminology employed by other disciplines. Here we give an overview of synthetic strategies towards molecular knots and outline the principles of knot, braid, and tangle theory appropriate to chemistry and molecular structure.

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

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

Trapping a Knot into Tight Conformations by Intra-Chain Repulsions

Knots can occur in biopolymers such as DNA and peptides. In our previous study, we systematically investigated the effects of intra-chain interactions on knots and found that long-range repulsions can surprisingly tighten knots. Here, we use this knowledge to trap a knot into tight conformations in Langevin dynamics simulations. By trapping, we mean that the free energy landscape with respect to the knot size exhibits a potential well around a small knot size in the presence of long-range repulsions, and this potential can well lead to long-lived tight knots when its depth is comparable to or larger than thermal energy. We tune the strength of intra-chain repulsion such that a knot is weakly trapped. Driven by thermal fluctuations, the knot can escape from the trap and is then re-trapped. We find that the knot switches between tight and loose conformations-referred to as "knot breathing". We use a Yukawa potential to model screened electrostatic interactions to explore the relevance of knot trapping and breathing in charged biopolymers. We determine the minimal screened length and the minimal strength of repulsion for knot trapping. We find that Coulomb-induced knot trapping is possible to occur in single-stranded DNA and peptides for normal ionic strengths.

Tuning knot abundance in semiflexible chains with crowders of different sizes: a Monte Carlo study of DNA chains

We use stochastic simulation techniques to sample the conformational space of linear semiflexible polymers in a crowded medium and study how the knotting properties depend on the crowder size and concentration. The abundance of physical knots in the chains, which for definiteness we model on 10 kb long DNA filaments, is shown to have a non-monotonic, unimodal dependence on the colloid diameter, dc. The maximum incidence of knots occurs when dc is about equal to half of the gyration radius of the isolated chain. The degree of enhancement of knots grows rapidly with the solution density and can be very conspicuous relative to the case of isolated chains with no crowders. For instance, at 30% volume fraction the relative increase is more than fourfold. This dramatic enhancement is shown to originate from the depletion-induced chain compaction over multiple and concurring length scales. The same effect accounts for the variations of the knot length that accompany the changes in knotting probability. The findings suggest that crowded media could be viably used as a passive physical means for controlling and modulating the incidence and length of knots in DNA and other types of semiflexible polymers.

Stretching Response of Knotted and Unknotted Polymer Chains

Recent theoretical and experimental advances have clarified the major effects of knotting on the properties of stretched chains. Yet, how knotted chains respond to weak mechanical stretching and how this behavior differs from the unknotted case are still open questions and we address them here by profiling the complete stretching response of chains of hundreds of monomers and different topology. We find that the ratio of the knotted and unknotted chain extensions varies nonmonotonically with the applied force. This surprising feature is shown to be a signature of the crossover between the well-known high-force stretching regime and the previously uncharacterized low-force one. The observed differences of knotted and unknotted chain response increases with knot complexity and are sufficiently marked that they could be harnessed in single-molecule contexts to infer the presence and complexity of physical knots in micron-long biomolecules.

An Anisotropic Effective Model for the Simulation of Semiflexible Ring Polymers

We derive and introduce anisotropic effective pair potentials to coarse-grain solutions of semiflexible ring polymers of various lengths. The system has been recently investigated by means of full monomer-resolved computer simulations, revealing a host of unusual features and structure formation, which, however, cannot be captured by a rotationally averaged effective pair potential between the rings' centers of mass [Bernabei M.; Soft Matter2013, 9, 1287]. Our new coarse-graining strategy is to picture each ring as a soft, penetrable disk. We demonstrate that for the short- and intermediate-length rings the new model is quite capable of capturing the physics in a quantitative fashion, whereas for the largest rings, which resemble flexible ones, it fails at high densities. Our work opens the way for the physical justification of general, anisotropic penetrable interaction potentials.

Origin of metastable knots in single flexible chains

Recent theoretical progress has explained the physics of knotting of semiflexible polymers, yet knotting of flexible polymers is relatively unexplored. We herein develop a new theory for the size distribution of knots on a flexible polymer and the existence of metastable knots. We show the free energy of a flexible molecule in a tube can be mapped to quantitatively reproduce the free energy distribution of a knot on a flexible chain. The size distribution of knots on flexible chains is expected to be universal and might be observed at a macroscopic scale, such as a string of hard balls.

Multi-blob coarse graining for ring polymer solutions

We present a multi-scale molecular modeling of concentrated solutions of unknotted and non-concatenated ring polymers under good solvent conditions. The approach is based on a multi-blob representation of each ring polymer, which is capable of overcoming the shortcomings of single-blob approaches that lose their validity at concentrations exceeding the overlap density of the solution [A. Narros, A. J. Moreno, and C. N. Likos, Soft Matter, 2010, 6, 2435]. By means of a first principles coarse-graining strategy based on analytically determined effective pair potentials between the blobs, computed at zero density, we quantitatively reproduce the single molecule and solution properties of a system with well-defined topological constraints. Detailed comparisons with the underlying, monomer-resolved model demonstrate the validity of our approach, which employs fully transferable pair potentials between connected and unconnected blobs. We demonstrate that the pair structure between the centers of mass of the rings is accurately reproduced by the multi-blob approach, thus opening the way for simulation of arbitrarily long polymers. Finally, we show the importance of the topological constraint of non-concatenation on the structure of the concentrated solution and in particular on the size of the correlation hole and the shrinkage of the rings as melt concentrations are approached.

Untying Knotted DNA with Elongational Flows

We present Brownian dynamics simulations of initially knotted double-stranded DNA molecules untying in elongational flows. We show that the motions of the knots are governed by a diffusion-convection equation by deriving scalings that collapse the simulation data. When being convected, all knots displace nonaffinely, and their rates of translation along the chain are topologically dictated. We discover that torus knots "corkscrew" when driven by flow, whereas nontorus knots do not. We show that a simple mechanism can explain a coupling between this rotation and the translation of a knot, explaining observed differences in knot translation rates. These types of knots are encountered in nanoscale manipulation of DNA, occur in biology at multiple length scales (DNA to umbilical cords), and are ubiquitous in daily life (e.g., hair). These results may have a broad impact on manipulations of such knots via flows, with applications to genomic sequencing and polymer processing.

Knotting and Unknotting Dynamics of DNA Strands in Nanochannels

The self-knotting dynamics of DNA strands confined in nanochannels is studied with Brownian simulations. The model DNA chains are several microns long and placed inside channels that are 50-300 nm wide. This width range covers the transition between different metric scaling regimes and the concomitant drop of DNA knotting probability for channel widths below ∼75 nm. We find that knots typically originate from deep looping and backfoldings of the chain ends. Upon lowering the channel width, backfoldings become shallower and rarer and the lifetime of knots decreases while that of unknots increases. This lifetimes interplay causes the dramatic reduction of knots incidence for increasing confinement. The results can aid the design of nanochannels capable of harnessing the self-knotting dynamics to quench or relax the DNA topological state as desired.

Driving knots on DNA with AC/DC electric fields: topological friction and memory effects

The dynamical properties of entangled polyelectrolytes are investigated theoretically and computationally for a proposed novel micromanipulation setup. Specifically, we investigate the effects of DC and AC electric fields acting longitudinally on knotted DNA chains, modelled as semiflexible chains of charged beads, under mechanical tension. We consider various experimentally accessible values of the field amplitude and frequency as well as several of the simplest knot types. In particular, we consider both torus and twist knots because they are respectively known to be able or unable to slide along macroscopic threads and ropes. Strikingly, this qualitative distinction disappears in this microscopic context because all the considered knot types acquire a systematic drift in the direction of the electric force. Notably, the knot drift velocity and diffusion coefficient in zero field (both measurable also experimentally) can be used to define a characteristic "frictional" lengthscale for the various knot types. This previously unexplored length provides valuable information on the extent of self-interactions in the nominal knotted region. It is finally observed that the motion of a knot can effectively follow the AC field only if the driving period is larger than the knot relaxation time (for which the self-diffusion time provides an upper bound). These results suggest that salient aspects of the intrinsic dynamics of knots in DNA chains could be probed experimentally by means of external, time-dependent electric fields.

Influence of Rigidity and Knot Complexity on the Knotting of Confined Polymers

We employ computer simulations and thermodynamic integration to analyze the effects of bending rigidity and slit confinement on the free energy cost of tying knots, ΔF_knotting, on polymer chains under tension. A tension-dependent, nonzero optimal stiffness κ_min exists, for which ΔF_knotting is minimal. For a polymer chain with several stiffness domains, each containing a large amount of monomers, the domain with stiffness κ_min will be preferred by the knot. A local analysis of the bending in the interior of the knot reveals that local stretching of chains at the braid region is responsible for the fact that the tension-dependent optimal stiffness has a nonzero value. The reduction in ΔF_knotting for a chain with optimal stiffness relative to the flexible chain can be enhanced by tuning the slit width of the 2D confinement and increasing the knot complexity. The optimal stiffness itself is independent of the knot types we considered, while confinement shifts it toward lower values.