Unifying model of driven polymer translocation

Na+-Mg2+ ion effects on conformation and translocation dynamics of single-stranded RNA: Cooperation and competition

The nanopore-based translocation of a single-stranded RNA (ssRNA) in mixed salt solution has garnered increasing interest for its biological and technological significance. However, it is challenging to comprehensively understand the effects of the mixed ion species on the translocation dynamics due to their cooperation and competition, which can be directly reflected by the ion screening and neutralizing effects, respectively. In this study, Langevin dynamics simulation is employed to investigate the properties of ssRNA conformation and translocation in mixed Na+-Mg2+ ion environments. Simulation results reveal that the ion screening effect dominates the change in the ssRNA conformational size, the ion neutralizing effect controls the capture rate of the ssRNA by the nanopore, and both of them take charge of the different changes in translocation time of the ssRNA under various mixed ion environments. Under high Na+ ion concentration, as Mg2+ concentration increases, the ion neutralizing effect strengthens, weakening the driving force inside the nanopore, leading to longer translocation time. Conversely, at low Na+ concentration, an increase in Mg2+ concentration enhances the ion screening effect, aiding in faster translocation. Furthermore, these simulation results will be explained by quantitative analysis, advancing a deeper understanding of the complicated effects of the mixed Na+-Mg2+ ions.

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.

Translocation Dynamics of an Active Filament through a Long-Length Scale Channel

Active filament translocation through a confined space is crucial for diverse biological processes. By using Langevin dynamics simulations, we investigate the translocation dynamics of an axially self-propelled chain through a channel. First, results show a suggestive reciprocal scaling of translocation time versus active force. Second, in the case of a long channel, we demonstrate a very intriguing nonmonotonic change of translocation time with increasing channel width. The driving force shows a similar trend, providing a consistent picture to understand the unexpected channel width effect. In particular, in a moderately broad channel, the disordered chain conformation results in a loss of driving force and thus inhibits translocation dynamics. Chain adsorption might occur in a wide channel, which accounts for a facilitated translocation. Lastly, we connect the translocation process to tension propagation (TP). A modified TP picture is proposed to interpret the waiting time distribution. Our work highlights the new phenomenology owing to the crucial interplay of activity and spacial confinement, which drives the translocation dynamics, going beyond the traditional entropic barrier scenario.

Ratcheting Charged Polymers through Symmetric Nanopores Using Pulsed Fields: Designing a Low Pass Filter for Concentrating Polyelectrolytes

We present a new concept for the separation of DNA molecules by contour length that combines a nanofluidic ratchet, nanopore translocation, and pulsed fields. Using Langevin dynamics simulations, we show that it is possible to design pulsed field sequences to ratchet captured semiflexible molecules in such a way that only short chains successfully translocate, effectively transforming the nanopore process into a low pass molecular filter. We also show that asymmetric pulses can significantly enhance the device efficiency. The process itself can be performed with many pores in parallel, and it should be possible to integrate it directly into nanopore sequencing devices, increasing its potential utility.

Nanopore translocation of topologically linked DNA catenanes

The electrical signal associated with a biopolymer translocating through a nanoscale pore depends on the size, topology, and configuration of each molecule. Building upon recent interest in using solid-state nanopores for studying the topology of knotted and supercoiled DNA, we present experimental observations of topologically linked catenanes translocating through a solid-state nanopore. Using restriction enzymes, linked circular molecules were isolated from the mitochondrial DNA of Crithidia fasciculata, a structure known as a kinetoplast that comprises thousands of topologically interlocked minicircles. Digested kinetoplasts produce a spectrum of catenane topologies, which are identified from their nanopore translocation signals by spikes in the blockade current associated with the topological linkages. We attribute the different patterns of the measured electrical signals to 2-catenanes, linear and triangular 3-catenanes, and several types of 4- and 5-catenanes as well as more complex structures. Measurements of the translocation time of signals consistent with 2- and 3-catenanes suggest that topological friction between the linkages and the pore slows the translocation time of these structures, as predicted in recent simulations.

Nonequilibrium Thermodynamics of DNA Nanopore Unzipping

Using theory and simulations, we carried out a first systematic characterization of DNA unzipping via nanopore translocation. Starting from partially unzipped states, we found three dynamical regimes depending on the applied force f: (i) heterogeneous DNA retraction and rezipping (f<17  pN), (ii) normal (17  pN<f<60  pN), and (iii) anomalous (f>60  pN) drift-diffusive behavior. We show that the normal drift-diffusion regime can be effectively modeled as a one-dimensional stochastic process in a tilted periodic potential. We use the theory of stochastic processes to recover the potential from nonequilibrium unzipping trajectories and show that it corresponds to the free-energy landscape for single-base-pair unzipping. Applying this general approach to other single-molecule systems with periodic potentials ought to yield detailed free-energy landscapes from out-of-equilibrium trajectories.

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.

