Flow-induced translocation of polymers through a fluidic channel: a dissipative particle dynamics simulation study

The dynamics of flow-induced translocation of polymers through a fluidic channel has been studied by dissipative particle dynamics (DPD) approach. Unlike implicit solvent models, the many-body energetic and hydrodynamic interactions are preserved naturally by incorporating explicit solvent particles in this approach. The no-slip wall boundary and the adaptive boundary conditions have been implemented in the modified DPD approach to model the hydrodynamic flow within a specific wall structure of fluidic channel and control the particles' density fluctuations. The results show that the average translocation time versus polymer chain length satisfies a power-law scaling of τ ～N(1.152). The conformational changes and translocation dynamics of polymers through the fluidic channel have also been investigated in our simulations, and two different translocation processes, i.e., the single-file and double-folded translocation events, have been observed in detail. These findings may be helpful in understanding the conformational and dynamic behaviors of such polymer and/or DNA molecules during the translocation processes.     

Effect of Hydrodynamic Interaction on Flow-induced Polymer Translocation Through a Nanotube

We investigated the effect of hydrodynamic interaction(HI) on flow-induced polymer translocation through a nanotube by Brownian dynamics simulations. Whether there is HI in the simulation system is separately controlled by using different diffusion tensors. It is found that HI has no effect on critical velocity flux for long polymer chains due to the competition between more drag force and the hindrance of chain stretching from HI, however, HI broadens the transition interval. In addition, for flow-induced polymer translocation with HI, the critical velocity flux firstly slowly decreases with the increase of chain length and then becomes identical to that of it without HI, that is, the critical velocity flux is independent of chain length. At the same time, HI also accelerates the translocation process and makes the relative variation amplitude of single bead translocation time smaller. In fact, HI can enhance the intrachain cooperativity to make the whole chain obtain more drag force from fluid field and hinder chain stretching, both of which play an important role in translocation process.     

Effect of charge distribution on the translocation of an inhomogeneously charged polymer through a nanopore

We investigate the voltage-driven translocation of an inhomogeneously charged polymer through a nanopore by utilizing discrete and continuous stochastic models. As a simplified illustration of the effect of charge distribution on translocation, we consider the translocation of a polymer with a single charged site in the presence and absence of interactions between the charge and the pore. We find that the position of the charge that minimizes the translocation time in the absence of pore-polymer interactions is determined by the entropic cost of translocation, with the optimum charge position being at the midpoint of the chain for a rodlike polymer and close to the leading chain end for an ideal chain. The presence of attractive and repulsive pore-charge interactions yields a shift in the optimum charge position toward the trailing end and the leading end of the chain, respectively. Moreover, our results show that strong attractive or repulsive interactions between the charge and the pore lengthen the translocation time relative to translocation through an inert pore. We generalize our results to accommodate the presence of multiple charged sites on the polymer. Our results provide insight into the effect of charge inhomogeneity on protein translocation through biological membranes.     

Effect of Interaction upon Translocation of Confined Polymer Chain Through Nanopore

The effect of the interaction between nanopore and chain monomer on the translocation of a single polymer chain confined in a finite size square through an interacting nanopore to a large space has been studied by two-dimensional bond fluctuation model with Monte Carlo simulation. Results indicate that the free energy barrier before the successful translocation of the chain depends linearly on the chain length as well as the nanopore length for different pore-polymer interaction, and the attractive interaction reduces the free energy barrier, leading to the reduction of the average trapping time.     

Simulation Study on Translocation of Confined Chain Through Interacting Nanopore

The translocation of a confined polymer chain through an interacting nanopore has been studied using two-dimensional bond fluctuation model with Monte Carlo dynamics. For different pore-polymer interaction, the average escaping time〈Τesc〉of the polymer chain through the nanopore, increases roughly linearly with the chain length and the nanopore length, respectively. However, the large repulsive and attractive pore-polymer interaction adds the difficulty of the monomers of the chain entering and leaving the nanopore, respec-tively, leading to the nonmonotonical dependence of〈Τesc〉on the pore-polymer interaction. The detailed translocation dynamics of the chain through the interacting nanopore is inves-tigated too.     

Numerical and theoretical study on the mechanism of biopolymer translocation process through a nano-pore

