Numerical Study of Linear and Circular Model DNA Chains Confined in a Slit: Metric and Topological Properties

Unraveling the Influence of Topology and Spatial Confinement on Equilibrium and Relaxation Properties of Interlocked Ring Polymers

We use Langevin dynamics simulations to study linked ring polymers in channel confinement. We address the in- and out-of-equilibrium behavior of the systems for varying degrees of confinement and increasing topological and geometrical complexity of the interlocking. The main findings are three. First, metric observables of different link topologies collapse onto the same master curve when plotted against the crossing number, revealing a universal response to confinement. Second, the relaxation process from initially stretched states is faster for more complex links. We ascribe these properties to the interplay of several effects, including the dependence of topological friction on the link complexity. Finally, we show that transient forms of geometrical entanglement purposely added to the initial stressed state can leave distinctive signatures in force-spectroscopy curves. The insight provided by the findings could be leveraged in single-molecule nanochannel experiments to identify geometric entanglement within topologically linked rings.

Spatial Organization of Slit-Confined Melts of Ring Polymers with Nonconserved Topology: A Lattice Monte Carlo Study

We present Monte Carlo computer simulations for melts of semiflexible randomly knotted and randomly concatenated ring polymers on the fcc lattice and in slit confinement. Through systematic variation of the slit width at fixed melt density, we explore the influence of confinement on single-chain conformations and interchain interactions. We demonstrate that confinement makes chains globally larger and more elongated while enhancing both contacts and knottedness propensities. As for multichain properties, we show that ring-ring contacts decrease with the confinement, yet neighboring rings overlap more as confinement grows. These aspects are accompanied by a marked decrease in the links formed between pairs of neighboring rings. In connection with the quantitative relation between links and entanglements in polymer melts recently established by us [Ubertini M. A.; Rosa A.Macromolecules2023, 56, 3354-3362], we propose that confinement can be used to set polymer networks that act softer under mechanical stress and suggest a viable experimental setup to validate our results.

Multi-feature clustering of CTCF binding creates robustness for loop extrusion blocking and Topologically Associating Domain boundaries

Topologically Associating Domains (TADs) separate vertebrate genomes into insulated regulatory neighborhoods that focus genome-associated processes. TADs are formed by Cohesin-mediated loop extrusion, with many TAD boundaries consisting of clustered binding sites of the CTCF insulator protein. Here we determine how this clustering of CTCF binding contributes to the blocking of loop extrusion and the insulation between TADs. We identify enrichment of three features of CTCF binding at strong TAD boundaries, consisting of strongly bound and closely spaced CTCF binding peaks, with a further enrichment of DNA-binding motifs within these peaks. Using multi-contact Nano-C analysis in cells with normal and perturbed CTCF binding, we establish that individual CTCF binding sites contribute to the blocking of loop extrusion, but in an incomplete manner. When clustered, individual CTCF binding sites thus create a stepwise insulation between neighboring TADs. Based on these results, we propose a model whereby multiple instances of temporal loop extrusion blocking create strong insulation between TADs.

Diffusion of knots in nanochannel-confined DNA molecules

We used Langevin dynamics simulations without hydrodynamic interactions to probe knot diffusion mechanisms and the time scales governing the evolution and the spontaneous untying of trefoil knots in nanochannel-confined DNA molecules in the extended de Gennes regime. The knot untying follows an "opening up process," wherein the initially tight knot continues growing and fluctuating in size as it moves toward the end of the DNA molecule before its annihilation at the chain end. The mean knot size increases significantly and sub-linearly with increasing chain contour length. The knot diffusion in nanochannel-confined DNA molecules is subdiffusive, with the unknotting time scaling with chain contour length with an exponent of 2.64 ± 0.23 to within a 95% confidence interval. The scaling exponent for the mean unknotting time vs chain contour length, along with visual inspection of the knot conformations, suggests that the knot diffusion mechanism is a combination of self-reptation and knot region breathing for the simulated parameters.

Influence of ionic conditions on knotting in a coarse-grained model for DNA

