Role of non-equilibrium conformations on driven polymer translocation

Homopolymer and heteropolymer translocation through patterned pores under fluctuating forces

We investigate the translocation of a semiflexible polymer through extended patterned pores using Langevin dynamics simulations, specifically focusing on the influence of a time-dependent driving force. Our findings reveal that, akin to its flexible counterpart, a rigid chain-like molecule translocates faster when subjected to an oscillating force than a constant force of equivalent average magnitude. The enhanced translocation is strongly correlated with the stiffness of the polymer and the stickiness of the pores. The arrangement of the pores plays a pivotal role in translocation dynamics, deeply influenced by the interplay between polymer stiffness and pore-polymer interactions. For heterogeneous polymers with periodically varying stiffness, the oscillating force introduces significant variations in the translocation time distributions based on segment sizes and orientations. On the basis of these insights, we propose a sequencing approach that harnesses distinct pore surface properties that are capable of accurately predicting sequences in heteropolymers with diverse bending rigidities.

Gain reversal in the translocation dynamics of a semiflexible polymer through a flickering pore

We study the driven translocation of a semiflexible polymer through an attractive extended pore with a periodically oscillating width. Similar to its flexible counterpart, a stiff polymer translocates through an oscillating pore more quickly than a static pore whose width is equal to the oscillating pore's mean width. This efficiency quantified as a gain in the translocation time, highlights a considerable dependence of the translocation dynamics on the stiffness of the polymer and the attractive nature of the pore. The gain characteristics for various polymer stiffness exhibit a trend reversal when the stickiness of the pore is changed. The gain reduces with increasing stiffness for a lower attractive strength of the pore, whereas it increases with increasing stiffness for higher attractive strengths. Such a dependence leads to the possibility of a high degree of robust selectivity in the translocation process.

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.

Driven translocation of a semiflexible polymer through a conical channel in the presence of attractive surface interactions

We study the translocation of a semiflexible polymer through a conical channel with attractive surface interactions and a driving force which varies spatially inside the channel. Using the results of the translocation dynamics of a flexible polymer through an extended channel as control, we first show that the asymmetric shape of the channel gives rise to non-monotonic features in the total translocation time as a function of the apex angle of the channel. The waiting time distributions of individual monomer beads inside the channel show unique features strongly dependent on the driving force and the surface interactions. Polymer stiffness results in longer translocation times for all angles of the channel. Further, non-monotonic features in the translocation time as a function of the channel angle changes substantially as the polymer becomes stiffer, which is reflected in the changing features of the waiting time distributions. We construct a free energy description of the system incorporating entropic and energetic contributions in the low force regime to explain the simulation results.

How capture affects polymer translocation in a solitary nanopore

DNA capture with high fidelity is an essential part of nanopore translocation. We report several important aspects of the capture process and subsequent translocation of a model DNA polymer through a solid-state nanopore in the presence of an extended electric field using the Brownian dynamics simulation that enables us to record statistics of the conformations at every stage of the translocation process. By releasing the equilibrated DNAs from different equipotentials, we observe that the capture time distribution depends on the initial starting point and follows a Poisson process. The field gradient elongates the DNA on its way toward the nanopore and favors a successful translocation even after multiple failed threading attempts. Even in the limit of an extremely narrow pore, a fully flexible chain has a finite probability of hairpin-loop capture, while this probability decreases for a stiffer chain and promotes single file translocation. Our in silico studies identify and differentiate characteristic distributions of the mean first passage time due to single file translocation from those due to translocation of different types of folds and provide direct evidence of the interpretation of the experimentally observed folds [M. Gershow and J. A. Golovchenko, Nat. Nanotechnol. 2, 775 (2007) and Mihovilovic et al., Phys. Rev. Lett. 110, 028102 (2013)] in a solitary nanopore. Finally, we show a new finding—that a charged tag attached at the 5′ end of the DNA enhances both the multi-scan rate and the uni-directional translocation (5′ → 3′) probability that would benefit the genomic barcoding and sequencing experiments.

Nonequilibrium free energy during polymer chain growth