We conducted a numerical study on the translocation of a biopolymer from the cis side to the trans side of a membrane through a synthetic nano-pore driven by an external electric field in the presence of hydrodynamic interactions (HIs). The motion of the polymer is simulated by 3D Langevin dynamics technique using a worm-like chain model of N identical beads, while HI between the polymer and fluid are incorporated by the lattice Boltzmann equation. The translocation process is induced by electrophoretic force, which sequentially straightens out the folds of the initial random configuration of the polymer chain on the cis side. Our simulation results on translocation time and velocity are in good quantitative agreement with the corresponding experimental ones when the surface charge on the nano-pore and the HI effect are considered explicitly. We found that the translocation velocity of each bead inside the nano-pore mainly depends upon the length of the straightened portion of the polymer in forced motion near the pore. We confirmed this by a theoretical formula. After performing simulations with different pore lengths, we observed that translocation velocity mainly depends upon the applied potential difference rather than upon the electric field inside the nano-pore.     

Dynamics of polymer translocation into a circular nanocontainer through a nanopore

Using Langevin dynamics simulations, we investigate the dynamics of polymer translocation into a circular nanocontainer through a nanopore under a driving force F. We observe that the translocation probability initially increases and then saturates with increasing F, independent of φ, which is the average density of the whole chain in the nanocontainer. The translocation time distribution undergoes a transition from a Gaussian distribution to an asymmetric distribution with increasing φ. Moreover, we find a nonuniversal scaling exponent of the translocation time as chain length, depending on φ and F. These results are interpreted by the conformation of the translocated chain in the nanocontainer and the time of an individual segment passing through the pore during translocation.     

First passage time distribution of chaperone driven polymer translocation through a nanopore: homopolymer and heteropolymer cases

Combining the advection-diffusion equation approach with Monte Carlo simulations we study chaperone driven polymer translocation of a stiff polymer through a nanopore. We demonstrate that the probability density function of first passage times across the pore depends solely on the Pe?clet number, a dimensionless parameter comparing drift strength and diffusivity. Moreover it is shown that the characteristic exponent in the power-law dependence of the translocation time on the chain length, a function of the chaperone-polymer binding energy, the chaperone concentration, and the chain length, is also effectively determined by the Pe?clet number. We investigate the effect of the chaperone size on the translocation process. In particular, for large chaperone size, the translocation progress and the mean waiting time as function of the reaction coordinate exhibit pronounced sawtooth-shapes. The effects of a heterogeneous polymer sequence on the translocation dynamics is studied in terms of the translocation velocity, the probability distribution for the translocation progress, and the monomer waiting times.     

Chaperone-assisted translocation of a polymer through a nanopore

Using Langevin dynamics simulations, we investigate the dynamics of chaperone-assisted translocation of a flexible polymer through a nanopore. We find that increasing the binding energy ε between the chaperone and the chain and the chaperone concentration N(c) can greatly improve the translocation probability. Particularly, with increasing the chaperone concentration a maximum translocation probability is observed for weak binding. For a fixed chaperone concentration, the histogram of translocation time τ has a transition from a long-tailed distribution to a gaussian distribution with increasing ε. τ rapidly decreases and then almost saturates with increasing binding energy for a short chain; however, it has a minimum for longer chains at a lower chaperone concentration. We also show that τ has a minimum as a function of the chaperone concentration. For different ε, a nonuniversal dependence of τ on the chain length N is also observed. These results can be interpreted by characteristic entropic effects for flexible polymers induced by either the crowding effect from a high chaperone concentration or the intersegmental binding for the high binding energy.     

流场驱动高分子链迁移穿过微通道的耗散粒子动力学模拟

利用耗散粒子动力学模拟方法研究了高分子链在流场驱动作用下迁移穿过微通道过程中的链构象变化和动力学行为.在足够大的流场力驱动作用下,高分子链在沿着流场方向逐渐被拉伸,从而能够穿过管径小于其自身尺寸的微通道.耗散粒子动力学模拟结果表明高分子链的迁移过程主要分为3个步骤:(1)在流场驱动作用下,高分子链漂移并逐渐靠近微通道入口;(2)高分子链逐渐调整自身构象,并使其部分进入微通道;(3)高分子链成功穿过微通道.同时,模拟还发现当高分子链尺寸大于微通道细管道管径时,高分子链穿过微通道所需的平均迁移时间随着流量的增加而逐渐减小.此外,为了研究高分子链刚性对高分子链穿过微通道的影响,模型中还引入了蠕虫状高分子链模型.模拟结果发现,高分子链的链刚性越强,其迁移穿过微通道的时间越长.     

