Metastable Tight Knots in Semiflexible Chains

Effect of Counterion Size on Knotted Polyelectrolyte Conformations

Using Langevin dynamics simulations and a coarse-grained primitive model of electrolytes, we show that the behavior of knotted circular strong polyelectrolytes (PEs) in diluted aqueous solution is largely affected by the diameter of the counterions (CIs), σCI. Indeed, we observe that both gyration radius and knot length vary nonmonotonically with σCI, with both small and bulky CIs favoring knot localization, while medium-sized ones promote delocalized knots. We also show that the conformational change from delocalized to tight knots occurs via the progressive coalescence of the knot's essential crossings. The emerging conformers correspond to the minima of the free energy landscape profiled as a function of the knot length or PE size. We demonstrate that different conformational states can coexist, the transition between them appearing first-order-like and controlled by the enthalpic and entropic trade-off of the amount of CIs condensed on the PE. Such balance can be further altered by varying CI concentrations, thus providing an additional and more convenient tuning parameter for the system properties. Our results lay the foundation for achieving broader and more precise external adjustability of knotted PE size and shape by choosing the nature of its CIs. Thus, they offer new intriguing possibilities for designing novel PE-based materials that are capable of responding to changes in ionic solution properties.

Expansion of Single Chains Released from a Spherical Cavity

A two-stage model is developed to explain the phenomena of chain expansion, released from a confining cavity. In the first stage, the chain is assumed to expand as a sphere, while in the second stage it expands like a coil. The kinetic equations for the variation of chain size are derived in the two stages by balancing the rate of the free energy change with the rate of the energy dissipation. Langevin dynamics simulations are then performed to examine the theory. We find that the expansion process is dominated by the second stage and the evolution of chain size follows, mainly, the predicted curve for coil expansion, which depends on the chain length and is not sensitive to the confining volume fraction. It permits to define the expansion time for the process. Further study reveals that the chain does undergo a spherical expansion in the first stage with the characteristic time much shorter than the one for the second stage. As a consequence, the first-stage variation of chain size can be regarded as an add-on to the principal curve of expansion designated by the second stage. The scaling behaviors and the associated scaling exponents are analyzed in details. The simulation results well support the theory.

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.

Molecular weaving

Historically, the interlacing of strands at the molecular level has mainly been limited to coordination polymers and DNA. Despite being proposed on a number of occasions, the direct, bottom-up assembly of molecular building blocks into woven organic polymers remained an aspirational, but elusive, target for several decades. However, recent successes in two-dimensional and three-dimensional molecular-level weaving now offer new opportunities and research directions at the interface of polymer science and molecular nanotopology. This Perspective provides an overview of the features and potential of the periodic nanoscale weaving of polymer chains, distinguishing it from randomly entangled polymer networks and rigid crystalline frameworks. We review the background and experimental progress so far, and conclude by considering the potential of molecular weaving and outline some of the current and future challenges in this emerging field.

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.

DNA Knot Malleability in Single-Digit Nanopores

Knots in long DNA molecules are prevalent in biological systems and serve as a model system for investigating static and dynamic properties of biopolymers. We explore the dynamics of knots in double-stranded DNA in a new regime of nanometer-scale confinement, large forces, and short time scales, using solid-state nanopores. We show that DNA knots undergo isomorphic translocation through a nanopore, retaining their equilibrium morphology by swiftly compressing in a lateral direction to fit the constriction. We observe no evidence of knot tightening or jamming, even for single-digit nanopores. We explain the observations as the malleability of DNA, characterized by sharp buckling of the DNA in nanopores, driven by the transient disruption of base pairing. Our molecular dynamics simulations support the model. These results are relevant not only for the understanding of DNA packing and manipulation in living cells but also for the polymer physics of DNA and the development of nanopore-based sequencing technologies.

Comparing equilibration schemes of high-molecular-weight polymer melts with topological indicators

