Molecular dynamics simulations of DNA within a nanopore: arginine-phosphate tethering and a binding/sliding mechanism for translocation

Protein nanopores show great potential as low-cost detectors in DNA sequencing devices. To date, research has largely focused on the staphylococcal pore α-hemolysin (αHL). In the present study, we have developed simplified models of the wild-type αHL pore and various mutants in order to study the translocation dynamics of single-stranded DNA under the influence of an applied electric field. The model nanopores reflect the experimentally measured conductance values in planar lipid bilayers. We show that interactions between rings of cationic amino acids and DNA backbone phosphates result in metastable tethering of nucleic acid molecules within the pore, leading us to propose a "binding and sliding" mechanism for translocation. We also observe folding of DNA into nonlinear conformational intermediates during passage through the confined nanopore environment. Despite adopting nonlinear conformations, the DNA hexamer always exits the pore in the same orientation as it enters (3' to 5') in our simulations. The observations from our simulations help to rationalize experimentally determined trends in residual current and translocation efficiency for αHL and its mutants.     

Urea Facilitates the Translocation of Single-Stranded DNA and RNA Through the α-Hemolysin Nanopore

The staphylococcal α-hemolysin (αHL) protein nanopore is under investigation as a fast, cheap detector for nucleic acid analysis and sequencing. Although discrimination of all four bases of DNA by the αHL pore has been demonstrated, analysis of single-stranded DNAs and RNAs containing secondary structure mediated by basepairing is prevented because these nucleic acids cannot be translocated through the pore. Here, we show that a structured 95-nucleotide single-stranded DNA and its RNA equivalent are translocated through the αHL pore in the presence of 4 M urea, a concentration that denatures the secondary structure of the polynucleotides. The αHL pore is functional even in 7 M urea, and therefore it is easily stable enough for analyses of challenging DNA and RNA species.     

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.     

DNA Translocation and Unzipping through a Nanopore: Some Geometrical Effects

This article explores the role of some geometrical factors on the electrophoretically driven translocations of macromolecules through nanopores. In the case of asymmetric pores, we show how the entry requirements and the direction of translocation can modify the information content of the blocked ionic current as well as the transduction of the electrophoretic drive into a mechanical force. To address these effects we studied the translocation of single-stranded DNA through an asymmetric α-hemolysin pore. Depending on the direction of the translocation, we measure the capacity of the pore to discriminate between both DNA orientations. By unzipping DNA hairpins from both sides of the pores we show that the presence of single strand or double strand in the pore can be discriminated based on ionic current levels. We also show that the transduction of the electrophoretic drive into a denaturing mechanical force depends on the local geometry of the pore entrance. Eventually we discuss the application of this work to the measurement of energy barriers for DNA unzipping as well as for protein binding and unfolding.     

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.     

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.     

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.     

Nucleotide discrimination with DNA immobilized in the MspA nanopore

Nanopore sequencing has the potential to become a fast and low-cost DNA sequencing platform. An ionic current passing through a small pore would directly map the sequence of single stranded DNA (ssDNA) driven through the constriction. The pore protein, MspA, derived from Mycobacterium smegmatis, has a short and narrow channel constriction ideally suited for nanopore sequencing. To study MspA's ability to resolve nucleotides, we held ssDNA within the pore using a biotin-NeutrAvidin complex. We show that homopolymers of adenine, cytosine, thymine, and guanine in MspA exhibit much larger current differences than in α-hemolysin. Additionally, methylated cytosine is distinguishable from unmethylated cytosine. We establish that single nucleotide substitutions within homopolymer ssDNA can be detected when held in MspA's constriction. Using genomic single nucleotide polymorphisms, we demonstrate that single nucleotides within random DNA can be identified. Our results indicate that MspA has high signal-to-noise ratio and the single nucleotide sensitivity desired for nanopore sequencing devices.     

How Nanopore Translocation Experiments Can Measure RNA Unfolding

Electrokinetic translocation of biomolecules through solid-state nanopores represents a label-free single-molecule technique that may be used to measure biomolecular structure and dynamics. Recent investigations have attempted to distinguish individual transfer RNA (tRNA) species based on the associated pore translocation times, ion-current noise, and blockage currents. By manufacturing sufficiently smaller pores, each tRNA is required to undergo a deformation to translocate. Accordingly, differences in nanopore translocation times and distributions may be used to infer the mechanical properties of individual tRNA molecules. To bridge our understanding of tRNA structural dynamics and nanopore measurements, we apply molecular dynamics simulations using a simplified “structure-based” energetic model. Calculating the free-energy landscape for distinct tRNA species implicates transient unfolding of the terminal RNA helix during nanopore translocation. This provides a structural and energetic framework for interpreting current experiments, which can aid the design of methods for identifying macromolecules using nanopores.     

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.     

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.     

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.     

Peptide Nanopores and Lipid Bilayers: Interactions by Coarse-Grained Molecular-Dynamics Simulations

A set of 49 protein nanopore-lipid bilayer systems was explored by means of coarse-grained molecular-dynamics simulations to study the interactions between nanopores and the lipid bilayers in which they are embedded. The seven nanopore species investigated represent the two main structural classes of membrane proteins (α-helical and β-barrel), and the seven different bilayer systems range in thickness from ∼28 to ∼43 Å. The study focuses on the local effects of hydrophobic mismatch between the nanopore and the lipid bilayer. The effects of nanopore insertion on lipid bilayer thickness, the dependence between hydrophobic thickness and the observed nanopore tilt angle, and the local distribution of lipid types around a nanopore in mixed-lipid bilayers are all analyzed. Different behavior for nanopores of similar hydrophobic length but different geometry is observed. The local lipid bilayer perturbation caused by the inserted nanopores suggests possible mechanisms for both lipid bilayer-induced protein sorting and protein-induced lipid sorting. A correlation between smaller lipid bilayer thickness (larger hydrophobic mismatch) and larger nanopore tilt angle is observed and, in the case of larger hydrophobic mismatches, the simulated tilt angle distribution seems to broaden. Furthermore, both nanopore size and key residue types (e.g., tryptophan) seem to influence the level of protein tilt, emphasizing the reciprocal nature of nanopore-lipid bilayer interactions.     

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.     

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.     

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.     

How B-DNA Dynamics Decipher Sequence-Selective Protein Recognition

The rules governing sequence-specific DNA–protein recognition are under a long-standing debate regarding the prevalence of base versus shape readout mechanisms to explain sequence specificity and of the conformational selection versus induced fit binding paradigms to explain binding-related conformational changes in DNA. Using a combination of atomistic simulations on a subset of representative sequences and mesoscopic simulations at the protein–DNA interactome level, we demonstrate the prevalence of the shape readout model in determining sequence-specificity and of the conformational selection paradigm in defining the general mechanism for binding-related conformational changes in DNA. Our results suggest that the DNA uses a double mechanism to adapt its structure to the protein: it moves along the easiest deformation modes to approach the bioactive conformation, while final adjustments require localized rearrangements at the base-pair step and backbone level. Our study highlights the large impact of B-DNA dynamics in modulating DNA–protein binding.