# Patent

## Title
Nanopore-facilitated single molecule detection of nucleic acid

## Abstract

The present invention provides a new and improved oligonucleotide detection method based on the nanopore technology with a probe containing a complementary sequence to the target oligonucleotide and a terminal extension at the probe's 3′ terminus, 5′ terminus, or both termini. The improved nanopore sensor with the probe enables sensitive, selective, and direct detection, differentiation and quantification of target oligonucleotides such as miRNAs. The inventive detection method may also be employed as a non-invasive and cost-effective diagnostic method for cancer detection based on miRNA levels in the patient's blood sample.

## Background
RELATED APPLICATIONS
This application is a National Stage application of International Patent Application No PCT/US2011/044082, filed Jul. 14, 2011, which is incorporated herein by reference in its entirety and which claims the benefit of U.S. Provisional Patent Application No. 61/399,578, filed Jul. 14, 2010, which is incorporated herein by reference in its entirety.


GRANT STATEMENT
This invention was made with government support under Grant No. GM079613 awarded by the National Institute for Health. The government has certain rights in the invention.     


FIELD OF INVENTION
This product relates to a method/apparatus of single-molecule detection, more specifically, to a method/system for quantitative detection of single strand nucleic acids, such as microRNAs, employing an ultrasensitive, low noise nanopore-based single-molecule technology.
SEQUENCE LISTING STATEMENT
The sequence listing that is contained in the file named “52553_110702_ST25.txt ”, which is 12,350 bytes in size (measured in operating system MS-Windows), created on Jan. 14, 2013, is filed herewith by electronic submission and incorporated herein by reference in its entirety.
BACKGROUND OF INVENTION
MicroRNAs.
MicroRNAs (miRNAs) are a class of short ( ̃18-24 nucleotides) noncoding RNAs that regulate gene expression at the post-transcriptional level2. Depending on the degree of homology to their target sequences, miRNA binding induces either translational repression or cleavage of target mRNAs2. As powerful gene regulators, miRNAs play important roles in development, cell differentiation, and regulation of cell cycle, apoptosis and signaling pathways2,3. Aberrant expression of miRNAs has been found in all types of tumors4,5; the different cancer types have distinct miRNA expression profiles6. Specific miRNAs including some miRNA families containing a few nucleotide differences are constantly released from the primary tumor into blood stream and present in an incredible stable form7. Recent studies demonstrated that circulating miRNAs are enveloped inside exosomal vesicles and can be transferable and functional in the recipient cells8,9. Thus, detection of tumor specific circulating miRNAs provides a powerful tool for early diagnosis, staging, and monitoring of cancer cells10.
MiRNA Detection.
Several technologies including reverse transcription real-time PCR (RT-qPCR) and microarray for miRNA detection have been developed11-13. Each technology has its own advantages, but limitations include requiring enzymatic amplification and semi-quantitative results14. In particular the short miRNA sequences make it difficult to selectively design the primers or probes, resulting in cross-hybridization and low selectivity. This is especially true when the miRNAs contain a few or a single nucleotide difference in a miRNA family. Emerging techniques based on colorimetry, bioluminescence, enzyme turnovers and electrochemistry have been proposed, and nanoparticles and molecular beacon have been applied to miRNA detection with high sensitivity and selectivity (review14). But the intrinsic versatility needs to be improved. Recently, the integration of single-molecule fluorescence and lock-nucleic acids (LNA)15 probes provided a sensitive method for miRNA profiling in tissue samples16, though this method requires expensive instrument.
Nanopore Single Molecule Detection.
In a nanometer-scaled pore structure, the ion current becomes very sensitive to the presence, location and conformation of single target molecules occupying the ion pathway17. This sensitivity allows “visualizing” single molecules, elucidating their kinetics from characteristic change in the pore conductance, and quantifying the target from the occurrence of single molecule signature events. Nanopores have been developed as receptive single molecule detectors for broad biotechnological applications (reviews17-19). The nanopore is also recognized as one of the next generations of DNA sequencing technologies20,21. For example, the 2-nm nanopore, α-hemolysin transmembrane protein pore, allows rapid translocation of single-stranded oligonucleotide, which has been well characterized for DNA sequencmg22-27. However, the molecular translocation-based sensing mode is not suitable for miRNA detection because the sequences of all mature miRNAs are short (18-24 nt), and when traversing the nanopore, the current signals by different miRNAs are indistinguishable.
Therefore, there is a need to provide a new miRNA detection method based on nano-scale pore structure with improved sensitivity, speedy process, and cost efficiency.
SUMMARY OF INVENTION
In one aspect of the invention, a new and improved nanopore-based sensing system for detection and differentiation of single strand oligonucleotides, such as miRNAs, is described. The inventive system for detecting a target single strand oligonucleotide comprises 1) a nanopore, 2) a power source providing sufficient voltage to induce unzipping, 3) a probe with its center domain complementary to the target oligonucleotide, whereas the unzipping of the hybrid of target oligonucleotide and the probe in the nanopore induces certain identifiable current signal changes, and 4) means for detecting the current signal changes. The inventive probe further comprises at least one signal tag at its 3′ or 5′ terminal (or both). The signal tag may be of any charged single chain molecule with sufficient length to assist the unzipping translocation through the nanopore driven by the voltage. For example, the signal tag may be oligonucleotides such as poly(dC)n, poly(dA)n, and poly(dT)n, or charged polypeptides.
In another aspect of the invention, a new and improved method based on nanopore technology for detecting and differentiating single strand oligonucleotide is described. The inventive method detects the current changes triggered by the unzipping of the hybrid of the target oligonucleotide and its probe in a nanopore. The inventive method includes the step of 1) mixing the target oligonucleotide with a pre-designed probe, which has its central domain matching the target sequence and a charged single chain molecule tagged to at least one of its 3′ and 5′ terminals, to produce a sample mixture, 2) loading the mixture into the cis chamber of a nanopore system, and a voltage is added from the trans chamber, and 3) recording current output for a pre-determined time period.
In yet another aspect of the invention, a new and improved method for detecting and monitoring cancer-related miRNAs in patients' blood sample is described. The inventive method includes the steps of 1) mixing the total plasma RNAs extracted from a patient's blood sample with the miRNA probe that contains the complementary sequence to the targeting miRNA and a signal tag at the probe's 3 ‘-terminal, 5’-terminal, or both, 2) adding the mixture into a nanopore chamber with a preselected voltage, and 3) monitoring and analyzing the signature events in the output current traces that serves as an electrical marker for single miRNA molecule recognition.
In certain embodiments, a probe molecule for detecting of a single strand oligonucleotide, such as miRNA, using a nanopore comprising: 1) a center domain with a complementary sequence to the target oligonucleotide, and 2) a terminal extension tagged to at least one of the center domain's 3′ or 5′ terminals is provided. In certain embodiments, the terminal extension is a charged chain molecule. In certain embodiments, the terminal extension is a charged polypeptide. In certain embodiments, terminal extension is a charged polymer.
In certain embodiments, the invention provides a method of detecting single strand oligonucleotide with a dual-compartment nanopore system, whereas the system includes a cis compartment and a trans compartment divided by a partition with an opening at its center region, recording solution filling both compartments and a lipid bilayer formed at the opening on the partition, a nanopore plugging through the lipid bilayer bridging the cis and trans chamber, a voltage loaded upon the system via a pair of electrodes each extruding from the cis or trans compartment, and a current detector monitoring the current changes, includes the steps of 1) mixing the target oligonucleotide with a pre-designed probe, which has its central domain matching the target sequence and a charged single chain molecule tagged to at least one of its 3′ and 5′ terminals, to produce a sample mixture, 2) loading the mixture into the cis compartment, 3) providing the system with a pre-determined voltage, and 4) recording current output for a pre-determined time period. In certain embodiments of the methods, the recording step can further comprise the step of analyzing the current change induced by the hybrid of the target oligonucleotide and the probe undergoes unzipping in the nanopore.
The instant invention also includes probes, nanopores, kits comprising the probes and nanopores, and associated methods of use described in the following portions of the specification, drawings, and claims provided herewith.



