Next Article in Journal
Multi-Physics Fields Based Nonlinear Dynamic Behavior Analysis of Air Bearing Motorized Spindle
Next Article in Special Issue
A Numerical Investigation of Enhancing Microfluidic Heterogeneous Immunoassay on Bipolar Electrodes Driven by Induced-Charge Electroosmosis in Rotating Electric Fields
Previous Article in Journal
Numerical Investigation on the Optimum Thermal Design of the Shape and Geometric Parameters of Microchannel Heat Exchangers with Cavities
Previous Article in Special Issue
Dielectrophoretic Separation of Particles Using Microfluidic Chip with Composite Three-Dimensional Electrode
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Influence of Electric Field Intensity and Particle Length on the Electrokinetic Transport of Cylindrical Particles Passing through Nanopore

1
Mechanical and Electrical Engineering College, Hainan University, Haikou 570228, China
2
Institute of Biomedical Engineering, Chinese Academy of Medical Sciences and Peking Union Medical College, Tianjin 300192, China
3
School of Mechanical Engineering, Yeungnam University, Gyongsan 38541, Korea
*
Authors to whom correspondence should be addressed.
Micromachines 2020, 11(8), 722; https://doi.org/10.3390/mi11080722
Submission received: 27 May 2020 / Revised: 22 July 2020 / Accepted: 22 July 2020 / Published: 25 July 2020
(This article belongs to the Special Issue Electrokinetics in Micro-/nanofluidic Devices)

Abstract

:
The electric transport of nanoparticles passing through nanopores leads to a change in the ion current, which is essential for the detection technology of DNA sequencing and protein determination. In order to further illustrate the electrokinetic transport mechanism of particles passing through nanopores, a fully coupled continuum model is constructed by using the arbitrary Lagrangian–Eulerian (ALE) method. The model consists of the electric field described by the Poisson equation, the concentration field described by Nernst–Planck equation, and the flow field described by the Navier–Stokes equation. Based on this model, the influence of imposed electric field and particle length on the electrokinetic transport of cylindrical particles is investigated. It is found firstly the translation velocities for the longer particles remain constant when they locate inside the nanopore. Both the ion current blockade effect and the ion current enhancement effect occur when cylindrical particles enter and exit the nanopore, respectively, for the experimental parameters employed in this research. Moreover, the particle translation velocity and current fluctuation amplitude are dominated by the electric field intensity, which can be used to adjust the particle transmission efficiency and the ion current detectability. In addition, the increase in particle length changes the particle position corresponding to the peak value of the ion current, which contributes to distinguishing particles with different lengths as well.

1. Introduction

With the completion of the Human Genome Project, the era of post genomics is coming [1]. An urgent demand has been proposed for cheaper and faster DNA sequencing, and thus a new generation of DNA sequencing technology based on the nanopore has been developed rapidly [2,3,4]. Currently, there are two kinds of DNA sequencing technology based on nanopores: optical detection and electrical detection. Optical detection is to hybridize nucleotide and marker, and the marker can be detected when DNA is translocated so as to achieve sequencing [5,6]. Electrical detection mainly includes three detection methods, namely ion current [7,8], tunnel current [9,10], and potential difference [11,12]. Their basic principle is that when DNA molecules pass through the nanopore, they will temporarily affect the conductivity of the nanopore. The sequence of the base pairs can be distinguished based on the different influence of the base pairs on the ion current, the tunnel current and the potential difference between the two ends of the nanopore.
Generally, two kinds of nanopores are employed in the aforementioned analysis and detection technology. One is a biological nanopore, and the other is a solid nanopore. In a pioneering study of biological nanopores [13], α-hemolysin was employed to explore the movement of single-stranded DNA in 1996. It was found that the current blockade time was proportional to the length of the polymerization chain when single-stranded DNA molecules passed through the double-layer membrane structure with a diameter of 2.6 nm, and thus the feasibility of DNA sequencing using nanopores was proven. Subsequently, the movement of DNA in nanopores was studied using phage Phi29 connector [14] and Aerolysin [15]. Due to some intrinsic problems of biological nanopores, associated with difficulty in preparation, unstable chemical properties and poor repeatability, more and more researchers turned to solid nanopores. The first solid nanopore [16] was synthesized solid on silicon nitride thin films for DNA detection by means of an ion beam in 2001. On the basis of successfully preparing nanopores, other researchers further investigated the movement of DNA using nanopores made of SiO2 [17], alumina film [18], graphene nanopore [19], molybdenum disulfide [20], and glass capillary [21,22]. To complement experimental investigations demanding complex preparation processes of solid nanopores and expensive detection equipment, a large number of researchers have applied the continuum model [23,24,25,26,27], molecular dynamics [28,29,30], atomic level Brownian dynamics [31], Brownian kinematics [32] to carry out numerical simulation on particle motion and ion current change in nanopores.
Despite many enlightening reports on particle dynamics through nanopores, theoretical investigations based on numerical simulations deserve further attention. Some key parameters, including the particle length, have not been exploited, and the details of the particle transport mechanism have not been fully elucidated. Herein, a fully coupled continuous model is constructed by the arbitrary Lagrangian−Eulerian (ALE) method. The model is composed of an electric field described by the Poisson equation, a concentration field described by the Nernst−Planck equation, and a flow field described by the Navier−Stokes equation. Based on the model, the electrokinetic mechanism of cylindrical particles is investigated in detail to elucidate the movement of a DNA strand passing through the nanopore. In particular, the effects of the applied electric field intensity and the particle length on the translation velocity in the nanopore and the change in ion current are studied.

