Dynamics of a polymer chain translocating through varying cone-shaped channels

Polymer translocation: effects of periodically driven confinement

We study the influence of confinement on the dynamics of translocation of a linear polymer chain in a good solvent through a cone-shaped pore. Using the Langevin dynamics simulations, we calculate both the first attempt time and translocation time as a function of the position of the back wall and apex angle α. As the in vivo confining environment is inherently dynamic, we extended the present study to explore the consequences of a periodically driven back wall and apex angles on the translocation dynamics. Our findings reveal that the translocation time initially decreases as the driving frequency increases, but increases after a certain frequency. The frequency at which the translocation time is found to be minimum is referred to as the resonance activation. Analyzing the distribution of translocation times around this frequency renders interesting information about the translocation process. We further explore the translocation dynamics by calculating the residence time of individual monomers, shedding light on the microscopic aspects of the process.

Polymer translocation: Effects of confinement

We investigate the influence of varying confinement on the dynamics of polymer translocation through a cone-shaped channel. For this, a linear polymer chain is modeled using self-avoiding walks on a square lattice. The cis side of a cone-shaped channel has a finite volume, while the trans side has a semi-infinite space. The confining environment is varied either by changing the position of the back wall while keeping the apex angle fixed or altering the apex angle while keeping the position of the back wall fixed. In both cases, the effective space ϕ, which represents the number of monomers in a chain relative to the total number of accessible sites within the cone, is reduced due to the imposed confinement. Consequently, the translocation dynamics are affected. We analyze the entropy of the confined system as a function of ϕ, which exhibits nonmonotonic behavior. We also calculate the free energy associated with the confinement as a function of a virtual coordinate for different positions of the back wall (base of the cone) along the conical axis for various apex angles. Employing the Fokker-Planck equation, we calculate the translocation time as a function of ϕ for different solvent conditions across the channel. Our findings indicate that the translocation time decreases as ϕ increases, but it eventually reaches a saturation point at a certain value of ϕ. Moreover, we highlight the possibility of controlling the translocation dynamics by manipulating the solvent quality across the channel. Furthermore, our investigation delves into the intricacies of polymer translocation through a cone-shaped channel, considering both repulsive and neutral interactions with the channel wall. This exploration unveils nuanced dynamics and sheds light on the factors that significantly impact translocation within confined channels.

Polymer Translocation and Nanopore Sequencing: A Review of Advances and Challenges

Various biological processes involve the translocation of macromolecules across nanopores; these pores are basically protein channels embedded in membranes. Understanding the mechanism of translocation is crucial to a range of technological applications, including DNA sequencing, single molecule detection, and controlled drug delivery. In this spirit, numerous efforts have been made to develop polymer translocation-based sequencing devices, these efforts include findings and insights from theoretical modeling, simulations, and experimental studies. As much as the past and ongoing studies have added to the knowledge, the practical realization of low-cost, high-throughput sequencing devices, however, has still not been realized. There are challenges, the foremost of which is controlling the speed of translocation at the single monomer level, which remain to be addressed in order to use polymer translocation-based methods for sensing applications. In this article, we review the recent studies aimed at developing control over the dynamics of polymer translocation through nanopores.

Delayed collapse transitions in a pinned polymer system

Employing Langevin dynamics simulations, we investigated the kinetics of the collapse transition for a polymer of length N when a particular monomer at a position 1≤X≤N is pinned. The results are compared with the kinetics of a free polymer. The equilibrium θ-point separating the coil from the globule phase is located by a crossover in 〈R_{g}^{2}〉/N plots of different chain lengths. Our simulation supports a three-stage mechanism for free and pinned polymer collapse: the formation of pearls, the coarsening of pearls, and the formation of a compact globule. Pinning the central monomer has negligible effects on the kinetics as it does not break the symmetry. However, pinning a monomer elsewhere causes the process to be delayed by a constant factor ϕ_{X} depending linearly upon X. The total collapse time scales with N as τ_{c}∼ϕ_{X}N^{1.60±0.03}, which implies τ_{c} is maximum when an end monomer is pinned (X=1 or N), while when pinning the central monomer (X=N/2) it is minimum and identical to that of a free polymer. The average cluster size N_{c}(t) grows in time as t^{z}, where z=1.00±0.04 for a free particle, whereas we identify two time regimes separated by a plateau for pinned polymers. At longer times, z=1.00±0.04, while it deviates in early time regimes significantly, depending on the value of X.

