Abstract
Owing to the potential use for real personalized genome sequencing, DNA sequencing with solid-state nanopores has been investigated intensively in recent time. However, the area still confronts problems and challenges. In this work, we present a brief overview of computational studies of key issues in DNA sequencing with solid-state nanopores by addressing the progress made in the last few years. We also highlight future challenges and prospects for DNA sequencing using this technology.
Introduction
Sequencing the human genome allows us to better understand the relationships among diseases, inheritance, and individuality. Although the increasing need for cheap and fast genome sequencing has promoted the development of new sequencing technologies, the cost of genomic sequencing is still far from the ideal price point of “the $1000 genome” (Porcu et al., ). This means that even cheaper and faster sequencing methods need to be developed. For this purpose, “single molecule” sequencing has been considered to be a promising technology (Harding and Keller, ; Pushkarev et al., ). In particular, “single molecule” sequencing with nanopores has been recommended as the next generation platform for DNA sequencing (Rhee and Burns, ; Maitra et al., ; Ku and Roukos, ).
The basic principle of DNA sequencing with a nanopore is illustrated in Figure 1. A DNA molecule, either double strand (dsDNA) or single-strand (ssDNA), dispersed in a salt solution (such as a KCl solution) is driven by an applied electric field to pass through a nanopore for sequencing. As the DNA translocates through the nanopore, the flow of ions is interrupted and the ionic current is blocked as a function of time. This phenomenon was first observed by Kasianowicz et al. in a study of a DNA molecule translocation through the α-haemolysin membrane pores (Kasianowicz et al., ). Later, it was found that the ionic current was associated with the particular nucleotides [adenine (A), cytosine (C), guanine (G), and thymine (T)] passing through the nanopore (Clarke et al., ). This makes it possible for nucleotides to be distinguished in terms of the detected current. DNA sequencing with a nanopore is thus to identify the composition of a DNA molecule from the pattern of the ionic current as the DNA passes through the nanopore in a salt solution (Purnell and Schmidt, ).
Figure 1
Both biological and solid-state nanopores can be used for DNA sequencing. It has been reported that sequence information can be obtained with biological nanopores (Cherf et al., ; Manrao et al., ). This reflects that DNA sequencing with biological nanopores has achieved a significant progress. However, the fact that the lipid membrane used to fix biological nanopores is delicately sensitive to the temperature, pH value and salt concentration makes it difficult to control their stability. For this reason, solid-state nanopores have become the promising alternative to biological nanopores (Merchant et al., ; Schneider et al., ; Garaj et al., ).
DNA sequencing with solid-state nanopores has experienced fast development in recent years (see for example, Venkatesan and Bashir, ,; Yang et al., for reviews). Yet, it still confronts problems, such as that the error rate is too large (Schadt et al., ) and the DNA translocation speed is too high. Theoretical studies can help us to understand the molecular details of DNA translocation through a nanopore. In this mini-review, we focus on the recent progress in computational studies on key issues in DNA sequencing with solid-state nanopores, such as control of the DNA translocation speed and the electronic signature of the identity of bases. We also address future challenges and prospects of DNA sequencing with this technology.
Control of the DNA translocation speed in solid-state nanopores
One of the main challenges in DNA sequencing with a solid-state nanopore is to reduce the speed of the DNA molecule passing through the nanopore (Kowalczyk and Dekker, ). The DNA translocation speed is still too high using the currently available solid-state nanopores. However, the development in the field indicates that ideal speeds are expected to be achieved around 2015 (Venkatesan and Bashir, ). DNA translocation through a nanopore is affected by many factors, such as the solvent viscosity, ion concentration, surface charge, and diameter of the pore. In theory, assuming that the end effects of a nanopore are negligible and that a DNA molecule can only move along the center axis of the pore, we can determine the translocation speed of the DNA in an applied electric field by solving a set of coupled differential equations (Ghosal, ; Kejian et al., ), i.e., the Stokes equation for the flow of the DNA to the pore surface and the Poission equation for the distribution of the electric potential in the pore. The resulting translocation speed is simply given as (Luan et al., ):
where ε is the dielectric constant of the electrolyte, ζd and ζw are the zeta potentials on the DNA and pore surfaces, respectively, η is the viscosity of the electrolyte, and E is the bias electric field. According to Equation (1), we can tune the electric field, solvent viscosity, etc. to control the speed. Experiments as well as theoretical calculations have been carried out to study the translocation time by changing the electric field or solvent viscosity (Japrung et al., ; Luan et al., ).
By lowering the temperature of a graphene nanopore, Fologea et al. decreased the translocation speed by a factor of ~2 (Fologea et al., ). Luan et al. demonstrated that the translocation speed of an ssDNA molecule could be reduced by 10 times in a SiO2 nanopore in a glycerol solvent (Luan et al., ,). By decreasing the size of the counter ions, for example from K+ to Na+ to Li+, Dekker et al. revealed that the translocation time of a DNA molecule through a solid-state nanopore increased significantly (Kowalczyk et al., ). Molecular dynamics (MD) simulations indicated that Li+ binds to DNA more strongly than K+ does. The stronger interaction of the ions with the DNA thus creates a much stronger drag because the movements of the ions and the DNA occur in the opposite directions (Kowalczyk et al., ).
In the above mentioned studies, the DNA translocation speed in the solid-state nanopores was indeed decreased. However, the variations in the translocation dynamics due to the DNA-pore interactions (Aksimentiev et al., ; Branton et al., ; Gierhart et al., ) were not reduced and accordingly the identification of different bases could not be enhanced. Computer modeling methods have been used to study the change of translocation dynamics of DNA passing through nanopores. Employing steered molecular dynamics (SMD) simulations, Luan et al. found that the translocation time was increased by chemical modification of the nanopore surface (Luan et al., ). The work of Luan et al. has enhanced our understanding of factors in decreasing the translocation speed in solid-state nanopores. By using a continuum-based model, Zhang et al. indicated that the local permittivity environment could reduce the DNA translocation speed (Zhang et al., ). Theoretical modeling has also been used to study the translocation time. Using the Lubensky-Nelson model (Lubensky and Nelson, ), Reimann et al. showed that the time distribution derived from the model matches the experimental data very well (Reimann et al., ). From Monte Carlo simulations, Polson et al. found that the translocation time is proportional to (N − NP)2, i.e., 〈τ〉 ∝ (N − NP)2 for a sufficiently narrow pore, where N means the polymer length and Np is the average number of monomers in the nanopore (Polson and Mccaffrey, ). In addition, a theoretical study reported by He et al. indicated that by utilizing the cross-pore thermal gradient, the DNA translocation speeds could be orders of magnitude slower than the electrophoretic counterpart (He et al., ).
Base identification with solid-state nanopores
There have been many theoretical studies on detecting nucleotide differences by solid-state nanopores (Postma, ; Wells et al., ). As indicated from computer simulations, A-T and G-C base pairs could be discriminated using graphene nanopores (Sathe et al., ). Through theoretical modeling, Liang et al. illustrated that the DNA translocation time could be extended and the discrimination of the bases could be improved by narrowing the nanopores (Liang et al., ). Aksimentiev et al. showed that the hydrophobic interactions of the nucleotides with the graphene membrane led to a dramatic reduction in the conformational fluctuations of the nucleotides in the pore and that the resulting ionic current blockades were different for different DNA nucleotides (Wells et al., ). The results from SMD simulations indicated that each nucleotide, except for cytosine and thymine, in an ssDNA could be identified and characterized by the sub-2-nanometer nanopore (Qiu and Guo, ). Garaj et al. pointed out that a molecule-hugging graphene nanopore could have a resolution higher than 0.6 nm along the length of the molecule (Garaj et al., ). In addition, theoretical calculations have also shown that different nucleotides can be discriminated by the difference in the conductance spectra (Nelson et al., ; Avdoshenko et al., ). In these studies, the read lengths of the DNA molecules were very short (shorter than 100 base pairs) and the data may therefore not have reflected the complexity in experiment where the number of base pairs is often more than 1000. However, theoretical studies significantly enhance our understanding of base identification for DNA sequencing with solid-state nanopores, especially with graphene nanopores, and can lay a solid theoretical foundation for DNA sequencing with nanopores in the future.
Future challenges and outlook
In using solid-state nanopore-based devices for DNA sequencing, the major concerns are the range of DNA lengths, the sequencing speed, and the error rate in the sequencing. To decrease the DNA translocation speed we can adopt some measures such as using more sticky solvents or decreasing the temperature. However, by these measures, the DNA translocation dynamics, which originates from the DNA-pore interactions, cannot be changed. The most promising approach to change the DNA translocation dynamics is to change the nanopore material because we can evidently not change the properties of the DNA. A few studied ways to deal with this issue include those changing the pore diameter (Liang et al., ), pore geometry (Aksimentiev et al., ; Venkatesan et al., ), the charge on the pore surface (Luan and Aksimentiev, ), and using surfaces modified by SAM (Luan et al., ).
Normally, solid-state nanopores are made of silicon nitride, aluminum oxide, or silicon oxide. In recent years, graphene nanopores have been suggested as a potential material for sequencing DNA because of their unique electric and mechanical properties. However, our understanding of DNA sequencing with graphene nanopores is still rather limited, necessitating further studies in this area.
Simulations can greatly enhance our understanding of the atomic details of DNA translocation dynamics and DNA-pore interactions. We have addressed the use of simulations to study the DNA translocation in solid-state nanopores. Many simulation studies have captured important physical properties specific to DNA sequencing with nanopores. However, the underlying simulation models are relatively simple, and may not reflect the complexity of the real systems. Compared to the experimental conditions used, the applied voltages in the simulations are often too large, leading to too high translocation speeds. This may produce unpredictable influence on the translocation dynamics. In a real system, the charge on the nanopore surface can also change with the pH value, solvent, and ionic concentration. Furthermore, the nanopore geometry can greatly modulate the drag forces, something that is not considered in most simulations. In addition, the error rate of DNA sequencing has not been assessed in most computational studies since the read lengths of the DNA molecules used in these studies have been very short. All these issues need to be attended in the future. Given so many new opportunities and given the considerable progress we have witnessed in this field, it can be foreseen that studies of DNA sequencing with solid-state nanopores will remain active for years to come.
Conflict of interest statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Statements
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
References
1
AksimentievA.HengJ. B.TimpG.SchultenK. (2004). Microscopic kinetics of DNA translocation through synthetic nanopores. Biophys. J. 87, 2086–2097. 10.1529/biophysj.104.042960
2
AvdoshenkoS. M.NozakiD.Gomes Da RochaC.GonzálezJ. W.LeeM. H.GutierrezR.et al. (2013). Dynamic and electronic transport properties of DNA translocation through graphene nanopores. Nano Lett. 13, 1969–1976. 10.1021/nl304735k
3
BrantonD.DeamerD. W.MarzialiA.BayleyH.BennerS. A.ButlerT.et al. (2008). The potential and challenges of nanopore sequencing. Nat. Biotechnol. 26, 1146–1153. 10.1038/nbt.1495
4
CherfG. M.LiebermanK. R.RashidH.LamC. E.KarplusK.AkesonM. (2012). Automated forward and reverse ratcheting of DNA in a nanopore at 5-A precision. Nat. Biotechnol. 30, 344–348. 10.1038/nbt.2147
5
ClarkeJ.WuH. C.JayasingheL.PatelA.ReidS.BayleyH. (2009). Continuous base identification for single-molecule nanopore DNA sequencing. Nat. Nanotechnol. 4, 265–270. 10.1038/nnano.2009.12
6
FologeaD.UplingerJ.ThomasB.McNabbD. S.LiJ. L. (2005). Slowing DNA translocation in a solid-state nanopore. Nano Lett. 5, 1734–1737. 10.1021/nl051063o
7
GarajS.LiuS.GolovchenkoJ. A.BrantonD. (2013). Molecule-hugging graphene nanopores. Proc. Natl. Acad. Sci. U.S.A. 110, 12192–12196. 10.1073/pnas.1220012110
8
GhosalS. (2007). Effect of salt concentration on the electrophoretic speed of a polyelectrolyte through a nanopore. Phys. Rev. Lett. 98:238104. 10.1103/PhysRevLett.98.238104
9
GierhartB. C.HowittD. G.ChenS. J.ZhuZ.KoteckiD. E.SmithR. L.et al. (2008). Nanopore with transverse nanoelectrodes for electrical characterization and sequencing of DNA. Sens. Actuators B Chem. 132, 593–600. 10.1016/j.snb.2007.11.054
10
HardingJ. D.KellerR. A. (1992). Single-molecule detection as an approach to rapid DNA sequencing. Trends Biotechnol. 10, 55–57. 10.1016/0167-7799(92)90170-Z
11
HeY.TsutsuiM.ScheicherR. H.BaiF.TaniguchiM.KawaiT. (2012). Thermophoretic Manipulation of DNA translocation through nanopores. ACS Nano7, 538–546. 10.1021/nn304914j
12
JaprungD.HenricusM.LiQ. H.MagliaG.BayleyH. (2010). Urea facilitates the translocation of single-stranded DNA and RNA through the alpha-hemolysin nanopore. Biophys. J. 98, 1856–1863. 10.1016/j.bpj.2009.12.4333
13
KasianowiczJ. J.BrandinE.BrantonD.DeamerD. W. (1996). Characterization of individual polynucleotide molecules using a membrane channel. Proc. Natl. Acad. Sci. U.S.A. 93, 13770–13773. 10.1073/pnas.93.24.13770
14
KejianD.SunW. M.ZhangH. Y.PengX. L.HuH. G. (2009). Dependence of zeta potential on polyelectrolyte moving through a solid-state nanopore. Appl. Phys. Lett. 94, 014101–014104-3. 10.1063/1.3065029
15
KowalczykS. W.DekkerC. (2012). Measurement of the docking time of a DNA molecule onto a solid-state nanopore. Nano Lett. 12, 4159–4163. 10.1021/nl301719a
16
KowalczykS. W.WellsD. B.AksimentievA.DekkerC. (2012). Slowing down DNA translocation through a nanopore in lithium chloride. Nano Lett. 12, 1038–1044. 10.1021/nl204273h
17
KuC. S.RoukosD. H. (2013). From next-generation sequencing to nanopore sequencing technology: paving the way to personalized genomic medicine. Expert Rev. Med. Devices10, 1–6. 10.1586/erd.12.63
18
LiangL.CuiP.WangQ.WuT.AgrenH.TuY. (2013). Theoretical study on key factors in DNA sequencing with graphene nanopores. RSC Adv. 3, 2445–2453. 10.1039/c2ra22109h
19
LuanB.AfzaliA.HarrerS.PengH.WaggonerP.PolonskyS.et al. (2010). Tribological effects on DNA translocation in a nanochannel coated with a self-assembled monolayer. J. Phys. Chem. B114, 17172–17176. 10.1021/jp108865q
20
LuanB.AksimentievA. (2008). Electro-osmotic screening of the DNA charge in a nanopore. Phys. Rev. E Stat. Nonlin. Soft Matter Phys. 78:021912. 10.1103/PhysRevE.78.021912
21
LuanB.WangD.ZhouR.HarrerS.PengH.StolovitzkyG. (2012b). Dynamics of DNA translocation in a solid-state nanopore immersed in aqueous glycerol. Nanotechnology23:455102. 10.1088/0957-4484/23/45/455102
22
LuanB. Q.StolovitzkyG.MartynaG. (2012a). Slowing and controlling the translocation of DNA in a solid-state nanopore. Nanoscale4, 1068–1077. 10.1039/c1nr11201e
23
LubenskyD. K.NelsonD. R. (1999). Driven polymer translocation through a narrow pore. Biophys. J. 77, 1824–1838. 10.1016/S0006-3495(99)77027-X
24
MaitraR. D.KimJ.DunbarW. B. (2012). Recent advances in nanopore sequencing. Electrophoresis33, 3418–3428. 10.1002/elps.201200272
25
ManraoE. A.DerringtonI. M.LaszloA. H.LangfordK. W.HopperM. K.GillgrenN.et al. (2012). Reading DNA at single-nucleotide resolution with a mutant MspA nanopore and phi29 DNA polymerase. Nat. Biotechnol. 30, 349–353. 10.1038/nbt.2171
26
MerchantC. A.HealyK.WanunuM.RayV.PetermanN.BartelJ.et al. (2010). DNA Translocation through Graphene Nanopores. Nano Lett. 10, 2915–2921. 10.1021/nl101046t
27
NelsonT.ZhangB.PrezhdoO. V. (2010). Detection of nucleic acids with graphene nanopores: ab initio characterization of a novel sequencing device. Nano Lett. 10, 3237–3242. 10.1021/nl9035934
28
PolsonJ. M.MccaffreyA. C. M. (2013). Polymer translocation dynamics in the quasi-static limit. J. Chem. Phys. 138, 174902–174912. 10.1063/1.4803022
29
PorcuE.SannaS.FuchsbergerC.FritscheL. G. (2013). Genotype imputation in genome-wide association studies. Curr. Protoc. Hum. Genet. 78, 1.25.1–1.25.14. 10.1002/0471142905.hg0125s78
30
PostmaH. W. (2010). Rapid sequencing of individual DNA molecules in graphene nanogaps. Nano Lett. 10, 420–425. 10.1021/nl9029237
31
PurnellR. F.SchmidtJ. J. (2009). Discrimination of single base substitutions in a DNA strand immobilized in a biological nanopore. ACS Nano3, 2533–2538. 10.1021/nn900441x
32
PushkarevD.NeffN. F.QuakeS. R. (2009). Single-molecule sequencing of an individual human genome. Nat. Biotechnol. 27, 847–850. 10.1038/nbt.1561
33
QiuH.GuoW. (2012). Detecting ssDNA at single-nucleotide resolution by sub-2-nanometer pore in monoatomic graphene: a molecular dynamics study. Appl. Phys. Lett. 100, 083106–083104. 10.1063/1.3686921
34
ReimannP.MeyerA.GetfertS. (2012). On the lubensky-nelson model of polymer translocation through nanopores. Biophys. J. 103, 889–897. 10.1016/j.bpj.2012.07.036
35
RheeM.BurnsM. A. (2006). Nanopore sequencing technology: research trends and applications. Trends Biotechnol. 24, 580–586. 10.1016/j.tibtech.2006.10.005
36
SatheC.ZouX. Q.LeburtonJ. P.SchultenK. (2011). Computational investigation of DNA detection using graphene nanopores. ACS Nano5, 8842–8851. 10.1021/nn202989w
37
SchadtE. E.TurnerS.KasarskisA. (2010). A window into third-generation sequencing. Hum. Mol. Genet. 19, R227–R240. 10.1093/hmg/ddq416
38
SchneiderG. F.KowalczykS. W.CaladoV. E.PandraudG.ZandbergenH. W.VandersypenL. M. K.et al. (2010). DNA translocation through graphene nanopores. Nano Lett. 10, 3163–3167. 10.1021/nl102069z
39
VenkatesanB. M.BashirR. (2011a). Nanopore sensors for nucleic acid analysis. Nat. Nanotechnol. 6, 615–624. 10.1038/nnano.2011.129
40
VenkatesanB. M.BashirR. (2011b). Solid-state nanopore sensors for nucleic acid analysis. Nanopores Sens. Fund. Biol. Interact. 1–33. 10.1007/978-1-4419-8252-0_1
41
VenkatesanB. M.DorvelB.YemeniciogluS.WatkinsN.PetrovI.BashirR. (2009). Highly Sensitive, Mechanically Stable Nanopore Sensors for DNA Analysis. Adv. Mater. 21, 2771. 10.1002/adma.200803786
42
WellsD. B.BelkinM.ComerJ.AksimentievA. (2012). Assessing graphene nanopores for sequencing DNA. Nano Lett. 12, 4117–4123. 10.1021/nl301655d
43
YangY.LiuR.XieH.HuiY.JiaoR.GongY.et al. (2013). Advances in nanopore sequencing technology. J. Nanosci. Nanotechnol. 13, 4521–4538. 10.1166/jnn.2013.7756
44
ZhangM.YehL.-H.QianS.HsuJ.-P.JooS. W. (2012). DNA Electrokinetic translocation through a nanopore: local permittivity environment effect. J. Phys. Chem. C116, 4793–4801. 10.1021/jp211798x
Summary
Keywords
solid-state nanopore, graphene nanopore, computer modeling, bio-nanotechnology, DNA sequencing, theory, molecular simulation
Citation
Liang L, Wang Q, Ågren H and Tu Y (2014) Computational studies of DNA sequencing with solid-state nanopores: key issues and future prospects. Front. Chem. 2:5. doi: 10.3389/fchem.2014.00005
Received
26 August 2013
Accepted
05 February 2014
Published
21 February 2014
Volume
2 - 2014
Edited by
Alan Aspuru-Guzik, Harvard University, USA
Reviewed by
Joshua Schrier, Haverford College, USA; Ryan Rollings, Harvard University, USA
Copyright
© 2014 Liang, Wang, Ågren and Tu.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yaoquan Tu, Division of Theoretical Chemistry and Biology, School of Biotechnology, KTH Royal Institute of Technology, Roslagstullsbacken 15, SE-10691 Stockholm, Sweden e-mail: tu@theochem.kth.se
This article was submitted to The oretical and Computational Chemistry, a section of the journal Frontiers in Chemistry.
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.