2. Mathematical Model

We studied a membrane structure with pore diameter b and thickness h, connecting two nanounits with uniform height H and width W to form a nanopore channel, as shown in Figure 1. The channel was filled with KCl aqueous solution with concentration C0, dynamic viscosity μ , the density ρ f , and permittivity ε f . A Cartesian coordinate system (x, y) was established at the center of the channel. An electric potential difference was imposed between the upper and lower walls of the nanopore channel, so the ions’ motion in the channel generated current under the action of the DC electric field. A negatively charged cylindrical particle with its length Lp and the cap radius a was put into the nanochannel, the cylindrical particle was defined as a rigid domain, its density was set equal to that of the aqueous solution, the gravity and buoyancy on the particle could be regarded as offsetting each other accordingly, then the particle would move continuously under the action of the electrophoretic force and pass through the nanopore [26]. Since the particle displaced a part of the fluid as it passed through the nanopore, the ions passing through the nanopore decreased, and the ionic current was blocked, inducing a detectable ionic current deviation. The fluctuation magnitude of the current was related to the particle size and the applied electrical field strength.
The electric field distribution in the nanopore is described by the Poisson equation
ε f 2 ϕ = F ( z 1 c 1 + z 2 c 2 )
where ε f is the permittivity of fluid in the nanopore, ϕ is the electric potential, F is the Faraday constant, z i is the valence of the i th ionic species, c i is the concentration of the i th ionic species ( i = 1,2), and the right side of the equation represents the bulk charge density of the solution. The potential boundary conditions on the upper and lower walls of the nanopore channel are given as
ϕ ( x , ( H + h / 2 ) ) = 0
ϕ ( x , ( H + h / 2 ) ) = ϕ 0
The particle is charged with surface charge density σ p , which is imposed as a boundary condition on the moving particle surface, where n is the outward unit normal vector on the particle surface. Other rigid boundaries are set to be electrically insulated.
n ϕ =   σ p
The ion distribution in the nanopore is described by the Nernst−Planck equation, expressed as
c i t + N i = 0 , i = 1   ,   2
N i = u f c i D i c i z i D i R T F c i ϕ , i = 1   ,   2
where N i is the flux of the i th ionic species, u f is the flow velocity in the liquid field, D i is the diffusion coefficient of the i th ionic species, R is the universal gas constant, and T is the absolute temperature of the medium. Each term on the right-hand side of Equation (6) represents the convection, the diffusion, and the electro migration of ionic species, respectively. The ion boundary conditions at the upper and lower walls of the nanopore channel are given as
c i ( x , ± ( H + h / 2 ) ) = c 0 , i = 1 , 2
The convective flux dominates on the moving particle surface, and thus the normal ionic flux of particle surface is set as
n N i = n ( u f c i )
There is no ion penetrating the stationary nanopore wall, where the normal flux is set to zero. The upper and lower edges of the two nanounits are set as open boundaries that describe boundaries in contact with a large volume of fluid; fluid can both enter and leave the domain on these boundaries. On the other boundaries, indicated by the dotted lines in Figure 1, no penetration and no shear-stress conditions are prescribed:
u f n = 0
[ p Ι + μ ( u f + u f T ) ] n = 0
The Navier−Stokes equation is applied to describe the fluid motion in the nanopore. Since the Reynolds number (Re) in the nanopore is extremely small, we can ignore the inertia, and so for an incompressible flow the conservation of momentum and mass are expressed by
ρ f u f t = [ p Ι + μ ( u f + u f T ) ] + f
u f = 0
where ρ f is the density of the fluid, u f is the velocity of the flow field, p is the fluid pressure, μ is the dynamic viscosity of fluid, and f is the volume force per unit fluid, f = F ( z 1 c 1 + z 2 c 2 ) ϕ , representing the interactive force between the electric field and the net charge in the fluid. The boundary conditions of the upper and lower walls of the nanopore channel are
p ( x , ± ( H + h / 2 ) ) = 0
The particle in the nanopore is forced by both the electric and the flow field, F p = F E + F f , where the electric force F E and the hydrodynamic force F f are, respectively, obtained by integrating the Maxwell tensor σ E and the hydrodynamic tensor σ f on the particle surface Γ p as follows:
F E = σ E n d Γ p = ε p [ E E 1 2 ( E E ) I ] n d Γ p
F f = σ f n d Γ p = [ p I + μ ( u + u T ) ] n d Γ p
where E represents the electric field intensity determined by E = ϕ , I denotes unit tensor, and ε p represents the permittivity of the particle. The superscript T represents the transpose of the matrix. The particle surface and nanopore wall are set as no slip boundaries, and so the particle surface velocity u is identical to the fluid velocity, u = u f . The other boundaries are set as symmetrical boundaries.
In addition, the ion current passing through the cross section S of the nanopore is obtained by
I = S F ( z i N i ) n d S
where I is the ion current, and N i is the ionic flux corresponding to the cross section S.
Since the magnitude of ion current I passing through the nanopore is very small, the ionic current when the particle is far away from the nanopore is taken as the reference current I0 (base ionic current), and the dimensionless ion current deviation I* is used to express the difference between them as
I = ( I I 0 ) / I 0
I* > 0 means that the ion current passing through the nanopore is greater than the base ionic current, resulting in the ion current enhancement effect, while I* < 0 means that the ion current passing through the nanopore is smaller than the base ionic current, resulting in the ion current blockade effect.
The commercial finite element software COMSOL (version 5.3a, COMSOL Inc., Stockholm, Sweden) was employed to solve the fully coupled electric field, ion transport and particle transport based on the ALE method, which had been proven to be effective in simulating moving boundaries. When the particle started to move (passed the center of the nanopore), the minimum total number of grid and edge elements required were 27,164 (56,900) and 774 (1,357) units, respectively, for convergence, which was also consistent with reference [26]. The model adopted the built-in directive PARDISO solver of COMSOL, with a minimum time step of 8 × 107 [s] and relative tolerance of 0.005 to ensure the convergence of the solution process and the independence of the solution results. In previous reports the integrity of the computational model used was verified by investigating the electrodynamics of deformable or rigid particles inside micro- and nanochannels [33,34,35,36]. Figure 2 additionally shows the high consistency between the numerical simulation used in this study and that in reference [26]. Here we applied this model to further explore the electrokinetic transport mechanism of cylindrical particles in nanopores.