Recent theoretical studies have demonstrated that the behaviour of molecular knots is a sensitive indicator of polymer structure. Here, we use knots to verify the ability of two state-of-the-art algorithms-configuration assembly and hierarchical backmapping-to equilibrate high-molecular-weight (MW) polymer melts. Specifically, we consider melts with MWs equivalent to several tens of entanglement lengths and various chain flexibilities, generated with both strategies. We compare their unknotting probability, unknotting length, knot spectra, and knot length distributions. The excellent agreement between the two independent methods with respect to knotting properties provides an additional strong validation of their ability to equilibrate dense high-MW polymeric liquids. By demonstrating this consistency of knotting behaviour, our study opens the way for studying topological properties of polymer melts beyond time and length scales accessible to brute-force molecular dynamics simulations.

Can Soft Models Describe Polymer Knots?

Similar to macroscopic ropes and cables, long polymers create knots. We address the fundamental question whether and under which conditions it is possible to describe these intriguing objects with crude models that capture only mesoscale polymer properties. We focus on melts of long polymers which we describe by a model typical for mesoscopic simulations. A worm-like chain model defines the polymer architecture. To describe nonbonded interactions, we deliberately choose a generic "soft" repulsive potential that leads to strongly overlapping monomers and coarse local liquid structure. The soft model is parametrized to accurately reproduce mesoscopic structure and conformations of reference polymer melts described by a microscopic model. The microscopically resolved samples retain all generic features affecting polymer topology and provide, therefore, reliable reference data on knots. We compare characteristic knotting properties in mesoscopic and microscopically resolved melts for different cases of chain stiffness. We conclude that mesoscopic models can reliably describe knots in those melts, where the length scale characterizing polymer stiffness is substantially larger than the size of monomer-monomer excluded volume. In this case, simplified local liquid structure influences knotting properties only marginally. In contrast, mesoscopic models perform poorly in melts with flexible chains. We qualitatively explain our findings through a free energy model of simple knots available in the literature.

Topological Disentanglement of Linear Polymers under Tension

We develop a theoretical description of the topological disentanglement occurring when torus knots reach the ends of a semiflexible polymer under tension. These include decays into simpler knots and total unknotting. The minimal number of crossings and the minimal knot contour length are the topological invariants playing a key role in the model. The crossings behave as particles diffusing along the chain and the application of appropriate boundary conditions at the ends of the chain accounts for the knot disentanglement. Starting from the number of particles and their positions, suitable rules allow reconstructing the type and location of the knot moving on the chain Our theory is extensively benchmarked with corresponding molecular dynamics simulations and the results show a remarkable agreement between the simulations and the theoretical predictions of the model.

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.

Identifying knot types of polymer conformations by machine learning

We investigate the use of artificial neural networks (NNs) as an alternative tool to current analytical methods for recognizing knots in a given polymer conformation. The motivation is twofold. First, it is of interest to examine whether NNs are effective at learning the global and sequential properties that uniquely define a knot. Second, knot classification is an important and unsolved problem in mathematical and physical sciences, and NNs may provide insights into this problem. Motivated by these points, we generate millions of polymer conformations for five knot types: 0, 3_{1}, 4_{1}, 5_{1}, and 5_{2}, and we design various NN models for classification. Our best model achieves a five-class classification accuracy of above 99% on a polymer of 100 monomers. We find that the sequential modeling ability of recurrent NNs is crucial for this result, as it outperforms feed-forward NNs and successfully generalizes to differently sized conformations as well. We present our methods and suggest that deep learning may be used in specific applications of knot detection where some error is permissible. Hopefully, with further development, NNs can offer an alternative computational method for knot identification and facilitate knot research in mathematical and physical sciences.

Complex DNA knots detected with a nanopore sensor

Equilibrium knots are common in biological polymers-their prevalence, size distribution, structure, and dynamics have been extensively studied, with implications to fundamental biological processes and DNA sequencing technologies. Nanopore microscopy is a high-throughput single-molecule technique capable of detecting the shape of biopolymers, including DNA knots. Here we demonstrate nanopore sensors that map the equilibrium structure of DNA knots, without spurious knot tightening and sliding. We show the occurrence of both tight and loose knots, reconciling previous contradictory results from different experimental techniques. We evidence the occurrence of two quantitatively different modes of knot translocation through the nanopores, involving very different tension forces. With large statistics, we explore the complex knots and, for the first time, reveal the existence of rare composite knots. We use parametrized complexity, in concert with simulations, to test the theoretical assumptions of the models, further asserting the relevance of nanopores in future investigation of knots.

