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.     

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.     

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.     

Rapid discrimination among individual DNA hairpin molecules at single-nucleotide resolution using an ion channel

RNA and DNA strands produce ionic current signatures when driven through an alpha-hemolysin channel by an applied voltage. Here we combine this nanopore detector with a support vector machine (SVM) to analyze DNA hairpin molecules on the millisecond time scale. Measurable properties include duplex stem length, base pair mismatches, and loop length. This nanopore instrument can discriminate between individual DNA hairpins that differ by one base pair or by one nucleotide.     

A toolbox for generating single‐stranded DNA in optical tweezers experiments

Essential genomic transactions such as DNA‐damage repair and DNA replication take place on single‐stranded DNA (ssDNA) or require specific single‐stranded/double‐stranded DNA (ssDNA/dsDNA) junctions (SDSJ). A significant challenge in single‐molecule studies of DNA–protein interactions using optical trapping is the design and generation of appropriate DNA templates. In contrast to dsDNA, only a limited toolbox is available for the generation of ssDNA constructs for optical tweezers experiments. Here, we present several kinds of DNA templates suitable for single‐molecule experiments requiring segments of ssDNA of several kilobases in length. These different biotinylated dsDNA templates can be tethered between optically trapped microspheres and can, by the subsequent use of force‐induced DNA melting, be converted into partial or complete ssDNA molecules. We systematically investigated the time scale and efficiency of force‐induced melting at different ionic strengths for DNA molecules of different sequences and lengths. Furthermore, we quantified the impact of microspheres of different sizes on the lifetime of ssDNA tethers in optical tweezers experiments. Together, these experiments provide deeper insights into the variables that impact the production of ssDNA for single molecules studies and represent a starting point for further optimization of DNA templates that permit the investigation of protein binding and kinetics on ssDNA. © 2013 Wiley Periodicals, Inc. Biopolymers 99:611–620, 2013.     

Sequence-dependent gating of an ion channel by DNA hairpin molecules

DNA hairpins produce ionic current signatures when captured by the alpha-hemolysin nano-scale pore under conditions of single molecule electrophoresis. Gating patterns produced by individual DNA hairpins when captured can be used to distinguish differences of a single base pair or even a single nucleotide [Vercoutere,W.A. et al. (2003) Nucleic Acids Res., 31, 1311–1318]. Here we investigate the mechanism(s) that may account for the ionic current gating signatures. The ionic current resistance profile of conductance states produced by DNA hairpin molecules with 3–12 bp stems showed a plateau in resistance between 10 and 12 bp, suggesting that hairpins with 10–12 bp stems span the pore vestibule. DNA hairpins with 9–12 bp stems produced gating signatures with the same relative conductance states. Systematic comparison of the conductance state dwell times and apparent activation energies for a series of 9–10 bp DNA hairpins suggest that the 3′ and 5′ ends interact at or near the limiting aperture within the vestibule of the alpha-hemolysin pore. The model presented may be useful in predicting and interpreting DNA detection using nanopore detectors. In addition, this well-defined molecular system may prove useful for investigating models of ligand-gated channels in biological membranes.     

Single-molecule studies of the stringency factors and rates governing the polymerization of RecA on double-stranded DNA

RecA is a key protein in homologous recombination. During recombination, one single-stranded DNA (ssDNA) bound to site I in RecA exchanges Watson-Crick pairing with a sequence-matched ssDNA that was part of a double-stranded DNA molecule (dsDNA) bound to site II in RecA. After strand exchange, heteroduplex dsDNA is bound to site I. In vivo, direct polymerization of RecA on dsDNA through site I does not occur, though it does in vitro. The mechanisms underlying the difference have been unclear. We use single-molecule experiments to decouple the two steps involved in polymerization: nucleation and elongation. We find that elongation is governed by a fundamental clock that is insensitive to force and RecA concentration from 0.2 and 6 μM, though rates depend on ionic conditions. Thus, we can probe nucleation site stability by creating nucleation sites at high force and then measuring elongation as a function of applied force. We find that in the presence of ATP hydrolysis a minimum force is required for polymerization. The minimum force decreases with increasing RecA or ATP concentrations. We propose that force reduces the off-rate for nucleation site binding and that nucleation site stability is the stringency factor that prevents in vivo polymerization.     

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.     

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.     

