Microscopic Kinetics of DNA Translocation through synthetic nanopores

We have previously demonstrated that a nanometer-diameter pore in a nanometer-thick metal-oxide-semiconductor-compatible membrane can be used as a molecular sensor for detecting DNA. The prospects for using this type of device for sequencing DNA are avidly being pursued. The key attribute of the sensor is the electric field-induced (voltage-driven) translocation of the DNA molecule in an electrolytic solution across the membrane through the nanopore. To complement ongoing experimental studies developing such pores and measuring signals in response to the presence of DNA, we conducted molecular dynamics simulations of DNA translocation through the nanopore. A typical simulated system included a patch of a silicon nitride membrane dividing water solution of potassium chloride into two compartments connected by the nanopore. External electrical fields induced capturing of the DNA molecules by the pore from the solution and subsequent translocation. Molecular dynamics simulations suggest that 20-basepair segments of double-stranded DNA can transit a nanopore of 2.2 x 2.6 nm(2) cross section in a few microseconds at typical electrical fields. Hydrophobic interactions between DNA bases and the pore surface can slow down translocation of single-stranded DNA and might favor unzipping of double-stranded DNA inside the pore. DNA occluding the pore mouth blocks the electrolytic current through the pore; these current blockades were found to have the same magnitude as the blockade observed when DNA transits the pore. The feasibility of using molecular dynamics simulations to relate the level of the blocked ionic current to the sequence of DNA was investigated.     

An ON/OFF biosensor based on blockade of ionic current passing through a solid-state nanopore

Single nanopores have attracted interest for their use as biosensing devices. In general, methods involve measuring ionic current blockades associated with translocation of analytes through the nanopore, but the detection of such short time lasting events requires complex equipment and setup that are critical for convenient routine biosensing. Here we present a novel biosensing concept based on a single nanopore in a silicon nitride membrane and two anchor-linked DNA species that forms trans-pore hybrids, realizing a stable blockade of ionic current through the pore. Molecular recognition events affecting the DNA hybrids cause a pore opening and the consequent establishment of an ionic current. In the present implementation of the device, we constructed a magnetic bead/streptavidin/biotin-DNA1/DNA2-biotin/streptavidin/Quantumdot-cluster complex (where DNA1 is a mismatched reverse complement of DNA2) through a sub-micrometric pore and monitored DNA strand displacement events occurring after addition of an oligonucleotide complementary to DNA2. The electric and mechanical aspects of the novel device, as well as its potential in biosensing are discussed.     

Microscopic Mechanics of Hairpin DNA Translocation through Synthetic Nanopores

Nanoscale pores have proved useful as a means to assay DNA and are actively being developed as the basis of genome sequencing methods. Hairpin DNA (hpDNA), having both double-helical and overhanging coil portions, can be trapped in a nanopore, giving ample time to execute a sequence measurement. In this article, we provide a detailed account of hpDNA interaction with a synthetic nanopore obtained through extensive all-atom molecular dynamics simulations. For synthetic pores with minimum diameters from 1.3 to 2.2 nm, we find that hpDNA can translocate by three modes: unzipping of the double helix and—in two distinct orientations—stretching/distortion of the double helix. Furthermore, each of these modes can be selected by an appropriate choice of the pore size and voltage applied transverse to the membrane. We demonstrate that the presence of hpDNA can dramatically alter the distribution of ions within the pore, substantially affecting the ionic current through it. In experiments and simulations, the ionic current relative to that in the absence of DNA can drop below 10% and rise beyond 200%. Simulations associate the former with the double helix occupying the constriction and the latter with accumulation of DNA that has passed through the constriction.     

Translocation of structured polynucleotides through nanopores

We investigate theoretically the translocation of structured RNA/DNA molecules through narrow pores which allow single but not double strands to pass. The unzipping of basepaired regions within the molecules presents significant kinetic barriers for the translocation process. We show that this circumstance may be exploited to determine the full basepairing pattern of polynucleotides, including RNA pseudoknots. The crucial requirement is that the translocation dynamics (i.e. the length of the translocated molecular segment) needs to be recorded as a function of time with a spatial resolution of a few nucleotides. This could be achieved, for instance, by applying a mechanical driving force for translocation and recording force-extension curves (FECs) with a device such as an atomic force microscope or optical tweezers. Our analysis suggests that, with this added spatial resolution, nanopores could be transformed into a powerful experimental tool to study the folding of nucleic acids.     