DESCRIPTION OF DRAWINGS
FIG. 1 is a schematic diagram of an exemplary nanopore sensing system, according to one embodiment of the invention.
FIG. 2 is a schematic amplified illustration of an exemplary nanopore plugged through the lipid bilayer.
FIG. 3 includes an exemplary current trace, an amplified electrical mark of the signature event, and a schematic illustration of the unzipping and translocation event.
FIGS. 4(A) to (D) illustrate an exemplary detection of miR-155 with the inventive probe, P155: FIG. 4(A) is a current trace showing different types of blocks. The block profiles and corresponding molecular processes is depicted in panel B, C and D; FIG. 4(b)/(b′) is an exemplary spike-like short block produced by free miR-155 or P155 molecules translocating through the pore; FIG. 4(c)/(c′) is an amplified long block with multiple conductance levels which are sequentially generated by unzipping of a miR-155·P155 hybrid, confinement of miR-155 in the nanocavity and translocation of miR-155; FIG. 4(d)/(d′) is an exemplary long blocks with a single conductance level, produced by a trapped mir-155·P155 hybrid that exits from the cis entrance without unzipping.
FIGS. 5(a) and (b) illustrate the invention employed in the quantification of miRNAs using the nanopore sensor. a. Current traces in the presence of 100 nM P155 and different concentrations of miR-155. The signature events for miR-155·P155 hybrid interacting with the pore are marked with red arrows. b. Correlation between the miR-155 concentration [miR-155] and the frequency of signature events (fsig). Significance (p<0.01) is valid between detections in any two miR-155 concentrations. The fsig-[miR-155] curve is fitted using Eq. 1.
FIGS. 6(a) and (b) illustrate the invention employed in the differentiation of miRNAs with similar sequences. a. Current traces for detections of let-7a and let-7b using the probe Pa (+120 mV). b. Durations of signature events (τsig) for let-7a·Pa, let-7b·Pa, let-7a·Pb and let-7b·Pb.
FIGS. 7(a) to (f) illustrate the invention employed in the detection of miR-155 in lung cancer patients' plasma. a through d, the signature events were found in current traces for total plasma RNAs from normal volunteers (a) and lung cancer patients (b) in the presence of 100 nM P155 probe, but not observed in the absence of P155 (c and d). The traces were recorded in 1 M KCl at +100 mV. e. Frequencies of miR-155 signature events (fsig) from six normal individuals (N1 to N6) and six patients with lung cancer (P1 to P6). Each sample was measured n times (n≧4) with independent nanopores. The data was given as mean±SD. The patient conditions are, P1, metastatic squamous lung carcinoma; P2, recurrent small-cell cancer; P3, early stage of small-cell carcinoma, status post chemotherapy and radiation; P4, early stage of small-cell cancer, status post chemotherapy; P5, late stage non-small cell carcinoma, statue post resection and chemotherapy; P6, late stage adenocarcinoma, status post chemotherapy. f. Relative miR-155 levels in normal group (left) and lung cancer patient group (right), measured with the nanopore sensor and qRT-PCR. The means with SD were shown.
FIG. 8 A-G shows the capture of single miRNA molecules with a specific probe in the nanopore. A. Molecular diagram of the miRNA·probe hybrid; B. Sequentially-occurred nanopore current blocks in the presence of 100 nM miR-155/P155 mixture in the cis solution. The recording solution contained 1 M KCl buffered with 10 mM Tris (pH 8.0). Traces were recorded at +100 mV. The identified current patterns and corresponding molecular mechanisms were depicted in panel c, f and g. The framed blocks demonstrated the multi-level current pattern depicted in panel c and d; c. Multi-level long block at +100 mV, generated by the miR-155·P155 hybrid that was trapped in the pore, unzipped, followed by sequentially translocation of unzipped P155 and miR-155 through the pore; d. Characteristic multi-level long blocks at +150 mV and +180 mV; e. qRT-PCR-detected miR-155 levels in cis and trans solutions after  ̃6 hours electrical recording at different concentrations of miR-155 and P155 presented in the cis solution (Text in Supplementary Information); f. Single-level current pattern generated by a trapped mir-155·P155 hybrid that exited the pore from the cis entrance without unzipping; g. Spike-like short block generated by translocation of un-hybridized miR-155 or P155 from the cis solution.
FIG. 9 A-B shows enhancing detection sensitivity by optimizing the probe sequence. A. Left, current traces showing the frequency of signature events for miR-155 hybridized with the probes P5′-C30 (top), P3′-C30 (middle) and P155 (bottom), monitored at +100 mV in 1 M KCl. Right, the occurrence rate constant of signature events for miR-155 detection with different probes (Table 5). Significance (p<0.005) is valid between results with any two probes; B. Left, [miR-155]-f155 correlation for target concentration ranging between 10-100 nM. Right, [miR-155]-f155 correlation measured in 0.5 M/3 M (cis/trans) KCl asymmetrical solutions for much lower target concentration between 0.1-100 pM. Significance (p<0.01) was valid between detections at any two miR-155 concentrations.
FIG. 10 A-D shows differentiation of let-7 miRNAs containing one or two different nucleotides. The sequence of let-7a, -7b, and -7c were given in Table 3). a. Detections of let-7a and let-7b using the probe Pa or Pb at +120 mV. Left, current traces, and right, comparison of signature event duration (τsig); b. Detections of let-7a and -7c using the probe Pa or Pc at +100 mV. Left, current traces, and right, comparison of signature event duration (τsig). Data was shown in Table 6; c. Receiver Operating Characteristic (ROC) curves for discrimination of events for miRNA·probe hybrids without (positive) and with mismatches (negative). □: let-7a·Pa/let-7b·Pa, ◯: let-7b·Pb/let-7a·Pb, ▪: let-7a·Pa/let-7c·Pa, and ●: let-7c·Pc/let-7a·Pc; d. Correlation between the areas under the ROC curves (AUC) and the duration ratio between fully-match events and mismatch events. ●: AUC measured from the ROC curves in panel c (Table 7), ◯: AUC calculated from ROC analysis based on simulated datasets (FIG. 16 and Table 8). The events were generated with an exponentially distributed duration. The duration ratios of fully-match events (positive) and mismatch events (negative) were 1, 2, 3, 4, 5 and 10 respectively.
FIG. 11 A-H shows detection of miR-155 in lung cancer patients' plasma. a through d. Signature events found in current traces for total plasma RNAs from normal volunteers (a) and lung cancer patients (b) in the presence of 100 nM P155 probe, no signature events observed in the absence of P155 (c and d). The traces were recorded in 1 M KCl at +100 mV. e. Frequencies of miR-155 signature events (f155) from six normal individuals (#1 to #6) and six patients with lung cancer (#7 to #12) in the presence of spiked-in synthetic miR-39. f. Frequencies of miR-39 signature events detected by P39 (see the sequence in Table 3) from all sample used in e. Each sample was measured n times (n≧4) with independent nanopores. The data was given as mean±SD. The patient conditions were, #7, metastatic squamous lung carcinoma; #8, recurrent small-cell cancer; #9, early stage of small-cell carcinoma, status post chemotherapy and radiation; #10, early stage of small-cell cancer, status post chemotherapy; #11, late stage non-small cell carcinoma, statue post resection and chemotherapy; #12, late stage adenocarcinoma, status post chemotherapy. g. f155/f39 calculated from panel e and f. h. Box and Whiskers plot of relative miR-155 level in normal and lung cancer groups, measured with the nanopore sensor and qRT-PCR. The boxes mark the interval between 25th and 75th percentiles. The black lines inside the boxes denote the medians. The whiskers denote the interval between the 5th and 95th percentiles. Filled circles indicate data points outside the 5th and 95th percentiles. Data were given in Table 9.
FIG. 12 A-D shows histograms of block durations. a. Signature blocks generated by the mir-155·P155 hybrid. b. The short Level 3 state in the signature block. c. and d. Short blocks by translocation of miR-155 (c) and P155 (d) alone. Data was obtained from current traces recorded in 1 M KCl at +100 mV.
FIG. 13 shows voltage-dependent frequency of mir-155·P155 signature events. Data was obtained from current traces recorded in 1 M KCl with 10 (Δ) and 25 (□) nM mir-155 in the presence of 100 nM P155, and 10 pM mir-155 in the presence of 5 pM P155 (◯)
FIG. 14 shows the frequency of miR-155 signature events detected using P155 (100 nM) in the presence of other synthetic miRNA components. The three bars represented miR-155 alone (50 nM), miR-155 in the presence of Let-7a (50 nM), and that in the presence of both Let-7a (50 nM) and -7b (50 nM). Data was obtained from current traces recorded in 1 M KCl at +100 mV.
FIG. 15 shows duration histograms of signature events formed by Let-7a·Pa and Let-7b·Pa hybrids. Data was obtained from current traces recorded in 1 M KCl at +120 mV. a. Let-7a·Pa. b. Let-7b·Pa. Concentrations of all RNA and DNA components were 100 nM.
FIG. 16 A-B shows simulation on separation of fully-match (positive) and mismatch (negative) events based on event duration. a. ROC curves at various duration ratios. There were 400 events of both types participating in the analysis; b. ROC curves at various event number ratios of the two type of events. The duration ratio τP/τN=3.
FIG. 17 shows translocation frequencies of miR-155 and P155. Data was obtained from current traces recorded in 1 M KCl at +100 mV. The concentrations of both oligos were 100 nM.
FIG. 18 A-F shows: A) the diagram of the miRNA/probe complex. FIG. 18 B shows events for translocation of the peptide-PNA probe, P7b. The characteristic events last for 3 ms and reduce the current to 10 pA at +180 mV. FIG. 18C shows that no block events can be observed with free miRNA let-7b (without probe) in the solution at +180 mV. FIG. 18D shows signature events for the trapping of the let-7b/P7b complex. FIG. 18E shows that Let-7c, which has two different nucleotides from Let-7b, cannot bind to PNA of the probe P7b, therefore does not generate signature events as in FIG. A1c. Almost all observed events are due to the probe itself. FIG. 18F compares the duration-amplitude property for P7b binding to Let-7b (fully match, two separate clusters without overlay) and Let-7c (2 mismatches, two clusters fully overlay).
FIG. 19 A-C shows: FIG. 19A, when employing HP-C30 with a hairpin at the 3′-end of short strand, we observed a novel type of three-level current pattern. FIG. 19B, when using SA-C30 attached with a streptavidin at the 3′-end of the short strand, we also observed a new multi-level current pattern. FIG. 19C shows when using a short oligonucleotide to link two DNAs, the complex can be sequentially unzipped in the nanopore in two steps. The two unzipping can be clearly revealed by the two Level 2 states.



