Sequence dependence of DNA translocation through a nanopore

We investigate the dynamics of DNA translocation through a nanopore using 2D Langevin dynamics simulations, focusing on the dependence of the translocation dynamics on the details of DNA sequences. The DNA molecules studied in this work are built from two types of bases A and C, which have been shown previously to have different interactions with the pore. We study DNA with repeating blocks A(n)C(n) for various values of n and find that the translocation time depends strongly on the block length 2n as well as on the orientation of which base enters the pore first. Thus, we demonstrate that the measurement of translocation dynamics of DNA through a nanopore can yield detailed information about its structure. We have also found that the periodicity of the block sequences is contained in the periodicity of the residence time of the individual nucleotides inside the pore.     

Effect of Interaction upon Translocation of Confined Polymer Chain Through Nanopore

运用二维的键长涨落模型和蒙特卡洛方法研究高分子链从一个受限空间到自由空间穿孔过程中,链单体与纳米孔之间的相互作用.结果表明,在不同的链长和纳米孔交互作用下,高分子链成功穿越自由能能垒取决于链长和纳米孔长度,并且由于交互作用降低了自由能能垒,导致高分子链在纳米管的平均捕获时间缩短.     

Free energy and scalings for polymer translocation through a nanopore: A molecular dynamics simulation study combined with milestoning

Coarse-grained molecular dynamics simulations combined with milestoning method are used to study the stochastic process of polymer chain translocation though a nanopore. We find that the scalings for polymer translocation process (the chain is initialized with the first monomer in the nanopore) and for polymer escape process (the chain is initialized with the middle monomer in the nanopore) are different. The translocation process is mainly controlled by the entropic barrier, while the polymer escape process is driven by the effective force due to free energy difference.     

纳米孔壁面作用对蛋白质过孔影响的粗粒化分子动力学模拟

基于粗粒化分子动力学方法模拟电驱动蛋白质过孔过程，研究纳米孔-水/纳米孔-蛋白质相互作用对电泳迁移率的影响；用操控式分子动力学模拟分析蛋白质在不同相互作用下过孔摩擦系数和摩擦阻力.研究发现：蛋白质黏附纳米孔壁面对其过孔特性影响并不明显，而纳米孔-水相互作用对蛋白质过孔电泳迁移率和摩擦系数影响较大.随纳米孔-水相互作用增强，纳米孔壁面与蛋白质附近水分子运动差异显现，蛋白质过孔摩擦阻力显著增大，过孔摩擦系数随之增大，进而影响蛋白质过孔电泳迁移率.所得结果可为纳米孔材料设计提供理论指导.     

Translocation of Single Polymer Chain from Nanopore on a Membrane： Solvent Effect

We investigate the translocation of single polymer chain through a nanopore located on a membrane with different solvents in the two sides of the membrane. For the case under study, the effect of solvents on the translocation dynamics is significant, and as a result, the mean first passage time shortens remarkably compared with that calculated in the case of good solvents on both the sides of the membrane. In addition, we also discuss the condition such that the present result holds true.     

脱氧核糖核酸单链序列的预测

利用Langevin动力学方法模拟了脱氧核糖核酸(DNA)单链在电场力作用下穿越纳米孔道的动力学过程.研究表明,不同种类的单体对应着不同的居留时间,相邻单体的居留时间随着孔道长度的增大而减小.在简化模型的基础上,可以从居留时间图中一次性地推测出一条DNA链的嘌呤和嘧啶的分布.应用该方法对17条不同序列的DNA链进行了预测,平均准确率为951%.在此方法的基础上做一些改进,可以为DNA链的测序提供一种高效的低成本方法. 关键词： Langevin动力学 脱氧核糖核酸单链 序列预测     

Influence of polymer-pore interactions on translocation

We investigate the influence of polymer-pore interactions on the translocation dynamics using Langevin dynamics simulations. An attractive interaction can greatly improve the translocation probability. At the same time, it also increases the translocation time slowly for a weak attraction while an exponential dependence is observed for a strong attraction. For fixed driving force and chain length the histogram of translocation time has a transition from Gaussian distribution to long-tailed distribution with increasing attraction. Under a weak driving force and a strong attractive force, both the translocation time and the residence time in the pore show a nonmonotonic behavior as a function of the chain length. Our simulations results are in good agreement with recent experimental data.     

Monte Carlo analysis of polymer translocation with deterministic and noisy electric fields

Polymer translocation through the nanochannel is studied by means of a Monte Carlo approach, in the presence of a static or oscillating external electric voltage. The polymer is described as a chain molecule according to the two-dimensional “bond fluctuation model”. It moves through a piecewise linear channel, which mimics a nanopore in a biological membrane. The monomers of the chain interact with the walls of the channel, modelled as a reflecting barrier. We analyze the polymer dynamics, concentrating on the translocation time through the channel, when an external electric field is applied. By introducing a source of coloured noise, we analyze the effect of correlated random fluctuations on the polymer translocation dynamics.     