3. Results and Discussion

Among many cases studied we presented representative cases with nanounit size W = 100 nm and H = 200 nm, the membrane structure h = 5 nm and b = 5 nm, and the radius of the particle cap a = 1 nm. The concentration of KCl solution was C0 = 10 mol/m3, while the fluid permittivity, ε f = 7.08 × 10 10   F / m , the permittivity of particle ε p was set identical to ε f , the solution density ρ = 1 × 10 3   kg   / m 3 , the fluid viscosity μ = 1 × 10 3   Pa s , the diffusivity of K + , D 1 = 1.95 × 10 9   m 2 / s , and the diffusivity of C l , D 2 = 2.03 × 10 9   m 2 / s , the solution temperature T = 300 K , and the particle surface charge density σ p = 0.01 C / m 2 . Considering that the pitch of the double-helix structure of the common B-DNA molecule is 3.4 nm and 10 base pairs (bp) in one spiral period, we set 1 bp = 0.34 nm. The coordinate variable of particle center a (Xp, Yp) had its initial position (0, −300 bp), below the nanopore 300 bp away from the nanopore center. It should be noted that although the randomly distributed cylindrical particles generally experienced a process of rotating and lateral offset to a position parallel to the electric field prior to entering the nanopore, the orientation of cylindrical particles generally remained unchanged when they approached and entered the nanopore [23,26]. We thus regarded the motion of the cylindrical particles near the nanopore as translational motion, and the velocity u on the particle surface equaled the velocity of the particle as a whole, which simplified the simulation and subsequent analysis. Here we elucidate the translation mechanism of the particle and current fluctuation in the nanopore with focus on the electric field intensity and the length of DNA strand as parameters.

3.1. The Effect of Imposed Electric Field

