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.     

Direct Detection of Unnatural DNA Nucleotides dNaM and d5SICS using the MspA Nanopore

Malyshev et al. showed that the four-letter genetic code within a living organism could be expanded to include the unnatural DNA bases dNaM and d5SICS. However, verification and detection of these unnatural bases in DNA requires new sequencing techniques. Here we provide proof of concept detection of dNaM and d5SICS in DNA oligomers via nanopore sequencing using the nanopore MspA. We find that both phi29 DNA polymerase and Hel308 helicase are capable of controlling the motion of DNA containing dNaM and d5SICS through the pore and that single reads are sufficient to detect the presence and location of dNaM and d5SICS within single molecules.     

MspA nanopores from subunit dimers

Mycobacterium smegmatis porin A (MspA) forms an octameric channel and represents the founding member of a new family of pore proteins. Control of subunit stoichiometry is important to tailor MspA for nanotechnological applications. In this study, two MspA monomers were connected by linkers ranging from 17 to 62 amino acids in length. The oligomeric pore proteins were purified from M. smegmatis and were shown to form functional channels in lipid bilayer experiments. These results indicated that the peptide linkers did not prohibit correct folding and localization of MspA. However, expression levels were reduced by 10-fold compared to wild-type MspA. MspA is ideal for nanopore sequencing due to its unique pore geometry and its robustness. To assess the usefulness of MspA made from dimeric subunits for DNA sequencing, we linked two M1-MspA monomers, whose constriction zones were modified to enable DNA translocation. Lipid bilayer experiments demonstrated that this construct also formed functional channels. Voltage gating of MspA pores made from M1 monomers and M1-M1 dimers was identical indicating similar structural and dynamic channel properties. Glucose uptake in M. smegmatis cells lacking porins was restored by expressing the dimeric mspA M1 gene indicating correct folding and localization of M1-M1 pores in their native membrane. Single-stranded DNA hairpins produced identical ionic current blockades in pores made from monomers and subunit dimers demonstrating that M1-M1 pores are suitable for DNA sequencing. This study provides the proof of principle that production of single-chain MspA pores in M. smegmatis is feasible and paves the way for generating MspA pores with altered stoichiometries. Subunit dimers enable better control of the chemical and physical properties of the constriction zone of MspA. This approach will be valuable both in understanding transport across the outer membrane in mycobacteria and in tailoring MspA for nanopore sequencing of DNA.     

Control of subunit stoichiometry in single-chain MspA nanopores

Transmembrane protein channels enable fast and highly sensitive detection of single molecules. Nanopore sequencing of DNA was achieved using an engineered Mycobacterium smegmatis porin A (MspA) in combination with a motor enzyme. Due to its favorable channel geometry, the octameric MspA pore exhibits the highest current level compared with other pore proteins. To date, MspA is the only protein nanopore with a published record of DNA sequencing. While widely used in commercial devices, nanopore sequencing of DNA suffers from significant base-calling errors due to stochastic events of the complex DNA-motor-pore combination and the contribution of up to five nucleotides to the signal at each position. Different mutations in specific subunits of a pore protein offer an enormous potential to improve nucleotide resolution and sequencing accuracy. However, individual subunits of MspA and other oligomeric protein pores are randomly assembled in vivo and in vitro, preventing the efficient production of designed pores with different subunit mutations. In this study, we converted octameric MspA into a single-chain pore by connecting eight subunits using peptide linkers. Lipid bilayer experiments demonstrated that single-chain MspA formed membrane-spanning channels and discriminated all four nucleotides identical to MspA produced from monomers in DNA hairpin experiments. Single-chain constructs comprising three, five, six, and seven connected subunits assembled to functional channels, demonstrating a remarkable plasticity of MspA to different subunit stoichiometries. Thus, single-chain MspA constitutes a new milestone in the optimization of MspA as a biosensor for DNA sequencing and many other applications by enabling the production of pores with distinct subunit mutations and pore diameters.     

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.     

Structures of phi29 DNA polymerase complexed with substrate: the mechanism of translocation in B-family polymerases

Replicative DNA polymerases (DNAPs) move along template DNA in a processive manner. The structural basis of the mechanism of translocation has been better studied in the A-family of polymerases than in the B-family of replicative polymerases. To address this issue, we have determined the X-ray crystal structures of phi29 DNAP, a member of the protein-primed subgroup of the B-family of polymerases, complexed with primer-template DNA in the presence or absence of the incoming nucleoside triphosphate, the pre- and post-translocated states, respectively. Comparison of these structures reveals a mechanism of translocation that appears to be facilitated by the coordinated movement of two conserved tyrosine residues into the insertion site. This differs from the mechanism employed by the A-family polymerases, in which a conserved tyrosine moves into the templating and insertion sites during the translocation step. Polymerases from the two families also interact with downstream single-stranded template DNA in very different ways.     

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.     

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.     

