Polymer translocation through a nanopore. II. Excluded volume effect

Nonmonotonic polymer translocation kinetics through nanopores under changing surface-polymer interactions

Understanding the dynamics of polymers in confined environments is pivotal for diverse applications ranging from polymer upcycling to bioseparations. In this study, we develop an entropic barrier model using self-consistent field theory that considers the effect of attractive surface interactions, solvation, and confinement on polymer kinetics. In this model, we consider the translocation of a polymer from one cavity into a second cavity through a single-segment-width nanopore. We find that, for a polymer in a good solvent (i.e., excluded volume, u0 > 0), there is a nonmonotonic dependence of mean translocation time (τ) on surface interaction strength, ɛ. At low ɛ, excluded volume interactions lead to an energetic penalty and longer translocation times. As ɛ increases, the surface interactions counteract the energetic penalty imposed by excluded volume and the polymer translocates faster through the nanopore. However, as ɛ continues to increase, an adsorption transition occurs, which leads to significantly slower kinetics due to the penalty of desorption from the first cavity. The ɛ at which this adsorption transition occurs is a function of the excluded volume, with higher u0 leading to an adsorption transition at higher ɛ. Finally, we consider the effect of translocation across different size cavities. We find that the kinetics for translocation into a smaller cavity speeds up while translocation to a larger cavity slows down with increasing ɛ due to higher surface contact under stronger confinement.

A comparative study of the target search of end monomers of real and Rouse chains under spherical confinement

We study the dynamics of the end monomers of a real chain confined in a spherical cavity to search for a small target on the cavity surface using Langevin dynamics simulation. The results are compared and contrasted with those of a Rouse chain to understand the influence of excluded volume interactions on the search dynamics, as characterized by the first passage time (FPT). We analyze how the mean FPT depends on the cavity size Rb, the target size a, and the degree of confinement quantified by Rg/Rb, with Rg being the polymer radius of gyration in free space. As a basic finding, the equilibrium distribution of the end monomers of a real chain in a closed spherical cavity differs from that of a Rouse chain at a given Rg/Rb, which leads to the differences between the mean FPTs of real and Rouse chains. Fitting the survival probability S(t) by a multi-exponential form, we show that the S(t) of real chains exhibits multiple characteristic times at large Rg/Rb. Our simulation results indicate that the search dynamics of a real chain exhibit three characteristic regimes as a function of Rg/Rb, including the transition from the Markovian to non-Markovian process at Rg/Rb ≈ 0.39, along with two distinct regimes at 0.39 < Rg/Rb < 1.0 and Rg/Rb > 1.0, respectively, where S(t) exhibits a single characteristic time and multiple characteristic times.

Forced and spontaneous translocation dynamics of a semiflexible active polymer in two dimensions

Polymer translocation is a fundamental topic in non-equilibrium physics and is crucially important to many biological processes in life. In the present work, we adopt two-dimensional Langevin dynamics simulations to study the forced and spontaneous translocation dynamics of an active filament. The influence of polymer stiffness on the underlying dynamics is explicitly analyzed. For the forced translocation, the results show a robust stiffness-induced inhibition, and the translocation time exhibits a dual-exponent scaling relationship with the bending modulus. Tension propagation (TP) is also examined, where we find prominent modifications in terms of both activity and stiffness. For spontaneous translocation into a pure solvent, the translocation time is almost independent of the polymer stiffness. However, when the polymer is translocated into a porous medium, an intriguing non-monotonic alteration of translocation time with increasing chain stiffness is demonstrated. The semiflexible chain is beneficial for translocation while the rigid chain is not conducive. Stiffness regulation on the diffusion dynamics of the polymer in porous media shows a consistent scenario. The interplay of activity, stiffness, and porous crowding provides a new mechanism for understanding the non-trivial translocation dynamics of an active filament in complex environments.

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.

Constraints on Knot Insertion, Not Internal Jamming, Control Polycatenane Translocation Dynamics through Crystalline Pores