The electric field imposed is a dominant parameter dictating the dynamics of particles by electrophoresis, and so it is necessary to study its effect on particle transport in detail. Since the original position of the cylindrical particle is arranged below the nanopore with no lateral offset with respect to the y-axis, the particles always translate along the central axis of the nanopore(see Video S1 in the Supplementary Materials). Figure 3 shows the particle translation velocity depending on the particle position for a particle with length Lp = 200 bp under the applied electric field strength of E = 1, E = 2 and E = 3 MV/m. Before the particle reached the nanopore, its velocity stayed almost constant. As the particle entered the nanopore, the velocity increased rapidly until it reached a peak value near the position yp = −75 bp, after which it decreased conspicuously until the particle center moved to the nanopore center (yp = 0). As the particle center moved beyond the position yp = 75 bp, its velocity decreased monotonically until it exited the nanopore completely and translated away, as shown in the figure through the position yp = 300 bp. With the increase in the electric field strength, the translation velocity increased consistently during the entire transport process of the particle, and this trend was more pronounced when the particle was located in the mid-section (−75bp < yp < 75bp). The peak values calculated for the applied electric field strength E = 1, E = 2, and E = 3 MV/m were approximately 50, 125, and 210 mm/s, respectively.
The above velocity profile can be explained by examining the electrophoretic force acting on the particle and electrostatic interaction between the particles and the nanopore. Figure 4 shows the electric field around the nanopore when the particle center was located at yp = −75, yp = 0, yp = 75bp. As shown, strong electric fields were induced around the particle and the nanopore, which was attributed to the formation of the electric double layer (EDL) adjacent to their surfaces. On the one hand, the imposed electric field generated the electrophoretic force compelling the negative charged particles to pass through the nanopore (when the electric polarity of particle or electric field reverses, the particle will move downward rather than through the nanopore, see Video S2 in the Supplementary Materials). On the other hand, the electrostatic interaction was generated when the electric fields around the particles and the nanopore overlapped. When the particle approached the nanopore from below, the attractive electrostatic force between the particles and the nanopore promoted the translation of the particle due to the negative particle surface charge density. Then, when the particle entirely exited the nanopore, the attractive electrostatic force hindered the particle translation, slowing its transport velocity. When the particle was located in the midsection of the nanopore, the distribution of electric field around the cylindrical particle was basically symmetrical, and the electrostatic force on the particles was offset, so that the electrostatic contribution was less pronounced. It is worth mentioning that even if the particles themselves are not charged, if an initial angle between particle axis and the electric field is preset, the particles will rotate rapidly to the direction consistent with the electric field before passing through the nanopore, which can be attributed to the action of the dielectrophoresis force (see Videos S3 and S4 in the Supplementary Materials).
Figure 5 shows the concentration distribution and the streamlines when the particle was located at yp = −75, yp = 0, yp = 75 bp. The color variance in this figure implies the difference between the cation and anion concentration, c1–c2. From Figure 5, it can be seen that c1–c2 was almost close to 0 in the nanocavity solution, but there were a lot of cations gathering on the surface of the cylindrical particles due to the predetermined negative charge density on the particle surface. At the same time, the fluid around the particles moved upwards, and when the particles passed through the nanopore, symmetrical vortices were formed on both sides around the particles due to the rapid translation of the particles (see Video S5 in the Supplementary Materials).
Figure 6 shows the ionic current deviation I* passing through the nanopore depending on the particle position for E = 1, E = 2, and E = 3 MV/m for a particle length L = 200 bp. As the particle started to move, the ionic current deviation I* was almost absent. Considering that dimensionless I* is defined as the difference between the ionic current and the base ionic current in the nanopore, I* = 0 implies that the current in the nanopore equals the base ion current. It then decreased to be negative as the particle approached the nanopore. When the particle moved to yp = −100 bp, the ionic current deviation I* reached its minimum value, which meant that the current in the nanopore was far less than the base ion current. In other words, the current blockade effect occurred. With the continuous translation of the particle in the nanopore, the ionic current deviation I* increased gradually. When the particle moved to yp = 117 bp, I* reached its maximum value, which meant that the current in the nanopore exceeded the base ion current. The current enhancement effect occurred, after which the ionic current decreased gradually. When the particles were far away from the nanopore (yp = 300 bp), the ionic current deviation I* tended to recover zero value, and the current in the nanopore basically recovered to the base ion current. With the increase in the applied electric field, both the ionic current blockade and the enhancement effect were amplified significantly.
When the particle began to enter the nanopore, the ion transport was hindered due to the presence of particles, resulting in the ionic current decrease in the nanopore from the base ionic current, thus the ion current blockade effect. When the particle center arrived at yp = −100 bp, about half of the particle in the nanopore, the degree of the ion blockage reached its maximum. However, with the particle passing through the nanopore, the ionic current in the nanopore was increasing. Because a large number of cations carried on the particle surface flowed out of the nanopore, the ionic current recovered rapidly and even exceeded the base ion current, resulting in the current enhancement effect. In addition, the increase in fluid velocity due to particles passing through the nanopore also increased the ionic current to some extent (see Videos S5 and S6 in the Supplementary Materials). When yp > 117 bp, the entire particle exited the nanopore, the ionic current in the nanopore gradually recovered the base ionic current.
With the increase in the applied electric field, stronger electric fields were induced around the particle and nanopore, which prevented more ions in solution from passing through the nanopores. As a result, the ion blockage effect was promoted. At the same time, more cations were enriched on the surface of the particles under the stronger electric field. As the particles exited the nanopore, the enhancement effect of the ion current was also magnified accordingly.
In summary, both the ionic current blockade and the enhancement effects occurred for the experimental parameters employed in this research. The electric field intensity is an important parameter for adjusting the particle transport velocity and the current fluctuation amplitude. Therefore, if the electric potential or electric field intensity is increased in the experiment, the perforation time of particles could be decreased, thus the transmission efficiency is elevated, which is conducive to the realization of high-throughput sequencing of particles. The increase in the electric field intensity could also significantly increase the fluctuation range of the ionic current, which would reduce the requirements for the sensitivity of current detection instruments, accommodating improvements in the detectability of ionic current.

3.2. The Effect of Cylindrical Particle Length