Long dwell-time passage of DNA through nanometer-scale pores: kinetics and sequence dependence of motion

A detailed understanding of the kinetics of DNA motion though nanometer-scale pores is important for the successful development of many of the proposed next-generation rapid DNA sequencing and analysis methods. Many of these approaches require DNA motion through nanopores to be slowed by several orders of magnitude from its native translocation velocity so that the translocation times for individual nucleotides fall within practical timescales for detection. With the increased dwell time of DNA in the pore, DNA-pore interactions begin to play an increasingly important role in translocation kinetics. In previous work, we and others observed that when the DNA dwell time in the pore is substantial (>1 ms), DNA motion in α-hemolysin (α-HL) pores leads to nonexponential kinetics in the escape of DNA out of the pore. Here we show that a three-state model for DNA escape, involving stochastic binding interactions of DNA with the pore, accurately reproduces the experimental data. In addition, we investigate the sequence dependence of the DNA escape process and show that the interaction strength of adenine with α-HL is substantially lower relative to cytosine. Our results indicate a difference in the process by which DNA moves through an α-HL nanopore when the motion is fast (microsecond timescale) as compared with when it is slow (millisecond timescale) and strongly influenced by DNA-pore interactions of the kind reported here. We also show the ability of wild-type α-HL to detect and distinguish between 5-methylcytosine and cytosine based on differences in the absolute ionic current through the pore in the presence of these two nucleotides. The results we present here regarding sequence-dependent (and dwell-time-dependent) DNA-pore interaction kinetics will have important implications for the design of methods for DNA analysis through reduced-velocity motion in nanopores.     

Reading DNA at single-nucleotide resolution with a mutant MspA nanopore and phi29 DNA polymerase

Nanopore technologies are being developed for fast and direct sequencing of single DNA molecules through detection of ionic current modulations as DNA passes through a pore's constriction. Here we demonstrate the ability to resolve changes in current that correspond to a known DNA sequence by combining the high sensitivity of a mutated form of the protein pore Mycobacterium smegmatis porin A (MspA) with phi29 DNA polymerase (DNAP), which controls the rate of DNA translocation through the pore. As phi29 DNAP synthesizes DNA and functions like a motor to pull a single-stranded template through MspA, we observe well-resolved and reproducible ionic current levels with median durations of ～28 ms and ionic current differences of up to 40 pA. Using six different DNA sequences with readable regions 42-53 nucleotides long, we record current traces that map to the known DNA sequences. With single-nucleotide resolution and DNA translocation control, this system integrates solutions to two long-standing hurdles to nanopore sequencing.     

DNA translocation governed by interactions with solid-state nanopores

We investigate the voltage-driven translocation dynamics of individual DNA molecules through solid-state nanopores in the diameter range 2.7-5 nm. Our studies reveal an order of magnitude increase in the translocation times when the pore diameter is decreased from 5 to 2.7 nm, and steep temperature dependence, nearly threefold larger than would be expected if the dynamics were governed by viscous drag. As previously predicted for an interaction-dominated translocation process, we observe exponential voltage dependence on translocation times. Mean translocation times scale with DNA length by two power laws: for short DNA molecules, in the range 150-3500 bp, we find an exponent of 1.40, whereas for longer molecules, an exponent of 2.28 dominates. Surprisingly, we find a transition in the fraction of ion current blocked by DNA, from a length-independent regime for short DNA molecules to a regime where the longer the DNA, the more current is blocked. Temperature dependence studies reveal that for increasing DNA lengths, additional interactions are responsible for the slower DNA dynamics. Our results can be rationalized by considering DNA/pore interactions as the predominant factor determining DNA translocation dynamics in small pores. These interactions markedly slow down the translocation rate, enabling higher temporal resolution than observed with larger pores. These findings shed light on the transport properties of DNA in small pores, relevant for future nanopore applications, such as DNA sequencing and genotyping.     

The electromechanics of DNA in a synthetic nanopore

