Simulations of knotting of DNA during genome mapping

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.

Interactions between two knots in nanochannel-confined DNA molecules

Experimental data on the interaction between two knots in deoxyribonucleic acid (DNA) confined in nanochannels produced two particular behaviors of knot pairs along the DNA molecules: (i) widely separated knots experience an attractive interaction but only remain in close proximity for several seconds and (ii) knots tend to remain separated until one of the knots unravels at the chain end. The associated free energy profile of the knot-knot separation distance for an ensemble of DNA knots exhibits a global minimum when knots are separated, indicating that the separated knot state is more stable than the intertwined knot state, with dynamics in the separated knot state that are consistent with independent diffusion. The experimental observations of knot-knot interactions under nanochannel confinement are inconsistent with previous simulation-based and experimental results for stretched polymers under tension wherein the knots attract and then stay close to each other. This inconsistency is postulated to result from a weaker fluctuation-induced attractive force between knots under confinement when compared to the knots under tension, the latter of which experience larger fluctuations in transverse directions.

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.

Separation of Geometrical and Topological Entanglement in Confined Polymers Driven out of Equilibrium

We use Brownian dynamics simulations and advanced topological profiling methods to characterize the out-of-equilibrium evolution of self-entanglement in linear polymers confined into nanochannels and under periodic compression. By introducing suitable observables, we can distinguish two main forms of entanglement that we term geometrical and topological. The latter is measured by the number of (essential) crossings of the physical knot detected after a suitable bridging of the chain termini. The former is instead measured as the average number of times a linear chain appears to cross itself when viewed under all projections and is irrespective of the physical knotted state. The key discovery of our work is that these two forms of entanglement are uncoupled and evolve with distinct dynamics. While geometrical entanglement is typically in phase with the compression-elongation cycles and it is primarily sensitive to its force f, the topological measure is mildly sensitive to cyclic modulation but strongly depends on both compression force f and duration k. The findings could assist the interpretation of experiments using fluorescence molecular tracers to track physical knots in polymers. Furthermore, we identify optimal regions in the experimentally controllable parameter space in which to obtain more/less topological and geometrical entanglement; this may help designing polymers with targeted topology.

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.

Odijk excluded volume interactions during the unfolding of DNA confined in a nanochannel

We report experimental data on the unfolding of human and E. coli genomic DNA molecules shortly after injection into a 45 nm nanochannel. The unfolding dynamics are deterministic, consistent with previous experiments and modeling in larger channels, and do not depend on the biological origin of the DNA. The measured entropic unfolding force per friction per unit contour length agrees with that predicted by combining the Odijk excluded volume with numerical calculations of the Kirkwood diffusivity of confined DNA. The time scale emerging from our analysis has implications for genome mapping in nanochannels, especially as the technology moves towards longer DNA, by setting a lower bound for the delay time before making a measurement.