Since the electrostatic interaction between the nanopore and the particle is significant only when the particle is located inside the nanopore, the particle length Lp is another parameter of importance to study in order to further explain the detailed mechanism of the particle dynamics. Figure 7 shows the relationship between transport velocity and the particle position for particle length Lp = 50, Lp = 100, Lp = 150, Lp = 200 bp under the applied electric field intensity of E = 2 MV/m. It can be seen that when the particle was relatively short (Lp = 50 and Lp = 100 bp), the region for constant transport velocity before the nanopore was further extended toward the nanopore. The rapid velocity increase at the entrance of the nanopore culminated in a peak value at a location closer to the midsection near yp = 0 as opposed to the premature peak location shown in Figure 3. The peak value increased with the decrease in the particle length. The velocity decay beyond the midsection was more rapid, and so the recovery of the initial velocity seemed to occur nearer to the nanopore. With the increase in the particle length (Lp = 150 bp and Lp = 200 bp), the maximum transport velocity of particles was reduced, and the particle position corresponding to the maximum transport velocity shifted toward the nanopore entrance. The monotonic acceleration and deceleration shown for shorter particles were replaced by more complicated fluctuations in the particle velocity inside the nanopore.
The above velocity change was also attributed to the electrostatic interaction between the nanopore and the particle for different lengths. When the particle is shorter, a larger part of it is dominated by the electric field around the nanopore. When the particle length is longer, only a small part of it is affected by the electric field around the nanopore. As stated before, when the particle is located in the midsection of the nanopore, the electric field distribution is symmetrical with respect to the cylindrical particle, and the electrostatic force acting on the particles is basically offset, so that the particle’s velocity remains basically unchanged. With the increase in the particle length, its center position as the particle enters the nanopore is pushed downward leading to that fact the particle’s position corresponding to the maximum transport velocity is also pushed downward and the acceleration distance and the maximum transport velocity of particles is reduced.
Figure 8 shows the ionic current deviation I* passing through the nanopore depending on the particle position for electric field E = 2 MV/m and particle length Lp = 50, Lp = 100, Lp = 150, and Lp = 200 bp. It can be seen that both the ionic current blockage and enhancement effects occurred when these cylindrical particles passed through the nanopore, which was consistent with previous research [24,26]. With the increase in the particle length, the phenomenon of the ionic current blockage was amplified, whereas the ionic current enhancement effect had little change. At the same time, the particle position corresponding to the ion current blockage was pushed downward, while the particle position corresponding to the ion current enhancement was pushed upward. In addition, when particles with different lengths reached the center of the nanopore yp = 0, almost identical ion currents passing through the nanopore were obtained.
With the increase in the particle length, more space was occupied by the particle in the nanopore, and the degree of blockade of ion transport was aggravated. The effect of the ionic blockade thus was magnified. Moreover, the particle center position yp shifted downward, as the upper end of the particle entered the nanopore, which led to the fact that the yp corresponding to ion current blockage effect was pushed downward accordingly. Similarly, the yp shifted upward, as the lower end of the particle exited the nanopore, and the yp corresponding to the ion current enhancement effect was pushed upward. However, when the particle was located in the center of the nanopore, the degree of blockade of ion transport had little difference for particles with different lengths, which made the ion currents passing through the nanopore remain almost the same.
In summary, the increase in the particle length gives rise to a decrease in the peak value of particle translation velocity. The positions of the ion current peak values corresponding to particles with different lengths are also different, which is helpful in distinguishing particles with different lengths by detecting the ionic current peak values and their corresponding residence time.

4. Conclusions

In this study, fully coupled Poisson−Nernst−Planck and Navier−Stokes equations are solved by an arbitrary Lagrangian−Eulerian (ALE) method to investigate the transport mechanism of cylindrical particles in a nanopore. In particular, the influence of an applied electric field and the particle length on the transport of cylindrical particles in nanopores is studied. The results show that the acceleration and deceleration of the cylindrical particles inside the entrance and exit regions of the nanopore become less pronounced with the increase in the particle length. The translation velocity for long particles is maintained relatively uniform inside the nanopore. The ionic current blockade effect occurs as particles enter the nanopore, and the ionic current enhancement effect occurs as they exit the nanopore for the parameter ranges studied. As an important parameter to control the particle transport velocity and current fluctuation amplitude, the electric field intensity can be used to adjust the particle transmission efficiency and ion current detectability. The increase in particle length reduces the peak value of particle velocity and changes the particle position corresponding to the peak value of the ion current, which contributes to distinguishing particles with different lengths. The results of this study deepen the understanding of the transport mechanism of long cylindrical particles in nanopores, with potential applications to bio-macromolecule detection technology.

Supplementary Materials

Supplementary materials can be found at https://www.mdpi.com/2072-666X/11/8/722/s1. Video S1: The electric field fluctuation as the particle passes through nanopore; Video S2: The particle moves downward at the case of the reversed electric field; Video S3: The particle rotates clockwise to the direction of the electric field (It should be noted that the particle itself is not charged, the round cap of the particle is removed to facilitate the drawing); Video S4: The particle rotates counter-clockwise to the direction of the electric field. (It should be noted that the particle itself is not charged, the round cap of the particle is removed to facilitate the drawing); Video S5: The fluid velocity field fluctuation as the particle passes through nanopore; Video S6: The fluid pressure field fluctuation as the particle passes through nanopore.

Author Contributions

J.G. contributed to mathematical model construction and validation, X.H. visualized the results, L.S. and T.Z. formally analyzed the results and wrote the manuscript. T.L. and S.W.J. supervised and edited the paper. All authors have read and agreed to the published version of the manuscript.

Funding

