Next Article in Journal
Rapamycin-Induced Apoptosis in HGF-Stimulated Lens Epithelial Cells by AKT/mTOR, ERK and JAK2/STAT3 Pathways
Next Article in Special Issue
What Controls DNA Looping?
Previous Article in Journal
The Inflammation-Related Gene S100A12 Is Positively Regulated by C/EBPβ and AP-1 in Pigs
Previous Article in Special Issue
DNA Break Mapping Reveals Topoisomerase II Activity Genome-Wide
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Theoretical Study of the Transpore Velocity Control of Single-Stranded DNA

Department of Mechanical Science and Bioengineering, Graduate School of Engineering Science, Osaka University, 1-3 Machikaneyama, Toyonaka, Osaka 560-8531, Japan
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2014, 15(8), 13817-13832; https://doi.org/10.3390/ijms150813817
Submission received: 10 June 2014 / Revised: 15 July 2014 / Accepted: 22 July 2014 / Published: 11 August 2014
(This article belongs to the Special Issue Identification and Roles of the Structure of DNA)

Abstract

:
The electrokinetic transport dynamics of deoxyribonucleic acid (DNA) molecules have recently attracted significant attention in various fields of research. Our group is interested in the detailed examination of the behavior of DNA when confined in micro/nanofluidic channels. In the present study, the translocation mechanism of a DNA-like polymer chain in a nanofluidic channel was investigated using Langevin dynamics simulations. A coarse-grained bead-spring model was developed to simulate the dynamics of a long polymer chain passing through a rectangular cross-section nanopore embedded in a nanochannel, under the influence of a nonuniform electric field. Varying the cross-sectional area of the nanopore was found to allow optimization of the translocation process through modification of the electric field in the flow channel, since a drastic drop in the electric potential at the nanopore was induced by changing the cross-section. Furthermore, the configuration of the polymer chain in the nanopore was observed to determine its translocation velocity. The competition between the strength of the electric field and confinement in the small pore produces various transport mechanisms and the results of this study thus represent a means of optimizing the design of nanofluidic devices for single molecule detection.

1. Introduction

The high-speed reading of deoxyribonucleic acid (DNA) sequences is an important means of elucidating complete genetic sequences, and may enable the development of new medical treatments [1,2]. Recently, novel DNA and ribonucleic acid (RNA) sequencing technologies have been developed. Among these, nanopore sequencing devices are one of the most significant issues and represent an emerging non-optical process for high-throughput single-molecule detection [1,2,3,4], in which individual nucleobases are identified by measuring transpore ionic current blockade [5,6,7] or transverse tunneling current [8,9,10,11] during the transport of single-stranded DNA (ssDNA) through a nanometer-sized gap. Understanding biological polymer transport phenomena is a crucial issue in the development of DNA sequencing techniques, as well as in the study of many of the physical properties of polymers [12], and both the theoretical [13,14,15,16,17] and experimental [18,19,20,21,22,23,24,25,26,27,28] aspects of polymer translocation through nanopores have been widely studied. Computational studies have provided particularly valuable insights into the physics of transport within confined micro/nanochannels and previous works have examined the variation of translocation time with polymer chain length [13,14,15,19,21,29,30,31,15,19,21,29], pore dimensions [31], driving force [15,19,21,31], sequences and secondary structures [21,22,32], polymer–pore interactions [21,22,23,33], and polymer configurations [15,34]. Table 1 lists the various nanopore devices and polymers used in the pioneering research studies investigating these subjects with the aim of achieving an advanced DNA sequencer.
Table 1. Nanopore devices and deoxyribonucleic acid (DNA)/ribonucleic acid (RNA) samples used in published experimental studies.
Table 1. Nanopore devices and deoxyribonucleic acid (DNA)/ribonucleic acid (RNA) samples used in published experimental studies.
GroupPore TypeDiameter (nm)Length (nm)Voltage (mV)Voltage/Length (×106 V/m)Polymer Length (bp or nt)Polymer Type
Kasianowicz et al. [5]α-HL1.35.212023150ssDNA, ssRNA *
Meller et al. [18,19,20]α-HL1.35.250–3009.6–585–100ssDNA
Butler et al. [35]MspA110140, 18014, 1850ssDNA
Wendell et al. [36]Phi293.37.540, 755.3, 105.5 kdsDNA
Franceschini et al. [37]ClyA7.8131007.7290 bp, 51 ntdsDNA, ssDNA
Li et al. [21]SiN3, 105–1060, 1206–243–10 kdsDNA
Storm et al. [22]SiN1020100–6005–3010–97 kdsDNA
Skinner et al. [24]SiN1020100–6005–3010–30 kdsDNA, dsRNA, ssRNA *
Tsutsui et al. [26]SiN502001000548.5 kdsDNA
Fologea et al. [27,28]SiN1010, 2801200.43, 123 kdsDNA, ssDNA
Schneider et al. [25]Graphene220.320067048.5 kdsDNA
* ssRNA denotes poly(A), poly(C), and poly(U).
Sung and Park [13] and Muthukumar [14] studied the passage of single polymer molecules through the pores of a membrane during diffusion across a free energy barrier due to chemical potential differences. Both groups modeled the stochastic processes associated with the transport of long polymers based on the Fokker–Planck equation and were able to predict a scaling law describing translocation time, τ, as a function of polymer length, N. Storm et al. [22] and Skinner et al. [24] investigated the translocation of double-stranded DNA (dsDNA) through silicon nitride (SiN) nanopores that were 10 nm in diameter and 30 nm thick. They also identified that a power law best described the relationship between τ and the polymer length, such that τ~N1.27. The use of ultrathin nanopores (0.3 nm thick) fabricated within a graphene monolayer [25] is known to result in a slightly larger τ value than that obtained using SiN nanopores [21,22] and it has been suggested that these small pores as well as interactions with the graphene result in the slower translocation. This phenomenon has also been investigated on the basis of Langevin dynamics simulations [31,34]. In other works, Meller et al. [19] studied the translocation of ssDNA through a biological α-hemolysin (α-HL) nanopore and determined that the translocation velocity of short polymers exhibited a significant dependence on the length of the polymer, whereas there was no dependence in the case of long polymers. The engineered Mycobacterium smegmatis porin A (MspA) [35] and phi29 [36] protein nanopores were found to allow the translocation of ssDNA and dsDNA with remarkable stability against environmental stresses. It is indicated that an engineered DNA transporter is able to recognize and chaperone the specific DNA molecule across a biological membrane, making a further step for the application of a nanofluidic platform [37]. It was also found that, during forced translocation in narrow pores, the scaling exponents depended on the driving force, F, based on the relationship τ~F−1 [15,30,31]. Although the hydrodynamic effects on polymer chains appear to account for part of the force counteracting external forces [13,29], these effects seem to make only a minor contribution to the transport of DNA, since it has a large number of charges and small surface areas in comparison to other polymer particles [38,39,40]. In particular, the electrokinetic transport of DNA passing through very narrow spaces is predominantly affected by collisions with channel walls [41].
We are interested in a long polymer translocation mechanism in micro/nanochannels and nanopores [42,43,44] under the effects of nonuniform electric fields, since such mechanisms have not yet been sufficiently elucidated. In the present study, we attempt to gain a better understanding of the translocation mechanism of a DNA-like polymer chain, equivalent to 48 × 103 nucleotides (48 knt), penetrating a solid-sate nanopore in the presence of nonuniform electric fields, as illustrated in Figure 1. The cross-section of the nanopore is expected to play an important role in terms of controlling the translocation process. While nanopores embedded in nanochannels were supposed to be effective to slow down the transpore velocity of ssDNA [42,44], the mechanism has remained to be clarified. Herein, focusing on the multiply-connected nanofluidic channels, the retardation process and its advantage are discussed from a theoretical point of view. As part of this work, we develop a coarse-grained ssDNA model [16,45,46] and perform Langevin dynamics simulations of ssDNA transport under nonuniform electric fields in a rectangular nanochannel containing a nanopore with various cross-sections [47,48,49,50,51], where the electric fields are calculated for the cross-sections, ranging from 20 × 20 to 50 × 50 nm2. The results allow a visual analysis of the electrokinetic transport dynamics of ssDNA chains and allow us to determine the most suitable morphology for nanofluidic flow channels for single molecule detection. Furthermore, the simulation results are clearly understood by a theoretical model in the framework of the Langevin equation. Consequently, a relationship among the electrokinetic transport of ssDNA, pore dimensions, and multiply-connected structures of the nanofluidic channel are clarified and a desirable design to control the translocation velocity is concluded.
Figure 1. (a) Schematic illustration of a rectangular nanochannel used in Langevin dynamics simulations, in which a nanopore is embedded in the nanochannel. In the three-dimensional model, the width (w) and height (h) of a 200 nm long nanopore are varied as simulation parameters in a nanochannel of 2200 nm in length, 500 nm in width, and h in height. The center of mass of single-stranded DNA (ssDNA) is initially located at a distance of the radius of gyration (Rg = 300 nm) away from the nanochannel entrance where the x- and y-coordinate of the mass center are in coincidence with the center of nanochannel; and (b) An illustration showing the coarse-grained bead-spring model of ssDNA in the simulations to assess the electrokinetic transport dynamics and to optimize the structure of the nanofluidic channel for single molecule sensing.
Figure 1. (a) Schematic illustration of a rectangular nanochannel used in Langevin dynamics simulations, in which a nanopore is embedded in the nanochannel. In the three-dimensional model, the width (w) and height (h) of a 200 nm long nanopore are varied as simulation parameters in a nanochannel of 2200 nm in length, 500 nm in width, and h in height. The center of mass of single-stranded DNA (ssDNA) is initially located at a distance of the radius of gyration (Rg = 300 nm) away from the nanochannel entrance where the x- and y-coordinate of the mass center are in coincidence with the center of nanochannel; and (b) An illustration showing the coarse-grained bead-spring model of ssDNA in the simulations to assess the electrokinetic transport dynamics and to optimize the structure of the nanofluidic channel for single molecule sensing.
Ijms 15 13817 g001

2. Results and Discussion

2.1. Validation of the Coarse-Grained Single-Stranded DNA (ssDNA) Model

As a result of the Langevin dynamics simulation, Figure 2a shows the mean square displacement of the center of mass of the ssDNA model as a function of time. This plot represents the average of results from 90 simulations at each data point and is clearly linear. The associated diffusion coefficient, D, can be obtained according to the Einstein relation and is calculated to be 2.25 × 10−12 m2/s. Figure 2b presents a plot of the distance from the origin to the center of mass under an external electric field of 1 × 106 V/m, in which each data point is the average of 270 simulations. The electrophoretic mobility calculated from the ratio of the velocity values to the applied electric field strengths is μ = 2.87 × 10−8 m2/Vs. These two results are in close agreement with the experimental data described in the methodology (Section 3) [52]. From the viewpoint of diffusivity and electrophoretic mobility, the present parameter set is therefore acceptable when assessing the electrokinetic transport of ssDNA.
Figure 2. (a) Mean square displacement of the bead-spring model resulting from 90 simulation runs and (b) displacement of the mass center of the bead-spring chain under an electric field of 1 × 106 V/m obtained from 270 simulation runs. Each result is well fitted with straight lines by the least-squares method. The slope of plot (a) corresponds to the diffusion coefficient while that of (b) is the velocity that translates to the electrophoretic mobility.
Figure 2. (a) Mean square displacement of the bead-spring model resulting from 90 simulation runs and (b) displacement of the mass center of the bead-spring chain under an electric field of 1 × 106 V/m obtained from 270 simulation runs. Each result is well fitted with straight lines by the least-squares method. The slope of plot (a) corresponds to the diffusion coefficient while that of (b) is the velocity that translates to the electrophoretic mobility.
Ijms 15 13817 g002

2.2. Langevin Dynamics Simulations

As shown in Figure 3, electrostatic potentials across the microchannel, nanochannel, and nanopore are determined from the finite element method (FEM) analysis [53], in which the potential curves extracted along the central axis are presented for several nanopore cross-sections. It is found that the slope of the electrostatic potential becomes steeper in the narrower channels, as shown in Figure 3a,b. A large drop in the potential at the nanopore suppresses the potential difference outside the nanopore. As can be seen from Figure 3c, the electric field strength increases as the cross-sectional area of the nanopore is reduced. The electric field strengths calculated along the central axis of the nanochannel and nanopore with various cross-sections are also summarized in Table 2. In the previous experimental studies as summarized in Table 1, as well as in numerical analyses, other researchers have also found that strong electric fields are associated with nanopores [54]. Some publications note that such a strong electric field tends to be proportional to the value of (dpore/deff)2, where dpore is the diameter of the pore and deff is the effective diameter outside the pore [54,55]. Our computational results also agree with the potential drop resulting from variations in the nanopore cross-section.
Table 2. Electric field and ssDNA transpore properties in the nanochannel and nanopore.
Table 2. Electric field and ssDNA transpore properties in the nanochannel and nanopore.
Pore Size (nm2)Echannel (V/m)Epore (V/m)NporeΔNx (m−1)vchannel (mm/s)vpore (mm/s)
30 × 308.2 × 1042.7 × 1066.481.30 × 1092.31.3
40 × 401.2 × 1052.2 × 1066.819.64 × 1083.31.1
50 × 501.5 × 1051.9 × 1069.297.45 × 1084.31.2
Figure 3. (a) Electrostatic potential resulting from FEM analysis for the system including microchannel, nanochannel, and nanopore whose cross-section is 50 × 50 nm2; (b) Magnified view of (a) focusing onnear the nanopore;and (c) Electrostatic potentials, ϕ, for nanopores of cross-section 20 × 20, 30 × 30, 40 × 40, and 50 × 50 nm2, resulting from three-dimensional FEM analyses. The entrance and exit of the nanopores are indicated by the dotted lines.
Figure 3. (a) Electrostatic potential resulting from FEM analysis for the system including microchannel, nanochannel, and nanopore whose cross-section is 50 × 50 nm2; (b) Magnified view of (a) focusing onnear the nanopore;and (c) Electrostatic potentials, ϕ, for nanopores of cross-section 20 × 20, 30 × 30, 40 × 40, and 50 × 50 nm2, resulting from three-dimensional FEM analyses. The entrance and exit of the nanopores are indicated by the dotted lines.
Ijms 15 13817 g003
As shown in Figure 4a, we also ascertained the number of beads along a 200 nm long nanopore during the simulations. In this figure, the entire data set resulting from 20 simulation runs for a 30 × 30 nm2 cross-section nanopore is presented. At t = 0 s, a leading bead entered the pore, at which point the elapsed time was tracked until the end of the chain left the pore. The distribution of bead numbers seems to be discretized at several specific numerical values. As can be seen in the insets to this figure, which show illustrations of the nanopore, the discretized numbers correspond to specific folded structures of the polymer chain. Sufficiently uncoiled ssDNA chains tend to pass through the nanopore in an unfolded form and therefore, the translocation time is relatively long. In contrast, coiled chains adopt folded forms in the pore, resulting in shorter translocation periods. Figure 4b presents a summary of the data in Figure 4a in the form of a histogram. The highest peak in this plot corresponds to an unfolded structure, while the second and third highest peaks equate to 1- and 2-fold forms, respectively. More detailed illustrations of the unfolded, 1-fold, and 2-fold forms at the nanopore are also presented in Figure 4c–e, respectively. Figure 5a,b show the results for 40 × 40 and 50 × 50 nm2 cross-section nanopores, respectively, where four typical samples are presented by color variations. The time series data in the plots apparently fluctuate with increasing cross-sectional areas. The larger the cross-sections become, the more frequently the ssDNA will change its configuration, thus producing multifold forms in the pore. As a result, the distribution of multifold-structures increases as the cross-sectional area increases. The fitted distributions obtained from 20 simulation runs for each condition are summarized in Figure 5c. The concentrations of electric charges resulting from the folded structures increase the translocation speed due to the associated strong electric force. Although we also performed simulations for a 20 × 20 nm2 cross-section nanopore, the electric field outside a nanopore of this size was evidently too weak to introduce the ssDNA into the pore. This result implies that an excessively small pore will require a long period of time to attract charged molecules to it. A weak electric field outside the pore, as is produced in the case of an overly small pore, is therefore disadvantageous for the polymer chain to overcome the entropic barriers at interconnections in the channel, because the large difference of cross-sections requires strong force to uncoil a coiled structure to introduce it into the nanopore [56,57]. Figure 5d presents a normalized version of the distribution data in Figure 5c. Comparing the three cross-sections, it is evident that multifold-structures become prominent as the cross-sectional area increases. With respect to single-molecule detection, it is desirable to maintain unfolded configurations for as long as possible to slow down the translocation speed. Thus, the 30 × 30 nm2 cross-section pore is suggested to be the most suited to the sequential transport of ssDNA molecules.
Figure 4. (a) Time series data, where each run is presented by color variation; (b) histogram of the number of beads in a 30 × 30 nm2 cross-section nanopore, obtained from 20 runs of the Langevin dynamics simulation; A leading head enters the pore at t = 0 sand the time elapsed is recorded until the end bead leaves the pore. Illustrations of the entire ssDNA chain are also presented, showing (c) unfolded; (d) 1-fold; and (e) 2-fold forms.
Figure 4. (a) Time series data, where each run is presented by color variation; (b) histogram of the number of beads in a 30 × 30 nm2 cross-section nanopore, obtained from 20 runs of the Langevin dynamics simulation; A leading head enters the pore at t = 0 sand the time elapsed is recorded until the end bead leaves the pore. Illustrations of the entire ssDNA chain are also presented, showing (c) unfolded; (d) 1-fold; and (e) 2-fold forms.
Ijms 15 13817 g004
Figure 5. Time series data indicating the number of beads in nanopores of cross-section (a) 40 × 40 and (b) 50 × 50 nm2, in which only four typical data are presented by color variation in each case; (c) fitted distribution resulting from the complete data acquired from 20 simulations for each condition; and (d) the normalized distribution of (c).
Figure 5. Time series data indicating the number of beads in nanopores of cross-section (a) 40 × 40 and (b) 50 × 50 nm2, in which only four typical data are presented by color variation in each case; (c) fitted distribution resulting from the complete data acquired from 20 simulations for each condition; and (d) the normalized distribution of (c).
Ijms 15 13817 g005
For deep understanding of the electrokinetic transport phenomena in the nanofluidic device, the simulation results are analyzed by a theoretical model of the Langevin equation as also described in detail in the methodology section. Particularly, a relationship between the translocation velocity and the pore size attracts most of our interests. Figure 6 shows velocity profiles of the mass centers of the ssDNA for the three cases presented in Figure 5, in which xG denotes the position of mass center along the x-axis measured from the nanochannel entrance (2200 nm in total), and the nanopore is located from xG = 1000 to 1200 nm. In overdamped Langevin dynamics simulations, the velocity of a particle is directly proportional to the force on it, as theoretically described in the next section. For each cross-section, the velocity linearly increased until the mass center reached an xG value of approximately 500 nm, at which point the leading bead moved into the stronger field while the remainder of the chain was still in front of the nanochannel entrance. Therefore, the number of beads in the nanochannel increased in a stepwise fashion over time. According to the Langevin equation, the equation of motion of the mass center along the pore axis may be roughly expressed by:
Ijms 15 13817 i001
where, assuming conservative force and thermal fluctuations, we can apply:
Ijms 15 13817 i002
Ijms 15 13817 i003
The term ΔF/(NΔx) represents the ratio of the external force on the mass center to the displacement and ν0 is the initial velocity at the entrance. ζ is the friction coefficient fitted to represent the property of ssDNA and results in 4.68 × 10−12 kg/s. In Figure 6, the slopes of the plotted data in the initial portion of each graph are respectively 9.09 × 103, 9.10 × 103, and 9.80 × 103 s−1 for the 30 × 30, 40 × 40, and 50 × 50 nm2 nanopores, giving an average value of 9.33 × 103 s−1. In this region, the increment in which beads enter the nanochannel is almost constant despite the different channel cross-sections. In addition, when ΔFx is primarily due to the electric force in the nanochannel, we can write ΔF = QEchannelΔN, meaning that the change in the force is governed by the increase in the number of beads entering the channel under the almost uniform electric field. Equation (1) can then be replaced by:
Ijms 15 13817 i004
using Echannel approximated by the electric field at the center of the nanochannel as listed in Table 2. For the three cases, the values of ΔNx are 1.30 × 109, 9.64 × 108, and 7.45 × 108 m−1for the 30 × 30, 40 × 40, and 50 × 50 nm2 nanopores, respectively (Table 2). A charged bead in an Echannel field generates QEchannel such that the displacement, Δx, of the mass center related to each bead increment is proportional to Echannel and this explains why dvG/dxG is almost constant for all three cross-sections. In the following region, when the mass center approaches xG = 500 nm, there are obvious differences in velocity between the three cases. At this point, some beads are already in the nanopore. Subsequently, the velocity shows a moderate increase and appears to reach a terminal velocity when the center of mass passes through the nanopore. At this stage, the beads in the nanopore are driven forward due to the strong electric field and simultaneously experience counteracting force, being pushed back by the leading portion of the chain and pulled by the following portion. Since small nanopore cross-sections produce a strong driving force, the confinement in this region also gives rise to the counteraction including the entropic force and polymer–wall interactions. In this region, νG can be represented as:
Ijms 15 13817 i005
where Npore and Nchannel are the number of beads in the nanopore and nanochannel, respectively, and Epore is the electric field strength in the pore. Here, the ssDNA chain is usually stretched and rarely collides with the wall as it passes through the interface between the nanochannel and the nanopore, and so the counteracting force is negligibly small compared to the other terms. As shown in Figure 5d and Table 2, the average number Npore of beads in the nanopores is determined from the distributions. The remaining beads are in the nanochannel, such that Nchannel = NNpore. Using the electric field Epore at the center of the nanopore and Echannel, vG in Equation (5) results in 3.54, 4.41, and 5.49 mm/s for the 30 × 30, 40 × 40, and 50 × 50 nm2 cross-section nanopores, respectively. These theoretically derived values are in good agreement with the simulations shown in Figure 6. Particularly, in the 30 × 30 nm2pore, the rapid change in curvature of the plot occurs at an xG value of approximately 500 nm, indicating that the translocation process immediately reaches a steady state condition as the nanopore works to pump beads into the trans channel. In contrast, in the other pores, more moderate transitions of the velocity are observed and apparent transition points cannot be determined. As a result, the terminal velocities approach the theoretical values. Our data indicate that folded configurations of ssDNA chains in large cross-section pores cause moderate increases in the velocity of the mass center, and this results in high terminal velocities. In other words, our results explain why electrophoretic mobility decreases during transport through a confined space embedded in the fluidic channel [17,26,58]. From the viewpoint of molecular sequencing, increased knowledge of changes in the velocity and suppression of excessive increases in this velocity are desirable when attempting to ascertain details concerning the configuration changes of polymer molecules.
Figure 6. Velocity profile of the centers of mass of ssDNA chains passing through nanopores of cross-section 30 × 30, 40 × 40, and 50 × 50 nm2, in which xG is the position of mass center measured from the nanochannel entrance along the x-axis. The results of theoretical calculations using Equations (4) and (5) are shown as solid lines. The start and end of the nanopore are situated at xG = 1.0 and 1.2 μm and the end of the nanochannel is at xG = 2.2 μm, all of which are indicated by dashed lines.
Figure 6. Velocity profile of the centers of mass of ssDNA chains passing through nanopores of cross-section 30 × 30, 40 × 40, and 50 × 50 nm2, in which xG is the position of mass center measured from the nanochannel entrance along the x-axis. The results of theoretical calculations using Equations (4) and (5) are shown as solid lines. The start and end of the nanopore are situated at xG = 1.0 and 1.2 μm and the end of the nanochannel is at xG = 2.2 μm, all of which are indicated by dashed lines.
Ijms 15 13817 g006

3. Langevin Dynamics Simulations of Polymer Chain Motion

A Langevin dynamics simulation was applied to investigate the behavior of a polymer chain passing through a three-dimensional nanopore embedded in a nanochannel, where the presence of solvent molecules could effectively be treated as a random force acting on the coarse-grained polymer molecule [16,45,46]. In the present model, strong effects of intramolecular interactions on the inertial force were coarse-grained and the kinetics of ssDNA were mainly affected by external electric fields. In such a case, the behavior of a particle can be expressed by an over-damped Langevin equation [16,45,46]:
ζivi = −Ui + Fi + Ri
where ζi is the friction coefficient of the ith particle, −∇Ui is the conservative force, including interactions between particles, and Fi denotes the external electrostatic force, such that Fi = −Qiϕ, where Qi is the electric charge on the polymer molecule. For the purposes of a three-dimensional simulation, the electric potential, ϕ, in a rectangular nanofluidic channel was analyzed by solving for the Laplace equation 2ϕ = 0 with Neumann boundary conditions n·ϕ = 0 at the sidewall surfaces, where n was the surface normal vector, and with constant electric potentials at both ends of the channel. The FEM (Femtet®, Murata Software Co., Ltd., Tokyo, Japan) [53] was employed to solve for the electric potential. Fi was calculated by averaging the gradient of ϕ around each position [46]. In Equation (6), the random force Ri satisfied the fluctuation-dissipation theorem, such that:
Ijms 15 13817 i006
where kB was the Boltzmann constant, T was temperature, δij was Kronecker’s delta and δ(tt') was the Dirac delta function where t and t' were time. In this study, we focused on ssDNA and developed a bead-spring model for use in the Langevin dynamics simulations. Details of our coarse-grained model were also described in previous studies [16,45,46]. In order to model a ssDNA consisting of 48,000 nucleotides (48 knt), neighboring beads were connected with a harmonic spring [46]:
Ijms 15 13817 i007
where k was the spring constant and rij was the distance between the ith and jth particles. The equilibrium distance, req, between each connected pair of beads was defined as req = αrntNnt/N, where α was a variable parameter, Nnt was the number of nucleotides, and N was the number of beads. The equilibrium distance between the nucleotides in ssDNA is known to be rnt = 0.43 nm [59] and so, applying an α value of 0.847 [16], Nnt = 48,000, and N = 400, we obtained req = 43.7 nm. The above value for the parameter α was selected so as to properly replicate the radius of gyration (Rg) [60] of ssDNA, as well as the diffusion coefficient and electrophoretic mobility. The harmonic spring constant was calculated as k = kBT2, where T was set to 300 K and a χ value of 0.1σ was applied for thermal fluctuations based on previous studies [45], where σ was a Lennard–Jones parameter described below. Interactions between two beads, or between a bead and a channel wall, were represented by the Lennard–Jones potential, ULJ, taking into account the volume exclusion effect [46]:
Ijms 15 13817 i008
where σ was the characteristic length of ssDNA, and ω was the energy well-depth. A mirror reflection was assumed, meaning that the repulsive force from the wall effectively worked only along the direction perpendicular to the surface. The length parameter, σ, was determined from the persistence length of ssDNA necessary to reproduce the volume effect, such that σ = 5 nm [59]. ULJ was applied to non-adjacent beads and ω was set to kBT [16,46]. For the purposes of volume exclusion, the potential was truncated at to allow for purely repulsive interactions between the beads. The term ζi in Equations (6) and (7) was evaluated based on experimental measurements of the ssDNA diffusion coefficient, Di, according to NζiDi = kBT [52]. Applying N = 400, Di = 2.21 × 10−12 m2/s, and T = 300 K, ζi was determined to be 4.68 × 10−12 kg/s for each bead. Considering the existence of counterions around the ssDNA, the effective charge of an individual bead could be calculated according to Qi = ζiμi = μikBT/NDi [52]. Thus, based on the experimental value of μ = 2.84 × 10−8 m2/Vs [52], a Qi value resulted in −0.83e per bead (consisting 120 nt), where e is the elementary charge. This value was determined in terms of electrophoretic mobility of the coarse-grained ssDNA including counterions and thus, it might underestimate the monomer charge previously known [58,61]. In order to verify this quantity, we performed Langevin dynamics simulations for the ssDNA model in free solution.
The overall structure of the fluidic channel, including the reservoirs outside the nanochannel, was taken into account in the preliminary analysis, as shown in Figure 3a. There was a reservoir of 2.5 × 4.0 × 0.5 μm (length × width × height) on either side of the nanochannel and the electrodes were 2.5 μm from the nanochannel entrance. The electric potentials at the electrodes were set to −0.400 and 0.400 V at the cis and trans sides, respectively, based on the experimental conditions summarized in Table 1. Additionally, the Laplace equation was solved in the nanochannel and nanopore with a fine resolution of 10 nm.
At equilibrium, Rg was maintained in the vicinity of 300 nm, such that Rg2 was approximately equal to the product of the persistence length and the contour length [59,62]. Stable configurations such as this were employed as initial conditions for the simulations. The center of mass of the ssDNA was initially placed at a distance equivalent to Rg from the entrance of the nanochannel, as presented schematically in Figure 1a. In the next stage, the nonuniform electric field resulting from the FEM analysis was applied and the trajectories of the ssDNA were tracked. Equation (5) was integrated using the Euler algorithm with time steps of Δt = 10 ps [46].

4. Conclusions

In this study, we investigated the electrostatic potentials in nanopores embedded in a rectangular nanochannel. We obtained considerable agreement in the electric field strengths on the order of 106 V/m compared with previously published data [54]. Induction of strong electric fields in the narrowest space due to the connections of different-sized channels was confirmed [55]. Using such electric fields, we performed Langevin dynamics simulations by applying a coarse-grained model of ssDNA. The present model replicated the diffusion coefficient and electrophoretic mobility of long ssDNA, which allowed us to treat electrokinetic transport phenomena in the actual time and spatial scales. It was found that a nanoscale cross-sectional area was important with respect to uncoiling long-chained ssDNA molecules in a strong electric field and, as a result, reducing the translocation speed of the molecules. By adjusting the nanopore size, the quantity of ssDNA chains in the pore region can be constrained at a constant number, effectively producing a terminal velocity. With regard to the aim of obtaining single-molecule detection, this study suggests a preferred structure for nanofluidic channels.

Acknowledgments

The authors thank M. Tsutsui, M. Taniguchi, and T. Kawai, Osaka University, for fruitful comments on the electrokinetic transport of DNA in nanopores.

Author Contributions

W.Q. and K.M. mainly developed the simulation model and performed computations. W.Q., K.D., S.U., and S.K. analyzed the simulation results by developing the theoretical model. This article was written by W.Q., K.D., S.U., and S.K. All authors read and approved the final version of the manuscript.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Dekker, C. Solid-state nanopores. Nat. Nanotechnol. 2007, 2, 209–215. [Google Scholar] [CrossRef]
  2. Branton, D.; Deamer, D.W.; Marziali, A.; Bayley, H.; Benner, S.A.; Butler, T.; di Ventra, M.; Garaj, S.; Hibbs, A.; Huang, X.; et al. The potential and challenges of nanopore sequencing. Nat. Biotechnol. 2008, 26, 1146–1153. [Google Scholar] [CrossRef]
  3. Zwolak, M.; di Ventra, M. Colloquium: Physical approaches to DNA sequencing and detection. Rev. Mod. Phys. 2008, 80, 141–165. [Google Scholar]
  4. Venkatesan, B.M.; Bashir, R. Nanopore sensors for nucleic acid analysis. Nat. Nanotechnol. 2011, 6, 615–624. [Google Scholar]
  5. Kasianowicz, J.J.; Brandin, E.; Branton, D.; Deamer, D.W. Characterization of individual polynucleotide molecules using a membrane channel. Proc. Natl. Acad. Sci. USA 1996, 93, 13770–13773. [Google Scholar] [CrossRef]
  6. Clarke, J.; Wu, H.-C.; Jayasinghe, L.; Patel, A.; Reid, S.; Bayley, H. Continuous base identification for single-molecule nanopore DNA sequencing. Nat. Nanotechnol. 2009, 4, 265–270. [Google Scholar] [CrossRef]
  7. Lieberman, K.R.; Cherf, G.M.; Doody, M.J.; Olasagasti, F.; Kolodji, Y.; Akeson, M. Processive replication of single DNA molecules in a nanopore catalyzed by phi29 DNA polymerase. J. Am. Chem. Soc. 2010, 132, 17961–17972. [Google Scholar]
  8. Lagerqvist, J.; Zwolak, M.; di Ventra, M. Fast DNA sequencing via transverse electronic transport. Nano Lett. 2006, 6, 779–782. [Google Scholar]
  9. Liang, X.; Chou, S.Y. Nanogap detector inside nanofluidic channel for fast real-time label-free DNA analysis. Nano Lett. 2008, 8, 1472–1476. [Google Scholar] [CrossRef]
  10. Tsutsui, M.; Taniguchi, M.; Yokota, K.; Kawai, T. Identifying single nucleotides by tunnelling current. Nat. Nanotechnol. 2010, 5, 286–290. [Google Scholar] [CrossRef]
  11. Szarek, P.; Suwannawong, S.; Doi, K.; Kawano, S. Theoretical study on physicochemical aspects of a single molecular junction: Application to the bases of ssDNA. J. Phys. Chem. C 2013, 117, 109809:1–109809:9. [Google Scholar]
  12. Dorfman, K.D. DNA electrophoresis in microfabricated devices. Rev. Mod. Phys. 2010, 82, 2903–2947. [Google Scholar] [CrossRef]
  13. Sung, W.; Park, P.J. Polymer translocation through a pore in a membrane. Phys. Rev. Lett. 1996, 77, 783–786. [Google Scholar] [CrossRef]
  14. Muthukumar, M. Polymer translocation through a hole. J. Chem. Phys. 1999, 111, 10371–10374. [Google Scholar]
  15. Huopaniemi, I.; Luo, K.; Ala-Nissila, T.; Ying, S.-C. Langevin dynamics simulations of polymer translocation through nanopores. J. Chem. Phys. 2006, 125, 124901:1–124901:8. [Google Scholar]
  16. Nagahiro, S.; Kawano, S.; Kotera, H. Separation of long DNA chains using a nonuniform electric field: A numerical study. Phys. Rev. E 2007, 75, 011902:1–011902:5. [Google Scholar]
  17. Uehara, S.; Tsutsui, M.; Doi, K.; Taniguchi, M.; Kawano, S.; Kawai, T. Fluid dynamics and electrical detection of λDNA in electrode-embedded nanochannels. J. Biomech. Sci. Eng. 2013, 8, 244–256. [Google Scholar] [CrossRef]
  18. Meller, A.; Nivon, L.; Brandin, E.; Golovchenko, J.; Branton, D. Rapid nanopore discrimination between single polynucleotide molecules. Proc. Natl. Acad. Sci. USA 2000, 97, 1079–1084. [Google Scholar]
  19. Meller, A.; Nivon, L.; Branton, D. Voltage-driven DNA translocations through a nanopore. Phys. Rev. Lett. 2001, 86, 3435–3438. [Google Scholar] [CrossRef]
  20. Meller, A.; Branton, D. Single molecule measurements of DNA transport through a nanopore. Electrophoresis 2002, 23, 2583–2591. [Google Scholar] [CrossRef]
  21. Li, J.; Gershow, M.; Stein, D.; Brandin, E.; Golovchenko, J.A. DNA molecules and configurations in a solid-state nanopore microscope. Nat. Mater. 2003, 2, 611–615. [Google Scholar]
  22. Storm, A.J.; Storm, C.; Chen, J.; Zandbergen, H.; Joanny, J.-F.; Dekker, C. Fast DNA translocation through a solid-state nanopore. Nano Lett. 2005, 5, 1193–1197. [Google Scholar] [CrossRef]
  23. Krasilnikov, O.V.; Rodrigues, C.G.; Bezrukov, S.M. Single polymer molecules in a protein nanopore in the limit of a strong polymer-pore attraction. Phys. Rev. Lett. 2006, 97, 018301:1–018301:4. [Google Scholar]
  24. Skinner, G.M.; van den Hout, M.; Broekmans, O.; Dekker, C.; Dekker, N.H. Distinguishing single- and double-stranded nucleic acid molecules using solid-state nanopores. Nano Lett. 2009, 9, 2953–2960. [Google Scholar] [CrossRef]
  25. 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]
  26. Tsutsui, M.; He, Y.; Furuhashi, M.; Rahong, S.; Taniguchi, M.; Kawai, T. Transverse electric field dragging of DNA in a nanochannel. Sci. Rep. 2012, 2, 394. [Google Scholar]
  27. Fologea, D.; Uplinger, J.; Thomas, B.; McNabb, D.S.; Li, J. Slowing DNA translocation in a solid-state nanopore. Nano Lett. 2005, 5, 1734–1737. [Google Scholar] [CrossRef]
  28. 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]
  29. Lehtola, V.V.; Linna, R.P.; Kaski, K. Dynamics of forced biopolymer translocation. EPL 2009, 85, 58006:1–58006:6. [Google Scholar]
  30. Yong, H.; Wang, Y.; Yuan, S.; Xu, B.; Luo, K. Driven polymer translocation through a cylindrical nanochannel: Interplay between the channel length and the chain length. Soft Matter 2012, 8, 2769–2774. [Google Scholar] [CrossRef]
  31. Edmonds, C.M.; Hudiono, Y.C.; Ahmadi, A.G.; Hesketh, P.J.; Nair, S. Polymer translocation in solid-state nanopores: Dependence of scaling behavior on pore dimensions and applied voltage. AIP 2012, 136, 065105:1–065105:10. [Google Scholar]
  32. Luo, K.; Ala-Nissila, T.; Ying, S.-C.; Bhattacharya, A. Sequence dependence of DNA translocation through a nanopore. Phys. Rev. Lett. 2008, 100, 058101:1–058101:4. [Google Scholar]
  33. Ramachandran, A.; Guo, Q.; Iqbal, S.M.; Liu, Y. Coarse-grained molecular dynamics simulation of DNA translocation in chemically modified nanopores. J. Phys. Chem. B 2011, 115, 6138–6148. [Google Scholar]
  34. Forrey, C.; Muthukumar, M. Langevin dynamics simulations of dsDNA translocation through synthetic nanopores. J. Chem. Phys. 2007, 127, 015102:1–015102:10. [Google Scholar]
  35. Butler, T.Z.; Pavlenok, M.; Derrington, I.M.; Niederweis, M.; Gundlach, J.H. Single-molecule DNA detection with an engineered MspA protein nanopore. Proc. Natl. Acad. Sci. USA 2008, 105, 20647–20652. [Google Scholar]
  36. Wendell, D.; Jing, P.; Geng, J.; Subramaniam, V.; Lee, T.J.; Montemagnom, C.; Guo, P. Translocation of double stranded DNA through membrane adapted phi29 motor protein nanopore. Nat. Nanotechnol. 2009, 4, 765–772. [Google Scholar] [CrossRef]
  37. Franceschini, L.; Soskine, M.; Biesemans, A.; Maglia, G. A nanopore machine promotes the vectorial transport of DNA cross membranes. Nat. Commun. 2013, 4, 2415. [Google Scholar]
  38. Gauthier, M.G.; Slater, G.W. Molecular dynamics simulation of a polymer chain translocating through a nanoscopic pore. Eur. Phys. J. E 2008, 25, 17–23. [Google Scholar] [CrossRef]
  39. Izmitli, A.; Schwartz, D.C.; Graham, M.D.; de Pablo, J.J. The effect of hydrodynamic interactions on the dynamics of DNA translocation through pores. J. Chem. Phys. 2008, 128, 085102:1–085102:7. [Google Scholar]
  40. Fyta, M.; Melchionna, S.; Succi, S.; Kaxiras, E. Hydrodynamic correlations in the translocation of a biopolymer through a nanopore: Theory and multiscale simulations. Phys. Rev. E 2008, 78, 036704:1–036704:7. [Google Scholar]
  41. Luo, K.; Metzler, R. The chain sucker: Translocation dynamics of a polymer chain into a long narrow channel driven by longitudinal flow. J. Chem. Phys. 2011, 134, 135102:1–135102:8. [Google Scholar]
  42. Kawaguchi, C.; Noda, T.; Tsutsui, M.; Taniguchi, M.; Kawano, S.; Kawai, T. Electrical detection of single pollen allergen particles using electrode-embedded microchannels. J. Phys. 2012, 24, 164202:1–164202:6. [Google Scholar]
  43. Yasui., T.; Rahong, S.; Motoyama, K.; Yanagida, T.; Wu, Q.; Kaji, N.; Kanai, M.; Doi, K.; Nagashima, K.; Tokeshi, M.; et al. DNA manipulation and separation in sublithographic-scale nanowire array. ACS Nano 2013, 7, 3029–3035. [Google Scholar] [CrossRef]
  44. He, Y.; Tsutsui, M.; Fan, C.; Taniguchi, M.; Kawai, T. Controlling DNA translocation through gate modulation of nanopore wall surface charges. ACS Nano 2011, 5, 5509–5518. [Google Scholar] [CrossRef]
  45. Doi, K.; Haga, T.; Shintaku, H.; Kawano, S. Development of coarse-graining DNA models for single-nucleotide resolution analysis. Philos. Trans. R. Soc. A 2010, 368, 2615–2628. [Google Scholar] [CrossRef]
  46. Doi, K.; Qian, W.; Uehara, S.; Tsutsui, M.; Taniguchi, M.; Kawai, T.; Kawano, S. Langevin dynamics study on electrokinetic transport of long-chained DNA through nanogap embedded in nanochannel. Int. J. Emerg. Multidiscip. Fluid Sci. 2014, in press. [Google Scholar]
  47. Hanasaki, I.; Takahashi, H.; Sazaki, G.; Nakajima, K.; Kawano, S. Single-molecule measurements and dynamical simulations of protein molecules near silicon substrates. J. Phys. D 2008, 41, 095301:1–095301:9. [Google Scholar]
  48. Hanasaki, I.; Shintaku, H.; Matsunami, S.; Kawano, S. Structural and tensile properties of self-assembled DNA network on mica surface. Comput. Model. Eng. Sci. 2009, 46, 191–207. [Google Scholar]
  49. Doi, K.; Uemura, T.; Kawano, S. Molecular dynamics study of solvation effect on diffusivity changes of DNA fragments. J. Mol. Model. 2011, 17, 1457–1465. [Google Scholar]
  50. Doi, K.; Toyokita, Y.; Akamatsu, S.; Kawano, S. Reaction–diffusion wave model for self-assembled network formation of poly(dA)·poly(dT) DNA on mica and HOPG surfaces. Compt. Methods Biomech. Biomed. Eng. 2014, 17, 661–677. [Google Scholar] [CrossRef]
  51. Doi, K.; Takeuchi, H.; Nii, R.; Akamatsu, S.; Kakizaki, T.; Kawano, S. Self-assembly of 50 bp poly(dA)·poly(dT) DNA on highly oriented pyrolytic graphite via atomic force microscopy observation and molecular dynamics simulation. J. Chem. Phys. 2013, 139, 085102:1–085102:9. [Google Scholar]
  52. Stellwagen, E.; Lu, Y.; Stellwagen, N.C. Unified description of electrophoresis and diffusion for DNA and other polyions. Biochemistry 2003, 42, 11745–11750. [Google Scholar] [CrossRef]
  53. Femtet®; Murata Software Co., Ltd.: Tokyo, Japan, 2013.
  54. He, Y.; Tsutsui, M.; Fan, C.; Taniguchi, M.; Kawai, T. Gate manipulation of DNA capture into nanopores. ACS Nano 2011, 5, 8391–8397. [Google Scholar] [CrossRef]
  55. Trepargnier, E.H.; Radenovic, A.; Sivak, D.; Geissler, P.; Liphardt, J. Controlling DNA capture and propagation through artificial nanopores. Nano Lett. 2007, 7, 2824–2830. [Google Scholar] [CrossRef]
  56. Mannion, J.T.; Reccius, C.H.; Cross, J.D.; Craighead, H.G. Conformational analysis of single DNA molecules undergoing entropically induced motion in nanochannels. Biophys. J. 2006, 90, 4538–4545. [Google Scholar] [CrossRef]
  57. Salieb-Beugelaar, G.B.; Dorfman, K.D.; van den Berg, A.; Eijkel, J.C.T. Electrophoretic separation of DNA in gels and nanostructures. Lab. Chip 2009, 9, 2508–2523. [Google Scholar] [CrossRef]
  58. Uehara, S.; Shintaku, H.; Kawano, S. Electrokinetic flow dynamics of weakly aggregated λDNA confined in nanochannels. Trans. ASME J. Fluids Eng. 2011, 133, 121203:1–121203:8. [Google Scholar]
  59. Tinland, B.; Pluen, A.; Sturm, J.; Weill, G. Persistence length of single-stranded DNA. Macromolecules 1997, 30, 5763–5765. [Google Scholar]
  60. Kuszewski, J.; Gronenborn, A.M.; Clore, G.M. Improving the packing and accuracy of NMR structures with a pseudopotential for the radius of gyration. J. Am. Chem. Soc. 1999, 121, 2337–2338. [Google Scholar] [CrossRef]
  61. Manning, G.S. Limiting laws and counterion condensation in polyelectrolyte solutions. 7. Electrophoretic mobility and conductance. J. Phys. Chem. 1981, 85, 1506–1515. [Google Scholar]
  62. Rechendorff, K.; Witz, G.; Adamcik, J.; Dietler, G. Persistence length and scaling properties of single-stranded DNA adsorbed on modified graphite. J. Chem. Phys. 2009, 131, 095103:1–095103:6. [Google Scholar]

Share and Cite

MDPI and ACS Style

Qian, W.; Doi, K.; Uehara, S.; Morita, K.; Kawano, S. Theoretical Study of the Transpore Velocity Control of Single-Stranded DNA. Int. J. Mol. Sci. 2014, 15, 13817-13832. https://doi.org/10.3390/ijms150813817

AMA Style

Qian W, Doi K, Uehara S, Morita K, Kawano S. Theoretical Study of the Transpore Velocity Control of Single-Stranded DNA. International Journal of Molecular Sciences. 2014; 15(8):13817-13832. https://doi.org/10.3390/ijms150813817

Chicago/Turabian Style

Qian, Weixin, Kentaro Doi, Satoshi Uehara, Kaito Morita, and Satoyuki Kawano. 2014. "Theoretical Study of the Transpore Velocity Control of Single-Stranded DNA" International Journal of Molecular Sciences 15, no. 8: 13817-13832. https://doi.org/10.3390/ijms150813817

APA Style

Qian, W., Doi, K., Uehara, S., Morita, K., & Kawano, S. (2014). Theoretical Study of the Transpore Velocity Control of Single-Stranded DNA. International Journal of Molecular Sciences, 15(8), 13817-13832. https://doi.org/10.3390/ijms150813817

Article Metrics

Back to TopTop