We have explored the electromechanical properties of DNA on a nanometer-length scale using an electric field to force single molecules through synthetic nanopores in ultrathin silicon nitride membranes. At low electric fields, E < 200 mV/10 nm, we observed that single-stranded DNA can permeate pores with a diameter >/=1.0 nm, whereas double-stranded DNA only permeates pores with a diameter >/=3 nm. For pores <3.0 nm diameter, we find a threshold for permeation of double-stranded DNA that depends on the electric field and pH. For a 2 nm diameter pore, the electric field threshold is approximately 3.1 V/10 nm at pH = 8.5; the threshold decreases as pH becomes more acidic or the diameter increases. Molecular dynamics indicates that the field threshold originates from a stretching transition in DNA that occurs under the force gradient in a nanopore. Lowering pH destabilizes the double helix, facilitating DNA translocation at lower fields.     

Theoretical study of sequence-dependent nanopore unzipping of DNA

We theoretically investigate the unzipping of DNA electrically driven through a nanometer-size pore. Taking the DNA base sequence explicitly into account, the unpairing and translocation process is described by a biased random walk in a one-dimensional energy landscape determined by the sequential basepair opening. Distributions of translocation times are numerically calculated as a function of applied voltage and temperature. We show that varying these two parameters changes the dynamics from a predominantly diffusive behavior to a dynamics governed by jumps over local energy barriers. The work suggests experimentally studying sequence effects, by comparing the average value and standard deviation of the statistical distribution of translocation times.     

DnaB helicase activity is modulated by DNA geometry and force

The replicative helicase for Escherichia coli is DnaB, a hexameric, ring-shaped motor protein that encircles and translocates along ssDNA, unwinding dsDNA in advance of its motion. The microscopic mechanisms of DnaB are unknown; further, prior work has found that DnaB's activity is modified by other replication proteins, indicating some mechanistic flexibility. To investigate these issues, we quantified translocation and unwinding by single DnaB molecules in three tethered DNA geometries held under tension. Our data support the following conclusions: 1), Unwinding by DnaB is enhanced by force-induced destabilization of dsDNA. 2), The magnitude of this stimulation varies with the geometry of the tension applied to the DNA substrate, possibly due to interactions between the helicase and the occluded ssDNA strand. 3), DnaB unwinding and (to a lesser extent) translocation are interrupted by pauses, which are also dependent on force and DNA geometry. 4), DnaB moves slower when a large tension is applied to the helicase-bound strand, indicating that it must perform mechanical work to compact the strand against the applied force. Our results have implications for the molecular mechanisms of translocation and unwinding by DnaB and for the means of modulating DnaB activity.     

Dynamics of the DNA duplex formation studied by single molecule force measurements

DNA is partly denatured in vitro by applying a force that mechanically separates the two strands of the double helix. Sudden reduction of the imposed displacement triggers spontaneous reannealing of the molecule. The corresponding force signals are measured by optical trapping interferometry for backward steps of various amplitudes and base sequence intervals. The measured signals frequently show plateaus of varying duration at discrete values that are dependent on the base sequence. Additional measurements are performed with proteins bound to the double helix. When the opening fork encounters such a protein during mechanical unzipping, force increases until the protein is ejected. This ejection induces fast release of tension and fast unzipping. Comparing our different measurements, we find that both DNA unzipping and the relaxation of tension in DNA are faster than the formation of the double helix.     

Pore Direction in Relation to Anisotropy of Mechanical Strength in a Cubic Starch Compact

The purpose of this research was to evaluate the relation between preferential direction of pores and mechanical strength of cubic starch compacts. The preferential pore direction was quantified in SEM images of cross sections of starch compacts using a previously described algorithm for determination of the quotient of transitions (Q). This parameter and the mechanical strength were evaluated in compacts of different porosities. Starch was chosen as a model compound for materials with ductile behaviour of which tablets with low porosities can be made and which shows some elastic recovery after compaction. At medium and high porosity Q was significantly higher in the images providing a side view of the compact than in the images providing a top view (0.973 vs. 0.927 and 0.958 vs. 0.874 at 0 mm from the side of the compact and 0.956 vs. 0.854 and 0.951 vs. 0.862 at 3.5 mm), indicating that the pores were mainly oriented in the direction perpendicular to the direction of compression. This was accompanied by a lower crushing force in this direction. This could be explained by considering the pores as cracks which propagate through the sample during crushing. For both directions the crushing force decreased with increasing porosity. The yield strength of the compacts also decreased with increasing porosity, but this parameter was not dependent on the direction of crushing when the porosity was below 10%. The results show that pore direction significantly influences the crushing force but does not influence the yield strength, at porosities below 10%.     

Probing DNA base pairing energy profiles using a nanopore

We experimentally show that the voltage driven unzipping of long DNA duplexes by an α-hemolysin pore is sensitive to the shape of the base pairing energy landscape. Two sequences of equal global stability were investigated. The sequence with an homogeneous base pairing profile translocates faster than the one with alternative weak and strong regions. We could qualitatively account for theses observations by theoretically describing the voltage driven translocation as a biased random walk of the unzipping fork in the sequence dependent energy landscape. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Smooth DNA Transport through a Narrowed Pore Geometry

Voltage-driven transport of double-stranded DNA through nanoscale pores holds much potential for applications in quantitative molecular biology and biotechnology, yet the microscopic details of translocation have proven to be challenging to decipher. Earlier experiments showed strong dependence of transport kinetics on pore size: fast regular transport in large pores (> 5 nm diameter), and slower yet heterogeneous transport time distributions in sub-5 nm pores, which imply a large positional uncertainty of the DNA in the pore as a function of the translocation time. In this work, we show that this anomalous transport is a result of DNA self-interaction, a phenomenon that is strictly pore-diameter dependent. We identify a regime in which DNA transport is regular, producing narrow and well-behaved dwell-time distributions that fit a simple drift-diffusion theory. Furthermore, a systematic study of the dependence of dwell time on DNA length reveals a single power-law scaling of 1.37 in the range of 35–20,000 bp. We highlight the resolution of our nanopore device by discriminating via single pulses 100 and 500 bp fragments in a mixture with >98% accuracy. When coupled to an appropriate sequence labeling method, our observation of smooth DNA translocation can pave the way for high-resolution DNA mapping and sizing applications in genomics.     

Quantitative analysis of the nanopore translocation dynamics of simple structured polynucleotides

Nanopore translocation experiments are increasingly applied to probe the secondary structures of RNA and DNA molecules. Here, we report two vital steps toward establishing nanopore translocation as a tool for the systematic and quantitative analysis of polynucleotide folding: 1), Using α-hemolysin pores and a diverse set of different DNA hairpins, we demonstrate that backward nanopore force spectroscopy is particularly well suited for quantitative analysis. In contrast to forward translocation from the vestibule side of the pore, backward translocation times do not appear to be significantly affected by pore-DNA interactions. 2), We develop and verify experimentally a versatile mesoscopic theoretical framework for the quantitative analysis of translocation experiments with structured polynucleotides. The underlying model is based on sequence-dependent free energy landscapes constructed using the known thermodynamic parameters for polynucleotide basepairing. This approach limits the adjustable parameters to a small set of sequence-independent parameters. After parameter calibration, the theoretical model predicts the translocation dynamics of new sequences. These predictions can be leveraged to generate a baseline expectation even for more complicated structures where the assumptions underlying the one-dimensional free energy landscape may no longer be satisfied. Taken together, backward translocation through α-hemolysin pores combined with mesoscopic theoretical modeling is a promising approach for label-free single-molecule analysis of DNA and RNA folding.     

Nanopore unzipping of individual DNA hairpin molecules

We have used the nanometer scale alpha-Hemolysin pore to study the unzipping kinetics of individual DNA hairpins under constant force or constant loading rate. Using a dynamic voltage control method, the entry rate of polynucleotides into the pore and the voltage pattern applied to induce hairpin unzipping are independently set. Thus, hundreds of unzipping events can be tested in a short period of time (few minutes), independently of the unzipping voltage amplitude. Because our method does not entail the physical coupling of the molecules under test to a force transducer, very high throughput can be achieved. We used our method to study DNA unzipping kinetics at small forces, which have not been accessed before. We find that in this regime the static unzipping times decrease exponentially with voltage with a characteristic slope that is independent of the duplex region sequence, and that the intercept depends strongly on the duplex region energy. We also present the first nanopore dynamic force measurements (time varying force). Our results are in agreement with the approximately logV dependence at high V (where V is the loading rate) observed by other methods. The extension of these measurements to lower loading rates reveals a much weaker dependence on V.     

Determination of RNA orientation during translocation through a biological nanopore

We investigate single-molecule electrophoretic translocation of A(50), C(50), A(25)C(50), and C(50)A(25) RNA molecules through the alpha-hemolysin transmembrane protein pore. We observe pronounced bilevel current blockages during translocation of A(25)C(50) and C(50)A(25) molecules. The two current levels observed during these bilevel blockages are very similar to the characteristic current levels observed during A(50) and C(50) translocation. From the temporal ordering of the two levels within the bilevel current blockages, we infer whether individual A(25)C(50) and C(50)A(25) molecules pass through the pore in a 3'-->5' or 5'-->3' orientation. Correlation between the level of current obstruction and the inferred A(25)C(50) or C(50)A(25) orientation indicates that 3'-->5' translocation of a poly C segment causes a significantly deeper current obstruction than 5'-->3' translocation. Our analysis also suggests that the 3' ends of C(50) and A(25)C(50) RNA molecules are more likely to initiate translocation than the 5' ends. Orientation dependent differences in a smaller current blockage that immediately precedes many translocation events suggest that this blockage also contains information about RNA orientation during translocation. These findings emphasize that the directionality of polynucleotide molecules is an important factor in translocation and demonstrate how structure within ionic current signals can give new insights into the translocation process.     

Single-Stranded DNA within Nanopores: Conformational Dynamics and Implications for Sequencing; a Molecular Dynamics Simulation Study

Engineered protein nanopores, such as those based on α-hemolysin from Staphylococcus aureus have shown great promise as components of next-generation DNA sequencing devices. However, before such protein nanopores can be used to their full potential, the conformational dynamics and translocation pathway of the DNA within them must be characterized at the individual molecule level. Here, we employ atomistic molecular dynamics simulations of single-stranded DNA movement through a model α-hemolysin pore under an applied electric field. The simulations enable characterization of the conformations adopted by single-stranded DNA, and allow exploration of how the conformations may impact on translocation within the wild-type model pore and a number of mutants. Our results show that specific interactions between the protein nanopore and the DNA can have a significant impact on the DNA conformation often leading to localized coiling, which in turn, can alter the order in which the DNA bases exit the nanopore. Thus, our simulations show that strategies to control the conformation of DNA within a protein nanopore would be a distinct advantage for the purposes of DNA sequencing.     

Processive translocation mechanism of the human Bloom’s syndrome helicase along single-stranded DNA

BLM, one of the human RecQ helicases, plays a fundamental role in homologous recombination-based error-free DNA repair pathways, which require its translocation and DNA unwinding activities. Although translocation is essential in vivo during DNA repair processes and it provides a framework for more complex activities of helicases, including strand separation and nucleoprotein displacement, its mechanism has not been resolved for any human DNA helicase. Here, we present a quantitative model for the translocation of a monomeric form of BLM along ssDNA. We show that BLM performs translocation at a low adenosine triphosphate (ATP) coupling ratio (1 ATP consumed/1 nucleotide traveled) and moderate processivity (with a mean number of 50 nucleotides traveled in a single run). We also show that the rate-limiting step of the translocation cycle is a transition between two ADP-bound enzyme states. Via opening of the helicase core, this structural change may drive the stepping of BLM along the DNA track by a directed inchworm mechanism. The data also support the conclusion that BLM performs double-stranded DNA unwinding by fully active duplex destabilization.     

Osmium-Based Pyrimidine Contrast Tags for Enhanced Nanopore-Based DNA Base Discrimination

Nanopores are a promising platform in next generation DNA sequencing. In this platform, an individual DNA strand is threaded into nanopore using an electric field, and enzyme-based ratcheting is used to move the strand through the detector. During this process the residual ion current through the pore is measured, which exhibits unique levels for different base combinations inside the pore. While this approach has shown great promise, accuracy is not optimal because the four bases are chemically comparable to one another, leading to small differences in current obstruction. Nucleobase-specific chemical tagging can be a viable approach to enhancing the contrast between different bases in the sequence. Herein we show that covalent modification of one or both of the pyrimidine bases by an osmium bipyridine complex leads to measureable differences in the blockade amplitudes of DNA molecules. We qualitatively determine the degree of osmylation of a DNA strand by passing it through a solid-state nanopore, and are thus able to gauge T and C base content. In addition, we show that osmium bipyridine reacts with dsDNA, leading to substantially different current blockade levels than exhibited for bare dsDNA. This work serves as a proof of principle for nanopore sequencing and mapping via base-specific DNA osmylation.