The translocation of polymers through pores and channels is an archetypal process in biology and is widely studied and exploited for applications in bio- and nanotechnology. In recent times, the translocation of polymers of various different topologies has been studied both experimentally and by computer simulation. However, in some cases, a clear understanding of the precise mechanisms that drive their translocation dynamics can be challenging to derive. Experimental methods are able to provide statistical details of polymer translocation, but computer simulations are uniquely placed to uncover a finer level of mechanistic understanding. In this work, we use high-throughput molecular simulations to reveal the importance that knot insertion rates play in controlling translocation dynamics in the small pore limit, where unexpected nonpower law behavior emerges. This work both provides new predictive understanding of polycatenane translocation and shows the importance of carefully considering the role of the definition of translocation itself.

A general theory of polymer ejection tested in a quasi two-dimensional space

A general ejection theory of polymer is developed in a two- and three-dimensional space. A polymer is confined initially in a cavity and ejects spontaneously to the outer space through a nanopore channel without the help of any external stimulus. A reflective wall boundary is set at the pore entrance to prevent the falling of the head monomer of chain into the cavity. Three stages are distinguished in a process: (1) an entering stage, in which the head monomer enters the pore to search for a way to traverse the pore channel, (2) a main ejection stage, in which the chain body is transported from the cavity to the outer space, (3) a leaving stage, in which the tail monomer passes through and leaves the pore channel. Depending on the number of the monomers remaining in the cavity, the main ejection stage can be divided into the confined and the non-confined stages. The non-confined stage can be further split into the thermal escape and the entropic pulling stages. The Onsager’s variational principle is applied to derive the kinetics equation of ejection. The escape time is calculated from the corresponding Kramers’ escape problem. Extensive molecular dynamics simulations are then performed in a quasi two-dimensional space to verify the theory. The variation of the ejection speed is carefully examined. The decreasing behavior of the number of monomers in the cavity is studied in details. The scaling properties of the spending time at each processing stage are investigated systematically by varying the chain length, the cavity diameter, and the initial volume fraction of chain. The results of simulation support firmly the predictions of the theory, cross-checked in the studies of various topics. In combining with the previous investigations in the three-dimensional space, the generalized theory is very robust, able to explain the two seemly different phenomena, polymer ejection and polymer translocation, together under the same theoretical framework in the two space dimensions.

Synchronized resistive-pulse analysis with flow visualization for single micro- and nanoscale objects driven by optical vortex in double orifice

Resistive-pulse analysis is a powerful tool for identifying micro- and nanoscale objects. For low-concentration specimens, the pulse responses are rare, and it is difficult to obtain a sufficient number of electrical waveforms to clearly characterize the targets and reduce noise. In this study, we conducted a periodic resistive-pulse analysis using an optical vortex and a double orifice, which repetitively senses a single micro- or nanoscale target particle with a diameter ranging from 700 nm to 2 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu$$\end{document}μm. The periodic motion results in the accumulation of a sufficient number of waveforms within a short period. Acquired pulses show periodic ionic-current drops associated with the translocation events through each orifice. Furthermore, a transparent fluidic device allows us to synchronously average the waveforms by the microscopic observation of the translocation events and improve the signal-to-noise ratio. By this method, we succeed in distinguishing single particle diameters. Additionally, the results of measured signals and the simultaneous high-speed observations are used to quantitatively and systematically discuss the effect of the complex fluid flow in the orifices on the amplitude of the resistive pulse. The synchronized resistive-pulse analysis by the optical vortex with the flow visualization improves the pulse-acquisition rate for a single specific particle and accuracy of the analysis, refining the micro- and nanoscale object identification.

Digital Pathology Platform for Respiratory Tract Infection Diagnosis via Multiplex Single-Particle Detections

The variability of bioparticles remains a key barrier to realizing the competent potential of nanoscale detection into a digital diagnosis of an extraneous object that causes an infectious disease. Here, we report label-free virus identification based on machine-learning classification. Single virus particles were detected using nanopores, and resistive-pulse waveforms were analyzed multilaterally using artificial intelligence. In the discrimination, over 99% accuracy for five different virus species was demonstrated. This advance is accessed through the classification of virus-derived ionic current signal patterns reflecting their intrinsic physical properties in a high-dimensional feature space. Moreover, consideration of viral similarity based on the accuracies indicates the contributing factors in the recognitions. The present findings offer the prospect of a novel surveillance system applicable to detection of multiple viruses including new strains.

All-or-none folding of a flexible polymer chain in cylindrical nanoconfinement