During fast diffusion-influenced polymerization, nonequilibrium behavior of the polymer chains and the surrounding reactive monomers has been reported recently. Based on the laws of thermodynamics, the emerging nonequilibrium structures should be characterizable by some “extra free energy” (excess over the equilibrium Helmholtz free energy). Here, we study the nonequilibrium thermodynamics of chain-growth polymerization of ideal chains in a dispersion of free reactive monomers, using off-lattice, reactive Brownian dynamics computer simulations in conjunction with approximative statistical mechanics and relative entropy (Gibbs–Shannon and Kullback–Leibler) concepts. In the case of fast growing polymers, we indeed report increased nonequilibrium free energies Δ Fneq of several kB T compared to equilibrium and near-equilibrium, slowly growing chains. Interestingly, Δ Fneq is a non-monotonic function of the degree of polymerization and thus also of time. Our decomposition of the thermodynamic contributions shows that the initial dominant extra free energy is stored in the nonequilibrium inhomogeneous density profiles of the free monomer gas (showing density depletion and wakes) in the vicinity of the active center at the propagating polymer end. At later stages of the polymerization process, we report significant extra contributions stored in the nonequilibrium polymer conformations. Finally, our study implies a nontrivial relaxation kinetics and “restoring” of the extra free energy during the equilibration process after polymerization.

Quantifying and controlling bond multivalency for advanced nanoparticle targeting to cells

Nanoparticles have drawn intense interest as delivery agents for diagnosing and treating various cancers. Much of the early success was driven by passive targeting mechanisms such as the enhanced permeability and retention (EPR) effect, but this has failed to lead to the expected clinical successes. Active targeting involves binding interactions between the nanoparticle and cancer cells, which promotes tumor cell-specific accumulation and internalization. Furthermore, nanoparticles are large enough to facilitate multiple bond formation, which can improve adhesive properties substantially in comparison to the single bond case. While multivalent binding is universally believed to be an attribute of nanoparticles, it is a complex process that is still poorly understood and difficult to control. In this review, we will first discuss experimental studies that have elucidated roles for parameters such as nanoparticle size and shape, targeting ligand and target receptor densities, and monovalent binding kinetics on multivalent nanoparticle adhesion efficiency and cellular internalization. Although such experimental studies are very insightful, information is limited and confounded by numerous differences across experimental systems. Thus, we focus the second part of the review on theoretical aspects of binding, including kinetics, biomechanics, and transport physics. Finally, we discuss various computational and simulation studies of nanoparticle adhesion, including advanced treatments that compare directly to experimental results. Future work will ideally continue to combine experimental data and advanced computational studies to extend our knowledge of multivalent adhesion, as well as design the most powerful nanoparticle-based agents to treat cancer.

Controlling solvent quality by time: Self-avoiding sprints in nonequilibrium polymerization

A fundamental paradigm in polymer physics is that macromolecular conformations in equilibrium can be described by universal scaling laws, being key for structure, dynamics, and function of soft (biological) matter and in the materials sciences. Here we reveal that during diffusion-influenced, nonequilibrium chain-growth polymerization, scaling laws change qualitatively, in particular, the growing polymers exhibit a surprising self-avoiding walk behavior in poor and θ solvents. Our analysis, based on monomer-resolved, off-lattice reaction-diffusion computer simulations, demonstrates that this phenomenon is a result of (i) nonequilibrium monomer density depletion correlations around the active polymerization site, leading to a locally directed and self-avoiding growth, in conjunction with (ii) chain (Rouse) relaxation times larger than the competing polymerization reaction time. These intrinsic nonequilibrium mechanisms are facilitated by fast and persistent reaction-driven diffusion ("sprints") of the active site, with analogies to pseudochemotactic active Brownian particles. Our findings have implications for time-controlled structure formation in polymer processing, as in, e.g., reactive self-assembly, photocrosslinking, and three-dimensional printing.

Characteristic time for the end monomers of a spherically confined polymer to find a nano-pore

Translocation of a polymer through a nano-pore is relevant in a variety of contexts such as passage of RNAs through a nuclear pore and transportation of proteins across a membrane. An essential step in polymer translocation is for the end monomers to search the pore. This process requires a characteristic time, referred to as the "attempt time" in this work. Here, we study the attempt time τ of a confined polymer inside a spherical surface by combining a scaling approach and Langevin dynamics simulations. For a moderately to strongly confined polymer, our results suggest that τ ∼ R^3.67 for R > P and τ ∼ R^2.67 for R < P, where R is the radius of the spherical surface and P is the persistence length of the polymer. All simulation data obtained for an intermediate range of the volume fraction of monomers ϕ(≲ 0.2) tend to collapse onto each other. This implies that τ does not explicitly depend on ϕ, in agreement with the theoretical predictions. These results will be useful for interpreting translocation as a two-step process: the initial attempt to find the pore and eventual pore crossing.

