Evolution of Nested Folding States in Compression of a Strongly Confined Semiflexible Chain

Transient physics in the compression and mixing dynamics of two nanochannel-confined polymer chains

We use molecular dynamics (MD) simulation and nanofluidic experiments to probe the non-equilibrium transient physics of two nanochannel-confined polymers driven against a permeable barrier in a flow field. For chains with a persistence length P smaller than the channel diameter D, both simulation and experiment with dsDNA reveal nonuniform mixing of the two chains, with one chain dominating locally in what we term "aggregates." Aggregates undergo stochastic dynamics, persisting for a limited time, then disappearing and reforming. Whereas aggregate-prone mixing occurs immediately at sufficiently high flow speeds, chains stay segregated at intermediate flow for some time, often attempting to mix multiple times, before suddenly successfully mixing. Observation of successful mixing nucleation events in nanofluidic experiments reveal that they arise through a peculiar "back-propagation" mechanism whereby the upstream chain, closest to the barrier, penetrates and passes through the downstream chain (farthest from the barrier) moving against the flow direction. Simulations suggest that the observed back-propagation nucleation mechanism is favored at intermediate flow speeds and arises from a special configuration where the upstream chain exhibits one or more folds facing the downstream chain, while the downstream chain has an unfolded chain end facing upstream.

Knot Formation on DNA Pushed Inside Chiral Nanochannels

We performed coarse-grained molecular dynamics simulations of DNA polymers pushed inside infinite open chiral and achiral channels. We investigated the behavior of the polymer metrics in terms of span, monomer distributions and changes of topological state of the polymer in the channels. We also compared the regime of pushing a polymer inside the infinite channel to the case of polymer compression in finite channels of knot factories investigated in earlier works. We observed that the compression in the open channels affects the polymer metrics to different extents in chiral and achiral channels. We also observed that the chiral channels give rise to the formation of equichiral knots with the same handedness as the handedness of the chiral channels.

Mixing and demixing arising from compression of two semiflexible polymer chains in nanochannels

We use molecular dynamics simulation to probe the non-equilibrium physics of two nanochannel-confined semiflexible polymers in a homogeneous flow field. We find that for sufficiently stiff chains the internal organization of the two chains takes the form of interwoven folds and circular coils. This organization can lead to mixing or demixing depending on chain stiffness and flow speed. At low and intermediate flow, the two chains adopt a folded configuration, which favours mixing. At high flow, the two chains adopt a predominantly coiled configuration. The coiled configuration results in demixing when the chains are compressed from an initially demixed condition and mixing when the chains are compressed from an initially mixed condition. We find that the mixing/demixing behaviour is governed by the ratio of the number of folded segments of one chain relative to the other at low flow and by the degree of coiling in both chains at high flow. For decreasing stiffness, the chains start to aggregate locally instead of mixing smoothly at low and intermediate flow. In the limit of completely flexible chains, the two chains either completely segregate at low flow, or adopt a locally demixed configuration consisting of large aggregates of one chain relative to the other that undergo complex stochastic dynamics, diffusing, disintegrating, and reforming at intermediate flow. The transition from complete segregation to the aggregate-dominated configuration occurs when the linear intra-chain ordering breaks down.

Organized states arising from compression of single semiflexible polymer chains in nanochannels

We use molecular dynamics simulation to probe the nonequilibrium physics of single nanochannel-confined semiflexible polymers in a homogeneous flow field. The flow field compresses the polymer against the end of the nanochannel, simulating an experiment of a nanochannel confined chain compressed against a slit barrier. The flow-based compression gives rise to a packing of the chain against the channel end that possesses a striking organization, consisting of interweaving of folds and circular coils. For stiff chains at low flow, we find that the organization is dominated by repeated hairpin folds. For stiff chains at higher flow, we observe that circular coils arise along with the folds, with folding and coiling domains becoming interwoven at the highest flow speeds. Chain organization is retained even when the chain persistence length is on order of the channel width. We show that the global polymer organization, consisting of a number of defined folds and coiled loops, arises from the minimization of the total chain free energy.

Knot formation of dsDNA pushed inside a nanochannel