Cheminformatics methods for novel nanopore analysis of HIV DNA termini

BACKGROUND: Channel current feature extraction methods, using Hidden Markov Models (HMMs) have been designed for tracking individual-molecule conformational changes. This information is derived from observation of changes in ionic channel current blockade "signal" upon that molecule's interaction with (and occlusion of) a single nanometer-scale channel in a "nanopore detector". In effect, a nanopore detector transduces single molecule events into channel current blockades. HMM analysis tools described are used to help systematically explore DNA dinucleotide flexibility, with particular focus on HIV's highly conserved (and highly flexible/reactive) viral DNA termini. One of the most critical stages in HIV's attack is the binding between viral DNA and the retroviral integrase, which is influenced by the dynamic-coupling induced high flexibility of a CA/TG dinucleotide positioned precisely two base-pairs from the blunt terminus of the duplex viral DNA. This suggests the study of a family of such CA/TG dinucleotide molecules via nanopore measurement and cheminformatics analysis. RESULTS: HMMs are used for level identification on the current blockades, HMM/EM with boosted variance emissions are used for level projection pre-processing, and time-domain FSAs are used to parse the level-projected waveform for kinetic information. The observed state kinetics of the DNA hairpins containing the CA/TG dinucleotide provides clear evidence for HIV's selection of a peculiarly flexible/interactive DNA terminus.     

纳米孔测序信号处理及其在DNA数据存储的应用

随着高通量测序技术的不断更新,可以在单个分子水平读取核苷酸序列的第三代测序技术迅速发展,纳米孔测序技术是其具有代表性的单分子测序技术,该技术通过检测DNA单链分子穿过纳米孔时引起的跨膜电流信号的变化,实现碱基识别.纳米孔测序仪在便携性、碱基读取速度、测序读段长度等方面较传统的第一代与第二代测序技术都有明显优势.随着纳米...     

Nanopore cheminformatics

A cheminformatics method is described for classification, and biophysical examination, of individual molecules. A novel molecular detector is used--one based on current blockade measurements through a nanometer-scale ion channel (alpha-hemolysin). Classification results are described for blockades caused by DNA molecules in the alpha-hemolysin nanopore detector, with signal analysis and pattern recognition performed using a combination of methods from bioinformatics and machine learning. Due to the size of the alpha-hemolysin protein channel, the blockade events report on one DNA molecule at a time, which enables a variety of reproducible, single-molecule biophysical experiments. To capture the full sensitivity of the nanopore detector's blockade signal, Hidden Markov Models (HMMs) were used with Expectation/Maximization for denoising and for associating a feature vector with the ionic current blockade of each captured DNA molecule. Support Vector Machines (SVMs) that employ novel kernel designs were then used as discriminators. With SVM training performed off-line, and economical HMM processing on-line, blockade classification was possible during capture. HMMs were also used in conjunction with a time-domain finite state automaton (off-line) for feature discovery and kinetics analysis. Analysis of the DNA data indicates a variety of binding (DNA-protein), fraying, and conformational shifts that are consistent with data obtained from thermodynamic analyses (melting curves), X-ray crystallography, and NMR studies. The software tools are designed for analysis of generic blockades in ionic channels, including those in other biological pore-forming toxins, other biological channels in general, and semiconductor-based channels.     

Direct observation of translocation in individual DNA polymerase complexes