Conformational State Hopping of Knots in Tensioned Polymer Chains

We use Brownian dynamics simulations to study the conformational states of knots on tensioned chains. Focusing specifically on the 8₁ knot, we observe knot conformational state hopping and show that the process can be described by a two-state kinetic model in the presence of an external force. The distribution of knot conformational states depends on the applied chain tension, which leads to a force-dependent distribution of knot untying pathways. We generalize our findings by considering the untying pathways of other knots and find that the way knots untie is generally governed by the force applied to the chain. From a broader perspective, being able to influence how a knot unties via external force can potentially be useful for applications of single-molecule techniques in which knots are unwanted.

Topological Disentanglement Dynamics of Torus Knots on Open Linear Polymers

We simulate and study the topological disentanglement occurring when torus knots reach the ends of a semiflexible open polymer (decay into simpler knots or unknotting). Through a rescaling procedure and the application of appropriate boundary conditions, we show that the full unknotting process can be understood in terms of point-like particles representing essential crossings, diffusing on the support [0, 1]. We address the bending and configurational free energy drives on the diffusion process, together with the scaling properties of the effective diffusion and friction coefficients. Agreement with simulations suggests universal features for these two model parameters.

The folding pathways and thermodynamics of semiflexible polymers

Inspired by the protein folding and DNA packing, we have systematically studied the thermodynamic and kinetic behaviors of single semiflexible homopolymers by Langevin dynamics simulations. In line with experiments, a rich variety of folding products, such as rod-like bundles, hairpins, toroids, and a mixture of them, are observed in the complete diagram of states. Moreover, knotted structures with a significant population are found in a certain range of bending stiffness in thermal equilibrium. As the solvent quality becomes poorer, the population of the intermediate occurring in the folding process increases, which leads to a severe chevron rollover for the folding arm. However, the population of the intermediates in the unfolding process is very low, insufficient to induce unfolding arm rollover. The total types of folding pathways from the coil state to the toroidal state for a semiflexible polymer chain remain unchanged by varying the solvent quality or temperature, whereas the kinetic partitioning into different folding events can be tuned significantly. In the process of knotting, three types of mechanisms, namely, plugging, slipknotting, and sliding, are discovered. Along the folding evolution, a semiflexible homopolymer chain can knot at any stage of folding upon leaving the extended coil state, and the probability to find a knot increases with chain compactness. In addition, we find rich types of knotted topologies during the folding of a semiflexible homopolymer chain. This study should be helpful in gaining insight into the general principles of biopolymer folding.

Motion of Knots in DNA Stretched by Elongational Fields

Knots in DNA occur in biological systems, serve as a model system for polymer entanglement, and affect the efficacy of modern genomics technologies. We study the motion of complex knots in DNA by stretching molecules with a divergent electric field that provides an elongational force. We demonstrate that the motion of knots is nonisotropic and driven towards the closest end of the molecule. We show for the first time experimentally that knots can go from a mobile to a jammed state by varying an applied strain rate, and that this jamming is reversible. We measure the mobility of knots as a function of strain rate, demonstrating the conditions under which knots can be driven towards the ends of the molecule and untied.

A nanofluidic knot factory based on compression of single DNA in nanochannels

Knots form when polymers self-entangle, a process enhanced by compaction with important implications in biological and artificial systems involving chain confinement. In particular, new experimental tools are needed to assess the impact of multiple variables influencing knotting probability. Here, we introduce a nanofluidic knot factory for efficient knot formation and detection. Knots are produced during hydrodynamic compression of single DNA molecules against barriers in a nanochannel; subsequent extension of the chain enables direct assessment of the number of independently evolving knots. Knotting probability increases with chain compression as well as with waiting time in the compressed state. Using a free energy derived from scaling arguments, we develop a knot-formation model that can quantify the effect of interactions and the breakdown of Poisson statistics at high compression. Our model suggests that highly compressed knotted states are stabilized by a decreased free energy as knotted contour contributes a lower self-exclusion derived free energy.

Knots modify the coil-stretch transition in linear DNA polymers