DETAILED DESCRIPTION OF INVENTION
The invention provides a robust nanopore sensing system that enables sensitive, selective and direct detection, differentiation and quantification of single strand oligonucleotide, such as miRNAs. Additionally, the inventive sensing technology can also be employed to discriminate single nucleotide differences in miRNA family members. Furthermore, the inventive technology has the potential for non-invasive and cost-effective early diagnosis and continuous monitoring of cancer markers in patients' blood samples.
The inventive nanopore sensing system for detecting a target single-strand oligonucleotide, such as a miRNA, includes 1) a nanopore allowing rapid translocation of single-stranded oligonucleotide, 2) a power source providing a pre-determined voltage as driving force to induce unzipping of a double-stranded oligonucleotide, 3) a probe molecule to be mixed with the target oligonucleotide and loaded into the nanopore, and the unzipping of the hybrid of target oligonucleotide and the probe in the pore produces certain identifiable current signal changes, and 4) a means for detecting current changes.
Refer to FIG. 1, which is a schematic illustration of an exemplary nanopore sensing system. As shown in FIG. 1, the sensing chamber, 1, includes a cis compartment, 2, and a trans compartment, 3, which are divided by a partition, 4. Both compartments are filled with a pre-selected recording solution such as 1 M KCl. The partition, 4, has an opening, 5, in its center region, over which a lipid bilayer is formed, and the nanopore, 6, is plugged through the lipid bilayer. The power, 7, provides a voltage that is loaded through a pair of electrodes in the two compartments; the current detector, such as a pico-Ampere amplifier, 8, is connected to monitor the current changes. Upon the testing, a mixture sample of the target oligonucleotide, 9, and its complementary probe, 10, is loaded into the cis compartment, 2.
Refer to FIG. 2, which is a schematic amplified illustration of the nanopore, 6. As shown in FIG. 2, the nanopore, 6, is in conical or funnel shape with two openings, the cis opening, 11, at the wide end and the trans opening, 12, down the narrow end. During the detection, the paired oligonucleotides, 9/10, is captured into the nanocavity, 13. The voltage then drive the oligonucleotides, 9/10, to unzip at the constriction, 14, with the probe, 10, first traversing through the β-barrel, 15, and out off the trans opening, 12, and followed by the traversing of the target oligonucleotide, 9.
The nanopore may be any ion channel of cone-shape or any asymmetrical shape with a wide and a narrow opening plugged into the planar lipid bilayer that has a wider cavity followed by a narrow channel that can facilitate unzipping translocation events. The nanopore may be any existing protein ion channels, such as the α-hemolysin transmembrane protein pore adopted in the examples below, or various synthetic pores fabricated using fashion nanotechnologies with abiotic materials such as silicon.
The inventive probe is a multi-domain single strand molecule, which comprises a central domain fully complementary to the target oligonucleotide and at least one terminal extension, i.e., signal tag, at its 3′ or 5′ terminal, with signal tags at both terminals as preferred. The invention suggests the 3′-tagged probe is preferred over the 5′-tagged probe. The probe directionality-dependence of the capture rate is possibly due to that the bases of ssDNA tilt collectively toward the 5′ end of the strand38, and this asymmetric base orientation makes DNA move more easily from 3′-end than 5′-end.
The terminal extension (signal tag) may be of any charged single chain molecule with sufficient length to assist the unzipping translocation through the nanopore driven by the voltage. The signal tag may be a charged polymer chain, which can be an oligonucleotide such as poly(dC)n, poly(dA)n, and or poly(dT)n, or a charged polypeptide. For example, when α-hemolysin transmembrane protein pore is employed as the nanopore, the poly(dC) tag is more preferred over poly(dA) or poly(dT) tags; furthermore, the poly(dC)30 is much more efficient in generating signature events (discussed below) than that with a shorter tag such as poly(dC)8. The capture rate can be further enhanced once combined with other effective approaches, including detection at high voltage, use of engineered pores with designed charge profile in the lumen33, and detection in asymmetrical salt concentrations between both sides of the pore39.
The invention also provides a method of detecting and differentiating single strand oligonucleotides by monitoring the current changes induced by the unzipping and translocation of the oligonucleotides through a nanopore. The inventive method of detecting single strand oligonucleotide with a dual-compartment nanopore system, as the one illustrated in FIG. 1, includes the steps of 1) mixing the target oligonucleotide with a pre-designed probe, which has its central domain matching the target sequence and a charged single chain molecule tagged to at least one of its 3′ and 5′ terminals, to produce a sample mixture, 2) loading the mixture into the cis compartment, 3) providing the system with a pre-determined voltage, and 4) recording current output for a pre-determined time period. The current change induced by the unzipping and translocation of the hybrid of the target oligonucleotide and its complementary probe through the nanopore is a unique signature event, which is used to detect and differentiate the target oligonucleotide. Refer to FIG. 3, which includes an exemplary current trace recorded during an exemplary detection, an amplified electrical mark of the signature event, and a schematic illustration of the unzipping-translocation event.
Further Description
DEFINITIONS
As used herein, the term “ROC curve” refers to a Receiver Operating Characteristic Curve. An ROC curve used to analyze the relationship between selectivity and sensitivity. An ROC curve separates the plot into up and lower regions.
As used herein, the term “AUC” refers to the Area under the ROC curve. An AUC can range between 0.5-1.0. The higher the AUC value, the better the separation result.
As used herein, the term “OCP” refers to an Optimized Cutoff Point. In certain embodiments, an OCP can be calculated from ROC curves. In certain embodiments, an OCP is a cutoff duration at the maximal value of a Youden index.
As used herein, the phrase “Youden index” is defined as {sensitivity+selectivity−1}. A Youden index is calculated from the ROC curve, and can range between 0 and 1. A cutoff duration leading to complete separation of long and short duration distribution results in Youden index=1, whereas complete overlap gives Youden index=0. In certain embodiments, a cutoff duration value that returns the maximum of Youden index, i.e. “optimal” cutoff point (OCP) (Greiner et al., 2000 Preventive Veterinary Medicine 45, 23-41) gives the most accurate separation.
Description of the Probes, Nanopores, Kits Comprising the Probes and Nanopores, and Associated Methods of Use
In one broad aspect, the instant invention is directed to probes, nanopores, kits comprising the probes and nanopores, and associated methods of use, that provide for “signature” current blockage events that distinguish those events arising from interactions with the probe and target from other events. In this context, the other events are referred to as “background” events. Background events include, but are not limited to, interactions of a probe with nucleic acid that is not a target, interactions of a probe with other components present in a nanopore detection system, free nucleic acids present in the nanopore detection system, and the like. Such features of such signature events include, but are not limited to, at least one of a: i) a current block of different duration than a background current block; ii) a different number of distinct current blockade levels than a background current block; iii) a different order of occurrence of current blockade levels than a background current block; iv) a different current amplitude at a blockade level than a background current block; v) a different current amplitude of each blockade level than a background current block; or any combination of (i), (ii), (iii), (iv), or (v). In certain embodiments, a signature blockage event can be distinguished from a background blockage event by differences in a characteristic background noise of each blockage event. In certain embodiments, the distinct durations, numbers, or amplitude(s) in the signature event are greater than those observed in the background event. In certain embodiments, the distinct durations, numbers, or amplitude(s) in the signature event are less than those observed in the background event. In certain embodiments, the distinct durations, numbers, orders, or amplitude(s) in a signature event are statistically distinguishable from those of a background event. In certain embodiments, the signature events are provided in nanopore systems comprising a protein nanopore formed by alpha-hemolysin (αHL) or engineered variants thereof in a planar lipid bilayer system. In certain embodiments, the signature events can be provided in a biochip formed by hydrogel-encapsulated lipid bilayer with a single protein nanopore embedded therein or a micro-droplet bilayer system. Biochips and micro-droplet bilayer systems have been described (Shim and Gu; Stochastic Sensing on a Modular Chip Containing a Single-Ion Channel Anal. Chem. 2007, 79, 2207-2213; Bayley, H. et al. Droplet interface bilayers. Mol. Biosyst. 4, 1191-1208 (2008).
In certain embodiments, the signature events can be provided in a synthetic nanopore. Synthetic nanopores include, but are not limited to, nanopores comprising silicon nitride or graphene.
Probe molecules provided herein comprise terminal extensions at one or both of their 5′ and/or 3′ termini. Without seeking to be limited by theory, it is believed that these terminal extensions provide useful functions that include, but are not limited to, trapping of the probe/target complex into the nanopore at a high rate (i.e. the number of signature events per unit target concentration per unit recording time). The trapping rate directly determines the sensitivity. In the same target concentration and the same recording time, a higher trapping rate gives a more precise sensing result. Without seeking to be limited by theory, it is also believed that these terminal extensions provide useful functions that include, but are not limited to, inducing the voltage-driven dissociation of the probe/target complex. This dissociation function generates a signature event that can be used to discriminate interactions of the probe with the target from other components in the mixture, thereby ensuring the selectivity or specificity.
Probe terminal extensions can comprise a charged polymer of any length. In certain embodiments, the polymer can be a negatively charged single-stranded nucleic acid. Advantages of such nucleic acid terminal extensions include, but are not limited to, extremely low cost of synthesis and controllable charge by pH, salt concentration and temperature. Such nucleic acid extensions can comprise homopolymers, heteropolymers, copolymers or combinations thereof. In certain embodiments, the lengths of such nucleic acid terminal extensions can range from about 1 or 2 nucleotides to about 50 nucleotides. In still other embodiments, the nucleic acid extensions can range in length from about 5 to about 40 nucleotides, about 15 to about 35 nucleotides, or from about 20 to about 35 nucleotides. An exemplary terminal extension provided herewith is homopolymer poly(dC)30. However, a heteropolymeric sequence, including but not limited to, di- or tri-nucleotide heteropolymers such as CTCTCTCT . . . , or CATCATCAT . . . , can also be used. In certain embodiments, co-polymers comprising abases or polyethylene glycol (PEG) can be used in the terminal extension. These co-polymers, or domains thereof in a terminal extension, can confer new functions on the terminal extension of the probe. An abase is a nucleotide without the base, but carries a negative charge provided by the phosphate. As the dimension of abase is narrower than normal nucleotides, it may generate a signature event signal different from that formed by the neighbor nucleotides. PEG is not charged. Without seeking to be limited by theory, it is believed that when the PEG domain in a nucleic acid sequence is trapped in the pore, it can reduce the driving force, thus precisely regulating the dissociation of the probe/target complex.
Probe terminal extensions can also comprise a polypeptide. The richer choice of amino acids makes the sequence and functionality of the polypeptide terminal extension more programmable than an oligonucleotide terminal extension. For example, polypeptide terminal extensions allow insertion of charged amino acids in the optimized positions to generate more distinguishable probe/target signature events. While not seeking to be limited by theory, it is believed that the probe/target complex can be selectively trapped using a probe comprising a positively charged polypeptide terminal extension under an appropriate voltage while all other negatively charged non-target oligonucleotides in the mixture are prevented from entering into the pore, resulting in ultra-selective detection. In certain embodiments, the polypeptide terminal extensions can comprise two, three, four, or more amino acid residues that can carry a positive charge (i.e. lysine and/or arginine and/or histidine). In certain embodiments, sufficient numbers of positively charged residues are included in the polypeptide terminal extension to provide a net positive charge when said probe is hybridized to a target oligonucleotide. In certain embodiments where probes comprising terminal extensions with positive charges conferred by residues such as lysine, arginine or histidine, performance of the associated nanopore based detection methods can be enhanced under acidic conditions (i.e. when the pH value is less than 7) or conditions where the residue will be protonated. Thus, the use of such probes at pH values of about 1 to about 6.9, 1 to about 6.0, about 1 to about 5.5, about 3 to about 5.5, and the like. In certain embodiments, the lengths of such polypeptide terminal extensions can range from about 1 or 2 residues to about 30 residues. In still other embodiments, the polypeptide extensions can range in length from about 5 to about 20 residues, about 8 to about 20 residues, or from about 8 to about 15 residues. In an exemplary embodiment, an HIV-TAT polypeptide comprising positively charged arginine and lysine residues can be used as the terminal extension. In certain embodiments, the center domain of the probe that is complementary to the target oligonucleotide can comprise a peptide nucleic acid that is covalently linked to a terminal extension comprising amino acids that carry a positive charge. In certain embodiments, a center domain comprising a peptide nucleic acid is used in conjunction with a terminal extension comprising amino acids that carry a positive charge to provide a net positive charge when said probe is hybridized to a target oligonucleotide. In certain embodiments, polypeptide terminal extensions comprising amino acids with aromatic side chains including, but not limited to, phenylalanine, tryptophan, tyrosine, thyroxine, and the like, can be incorporated into the polypeptide terminal extensions. While not seeking to be limited by theory, it is believed that such aromatic amino acids can interact with the pore through aromatic stacking and provide for useful changes in the signature obtained in nanopore based detection methods.
Without seeking to be limited by theory, it is believed that if there are a sufficient number of positively-charged amino acids in the polypeptide terminal extension such that the net charges of the target oligonucleotide/probe complex are still positive when the probe comprising the terminal extension is hybridized to the target oligonucleotide, then the entire target oligonucleotide/probe complex will form a strong dipole molecule. It is believed that the positively-charged peptide domain of the probe dipole will be both pushed by the positive voltage (cis grounded) and attracted by the negative ring at the trans opening, guiding the trapping of the oligonucleotide/probe complex into the pore. At the positive voltage, any other free nucleic acids components will be repulsed from entering the pore due to the negative charge that is carried by the free, unhybridized nucleic acids. This significantly reduces signals by free nucleic acid components, such that the majority of the observed current blockage events are either due to the trapping of the oligonucleotide/probe complex or to the translocation of the probe.
In certain embodiments, the oligonucleotide/probe complexes with a net positive charge can be directed to a nanopore with a negatively-charged ring at the trans-opening of the pore. In this context, a trans opening of a pore is understood to be that portion of the pore from which a molecule would emerge whereas a cis opening of a pore from which a molecule would enter. In these embodiments, it is understood that a negative charged ring at the trans-opening of the pore can be obtained by using any type of nanopore that has been suitably synthesize and/or derivatized so as to have a negative charged ring at the trans-opening of the pore. Such nanopores with a negatively charged ring at the trans opening of the pore include, but are not limited to, protein nanopores and synthetic nanopores. Protein nanopores with a negatively charged ring at the trans opening of the pore include, but are not limited to, engineered variants of an alpha-hemolysin protein. In certain embodiments, the engineered alpha hemolysin variant can comprise a Staphylococcus aureus alpha hemolysin containing a K131D, a K131E, or a K131H amino acid substitution. Exemplary and non-limiting Staphylococcus aureus alpha hemolysin wild type sequences are provided herein (SEQ ID NO:20, nucleic acid coding region; SEQ ID NO:21: protein coding region) and available elsewhere (National Center for Bioinformatics or GenBank Accession Numbers M90536 and AAA26598). An exemplary and non-limiting Staphylococcus aureus alpha hemolysin variant comprising a K131D substitution is provided as SEQ ID NO:22. In certain embodiments, the engineered alpha hemolysin variant can comprise a suitably derivatized variant that is derivatized with moieties that provide for a negatively charged ring at the trans opening of the pore. An exemplary wild type S. aureus alpha hemolysin protein that can be substituted or derivatized to provide for a protein nanopore with a negative charged ring at the trans-opening of the pore is provided herewith as SEQ ID NO: 21. However, variants of other hemolysins capable of forming pores can be substituted or derivatized to provide for a protein nanopore with a negative charged ring at the trans-opening of the pore. Synthetic nanopores with a negatively charged ring at the trans opening of the pore are also provided. In certain embodiments, such synthetic nanopores with a negatively charged ring at the trans opening of the pore include, but are not limited to, silicon nitride or graphene nanopores that have been suitably derivatized with moieties that provide for a negatively charged ring at the trans opening of the pore.
The center domain of probes provided herein is used to capture the target molecule. In certain embodiments, the center domain can be fully complementary or partially complementary to the target sequence. In certain embodiments, a center domain can comprise an oligonucleotide comprising natural nucleotides (A, T, G, C (DNA) or a, u, g, c (RNA)), and/or artificial nucleotides including, but not limited to, nucleosides such as inosine, xanthosine, 7-methylguanosine, Dihydrouridine, and 5-methylcytidine. In certain embodiments, the center domain can comprise a locked nucleic acid (LNA) or a peptide nucleic acid (PNA). Locked nucleic acids comprise RNA derivatives where the ribose ring contains a methylene linkage between the 2′-oxygen and the 4′-carbon. Peptide nucleic acids (PNA) comprise a peptide backbone with nucleobase side chains. In certain embodiments, a LNA or a PNA center domain can comprise natural nucleobases (adenine, guanine, thymine, cytosine or uracil) and/or artificial nucleobases including, but not limited to, hypoxanthine, xanthosine, 7-methylguanine, 5,6-dihydrouracil, and 5-methyl cytosine. In certain embodiments, probe center domains comprising co-polymers of oligonucleotides, LNA, or PNA are provided. In certain embodiments, a center domain of a probe will have at least about 4, 6, 8, 10, 12, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, or 25 nucleotide or nucleobase residues that are complementary to the target nucleic acid. In certain embodiments, a central region of a probe will have at least about 4, 6, 8, 10, 12, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, or 25 to any of about 30, 35, 40, or 50 nucleotide or nucleobase residues that are complementary to the target nucleic acid. In certain embodiments, synthetic nucleotides or nucleobases inserted in the sequence can precisely adjust the hybridization energy with the target, such that one can distinguish the characters of targets such as single-nucleotide polymorphism, methylation, or interaction between miRNA and its target messenger RNA.
A variety or target nucleic acids or oligonucleotides that can be detected and distinguished from non-target nucleic acids by the probes, nanopores, kits comprising the probes and nanopores, and associated methods of use probes, provided herein. In certain embodiments, the target can be a nucleic acid or a fragment thereof from cells, body fluid, tissues, bacteria, or a virus. In certain embodiments, the target can be a PCR products or a synthetic oligonucleotide. In certain embodiments, a target can comprise a genomic DNA, an mRNA, a pre-mature or mature miRNA, an artificial miRNA, non-coding DNA or RNA, a nucleic acid biomarker, or a synthetic aptamer. In certain embodiments, a miRNA targets may come from the RNA extraction from bio-fluid from any tissues such as plasma and formalin-fixed and paraffin-embedded tissues. In certain embodiments, a target nucleic acid can comprise be a nucleic acid fragment complexed with any of a nucleic acid binding protein, an antibody, or an aptamer bound with a target protein. In certain embodiments, a target nucleic acid can comprise be a nucleic acid fragment complexed with a low molecule weight compound, including, but not limited to, a drug. In certain embodiments, targets can include sequences with mutations, with single-nucleotide polymorphisms, or with chemical modifications such as methylation and phosphorylation.
EXAMPLES
Example 1
Detection of miR-155, a Lung Cancer Biomarker
The invention further provides an exemplary nanopore sensing system for detection of miR-155, a lung cancer biomarker. The nanopore sensing system includes an α-hemolysin transmembrane protein pore and a pre-designed probe for miR-155. The probe is a DNA multiple-block copolymer with its central domain complementary to the target miR-155, and at least one poly(dC)30 extension at 3′-, 5′-, or both terminals functioning as signal tags. Table 1 lists the sequences of miR-155 and the exemplary probes with the tri-block copolymer, P155 as preferred.







TABLE 1







Sequences of miRNAs and their probes









miRNA
Probe
Sequence





mir-155

5′-UUAAUGCUAAUCGUGAUAGGGG-3′




(SEQ ID NO: 1)


 



Pnt
5′-CCCCTATCACGATTAGCATTAA-3′




(SEQ ID NO: 2)


 



P5′-C30
5′-C30-CCCCTATCACGATTAGCATTAA-3′




(SEQ ID NO: 3)


 



P3′-C30
5′-CCCCTATCACGATTAGCATTAA-C30-3′




(SEQ ID NO: 4)


 



P155
5′-C30-CCCCTATCACGATTAGCATTAA-C30-3′




(SEQ ID NO: 5)









During an exemplary sensing process, a mixture of miR-155/P155 is added to the cis side of the pore, a current trace with a series of short- and long-lived current blocks can be recorded, as shown in FIG. 4a, while being monitored in 1 M KCl at +100 mV. The spike-like short blocks, labeled as b in FIG. 4a and also shown in FIG. 4b, have duration of 220±21 μs and almost fully reduce the pore conductance with a relative conductance g/g0=0.16, where g and g0 are conductance of blocks and unoccupied nanopore. Both the duration and conductance are similar to that of blocks by miR-155 and P155 alone, thus the short blocks in the mixture are associated with the rapid passage of single-stranded free miR-155 or P155 through the pore, as illustrated in FIG. 4b′. 
In contrast to short blocks, the long blocks, labeled as c and d in FIG. 4a, in the recording persist for 250±58 ms. One type of long blocks, labeled as c in FIG. 4a, features three sequential conductance levels, Level 1→Level 2→Level 3 (expended in FIG. 4c). This type of blocks is not observed when either miR-155 or P155 alone is presented, indicating that it is originated from the miR-155/P155 hybrid (miR-155·P155). Level 1 almost fully reduces the conductance to g/g0=0.15. This conductance level is consistent with a configuration that miR-155·P155 is trapped in the pore from the wider opening (cis), with either 3′- or 5′-signal tag of P155 occupying the narrowest β-barrel (FIG. 4c′ level 1). The signal tag in the β-barrel can induce unzipping of miR-155·P155, driven by voltage. The unzipping time, or the duration of Level 1, is comparable to previously reported time scales for DNA unzipping in the pore, e.g.  ̃435 ms for unzipping a 50 bps dsDNA at +140 mV, and  ̃40 ms for a 10 bps hairpin DNA at +90 mV. The unzipping process was further evidenced by the discrete transition from Level 1 to Level 2. Level 2 lasts 410±20 μs and its conductance significantly increases to g/g0=0.42 (FIG. 1c′ level 2). This partial block can not be interpreted as an oligonucleotide occupying in the β-barrel. It is very likely that, after unzipping of miR-155·P155 followed by translocation of P155, mir-155 can be temporarily confined in the nanocavity of the pore. It has been verified that a single-stranded oligonucleotide residing in the nanocavity can generate such a partial block33,34. The miR-155 molecule in the nanocavity finally traversed the β-barrel, generating Level 3 which fully reduces the conductance to g/g0=0.08 (FIG. 4c). The duration of Level 3 is 270±30 μs, close to the 220 μs for short blocks by mir-155 alone, and consistent with the time scale of  ̃400 μs for translocation of a 75 bases RNA at +120 mV,35 and 800 μs for a 210 bases RNA at +120 mV.36 The duration of Level 3 becomes shortened as the voltage increases, further supporting the translocation of a single-stranded oligonucleotide for this conductance level.
In addition to the multi-conductance long blocks, the long blocks with single conductance level at g/g0=0.15 (labeled as d in FIG. 4a and expended in FIG. 4d) have also been observed. This type of long blocks may occur when the arrested miR-155·P155 exits the pore from the cis entry without unzipping.
The invention teaches that the characteristic long blocks can serve as signature events for identifying single molecules of target miRNAs. From the frequency of signature events (fsig), the target miRNA can be quantified using Eq. 1 (Methods in Supplementary Materials, provided below in Example 2),




  
    
      
        
          f
          sig
        
        =
        
          
            k
            on
          
          ⁢
          
            
              
                
                  
                    
                      (
                      
                        
                          
                            [
                            miR
                            ]
                          
                          0
                        
                        +
                        
                          
                            [
                            P
                            ]
                          
                          0
                        
                        +
                        
                          K
                          d
                        
                      
                      )
                    
                    -
                  
                
              
              
                
                  
                    
                      
                        
                          (
                          
                            
                              
                                [
                                miR
                                ]
                              
                              0
                            
                            +
                            
                              
                                [
                                P
                                ]
                              
                              0
                            
                            +
                            
                              K
                              d
                            
                          
                          )
                        
                        2
                      
                      -
                      
                        
                          
                            
                              4
                              ⁡
                              
                                [
                                miR
                                ]
                              
                            
                            0
                          
                          ⁡
                          
                            [
                            P
                            ]
                          
                        
                        0
                      
                    
                  
                
              
            
            2
          
        
      
    
    
      
        (
        1
        )
      
    
  




In Eq. 1, [miR]0 and [P]0 are the initial concentrations of miRNA and the probe, kon is the occurrence rate of signature events and Kd is the dissociation constant for miR·P in the solution. When [miR]0<<[P]0, fsig≈kon[miR]0. The current traces indeed show more frequent miR-155·P155 signature events as the miR-155 concentration increases (FIG. 5a). The fsig-[miR-155] data can be best fitted using Eq. 1, with Kd=30 nM and kon=3.6×106 M−1s−1 (FIG. 5b). According to the literature8, the mean concentrations of circulating miRNAs were 158.6 ng/mL ( ̃25 nM) for the lung cancer group versus 68.1 ng/mL ( ̃10 nM) for the control group. Therefore we compared the fsig values at 10 nM and 25 nM mir-155 (FIG. 5b). Analysis indicates that the two levels of miRNA concentration can be separated (p<0.005), suggesting that the inventive method has the potential to differentially detect miRNA levels in lung cancer patients.

The invention also provides an exemplary process employing the inventive nanopore sensing system to differentiate highly similar miRNA sequences, let-7a and let-7b. let-7a and let-7b are members of the Let-7 tumor suppressing miRNA family4-6; and the two Let-7 members only contain different nucleotides at the position 17 and 19, which are adenines in let-7a and guanines in let-7b. The inventive probes Pa and Pb are designed for let-7a and let-7b respectively with sequences listed in Table 2.







TABLE 2







Sequences of let-7a and let-7b and their probes









miRNA
Probe
Sequence





Let-7a

5′-UGAGGUAGUAGGUUGUAUAGUU-3′




(SEQ ID NO: 6)


 



Pa
5′-C30-AACTATACAACCTACTACCTCA-C30-3′




(SEQ ID NO: 7)


 


Let-7b

5′-UGAGGUAGUAGGUUGUGUGGUU-3′




(SEQ ID NO: 8)


 



Pb
5′-C30-AACCACACAACCTACTACCTCA-C30-3′




(SEQ ID NO: 9)









As shown in FIG. 6a and FIG. 6b, when using Pa to detect each miRNA, the duration of signature events (τsig) for let-7a·Pa (without mismatch) is 155±28 ms, whereas τsig for let-7b·Pa (with 2-nt mismatches) is significantly shortened to 48±11 ms (p<0.005). Similarly, when using Pb to detect both miRNAs, τsig for let-7b·Pb (without mismatch) is 165±47 ms, significantly longer than the 24±2 ms for let-7a·Pb (with 2-nt mismatches) (p<0.005) as shown in FIG. 6b. The significant differences in durations can be interpreted that the mismatches significantly weaken the hybridization interaction between miRNA and the probe. When placed in the identical electrical field, the hybrid containing mismatches needs lower energy than fully matched hybrid to unzip. Thus, the inventive nanopore sensing system is able to differentiate single mismatches based on the unzipping time, thus demonstrating the potential to detect miRNAs with similar sequences and SNPs.
The invention further provides an exemplary process of detecting plasma miR-155 in lung cancer patients with the inventive nanopore sensing system. During the exemplary process, the peripheral blood samples were obtained from six lung cancer patients and six normal volunteers with a local IRB approval. Total plasma RNAs containing miRNAs were extracted from 350 μl of each plasma sample using miRVana PARIS Kit (Ambion), with a final elution volume of 100 μl, which were than divided into two aliquots (50 μl) for the nanopore and RT-PCR assay44. One aliquot was pre-mixed with P155 and directly added to the 2-ml recording solution in the nanopore chamber. The nanopore current retain a low level of noise even in the presence of plasma samples, and distinct short and long blocks (marked with red arrows) can be indentified in both the control group (FIG. 7a) and lung cancer group (FIG. 7b). The characteristic long blocks, including both with multiple conductance and single conductance, features the same conductance profiles and similar properties to that for synthetic miR-155 RNA in FIG. 1a. In the absence of P155, no such types of long blocks can be observed (FIGS. 7c and d), but short blocks were found for translocation of single-stranded oligonucleotides such as free miRNAs (FIGS. 7c and d). Overall, the characteristic long blocks could be attributed to miR-155·P155 hybrids and serve as signature events for single miRNA molecules detection.
The frequency of miR-155 signature events fsig for all samples in the lung cancer patient group varies between 1.15-1.51 min−1, with a mean of 1.40±0.16 min−1 (FIG. 7e). This level was significantly higher than fsig in the control group that ranges between 0.32-0.70 min−1 with a mean of 0.48±0.14 min−1 (FIG. 7e). Since all samples were prepared following a standard procedure (Methods in Supplementary Materials), it should be valid to compare relative miRNA levels in two groups. When the mean fsig value in normal plasma was set as 1, the folds of miR-155 in lung cancer plasma were compared with the two methods. FIG. 7f showed that the relative mir-155 level in lung cancer patients was 2.79 with the nanopore sensor (p<0.001). By comparison, the relative miR-155 level was 4.72 with RT-PCR method (p<0.02) with greater variability. Therefore, both nanopore and RT-PCR assay indicated a significant elevation of miR-155 in lung cancer patient plasma although there is a 1.69 fold difference. As the nanopore method does not require labeling and amplification, this may be a reason for smaller variability in the nanopore assay (FIG. 7f). Overall the nanopore sensor with engineered probes demonstrates the ability to detect circulating miRNAs in clinical lung cancer patients, which is verified by the independent RT-PCR method.
Example 2
Supplementary Information
Materials. Oligonucleotides including miRNAs and DNA probes were synthesized and electrophoresis-purified by Integrated DNA Technologies (Coralville, Iowa). Before testing, the mixtures of miRNA and DNA probe were heated to 90° C. for 5 minutes, then gradually cooled down to room temperature and stored at 4° C. The RNase-free water was used to prepare RNA solution.
Setup and method of nanopore detection. This section has been well-documented earlier (Shim, J. W., Tan, Q., & Gu, L. Q. Single-molecule detection of folding and unfolding of a single G-quadruplex aptamer in a nanopore nanocavity. Nucleic Acids Res. 37, 972-982 (2009)). Briefly, the recording apparatus was composed of two chambers (cis and trans) that were partitioned with a Teflon film. The planar lipid bilayer of 1,2-diphytanoyl-sn-glycerophosphatidylcholine (Avanti Polar Lipids) was formed spanning a 100-150 nm hole in the center of the partition. Both cis and trans chambers were filled with symmetrical 1 M salt solutions (KCl) buffered with 10 mM Tris and titrated to pH 8.0. All solutions are filtered before use. Single α-hemolysin proteins were inserted into the bilayer from the cis side to form molecular pores in the membrane. All the oligonucleotides including miRNAs and DNA probes and clinical RNA samples were also added to the cis solution. To record the pore current, the cis solution was grounded and the voltage was given from the trans solution. In this convention, a positive voltage can drive the translocation of a negatively charged DNA through the pore from cis to trans. Single-channel currents were recorded with an Axopatch 200A amplifier (Molecular Device Inc. Sunnyvale, Calif.), filtered with a built-in 4-pole low-pass Bessel Filter at 5 kHz, and acquired with Clampex 9.0 software (Molecular Device Inc.) through a Digidata 1332 A/D converter (Molecular Device Inc.) at a sampling rate of 20 kHz. The data were analyzed using Clampfit 9.0 (Molecular Device Inc.), Excel (MicroSoft) and SigmaPlot (SPSS) software.
The translocation of free miRNA or probe through the pore generated very short current block ( ̃101-102 μs). Some short blocks showed partially reduced pore conductance, which may be due to the filtering of the recording at 5 kHz, or formed by the trapped oligonucleotide returning back to the cis solution (Maglia, G., Restrepo, M. R., Mikhailova, E., & Bayley, H. Enhanced translocation of single DNA molecules through +/− hemolysin nanopores by manipulation of internal charge. Proc. Natl. Acad. Sci. U.S.A. 105, 19720-19725 (2008). In our experiments, all short blocks including these partial blocks were collected for histogram construction. Since the translocation events ( ̃101-102 μs) are well distinguished from the signature events ( ̃101-103 ms), we simply used 1 ms as the boundary for separation of short and long events. Data were given as the mean±SD, based on at least three separate experiments (n>3). In the t-test, p<0.05 was considered as a significant difference between two groups. The electrophysiology experiments were conducted at 22±2° C.
Total RNA Extraction From Plasma and miRNA Quantification by qRT-PCR.
Peripheral blood samples were obtained at the University of Missouri Ellis Fischel Cancer Center with an IRB approval. Whole blood with EDTA preservative was centrifuged at 1,600 g for 10 min at room temperature and the plasma was transferred to new tubes. Total RNAs containing miRNAs was extracted from 350 μl of plasma using miRVana PARIS Kit (Ambion, Austin, Tex., USA) according to the manufacturer's protocol. The final elution volume was 100 μl.
A SYBR green-based quantitative RT-PCR assay was employed for miRNA quantification. In brief, 10 μl of total RNA sample containing miRNAs was polyadenylated by poly(A) polymerase (Ambion) and reverse transcribed to cDNA using SuperScript III Reverse Transcriptase (Invitrogen) according to the manufacturer's instructions with a poly(T) adapter primer (5′-GCGAGCACAGAATTAATACGACTCACTATAGGTTTTTTTTTTTTTTTVN-3′; SEQ ID NO: 10)). Real-time PCR was performed using iQ SYBR Green Supermix (Bio-Rad, Hercules, Calif., USA) with the miR-155 specific forward primer (5′-TTAATGCTAATCGTGATAGGGGT-3′; SEQ ID NO:11) and the sequence complementary to the poly(T) adapter as the reverse primer (5′-GCGAGCACAGAATTAATACGAC-3′; SEQ ID NO:12) in iQ5 Real-time PCR system (Bio-Rad, USA). The PCR was carried out as follows: after initial denaturation at 95° C. for 3 min, 40 cycles of 95° C. for 15 s and 60° C. for 1 min were followed. The relative level of miR-155 was calculated using 2−delta Ct method where the level of normal plasma was normalized as 1. Data was presented as mean±SD of three independent experiments, and the differences were considered statistically significant at p<0.05 by using the Student's t-test.
Normalization of the nanopore and qRT-PCR data using spiked-in C. elegans miRNA miR-39 as control. We introduced spiked-in synthetic miRNA as control, to convincingly validate the nanopore sensor's capability of miRNAs detection in human samples. The spiked-in RNA oligonucleotide in the detection matches the sequence of C. elegans miR-39, a miRNA that is absent in the human genomes. 3.5 μL, of 1 nM synthetic miR-39 solution was introduced to each 350 μL plasma sample after addition of the 2× Denaturing Solution (miRVana PARIS Kit) to the plasma, thus the miR-39 concentration in plasma was 10 pM. The Denaturing Solution prevents RNAs from undergoing degradation by inhibiting endogenous plasma RNAases. For each sample, both miR-155 and spiked-in miR-39 were measured using the nanopore sensor and SYBR green-based qRT-PCR. The nanopore data and normalization result were shown in Table 9. In the nanopore detection, the probes for miR-155 and miR-39 were P155 and P39. We first measured the signature event frequencies, f155 and f39, of the hybrids miR-155·P155 and miR-39·P39 respectively. The variability of f39 reflected the difference in miR-39 concentrations among samples after RNA extraction. Therefore the ratio of the two frequencies, f155/f39, should principally eliminate this variability. Finally, we used the Mean f155/f39 of six normal samples (Anormal) as the standard, and calculated each sample's relative miR-155 level by normalizing f155/f39 to Anormal, i.e. f55/f39/Anormal.
Correlation between miRNA concentration and frequency of signature events. In the deduction, “miR” represents miRNA; “P”, probe; Kd, equilibrium dissociation constant for miR·P; kon, occurrence rate constant of miR·P signature events; and fsig, frequency of signature events. In the mixture of miRNA and probe, the equilibrium can be established the reactors miR and P and the product miR·P,




  
    
      
        
          miR
          +
          P
        
        ⁢
        
          
            →
            
              K
              d
            
          
          ←
        
        ⁢
        
          miR
          ·
          P
        
      
    
    
      
        (
        
          Scheme
          ⁢
          
            
          
          ⁢
          1
        
        )
      
    
  




Kd is determined by





  
    
      
        
          K
          d
        
        =
        
          
            
              (
              
                
                  
                    [
                    miR
                    ]
                  
                  0
                
                -
                
                  [
                  
                    miR
                    ·
                    P
                  
                  ]
                
              
              )
            
            ⁢
            
              (
              
                
                  
                    [
                    P
                    ]
                  
                  0
                
                -
                
                  [
                  
                    miR
                    ·
                    P
                  
                  ]
                
              
              )
            
          
          
            [
            
              miR
              ·
              P
            
            ]
          
        
      
    
    
      
        (
        S1
        )
      
    
  




where [miR]0 and [P]0 are total concentrations of miRNA and the probe, and [miR·P] is the concentration of miR·P. Thus the relationship of [miR·] and [miR]0 is





  
    
      
        
          [
          
            miR
            ·
            P
          
          ]
        
        =
        
          
            
              
                
                  
                    (
                    
                      
                        
                          [
                          miR
                          ]
                        
                        0
                      
                      +
                      
                        
                          [
                          P
                          ]
                        
                        0
                      
                      +
                      
                        K
                        d
                      
                    
                    )
                  
                  -
                
              
            
            
              
                
                  
                    
                      
                        (
                        
                          
                            
                              [
                              miR
                              ]
                            
                            0
                          
                          +
                          
                            
                              [
                              P
                              ]
                            
                            0
                          
                          +
                          
                            K
                            d
                          
                        
                        )
                      
                      2
                    
                    -
                    
                      
                        4
                        ⁡
                        
                          [
                          miR
                          ]
                        
                      
                      0
                    
                    +
                    
                      
                        [
                        P
                        ]
                      
                      0
                    
                  
                
              
            
          
          2
        
      
    
    
      
        (
        S2
        )
      
    
  




The kinetics for trapping and unzipping of miR·P in the nanopore is





Because fsig is linearly related to [miR·P],

fsig=kon[miR·P]  (S3)

from Eq. S2 and Eq. S3,





  
    
      
        
          f
          sig
        
        =
        
          
            k
            on
          
          ⁢
          
            
              
                
                  
                    
                      (
                      
                        
                          
                            [
                            miR
                            ]
                          
                          0
                        
                        +
                        
                          
                            [
                            P
                            ]
                          
                          0
                        
                        +
                        
                          K
                          d
                        
                      
                      )
                    
                    +
                  
                
              
              
                
                  
                    
                      
                        
                          (
                          
                            
                              
                                [
                                miR
                                ]
                              
                              0
                            
                            +
                            
                              
                                [
                                P
                                ]
                              
                              0
                            
                            +
                            
                              K
                              d
                            
                          
                          )
                        
                        2
                      
                      -
                      
                        
                          
                            
                              4
                              ⁡
                              
                                [
                                miR
                                ]
                              
                            
                            0
                          
                          ⁡
                          
                            [
                            P
                            ]
                          
                        
                        0
                      
                    
                  
                
              
            
            2
          
        
      
    
    
      
        (
        S4
        )
      
    
  




Eq. S4 suggested that fsig is not in exact proportion to [miR]0, the total concentration of the target miRNA. However, when [miR]0 is considerably smaller than [P]0, which is the case in our miRNA detection, Eq. S4 can be simplified as

fsig≈kon[miR·P]  (S5)

In this condition, fsig is proportional to [miR]0. Eq. S4 also suggested that fsig will ultimately become saturated. This is because fsig measures the capture frequency of miR·P, and the maximal concentration of miR·P ([miR·P]) can not be higher than that of the probe ([P]0).

Identification of miR-155 in trans solutions using RT-PCR. As shown in the model (FIG. 8), with the miRNA·probe complex unzipped in the pore, the separated probe and miRNA can sequentially translocate through the β-barrel to the trans solution. To verify this model, we employed RT-PCR to detect the unzipped miRNAs in the trans solution. However, the PCR method cannot discriminate if the trans miRNAs are from the unzipped miRNAs or from the free miRNA (un-hybridized) that simply translocated from the cis solution to the trans solution. We therefore added a much higher concentration of the probe than miRNA in the cis solution, so that most of the miRNAs molecules are bound with the probe and there is little free miRNA left, eliminating the translocation of free miRNA to the trans solution that can interfere with the PCR result. Our target was miR-155 and the probe was P155. The target/probe concentrations were 0.1/1000, 1/1000 and 10/1000 (nM). After over 6 hours bilayer recording for many pores, 2 μL of both cis and trans solutions were subjected to polyadenylation, reverse transcription and RT-PCR to detect the concentration of miR-155 as indicated in the method described above. Meanwhile, a series of dilution of synthesized miR-155 were performed to construct standard curve for calibration. The cis and trans miR-155 were measured separately. In the case of 10/1000 (nM) miR-155/P155 concentrations, the peaks of melting curves for trans RNA samples were the same as synthesized miR-155. According to the standard curve of miR-155, the miR-155 concentrations in trans solutions were 14, 34 and 63 aM (10−8 M), indicating that a trace amount of miR-155 transported to the trans side of the pore.
Note S1. Precursor miRNAs (pre-miRNAs) are stem-loop RNAs of  ̃70 nucleotides bearing the 2 nucleotides 3′-overhang as a signature of RNase III-mediated cleavage (Lee, Y., Jeon, K., Lee, J. T., Kim, S., & Kim, V. N. MicroRNA maturation: stepwise processing and subcellular localization. EMBO J 21, 4663-4670 (2002)). It is not known whether the plasma total RNA extract contains pre-miRNAs. However, we have verified that the capture rate for a miR·P is very low if the signal tag in the probe is very short (FIG. 9, FIG. 16 and unpublished data). Therefore, we expected that, even though pre-miRNAs exist, the short overhang will prevent them from trapping in the nanopore.
Note S2. We also compared the time spent for analyzing miRNA via nanopore and qRT-PCR (in Supplementary Information). Our RT-PCR detection can test at most 16 samples at once, including triplicates for each sample with and without spike-in. PCR takes about 5-6 hours, including 1 hour PolyA reaction, 1 hour reverse transcription, 2.5 hours qPCR, plus 1 hour for sample addition. The nanopore method is label free, does not need amplification, and is selective for short nucleic acids fragments. But our current nanopore setup allows detecting one sample at once. The average recording time for miR-155 from human plasma sample was about 90 minutes, collecting  ̃100 events. Therefore, high through-put nanopore methods need developing. This is feasible because both the synthetic nanopore array [Advanced Materials 18, 3149-3153 (2006)] and the protein pore-synthetic pore hybridized system [Nat. Nanotechnol. 5, 874-877 (2010)] have been reported







TABLE 3





Sequences of studied miRNAs and their probes


















mir-155
5′-UUAAUGCUAAUCGUGAUAGGGG-3′




(SEQ ID NO: 1)



 



Pnt
5′-CCCCTATCACGATTAGCATTAA-3′




(SEQ ID NO: 2)



 



P5′-C30
5′-C30-CCCCTATCACGATTAGCATTAA-3′




(SEQ ID NO: 3)



 



P3′-C30
5′-CCCCTATCACGATTAGCATTAA-C30-3′




(SEQ ID NO: 4)



 



P155
5′-C30-CCCCTATCACGATTAGCATTAA-C30-3′




(SEQ ID NO: 5)



 



Let-7a
5′-UGAGGUAGUAGGUUGUAUAGUU-3′




(SEQ ID NO: 6)



 



Pa
5′-C30-AACTATACAACCTACTACCTCA-C30-3′




(SEQ ID NO: 7)



 



Let-7b
5′-UGAGGUAGUAGGUUGUGUGGUU-3′




(SEQ ID NO: 8)



 



Pb
5′-C30-AACCACACAACCTACTACCTCA-C30-3′




(SEQ ID NO: 9)



 



Let-7c
5′-UGAGGUAGUAGGUUGUAUGGUU-3′




(SEQ ID NO: 13)



 



Pc
5′-C30-AACCATACAACCTACTACCTCA-C30-3′




(SEQ ID NO: 14)



 



miR-39
5′-UCACCGGGUGUAAAUCAGCUUG-3′




(SEQ ID NO: 15)



 



P39
5′-C30-CAAGCTGATTTACACCCGGTGA-C30-3′




(SEQ ID NO: 16)

















TABLE 4







Conductance of signature events produced by the miR-155·P155 hybrid


and spike-like short blocks by the translocation of miR-155 or P155 alone










Signature block (miR-155·P155)













Level 1
Level 2
Level 3
Short event

















g (pS)
203
585
110
227



g/g0a
0.15
0.42
0.08
0.16







ag0, the conductance of unoccupied α-hemolysin pore. g0 = 1380 pS at +100 mV in 1M KCl (pH8.0).













TABLE 5







Frequencies and occurrence rate constants of signature events detected


with various probes











Probe
Pnt
P5′-C30
P3′-C30
P155





fsig (min−1)a
0.168 ± 0.042
0.408 ± 0.084
8.31 ± 2.1
11.9 ± 1.2



(n = 5)
(n = 4)
(n = 5)
(n = 6)


kon (M−1s−1)b
2.8 ± 0.6 × 104
6.8 ± 1.3 × 104
1.4 ± 0.3 ×
2.0 ±





106
0.2 × 106





aFrequencies in the presence of 100 nM miR-155 and the probe at +100 mV


bOccurrence rate constant from fsig≈kon [miR-155], where [miR-155] = 100 nM













TABLE 6





Durations of signature events for fully-matched


miRNA•probe hybrids and that with mismatches



















miRNA•Probe
let-7a•Pa
let-7b•Pa
let-7b•Pb
let-7a•Pb





τsig (ms)
155 ± 28
48 ± 11
165 ± 47
24 ± 2


at +120 mV
(n = 7)
(n = 7)
(n = 6)
(n = 5)









p-value
<0.005
<0.005














miRNA•Probe
let-7a•Pa
let-7c•Pa
let-7c•Pc
let-7a•Pc





τsig (ms)
303 ± 45
124 ± 39
343 ± 49
179 ± 38


at +100 mV
(n = 6)
(n = 4)
(n = 6)
(n = 4)









p-value
<0.005
<0.05
















TABLE 7







Areas under ROC curves (AUC) for separation of miRNAs with one


nucleotide difference (let-7a and let-7c) and with two nucleotide


difference (let-7a and let-7b)a













let-7a · Pa
let-7b · Pb

let-7a · Pa
let-7c · Pc
















let-7b · Pa
0.75
n.a.
let-7c · Pa
0.73
n.a.


let-7a · Pb
n.a.
0.83
let-7a · Pc
n.a.
0.71





aThe receiver operating characteristic (ROC) curve is a plot of the true positive rate (sensitivity) against the false positive rate (1-selectivity) for the different possible cutoff points that separate the entire duration distribution into the positive and negative components. In the miRNA detection, the events for fully matched miRNA.probe hybrids were denoted as “positive”, and that for mismatched hybrids as “negative”. The separation accuracy was measured by the area under the ROC curve (AUC). An AUC of 1 represents a perfect separation; an area of 0.5 represents no separation ability. AUC was analyzed online using free software on the world wide web (internet address) rad.jhmi.edu/jeng/javarad/roc/JROCFITi.html.













TABLE 8





Areas under ROC curves (AUC) and optimal cutoff point (OCP)


at various duration ratio and event number ratio a





















τP/τN (s/s) b
1/1
2/1
3/1
4/1
5/1
10/1





AUC
0.51
0.72
0.73
0.76
0.78
0.93


OCP
n.a.
1.33
1.74
1.88
1.98
2.18





NP/NN c
200:800
200:400
200:200
200:150
200:100
200:50





AUC
0.83
0.81
0.76
0.76
0.79
0.78


OCP
1.88
1.88
1.85
1.79
1.81
1.97





a: Both AUC and OCP were calculated from the ROC curves shown in FIG. 16A-B. OCP is a cutoff duration at the maximal value of Youden index. Youden index is defined as {sensitivity + selectivity − 1}, calculated from the ROC curve, and range between 0 and 1. A cutoff duration leading to complete separation of long and short duration distribution results in Youden index = 1, whereas complete overlap gives Youden index = 0. The cutoff duration value that returns the maximum of Youden index, i.e. “optimal” cutoff point (OCP) (Greiner et al., 2000 Preventive Veterinary Medicine 45, 23-41) gives the most accurate separation.


b τP/τN: The duration ratio of the “positive” and “negative” datasets. Each dataset contained 200 exponentially-distributed duration values. The dataset with a longer mean duration was denoted as “positive”; the shorter one as “negative”.


c NP/NN: The event number ratio in the “positive” (NP) versus the “negative” dataset. τP/τN was 5 in this simulation.













TABLE 9





Levels of miR-155 and spiked-in miR-39 in human plasma samples


detected by the nanopore sensor
















a Relative miR-155 level was obtained by normalizing each sample's f155/f39 to the mean f155/f39 of the normal samples 1-6, which was 0.239 as highlighted in the table.






Example 3
Peptide-Guided Selective Detection of MicroRNAs
Rationales. Our long-term goal is developing the nanopore single molecule sensor for accurate detection of target miRNAs in the total RNA extraction from plasma or tissues. The RNA extraction contains numerous and complicated nucleic acids components, at least including miRNAs (both pre-mature and mature miRNAs), mRNAs, tRNAs, other RNAs. All components commonly carry negative charges. Thus, if the target miRNA can be trapped in the pore at an applied voltage, any other component may also be driven to interact with the nanopore, generating non-specific current signals that interfere with the recognition of signature events generated by the target miRNA/probe complex. Here we provide a robust strategy for which only the target miRNA/probe complex can be trapped in the pore, but any other nucleic acids components are prevented from interacting with the pore, therefore greatly improve the both selectivity and sensitivity. This strategy is called peptide-guided selective detection of miRNAs.
Methods. We designed a peptide-PNA (peptide nucleic acid) co-polymer as the probe, as shown in FIG. 18 A. The PNA sequence (in “PNA” bracket in FIG. 18A) has a peptide backbone with side chain nucleobases that are complementary to the entire or partial sequence of the target miRNA (“miRNA (Let-7b, -7c)” in FIG. 18A bracket), and thus serves as the center domain for capturing the target miRNA in the solution. These sequences are provided as follows:
Probe
NH2-Tyr-Gly-Arg-Lys-Lys-Arg-Arg-Gln-Arg-Arg-Arg-AACCACACAA-COOH, where the entire probe has a peptide backbone (i.e. the AACCACACAA portion of the probe comprises a peptide backbone with the indicated AACCACACAA nucleobases); SEQ ID NO:17.
PNA (Center Domain) of Probe:
NH2-AACCACACAA-COOH, where the molecule comprises a peptide backbone with the indicated AACCACACAA nucleobases; (SEQ ID NO: 18).









HIV-TAT:



(SEQ ID NO: 19)



Tyr-Gly-Arg-Lys-Lys-Arg-Arg-Gln-Arg-Arg-Arg






As opposed to oligonucleotide probes, the reporter (or terminal extension) of the new probe is a peptide that carries a series of positively-charged amino acids (“Peptide reporter” in FIG. 18A bracket) and the center domain is a peptide nucleic acid comprising nucleotides that are complementary to the target nucleic acid. When there are a sufficient number of positively-charged amino acids in the reporter or terminal extension portion of the probe, the net charges of the miRNA/probe complex are still positive such that when the target miRNA binds to the PNA (peptide nucleic acid) domain of the probe the entire miRNA/probe complex forms a strong dipole molecule. We have engineered a nanopore with a negatively-charged residue ring at the trans opening of the pore (S. aureus alpha-hemolysin comprising a K131D mutation). The wild-type S. aureus alpha-hemolysin peptide sequence (National Center for Bioinformatics Accession NO. AAA26598.1); is:







(SEQ ID NO: 21)


ADSDINIKTGTTDIGSNTTVKTGDLVTYDKENGMHKKVFYSFIDDKNHN


 


KKLLVIRTKGTIAGQYRVYSEEGANKSGLAWPSAFKVQLQLPDNEVAQI


 


SDYYPRNSIDTKEYMSTLTYGFNGNVTGDDTGKIGGLIGANVSIGHTLK


 


YVQPDFKTILESPTDKKVGWKVIFNNMVNQNWGPYDRDSWNPVYGNQLF


 


MKTRNGSMKAADNFLDPNKASSLLSSGFSPDFATVITMDRKASKQQTNI


 


DVIYERVRDDYQLHWTSTNWKGTNTKDKWTDRSSERYKIDWEKEEMTN






The variant S. aureus K131D alpha-hemolysin peptide sequence is:







(SEQ ID NO: 22)


ADSDINIKTGTTDIGSNTTVKTGDLVTYDKENGMHKKVFYSFIDDKNHN


 


KKLLVIRTKGTIAGQYRVYSEEGANKSGLAWPSAFKVQLQLPDNEVAQI


 


SDYYPRNSIDTKEYMSTLTYGFNGNVTGDDTGDIGGLIGANVSIGHTLK


 


YVQPDFKTILESPTDKKVGWKVIFNNMVNQNWGPYDRDSWNPVYGNQLF


 


MKTRNGSMKAADNFLDPNKASSLLSSGFSPDFATVITMDRKASKQQTNI


 


DVIYERVRDDYQLHWTSTNWKGTNTKDKWTDRSSERYKIDWEKEEMTN






Therefore, the positively-charged peptide domain of the probe dipole will be both pushed by the positive voltage (cis grounded) and attracted by the negative ring at the trans opening, guiding the trapping of the miRNA/complex into the n-barrel of the pore. At the positive voltage, any other free nucleic acids components will be repulsed from entering the pore due to the negative charge carried. This significantly reduces signals by free RNA components, and most observed events are either due to the trapping of the miRNA/probe complex or the translocation of the probe. The use of peptide-PNA probe enables selective detection of the target miRNA.
Results. FIG. 18A shows the diagram of the miRNA/probe complex. The bracketed miRNA is target miRNA Let-7b. The bracketed probe P7b has a bracketed “Peptide Reporter” part and a bracketed “PNA” (peptide nucleic acid) part. The PNA is for capturing Let-7b, and the bracketed “Peptide reporter” is a positively-charged peptide corresponding the sequence of HIV-TAT, which contains +8e contributed by arginines and lysines. FIG. 18B shows events for translocation of the peptide-PNA probe, P7b. The characteristic events last for 3 ms and reduce the current to 10 pA at +180 mV. FIG. 18C shows no block events can be observed with free miRNA let-7b (without probe) in the solution at +180 mV.
FIG. 18D shows signature events for the trapping of the let-7b/P7b complex. These events characteristically last for 100 ms and reduce the current to 57 pA at +180 mV, completely different that for the probe. FIG. 18E shows that Let-7c, which has two different nucleotides from Let-7b, cannot bind to PNA of the probe P7b, therefore does not generate signature events as in FIG. 18C. Almost all observed events are due to the probe itself.
FIG. 18F compares the duration-amplitude property for P7b binding to Let-7b (fully match, two separate clusters without overlay) and Let-7c (2 mismatches, two clusters fully overlay). This suggests an accuracy of almost 100% in differentiating sequence-similar miRNAs with two different nucleotides.
Utilization of signature events to understand various molecular processes in the nanopore for biosensing applications is also provided. In FIG. 19A, we observed a novel type of three-level current pattern when employing HP-C30 with a hairpin at the 3′-end of short strand. Its Level 1 and Level 2 are consistent with the unzipping of HP-C30 and translocation of the unzipped short strand from the nanocavity to the β-barrel. However, the duration of Level 1′ was drastically prolonged by 80 folds to 15±1.9 ms, compared to the target without a hairpin. The prolonged Level 1′ is in agreement with the unzipping of hairpin prior to threading in the β-barrel. As many DNA or RNA structures such as aptamers contain hairpins, we can use this system and the signature events to study these structures and study their binding interaction with their protein targets. In FIG. 19B, we also demonstrated a new multi-level current pattern when using SA-C30 attached with a streptavidin at the 3′-end of the short strand. Again, both the fully-blocked Level 1 and partially-blocked level 2 are consistent with the unzipping of SA-C30 and translocation of the unzipped short strand from the nanocavity to the β-barrel. The current stayed at Level 1′ for minutes until it was forced to recover by a negative voltage. The long term Level 1′ can be interpreted by that although the short strand of SA-C30 moves into the β-barrel after unzipping, its translocation is prevented by the attached large streptavidin. This result suggested the potential of using signature events for protein detection. In FIG. 19C, we demonstrated that the complex can be sequentially unzipped in the nanopore in two steps when using a short oligonucleotide to link two DNAs. The unzipping of the two DNAs can be clearly revealed by the two Level 2 states.
Conclusion. The peptide-PNA probe enables 1) selective detection of the target miRNA, 2) greatly enhanced accuracy in differentiating sequence-similar miRNAs.
While the invention has been described in connection with specific embodiments thereof, it will be understood that the inventive device is capable of further modifications. This patent application is intended to cover any variations, uses, or adaptations of the invention following, in general, the principles of the invention and including such departures from the present disclosure as come within known or customary practice within the art to which the invention pertains and as may be applied to the essential features herein before set forth and as follows in scope of the appended claims.
REFERENCES


    1. Silvestri, G. A., Alberg, A. J., & Ravenel, J. The changing epidemiology of lung cancer with a focus on screening. BMJ 339, (2009).
    2. Carthew, R. W. & Sontheimer, E. J. Origins and Mechanisms of miRNAs and siRNAs. Cell 136, 642-655 (2009).
    3. Inui, M., Martello, G., & Piccolo, S., MicroRNA control of signal transduction. Nat Rev Mol. Cell Biol 11, 252-263 (2010).
    4. Garzon, R., Calin, G. A., and Croce, C. M. MicroRNAs in cancer. 60, 167-179. 2009.
    5. Ortholan, C., Puissegur, M. P., Ilie, M., Barbry, P., Mari, B., & Hofman, P. MicroRNAs and lung cancer: New oncogenes and tumor suppressors, new prognostic factors and potential therapeutic targets. Curr. Med. Chem. 16, 1047-1061 (2009).
    6. Calin, G. A. & Croce, C. M. MicroRNA signatures in human cancers. Nat Rev Cancer 6, 857-866 (2006).
    7. Mitchell, P. S., Parkin, R. K., Kroh, E. M., Fritz, B. R., Wyman, S. K., Pogosova-Agadjanyan, E. L., Peterson, A., Noteboom, J., O'Briant, K. C., Allen, A., Lin, D. W., Urban, N., Drescher, C. W., Knudsen, B. S., Stirewalt, D. L., Gentleman, R., Vessella, R. L., Nelson, P. S., Martin, D. B., & Tewari, M. Circulating microRNAs as stable blood-based markers for cancer detection. Proc. Natl. Acad. Sci. U.S.A. 105, 10513-10518 (2008).
    8. Rabinowits, G., GerÃ§el-Taylor, C., Day, J. M., Taylor, D. D., & Kloecker, G. H. Exosomal microRNA: A diagnostic marker for lung cancer. Clin. Lung Cancer 10, 42-46 (2009).
    9. Kosaka, N., Iguchi, H., Yoshioka, Y., Takeshita, F., Matsuki, Y., & Ochiya, T. Secretory mechanisms and intercellular transfer of microRNAs in living cells. J Biol Chem (2010).
    10. Rosell, R., Wei, J., & Taron, M. Circulating MicroRNA Signatures of Tumor-Derived Exosomes for Early Diagnosis of Non-Small-Cell Lung Cancer. Clin. Lung Cancer 10, 8-9 (2009).
    11. Thomson, J. M., Parker, J., Perou, C. M., & Hammond, S. M. A custom microarray platform for analysis of microRNA gene expression. Nat Methods 1, 47-53 (2004).
    12. Chen, C., Ridzon, D. A., Broomer, A. J., Zhou, Z., Lee, D. H., Nguyen, J. T., Barbisin, M., Xu, N. L., Mahuvakar, V. R., Andersen, M. R., Lao, K. Q., Livak, K. J., & Guegler, K. J. Real-time quantification of microRNAs by stem-loop RT-PCR. Nucleic Acids Res. 33, (2005).
    13. Li, W. & Ruan, K. MicroRNA detection by microarray. Anal. Bioanal. Chem. 394, 1117-1124 (2009).
    14. Hunt, E. A., Goulding, A. M., & Deo, S. K. Direct detection and quantification of microRNAs. Anal. Biochem. 387, 1-12 (2009).
    15. Kloosterman, W. P., Wienholds, E., de Bruijn, E., Kauppinen, S., & Plasterk, R. H. In situ detection of miRNAs in animal embryos using LNA-modified oligonucleotide probes. Nat Methods 3, 27-29 (2006).
    16. Neely, L. A., Patel, S., Garver, J., Gallo, M., Hackett, M., McLaughlin, S., Nadel, M., Harris, J., Gullans, S., & Rooke, J. A single-molecule method for the quantitation of microRNA gene expression. Nat Methods 3, 41-46 (2006).
    17. Bayley, H. & Jayasinghe, L. Functional engineered channels and pores—(Review). Molecular Membrane Biology 21, 209-220 (2004).
    18. Ma, L. & Cockroft, S. L. Biological nanopores for single-molecule biophysics. Chembiochem 11, 25-34 (2010).
    19. Movileanu, L. Interrogating single proteins through nanopores: challenges and opportunities. Trends Biotechnol. 27, 333-341 (2009).
    20. Bayley, H. Sequencing single molecules of DNA. Current Opinion in Chemical Biology 10, 628-637 (2006).
    21. Branton, D., Deamer, D. W., Marziali, A., Bayley, H., Benner, S. A., Butler, T., Di Ventra, M., Garaj, S., Hibbs, A., Huang, X., Jovanovich, S. B., Krstic, P. S., Lindsay, S., Ling, X. S., Mastrangelo, C. H., Meller, A., Oliver, J. S., Pershin, Y. V., Ramsey, J. M., Riehn, R., Soni, G. V., Tabard-Cossa, V., Wanunu, M., Wiggin, M., & Schloss, J. A. The potential and challenges of nanopore sequencing. Nature Biotechnology 26, 1146-1153 (2008).
    22. Kasianowicz, J. J., Brandin, E., Branton, D., & Deamer, D. W. Characterization of individual polynucleotide molecules using a membrane channel. Proc. Natl. Acad. Sci. U.S.A. 93, 13770-13773 (1996).
    23. Akeson, M., Branton, D., Kasianowicz, J. J., Brandin, E., & Deamer, D. W. Microsecond time-scale discrimination among polycytidylic acid, polyadenylic acid, and polyuridylic acid as homopolymers or as segments within single RNA molecules. Biophys. J 77, 3227-3233 (1999).
    24. Meller, A., Nivon, L., & Branton, D. Voltage-driven DNA translocations through a nanopore. Phys Rev Lett 86, 3435-3438 (2001).
    25. Deamer, D. W. & Branton, D. Characterization of nucleic acids by nanopore analysis. Acc. Chem. Res. 35, 817-825 (2002).
    26. Nakane, J., Akeson, M., & Marziali, A. Nanopore sensors for nucleic acid analysis. J Phys Condens Matter 15, (2003).
    27. Mitchell, N. & Howorka, S. Chemical tags facilitate the sensing of individual DNA strands with nanopores. Angew. Chem. Int. Ed. 47, 5565-5568 (2008).
    28. Rabinowits, G., GerÃ§el-Taylor, C., Day, J. M., Taylor, D. D., & Kloecker, G. H. Exosomal microRNA: A diagnostic marker for lung cancer. Clin. Lung Cancer 10, 42-46 (2009).
    29. Rosell, R., Wei, J., & Taron, M. Circulating microRNA signatures of tumor-derived exosomes for early diagnosis of non-small-cell lung cancer. Clin. Lung Cancer 10, 8-9 (2009).
    30. Meller, A., Nivon, L., Brandin, E., Golovchenko, J., & Branton, D. Rapid nanopore discrimination between single polynucleotide molecules. Proc. Natl. Acad. Sci. U.S.A. 97, 1079-1084 (2000).
    31. Sauer-Budge, A. F., Nyamwanda, J. A., Lubensky, D. K., & Branton, D. Unzipping kinetics of double-stranded DNA in a nanopore. Phys Rev Lett 90, (2003).
    32. Mathés, J., Visram, H., Viasnoff, V., Rabin, Y., & Meller, A. Nanopore unzipping of individual DNA hairpin molecules. Biophys. J 87, 3205-3212 (2004).
    33. Maglia, G., Restrepo, M. R., Mikhailova, E., & Bayley, H. Enhanced translocation of single DNA molecules through +/−-hemolysin nanopores by manipulation of internal charge. Proc. Natl. Acad. Sci. U.S.A. 105, 19720-19725 (2008).
    34. Shim, J. W., Tan, Q., & Gu, L. Q. Single-molecule detection of folding and unfolding of a single G-quadruplex aptamer in a nanopore nanocavity. Nucleic Acids Res. 37, 972-982 (2009).
    35. Butler, T. Z., Gundlach, J. H., & Troll, M. A. Determination of RNA orientation during translocation through a biological nanopore. Biophys. J. 90, 190-199 (2006).
    36. Kasianowicz, J. J., Brandin, E., Branton, D., & Deamer, D. W. Characterization of individual polynucleotide molecules using a membrane channel. Proc. Natl. Acad. Sci. U.S.A. 93, 13770-13773 (1996).
    37. Mathés, J., Aksimentiev, A., Nelson, D. R., Schulten, K., & Meller, A. Orientation discrimination of single-stranded DNA inside the Î±-hemolysin membrane channel. Proc. Natl. Acad. Sci. U.S.A. 102, 12377-12382 (2005).
    38. Purnell, R. F., Mehta, K. K., & Schmidt, J. J. Nucleotide identification and orientation discrimination of DNA homopolymers immobilized in a protein nanopore. Nano Lett. 8, 3029-3034 (2008).
    39. Wanunu, M., Morrison, W., Rabin, Y., Grosberg, A. Y., & Meller, A. Electrostatic focusing of unlabelled DNA into nanoscale pores using a salt gradient. Nature Nanotechnology 5, 160-165 (2010).
    40. Cho, W. C. Role of miRNAs in lung cancer. Expert Rev Mol. Diagn. 9, 773-776 (2009).
    41. Landi, M. T., Zhao, Y., Rotunno, M., Koshiol, J., Liu, H., Bergen, A. W., Rubagotti, M., Goldstein, A. M., Linnoila, I., Marincola, F. M., Tucker, M. A., Bertazzi, P. A., Pesatori, A. C., Caporaso, N. E., McShane, L. M., & Wang, E. MicroRNA expression differentiates histology and predicts survival of lung cancer. Clin. Cancer Res 16, 430-441 (2010).
    42. Patnaik, S. K., Kannisto, E., Knudsen, S., & Yendamuri, S. Evaluation of microRNA expression profiles that may predict recurrence of localized stage I non-small cell lung cancer after surgical resection. Cancer Res 70, 36-45 (2010).
    43. Yanaihara, N., Caplen, N., Bowman, E., Seike, M., Kumamoto, K., Yi, M., Stephens, R. M., Okamoto, A., Yokota, J., Tanaka, T., Calin, G. A., Liu, C. G., Croce, C. M., & Harris, C. C. Unique microRNA molecular profiles in lung cancer diagnosis and prognosis. Cancer Cell 9, 189-198 (2006).
    44. Shi, R. & Chiang, V. L. Facile means for quantifying microRNA expression by real-time PCR. BioTechniques 39, 519-524 (2005).
    45. Kim, M. J., Wanunu, M., Bell, D. C., & Meller, A. Rapid fabrication of uniformly sized nanopores and nanopore arrays for parallel DNA analysis. Advanced Materials 18, 3149-3153 (2006).