Complexes of phi29 DNA polymerase and DNA fluctuate on the millisecond time scale between two ionic current amplitude states when captured atop the α-hemolysin nanopore in an applied field. The lower amplitude state is stabilized by complementary dNTP and thus corresponds to complexes in the post-translocation state. We have demonstrated that in the upper amplitude state, the DNA is displaced by a distance of one nucleotide from the post-translocation state. We propose that the upper amplitude state corresponds to complexes in the pre-translocation state. Force exerted on the template strand biases the complexes toward the pre-translocation state. Based on the results of voltage and dNTP titrations, we concluded through mathematical modeling that complementary dNTP binds only to the post-translocation state, and we estimated the binding affinity. The equilibrium between the two states is influenced by active site-proximal DNA sequences. Consistent with the assignment of the upper amplitude state as the pre-translocation state, a DNA substrate that favors the pre-translocation state in complexes on the nanopore is a superior substrate in bulk phase for pyrophosphorolysis. There is also a correlation between DNA sequences that bias complexes toward the pre-translocation state and the rate of exonucleolysis in bulk phase, suggesting that during DNA synthesis the pathway for transfer of the primer strand from the polymerase to exonuclease active site initiates in the pre-translocation state.     

Characterizing and controlling the motion of ssDNA in a solid-state nanopore

Sequencing DNA in a synthetic solid-state nanopore is potentially a low-cost and high-throughput method. Essential to the nanopore-based DNA sequencing method is the ability to control the motion of a single-stranded DNA (ssDNA) molecule at single-base resolution. Experimental studies showed that the average translocation speed of DNA driven by a biasing electric field can be affected by ionic concentration, solvent viscosity, or temperature. Even though it is possible to slow down the average translocation speed, instantaneous motion of DNA is too diffusive to allow each DNA base to stay in front of a sensor site for its measurement. Using extensive all-atom molecular dynamics simulations, we study the diffusion constant, friction coefficient, electrophoretic mobility, and effective charge of ssDNA in a solid-state nanopore. Simulation results show that the spatial fluctuation of ssDNA in 1 ns is comparable to the spacing between neighboring nucleotides in ssDNA, which makes the sensing of a DNA base very difficult. We demonstrate that the recently proposed DNA transistor could potentially solve this problem by electrically trapping ssDNA inside the DNA transistor and ratcheting ssDNA base-by-base in a biasing electric field. When increasing the biasing electric field, we observed that the translocation of ssDNA changes from ratcheting to steady-sliding. The simulated translocation of ssDNA in the DNA transistor was theoretically characterized using Fokker-Planck analysis.     

Single-molecule manipulation of double-stranded DNA using optical tweezers: interaction studies of DNA with RecA and YOYO-1.

By using optical tweezers and a specially designed flow cell with an integrated glass micropipette, we constructed a setup similar to that of Smith et al. (Science 271:795-799, 1996) in which an individual double-stranded DNA (dsDNA) molecule can be captured between two polystyrene beads. The first bead is immobilized by the optical tweezers and the second by the micropipette. Movement of the micropipette allows manipulation and stretching of the DNA molecule, and the force exerted on it can be monitored simultaneously with the optical tweezers. We used this setup to study elongation of dsDNA by RecA protein and YOYO-1 dye molecules. We found that the stability of the different DNA-ligand complexes and their binding kinetics were quite different. The length of the DNA molecule was extended by 45% when RecA protein was added. Interestingly, the speed of elongation was dependent on the external force applied to the DNA molecule. In experiments in which YOYO-1 was added, a 10-20% extension of the DNA molecule length was observed. Moreover, these experiments showed that a change in the applied external force results in a time-dependent structural change of the DNA-YOYO-1 complex, with a time constant of approximately 35 s (1/e2). Because the setup provides an oriented DNA molecule, we determined the orientation of the transition dipole moment of YOYO-1 within DNA by using fluorescence polarization. The angle of the transition dipole moment with respect to the helical axis of the DNA molecule was 69 degrees +/- 3.     

Effect of temperature and ionic strength on the dissociation kinetics and lifetime of PNA-DNA triplexes