Geometric confinement of a polymer chain results in a loss of conformational entropy. For a chain that can fold into a compact native state via a first-order-like transition, as is the case for many small proteins, confinement typically provides an entropic stabilization of the folded state, thereby shifting the location of the transition. This allows for the possibility of confinement (entropy) driven folding. Here, we investigate such confinement effects for a flexible square-well-sphere N-mer chain (monomer diameter σ) confined within a long cylindrical pore (diameter D) or a closed cylindrical box (height H = D). We carry out Wang-Landau simulations to construct the density of states, which provides access to the complete thermodynamics of the system. For a wide pore, an entropic stabilization of the folded state is observed. However, as the pore diameter approaches the size of the folded chain (D ∼ N^1/3σ), we find a destabilization effect. For pore diameters smaller than the native ground-state, the chain folds into a different, higher energy, ground state ensemble and the T vs D phase diagram displays non-monotonic behavior as the system is forced into different ground states for different ranges of D. In this regime, isothermal reduction of the confinement dimension can induce folding, unfolding, or crystallite restructuring. For the cylindrical box, we find a monotonic stabilization effect with decreasing D. Scaling laws for the confinement free energy in the athermal limit are also investigated.

Back-Side Polymer-Coated Solid-State Nanopore Sensors

We systematically investigated the influence of polymer coating on temporal resolution of solid-state nanopores. We fabricated a Si₃N₄ nanopore integrated with a polyimide sheet partially covering the substrate surface. Upon detecting the nanoparticles dispersed in an electrolyte buffer by ionic current measurements, we observed a larger resistive pulse height along with a faster current decay at the tails under larger coverage of the polymeric layer, thereby suggesting a prominent role of the water-touching Si₃N₄ thin film as a significant capacitor serving to retard the ionic current response to the ion blockade by fast translocation of particles through the nanopores. From this, we came up with back-side polymer-coated chip designs and demonstrated improved pore sensor temporal resolution by developing a nanopore with a thick polymethyl-methacrylate layer laminated on the bottom surface. The present findings may be useful in developing integrated solid-state nanopore sensors with embedded nanochannels and nanoelectrodes.

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.

Free-energy cost of localizing a single monomer of a confined polymer

We describe a simple Monte Carlo simulation method to calculate the free-energy cost of localizing a single monomer of a polymer confined to a cavity. The localization position is chosen to be on the inside surface of the confining cavity. The method is applied to a freely jointed hard-sphere polymer chain confined to cavities of spherical and cubic geometries. In the latter case, we consider localization at a corner and at the center of a face of the confining cube. We consider cases of end-monomer localization both with and without tethering of the other end monomer to a point on the surface. We also examine localization of monomers at arbitrary positions along the contour of the polymer. We characterize the dependence of the free energy on the cavity size and shape, the localization position, and the polymer length. The quantitative trends can be understood using standard scaling arguments and use of a simple theoretical model. The results are relevant to those theories of polymer translocation that focus on the importance of the free-energy barrier as the translocation process requires an initial localization of a monomer to the position of a nanopore.

Polymer translocation through a hairy channel mimicking the inner plug of a nuclear pore complex

A microscopic transport model of a polymer translocating through a nuclear pore complex (NPC) is presented based on self-consistent field theory (SCFT), with the NPC and its nucleoporins mimicked by a hairy channel. Multiple cell environment effects (electrolyte effect, excluded volume effect, NPC drag effect, and hydrophobic effect) are all considered in this hairy channel model. The influence of various parameters (polymer chain length, length of NPC, strength of hydrophobic effect, and excluded volume effect) on translocation time is studied through theoretical analysis and numerical calculation. Numerical simulation results show that an area of low nucleoporin number density exists in the NPC, which facilitates the translocation of the polymer. The results also show that the translocation time curves with increasing NPC length and polymer charge number are concave. In addition, there are critical values for NPC length and polymer charge number for which the translocation time has a minimal value. The translocation time decreases with the increasing strength of the hydrophobic effect and excluded volume effect.

Identifying Single Particles in Air Using a 3D-Integrated Solid-State Pore