This work is funded by Grant NRF-2018R1A2B3001246 of the National Research Foundation of Korea.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Dulbecco, R. A turning point in cancer research: Sequencing the human genome. Science 1986, 231, 1055–1056. [Google Scholar] [CrossRef]
  2. Storm, A.J.; Storm, C.; Chen, J.; Zandbergen, H.; Joanny, J.F.O.; Dekker, C. Fast DNA Translocation through a Solid-State Nanopore. Nano Lett. 2005, 5, 1193–1197. [Google Scholar] [CrossRef] [Green Version]
  3. Meller, A.; Nivon, L.; Branton, D. Voltage-Driven DNA Translocations through a Nanopore. Phys. Rev. Lett. 2001, 86, 3435–3438. [Google Scholar] [CrossRef] [Green Version]
  4. Dekker, C. Solid-state nanopores. Nat. Nanotechnol. 2007, 2, 209–215. [Google Scholar] [CrossRef]
  5. Vlassiouk, I.; Krasnoslobodtsev, A.; Smirnov, S.; Germann, M. “Direct” Detection and Separation of DNA Using Nanoporous Alumina Filters. Langmuir 2004, 20, 9913–9915. [Google Scholar] [CrossRef]
  6. McNally, B.; Singer, A.; Yu, Z.L.; Sun, Y.J.; Weng, Z.P.; Meller, A. Optical Recognition of Converted DNA Nucleotides for Single-Molecule DNA Sequencing Using Nanopore Arrays. Nano Lett. 2010, 10, 2237–2244. [Google Scholar] [CrossRef] [Green Version]
  7. Fologea, D.; Gershow, M.; Ledden, B.; McNabb, D.S.; Golovchenko, J.A.; Li, J. Detecting single stranded DNA with a solid state nanopore. Nano Lett. 2005, 5, 1905–1909. [Google Scholar] [CrossRef] [Green Version]
  8. Chang, H.; Kosari, F.; Andreadakis, G.; Alam, M.A.; Vasmatzis, G.; Bashir, R. DNA-Mediated Fluctuations in Ionic Current through Silicon Oxide Nanopore Channels. Nano Lett. 2004, 4, 1551–1556. [Google Scholar] [CrossRef] [Green Version]
  9. Tsutsui, M.; Taniguchi, M.; Yokota, K.; Kawai, T. Identifying single nucleotides by tunnelling current. Nat. Nanotechnol. 2010, 5, 286–290. [Google Scholar] [CrossRef] [PubMed]
  10. Ivanov, A.P.; Instuli, E.; McGilvery, C.M.; Baldwin, G.; McComb, D.W.; Albrecht, T.; Edel, J.B. DNA Tunneling Detector Embedded in a Nanopore. Nano Lett. 2011, 11, 279–285. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. Xie, P.; Xiong, Q.H.; Fang, Y.; Qing, Q.; Lieber, C.M. Local electrical potential detection of DNA by nanowire-nanopore sensors. Nat. Nanotechnol. 2012, 7, 119–125. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Feliciano, G.T.; Sanz-Navarro, C.; Coutinho-Neto, M.D.; Ordejon, P.; Scheicher, R.H.; Rocha, A.R. Capacitive DNA Detection Driven by Electronic Charge Fluctuations in a Graphene Nanopore. Phys. Rev. Appl. 2015, 3, 034003. [Google Scholar] [CrossRef]
  13. Kasianowicz, J.J.; Brandin, E.; Branton, D. Characterization of individual polynucleotide molecules using a membrane channel. Proc. Natl. Acad. Sci. USA 1996, 93, 13770–13773. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Guasch, A.; Pous, J.; Ibarra, B.; Gomis-Rüth, F.X.; Valpuesta, J.M.a.; Sousa, N.; Carrascosa, J.L.; Coll, M. Detailed architecture of a DNA translocating machine: The high-resolution structure of the bacteriophage φ29 connector particle11Edited by R. Huber. J. Mol. Biol. 2002, 315, 663–676. [Google Scholar] [CrossRef]
  15. Cao, C.; Ying, Y.L.; Hu, Z.L.; Liao, D.F.; Tian, H.; Long, Y.T. Discrimination of oligonucleotides of different lengths with a wild-type aerolysin nanopore. Nat. Nanotechnol. 2016, 11, 713–718. [Google Scholar] [CrossRef]
  16. Li, J.; Stein, D.; McMullan, C.; Branton, D.; Aziz, M.J.; Golovchenko, J.A. Ion-beam sculpting at nanometre length scales. Nature 2001, 412, 166–169. [Google Scholar] [CrossRef] [Green Version]
  17. Storm, A.J.; Chen, J.H.; Ling, X.S.; Zandbergen, H.W.; Dekker, C. Fabrication of solid-state nanopores with single-nanometre precision. Nat. Mater. 2003, 2, 537–540. [Google Scholar] [CrossRef]
  18. Venkatesan, B.M.; Dorvel, B.; Yemenicioglu, S.; Watkins, N.; Petrov, I.; Bashir, R. Highly Sensitive, Mechanically Stable Nanopore Sensors for DNA Analysis. Adv. Mater. 2009, 21, 2771. [Google Scholar] [CrossRef]
  19. Schneider, G.F.; Kowalczyk, S.W.; Calado, V.E.; Pandraud, G.; Zandbergen, H.W.; Vandersypen, L.M.K.; Dekker, C. DNA Translocation through Graphene Nanopores. Nano Lett. 2010, 10, 3163–3167. [Google Scholar] [CrossRef] [Green Version]
  20. Liu, K.; Feng, J.; Kis, A.; Radenovic, A. Atomically Thin Molybdenum Disulfide Nanopores with High Sensitivity for DNA Translocation. ACS Nano 2014, 8, 2504–2511. [Google Scholar] [CrossRef]
  21. Sze, J.Y.Y.; Ivanov, A.P.; Cass, A.E.G.; Edel, J.B. Single molecule multiplexed nanopore protein screening in human serum using aptamer modified DNA carriers. Nat. Commun. 2017, 8, 1552. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Steinbock, L.J.; Otto, O.; Chimerel, C.; Gornall, J.; Keyser, U.F. Detecting DNA Folding with Nanocapillaries. Nano Lett. 2010, 10, 2493–2497. [Google Scholar] [CrossRef]
  23. Qian, S.; Joo, S.W.; Hou, W.S.; Zhao, X. Electrophoretic motion of a spherical particle with a symmetric nonuniform surface charge distribution in a nanotube. Langmuir 2008, 24, 5332–5340. [Google Scholar] [CrossRef] [PubMed]
  24. Zhang, M.K.; Yeh, L.H.; Qian, S.Z.; Hsu, J.P.; Joo, S.W. DNA Electrokinetic Translocation through a Nanopore: Local Permittivity Environment Effect. J. Phys. Chem. C 2012, 116, 4793–4801. [Google Scholar] [CrossRef]
  25. Zhang, M.K.; Ai, Y.; Sharma, A.; Joo, S.W.; Kim, D.S.; Qian, S.Z. Electrokinetic particle translocation through a nanopore containing a floating electrode. Electrophoresis 2011, 32, 1864–1874. [Google Scholar] [CrossRef] [PubMed]
  26. Ai, Y.; Qian, S.Z. Electrokinetic particle translocation through a nanopore. Phys. Chem. Chem. Phys. 2011, 13, 4060–4071. [Google Scholar] [CrossRef]
  27. Ai, Y.; Qian, S. Direct numerical simulation of electrokinetic translocation of a cylindrical particle through a nanopore using a Poisson-Boltzmann approach. Electrophoresis 2011, 32, 996–1005. [Google Scholar] [CrossRef]
  28. Shi, C.C.; Kong, Z.; Sun, T.Y.; Liang, L.J.; Shen, J.W.; Zhao, Z.Y.; Wang, Q.; Kang, Z.Z.; Agren, H.; Tu, Y.Q. Molecular dynamics simulations indicate that DNA bases using graphene nanopores can be identified by their translocation times. RSC Adv. 2015, 5, 9389–9395. [Google Scholar] [CrossRef]
  29. Li, J.P.; Zhang, Y.; Yang, J.K.; Bi, K.D.; Ni, Z.H.; Li, D.Y.; Chen, Y.F. Molecular dynamics study of DNA translocation through graphene nanopores. Phys. Rev. E 2013, 87, 062707. [Google Scholar] [CrossRef]
  30. Bhattacharya, S.; Derrington, I.M.; Pavlenok, M.; Niederweis, M.; Gundlach, J.H.; Aksimentiev, A. Molecular dynamics study of MspA arginine mutants predicts slow DNA translocations and ion current blockades indicative of DNA sequence. ACS Nano 2012, 6, 6960–6968. [Google Scholar] [CrossRef] [Green Version]
  31. Comer, J.; Aksimentiev, A. Predicting the DNA Sequence Dependence of Nanopore Ion Current Using Atomic-Resolution Brownian Dynamics. J. Phys. Chem. C 2012, 116, 3376–3393. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Hulings, Z.K.; Melnikov, D.V.; Gracheva, M.E. Brownian dynamics simulations of the ionic current traces for a neutral nanoparticle translocating through a nanopore. Nanotechnology 2018, 29, 445204. [Google Scholar] [CrossRef] [PubMed]
  33. Zhou, T.; Xu, Y.F.; Liu, Z.Y.; Joo, S.W. An Enhanced One-Layer Passive Microfluidic Mixer with an Optimized Lateral Structure with the Dean Effect. J. Fluids Eng.—Trans. ASME 2015, 137, 091102. [Google Scholar] [CrossRef]
  34. Zhou, T.; Liu, T.; Deng, Y.B.; Chen, L.M.; Qian, S.Z.; Liu, Z.Y. Design of microfluidic channel networks with specified output flow rates using the CFD-based optimization method. Microfluid. Nanofluidics 2017, 21, 11. [Google Scholar] [CrossRef]
  35. Zhou, T.; Ji, X.; Shi, L.; Zhang, X.; Song, Y.; Joo, S.W. AC dielectrophoretic deformable particle-particle interactions and their relative motions. Electrophoresis 2019. [Google Scholar] [CrossRef] [PubMed]
  36. Zhou, T.; Ge, J.; Shi, L.; Fan, J.; Liu, Z.; Woo Joo, S. Dielectrophoretic choking phenomenon of a deformable particle in a converging-diverging microchannel. Electrophoresis 2018, 39, 590–596. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Schematics of a cylindrical particle passing through a nanopore.
