This is an open access article distributed under the terms of the Creative Commons Attribution License (
Aptamers are nucleic acids selected for their ability to bind to molecules of interest and may provide the basis for a whole new class of medicines. If the aptamer is simply a dsDNA molecule with a ssDNA overhang (a "sticky" end) then the segment of ssDNA that complements that overhang provides a known binding target with binding strength adjustable according to length of overhang.
Two bifunctional aptamers are examined using a nanopore detector. They are chosen to provide sensitive, highly modulated, blockade signals with their captured ends, while their un-captured regions are designed to have binding moieties for complementary ssDNA targets. The bifunctional aptamers are duplex DNA on their channel-captured portion, and single-stranded DNA on their portion with binding ability. For short ssDNA, the binding is merely to the complementary strand of DNA, which is what is studied here – for 5-base and 6-base overhangs.
A preliminary statistical analysis using hidden Markov models (HMMs) indicates a clear change in the blockade pattern upon binding by the single captured aptamer. This is also consistent with the hypothesis that significant conformational changes occur during the annealing binding event. In further work the objective is to simply extend this ssDNA portion to be a well-studied ~80 base ssDNA aptamer, joined to the same bifunctional aptamer molecular platform.
Our nanopore detector is biologically based and uses a protein, the α-hemolysin (α-HL) toxin produced by the bacterium
There are important distinctions in how a nanopore detector can function: direct vs. indirect measurement of static, stationary, dynamic (possibly modulated), or non-stationary channel blockades (see [
Left Panel: A lipid bilayer supports the alpha-hemolysin heptamer that creates the pore, or channel, used to collect the data, as shown left. The bilayer is established and supported across an aperture that typically provides 5–25 um in effective bilayer diameter, and generally greater than 1 um. Right Panel: The assembled alpha-Hemolysin pore shown to scale, with a captured dsDNA molecule. As shown, the double stranded form is too wide to pass through the pore, while a single strand may pass through. Bottom Panel: One-second blockade patterns of four DNA hairpins, part of a test set of nine base-pair hairpins, with 4 dT hairpin loops, that have been studied extensively. The molecules only differ in their terminal base-pairs, yet their channel current blockade signals, "signatures", are easily resolved [4].
Classification software scheme and data flowchart. In real-time processing, the data from the detector is processed by wavelet FSA and stationarity analyses. From the HMM profiling, the SVM is used to classify the data, performing clustering analysis. The SVM is able to discriminate among known signal patterns with up to 99.9% accuracy
The signal processing architecture (Fig.
Aptamers are essentially synthetically-derived, single stranded RNA or DNA molecules up to ~80 oligonucleotides in length with a high affinity towards bonding to specific targets. Consequently, aptamers may be ligands of proteins or other large biomolecules while having modified functional groups either during production or for added utility. In 1990, it was independently reported that a new method dubbed SELEX (Systematic Evolution of Ligands by EXponential Enrichment) provided a process of producing aptamers from random DNA or RNA libraries ([
In a short number of years, there has been a growing preference for the use of aptamers over antibodies in a variety of different uses. The key to the success is how well aptamers do the same job in a diagnostic setting, as well as how versatile aptamers are in their design and usage. There is little to no conflict between design and function due to its origins in DNA and RNA, yet the structure is completely artificial. Since 1999, it has been reported that aptamers are being made in vivo as opposed to strictly in vitro [
Aptamers may bind to multiple molecules simultaneously, thus allowing certain molecules or functional groups to bind more specifically, e.g., fluorescein, biotin bound to aptamers which target a specific protein. In fact, they can discriminate between small functional groups (i.e., methyl vs OH-, and D, L-enantiomers), and the specificity rivals that of antibodies and surpasses in some cases [
Properties of aptamers include the ability of a given structure to have multiple, stable, structural configurations (especially compared to peptides) due to Watson-Crick base pairing. Aptamers may bind with anything from dyes, drugs, peptides, proteins, metal ions, antibodies, and enzymes for example. The values of Kd range between ~pML-1 to ~nML-1, better than that of antibodies [
The production of DNA-based aptamers does not depend on the inclusion of particular functional groups [
Aptamers are now replacing antibodies as detection reagents due to having several advantages over antibodies: versatility, the creation of a lab-on-a-chip to process, low detection limits, simpler reactions to perform, diversity and specificity of aptamer-target binding properties [
Since aptamers have high binding affinity and target specificity, there have been several investigations as to their structure and kinetic properties using techniques such as NMR and UV-Vis spectroscopy [
Experiments are performed with a linear molecule with a bulge in the center. To one side of the bulge is the blunt-ended stem sequence like that used in one of our DNA hairpin controls, where the bulge is now in the position of the hairpin's loop. To the other side of the bulge is a cap-section of base-pairs followed by an overhang section of length five bases. A similar set of experiments is performed with the "Y-aptamer", a Y-shaped DNA complex with one arm of the Y with an overhang (6 Ts), while the other arm is capped with a 4 dT loop. The base of the Y is a stem of 10 base-pairs length, prior to the Y-nexus of the molecule. Here the Y-nexus is in the place of the bulge, or the hairpin loop. Nine or ten base-pairs is approximately the length in dsDNA from the mouth of the channel to its limiting aperture. The significance of this length in the modeling is due to its delicate placement of the end of the captured molecule over the high electrophoretic field strength zone near the limiting aperture of the channel, permitting operation in transduction mode [
The linear duplex DNA molecule, with bulge, and ssDNA overhang, is shown captured by the channel in Figure
Translocation of a single DNA strand. (a) A crude molecular schematic illustrating a captured pseudo-aptamer before binding to the 5 nucleotide overhang complement. If the two strands were to separate, the strands would pass through the vestibule one strand at a time, as an aptamer melting event and ssDNA translocation. (b) Raw data shows the ssDNA translocation signals (most involving molecules that are already single-stranded upon capture).
(a) Collection of toggle signals from the captured pseudo-aptamer. The molecule appears as the one in Fig. 3(a), unattached to the 5 nucleotide ssDNA. (b) Collection of toggle signals from the pseudo-aptamer + 5 nucleotide ssDNA. The bound complex appears as in Fig. 3(a), with the 5 nucleotide sequence joined to the overhang section. One distinctive change only observed in the blockade signals after binding agent introduced, aside from the level dwell-time changes, are the much higher frequency upward "spike" transitions, from the lower level to the upper level.
Level duration feature extraction for the 6 T pseudo-aptamer. First feature extraction for the pseudo-aptamer bound to its complementary 5 base pair sequence. The y-axis is the transition probability, while the x-axis is the HMM Profile: 150 feature vectors components are obtained from the 50-state HMM-EM/Viterbi implementation in [4] are: the 50 dwell percentage in the different blockade levels (from the Viterbi trace-back states), the 50 variances of the emission probability distributions associated with the different states, and the 50 merged transition probabilities from the primary and secondary blockade occupation levels (fits to two-state dominant modulatory blockade signals).
Using our HMM-based cheminformatics software, we can determine the median values of the dwell time per level. The dwell time corresponds to the interval of time before transition to another level. A capture signal has noticeably different current value (between ~20–70 pA) than the open channel current (~120 pA). A blockade (or toggle) will transition between two or more current values, given the difference between the levels is greater than the signal noise level. (Notation: the assignment of levels is determined by increasing mean current value. For example, the lowest level may have a mean current around 35 pA, the intermediate level around 40 pA, and the upper level around 55 pA.) The median value can be interpreted as the half-life in biophysics terms or as ~1/koff in biochemical terms. The binding constant Kd is proportional to kon/koff. Figure
(a) Dwell-time distributions in the dominant blockade levels for the pseudo-aptamer. The vertical dotted lines indicate the median dwell time for two levels. (b) Same analysis for the pseudo-aptamer bound to its corresponding 5 base pair sequence, with the emergence of third level not present in the unbound pseudo-aptamer. The two vertical dotted lines indicate the median dwell times for the first and second levels, while the third level median dwell time value extends to 9.54 ms (plot scale has been truncated for space).
For the unbound pseudo-aptamer, the median values of these distributions are 1.78 ms and 3.12 ms (lower, upper levels respectively – see Fig.
Standard 150-component HMM-based feature extraction for the T6-Y10T1-GC (mfold diagram for molecule shown top), and for the T6-Y10T1-GC +T6 Y-aptamer blockade signals. (See Fig. 5 Caption for explanation of HMM features.)
The dwell times at each level for the Y-aptamers. (a) Kinetic feature extraction associated with analysis of the signal shown in Figure 7 for the Y-aptamer T6-Y10T1-GC. The vertical dotted lines stand in place of the median dwell times. (b) Same analysis for the Y-aptamer bound to the additional T6 appendage, with the disappearance of third level present in the unbound Y-aptamer. Again, the vertical lines stand in place of the median dwell times for each level.
The blockade stationary statistics or phases thereof, reveal information about the kinetics of the biopolymer resulting from interactions with another molecule, from interactions with the channel, or from undergoing conformational changes. LabVIEW Automation software is used to help manage the feedback linkage between patch-clamp amplifier measurements and in-house cheminformatics software. This has been used to do demonstrate molecular identification in the first 100 ms of capture, with return of classification information to the control of the amplifier – for voltage-controlled sample ejection if desired [
It is important to remember that the statistics obtained are for a single molecule, as opposed to a large collection of molecules. While diagnostic methods which employ aptamers in a dye-quenching process involve thousands if not millions of molecules (in a small volume) to obtain statistically significant results, the signal acquired by the nanopore detector is dependent on the statistics of small numbers. The process of detecting aptamer-beacons takes only a short while due to the large numbers of molecules contributing to the overall signal. If the diagnostic were carried out with a single fluorescing molecule bound to an aptamer, the photons absorbed and emitted would be few and far between, resulting in data completely buried by noise. Therefore, the statistics determined by the analysis of a single molecule are a snapshot of the conditions present at the time of data collection, as opposed to averaging numerous signals collected over time together.
This brings up the discussion of the fundamental physics involved in the examination of a single molecule. Does a single biomolecule behave in a way that complies with classical or quantum mechanics? A collection of gas molecules can be classical, while a single gas molecule emits radiation generally according to allowed transitions in quantum mechanics. Classical or quantum can affect the method of data collection and resulting statistical analysis. Our analysis of a single aptamer (and larger biomolecules) so far remains in the classical regime, since the molecules are thoroughly decohered in their quantum phase, immersed as they are in a highly interactive, stationary, ionic flow. If the data of interest were a measure of change occurring when an attached dye emitted or absorbed a photon, however, the measurements utility might be limited by quantum statistical noise.
The time-scale for the level-switch, or molecular state-changes, is clearly different before and after introduction of the complementary 5ss DNA binding partner. Furthermore the dominant levels increase from two to three, indicating a large-scale change to the molecule, possibly due to conformational change in the molecule itself upon binding. Similar observations hold for the annealing experiments with the Y-aptamer and its 5'-TTTTTT-3' overhang. Refinements of single-molecule binding information like this will lead to a much more informed aptamer design process. For the moment, a prototype of this informed aptamer design process is simply in the selection for clearer observation of on/off binding: Right now we see clear examples of on-binding (with a variety of controls to reassure us that that is what we are seeing), but nothing conclusive has been done with observation of off-binding. Once we can get an annealing or other binding interaction with kon and koff statistics, we can pursue a number of fundamental experiments as to the nanopore detectors overall capabilities in this regard. The ability to observe binding, and select for binding strength of interest, could aid in the identification of tumor-specific molecular markers, for example, which is a powerful tool in cancer diagnostics, and targeting of tumor specific pathways is the best hope for developing non-toxic and efficient anti cancer therapies. Targeting of cancer cells relies on the development of molecules, aptamers with toxins in tow, perhaps, that are suited to in vivo applications and that have the required affinity, specificity and favorable pharmacokinetic properties.
Each experiment is conducted using one α-hemolysin channel inserted into a diphytanoyl-phosphatidylcholine/hexadecane bilayer across a 25-micron-diameter horizontal Teflon aperture, as described previously [
The five DNA hairpins studied in the nanopore prototype experiment [
The linear aptamer with bulge consists of annealing the following two ssDNA strands:
(1) 5'-GAGGCTTGG TTT CAATAGGTA-3'
(2) 5'-ATTG TTT CCAAGCCTC-3'
The overhang region has sequence (3) as part of (1). The complementary 5 nucleotide ssDNA sequence (4):
(3) 5'-AGGTA-3'
(4) 5'-TACCT-3'
Likewise, the Y-aptamer DNA molecule consists of a three-way DNA junction created by annealing two DNA molecules:
(1) 5'-CTCCGTCGAC GAGTTTATAGAC TTTTTT-3'
(2) 5'-GTCTATAAACTC GCAGTCATGCTTTTGCATGACTGC GTCGACGGAG-3'
For the resulting Y-aptamer, one of the junctions' arms terminate in a 4 dT-loop and the other arm has a 6 T overhang in place of a 4 dT-loop. The blunt ended arm has to be carefully designed such that when it is captured by the nanopore it produces a toggling blockade. One of the arms of the Y-shaped aptamer (Y-aptamer) has a TATA sequence, and is meant to be a binding target for TBP. In general, any transcription factor binding site could be studied (or verified) in this manner. Similarly, transcription factor could be verified, or the efficacy of a synthetic transcription factor could be examined.
The features extracted for each blockade signal are as follows: the first 50 components are the level occupation probabilities for the blockade signal in the range 20% to 70% of the open channel current (with bin size at 1%). The blockade signal entirely resides in this blockade range, which is minimized for faster HMM processing. The middle 50 components are the variances of the emission probabilities for each of the 50 blockade levels (normalized as an L1 vector). The last 50 components shows a compression of the 2500 transition probabilities to the 50 values obtained from merging the transition probabilities from the two strongest levels indicated in the analysis of the first 50 components (with weighting as indicated by their support in the level occupation probabilities). The software is accessible via a web interface at
The authors declare that they have no competing interests.
The initial submission was written by KT and SWH, with revisions by SWH. Most of the data was gathered by KT after being trained on nanopore detector operations by IA and EM. IA and EM also helped in gathering the data. EM helped in constructing the nanopore device. Experimental design and pattern recognition software contributed by SWH.
Federal funding was provided by NIH K-22 (PI, 5K22LM008794), NIH NNBM R-21 (co-PI), and NIH R-01 (sub-award). State funding was provided from a LaBOR Enhancement (PI), a LaBOR Research Competitiveness Subcontract (PI), and a LaBOR/NASA LaSPACE Grant (PI). Thanks to Alex Ortiz for help gathering data on the linear-bulge aptamer. Funding also provided by New Orleans Children's Hospital and the University of New Orleans Computer Science Department.
This article has been published as part of