We investigate knotting probabilities of long double-stranded DNA strands in a coarse-grained Kratky-Porod model for DNA with Monte Carlo simulations. Various ionic conditions are implemented by adjusting the effective diameter of monomers. We find that the occurrence of knots in DNA can be reinforced considerably by high salt conditions and confinement between plates. Likewise, knots can almost be dissolved completely in a low salt scenario. Comparisons with recent experiments confirm that the coarse-grained model is able to capture and quantitatively predict topological features of DNA and can be used for guiding future experiments on DNA knots.

Quantifying the effects of slit confinement on polymer knots using the tube model

Knots can spontaneously form in DNA, proteins, and other polymers and affect their properties. These knots often experience spatial confinement in biological systems and experiments. While confinement dramatically affects the knot behavior, the physical mechanisms underlying the confinement effects are not fully understood. In this work, we provide a simple physical picture of the polymer knots in slit confinement using the tube model. In the tube model, the polymer segments in the knot core are assumed to be confined in a virtual tube due to the topological restriction. We first perform Monte Carlo simulation of a flexible knotted chain confined in a slit. We find that with the decrease of the slit height from H=+∞ (the 3D case) to H=2a (the 2D case), the most probable knot size L_{knot}^{*} dramatically shrinks from (L_{knot}^{*})_{3D}≈140a to (L_{knot}^{*})_{2D}≈26a, where a is the monomer diameter of the flexible chain. Then we quantitatively explain the confinement-induced knot shrinking and knot deformation using the tube model. Our results for H=2a can be applied to a polymer knot on a surface, which resembles DNA knots measured by atomic force microscopy under the conditions that DNA molecules are weakly absorbed on the surface and reach equilibrium 2D conformations. This work demonstrates the effectiveness of the tube model in understanding polymer knots.

Topological and physical links in soft matter systems

Linking, or multicomponent topological entanglement, is ubiquitous in soft matter systems, from mixtures of polymers and DNA filaments packedin vivoto interlocked line defects in liquid crystals and intertwined synthetic molecules. Yet, it is only relatively recently that theoretical and experimental advancements have made it possible to probe such entanglements and elucidate their impact on the physical properties of the systems. Here, we review the state-of-the-art of this rapidly expanding subject and organize it as follows. First, we present the main concepts and notions, from topological linking to physical linking and then consider the salient manifestations of molecular linking, from synthetic to biological ones. We next cover the main physical models addressing mutual entanglements in mixtures of polymers, both linear and circular. Finally, we consider liquid crystals, fluids and other non-filamentous systems where topological or physical entanglements are observed in defect or flux lines. We conclude with a perspective on open challenges.

Investigation of Ring and Star Polymers in Confined Geometries: Theory and Simulations

The calculations of the dimensionless layer monomer density profiles for a dilute solution of phantom ideal ring polymer chains and star polymers with f=4 arms in a Θ-solvent confined in a slit geometry of two parallel walls with repulsive surfaces and for the mixed case of one repulsive and the other inert surface were performed. Furthermore, taking into account the Derjaguin approximation, the dimensionless layer monomer density profiles for phantom ideal ring polymer chains and star polymers immersed in a solution of big colloidal particles with different adsorbing or repelling properties with respect to polymers were calculated. The density-force relation for the above-mentioned cases was analyzed, and the universal amplitude ratio B was obtained. Taking into account the small sphere expansion allowed obtaining the monomer density profiles for a dilute solution of phantom ideal ring polymers immersed in a solution of small spherical particles, or nano-particles of finite size, which are much smaller than the polymer size and the other characteristic mesoscopic length of the system. We performed molecular dynamics simulations of a dilute solution of linear, ring, and star-shaped polymers with N=300, 300 (360), and 1201 (4 × 300 + 1-star polymer with four arms) beads accordingly. The obtained analytical and numerical results for phantom ring and star polymers are compared with the results for linear polymer chains in confined geometries.

Effects of topological constraints on linked ring polymers in solvents of varying quality