The translocation dynamics of the polymer through a conical pore: Non-stuck, weak-stuck, and strong-stuck modes

The external voltage-driven polymer translocation through a conical pore (with a large opening at the entry and a small tip at the exit) is studied by using the Langevin dynamics simulation in this paper. The entire translocation process is divided into an approaching stage and a threading stage. First, the approaching stage starts from the polymer entering the large opening and ends up at a terminal monomer reaching the pore tip. In this stage, the polymer will undergo the conformation adjustment to fit the narrowed cross-sectional area of the pore, leading to three approaching modes: the non-stuck mode with a terminal monomer arriving at the pore tip smoothly, the weak-stuck mode for the polymer stuck inside the pore for a short duration with minor conformational adjustments, and the strong-stuck mode with major conformational changes and a long duration. The approaching times (the duration of the approaching stage) of the three approaching modes show different behavior as a function of the pore apex angle. Second, the threading stage describes that the polymer threads through the pore tip with a linear fashion. In this stage, an increase in the apex angle causes the reduction of the threading time (the duration of the threading stage) due to the increase in the driving force with the apex angle at the tip. Moreover, we also find that with the increase in the apex angle or the polymer length, the polymer threading dynamics will change from the quasi-equilibrium state to the non-equilibrium state.

An efficient kinetic Monte Carlo to study analyte capture by a nanopore: transients, boundary conditions and time-dependent fields

To better understand the capture process by a nanopore, we introduce an efficient Kinetic Monte Carlo (KMC) algorithm that can simulate long times and large system sizes by mapping the dynamic of a point-like particle in a 3D spherically symmetric system onto the 1D biased random walk. Our algorithm recovers the steady-state analytical solution and allows us to study time-dependent processes such as transients. Simulation results show that the steady-state depletion zone near pore is barely larger than the pore radius and narrows at higher field intensities; as a result, the time to reach steady-state is much smaller than the time required to empty a zone of the size of the capture radius λe. When the sample reservoir has a finite size, a second depletion region propagates inward from the outer wall, and the capture rate starts decreasing when it reaches the capture radius λe. We also note that the flatness of the electric field near the pore, which is often neglected, induces a traffic jam that can increase the transient time by several orders of magnitude. Finally, we propose a new proof-of-concept scheme to separate two analytes of the same mobility but different diffusion coefficients using time-varying fields.

Pore shapes effects on polymer translocation

We translocated polymers through pores of different shapes and interaction patterns in three dimensions by Langevin molecular dynamics. There were four simple cylindrical pores of the same length but with different diameters. The results showed that even though decreasing the pore diameter would always decrease the translocation velocity, it was strongly dependent on the shape of the increased pore diameter. Although increasing the pore diameter made the translocation faster in simple cylindrical pores, it was complicated in different pore shapes, e.g. increasing the diameter in the middle decreased the translocation velocity. Investigating polymer shapes through the translocation process and comparing the shapes by the cumulative waiting time for different pore structures reveals the non-equilibrium properties of translocation. Moreover, polymer shape parameters such as gyration radius, polymer center of mass, and average aspect ratio help us to distinguish different pore shapes and/or different polymers.

Translocation through a narrow pore under a pulling force

We employ a three-dimensional molecular dynamics to simulate a driven polymer translocation through a nanopore by applying an external force, for four pore diameters and two external forces. To see the polymer and pore interaction effects on translocation time, we studied nine interaction energies. Moreover, to better understand the simulation results, we investigate polymer center of mass, shape factor and the monomer spatial distribution through the translocation process. Our results reveal that increasing the polymer-pore interaction energy is accompanied by an increase in the translocation time and decrease in the process rate. Furthermore, for pores with greater diameter, the translocation becomes faster. The shape analysis of the polymer indicates that the polymer shape is highly sensitive to the interaction energy. In great interactions, the monomers come close to the pore from both sides. As a result, the translocation becomes fast at first and slows down at last. Overall, it can be concluded that the external force does not play a major role in the shape and distribution of translocated monomers. However, the interaction energy between monomer and nanopore has a major effect especially on the distribution of translocated monomers on the trans side.

Polymer translocation into cavities: Effects of confinement geometry, crowding, and bending rigidity on the free energy