Recent experiments demonstrated that knots in single molecule dsDNA can be formed by compression in a nanochannel. In this manuscript, we further elucidate the underlying molecular mechanisms by carrying out a compression experiment in silico, where an equilibrated coarse-grained double-stranded DNA confined in a square channel is pushed by a piston. The probability of forming knots is a non-monotonic function of the persistence length and can be enhanced significantly by increasing the piston speed. Under compression knots are abundant and delocalized due to a backfolding mechanism from which chain-spanning loops emerge, while knots are less frequent and only weakly localized in equilibrium. Our in silico study thus provides insights into the formation, origin and control of DNA knots in nanopores.

Adsorption of Semiflexible Polymers in Cylindrical Tubes

Conformations of wormlike chains in cylindrical pores with attractive walls are explored for varying pore radius and strength of the attractive wall potential by molecular dynamics simulations of a coarse-grained model. Local quantities such as the fraction of monomeric units bound to the surface and the bond-orientational order parameter as well as the radial density distribution are studied, as well as the global chain extensions parallel to the cylinder axis and perpendicular to the cylinder surface. A nonmonotonic convergence of these properties to their counterparts for adsorption on a planar substrate is observed due to the conflict between pore surface curvature and chain stiffness. Also the interpretation of partially adsorbed chains in terms of trains, loops, and tails is discussed.

Tug of war in a double-nanopore system

We simulate a tug-of-war (TOW) scenario for a model double-stranded DNA threading through a double nanopore (DNP) system. The DNA, simultaneously captured at both pores, is subject to two equal and opposite forces -f[over ⃗]_{L}=f[over ⃗]_{R} (TOW), where f[over ⃗]_{L} and f[over ⃗]_{R} are the forces applied to the left and the right pore, respectively. Even though the net force on the DNA polymer Δf[over ⃗]_{LR}=f[over ⃗]_{L}+f[over ⃗]_{R}=0, the mean first passage time (MFPT) 〈τ〉 depends on the magnitude of the TOW forces |f_{L}|=|f_{R}|=f_{LR}. We qualitatively explain this dependence of 〈τ〉 on f_{LR} from the known results for the single-pore translocation of a triblock copolymer A-B-A with ℓ_{pB}>ℓ_{pA}, where ℓ_{pA} and ℓ_{pB} are the persistence length of the A and B segments, respectively. We demonstrate that the time of flight of a monomer with index m [〈τ_{LR}(m)〉] from one pore to the other exhibits quasiperiodic structure commensurate with the distance between the pores d_{LR}. Finally, we study the situation where we offset the TOW biases so that Δf[over ⃗]_{LR}=f[over ⃗]_{L}+f[over ⃗]_{R}≠0, and qualitatively reproduce the experimental result of the dependence of the MFPT on Δf[over ⃗]_{LR}. We demonstrate that, for a moderate bias, the MFPT for the DNP system for a chain length N follows the same scaling ansatz as that for the single nanopore, 〈τ〉=(AN^{1+ν}+η_{pore}N)(Δf_{LR})^{-1}, where η_{pore} is the pore friction, which enables us to estimate 〈τ〉 for a long chain. Our Brownian dynamics simulation studies provide fundamental insights and valuable information about the details of the translocation speed obtained from 〈τ_{LR}(m)〉, and accuracy of the translation of the data obtained in the time domain to units of genomic distances.

Compression and Stretching of Single DNA Molecules under Channel Confinement

We study the compression and extension response of single dsDNA (double-stranded DNA) molecules confined in cylindrical channels by means of Monte Carlo simulations. The elastic response of micrometer-sized DNA to the external force acting through the chain ends or through the piston is markedly affected by the size of the channel. The interpretation of the force (f)-displacement (R) functions under quasi-one-dimensional confinement is facilitated by resolving the overall change of displacement ΔR into the confinement contribution ΔR_D and the force contribution ΔR_f. The external stretching of confined DNA results in a characteristic pattern of f-R functions involving their shift to the larger extensions due to the channel-induced pre-stretching ΔR_D. A smooth end-chain compression into loop-like conformations observed in moderately confined DNA can be accounted for by the relationship valid for a Gaussian chain in bulk. In narrow channels, the considerably pre-stretched DNA molecules abruptly buckle on compression by the backfolding into hairpins. On the contrary, the piston compression of DNA is characterized by a gradual reduction of the chain span S and by smooth f-S functions in the whole spatial range from the 3d near to 1d limits. The observed discrepancy between the shape of the f-R and f-S functions from two compression methods can be important for designing nanopiston experiments of compaction and knotting of single DNA in nanochannels.