Polymer translocation through a cylindrical channel

A formalism of polymer translocation through a cylindrical channel of finite diameter and length between two spherical compartments is developed. Unlike previous simplified systems, the finite diameter of the channel allows the number of polymer segments inside the channel to be adjusted during translocation according to the free energy of possible conformations. The translocation process of a Gaussian chain without excluded volume and hydrodynamic interactions is studied using exact formulas of confinement free energy under this formalism. The free energy landscape for the translocation process, the distribution of the translocation time, and the average translocation time are presented. The complex dependencies of the average translocation time on the length and diameter of the channel, the sizes of the donor and receptor compartments, and the chain length are illustrated.     

Unforced translocation of a polymer chain through a nanopore: the solvent effect

The authors have performed the Langevin dynamics simulation to investigate the unforced polymer translocation through a narrow nanopore in an impermeable membrane. The effects of solvent quality controlled by the attraction strength lambda of the Lennard-Jones cosine potential between polymer beads and beads on two sides of the membrane on the translocation processes are extensively examined. For polymer translocation under the same solvent quality on both sides of the membrane, the two-dimensional and three-dimensional simulations confirm the scaling law of tautrans approximately N1+2upsilon for the translocation in the good solvent, where tautrans is the translocation time, N is the chain length, and upsilon is the Flory exponent. For the three-dimensional polymer translocation under different solvent qualities on two sides of the membrane, the translocation efficiency may be notably improved. The scaling law between tautrans and N varies from tautrans approximately N1+2upsilon to tautrans approximately N with the increase of the difference of solvent qualities, and the crossover occurs at the theta temperature point, where a scaling law of tautrans approximately N1.27 is found. The simulation results here also show that the translocation time changes from a wide and asymmetric distribution with a long tail to a narrow and symmetric distribution with the increase of the difference of the solvent qualities.     

单链高分子在纳米孔隙中的跨膜输运问题

针对单链高分子通过薄膜上纳米孔隙的输运这一现象分别首次实验和理论方面介绍了现今的研究进展情况.文中以理论和实验上共同关注的单链高分子跨膜输运的平均通过时间问题展开论述.具体地,作者考察了薄膜的厚度、薄膜的曲率、膜两侧的溶剂性质、单链高分子聚集态的变化以及链间相互作用等因素对输运过程的影响,并讨论了与这些因素相对应的实际应用.同时,分别对几个研究组的实验结果、解析以及模拟的理论结果进行了分析比较,旨在促进这一现象在相关领域的研究.     

Polymer translocation through a nanopore. II. Excluded volume effect

Following our previous study of a Gaussian chain translocation, we have investigated the transport of a self-avoiding chain from one sphere to another sphere through a narrow pore, using the self-consistent field theory formalism. The free energy landscape for polymer translocation is significantly modified by excluded volume interactions among monomers. The free energy barrier for the placement of one of the chain ends at the pore depends on the chain length N nonmonotonically, in contrast to the N-independence for Gaussian chains. This results in a nonmonotonic dependence of the average arrival time [tau0] on N for self-avoiding chains. When the polymer chain is partitioned between the donor and recipient spheres, a local free energy minimum develops, depending on the strength w of the excluded volume interaction and the relative sizes of the donor and recipient spheres. If the sizes of spheres are comparable, the average translocation time tau (the average time taken by the polymer, after the arrival at the pore, to convert from the donor to the recipient) increases with an increase in w for a fixed N value. On the other hand, for the highly asymmetric sizes of the donor and recipient spheres, tau decreases with an increase in w. As in the case of Gaussian chains, tau depends nonmonotonically on the pore length.     

单链高分子通过薄膜上纳米孔隙的输运过程: 链间相互作用的影响

从理论上研究了高分子链通过位于薄膜上的一个纳米孔隙的输运问题. 考虑在膜的一侧引入与传输高分子链有特殊相互作用的高分子链段, 研究链间相互作用致使高分子链构型的改变对其平均首次通过时间的影响, 进而在不同条件下进行了相应的数值模拟, 并讨论了其它一些相关问题.     

Monte Carlo simulation on polymer translocation in crowded environment