## Claims

1. A method of detecting an oligonucleotide in a sample with a nanopore system, comprising the steps of:
a) mixing a sample suspected of comprising a target oligonucleotide and a probe, said probe comprising a central domain with a complementary sequence to said target oligonucleotide and a terminal extension tagged to at least one of a 3′ terminal or a 5′ terminal of the central domain, to produce a sample mixture wherein said central domain of said probe hybridizes to said target oligonucleotide,
b) applying a voltage to said sample mixture in a cis compartment of a dual chamber nanopore system sufficient to drive translocation of said hybridized probe and target oligonucleotide through a nanopore of said system by an unzipping process, and
c) analyzing an electrical current pattern in said nanopore system over time,
wherein a presence of said target oligonucleotide in said sample is indicated by at least one signature electrical current block that is distinct from a background electrical current block that occurs with said sample alone or said probe alone.

2. The method of claim 1, wherein said signature electrical current block comprises at least one of: i) a current block of different duration than a background current block; ii) a different number of distinct current blockade levels than a background current block; iii) a different order of occurrence of current blockade levels than a background current block; iv) a blockade level having a different current amplitude than a background current block; v) a current amplitude of each blockade level different than a current amplitude of a background current block; or any combination of (i), (ii), (iii), (iv), or (v).

3. The method of claim 1, wherein at least one terminal extension is a charged polypeptide or an oligonucleotide.