We investigate the effects of topological constraints in catenanes composed of interlinked ring polymers on their size in a good solvent as well as on the location of their θ-point when the solvent quality is worsened. We mainly focus on poly[n]catenanes consisting of n ring polymers each of length m interlocked in a linear fashion. Using molecular dynamics simulations, we study the scaling of the poly[n]catenane's radius of gyration in a good solvent, assuming in general that Rg∼mμnν and we find that μ = 0.65 ± 0.02 and ν = 0.60 ± 0.01 for the range of n and m considered. These findings are further rationalized with the help of a mean-field Flory-like theory yielding the values of μ = 16/25 and ν = 3/5, consistent with the numerical results. We show that individual rings within catenanes feature a surplus swelling due to the presence of NL topological links. Furthermore, we consider poly[n]catenanes in solvents of varying quality and we demonstrate that the presence of topological links leads to an increase of its θ-temperature in comparison to isolated linear and ring chains with the following ordering: T > T > T. Finally, we show that the presence of links similarly raises the θ-temperature of a single linked ring in comparison to an unlinked one, bringing its θ-temperature close to the one of a poly[n]catenane.

Correlation anisotropy and stiffness of DNA molecules confined in nanochannels

The anisotropy of orientational correlations in DNA molecules confined in cylindrical channels is explored by Monte Carlo simulations using a coarse-grained model of double-stranded (ds) DNA. We find that the correlation function ⟨C(s)⟩_⊥ in the transverse (confined) dimension exhibits a region of negative values in the whole range of channel sizes. Such a clear-cut sign of the opposite orientation of chain segments represents a microscopic validation of the Odijk deflection mechanism in narrow channels. At moderate-to-weak confinement, the negative ⟨C(s)⟩_⊥ correlations imply a preference of DNA segments for transverse looping. The inclination for looping can explain a reduction of stiffness as well as the enhanced knotting of confined DNA relative to that detected earlier in bulk at some channel sizes. Furthermore, it is shown that the orientational persistence length P_or fails to convey the apparent stiffness of DNA molecules in channels. Instead, correlation lengths P_∥ and P_⊥ in the axial and transverse directions, respectively, encompass the channel-induced modifications of DNA stiffness.

KymoKnot: A web server and software package to identify and locate knots in trajectories of linear or circular polymers

The KymoKnot software package and web server identifies and locates physical knots or proper knots in a series of polymer conformations. It is mainly intended as an analysis tool for trajectories of linear or circular polymers, but it can be used on single instances too, e.g. protein structures in PDB format. A key element of the software package is the so-called minimally interfering chain closure algorithm that is used to detect physical knots in open chains and to locate the knotted region in both open and closed chains. The web server offers a user-friendly graphical interface that identifies the knot type and highlights the knotted region on each frame of the trajectory, which the user can visualize interactively from various viewpoints. The dynamical evolution of the knotted region along the chain contour is presented as a kymograph. All data can be downloaded in text format. The KymoKnot package is licensed under the BSD 3-Clause licence. The server is publicly available at http://kymoknot.sissa.it/kymoknot/interactive.php .

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.

Evidence for the extended de Gennes regime of a semiflexible polymer in slit confinement

We use off-lattice, pruned-enriched Rosenbluth method (PERM) simulations to compute the confinement free energy of a real wormlike chain of effective width w and persistence length lp in a slit of height H. For slit heights much larger than the persistence length of the polymer and much smaller than the thermal blob size, the excess free energy of the confined chain is consistent with a modified version of the scaling theory for the extended de Gennes regime in a channel that reflects the blob statistics in slit confinement. Explicitly, for channel sizes [Formula: see text], the difference between the confinement free energy of the real chain and that of an ideal chain scales like w/H.

Shape of a Stretched Polymer

The shape of a polymer plays an important role in its interactions with surrounding molecules. We characterize the shape and the orientational properties of a polymer chain under tension in a good solvent, a physical condition that is often realized both in single-molecule experiments and in vivo. Our findings reveal the existence of hitherto unobserved universal laws encompassing polymers with different rigidities and including the possible presence of excluded-volume effects, showing that both shape and orientation are solely determined by the force contribution to the free energy. In doing so, they also provide a simple way to retrieve these quantities from the knowledge of the force-versus-extension curve.

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.

Spatial confinement induces hairpins in nicked circular DNA