Figure 1. Schematics of a cylindrical particle passing through a nanopore.
Micromachines 11 00722 g001
Figure 2. The y-component translational velocity u* of the cylinder particle as a function of the particle’s location y*p under E = 20 KV/m, the ratio of the particle radius to the Debye length, ka = 1.03 (red line and circles), 0.46 (black lines and squares). Symbols and lines represent, respectively, the numerical solution of reference [26] and the numerical results from the present model.
Figure 2. The y-component translational velocity u* of the cylinder particle as a function of the particle’s location y*p under E = 20 KV/m, the ratio of the particle radius to the Debye length, ka = 1.03 (red line and circles), 0.46 (black lines and squares). Symbols and lines represent, respectively, the numerical solution of reference [26] and the numerical results from the present model.
Micromachines 11 00722 g002
Figure 3. The particle translation velocity depending on particle position under the applied electric field strength of E = 1 MV/m (black line with squares), E = 2 MV/m (red line with circles), and E = 3 MV/m (blue line with triangles) for particle length Lp = 200 bp. The vertical dotted lines indicate the locations yp = −75 and yp = 75 bp.
Figure 3. The particle translation velocity depending on particle position under the applied electric field strength of E = 1 MV/m (black line with squares), E = 2 MV/m (red line with circles), and E = 3 MV/m (blue line with triangles) for particle length Lp = 200 bp. The vertical dotted lines indicate the locations yp = −75 and yp = 75 bp.
Micromachines 11 00722 g003
Figure 4. The electric field distribution around the nanopore for particle length Lp = 200 bp and E = 2 MV/m when the particle center is located at (a) yp = −75, (b) yp = 0, (c) yp = 75 bp. The darker colors in the diagram represent a stronger electric field, and the arrows indicate the direction of the electric field.
Figure 4. The electric field distribution around the nanopore for particle length Lp = 200 bp and E = 2 MV/m when the particle center is located at (a) yp = −75, (b) yp = 0, (c) yp = 75 bp. The darker colors in the diagram represent a stronger electric field, and the arrows indicate the direction of the electric field.
Micromachines 11 00722 g004
Figure 5. The distribution of the concentration difference c1−c2 between cation and anion and streamline diagram when the particle length is Lp = 200 bp, the particle locates at (a) yp = −75, (b) yp = 0, (c) yp = 75 bp. The darker colors in the diagram mean the larger value of c1−c2, and the lines with arrows denote the streamlines of the flow field.
Figure 5. The distribution of the concentration difference c1−c2 between cation and anion and streamline diagram when the particle length is Lp = 200 bp, the particle locates at (a) yp = −75, (b) yp = 0, (c) yp = 75 bp. The darker colors in the diagram mean the larger value of c1−c2, and the lines with arrows denote the streamlines of the flow field.
Micromachines 11 00722 g005
Figure 6. The ion current deviation I* passing through the nanopore depending on the particle position under E = 1 (black line plus frame), E = 2 (red line plus circle), E = 3 MV/m (blue line plus triangle) when the particle length is Lp = 200 bp. Dotted lines indicate the particle position yp = −100 and yp = 117 bp, and are not meant to indicate locations for maximum or minimum current.
Figure 6. The ion current deviation I* passing through the nanopore depending on the particle position under E = 1 (black line plus frame), E = 2 (red line plus circle), E = 3 MV/m (blue line plus triangle) when the particle length is Lp = 200 bp. Dotted lines indicate the particle position yp = −100 and yp = 117 bp, and are not meant to indicate locations for maximum or minimum current.
Micromachines 11 00722 g006
Figure 7. The particle translation velocity vs. particle position for E = 2 MV/m and particle length Lp = 50 (black line plus frame), Lp = 100 (red line plus circle), Lp = 150 (blue line plus triangle), and Lp = 200 bp (green line plus diamond). The vertical dotted line indicates the center of the nanopore.
Figure 7. The particle translation velocity vs. particle position for E = 2 MV/m and particle length Lp = 50 (black line plus frame), Lp = 100 (red line plus circle), Lp = 150 (blue line plus triangle), and Lp = 200 bp (green line plus diamond). The vertical dotted line indicates the center of the nanopore.
Micromachines 11 00722 g007
Figure 8. The ion current deviation I* passing through the nanopore depending on the particle position when the applied electric field strength E = 2 MV/m, the particle length is Lp = 50 (black line plus frame), Lp = 100 (red line plus circle), Lp = 150 (blue line plus triangle), and Lp = 200 bp (green line plus diamond). The dotted line indicates that the particles reach the center of the nanopore, yp = 0.
Figure 8. The ion current deviation I* passing through the nanopore depending on the particle position when the applied electric field strength E = 2 MV/m, the particle length is Lp = 50 (black line plus frame), Lp = 100 (red line plus circle), Lp = 150 (blue line plus triangle), and Lp = 200 bp (green line plus diamond). The dotted line indicates that the particles reach the center of the nanopore, yp = 0.
Micromachines 11 00722 g008