4. The method of claim 1, wherein said terminal extension comprises a charged polypeptide containing at least two positively charged amino acids, at least two aromatic amino acid residues, or at least two positively charged amino acids and at least two aromatic amino acid residues.

5. The method of claim 4, wherein said charged polypeptide provides for a net positive charge when said probe is hybridized to said target oligonucleotide.

6. The method of claim 1, wherein said nanopore comprises a negatively charged residue ring at a trans opening of the nanopore.

7. The method of claim 1, wherein said nanopore comprises an alpha-hemolysin variant comprising a K131D, K131E, or K131H amino acid substitution.

8. The method of claim 1, wherein steps (a) and (b) or step (b) are performed at a pH value of less than 7.

9. The method of claim 1, wherein said target oligonucleotide is a DNA molecule or an RNA molecule.

10. The method of claim 9, wherein said RNA molecule is an miRNA or a fragment thereof, and wherein said fragment comprises at least 15 nucleotides of said RNA molecule.

11. The method of claim 1, wherein a frequency of said signature electrical current block is used to quantitate an amount of said target oligonucleotide in said sample.

12. The method of claim 1, wherein a non-target oligonucleotide is spiked into said sample and a frequency of electrical current blocks attributable to said non-target oligonucleotide is used to normalize a frequency of the signature electrical current block.