In living cells, DNA is highly confined in space with the help of condensing agents, DNA binding proteins and high levels of supercoiling. Due to challenges associated with experimentally studying DNA under confinement, little is known about the impact of spatial confinement on the local structure of the DNA. Here, we have used well characterized slits of different sizes to collect high resolution atomic force microscopy images of confined circular DNA with the aim of assessing the impact of the spatial confinement on global and local conformational properties of DNA. Our findings, supported by numerical simulations, indicate that confinement imposes a large mechanical stress on the DNA as evidenced by a pronounced anisotropy and tangent-tangent correlation function with respect to non-constrained DNA. For the strongest confinement we observed nanometer sized hairpins and interwound structures associated with the nicked sites in the DNA sequence. Based on these findings, we propose that spatial DNA confinement in vivo can promote the formation of localized defects at mechanically weak sites that could be co-opted for biological regulatory functions.

Hyperplectonemes: A Higher Order Compact and Dynamic DNA Self-Organization

Bacterial chromosome has a compact structure that dynamically changes its shape in response to bacterial growth rate and growth phase. Determining how chromatin remains accessible to DNA binding proteins, and transcription machinery is crucial to understand the link between genetic regulation, DNA structure, and topology. Here, we study very large supercoiled dsDNA using high-resolution characterization, theoretical modeling, and molecular dynamics calculations. We unveil a new type of highly ordered DNA organization forming in the presence of attractive DNA-DNA interactions, which we call hyperplectonemes. We demonstrate that their formation depends on DNA size, supercoiling, and bacterial physiology. We compare structural, nanomechanic, and dynamic properties of hyperplectonemes bound by three highly abundant nucleoid-associated proteins (FIS, H-NS, and HU). In all these cases, the negative supercoiling of DNA determines molecular dynamics, modulating their 3D shape. Overall, our findings provide a mechanistic insight into the critical role of DNA topology in genetic regulation.

Simulations of knotting of DNA during genome mapping

Genome mapping involves the confinement of long DNA molecules, in excess of 150 kilobase pairs, in nanochannels near the circa 50 nm persistence length of DNA. The fidelity of the map relies on the assumption that the DNA is linearized by channel confinement, which assumes the absence of knots. We have computed the probability of forming different knot types and the size of these knots for long chains (approximately 164 kilobase pairs) via pruned-enriched Rosenbluth method simulations of a discrete wormlike chain model of DNA in channel sizes ranging from 35 nm to 60 nm. Compared to prior simulations of short DNA in similar confinement, these long molecules exhibit both complex knots, with up to seven crossings, and multiple knots per chain. The knotting probability is a very strong function of channel size, ranging from 0.3% to 60%, and rationalized in the context of Odijk's theory for confined semiflexible chains. Overall, the knotting probability and knot size obtained from these equilibrium measurements are not consistent with experimental measurements of the properties of anomalously bright regions along the DNA backbone during genome mapping experiments. This result suggests that these events in experiments are either knots formed during the processing of the DNA prior to injection into the nanochannel or regions of locally high DNA concentration without a topological constraint. If so, knots during genome mapping are not an intrinsic problem for genome mapping technology.

Proteins analysed as virtual knots

Long, flexible physical filaments are naturally tangled and knotted, from macroscopic string down to long-chain molecules. The existence of knotting in a filament naturally affects its configuration and properties, and may be very stable or disappear rapidly under manipulation and interaction. Knotting has been previously identified in protein backbone chains, for which these mechanical constraints are of fundamental importance to their molecular functionality, despite their being open curves in which the knots are not mathematically well defined; knotting can only be identified by closing the termini of the chain somehow. We introduce a new method for resolving knotting in open curves using virtual knots, which are a wider class of topological objects that do not require a classical closure and so naturally capture the topological ambiguity inherent in open curves. We describe the results of analysing proteins in the Protein Data Bank by this new scheme, recovering and extending previous knotting results, and identifying topological interest in some new cases. The statistics of virtual knots in protein chains are compared with those of open random walks and Hamiltonian subchains on cubic lattices, identifying a regime of open curves in which the virtual knotting description is likely to be important.

Sorting ring polymers by knot type with modulated nanochannels

In this theoretical study we discuss a novel method for sorting ring polymers according to their topological, knotted state. The proposed approach harnesses the rich dynamical behaviour of polymers confined inside spatially-modulated nanochannels. The longitudinal mobility of the rings is shown to have two key properties that are ideally suited for knot sorting. First, at fixed topology, the mobility has an intriguing oscillatory dependence on chain length. Second, the mobility ranking of different knot types is inverted upon increasing the chain length. We show that this complex interplay of channel geometry, chain length and topology can be rationalised within a simple theoretical framework based on Fick-Jacobs's diffusive theory. The results and the interpretative scheme ought to be useful for designing microfluidic devices with optimal topological sorting capabilities.

Comparison of a stripe and slab confinement for ring and linear macromolecules in nanochannel

The combined effects of the channel asymmetry and the closed chain topology on the chain extension, structure factor, and the orientation correlations were studied using coarse-grained molecular dynamics simulations for moderate chain lengths. These effects are related to applications in linearization experiments with a DNA molecule in nanofluidic devices. According to the aspect ratio, the channels are classified as a stripe or slabs. The chain segments do not have any freedom to move in the direction of the narrowest stripe size, being approximately the same size as the segment size. The chains of both ring and linear topologies are extended more in a stripe than in a slab; this effect is strengthened for a ring. For a ring in a stripe, the extension-confinement strength dependence leads to effective Flory exponents even larger than 3/4, which is characteristic for a self-avoiding two-dimensional chain. While the chain extension-confinement strength dependence for both topologies conforms to the de Gennes regime in a stripe, a linear chain undergoes gradual transition to the pseudoideal regime as the slab height increases in the slab-like confinement. For a confined circle, the onset of the pseudoideal regime is shifted to larger slab heights. The structure factor confirms the absence of the pseudoideal and extended de Gennes regime in a stripe and the transition from the extended to the pseudoideal regime of a circular and linear chain upon increasing the slab heights.

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.

Interactions of complex polymers with nanoporous substrate

With the advance of polymer synthesis, polymers that possess unique architectures such as stars or cyclic chains, and unique chemical composition distributions such as block copolymers or statistical copolymers have become frequently encountered. Characterization of these complex polymer systems drives the development of interactive chromatography where the adsorption of polymers on the porous substrate in chromatography columns is finely tuned. Liquid Chromatography at the Critical Condition (LCCC) in particular makes use of the existence of the Critical Adsorption Point (CAP) of polymers on solid surfaces and has been successfully applied to characterization of complex polymer systems. Interpretation and understanding of chromatography behaviour of complex polymers in interactive chromatography motivates theoretical/computational studies on the CAP of polymers and partitioning of these complex polymers near the CAP. This review article covers the theoretical questions encountered in chromatographic studies of complex polymers.

The polymer physics of single DNA confined in nanochannels

In recent years, applications and experimental studies of DNA in nanochannels have stimulated the investigation of the polymer physics of DNA in confinement. Recent advances in the physics of confined polymers, using DNA as a model polymer, have moved beyond the classic Odijk theory for the strong confinement, and the classic blob theory for the weak confinement. In this review, we present the current understanding of the behaviors of confined polymers while briefly reviewing classic theories. Three aspects of confined DNA are presented: static, dynamic, and topological properties. The relevant simulation methods are also summarized. In addition, comparisons of confined DNA with DNA under tension and DNA in semidilute solution are made to emphasize universal behaviors. Finally, an outlook of the possible future research for confined DNA is given.

Dimensional reduction of duplex DNA under confinement to nanofluidic slits