Diameter Dependent Melting and Softening of dsDNA Under Cylindrical Confinement

Carbon nanotubes (CNTs) are considered promising candidates for biomolecular confinement, including DNA encapsulation for gene delivery. Threshold values of diameters have been reported for double-stranded DNA (dsDNA) encapsulation inside CNTs. We have performed all-atom molecular dynamics (MD) simulations of dsDNAs confined inside single-walled CNTs (SWCNTs) at the physiologically relevant temperature of 300 K. We found that the dsDNA can be confined without being denatured only when the diameter of the SWCNT exceeds a threshold value. Below this threshold diameter, the dsDNA gets denatured and melts even at the temperature of 300 K. Our simulations using SWCNTs with chirality indices (20,20) to (30,30) at 300 K found the critical diameter to be 3.25 nm (corresponding to (24,24) chirality). Analyses of the hydrogen bonds (H-bonds), Van der Walls (VdW) energy, and other inter-base interactions show drastic reduction in the number of H-bonds, VdW energy, and electrostatic energies between the bases of dsDNA when it is confined in narrower SWCNTs (up to diameter of 3.12 nm). On the other hand, the higher interaction energy between the dsDNA and the SWCNT surface in narrower SWCNTs assists in the melting of the dsDNA. Electrostatic mapping and hydration status analyses show that the dsDNA is not adequately hydrated and the counter ion distribution is not uniform below the critical diameter of the SWCNT. As properly hydrated counter ions provide stability to the dsDNA, we infer that the inappropriate hydration of counter ions and their non-uniform distribution around the dsDNA cause the melting of the dsDNA inside SWCNTs of diameter below the critical value of 3.25 nm. For confined dsDNAs that do not get denatured, we computed their elastic properties. The persistence length of dsDNA was found to increase by a factor of about two and the torsional stiffness by a factor of 1.5 for confinement inside SWCNTs of diameters up to 3.79 nm, the stretch modulus also following nearly the same trend. Interestingly, for higher diameters of SWCNT, 3.79 nm and above, the dsDNA becomes more flexible, demonstrating that the mechanical properties of the dsDNA under cylindrical confinement depend non-monotonically on the confinement diameter.

Melting of confined DNA: static and dynamic properties

We study dsDNA (double strand DNA) melting in detail within varying strip-like confinement in a two-dimensional lattice model. The interplay between reduced configurational entropy and attractive base-pairing energy results in a non-monotonic melting profile of DNA. Structural transitions associated with confined DNA melting reveal a stretched or extended state for very strong confinement. By using the exact enumeration method, we investigate the emergence of a local denatured zone e.g. bubbles during DNA melting. The survival time of a single bubble within varying strip width is studied from the Fokker-Planck formalism by considering the bubble size as a reaction co-ordinate. We show that a simple lattice model can capture the sequence heterogeneity effect on DNA melting and bubble dynamics within the strip. Different time scales of bubble zipping for different DNA sequences are found, which may have potential applications in denaturation mapping.

Effect of solvent gradient inside the entropic trap on polymer migration

By employing the exact enumeration technique on the lattice model of a polymer, we study the migration of the polymer chain across an entropic trap in a quasiequilibrium condition and explore the effect of solvent gradient present in the entropic trap which acts both parallel and perpendicular to the direction of migration. The Fokker-Planck formalism utilizes the free energy landscape of a polymer chain across the channel in the presence of the entropic trap to calculate the migration time. It is revealed that the migration is fast when the solvent gradient acts along the migration axis (i.e., x axis) inside the channel in comparison to the channel having the entropic trap. We report here for the first time that the entropic trap makes the migration faster at a certain value of solvent gradient. We also study the effect of transverse solvent gradient (along the y axis) inside the trap and investigate the structural changes of the polymer during migration through the channel. We observe the nonmonotonic dependence of migration time on the solvent gradient.