Solid-state micro- and nanopores are a versatile sensor platform capable of detecting single particles in electrolyte solution by cross-pore ionic current. Here we report on a use of this technology to identify airborne particulate matter. The detection concept lies in an electrophoretic control of air-floating particles captured in liquid to deliver them into a pore detector via microfluidic channels. We demonstrate resistive pulse measurements to machine-learning-based discriminations of intragranular contents of cypress and cedar pollens at a single-particle level. This all-electrical-sensor technique would pave a new venue toward real-time monitoring of single particles and molecules in air.

Silicon substrate effects on ionic current blockade in solid-state nanopores

We investigated the roles of silicon substrate material compositions in ionic current blockade in solid-state nanopores. When detecting single nanoparticles using an ionic current in a Si3N4 nanopore supported on a doped silicon wafer, resistive pulses were found to be blunted significantly via signal retardation due to predominant contributions of large capacitance at the ultrathin membrane. Unexpectedly, in contrast, changing the substrate material to non-doped silicon led to the sharpening of the spike-like signal feature, suggesting a better temporal resolution of the cross-channel ionic current measurements by virtue of the thick intrinsic semiconductor layer that served to diminish the net chip capacitance. The present results suggest the importance of the choice of Si compositions regarding the capacitance effects to attain better spatiotemporal resolution in solid-state nanopore sensors.

Particle Capture in Solid-State Multipores

Utilization of multiple-channel structure is a promising way of accomplishing high-throughput detections of analytes in solid-state pore sensors. Here we report on systematic investigation of particle capture efficiency in Si₃N₄ multipore systems of various array configurations. We demonstrated enhanced detection throughput with increasing numbers of pore channels in a membrane. Meanwhile, we also observed significant contributions of the interchannel crosstalk in closely integrated multipores that tended to deteriorate throughput performance by causing shrinkage of the absorption zone via the interference-derived weakening of the electric field around the pore orifice. At the same time, the interference-derived electric field distributions were also found to diminish the electroosmotic contributions to the particle capture efficiency. The present findings can be useful in designing pore arrays with optimal throughput performance.

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.

Anomalous packing and dynamics of a polymer chain confined in a static porous environment

Polymers in confined porous environments are ubiquitous throughout biology, physics, materials science, and engineering. Several experiments have suggested that in some porous environments, chain dynamics can become extremely slow. While phenomenological explanations exist, the exact mechanisms for these slow dynamics have not been fully characterized. In this work, we initiate a joint simulation-theory study to investigate chain packing and dynamics in a static porous environment. The main theoretical concept is the free energy of the chain partitioning into several chambers of the porous environment. Both the theoretical results and Langevin dynamics simulations show that chain packing in each of the chambers is predominantly independent of chain length; it is determined by the maximal packing of segments in each chamber. Dynamically, short chains (compared to the chamber size) become trapped in a single chamber and dynamics become extremely slow, characteristic of an Ogston sieving-like behavior. For longer chains, on the other hand, a hierarchy of slow dynamics is observed due to entropic trapping, characterized by sub-diffusive behavior and a temporary plateau in the mean square displacement. Due to the slow nature of the dynamics, the inevitable long-time diffusive behavior of the chains is not captured by our simulations. Theoretically, the slow dynamics are understood in terms of a free energy barrier required to thread the chain from one chamber to the next. There is overall qualitative and quantitative agreement between simulations and theory. This work provides foundations for a better understanding of how chain dynamics are affected by porous environments.

Effects of Nanopore Charge Decorations on the Translocation Dynamics of DNA

We have investigated the dynamics of single-stranded DNA as it translocates through charge-mutated protein nanopores. Translocation of DNA is a crucial step in nanopore-based sequencing platforms, where control over translocation speed remains one of the main challenges. Taking advantage of the interactions between negatively charged DNA and positively charged amino acid residues, the translocation speed of DNA can be manipulated by deliberate charge decorations inside the nanopore. We employed coarse-grained Langevin dynamics simulations to monitor the step-by-step movement of DNA through different mutations of α-hemolysin protein nanopores. We found that although the average translocation time per nucleotide is longer, in agreement with experiments, the DNA nucleotides do not translocate with a uniform speed. Furthermore, the location and spacing of the charge decorations can alter the translocation dynamics significantly, trapping DNA in some cases. Our findings can give insights when designing charge patterns in nanopores.

The Semiflexible Polymer Translocation into Laterally Unbounded Region between Two Parallel Flat Membranes