Dissociation kinetics of triplexes formed by molecules of peptide nucleic acid (PNA) and DNA have been studied. The complexes consisted of oligomeric PNA containing 10 thymine bases and the dA(10) target incorporated in single-stranded (ssDNA) or double-stranded DNA (dsDNA). Their dissociation was followed by means of the gel mobility shift assay at various temperatures and sodium ion concentrations. In all experiments, the dissociation kinetics of triplexes were exponential; the effective lifetime of a triplex, tau, depended on temperature in accordance with the Arrhenius law. The tau values for T(10) PNA complexes with ss- and dsDNA were equal within the accuracy of experiments. The activation energy, U, value for T(10) PNA-DNA complexes did not change when the NaCl concentration was increased from 50 to 200 or 600 mM. Conversely, the tau values decreased with the increase in NaCl concentration. The equal lifetimes of the T(10) PNA-DNA triplexes containing ss- and dsDNA suggest that the loop formed in dsDNA does not noticeably affect the triplex structure. The decrease in the triplex lifetime tau with an increase in ionic strength was accounted for by the fact that the PNA backbone is neutral. The lack of relationship between the activation energy of dissociation and salt concentration suggests that the dissociation enthalpy does not depend on the ionic strength. Thus, the effect of ionic strength on the lifetime is entropic by its nature. Contrary to this, for complexes of ssDNA with bis-PNA 1743, which also consists of 10 thymine bases but contains 2 additional positive charges inside the sequence in 1 of the PNA arms, an increase of the dissociation enthalpy at low salt concentration was observed. We suggest that this effect is a result of a direct electrostatic interaction of the positive charges of the PNA with the DNA backbone. Finally, our results allow an estimate of the lifetime of a 10-mer triplex invasion complex in dsDNA at 37 degrees C in excess of several hundred days.     

Measuring single-molecule DNA hybridization by active control of DNA in a nanopore

We present a novel application of active voltage control of DNA captured in a nanopore to regulate the amount of time the DNA is available to molecules in the bulk phase that bind to the DNA. In this work, the control method is used to measure hybridization between a single molecule of DNA captured in a nanopore and complementary oligonucleotides in the bulk phase. We examine the effect of oligonucleotide length on hybridization, and the effect of DNA length heterogeneity on the measurements. Using a mathematical model, we are able to deduce the binding rate of complementary oligonucleotides, even when DNA samples in experiments are affected by heterogeneity in length. We analyze the lifetime distribution of DNA duplexes that are formed in the bulk phase and then pulled against the pore by reversing the voltage. The lifetime distribution reveals several dissociation modes. It remains to be resolved whether these dissociation modes are due to DNA heterogeneity or correspond to different states of duplex DNA. The control method is unique in its ability to detect single-molecule complex assembly in the bulk phase, free from external force and with a broad (millisecond-to-second) temporal range.     

The interaction of the A and A* proteins of bacteriophage phi X174 with single-stranded and double-stranded phi X DNA in vitro

The binding of the bacteriophage phi X 174-coded A and A* proteins to single-stranded (ssDNA) and double-stranded (dsDNA ) phi X DNA was studied by electron microscopy. The interaction of the A* protein with ssDNA and dsDNA was also studied by sedimentation velocity centrifugation. It was shown that the binding of the A and A* proteins to ssDNA occurs in a non-cooperative manner and requires no or very little sequence specificity under the conditions used here. Both protein-ssDNA complexes have the same compact structure caused by intrastrand cross-linking through the interaction of protein molecules with separate parts of the ssDNA molecule. The A protein does not bind to phi X dsDNA in the absence of divalent cations. The A* protein does bind to dsDNA, although it has a strong preference for binding to ssDNA. The structure of the A* protein-dsDNA complexes is different from that of the A* protein-ssDNA complexes, as the former have a rosette-like structure caused by protein-protein interactions. High ionic strengths favour the formation of large condensed aggregates.