We perform single-molecule DNA experiments to investigate the relaxation dynamics of knotted polymers and examine the steady-state behavior of knotted polymers in elongational fields. The occurrence of a knot reduces the relaxation time of a molecule and leads to a shift in the molecule's coil-stretch transition to larger strain rates. We measure chain extension and extension fluctuations as a function of strain rate for unknotted and knotted molecules. The curves for knotted molecules can be collapsed onto the unknotted curves by defining an effective Weissenberg number based on the measured knotted relaxation time in the low extension regime, or a relaxation time based on Rouse/Zimm scaling theories in the high extension regime. Because a knot reduces a molecule's relaxation time, we observe that knot untying near the coil-stretch transition can result in dramatic changes in the molecule's conformation. For example, a knotted molecule at a given strain rate can experience a stretch-coil transition, followed by a coil-stretch transition, after the knot partially or fully unties.

Steady-State and Transient Behavior of Knotted Chains in Extensional Fields

Recently, there has been a push to understand how molecular topology alters the nonequilibrium dynamics of polymer systems. In this paper, we probe how knotted polymers evolve in planar extensional fields using Brownian dynamics simulations and single-molecule experiments. In the first part of the study, we quantify the extension versus strain-rate curves of polymers and find that knots shift these curves to larger strain-rates. These trends can be quantitatively explained by Rouse-like scaling theories. In the second half of the study, we examine the consequences of knot untying on the time-dependent conformations of polymers in these external fields. We find that knot untying creates significant, transient changes in chain extension. If the topology is complex, the chain undergoes a wide range of time-dependent conformations since knot untying proceeds through many different stages. We provide examples of such untying trajectories over time.

The energy cost of polypeptide knot formation and its folding consequences

Knots are natural topologies of chains. Yet, little is known about spontaneous knot formation in a polypeptide chain—an event that can potentially impair its folding—and about the effect of a knot on the stability and folding kinetics of a protein. Here we used optical tweezers to show that the free energy cost to form a trefoil knot in the denatured state of a polypeptide chain of 120 residues is 5.8 ± 1 kcal mol⁻¹. Monte Carlo dynamics of random chains predict this value, indicating that the free energy cost of knot formation is of entropic origin. This cost is predicted to remain above 3 kcal mol⁻¹ for denatured proteins as large as 900 residues. Therefore, we conclude that naturally knotted proteins cannot attain their knot randomly in the unfolded state but must pay the cost of knotting through contacts along their folding landscape.

The effect of knots on protein stability and folding kinetics is not well understood. Here the authors combine optical tweezer experiments and calculations to experimentally determine the energy cost for knot formation, which indicates that knotted proteins evolved specific folding pathways because knot formation in unfolded chains is unfavorable.

Pore translocation of knotted DNA rings

We use an accurate coarse-grained model for DNA and stochastic molecular dynamics simulations to study the pore translocation of 10-kbp-long DNA rings that are knotted. By monitoring various topological and physical observables we find that there is not one, as previously assumed, but rather two qualitatively different modes of knot translocation. For both modes the pore obstruction caused by knot passage has a brief duration and typically occurs at a late translocation stage. Both effects are well in agreement with experiments and can be rationalized with a transparent model based on the concurrent tensioning and sliding of the translocating knotted chains. We also observed that the duration of the pore obstruction event is more controlled by the knot translocation velocity than the knot size. These features should advance the interpretation and design of future experiments aimed at probing the spontaneous knotting of biopolymers.

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.

Direct observation of DNA knots using a solid-state nanopore

Long DNA molecules can self-entangle into knots. Experimental techniques for observing such DNA knots (primarily gel electrophoresis) are limited to bulk methods and circular molecules below 10 kilobase pairs in length. Here, we show that solid-state nanopores can be used to directly observe individual knots in both linear and circular single DNA molecules of arbitrary length. The DNA knots are observed as short spikes in the nanopore current traces of the traversing DNA molecules and their detection is dependent on a sufficiently high measurement resolution, which can be achieved using high-concentration LiCl buffers. We study the percentage of molecules with knots for DNA molecules of up to 166 kilobase pairs in length and find that the knotting occurrence rises with the length of the DNA molecule, consistent with a constant knotting probability per unit length. Our experimental data compare favourably with previous simulation-based predictions for long polymers. From the translocation time of the knot through the nanopore, we estimate that the majority of the DNA knots are tight, with remarkably small sizes below 100 nm. In the case of linear molecules, we also observe that knots are able to slide out on application of high driving forces (voltage).