Using the dynamic Monte Carlo method, we investigate dynamics of semiflexible polymer translocation through a nanopore into laterally unbounded region between two parallel flat membranes with separation R in presence of an electric field inside the pore. The average translocation time τ initially decreases rapidly with increase of R in the range of R < 10 and then almost keeps constant for R ≥ 10, and the decline range increases with increase of dimensionless bending stiffness κ. We mainly study the effect of chain length N, κ and electric field strength E on the translocation process for R = 5. The translocation dynamics is significantly altered in comparison to an unconfined environment. We find τ ~ N^α, where the exponent α increases with increase of E for small κ. α initially increases slowly with increase of E and then keeps constant for moderate κ. α decreases with increase of E for large κ. However, α decreases with increase of κ under various E. In addition, we find τ ~ κ^β. β decreases with increase of N under various E. These behaviors are interpreted in terms of the probability distribution of translocation time and the waiting time of an individual monomer segment passing through the pore during translocation.

Numerical investigation of molecular nano-array in potential-energy profile for a single dsDNA

A Rigorous numerical investigation on dsDNA translocation in quasi-2-dimensional nano-array filter is performed using Molecular Dynamics (MD) method. Various dsDNA molecules with different sizes are chosen in order to model Ogston sieving in a nano-array filter. The radius of gyration of dsDNA molecule is less than the characteristic length of the shallow region in nano-array. The dsDNA molecule is assumed to be in the 0.05M NaCl electrolyte. MD shows acceptable results for potential-energy profile for nano-array filter. According to the MD outcomes, the dsDNA electrophoretic mobility decreases almost linearly with dsDNA size and show the same trend as Ogston sieving for gel electrophoresis. In addition, different shapes for nano-array filter are studied for a unique dsDNA molecule. It is concluded that steeping the nano-array wall can cause the retardation of dsDNA translocation and decreases dsDNA electrophoretic mobility.

Fluctuating bottleneck model studies on kinetics of DNA escape from α-hemolysin nanopores

We have proposed a fluctuation bottleneck (FB) model to investigate the non-exponential kinetics of DNA escape from nanometer-scale pores. The basic idea is that the escape rate is proportional to the fluctuating cross-sectional area of DNA escape channel, the radius r of which undergoes a subdiffusion dynamics subjected to fractional Gaussian noise with power-law memory kernel. Such a FB model facilitates us to obtain the analytical result of the averaged survival probability as a function of time, which can be directly compared to experimental results. Particularly, we have applied our theory to address the escape kinetics of DNA through α-hemolysin nanopores. We find that our theoretical framework can reproduce the experimental results very well in the whole time range with quite reasonable estimation for the intrinsic parameters of the kinetics processes. We believe that FB model has caught some key features regarding the long time kinetics of DNA escape through a nanopore and it might provide a sound starting point to study much wider problems involving anomalous dynamics in confined fluctuating channels.

Evaluating the applicability of the Fokker-Planck equation in polymer translocation: a Brownian dynamics study

Brownian dynamics (BD) simulations are used to study the translocation dynamics of a coarse-grained polymer through a cylindrical nanopore. We consider the case of short polymers, with a polymer length, N, in the range N = 21-61. The rate of translocation is controlled by a tunable friction coefficient, γ0p, for monomers inside the nanopore. In the case of unforced translocation, the mean translocation time scales with polymer length as <τ1> ∼ (N - Np)(α), where Np is the average number of monomers in the nanopore. The exponent approaches the value α = 2 when the pore friction is sufficiently high, in accord with the prediction for the case of the quasi-static regime where pore friction dominates. In the case of forced translocation, the polymer chain is stretched and compressed on the cis and trans sides, respectively, for low γ0p. However, the chain approaches conformational quasi-equilibrium for sufficiently large γ0p. In this limit the observed scaling of <τ1> with driving force and chain length supports the Fokker-Planck (FP) prediction that <τ> ∝ N/fd for sufficiently strong driving force. Monte Carlo simulations are used to calculate translocation free energy functions for the system. The free energies are used with the FP equation to calculate translocation time distributions. At sufficiently high γ0p, the predicted distributions are in excellent agreement with those calculated from the BD simulations. Thus, the FP equation provides a valid description of translocation dynamics for sufficiently high pore friction for the range of polymer lengths considered here. Increasing N will require a corresponding increase in pore friction to maintain the validity of the FP approach. Outside the regime of low N and high pore friction, the polymer is out of equilibrium, and the FP approach is not valid.