Translocation of polymers in a lattice model

Voltage-driven polymer translocation is studied by means of a stochastic lattice model. The model incorporates voltage drop over the membrane as a bias in the hopping rate through the pore and exhibits the two main ingredients of the translocation process: driven motion through the pore and diffusive supply of chain length towards the pore on the cis-side and the drift away from the pore on the trans-side. The translocation time is either bias limited or diffusion limited. In the bias-limited regime the translocation time is inversely proportional to the voltage drop over the membrane. In the diffusion-limited regime the translocation time is independent of the applied voltage, but it is rather sensitive to the motion rules of the model. We find that the whole regime is well described by a single curve determined by the initial slope and the saturation value. The dependence of these parameters on the length of the chain, the motion rules and the repton statistics are established. Repulsion of reptons as well as the increase of chain length decrease the throughput of the polymer through the pore. As for free polymers, the inclusion of a mechanism for hernia creations/annihilations leads to the cross-over from Rouse-like behaviour to reptation. For the experimentally most relevant case (Rouse dynamics) the bimodal power law dependence of the translocation time on the chain length is found.     

Translocation of closed polymers through a nanopore under an applied external field

The dynamic behaviours of the translocations of closed circular polymers and closed knotted polymers through a nanopore, under the driving of an applied field, are studied by three-dimensional Langevin dynamics simulations. The power-law scaling of the translocation time τ with the chain length N and the distribution of translocation time are investigated separately. For closed circular polymers, a crossover scaling of translocation time with chain length is found to be τ～ N α , with the exponent α varying from α = 0.71 for relatively short chains to α = 1.29 for longer chains under driving force F = 5. The scaling behaviour for longer chains is in good agreement with experimental results, in which the exponent α = 1.27 for the translocation of double-strand DNA. The distribution of translocation time D(τ) is close to a Gaussian function for duration time τ < τ p and follows a falling exponential function for duration time τ > τ p . For closed knotted polymers, the scaling exponent α is 1.27 for small field force (F = 5) and 1.38 for large field force (F = 10). The distribution of translocation time D(τ) remarkably features two peaks appearing in the case of large driving force. The interesting result of multiple peaks can conduce to the understanding of the influence of the number of strands of polymers in the pore at the same time on translocation dynamic process and scaling property.     

Two-dimensional diffusion model for the biopolymers dynamics at nanometer scale

In this paper the driven transport of linear polymers through a nanopore is presented. Biopolymer physical behavior in an external electric field is modeled and its motion is simulated using the Langevin impulse integrator method. Within fairly large limits, the polymer translocation time is inversely proportional with the electric field intensity and directly proportional with the polymer chain length.      

Single long-polymer translocation through a long pore

This paper theoretically studies the free energy and conformational entropy of a long polymer threading a long nanopore (n₀/N \ge 0.1) on external electric field. The polymer expanded model is built in this paper, that is, a single long polymer chain with N monomers (each of size a) threading a pore with n₀ monomers can be regarded as polymer with N+n_{0} monomers translocating a 2-dimension hole embedded in membrane. A theoretical approach is presented which explicitly takes into account the nucleation theory. Our calculations imply that, the structure of polymer changes more acutely than other situation, while its leading monomer reaches the second vacuum and its end monomer escapes the first vacuum. And it is also shown that the length scale of polymer and pore play a very important role for polymer translocation dynamics. The present model predicts that the translocation time depends on the chemical potential gradient and the property of the solvent on sides of pore to some extent.     

Dynamics of polymer translocation through nanopores: theory meets experiment

The dynamics of translocation of polymer molecules through nanopores is investigated via molecular dynamics. We find that an off-lattice minimalist model of the system is sufficient to reproduce quantitatively all the experimentally observed trends and scaling behavior. Specifically, simulations show (i) two translocation regimes depending on the ratio of pore and polymer length, (ii) two different regimes for the probability of translocation depending on applied voltage, (iii) an exponential dependence of translocation velocity upon applied voltage, and (iv) an exponential decrease of the translocation time with temperature. We also propose a simple theoretical explanation of each of the observed trends within a free energy landscape framework.     

Polymer translocation in a double-force arrangement

Using Langevin dynamics simulations, we investigate the translocation dynamics of an externally driven polymer chain through a nanopore, where a pulling force F is exerted on the first monomer whilst there is an opposing force F _E < F within the pore. Such a double-force arrangement has been proposed recently to allow better dynamical control of the translocation process in order to sequence biopolymers. We find that in the double-force arrangement translocation becomes slower as compared to the case under a single monomer pulling force of magnitude F - F _E , but scaling of the translocation time as a function of the chain length ∼ N ² does not change. The waiting time (m) for monomer m to exit the pore is found to be a monotonically increasing function of the bead number almost until m N , which indicates relatively well-defined slowing down and control of the chain velocity during translocation. We also study the waiting time distributions for the beads in the chain, and characterize in detail fluctuations in the bead positions and their transverse position coordinates during translocation. These data should be useful in estimating position-dependent sequencing errors in double-force experiments.     