13. The method of claim 1, wherein said central domain of said probe comprises at least 15 nucleotides or nucleobases that are complementary to the target oligonucleotide.

14. A method for detecting two or more distinct target oligonucleotides that differ by at least one nucleotide in a sample with a nanopore system, comprising the steps of:
a) mixing a sample suspected of comprising said distinct target oligonucleotides and a probe, said probe comprising a central domain with a fully or partially complementary sequence to said target oligonucleotides and a terminal extension tagged to at least one of a 3′ terminal or a 5′ terminal of the central domain, to produce a sample mixture wherein said central domain of said probe hybridizes to said target oligonucleotides,
b) applying a voltage to said sample mixture in a cis compartment of a dual chamber nanopore system sufficient to drive translocation of hybridized probe and target oligonucleotides through a nanopore of said system by an unzipping process, and
c) analyzing an electrical current pattern in said nanopore system over time, wherein a presence of two or more distinct target oligonucleotides in said sample is indicated by two or more distinct signature electrical current blocks.

15. The method of claim 14, wherein said two or more distinct signature electrical current blocks comprise at least one of: i) different signature electrical current block durations (τsig); ii) a different number of distinct current blockade levels; iii) a different order of occurrence of current blockade levels; iv) a different current amplitude; v) a different current amplitude of each blockade level; or any combination of (i), (ii), (iii), (iv), or (v).