Measuring cation dependent DNA polymerase fidelity landscapes by deep sequencing

High-throughput recording of signals embedded within inaccessible micro-environments is a technological challenge. The ideal recording device would be a nanoscale machine capable of quantitatively transducing a wide range of variables into a molecular recording medium suitable for long-term storage and facile readout in the form of digital data. We have recently proposed such a device, in which cation concentrations modulate the misincorporation rate of a DNA polymerase (DNAP) on a known template, allowing DNA sequences to encode information about the local cation concentration. In this work we quantify the cation sensitivity of DNAP misincorporation rates, making possible the indirect readout of cation concentration by DNA sequencing. Using multiplexed deep sequencing, we quantify the misincorporation properties of two DNA polymerases - Dpo4 and Klenow exo(-) - obtaining the probability and base selectivity of misincorporation at all positions within the template. We find that Dpo4 acts as a DNA recording device for Mn(2+) with a misincorporation rate gain of ～2%/mM. This modulation of misincorporation rate is selective to the template base: the probability of misincorporation on template T by Dpo4 increases >50-fold over the range tested, while the other template bases are affected less strongly. Furthermore, cation concentrations act as scaling factors for misincorporation: on a given template base, Mn(2+) and Mg(2+) change the overall misincorporation rate but do not alter the relative frequencies of incoming misincorporated nucleotides. Characterization of the ion dependence of DNAP misincorporation serves as the first step towards repurposing it as a molecular recording device.     

Dynamics of Translocation and Substrate Binding in Individual Complexes Formed with Active Site Mutants of Φ29 DNA Polymerase

The Φ29 DNA polymerase (DNAP) is a processive B-family replicative DNAP. Fluctuations between the pre-translocation and post-translocation states can be quantified from ionic current traces, when individual Φ29 DNAP-DNA complexes are held atop a nanopore in an electric field. Based upon crystal structures of the Φ29 DNAP-DNA binary complex and the Φ29 DNAP-DNA-dNTP ternary complex, residues Tyr-226 and Tyr-390 in the polymerase active site were implicated in the structural basis of translocation. Here, we have examined the dynamics of translocation and substrate binding in complexes formed with the Y226F and Y390F mutants. The Y226F mutation diminished the forward and reverse rates of translocation, increased the affinity for dNTP in the post-translocation state by decreasing the dNTP dissociation rate, and increased the affinity for pyrophosphate in the pre-translocation state. The Y390F mutation significantly decreased the affinity for dNTP in the post-translocation state by decreasing the association rate ∼2-fold and increasing the dissociation rate ∼10-fold, implicating this as a mechanism by which this mutation impedes DNA synthesis. The Y390F dissociation rate increase is suppressed when complexes are examined in the presence of Mn²⁺ rather than Mg²⁺. The same effects of the Y226F or Y390F mutations were observed in the background of the D12A/D66A mutations, located in the exonuclease active site, ∼30 Å from the polymerase active site. Although translocation rates were unaffected in the D12A/D66A mutant, these exonuclease site mutations caused a decrease in the dNTP dissociation rate, suggesting that they perturb Φ29 DNAP interdomain architecture.     

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.     

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.     

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.     

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.     

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.     

DNA conformational change induced by the bacteriophage phi 29 connector.

Translocation of viral DNA inwards and outwards of the capsid of double-stranded DNA bacteriophages occurs through the connector, a key viral structure that is known to interact with DNA. It is shown here that phage phi 29 connector binds both linear and circular double-stranded DNA. However, DNA-mediated protection of phi 29 connectors against Staphylococcus aureus endoprotease V8 digestion suggests that binding to linear DNA is more stable than to circular DNA. Endoprotease V8-protection assays also suggest that the length of the linear DNA required to produce a stable phi 29 connector-DNA interaction is, at least, twice longer than the phi 29 connector channel. This result is confirmed by experiments of phi 29 connector-protection of DNA against DNase I digestion. Furthermore, DNA circularization assays indicate that phi 29 connectors restrain negative supercoiling when bound to linear DNA. This DNA conformational change is not observed upon binding to circular DNA and it could reflect the existence of some left-handed DNA coiling or DNA untwisting inside of the phi 29 connector channel.     

TempliPhi,phi29 DNA polymerase based rolling circle amplification of templates for DNA sequencing