Concentration polarization in translocation of DNA through nanopores and nanochannels

In this Letter we provide a theory to show that high-field electrokinetic translocation of DNA through nanopores or nanochannels causes large transient variations of the ionic concentrations in front and at the back of the DNA due to concentration polarization (CP). The CP causes strong local conductivity variations, which can successfully explain the nontrivial current transients and ionic distributions observed in molecular dynamics simulations of nanopore DNA translocations as well as the transient current dips and spikes measured for translocating hairpin DNA. Most importantly, as the future of sequencing of DNA by nanopore translocation will be based on time-varying electrical conductance, CP, must be considered in experimental design and interpretation--currently these studies are mostly based on the incomplete pore conductance models that ignore CP and transients in the electrical conductance.     

打结高分子链穿孔行为的研究

以三叶草型结(即3₁结)为例,采用分子动力学(MD)方法,研究打结高分子链在外场力作用下穿越微孔的动力学过程.模拟发现,在拉动打结高分子链的过程中,结的大小呈涨落变化,直至最后散结.定性讨论了结的存在对高分子链穿孔速率的影响.在外场力作用下,打结高分子链平均穿孔时间(τ)与链长(N)满足标度关系τ～N ^α,其中标度系数α随外场力f增大而增大.对于短链,外场力越大,平均穿孔时间越短     

Poisson-Fokker-Planck model for biomolecules translocation through nanopore driven by electroosmotic flow

A non-continuous electroosmotic flow model (PFP model) is built based on Poisson equation, Fokker-Planck equation and Navier-Stokse equation, and used to predict the DNA molecule translocation through nanopore. PFP model discards the continuum assumption of ion translocation and considers ions as discrete particles. In addition, this model includes the contributions of Coulomb electrostatic potential between ions, Brownian motion of ions and viscous friction to ion transportation. No ionic diffusion coefficient and other phenomenological parameters are needed in the PFP model. It is worth noting that the PFP model can describe non-equilibrium electroosmotic transportation of ions in a channel of a size comparable with the mean free path of ion. A modified clustering method is proposed for the numerical solution of PFP model, and ion current translocation through nanopore with a radius of 1 nm is simulated using the modified clustering method. The external electric field, wall charge density of nanopore, surface charge density of DNA, as well as ion average number density, influence the electroosmotic velocity profile of electrolyte solution, the velocity of DNA translocation through nanopore and ion current blockade. Results show that the ion average number density of electrolyte and surface charge density of nanopore have a significant effect on the translocation velocity of DNA and the ion current blockade. The translocation velocity of DNA is proportional to the surface charge density of nanopore, and is inversely proportional to ion average number density of electrolyte solution. Thus, the translocation velocity of DNAs can be controlled to improve the accuracy of sequencing by adjusting the external electric field, ion average number density of electrolyte and surface charge density of nanopore. Ion current decreases when the ion average number density is larger than the critical value and increases when the ion average number density is lower than the critical value. Our numerical simulation shows that the translocation velocity of DNA given by the PFP model agrees with the experimental, results better than that given by PNP model or PB model.     

DNA in nanopores: counterion condensation and coion depletion

Molecular dynamics simulations are used to study the equilibrium distribution of monovalent ions in a nanopore connecting two water reservoirs separated by a membrane, both for the empty pore and that with a single stranded DNA molecule inside. In the presence of DNA, the counterions condense on the stretched macromolecule effectively neutralizing it, and nearly complete depletion of coions from the pore is observed. The implications of our results for experiments on DNA translocation through alpha-hemolysin nanopores are discussed.     

Optical detection of DNA translocation through silicon nanopore by ultraviolet light

In this paper, we propose a new optical detection scheme for nanopore-based DNA sequencing with high resolution towards eventual base identification. We use ultraviolet light for excitation of a fluorescent probe attached to DNA and a nanopore in the silicon membrane that has a significantly large refractive index and an extinction coefficient at ultraviolet wavelengths. In this study, numerical electromagnetic simulation revealed that the z-polarization component (perpendicular to the membrane plane) of the electric field was dominant near the nanopore and generated a large electric field gradient at the nanopore exit, typically with a decay length of 2 nm for a nanopore with diameter of 7 nm. The large extinction coefficient contributed to reduction in background noise coming from fluorophore-labeled DNA strands that remain behind the membrane (the cis side of the membrane). We observed a high signal-to-noise ratio of single DNA translocation events under the application of an electric field.