Share and Cite

MDPI and ACS Style

Shi, L.; He, X.; Ge, J.; Zhou, T.; Li, T.; Joo, S.W. The Influence of Electric Field Intensity and Particle Length on the Electrokinetic Transport of Cylindrical Particles Passing through Nanopore. Micromachines 2020, 11, 722. https://doi.org/10.3390/mi11080722

AMA Style

Shi L, He X, Ge J, Zhou T, Li T, Joo SW. The Influence of Electric Field Intensity and Particle Length on the Electrokinetic Transport of Cylindrical Particles Passing through Nanopore. Micromachines. 2020; 11(8):722. https://doi.org/10.3390/mi11080722

Chicago/Turabian Style

Shi, Liuyong, Xiaohan He, Jian Ge, Teng Zhou, Ting Li, and Sang Woo Joo. 2020. "The Influence of Electric Field Intensity and Particle Length on the Electrokinetic Transport of Cylindrical Particles Passing through Nanopore" Micromachines 11, no. 8: 722. https://doi.org/10.3390/mi11080722

APA Style

Shi, L., He, X., Ge, J., Zhou, T., Li, T., & Joo, S. W. (2020). The Influence of Electric Field Intensity and Particle Length on the Electrokinetic Transport of Cylindrical Particles Passing through Nanopore. Micromachines, 11(8), 722. https://doi.org/10.3390/mi11080722

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop