The impact of the number of layers of the graphene nanopores and the electrical field on ssDNA translocation

Graphene-based nanopore devices hold great promise for the next generation DNA sequencing because graphene is atomically thin which is extremely important for single base recognition. To understand the fundamental details of DNA translocation through a graphene nanopore, in this work, molecular dynamics simulations of ssDNA translocation through the nanopore were performed to trace the nucleobase trajectories and to investigate the impact of the number of layers of the graphene membrane and the electrical field on ssDNA translocation. We found that the velocity of ssDNA translocation was speeded up with the higher bias voltage, and the two-layered and five-layered graphene membrane with 1.0-nm diameter circular nanopore could discern different DNA strand by the translocation time.     

A Brownian ratchet for protein translocation including dissociation of ratcheting sites

We study a model for the translocation of proteins across membranes through a nanopore using a ratcheting mechanism. When the protein enters the nanopore it diffuses in and out of the pore according to a Brownian motion. Moreover, it is bound by ratcheting molecules which hinder the diffusion of the protein out of the nanopore, i.e. the Brownian motion is reflected such that no ratcheting molecule exits the pore. New ratcheting molecules bind at rate γ. Extending our previous approach (Depperschmidt and Pfaffelhuber in Stoch Processes Appl 120:901–925, 2010) we allow the ratcheting molecules to dissociate (at rate δ) from the protein (Model I). We also provide an approximate model (Model II) which assumes a Poisson equilibrium of ratcheting molecules on one side of the current reflection boundary. Using analytical methods and simulations we show that the speeds of both models are approximately the same. Our analytical results on Model II give the speed of translocation by means of a solution of an ordinary differential equation. This speed gives an approximation for the time it takes to translocate a protein of given length.     

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.     

Mechanism of How Salt-Gradient-Induced Charges Affect the Translocation of DNA Molecules through a Nanopore

Experiments using nanopores demonstrated that a salt gradient enhances the capture rate of DNA and reduces its translocation speed. These two effects can help to enable electrical DNA sequencing with nanopores. Here, we provide a quantitative theoretical evaluation that shows the positive net charges, which accumulate around the pore entrance due to the salt gradient, are responsible for the two observed effects: they reinforce the electric capture field, resulting in promoted molecule capture rate; and they induce cationic electroosmotic flow through the nanopore, thus significantly retarding the motion of the anionic DNA through the nanopore. Our multiphysical simulation results show that, during the polymer trapping stage, the former effect plays the major role, thus resulting in promoted DNA capture rate, while during the nanopore-penetrating stage the latter effect dominates and consequently reduces the DNA translocation speed significantly. Quantitative agreement with experimental results has been reached by further taking nanopore wall surface charges into account.     

Mechanism of How Salt-Gradient-Induced Charges Affect the Translocation of DNA Molecules through a Nanopore

Experiments using nanopores demonstrated that a salt gradient enhances the capture rate of DNA and reduces its translocation speed. These two effects can help to enable electrical DNA sequencing with nanopores. Here, we provide a quantitative theoretical evaluation that shows the positive net charges, which accumulate around the pore entrance due to the salt gradient, are responsible for the two observed effects: they reinforce the electric capture field, resulting in promoted molecule capture rate; and they induce cationic electroosmotic flow through the nanopore, thus significantly retarding the motion of the anionic DNA through the nanopore. Our multiphysical simulation results show that, during the polymer trapping stage, the former effect plays the major role, thus resulting in promoted DNA capture rate, while during the nanopore-penetrating stage the latter effect dominates and consequently reduces the DNA translocation speed significantly. Quantitative agreement with experimental results has been reached by further taking nanopore wall surface charges into account.     

Trapping DNA near a Solid-State Nanopore

We demonstrate that voltage-biased solid-state nanopores can transiently localize DNA in an electrolyte solution. A double-stranded DNA (dsDNA) molecule is trapped when the electric field near the nanopore attracts and immobilizes a nonend segment of the molecule across the nanopore orifice without inducing a folded molecule translocation. In this demonstration of the phenomenon, the ionic current through the nanopore decreases when the dsDNA molecule is trapped by the nanopore. By contrast, a translocating dsDNA molecule under the same conditions causes an ionic current increase. We also present finite-element modeling results that predict this behavior for the conditions of the experiment.     

Automated forward and reverse ratcheting of DNA in a nanopore at 5-Å precision

An emerging DNA sequencing technique uses protein or solid-state pores to analyze individual strands as they are driven in single-file order past a nanoscale sensor. However, uncontrolled electrophoresis of DNA through these nanopores is too fast for accurate base reads. Here, we describe forward and reverse ratcheting of DNA templates through the α-hemolysin nanopore controlled by phi29 DNA polymerase without the need for active voltage control. DNA strands were ratcheted through the pore at median rates of 2.5-40 nucleotides per second and were examined at one nucleotide spatial precision in real time. Up to 500 molecules were processed at ～130 molecules per hour through one pore. The probability of a registry error (an insertion or deletion) at individual positions during one pass along the template strand ranged from 10% to 24.5% without optimization. This strategy facilitates multiple reads of individual strands and is transferable to other nanopore devices for implementation of DNA sequence analysis.     

A parallel graph decomposition algorithm for DNA sequencing with nanopores

MOTIVATION: With the potential availability of nanopore devices that can sense the bases of translocating single-stranded DNA (ssDNA), it is likely that 'reads' of length approximately 10(5) will be available in large numbers and at high speed. We address the problem of complete DNA sequencing using such reads.We assume that approximately 10(2) copies of a DNA sequence are split into single strands that break into randomly sized pieces as they translocate the nanopore in arbitrary orientations. The nanopore senses and reports each individual base that passes through, but all information about orientation and complementarity of the ssDNA subsequences is lost. Random errors (both biological and transduction) in the reads create further complications. RESULTS: We have developed an algorithm that addresses these issues. It can be considered an extreme variation of the well-known Eulerian path approach. It searches over a space of de Bruijn graphs until it finds one in which (a) the impact of errors is eliminated and (b) both possible orientations of the two ssDNA sequences can be identified separately and unambiguously.Our algorithm is able to correctly reconstruct real DNA sequences of the order of 10(6) bases (e.g. the bacterium Mycoplasma pneumoniae) from simulated erroneous reads on a modest workstation in about 1 h. We describe, and give measured timings of, a parallel implementation of this algorithm on the Cray Multithreaded Architecture (MTA-2) supercomputer, whose architecture is ideally suited to this 'unstructured' problem. Our parallel implementation is crucial to the problem of rapidly sequencing long DNA sequences and also to the situation where multiple nanopores are used to obtain a high-bandwidth stream of reads.     

Hydroxymethyluracil modifications enhance the flexibility and hydrophilicity of double-stranded DNA

Oxidation of a DNA thymine to 5-hydroxymethyluracil is one of several recently discovered epigenetic modifications. Here, we report the results of nanopore translocation experiments and molecular dynamics simulations that provide insight into the impact of this modification on the structure and dynamics of DNA. When transported through ultrathin solid-state nanopores, short DNA fragments containing thymine modifications were found to exhibit distinct, reproducible features in their transport characteristics that differentiate them from unmodified molecules. Molecular dynamics simulations suggest that 5-hydroxymethyluracil alters the flexibility and hydrophilicity of the DNA molecules, which may account for the differences observed in our nanopore translocation experiments. The altered physico-chemical properties of DNA produced by the thymine modifications may have implications for recognition and processing of such modifications by regulatory DNA-binding proteins.     

Thermal Motion of DNA in an MspA Pore

We report on an experiment and calculations that determine the thermal motion of a voltage-clamped single-stranded DNA-NeutrAvidin complex in a Mycobacterium smegmatis porin A nanopore. The electric force and diffusion constant of DNA inside a Mycobacterium smegmatis porin A pore were determined to evaluate the thermal position fluctuations of DNA. We show that an out-of-equilibrium state returns to equilibrium so quickly that experiments usually measure a weighted average over the equilibrium position distribution. Averaging over the equilibrium position distribution is consistent with results of state-of-the-art nanopore sequencing experiments. It is shown how a reduction in thermal position fluctuations can be achieved by increasing the electrophoretic force used in nanopore sequencing devices.     

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.     

Translocation of E. coli RecQ helicase on single-stranded DNA

A member of the SF2 family of helicases, Escherichia coli RecQ, is involved in the recombination and repair of double-stranded DNA breaks and single-stranded DNA (ssDNA) gaps. Although the unwinding activity of this helicase has been studied biochemically, the mechanism of translocation remains unclear. To this end, using ssDNA of varying lengths, the steady-state ATP hydrolysis activity of RecQ was analyzed. We find that the rate of ATP hydrolysis increases with DNA length, reaching a maximum specific activity of 38 ± 2 ATP/RecQ/s. Analysis of the rate of ATP hydrolysis as a function of DNA length implies that the helicase has a processivity of 19 ± 6 nucleotides on ssDNA and that RecQ requires a minimal translocation site size of 10 ± 1 nucleotides. Using the T4 phage encoded gene 32 protein (G32P), which binds ssDNA cooperatively, to decrease the lengths of ssDNA gaps available for translocation, we observe a decrease in the rate of ATP hydrolysis activity that is related to lattice occupancy. Analysis of the activity in terms of the average gap sizes available to RecQ on the ssDNA coated with G32P indicates that RecQ translocates on ssDNA on average 46 ± 11 nucleotides before dissociating. Moreover, when bound to ssDNA, RecQ hydrolyzes ATP in a cooperative fashion, with a Hill coefficient of 2.1 ± 0.6, suggesting that at least a dimer is required for translocation on ssDNA. We present a kinetic model for translocation by RecQ on ssDNA based on this characterization.     

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.     

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.     

Direct measurement of single-stranded DNA translocation by PcrA helicase using the fluorescent base analogue 2-aminopurine.