We have developed a novel, isothermal DNA amplification strategy that employs phi29 DNA polymerase and rolling circle amplification to generate high-quality templates for DNA sequencing reactions. The TempliPhi DNA amplification kits take advantage of the fact that cloned DNA is typically obtained in circular vectors, which are readily replicated in vitro using phi29 DNA polymerase by a rolling circle mechanism. This single subunit, proofreading DNA polymerase has excellent processivity and strand displacement properties for generation of multiple, tandem double-stranded copies of the circular DNA, generating as much as 10(7)-fold amplification. Large amounts of product (1-3 microg) can be obtained in as little as 4 hours. Input DNA can be as little as 0.01 ng of purified plasmid DNA, a single bacterial colony, or a 1 microL of a saturated overnight culture. Additionally, the presence of an associated proof reading function within the phi29 DNA polymerase ensures high-fidelity amplification. Once completed, the product DNA can be used directly in sequencing reactions. Additionally, the properties of phi29 DNA polymerase and its use in applications such as amplification ofhuman genomic DNA for genotyping studies is discussed.     

Amplification of uncultured single-stranded DNA viruses from rice paddy soil

Viruses are known to be the most numerous biological entities in soil; however, little is known about their diversity in this environment. In order to explore the genetic diversity of soil viruses, we isolated viruses by centrifugation and sequential filtration before performing a metagenomic investigation. We adopted multiple-displacement amplification (MDA), an isothermal whole-genome amplification method with phi29 polymerase and random hexamers, to amplify viral DNA and construct clone libraries for metagenome sequencing. By the MDA method, the diversity of both single-stranded DNA (ssDNA) viruses and double-stranded DNA viruses could be investigated at the same time. On the contrary, by eliminating the denaturing step in the MDA reaction, only ssDNA viral diversity could be explored selectively. Irrespective of the denaturing step, more than 60% of the soil metagenome sequences did not show significant hits (E-value criterion, 0.001) with previously reported viral sequences. Those hits that were considered to be significant were also distantly related to known ssDNA viruses (average amino acid similarity, approximately 34%). Phylogenetic analysis showed that replication-related proteins (which were the most frequently detected proteins) related to those of ssDNA viruses obtained from the metagenomic sequences were diverse and novel. Putative circular genome components of ssDNA viruses that are unrelated to known viruses were assembled from the metagenomic sequences. In conclusion, ssDNA viral diversity in soil is more complex than previously thought. Soil is therefore a rich pool of previously unknown ssDNA viruses.     

Chimeric Phi29 DNA polymerase with helix–hairpin–helix motifs shows enhanced salt tolerance and replication performance

Phi29 DNA polymerase (Phi29 Pol) has been successfully applied in DNA nanoball-based sequencing, real-time DNA sequencing from single polymerase molecules and nanopore sequencing employing the sequencing by synthesis (SBS) method. Among these, polymerase-assisted nanopore sequencing technology analyses nucleotide sequences as a function of changes in electrical current. This ionic, current-based sequencing technology requires polymerases to perform replication at high salt concentrations, for example 0.3 M KCl. Nonetheless, the salt tolerance of wild-type Phi29 Pol is relatively low. Here, we fused helix–hairpin–helix (HhH)₂ domains E-L (eight repeats in total) of topoisomerase V (Topo V) from the hyperthermophile Methanopyrus kandleri to the Phi29 Pol COOH terminus, designated Phi29EL DNA polymerase (Phi29EL Pol). Domain fusion increased the overall enzyme replication efficiency by fourfold. Phi29EL Pol catalysed rolling circle replication in a broader range of salt concentrations than did Phi29 Pol, extending the KCl concentration range for activity up to 0.3 M. In addition, the mutation of Glu³⁷⁵ to Ser or Gln increased Phi29EL Pol activity in the presence of KCl. In this work, we produced a salt-tolerant Phi29 Pol derivative by means of (HhH)₂ domain insertion. The multiple advantages of this insertion make it a good substitute for Phi29 Pol, especially for use in nanopore sequencing or other circumstances that require high salt concentrations.     

RNA-mediated specificity of DNA packaging into hybrid lambda/phi 29 proheads.

A small RNA (pRNA, 174 nt) is known to be essential for DNA packaging in bacteriophage phi 29. However, in an in vitro DNA packaging system based on hybrid lambda/phi 29 proheads (made up of head proteins from phage lambda and connectors from phage phi 29), the specificity of DNA packaging is lost, and different RNA molecules fulfil the requirements for DNA packaging, albeit with less efficiency than phi 29 pRNA. Competition assays with RNAs from different sources have shown that phi 29 connectors bind preferentially pRNA. An increase in the efficiency of phi 29 DNA packaging into hybrid proheads induced by phi 29 pRNA is observed because, when phi 29 pRNA is incubated with hybrid proheads, phi 29 DNA is packaged more efficiently than other DNAs of similar length. Furthermore, when hybrid proheads carrying phi 29 pRNA are incubated with a mixture of DNAs from different sources, phi 29 DNA is selectively packaged, thus indicating that phi 29 pRNA determines the specificity of DNA packaging.