Monte Carlo simulations are used to study the translocation of a polymer into a cavity. Modeling the polymer as a hard-sphere chain with a length up to N=601 monomers, we use a multiple-histogram method to measure the variation of the conformational free energy of the polymer with respect to the number of translocated monomers. The resulting free-energy functions are then used to obtain the confinement free energy for the translocated portion of the polymer. We characterize the confinement free energy for a flexible polymer in cavities with constant cross-sectional area A for various cavity shapes (cylindrical, rectangular, and triangular) as well as for tapered cavities with pyramidal and conical shape. The scaling of the free energy with cavity volume and translocated polymer subchain length is generally consistent with predictions from simple scaling arguments, with small deviations in the scaling exponents likely due to finite-size effects. The confinement free energy depends strongly on cavity shape anisometry and is a minimum for an isometric cavity shape with a length-to-width ratio of unity. Entropic depletion at the edges or vertices of the confining cavity are evident in the results for constant-A and pyramidal cavities. For translocation into infinitely long cones, the scaling of the free energy with taper angle is consistent with a theoretical prediction employing the blob model. We also examine the effects of polymer bending rigidity on the translocation free energy for cylindrical cavities. For isometric cavities, the observed scaling behavior is in partial agreement with theoretical predictions, with discrepancies arising from finite-size effects that prevent the emergence of well-defined scaling regimes. In addition, translocation into highly anisometric cylindrical cavities leads to a multistage folding process for stiff polymers. Finally, we examine the effects of crowding agents inside the cavity. We find that the confinement free energy increases with crowder density. At constant packing fraction the magnitude of this effect lessens with increasing crowder size for a crowder-to-monomer size ratio ≥1.

The breakdown of the local thermal equilibrium approximation for a polymer chain during packaging

The conformational relaxation of a polymer chain often slows down in various biological and engineering processes. The polymer, then, may stay in nonequilibrium states throughout the process such that one may not invoke the local thermal equilibrium (LTE) approximation, which has been usually employed to describe the kinetics of various processes. In this work, motivated by recent single-molecule experiments on DNA packaging into a viral capsid, we investigate how the nonequilibrium conformations and the LTE approximation would affect the packaging of a polymer chain into small confinement. We employ a simple but generic coarse-grained model and Langevin dynamics simulations to investigate the packaging kinetics. The polymer segments (both inside and outside the confinement) stay away from equilibrium under strong external force. We devise a simulation scheme to invoke the LTE approximation during packaging and find that the relaxation of nonequilibrium conformations plays a critical role in regulating the packaging rate.

Trapped and non-trapped polymer translocations through a spherical pore

The polymer translocation through a spherical pore is studied using the Langevin dynamics simulation. The translocation events are classified into two types: one is the trapped translocation in which the entire polymer is trapped in the pore and the other is the non-trapped translocation where the pore cannot hold the whole polymer. We find that the trapped translocation is favored at large spheres and small external voltages. However, the monomer-pore attraction would lead to the non-monotonic behavior of the trapped translocation possibility out of all translocation events. Moreover, both the trapped and non-trapped translocation times are dependent on the polymer length, pore size, external voltage, and the monomer-pore attraction. There exist two pathways for the polymer in the trapped translocation: an actively trapped pathway for the polymer trapped in the pore before the head monomer arrives at the pore exit, and a passively trapped pathway for the polymer trapped in the pore while the head monomer is struggling to move out of the pore. The studies of trapped pathways can provide a deep understanding of the polymer translocation behavior.

Effects of solvent quality and non-equilibrium conformations on polymer translocation

The conformation and its relaxation of a single polymer depend on solvent quality in a polymer solution: a polymer collapses into a globule in a poor solvent, while the polymer swells in a good solvent. When one translocates a polymer through a narrow pore, a drastic conformational change occurs such that the kinetics of the translocation is expected to depend on the solvent quality. However, the effects of solvent quality on the translocation kinetics have been controversial. In this study, we employ a coarse-grained model for a polymer and perform Langevin dynamics simulations for the driven translocation of a polymer in various types of solvents. We estimate the free energy of polymer translocation using steered molecular dynamics simulations and Jarzynski's equality and find that the free energy barrier for the translocation increases as the solvent quality becomes poorer. The conformational entropy contributes most to the free energy barrier of the translocation in a good solvent, while a balance between entropy and energy matters in a poor solvent. Interestingly, contrary to what is expected from the free energy profile, the translocation kinetics is a non-monotonic function of the solvent quality. We find that for any type of solvent, the polymer conformation stays far away from the equilibrium conformation during translocation due to an external force and tension propagation. However, the degree of tension propagation differs depending on the solvent quality as well as the magnitude of the external force: the tension propagation is more significant in a good solvent than in a poor solvent. We illustrate that such differences in tension propagation and non-equilibrium conformations between good and poor solvents are responsible for the complicated non-monotonic effects of solvent quality on the translocation kinetics.

Simulation study on the translocation of polyelectrolyte through conical nanopores

Experiments have suggested that the conical nanopore may be a reasonable sensor for the biopolymer analysis as it can provide high-resolution current signal. In this paper, we use Langevin dynamics simulation to study the translocation of charged polymer (polyelectrolyte) through three different conical nanopores, a single-conical nanopore with large entry and small exit (pore I), a single-conical nanopore with small entry and large exit (pore II), and a double-conical nanopore with the tip (narrowest place) at the middle (pore III). Simulation shows that the detailed translocation behaviors are of pore structure dependence. Pore I might be the most reasonable one for the polyelectrolyte analysis, especially at strong monomer-pore attraction, since it can efficiently reduce the polyelectrolyte speed at the tip. The simulation results are interpreted by the free energy profiles of the polyelectrolyte translocation through different pores and the time of individual monomer passing through the tips.

Translocation of Charged Polymers through a Nanopore in Monovalent and Divalent Salt Solutions: A Scaling Study Exploring over the Entire Driving Force Regimes

Langevin dynamics simulations are performed to study polyelectrolytes driven through a nanopore in monovalent and divalent salt solutions. The driving electric field E is applied inside the pore, and the strength is varied to cover the four characteristic force regimes depicted by a rederived scaling theory, namely the unbiased (UB) regime, the weakly-driven (WD) regime, the strongly-driven trumpet (SD(T)) regime and the strongly-driven isoflux (SD(I)) regime. By changing the chain length N, the mean translocation time is studied under the scaling form 〈 τ 〉 ∼ N α E - δ . The exponents α and δ are calculated in each force regime for the two studied salt cases. Both of them are found to vary with E and N and, hence, are not universal in the parameter's space. We further investigate the diffusion behavior of translocation. The subdiffusion exponent γ p is extracted. The three essential exponents ν s , q, z p are then obtained from the simulations. Together with γ p , the validness of the scaling theory is verified. Through a comparison with experiments, the location of a usual experimental condition on the scaling plot is pinpointed.

Theoretical study on the polymer translocation into an attractive sphere

We report a non-sampling model, combining the blob method with the standard lattice-based approximation, to calculate the free energy for the polymer translocation into an attractive sphere (i.e., spherical confined trans side) through a small pore. The translocation time is then calculated by the Fokker-Planck equation based on the free energy profile. There is a competition between the confinement effect of the sphere and the polymer-sphere attraction. The translocation time is increased due to the confinement effect of the sphere, whereas it is reduced by the polymer-sphere attraction. The two effects offset each other at a special polymer-sphere attraction which is dependent on the sphere size, the polymer length, and the driving force. Moreover, the entire translocation process can be divided into an uncrowded stage where the polymer does not experience the confinement effect of the sphere and a crowded stage where the polymer is confined by the sphere. At the critical sphere radius, the durations of the two (uncrowded and crowded) stages are the same. The critical sphere radius R^* has a scaling relation with the polymer length N as R^* ∼ N^β. The calculation results show that the current model can effectively treat the translocation of a three-dimensional self-avoiding polymer into the spherical confined trans side.

Dynamics of a polyelectrolyte through aerolysin channel as a function of applied voltage and concentration⋆

We describe the behaviour of a polyelectrolyte in confined geometry. The transport of a polyelectrolyte, dextran sulfate, through a recombinant protein channel, aerolysin, inserted into a planar lipid bilayer is studied as a function of applied voltage and polyelectrolyte concentration and chain length. The aerolysin pore has a weak geometry asymmetry, a high number of charged residues and the polyelectrolyte is strongly negatively charged. The resulting current blockades were characterized by short and long dwelling times. Their frequency varies exponentially as a function of applied voltage and linearly as a function of polyelectrolyte concentration. The long blockade duration decreases exponentially when the electrical force increases. The ratio of the population of short events to the one of long events decreases when the applied voltage increases and displays an exponential variation. The long residence time increases with the polyelectrolyte chain length. We measure a reduction of the effective charge of the polyelectrolyte at the pore entry and inside the channel. For a fixed applied voltage, + / - 100 mV, at both sides of the protein pore entrance, the events frequency is similar as a function of dextran sulfate concentration. The mean blockade durations are independent of polyelectrolyte concentration and are similar for both entrances of the pore and remain constant as a function of the electrical force.