Use of the fluorescent base analogue 2-aminopurine has provided a direct demonstration of the translocation of PcrA helicase toward the 5'-end of single-stranded DNA. Single 2-aminopurine bases are introduced into otherwise standard oligonucleotides and produce a fluorescence signal when PcrA reaches their position. We demonstrate that random binding of PcrA to ssDNA is followed by translocation in an ATP-dependent manner toward the 5'-terminus at 80 bases per second at 20 degrees C. The data also provide information on the kinetics of ssDNA binding to the helicase and of the protein dissociation from the 5'-end of ssDNA. A full kinetic model is presented for ATP-dependent DNA translocation by PcrA helicase.     

T7 DNA helicase: a molecular motor that processively and unidirectionally translocates along single-stranded DNA

DNA helicases are molecular motors that use the energy from NTP hydrolysis to drive the process of duplex DNA strand separation. Here, we measure the translocation and energy coupling efficiency of a replicative DNA helicase from bacteriophage T7 that is a member of a class of helicases that assembles into ring-shaped hexamers. Presteady state kinetics of DNA-stimulated dTTP hydrolysis activity of T7 helicase were measured using a real time assay as a function of ssDNA length, which provided evidence for unidirectional translocation of T7 helicase along ssDNA. Global fitting of the kinetic data provided an average translocation rate of 132 bases per second per hexamer at 18 degrees C. While translocating along ssDNA, T7 helicase hydrolyzes dTTP at a rate of 49 dTTP per second per hexamer, which indicates that the energy from hydrolysis of one dTTP drives unidirectional movement of T7 helicase along two to three bases of ssDNA. One of the features that distinguishes this ring helicase is its processivity, which was determined to be 0.99996, which indicated that T7 helicase travels on an average about 75kb of ssDNA before dissociating. We propose that the ability of T7 helicase to translocate unidirectionally along ssDNA in an efficient manner plays a crucial role in DNA unwinding.     

Distinct complexes of DNA polymerase I (Klenow fragment) for base and sugar discrimination during nucleotide substrate selection

During each catalytic cycle, DNA polymerases select deoxyribonucleoside triphosphate (dNTP) substrates complementary to a templating base with high fidelity from a pool that includes noncomplementary dNTPs and both complementary and noncomplementary ribonucleoside triphosphates (rNTPs). The Klenow fragment of Escherichia coli DNA polymerase I (KF) achieves this through a series of conformational transitions that precede the chemical step of phosphodiester bond formation. Kinetic evidence from fluorescence and FRET experiments indicates that discrimination of the base and sugar moieties of the incoming nucleotide occurs in distinct, sequential steps during the selection pathway. Here we show that KF-DNA complexes formed with complementary rNTPs or with noncomplementary nucleotides can be distinguished on the basis of their properties when captured in an electric field atop the α-hemolysin nanopore. The average nanopore dwell time of KF-DNA complexes increased as a function of complementary rNTP concentration. The increase was less than that promoted by complementary dNTP, indicating that the rNTP complexes are more stable than KF-DNA binary complexes but less stable than KF-DNA-dNTP ternary complexes. KF-DNA-rNTP complexes could also be distinguished from KF-DNA-dNTP complexes on the basis of ionic current amplitude. In contrast to complementary rNTPs, noncomplementary dNTPs and rNTPs diminished the average nanopore dwell time of KF-DNA complexes in a concentration-dependent manner, suggesting that binding of a noncomplementary nucleotide keeps the KF-DNA complex in a less stable state. These results imply that nucleotide selection proceeds through a series of complexes of increasing stability in which substrates with the correct moiety promote the forward transitions.     

Dda helicase tightly couples translocation on single-stranded DNA to unwinding of duplex DNA: Dda is an optimally active helicase

Helicases utilize the energy of ATP hydrolysis to unwind double-stranded DNA while translocating on the DNA. Mechanisms for melting the duplex have been characterized as active or passive, depending on whether the enzyme actively separates the base pairs or simply sequesters single-stranded DNA (ssDNA) that forms due to thermal fraying. Here, we show that Dda translocates unidirectionally on ssDNA at the same rate at which it unwinds double-stranded DNA in both ensemble and single-molecule experiments. Further, the unwinding rate is largely insensitive to the duplex stability and to the applied force. Thus, Dda transduces all of its translocase activity into DNA unwinding activity so that the rate of unwinding is limited by the rate of translocation and that the enzyme actively separates the duplex. Active and passive helicases have been characterized by dividing the velocity of DNA unwinding in base pairs per second (V_un) by the velocity of translocation on ssDNA in nucleotides per second (V_trans). If the resulting fraction is 0.25, then a helicase is considered to be at the lower end of the “active” range. In the case of Dda, the average DNA unwinding velocity was 257 ± 42 bp/s, and the average translocation velocity was 267 ± 15 nt/s. The V_un/V_trans value of 0.96 places Dda in a unique category of being an essentially “perfectly” active helicase.