Dielectric Manipulation of Polymer Translocation Dynamics in Engineered Membrane Nanopores

The alteration of the dielectric membrane properties by membrane engineering techniques such as carbon nanotube (CNT) coating opens the way to novel molecular transport strategies for biosensing purposes. In this article, we predict a macromolecular transport mechanism enabling the dielectric manipulation of the polymer translocation dynamics in dielectric membrane pores confining mixed electrolytes. In the giant permittivity regime of these engineered membranes governed by attractive polarization forces, multivalent ions adsorbed by the membrane nanopore trigger a monovalent ion separation and set an electroosmotic counterion flow. The drag force exerted by this flow is sufficiently strong to suppress and invert the electrophoretic velocity of anionic polymers and also to generate the mobility of neutral polymers whose speed and direction can be solely adjusted by the charge and concentration of the added multivalent ions. These features identify the dielectrically generated transport mechanism as an efficient means to drive overall neutral or weakly charged analytes that cannot be controlled by an external voltage. We also reveal that, in anionic polymer translocation, multivalent cation addition into the monovalent salt solution amplifies the electric current signal by several factors. The signal amplification is caused by the electrostatic many-body interactions replacing the monovalent polymer counterions by the multivalent cations of higher electric mobility. The strength of this electrokinetic charge discrimination points out the potential of multivalent ions as current amplifiers capable of providing boosted resolution in nanopore-based biosensing techniques.

Polymer Translocation Time

The force- and flow-induced translocation processes of linear and ring polymers are studied using a combination of multiparticle collision dynamics and molecular dynamics, focusing on the behavior of the polymer translocation time. We compare the force- and flow-induced translocations of linear and ring polymers. It is found that when the translocation time (τ*) is characterized by scaling exponents, δ, δ', and α, via the relations τ* ∼ f^δN^α and τ* ∼ J^δ'N^α, the scaling exponents are not constants. For long chains tested, α = 1.0 for both force- and flow-induced translocations. The difference between the force- and flow-induced translocations stems from different monomer crowding effects due to distinct flow patterns outside the channel. Furthermore, general relations for polymer translocation time are derived for these two translocation processes, which are in good agreement with the simulation results. Our results provide clear molecular pictures for the force- and flow-induced translocations, which shed light on the understanding of translocation dynamics and provide guidance for practical applications such as molecular sequencing and ultrafiltration.

Simulation of nano elastic polymer chain displacement under pressure gradient/electroosmotic flow with the target of less dispersion of transition

Since non-scattering transfer of polymer chain in nanochannel is one of the important issue in biology, in this research, the behavior study of a long polymer chain in the nanofluid in two modes of free motion and restricted motion (fixed two ends) under two different forces including constant force (pressure gradient (PG)) and variable force (electroosmotic force (EOF)) has been investigated using dissipative particle dynamics (DPD) method. Our aim is that displacement of polymer chain carries out with less dispersion. Initially, without the presence of polymer, the results have been validated in a nanochannel by analytical results for both cases (PG, EOF) with an error of less than 10%. Then, assuming 50 beads of polymer chain, the polymer chain motion in free motion and fixed two ends modes has been examined by different spring coefficients between beads and different forces including PG (0.01 DPD unite) and EOF (zeta potential =  − 25 mV, electric field = 250 V/mm, kh parameter = 8). The results show that in free polymer motion-PG mode, by increasing 1.6 times of spring coefficient of the polymer, a 40% reduction in transition of polymer is achieved, which high dispersion of polymer chain is resulted for this mode. In the EOF, the spring coefficient has a slight effect on transferring of polymer and also, EOF moves the polymer chain with extremely low polymer chain scattering. Also, for fixed two ends-PG mode, a 36% reduction in displacement is achieved and in the same way, in EOF almost 39% declining in displacement is resulted by enhancing the spring coefficients. The results have developed to 25 and 100 beads which less dispersion of polymer chain transfer for free polymer chain-EOF is reported again for both circumstances and for restricted polymer chain state in two PG and EOF modes, less differences are reported for two cases. The results show that the EOF has the benefit of low dispersion for free polymer chain transfer, also, almost equal displacement for restricted polymer chain mode is observed for both cases.

Unfolding and Translocation of Knotted Proteins by Clp Biological Nanomachines: Synergistic Contribution of Primary Sequence and Topology Revealed by Molecular Dynamics Simulations