The effect of crowded environment with static obstacles on the translocation of a three-dimensional self-avoiding polymer through a small pore is studied using dynamic Monte Carlo simulation. The translocation time τ is dependent on polymer-obstacle interaction and obstacle concentration. The influence of obstacles on the polymer translocation is explained qualitatively by the free energy landscape. There exists a special polymer-obstacle interaction at which the translocation time is roughly independent of the obstacle concentration at low obstacle concentration, and the strength of the special interaction is roughly independent of chain length N. Scaling relation τ ~ N(1.25) is observed for strong driving translocations. The diffusion property of polymer chain is also influenced by obstacles. Normal diffusion is only observed in dilute solution without obstacles or in a crowded environment with weak polymer-obstacle attraction. Otherwise, subdiffusion behavior of polymer is observed.     

Excluded volume effect on confined polymer translocation through a short nanochannel

We simulated the translocation process of a polymer chain from a source container to a drain container through a short nanochannel. We utilized the bond fluctuation model coupled with Monte Carlo dynamics in our simulations. The calculation results show that the excluded volume effect significantly affects the polymer's translocation time tau. This time depends nonmonotonically on the polymer length N. For a fixed nanochannel length, tau decreases when the polymer length increases. tau, however, increases when the polymer length exceeds a certain threshold. This observation differs from those predicated for a Gaussian chain. In this paper, we will further present our findings to explain this phenomenon. The knowledge we gain from this research can enhance the understanding of complex transport processes in many biological systems.     

Effect of orientation in translocation of polymers through nanopores

The motion of polymers with inhomogeneous structure through nanopores is discussed theoretically. Specifically, we consider the translocation dynamics of polymers consisting of double-stranded and single-stranded blocks. Since only the single-stranded chain can go through the nanopore the double-stranded segment has to unzip before the translocation. Utilizing a simple analytical model, translocation times are calculated explicitly for different polymer orientations, i.e., when the single-stranded block enters the pore first and when the double-stranded segment is a leading one. The dependence of the translocation dynamics on external fields, energy of interaction in the double-stranded segment, size of the polymer, and the fraction of double-stranded monomers is analyzed. It is found that the order of entrance into the pore has a significant effect on the translocation dynamics. The theoretical results are discussed using free-energy landscape arguments.     

Theoretical study on the translocation of partially charged polymers through nanopore

The translocation time τ of partially charged polymers through a neutral nanopore is calculated using Fokker–Planck equation with adsorbing–adsorbing boundary conditions. For the polymer with one charged monomer, we find that τ is dependent on the position of the charged monomer and on the magnitude of the driving force f inside the nanopore. When the charge is located at the front half of the polymer chain, τ is larger than that of neutral polymer and increases with f. When the charge is located at the back half, it is smaller than that of the neutral polymer and decreases with increasing f. We have also studied the behavior of a symmetrical polymer with two like charges located symmetrically in the chain and that of an asymmetrical polymer with two unlike charges. Moreover, we have calculated the translocation time for a general condition of polymer with two randomly distributed charges. All results show that τ is dependent on the positions of charges in the polymer chain and on the magnitude of the driving force. The results can be explained qualitatively by the free‐energy landscape of polymer translocation. © 2017 Wiley Periodicals, Inc. J. Polym. Sci., Part B: Polym. Phys. 2017 , 55, 1017–1025     

Coupling Effects of Electrostatic Interactions and Salt Concentration Gradient in Polymer Translocation through a Nanopore: A Coarse-Grained Molecular Dynamics Simulations Study

We study the influence of polymer pore interactions and focus on the role played by the concentration gradient of salt in the translocation of polyelectrolytes (PE) through nanopores explicitly using coarse-grained Langevin dynamics simulations. The mean translocation time is calculated by varying the applied voltage, the pH, and the salt concentration gradient. Changing the pH can alter the electrostatic interaction between the protein pore and the polyelectrolyte chain. The polymer pore interaction is weakened by the increase in the strength of the externally applied electric field that drives translocation. Additionally, the screening effect of the salt can reduce the strong charge-charge repulsion between the PE beads which can make translocation faster. The simulation results show there can be antagonistic or synergistic coupling between the salt concentration-induced screening effect and the drift force originating from the salt concentration gradient thereby affecting the translocation time. Our simulation results are explained qualitatively with free energy calculations.