16. The method of claim 14, wherein said target oligonucleotides are DNA molecules or RNA molecules.

17. The method of claim 16, wherein said RNA molecules are miRNAs or fragments thereof, and wherein said fragments comprise at least 15 nucleotides of said RNA molecules.

18. The method of claim 14, wherein the terminal extension is a charged polypeptide or an oligonucleotide.

19. The method of claim 18, wherein the charged polypeptide comprises at least two positively charged amino acids at least two aromatic amino acid residues, or at least two positively charged amino acids and at least two aromatic amino acid residues.

20. The method of claim 19, wherein said charged polypeptide provides for a net positive charge when said probe is hybridized to said target oligonucleotides.

21. The method of claim 19, wherein said nanopore comprises a negatively charged residue ring at a trans opening of the nanopore.

22. The method of claim 19, wherein said nanopore comprises an alpha-hemolysin variant comprising a K131D, a K131E, or a K131H amino acid substitution.

23. The method of claim 14, wherein said central domain comprises at least 15 nucleotides or nucleobases that are complementary to the target oligonucleotides.

24. A kit comprising, (a) a probe comprising a central domain with a fully or partially complementary sequence to a target oligonucleotide and a terminal extension tagged to at least one of a 3′ terminal or a 5′ terminal of the central domain, wherein said probe comprises at least one terminal extension that is a charged polypeptide, and wherein said central domain of said probe can hybridize to the target oligonucleotide; and (b) a nanopore.

25. The kit of claim 24, wherein said central domain comprises a peptide nucleic acid probe.

26. The kit of claim 24, wherein the charged polypeptide comprises at least two positively charged amino acids, at least two aromatic amino acid residues, or at least two positively charged amino acids and at least two aromatic amino acid residues.

27. The kit of claim 26, wherein said charged polypeptide provides for a net positive charge when said probe is hybridized to the target oligonucleotide.

28. The kit of claim 24, wherein said nanopore comprises a negatively charged residue ring at a trans opening of the nanopore.

29. The kit of claim 24, wherein said nanopore comprises an alpha-hemolysin variant comprising a K131D, a K131E, or a K131H amino acid substitution.