We use Langevin dynamics simulations to model, at an atomistic resolution, how various natively knotted proteins are unfolded in repeated allosteric translocating cycles of the ClpY ATPase. We consider proteins representative of different topologies, from the simplest knot (trefoil 3₁), to the three-twist 5₂ knot, to the most complex stevedore, 6₁, knot. We harness the atomistic detail of the simulations to address aspects that have so far remained largely unexplored, such as sequence-dependent effects on the ruggedness of the landscape traversed during knot sliding. Our simulations reveal the combined effect on translocation of the knotted protein structure, i.e., backbone topology and geometry, and primary sequence, i.e., side chain size and interactions, and show that the latter can dominate translocation hindrance. In addition, we observe that due to the interplay between the knotted topology and intramolecular contacts the transmission of tension along the polypeptide chain occurs very differently from that of homopolymers. Finally, by considering native and non-native interactions, we examine how the disruption or formation of such contacts can affect the translocation processivity and concomitantly create multiple unfolding pathways with very different activation barriers.

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.

RNA Pore Translocation with Static and Periodic Forces: Effect of Secondary and Tertiary Elements on Process Activation and Duration

We use MD simulations to study the pore translocation properties of a pseudoknotted viral RNA. We consider the 71-nucleotide-long xrRNA from the Zika virus and establish how it responds when driven through a narrow pore by static or periodic forces applied to either of the two termini. Unlike the case of fluctuating homopolymers, the onset of translocation is significantly delayed with respect to the application of static driving forces. Because of the peculiar xrRNA architecture, activation times can differ by orders of magnitude at the two ends. Instead, translocation duration is much smaller than activation times and occurs on time scales comparable at the two ends. Periodic forces amplify significantly the differences at the two ends, for both activation times and translocation duration. Finally, we use a waiting-times analysis to examine the systematic slowing downs in xrRNA translocations and associate them to the hindrance of specific secondary and tertiary elements of xrRNA. The findings provide a useful reference to interpret and design future theoretical and experimental studies of RNA translocation.

Pulling a folded polymer through a nanopore

We investigate the translocation dynamics of a folded linear polymer which is pulled through a nanopore by an external force. To this end, we generalize the iso-flux tension propagation theory for end-pulled polymer translocation to include the case of two segments of the folded polymer traversing simultaneously trough the pore. Our theory is extensively benchmarked with corresponding molecular dynamics (MD) simulations. The translocation process for a folded polymer can be divided into two main stages. In the first stage, both branches are traversing the pore and their dynamics is coupled. If the branches are not of equal length, there is a second stage where translocation of the shorter branch has been completed. Using the assumption of equal monomer flux of both branches confirmed by MD simulations, we analytically derive the equations of motion for both branches and characterize the translocation dynamics in detail from the average waiting time and its scaling form.

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.

The journey of a single polymer chain to a nanopore

For a polymer to successfully thread through a nanopore, it must first find the nanopore. This so-called capture process is typically considered as a two-stage operation consisting of the chain being delivered at the entrance of the nanopore and then insertion of one of the ends. Studying molecular dynamics-lattice Boltzmann simulations of the capture of a single polymer chain under pressure driven hydrodynamic flow, we observe that the insertion can be essentially automatic with no delay for the ends searching for the nanopore. The deformation of the chain within the converging flow area and also, the interplay between the chain elastic forces and the hydrodynamic drag play an important role in the capture of the chain by the nanopore. Along the journey to the nanopore, the chain may form folded shapes. The competition between the elastic and hydrodynamic forces results in unraveling of the folded conformations (hairpins) as the chain approaches the nanopore. Although the ends are not the only monomers that can thread into the nanopore, the unraveling process can result in much higher probability of threading by the ends.

Pulling a DNA molecule through a nanopore embedded in an anionic membrane: tension propagation coupled to electrostatics

We consider the influence of electrostatic forces on driven translocation dynamics of a flexible polyelectrolyte being pulled through a nanopore by an external force on the head monomer. To this end, we augment the iso-flux tension propagation theory with electrostatics for a negatively charged biopolymer pulled through a nanopore embedded in a similarly charged anionic membrane. We show that in the realistic case of a single-stranded DNA molecule, dilute salt conditions characterized by weak charge screening, and a negatively charged membrane, the translocation dynamics is unexpectedly accelerated despite the presence of large repulsive electrostatic interactions between the polymer coil on thecisside and the charged membrane. This is due to the rapid release of the electrostatic potential energy of the coil during translocation, leading to an effectively attractive force that assists end-driven translocation. The speedup results in non-monotonic polymer length and membrane charge dependence of the exponentαcharacterizing the translocation timeτ∝N0αof the polymer with lengthN0. In the regime of long polymersN0 ≳   500, the translocation exponent exceeds its upper limitα= 2 previously observed for the same system without electrostatic interactions.

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.

Investigation of Magneto Hydro-Dynamics Effects on a Polymer Chain Transfer in Micro-Channel Using Dissipative Particle Dynamics Method