Conformation-dependent translocation of a star polymer through a nanochannel

The translocation process of star polymers through a nanochannel is investigated by dissipative particle dynamics simulations. The translocation process is strongly influenced by the star arm arrangement as the polymer enters the channel, and a scaling relation between the translocation time [Formula: see text] and the total number of beads N tot is obtained. Qualitative agreements are found with predictions of the nucleation and growth model for linear block co-polymer translocation. In the intermediate stage where the center of the star polymer is at the channel entrance, the translocation time is found to have power law-dependence on the number of arms outside the channel and very weakly dependent on the number of arms in the channel. Increasing the total number of star arms also increases the star translocation time.

Polymer translocation dynamics in the quasi-static limit

Monte Carlo (MC) simulations are used to study the dynamics of polymer translocation through a nanopore in the limit where the translocation rate is sufficiently slow that the polymer maintains a state of conformational quasi-equilibrium. The system is modeled as a flexible hard-sphere chain that translocates through a cylindrical hole in a hard flat wall. In some calculations, the nanopore is connected at one end to a spherical cavity. Translocation times are measured directly using MC dynamics simulations. For sufficiently narrow pores, translocation is sufficiently slow that the mean translocation time scales with polymer length N according to <τ> ∝ (N - N(p))(2), where N(p) is the average number of monomers in the nanopore; this scaling is an indication of a quasi-static regime in which polymer-nanopore friction dominates. We use a multiple-histogram method to calculate the variation of the free energy with Q, a coordinate used to quantify the degree of translocation. The free energy functions are used with the Fokker-Planck formalism to calculate translocation time distributions in the quasi-static regime. These calculations also require a friction coefficient, characterized by a quantity N(eff), the effective number of monomers whose dynamics are affected by the confinement of the nanopore. This was determined by fixing the mean of the theoretical distribution to that of the distribution obtained from MC dynamics simulations. The theoretical distributions are in excellent quantitative agreement with the distributions obtained directly by the MC dynamics simulations for physically meaningful values of N(eff). The free energy functions for narrow-pore systems exhibit oscillations with an amplitude that is sensitive to the nanopore length. Generally, larger oscillation amplitudes correspond to longer translocation times.

Translocation dynamics of freely jointed Lennard-Jones chains into adsorbing pores

Polymer translocation into adsorbing nanopores is studied by using the Fokker-Planck equation of chain diffusion along the energy landscape calculated with Monte Carlo simulations using the incremental gauge cell method. The free energy profile of a translocating chain was found by combining two independent sub-chains, one free but tethered to a hard wall, and the other tethered inside an adsorbing pore. Translocation dynamics were revealed by application of the Fokker-Planck equation for normal diffusion. Adsorption of polymer chains into nanopores involves a competition of attractive adsorption and repulsive steric hindrance contributions to the free energy. Translocation times fell into two regimes depending on the strength of the adsorbing pore. In addition, we found a non-monotonic dependence of translocation times with increasing adsorption strength, with sharp peak associated with local free energy minima along the translocation coordinate.

Effects of an attractive wall on the translocation of polymer under driving

The effects of an attractive wall at the trans side on the translocation of an eight-site bond-fluctuation model (BFM) polymer through a pore in a membrane under driving are simulated by the dynamic Monte Carlo method. The attractive wall shows two contrary effects: its excluded volume effect reduces configuration entropy and thus hinders the translocation of the polymer, while its attraction decreases the energy and thus accelerates the translocation. At a critical polymer-wall interaction ε* ≈- 1, we find that the two effects compensate each other and the translocation time τ is roughly independent of the separation distance between the wall and the pore. The value ε* ≈- 1 is roughly equal to the critical adsorption point for the BFM polymer. Moreover, the value of the critical attraction is roughly independent of chain length N and chemical potential difference Δμ. At last, a scaling relation τ ∼ N(α) is observed for polymer translocation at a high value of NΔμ. Though the translocation time is highly dependent on the polymer-wall interaction and pore-wall separation distance, the exponent α is always about 1.30 ± 0.05 so long as NΔμ is large enough.

Single-polymer dynamics under constraints: scaling theory and computer experiment

The relaxation, diffusion and translocation dynamics of single linear polymer chains in confinement is briefly reviewed with emphasis on the comparison between theoretical scaling predictions and observations from experiment or, most frequently, from computer simulations. Besides cylindrical, spherical and slit-like constraints, related problems such as the chain dynamics in a random medium and the translocation dynamics through a nanopore are also considered. Another particular kind of confinement is imposed by polymer adsorption on attractive surfaces or selective interfaces--a short overview of single-chain dynamics is also contained in this survey. While both theory and numerical experiments consider predominantly coarse-grained models of self-avoiding linear chain molecules with typically Rouse dynamics, we also note some recent studies which examine the impact of hydrodynamic interactions on polymer dynamics in confinement. In all of the aforementioned cases we focus mainly on the consequences of imposed geometric restrictions on single-chain dynamics and try to check our degree of understanding by assessing the agreement between theoretical predictions and observations.

Polymer translocation through alpha-hemolysin pore with tunable polymer-pore electrostatic interaction

We have measured the ionic current blockages produced by single molecules of sodium poly(styrene sulfonate) passing through an alpha-hemolysin protein pore under an electric field. Most of the blockage events were composed of one or two blockage levels of ionic current. By analyzing the statistics of different event types for different polymer lengths, applied voltages, and pH conditions, we have identified the molecular mechanism behind the two-level blockages. Our analysis of the data shows that not all blockages are successful translocation events and the propensity of successful translocation can be tuned by pH gradients across the protein pore. We interpret our results as the change in protein-polymer interaction via protonation of charged amino acid residues of alpha-hemolysin pore. In addition, we have constructed a stochastic theory for polymer translocation through alpha-hemolysin pore with tunable polymer-pore interactions. The theoretical calculations capture many features observed in our experiments.

Cooperative translocation dynamics of biopolymer chains through nanopores in a membrane: Slow dynamics limit

The cooperative translocation dynamics of two complementary single-stranded DNA chains through two nanopores located in a membrane is investigated theoretically. The translocation process is considered to be quasi-equilibrium, and then under the limit of slow dynamics the average translocation times are numerically presented under different conditions. It is shown that the effects of the chemical potential gradient, the recombination energy and the distance between the two nanopores on the cooperative translocation are significant. The present model predicts that the cooperative translocation of such two chains can shorten the translocation time, reduce the backward motion and thus improve the translocation efficiency.

Theory of capture rate in polymer translocation

The translocation of macromolecules through a nanopore requires the impingement of the molecules at the pore followed by threading through the pore. While most of the discussion on the translocation phenomenon focused so far on the threading process, the phenomenology on the frequency of encounters between the polymer and the pore exhibits diverse features in terms of polymer length, solution conditions, driving force, and pore geometry. We derive a general theory for the capture rate of polyelectrolyte molecules and the probability of successful translocation through a nanopore, under an externally imposed electric field. By considering the roles of entropic barrier at the pore entrance and drift of the polyelectrolyte under the electric field, we delineate two regimes: (a) entropic barrier regime and (b) drift regime. In the first regime dominated by the entropic barrier for the polyelectrolyte, the capture rate is an increasing nonlinear function in the electric field and chain length. In the drift regime, where the electric field dwarfs the role of entropic barriers, the capture rate is independent of chain length and linear in electric field. An analytical formula is derived for the crossover behavior between these regimes, and the general results are consistent with various experimentally observed trends.

Polymer translocation into a fluidic channel through a nanopore

Using two-dimensional Langevin dynamics simulations, we investigate the dynamics of polymer translocation into a fluidic channel with diameter R through a nanopore under a driving force F . Due to the crowding effect induced by the partially translocated monomers, the translocation dynamics is significantly altered in comparison to an unconfined environment, namely, we observe a nonuniversal dependence of the translocation time τ on the chain length N . τ initially decreases rapidly and then saturates with increasing R , and a dependence of the scaling exponent α of τ with N on the channel width R is observed. The otherwise inverse linear scaling of τ with F breaks down and we observe a minimum of α as a function of F . These behaviors are interpreted in terms of the waiting time of an individual segment passing through the pore during translocation.