Role of Bending Energy and Knot Chirality in Knot Distribution and Their Effective Interaction along Stretched Semiflexible Polymers

Knots appear frequently in semiflexible (bio)polymers, including double-stranded DNA, and their presence can affect the polymer’s physical and functional properties. In particular, it is possible and indeed often the case that multiple knots appear on a single chain, with effects which have only come under scrutiny in the last few years. In this manuscript, we study the interaction of two knots on a stretched semiflexible polymer, expanding some recent results on the topic. Specifically, we consider an idealization of a typical optical tweezers experiment and show how the bending rigidity of the chain—And consequently its persistence length—Influences the distribution of the entanglements; possibly more importantly, we observe and report how the relative chirality of the otherwise identical knots substantially modifies their interaction. We analyze the free energy of the chain and extract the effective interactions between embedded knots, rationalizing some of their pertinent features by means of simple effective models. We believe the salient aspect of the knot–knot interactions emerging from our study will be present in a large number of semiflexible polymers under tension, with important consequences for the characterization and manipulation of these systems—Be they artificial or biologica in origin—And for their technological application.

A Monte Carlo Study of Knots in Long Double-Stranded DNA Chains

We determine knotting probabilities and typical sizes of knots in double-stranded DNA for chains of up to half a million base pairs with computer simulations of a coarse-grained bead-stick model: Single trefoil knots and composite knots which include at least one trefoil as a prime factor are shown to be common in DNA chains exceeding 250,000 base pairs, assuming physiologically relevant salt conditions. The analysis is motivated by the emergence of DNA nanopore sequencing technology, as knots are a potential cause of erroneous nucleotide reads in nanopore sequencing devices and may severely limit read lengths in the foreseeable future. Even though our coarse-grained model is only based on experimental knotting probabilities of short DNA strands, it reproduces the correct persistence length of DNA. This indicates that knots are not only a fine gauge for structural properties, but a promising tool for the design of polymer models.

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.

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.

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.

Topological events in single molecules of E. coli DNA confined in nanochannels

We present experimental data concerning potential topological events such as folds, internal backfolds, and/or knots within long molecules of double-stranded DNA when they are stretched by confinement in a nanochannel. Genomic DNA from E. coli was labeled near the 'GCTCTTC' sequence with a fluorescently labeled dUTP analog and stained with the DNA intercalator YOYO. Individual long molecules of DNA were then linearized and imaged using methods based on the NanoChannel Array technology (Irys® System) available from BioNano Genomics. Data were collected on 189 153 molecules of length greater than 50 kilobases. A custom code was developed to search for abnormal intensity spikes in the YOYO backbone profile along the length of individual molecules. By correlating the YOYO intensity spikes with the aligned barcode pattern to the reference, we were able to correlate the bright intensity regions of YOYO with abnormal stretching in the molecule, which suggests these events were either a knot or a region of internal backfolding within the DNA. We interpret the results of our experiments involving molecules exceeding 50 kilobases in the context of existing simulation data for relatively short DNA, typically several kilobases. The frequency of these events is lower than the predictions from simulations, while the size of the events is larger than simulation predictions and often exceeds the molecular weight of the simulated molecules. We also identified DNA molecules that exhibit large, single folds as they enter the nanochannels. Overall, topological events occur at a low frequency (∼7% of all molecules) and pose an easily surmountable obstacle for the practice of genome mapping in nanochannels.

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.

Confined polymers in the extended de Gennes regime

We show that the problem of describing the conformations of a semiflexible polymer confined to a channel can be mapped onto an exactly solvable model in the so-called extended de Gennes regime. This regime (where the polymer is neither weakly nor strongly confined) has recently been intensively studied experimentally and by means of computer simulations. The exact solution predicts precisely how the conformational fluctuations depend upon the channel width and upon the microscopic parameters characterizing the physical properties of the polymer.