In this paper, the effect of Magneto Hydro-Dynamics (MHD) on a polymer chain in the micro channel is studied by employing the Dissipative Particle Dynamics simulation (DPD) method. First, in a simple symmetric micro-channel, the results are evaluated and validated for different values of Hartmann (Ha) Number. The difference between the simulation and analytical solution is below 10%. Then, two types of polymer chain including short and long polymer chain are examined in the channel and the effective parameters such as Ha number, the harmony bond coefficient or spring constant (K), and the length of the polymer chain (N) are studied in the MHD flow. It is shown that by increasing harmony bond constant to 10 times with Ha = 20, the reduction of about 80% in radius of gyration squared, and half in polymer length compared to Ha = 1 would occur for both test cases. For short and long length of polymer, proper transfer of a polymer chain through MHD particles flow is observed with less perturbations (80%) and faster polymer transfer in the symmetric micro-channel.

Theoretical Modeling of Polymer Translocation: From the Electrohydrodynamics of Short Polymers to the Fluctuating Long Polymers

The theoretical formulation of driven polymer translocation through nanopores is complicated by the combination of the pore electrohydrodynamics and the nonequilibrium polymer dynamics originating from the conformational polymer fluctuations. In this review, we discuss the modeling of polymer translocation in the distinct regimes of short and long polymers where these two effects decouple. For the case of short polymers where polymer fluctuations are negligible, we present a stiff polymer model including the details of the electrohydrodynamic forces on the translocating molecule. We first show that the electrohydrodynamic theory can accurately characterize the hydrostatic pressure dependence of the polymer translocation velocity and time in pressure-voltage-driven polymer trapping experiments. Then, we discuss the electrostatic correlation mechanisms responsible for the experimentally observed DNA mobility inversion by added multivalent cations in solid-state pores, and the rapid growth of polymer capture rates by added monovalent salt in α -Hemolysin pores. In the opposite regime of long polymers where polymer fluctuations prevail, we review the iso-flux tension propagation (IFTP) theory, which can characterize the translocation dynamics at the level of single segments. The IFTP theory is valid for a variety of polymer translocation and pulling scenarios. We discuss the predictions of the theory for fully flexible and rodlike pore-driven and end-pulled translocation scenarios, where exact analytic results can be derived for the scaling of the translocation time with chain length and driving force.

Dielectric Trapping of Biopolymers Translocating through Insulating Membranes

Sensitive sequencing of biopolymers by nanopore-based translocation techniques requires an extension of the time spent by the molecule in the pore. We develop an electrostatic theory of polymer translocation to show that the translocation time can be extended via the dielectric trapping of the polymer. In dilute salt conditions, the dielectric contrast between the low permittivity membrane and large permittivity solvent gives rise to attractive interactions between the cis and trans portions of the polymer. This self-attraction acts as a dielectric trap that can enhance the translocation time by orders of magnitude. We also find that electrostatic interactions result in the piecewise scaling of the translocation time τ with the polymer length L. In the short polymer regime L≲10 nm where the external drift force dominates electrostatic polymer interactions, the translocation is characterized by the drift behavior τ∼L2. In the intermediate length regime 10nm≲L≲κb−1 where κb is the Debye–Hückel screening parameter, the dielectric trap takes over the drift force. As a result, increasing polymer length leads to quasi-exponential growth of the translocation time. Finally, in the regime of long polymers L≳κb−1 where salt screening leads to the saturation of the dielectric trap, the translocation time grows linearly as τ∼L. This strong departure from the drift behavior highlights the essential role played by electrostatic interactions in polymer translocation.

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.

Reducing the variance in the translocation times by prestretching the polymer

Langevin dynamics simulations of polymer translocation are performed where the polymer is stretched via two opposing forces applied on the first and last monomer before and during translocation. In this setup, polymer translocation is achieved by imposing a bias between the two pulling forces such that there is net displacement towards the trans side. Under the influence of stretching forces, the elongated polymer ensemble contains less variations in conformations compared to an unstretched ensemble. Simulations demonstrate that this reduced spread in initial conformations yields a reduced variation in translocation times relative to the mean translocation time. This effect is explored for different ratios of the amplitude of thermal fluctuations to driving forces to control for the relative influence of the thermal path sampled by the polymer. Since the variance in translocation times is due to contributions coming from sampling both thermal noise and initial conformations, our simulations offer independent control over the two main sources of noise and allow us to shed light on how they both contribute to translocation dynamics. Simulation parameter space corresponding to experimentally relevant conditions is highlighted and shown to correspond to a significant decrease in the spread of translocation times, thus indicating that stretching DNA prior to translocation could assist nanopore-based sequencing and sizing applications.

Dynamics of driven translocation of semiflexible polymers

We study translocation of semiflexible polymers driven by force f_{d} inside a nanometer-scale pore using our three-dimensional Langevin dynamics model. We show that the translocation time τ increases with increasing bending rigidity κ. Similarly, the exponent β for the scaling of τ with polymer length N,τ∼N^{β}, increases with increasing κ as well as with increasing f_{d}. By comparing waiting times between semiflexible and fully flexible polymers we show that for realistic f_{d} translocation dynamics is to a large extent, but not completely, determined by the polymer's elastic length measured in number of Kuhn segments N_{Kuhn}. Unlike in driven translocation of flexible polymers, friction related to the polymer segment on the trans side has a considerable effect on the resulting dynamics. This friction is intermittently reduced by buckling of the polymer segment in the vicinity of the pore opening on the trans side. We show that in the experimentally relevant regime where viscosity is higher than in computer simulation models, the probability for this buckling increases with increasing f_{d}, giving rise to a larger contribution to the trans side friction at small f_{d}. Similarly to flexible polymers, we find significant center-of-mass diffusion of the cis side polymer segment which speeds up translocation. This effect is larger for smaller f_{d}. However, this speedup is smaller than the slowing down due to the trans side friction. At large enough N_{Kuhn}, the roles can be seen to be reversed, and the dynamics of flexible polymers can be reached. However, for example, polymers used in translocation experiments of DNA are elastically so short that the finite-length dynamics outlined here applies.

Enhanced polymer capture speed and extended translocation time in pressure-solvation traps

The efficiency of nanopore-based biosequencing techniques requires fast anionic polymer capture by like-charged pores followed by a prolonged translocation process. We show that this condition can be achieved by setting a pressure-solvation trap. Polyvalent cation addition to the KCl solution triggers the like-charge polymer-pore attraction. The attraction speeds-up the pressure-driven polymer capture but also traps the molecule at the pore exit, reducing the polymer capture time and extending the polymer escape time by several orders of magnitude. By direct comparison with translocation experiments [D. P. Hoogerheide et al., ACS Nano 8, 7384 (2014)1936-085110.1021/nn5025829], we characterize as well the electrohydrodynamics of polymers transport in pressure-voltage traps. We derive scaling laws that can accurately reproduce the pressure dependence of the experimentally measured polymer translocation velocity and time. We also find that during polymer capture, the electrostatic barrier on the translocating molecule slows down the liquid flow. This prediction identifies the streaming current measurement as a potential way to probe electrostatic polymer-pore interactions.

Granular chain escape from a pore in a wall in the presence of particles on one side: a comparison to polymer translocation

Escape of a granular chain from a pore in a wall in the presence of diffusing granular particles on one side of the wall is studied experimentally. The escape time shows power-law behavior as a function of the chain length (τ ∝ Nα). A Langevin dynamics simulation of a polymer chain in a similar geometry is also performed and similar results to those for a granular system are obtained. A simple scaling argument and an energetic argument (based on the Onsager principle) are introduced which explain our results very well. Experiments (simulations) show that by increasing the number of particles on one side of the wall from zero, the exponent α decreases from 2.6 ± 0.1 (3.1 ± 0.1) to about 2. Both scaling and the Onsager principle argument predict α = 2 at high particle concentration, in agreement with the experiments and simulations. In the absence of particles, the scaling predicts τ = N2.5 (in agreement with the experimental result for the granular chain) and the Onsager principle predictsτ = N3 ln N, supporting the simulation result for the polymer chain. Experiments, simulations, scaling, and the Onsager principle confirm an inverse relation between τ and the density of particles on one side of the wall.

Facilitated polymer capture by charge inverted electroosmotic flow in voltage-driven polymer translocation

The optimal functioning of nanopore-based biosensing tools necessitates rapid polymer capture from the ion reservoir. We identify an ionic correlation-induced transport mechanism that provides this condition without the chemical modification of the polymer or the pore surface. In the typical experimental configuration where a negatively charged silicon-based pore confines a 1 : 1 electrolyte solution, anionic polymer capture is limited by electrostatic polymer-membrane repulsion and the electroosmotic (EO) flow. Added multivalent cations suppress the electrostatic barrier and reverse the pore charge, inverting the direction of the EO flow that drags the polymer to the trans side. This inverted EO flow can be used to speed up polymer capture from the reservoir and to transport weakly or non-uniformly charged polymers that cannot be controlled by electrophoresis.

Sequencing of semiflexible polymers of varying bending rigidity using patterned pores

We study the translocation of a semiflexible polymer through extended pores with patterned stickiness, using Langevin dynamics simulations. We find that the consequence of pore patterning on the translocation time dynamics is dramatic and depends strongly on the interplay of polymer stiffness and pore-polymer interactions. For heterogeneous polymers with periodically varying stiffness along their lengths, we find that variation of the block size of the sequences and the orientation results in large variations in the translocation time distributions. We show how this fact may be utilized to develop an effective sequencing strategy. This strategy involving multiple pores with patterned surface energetics can predict heteropolymer sequences having different bending rigidity to a high degree of accuracy.

DNA Translocations through Nanopores under Nanoscale Preconfinement

To reduce unwanted variation in the passage speed of DNA through solid-state nanopores, we demonstrate nanoscale preconfinement of translocating molecules using an ultrathin nanoporous silicon nitride membrane separated from a single sensing nanopore by a nanoscale cavity. We present comprehensive experimental and simulation results demonstrating that the presence of an integrated nanofilter within nanoscale distances of the sensing pore eliminates the dependence of molecular passage time distributions on pore size, revealing a global minimum in the coefficient of variation of the passage time. These results provide experimental verification that the inter- and intramolecular passage time variation depends on the conformational entropy of each molecule prior to translocation. Furthermore, we show that the observed consistently narrower passage time distributions enables a more reliable DNA length separation independent of pore size and stability. We also demonstrate that the composite nanofilter/nanopore devices can be configured to suppress the frequency of folded translocations, ensuring single-file passage of captured DNA molecules. By greatly increasing the rate at which usable data can be collected, these unique attributes will offer significant practical advantages to many solid-state nanopore-based sensing schemes, including sequencing, genomic mapping, and barcoded target detection.

Polymer translocation under a pulling force: Scaling arguments and threshold forces

DNA translocation through nanopores is one of the most promising strategies for next-generation sequencing technologies. Most experimental and numerical works have focused on polymer translocation biased by electrophoresis, where a pulling force acts on the polymer within the nanopore. An alternative strategy, however, is emerging, which uses optical or magnetic tweezers. In this case, the pulling force is exerted directly at one end of the polymer, which strongly modifies the translocation process. In this paper, we report numerical simulations of both linear and structured (mimicking DNA) polymer models, simple enough to allow for a statistical treatment of the pore structure effects on the translocation time probability distributions. Based on extremely extended computer simulation data, we (i) propose scaling arguments for an extension of the predicted translocation times τ∼N^{2}F^{-1} over the moderate forces range and (ii) analyze the effect of pore size and polymer structuration on translocation times τ.

Quantification of tension to explain bias dependence of driven polymer translocation dynamics

Motivated by identifying the origin of the bias dependence of tension propagation, we investigate methods for measuring tension propagation quantitatively in computer simulations of driven polymer translocation. Here, the motion of flexible polymer chains through a narrow pore is simulated using Langevin dynamics. We measure tension forces, bead velocities, bead distances, and bond angles along the polymer at all stages of translocation with unprecedented precision. Measurements are done at a standard temperature used in simulations and at zero temperature to pin down the effect of fluctuations. The measured quantities were found to give qualitatively similar characteristics, but the bias dependence could be determined only using tension force. We find that in the scaling relation τ∼N^{β}f_{d}^{α} for translocation time τ, the polymer length N, and the bias force f_{d}, the increase of the exponent β with bias is caused by center-of-mass diffusion of the polymer toward the pore on the cis side. We find that this diffusion also causes the exponent α to deviate from the ideal value -1. The bias dependence of β was found to result from combination of diffusion and pore friction and so be relevant for polymers that are too short to be considered asymptotically long. The effect is relevant in experiments all of which are made using polymers whose lengths are far below the asymptotic limit. Thereby, our results also corroborate the theoretical prediction by Sakaue's theory [Polymers 8, 424 (2016)2073-436010.3390/polym8120424] that there should not be bias dependence of β for asymptotically long polymers. By excluding fluctuations we also show that monomer crowding at the pore exit cannot have a measurable effect on translocation dynamics under realistic conditions.

Role of non-equilibrium conformations on driven polymer translocation

One of the major theoretical methods in understanding polymer translocation through a nanopore is the Fokker-Planck formalism based on the assumption of quasi-equilibrium of polymer conformations. The criterion for applicability of the quasi-equilibrium approximation for polymer translocation is that the average translocation time per Kuhn segment, ⟨τ⟩/NK, is longer than the relaxation time τ0 of the polymer. Toward an understanding of conditions that would satisfy this criterion, we have performed coarse-grained three dimensional Langevin dynamics and multi-particle collision dynamics simulations. We have studied the role of initial conformations of a polyelectrolyte chain (which were artificially generated with a flow field) on the kinetics of its translocation across a nanopore under the action of an externally applied transmembrane voltage V (in the absence of the initial flow field). Stretched (out-of-equilibrium) polyelectrolyte chain conformations are deliberately and systematically generated and used as initial conformations in translocation simulations. Independent simulations are performed to study the relaxation behavior of these stretched chains, and a comparison is made between the relaxation time scale and the mean translocation time (⟨τ⟩). For such artificially stretched initial states, ⟨τ⟩/NK < τ0, demonstrating the inapplicability of the quasi-equilibrium approximation. Nevertheless, we observe a scaling of ⟨τ⟩ ∼ 1/V over the entire range of chain stretching studied, in agreement with the predictions of the Fokker-Planck model. On the other hand, for realistic situations where the initial artificially imposed flow field is absent, a comparison of experimental data reported in the literature with the theory of polyelectrolyte dynamics reveals that the Zimm relaxation time (τZimm) is shorter than the mean translocation time for several polymers including single stranded DNA (ssDNA), double stranded DNA (dsDNA), and synthetic polymers. Even when these data are rescaled assuming a constant effective velocity of translocation, it is found that for flexible (ssDNA and synthetic) polymers with NK Kuhn segments, the condition ⟨τ⟩/NK < τZimm is satisfied. We predict that for flexible polymers such as ssDNA, a crossover from quasi-equilibrium to non-equilibrium behavior would occur at NK ∼ O(1000).

Asymmetric dynamics of DNA entering and exiting a strongly confining nanopore

In nanopore sensing, changes in ionic current are used to analyse single molecules in solution. The translocation dynamics of polyelectrolytes is of particular interest given potential applications such as DNA sequencing. In this paper, we determine how the dynamics of voltage driven DNA translocation can be affected by the nanopore geometry and hence the available configurational space for the DNA. Using the inherent geometrical asymmetry of a conically shaped nanopore, we examine how DNA dynamics depends on the directionality of transport. The total translocation time of DNA when exiting the extended conical confinement is significantly larger compared to the configuration where the DNA enters the pore from the open reservoir. By using specially designed DNA molecules with positional markers, we demonstrate that the translocation velocity progressively increases as the DNA exits from confinement. We show that a hydrodynamic model can account for these observations.

Translocation of a charged polymer through confined nanoenvironments is highly dependent on their geometrical parameters. Here, the authors investigate experimentally the translocation dynamics of DNA through conical nanopores and provide a quantitative model for the translocation into and out of confinement.

Driven translocation of a semi-flexible polymer through a nanopore

We study the driven translocation of a semi-flexible polymer through a nanopore by means of a modified version of the iso-flux tension propagation theory, and extensive molecular dynamics (MD) simulations. We show that in contrast to fully flexible chains, for semi-flexible polymers with a finite persistence length [Formula: see text] the trans side friction must be explicitly taken into account to properly describe the translocation process. In addition, the scaling of the end-to-end distance R N as a function of the chain length N must be known. To this end, we first derive a semi-analytic scaling form for R N, which reproduces the limits of a rod, an ideal chain, and an excluded volume chain in the appropriate limits. We then quantitatively characterize the nature of the trans side friction based on MD simulations. Augmented with these two factors, the theory shows that there are three main regimes for the scaling of the average translocation time τ ∝ N α . In the rod [Formula: see text], Gaussian [Formula: see text] and excluded volume chain [Formula: see text] ≫ 10 6 limits, α = 2, 3/2 and 1 + ν, respectively, where ν is the Flory exponent. Our results are in good agreement with available simulations and experimental data.

Active translocation of a semiflexible polymer assisted by an ATP-based molecular motor

In this work we study the assisted translocation of a polymer across a membrane nanopore, inside which a molecular motor exerts a force fuelled by the hydrolysis of ATP molecules. In our model the motor switches to its active state for a fixed amount of time, while it waits for an ATP molecule which triggers the motor, during an exponentially distributed time lapse. The polymer is modelled as a beads-springs chain with both excluded volume and bending contributions, and moves in a stochastic three dimensional environment modelled with a Langevin dynamics at a fixed temperature. The resulting dynamics shows a Michaelis-Menten translocation velocity that depends on the chain flexibility. The scaling behavior of the mean translocation time with the polymer length for different bending values is also investigated.

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.

Translocation Time through a Nanopore with an Internal Cavity Is Minimal for Polymers of Intermediate Length

The translocation of polymers through nanopores with large internal cavities bounded by two narrow pores is studied via Langevin dynamics simulations. The total translocation time is found to be a nonmonotonic function of polymer length, reaching a minimum at intermediate length, with both shorter and longer polymers taking longer to translocate. The location of the minimum is shown to shift with the magnitude of the applied force, indicating that the pore can be dynamically tuned to favor different polymer lengths. A simple model balancing the effects of entropic trapping within the cavity against the driving force is shown to agree well with simulations. Beyond the nonmonotonicity, detailed analysis of translocation uncovers rich dynamics in which factors such as going to a high force regime and the emergence of a tail for long polymers dramatically change the behavior of the system. These results suggest that nanopores with internal cavities can be used for applications such as selective extraction of polymers by length and filtering of polymer solutions, extending the uses of nanopores within emerging nanofluidic technologies.

Polymer translocation through nanopore into active bath

Polymer translocation through nanopores into a crowded environment is of ubiquitous importance in many biological processes. Here we investigate polymer translocation through a nanopore into an active bath of self-propelled particles in two-dimensional space using Langevin dynamics simulations. Interestingly, we find that the mean translocation time τ can show a bell-shape dependence on the particle activity F_a at a fixed volume fraction ϕ, indicating that the translocation process may become slower for small activity compared to the case of the passive media, and only when the particle activity becomes large enough can the translocation process be accelerated. In addition, we also find that τ can show a minimum as a function of ϕ if the particle activity is large enough, implying that an intermediate volume fraction of active particles is most favorable for the polymer translocation. Detailed analysis reveals that such nontrivial behaviors result from the two-fold effect of active bath: one that active particles tend to accumulate near the pore, providing an extra pressure hindering the translocation, and the other that they also aggregate along the polymer chain, generating an effective pulling force accelerating the translocation. Such results demonstrate that active bath plays rather subtle roles on the polymer translocation process.

Driven polymer translocation in good and bad solvent: Effects of hydrodynamics and tension propagation

We investigate the driven polymer translocation through a nanometer-scale pore in the presence and absence of hydrodynamics both in good and bad solvent. We present our results on tension propagating along the polymer segment on the cis side that is measured for the first time using our method that works also in the presence of hydrodynamics. For simulations we use stochastic rotation dynamics, also called multiparticle collision dynamics. We find that in the good solvent the tension propagates very similarly whether hydrodynamics is included or not. Only the tensed segment is by a constant factor shorter in the presence of hydrodynamics. The shorter tensed segment and the hydrodynamic interactions contribute to a smaller friction for the translocating polymer when hydrodynamics is included, which shows as smaller waiting times and a smaller exponent in the scaling of the translocation time with the polymer length. In the bad solvent hydrodynamics has a minimal effect on polymer translocation, in contrast to the good solvent, where it speeds up translocation. We find that under bad-solvent conditions tension does not spread appreciably along the polymer. Consequently, translocation time does not scale with the polymer length. By measuring the effective friction in a setup where a polymer in free solvent is pulled by a constant force at the end, we find that hydrodynamics does speed up collective polymer motion in the bad solvent even more effectively than in the good solvent. However, hydrodynamics has a negligible effect on the motion of individual monomers within the highly correlated globular conformation on the cis side and hence on the entire driven translocation under bad-solvent conditions.

Electrostatics of polymer translocation events in electrolyte solutions

We develop an analytical theory that accounts for the image and surface charge interactions between a charged dielectric membrane and a DNA molecule translocating through the membrane. Translocation events through neutral carbon-based membranes are driven by a competition between the repulsive DNA-image-charge interactions and the attractive coupling between the DNA segments on the trans and the cis sides of the membrane. The latter effect is induced by the reduction of the coupling by the dielectric membrane. In strong salt solutions where the repulsive image-charge effects dominate the attractive trans-cis coupling, the DNA molecule encounters a translocation barrier of ≈10 kBT. In dilute electrolytes, the trans-cis coupling takes over image-charge forces and the membrane becomes a metastable attraction point that can trap translocating polymers over long time intervals. This mechanism can be used in translocation experiments in order to control DNA motion by tuning the salt concentration of the solution.

Chaperone-assisted translocation of flexible polymers in three dimensions

Polymer translocation through a nanometer-scale pore assisted by chaperones binding to the polymer is a process encountered in vivo for proteins. Studying the relevant models by computer simulations is computationally demanding. Accordingly, previous studies are either for stiff polymers in three dimensions or flexible polymers in two dimensions. Here, we study chaperone-assisted translocation of flexible polymers in three dimensions using Langevin dynamics. We show that differences in binding mechanisms, more specifically, whether a chaperone can bind to a single site or multiple sites on the polymer, lead to substantial differences in translocation dynamics in three dimensions. We show that the single-binding mode leads to dynamics that is very much like that in the constant-force driven translocation and accordingly mainly determined by tension propagation on the cis side. We obtain β≈1.26 for the exponent for the scaling of the translocation time with polymer length. This fairly low value can be explained by the additional friction due to binding particles. The multiple-site binding leads to translocation the dynamics of which is mainly determined by the trans side. For this process we obtain β≈1.36. This value can be explained by our derivation of β=4/3 for constant-bias translocation, where translocated polymer segments form a globule on the trans side. Our results pave the way for understanding and utilizing chaperone-assisted translocation where variations in microscopic details lead to rich variations in the emerging dynamics.