There has been much interest in the dimensional properties of double-stranded DNA (dsDNA) confined to nanoscale environments as a problem of fundamental importance in both biological and technological fields. This has led to a series of measurements by fluorescence microscopy of single dsDNA molecules under confinement to nanofluidic slits. Despite the efforts expended on such experiments and the corresponding theory and simulations of confined polymers, a consistent description of changes of the radius of gyration of dsDNA under strong confinement has not yet emerged. Here, we perform molecular dynamics (MD) simulations to identify relevant factors that might account for this inconsistency. Our simulations indicate a significant amplification of excluded volume interactions under confinement at the nanoscale due to the reduction of the effective dimensionality of the system. Thus, any factor influencing the excluded volume interaction of dsDNA, such as ionic strength, solution chemistry, and even fluorescent labels, can greatly influence the dsDNA size under strong confinement. These factors, which are normally less important in bulk solutions of dsDNA at moderate ionic strengths because of the relative weakness of the excluded volume interaction, must therefore be under tight control to achieve reproducible measurements of dsDNA under conditions of dimensional reduction. By simulating semi-flexible polymers over a range of parameter values relevant to the experimental systems and exploiting past theoretical treatments of the dimensional variation of swelling exponents and prefactors, we have developed a novel predictive relationship for the in-plane radius of gyration of long semi-flexible polymers under slit-like confinement. Importantly, these analytic expressions allow us to estimate the properties of dsDNA for the experimentally and biologically relevant range of contour lengths that is not currently accessible by state-of-the-art MD simulations.

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.

Knotting dynamics of DNA chains of different length confined in nanochannels

Langevin dynamics simulations are used to characterize the typical mechanisms governing the spontaneous tying, untying and the dynamical evolution of knots in coarse-grained models of DNA chains confined in nanochannels. In particular we focus on how these mechanisms depend on the chain contour length, Lc, at a fixed channel width D = 56 nm corresponding to the onset of the Odijk scaling regime where chain backfoldings and hence knots are disfavoured but not suppressed altogether. We find that the lifetime of knots grows significantly with Lc, while that of unknots varies to a lesser extent. The underlying kinetic mechanisms are clarified by analysing the evolution of the knot position along the chain. At the considered confinement, in fact, knots are typically tied by local backfoldings of the chain termini where they are eventually untied after a stochastic motion along the chain. Consequently, the lifetime of unknots is mostly controlled by backfoldings events at the chain ends, which is largely independent of Lc. The lifetime of knots, instead, increases significantly with Lc because knots can, on average, travel farther along the chain before being untied. The observed interplay of knots and unknots lifetimes underpins the growth of the equilibrium knotting probability of longer and longer chains at fixed channel confinement.

Effects of excluded volume and hydrodynamic interaction on the deformation, orientation and motion of ring polymers in shear flow

A ring polymer is a classical model to explore the behaviors of biomacromolecules. Compared with its linear counterpart in shear flow, the ring polymer should be more sensitive to excluded volume and hydrodynamic interaction attributed to the absence of chain ends. We carried out multiparticle collision dynamics combined with molecular dynamics simulation to study the effects of excluded volume and hydrodynamic interaction on the behaviors of ring polymers in shear flow. The results show that in the absence of the strong excluded volume interaction, the ring polymer prefers a two-strand linear conformation with high deformation and orientation in the flow-gradient plane, and the tank-treading motion is nearly negligible. Ring polymers without excluded volume show no significant difference from linear polymers in the scaling exponents for the deformation, orientation and tumbling motion. We also observed that the hydrodynamic interaction could efficiently slow down the relaxation of ring polymers while the scaling exponents against the Weissenberg number have rarely been affected.

Stretching self-entangled DNA molecules in elongational fields

We present experiments of self-entangled DNA molecules stretching under a planar elongational field, and their stretching dynamics are compared to identical molecules without entanglements. Self-entangled molecules stretch in a stage-wise fashion, persisting in an "arrested" state for decades of strain prior to rapidly stretching, slowing down the stretching dynamics by an order of magnitude compared to unentangled molecules. Self-entangled molecules are shown to proceed through a transient state where one or two ends of the molecule are protruding from an entangled, knotted core. This phenomenon sharply contrasts with the wide array of transient configurations shown here and by others for stretching polymers without entanglements. The rate at which self-entangled molecules stretch through this transient state is demonstrably slower than unentangled molecules, providing the first direct experimental evidence of a topological friction. These experimental observations are shown to be qualitatively and semi-quantitatively reproduced by a dumbbell model with two fitting parameters, the values of which are reasonable in light of previous experiments of knotted DNA.

Stripe to slab confinement for the linearization of macromolecules in nanochannels

We investigated the recently suggested advantageous analysis of chain linearization experiments with macromolecules confined in a stripe-like channel (Huang and Battacharya, EPL, 2014, 106, 18004) using Monte Carlo simulations. The enhanced chain extension in a stripe, which is due to the significant excluded volume interactions between the monomers in two dimensions, weakens considerably on transition to an experimentally feasible slit-like channel. Based on the chain extension-confinement strength dependence and the structure factor behavior for a chain in a stripe, we infer the excluded volume regime (de Gennes regime) typical for two-dimensional systems. On widening of the stripe in a direction perpendicular to the stripe plane, i.e. on the transition to the slab geometry, the advantageous chain extension decreases and a Gaussian regime is observed for not very long semiflexible chains. The evidence for pseudo-ideality in confined chains is based on four indicators: the extension curves, variation of the extension with the persistence length P, estimated limits for the regimes in the investigated systems, and the structure factor behavior. The slab behavior can be observed when the two-dimensional stripe (originally of a one-monomer thickness) reaches a reduced thickness D larger than approximately D/P ≈ 0.2 in the third dimension. This maximum height of a slab at which the advantage of a stripe is retained is very low and has implications for DNA linearization experiments.

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.

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.

Computational study on the progressive factorization of composite polymer knots into separated prime components

Using Monte Carlo simulations and advanced knot localization methods, we analyze the length and distribution of prime components in composite knots tied on freely jointed rings. For increasing contour length, we observe the progressive factorization of composite knots into separated prime components. However, we observe that a complete factorization, equivalent to the "decorated ring" picture, is not obtained even for rings of contour lengths N ≃ 3 N(0), about tens of times the most probable length of the prime knots tied on the rings. The decorated ring hypothesis has been used in the literature to justify the factorization of composite knot probabilities into the knotting probabilities of their prime components. Following our results, we suggest that such a hypothesis may not be necessary to explain the factorization of the knotting probabilities, at least when polymers excluding volume is not relevant. We rationalize the behavior of the system through a simple one-dimensional model in which prime knots are replaced by slip links randomly placed on a circle, with the only constraint being that the length of the loops has the same distribution as that of the length of the corresponding prime knots.

Circular Kardar-Parisi-Zhang equation as an inflating, self-avoiding ring polymer

We consider the Kardar-Parisi-Zhang equation for a circular interface in two dimensions, unconstrained by the standard small-slope and no-overhang approximations. Numerical simulations using an adaptive scheme allow us to elucidate the complete time evolution as a crossover between a short-time regime with the interface fluctuations of a self-avoiding ring or two-dimensional vesicle, and a long-time regime governed by the Tracy-Widom distribution expected for this geometry. For small-noise amplitudes, scaling behavior is only of the latter type. Large noise is also seen to renormalize the bare physical parameters of the ring, akin to analogous parameter renormalization for equilibrium three-dimensional membranes. Our results bear particular importance on the relation between relevant universality classes of scale-invariant systems in two dimensions.

Knotting of linear DNA in nano-slits and nano-channels: a numerical study

The amount and type of self-entanglement of DNA filaments is significantly affected by spatial confinement, which is ubiquitous in biological systems. Motivated by recent advancements in single DNA molecule experiments based on nanofluidic devices and by the introduction of algorithms capable of detecting knots in open chains, we investigate numerically the entanglement of linear, open DNA chains confined inside nano-slits. The results regard the abundance, type, and length of occurring knots and are compared with recent findings for DNA inside nano-channels. In both cases, the width of the confining region, D, spans the 30 nm-1 μm range and the confined DNA chains are 1-4 μm long. It is found that the knotting probability is maximum for slit widths in the 70-100 nm range. However, over the considered DNA contour lengths, the maximum incidence of knots remains below 20%, while for channel confinement it tops 50%. Further differences of the entanglement are seen for the average contour length of the knotted region, which drops significantly below D ~100 nm for channel-confinement, while it stays approximately constant for slit-like confinement. These properties ought to reverberate in different kinetic properties of linear DNA depending on confinement and could be detectable experimentally or exploitable in nano-technological applications.