help button home button Biophys. J.
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
Originally published as Biophys J. BioFAST on August 3, 2007.
doi:10.1529/biophysj.107.107003
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
biophysj.107.107003v1
93/9/3229    most recent
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Butler, T. Z.
Right arrow Articles by Troll, M.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Butler, T. Z.
Right arrow Articles by Troll, M.
Biophysical Journal 93:3229-3240 (2007)
© 2007 The Biophysical Society

Ionic Current Blockades from DNA and RNA Molecules in the {alpha}-Hemolysin Nanopore

Tom Z. Butler *, Jens H. Gundlach * and Mark Troll {dagger}

* Department of Physics, and {dagger} Department of Microbiology, University of Washington, Seattle, Washington

Correspondence: Address reprint requests to T. Z. Butler, E-mail: twb2{at}u.washington.edu.


   ABSTRACT
TOP
 ABSTRACT
INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 CONCLUSIONS
 ACKNOWLEDGEMENTS
 REFERENCES
 
We characterize the substate structure of current blockades produced when single-stranded polynucleotide molecules were electrophoretically driven into the {alpha}-hemolysin protein pore. We frequently observe substates where the ionic current is reduced by ~50%. Most of these substates can be associated with a molecular configuration where a polymer occupies only the vestibule region of the pore, though a few appear related to a polymer occupying only the transmembrane β-barrel region of the pore. The duration of the vestibule configuration depends on polymer composition and on which end of the polymer, 3' or 5', subsequently threads into the narrowest constriction and initiates translocation. Below ~140 mV a polymer is more likely to escape from the vestibule against the applied voltage gradient, while at higher voltages a polymer is more likely to follow the voltage gradient by threading through the narrowest constriction and translocating through the pore. Increasing the applied voltage also increases the duration of the vestibule configuration. A semiquantitative model of these trends suggests that escape has stronger voltage dependence than threading, and that threading is sensitive to polymer orientation while escape is not. These results emphasize the utility of {alpha}-hemolysin as a model system to study biologically relevant physical and chemical processes at the single-molecule level.


   INTRODUCTION
TOP
 ABSTRACT
 INTRODUCTION
MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 CONCLUSIONS
 ACKNOWLEDGEMENTS
 REFERENCES
 
Individual single-stranded DNA or RNA molecules can be detected as they are driven electrophoretically through a nanometer-scale pore (1Go,2Go). During translocation, polynucleotide molecules impede the flow of small ions and cause a transient blockade of the ionic current through the pore. These blockade events contain a wealth of information about the composition and concentration of polynucleotide molecules in the vicinity of the pore (3Go–7Go). This has motivated a growing interest in nanopore-based molecular detection technologies (2Go,4Go–6Go,8Go–12Go). The original and most widely used pore for polynucleotide detection is the Staphylococcus aureus {alpha}-hemolysin ({alpha}HL) protein pore (13Go). The {alpha}HL system is a sensitive probe of DNA structure (14Go) and basepairing kinetics (15Go–18Go), DNA-protein interactions (19Go), the physics of confined polymers (20Go), and the physics of ion permeation (21Go). The blockade events produced during the interaction of DNA with {alpha}HL allow the rare combination of high-throughput, real time, single-molecule observations. The composition of the polynucleotide molecules can be easily manipulated, and the structure and composition of pore interior are known (13Go) and can be engineered (10Go,22Go,23Go). Finally, the temperature, applied voltage, and buffer composition can be controlled and characterized with a high degree of precision. These properties make the {alpha}HL system a particularly attractive platform for engineering nanopore-based biosensors and an excellent system for developing and testing theoretical models of nanoscale biological processes (24Go–29Go).

Individual current blockade events observed during the interaction between single-stranded polynucleotide molecules and {alpha}HL often demonstrate multiple substates of distinct conductance (6Go,14Go,30Go). A Deep state where the current is reduced to ~0–30% of the unobstructed Open state is generally associated with translocation of a polynucleotide molecule through the pore (1Go,5Go,9Go,10Go,19Go,30Go,31Go). The dependence of the duration, depth, and frequency of occurrence of the Deep state on polymer length, applied voltage, temperature, and polymer composition have been experimentally characterized (1Go,3Go,5Go,7Go,19Go,31Go). A growing amount of theoretical work has sought to develop a model that can encompass as many of these observations as possible (24Go–28Go,32Go,33Go). A Mid state where the current is reduced to ~30–85% of the Open state is also often observed (6Go,14Go,30Go). Previous studies with double-stranded DNA constructs (9Go,15Go) indicate that a Mid state can be produced when a polynucleotide molecule is resident only in the vestibule portion of the pore. The dynamics of single-stranded polynucleotide molecules in the {alpha}HL vestibule have been investigated in simulation studies (33Go,34Go), but the detailed characteristics of Mid states produced by these polymers have not yet been investigated experimentally. In this article, we analyze the substate structure of current blockade events produced by single-stranded polynucleotide molecules in the {alpha}HL pore. We focus primarily on the vestibule-associated Mid state, exploring how its characteristics depend on polymer composition, polymer orientation, and applied voltage.


   MATERIALS AND METHODS
TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
RESULTS AND DISCUSSION
 CONCLUSIONS
 ACKNOWLEDGEMENTS
 REFERENCES
 
Homogeneous ssDNA oligonucleotides dA50, dC50, and dT50 were synthesized by Integrated DNA Technologies (Coralville, IA) and homogeneous RNA oligonucleotides rA50, rC50, and rU50 were synthesized by Dharmacon RNA Technologies (Lafayette, CO). None of these synthetic oligonucleotide constructs were phosphorylated at either the 5' or 3' end. All samples were PAGE-purified by the suppliers. We resuspended the oligonucleotide samples in the experimental buffer and stored them at –20°C until immediately before use.

Horizontal diphytanoyl-PC/hexadecene bilayers (Avanti Polar Lipids, Alabaster, AL) were formed across a ~20 µm diameter aperture in Teflon. Monomeric {alpha}HL (EMD Biosciences, San Diego, CA) was added to one side of the bilayer (cis side) at a concentration of ~50 µg/ml. The cis side was grounded and positive voltage was applied to the trans side of the bilayer. With 120 mV applied across the bilayer, insertion of a single {alpha}HL pore was marked by an abrupt jump in the current from 0 to ~120 pA. RNA or ssDNA was then added to the cis compartment to a concentration of 1–2 µM. In this configuration, the {alpha}HL pore was oriented such that polynucleotide molecules were driven into the pore through its extramembranous vestibule region. All experiments were performed at 21 ± 2°C in 1 M KCl, 10 mM HEPES/KOH buffered at pH 8.

We used an Axopatch-1B patch-clamp amplifier (Axon Instruments, Union City, CA) to apply voltage across the bilayer and to measure the ionic current flowing through the pore. The analog signal was low-pass filtered at 50 kHz with a four-pole Bessel filter. The amplified, filtered signal was digitized at 250 kHz with a NI PCI-6014 DAQ card (National Instruments, Austin, TX). Data acquisition was controlled with custom software written in LabWindows/CVI (National Instruments). Data analysis was implemented in MatLab (The MathWorks, Natick, MA).


   RESULTS AND DISCUSSION
TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
CONCLUSIONS
 ACKNOWLEDGEMENTS
 REFERENCES
 
Blockades demonstrated rapid transitions between well-defined conductance states
Entry of individual polynucleotide molecules into the {alpha}HL nanopore caused current blockade events with two primary conductance states; a Mid state where the ionic current was between 30% and 85% of the unobstructed Open state, and a Deep state where the current was between 0% and 30% of the Open state (Fig. 1). The rise time of Open state to Mid state or Mid state to Deep state transitions was generally comparable to the 7-µs limit imposed by our 50 kHz antialiasing filter. We thus modeled the structure of each event as a sequence of Mid and Deep conductance states. Using typical threshold-based nanopore event detection algorithms (3Go,4Go,31Go), we defined an Open Threshold at IOS-2{sigma}OS and a Deep Threshold at 0.3*IOS, where IOS is the average Open state current and {sigma}OS is the RMS noise of the Open state (Fig. 1). Mid states were identified as periods when the ionic current was between the Open and Deep Thresholds, and Deep states were identified as periods when the ionic current was below the Deep Threshold. We parameterized each event by the durations and average currents of its constituent substates, and we required that a substate have a duration of 30 µs or longer to be included in our description of the event.


Figure 1
View larger version (11K):
[in this window]
[in a new window]

  		FIGURE 1  Example of ionic-current traces illustrating the substate structure of polynucleotide-induced current blockade events. These traces are from an experiment with dT50 at 120 mV. The horizontal lines show the zero current level and the thresholds used to subdivide events into constituent Mid and Deep states. We defined three kinds of Mid states: M0 states occur in the absence of a Deep state, M1 states directly precede a Deep state, and M2 states directly follow a Deep state. All Deep states were abbreviated as D. We classified each event by its sequence of states. Our names for the various blockade event types appear above the corresponding traces.

 
We focused our analysis on the five types of events shown in Fig. 1. We defined a Mid state that occurred in the absence of a Deep state as an M0 state, a Mid state that directly preceded a Deep state as an M1 state, and a Mid state that directly followed a Deep state as an M2 state. All Deep state signals were defined as D-states. We classified events according to their sequence of states. Table 1 lists the percentages of event types observed for different polymers in experiments with 120 mV applied across the bilayer. M0, M1-D, and D were the most common types of events for all six polymers we investigated. Events with M2 states occurred with small but appreciable frequency only for rA50, dT50, and rU50. Events that were not well described by any of the five event types shown in Fig. 1 were identified by one or more of the following characteristics:
  1. The Mid state had an excursion into the Deep state that was longer than 22 µs but shorter than 30 µs, too long to be attributed to electronic noise but too short for characterization as a Deep state;
  2. The RMS noise of one of the states was larger than twice the RMS noise of the Open state; and
  3. The event demonstrated more than one Deep state signal.


View this table:
[in this window]
[in a new window]

  		TABLE 1  Distributions of event types for the six polymers we investigated

 
Molecular interpretation of blockade substates
The known structure of the {alpha}HL pore (13Go) is central to a physical interpretation of the different types of ionic current blockades. Fig. 2 shows a schematic cross section of an {alpha}HL pore in a lipid bilayer. On the cis side of the bilayer the pore entrance narrows to a diameter of ~2.5 nm before opening into the internal vestibule region. The vestibule region widens to a maximum diameter of ~4.6 nm before narrowing down and joining a short, ~1.5-nm-diameter inner constriction. The inner constriction leads into a ~5-nm-long, ~2-nm-diameter β-barrel that spans the bilayer and opens into the trans compartment. All our experiments were configured such that ssDNA and RNA molecules were driven into {alpha}HL from the cis side of the bilayer.


Figure 2
View larger version (29K):
[in this window]
[in a new window]

  		FIGURE 2  Molecular interpretation of the structure of the blockade events shown in Fig. 1. We hypothesize that the vestibule configuration produces both M0 and M1 states and that the threaded configuration produces Deep states. We also suggest that either retraction from the threaded configuration back into the vestibule configuration or pausing of a polymer in the exit configuration could produce M2 states.

 
Knowledge of the structure of {alpha}HL combined with a crude volume-exclusion model for the obstruction of current through {alpha}HL leads to a simple physical interpretation (Fig. 2) of the various event types shown in Fig. 1 (35Go). When a polymer occupies only the vestibule region of the pore, it reduces the current to the Mid state. This vestibule configuration leads to M0 events if the polymer escapes back into the cis compartment. Alternatively, one end of the polymer may thread through the inner constriction of the pore and into the β-barrel, impeding the passage of ions enough to reduce the current to the Deep state. If the polymer continues threading through the β-barrel, it will eventually exit on the trans side of the bilayer. This translocation scenario accounts for M1-D type events if the vestibule configuration persists long enough to produce a resolvable Mid state (longer than our 30-µs temporal resolution). We interpret D-type events as translocations where the vestibule configuration persisted for <30 µs. It is also possible, in principle, for a polymer that has threaded into the inner constriction to retract back into the vestibule configuration. In the last two traces of Fig. 1, a well-resolved Mid state is observed after the Deep state. It is clear how retraction from the threaded configuration back into the vestibule configuration followed by escape into the cis compartment could account for these M2 signals. However, if the Deep state corresponds to translocation, then the subsequent M2 state could result if the polymer's progress through the pore is arrested in a configuration where its trailing end occupies part of the β-barrel, but is no longer threaded through the inner constriction. We refer to this hypothetical configuration as the exit-configuration.

The association of the vestibule configuration with reduction of the current to the Mid state and the threaded configuration with further reduction of the current to the Deep state is supported by results from previous investigations. Howorka et al. investigated the binding of free oligonucleotides to complementary strands that were tethered to the interior of the {alpha}HL vestibule (36Go). When the transient, double-stranded DNA constructs had sufficiently long single-stranded overhangs, they produced current blockades that fluctuated between Mid and Deep conductance states. Complexes with shorter overhangs produced blockades that only demonstrated Mid states. Similar data and interpretations have been given for experiments with hairpin DNA having variable length duplex regions (9Go,37Go). These results strongly support the association of the Mid state with the vestibule configuration and the association of the Deep state with the threaded configuration. Further, a recent simulation study that simultaneously computed the time evolution of a single-stranded polymer interacting with the {alpha}HL pore and the corresponding modulation of the ionic current was able to reproduce this relationship between molecular configurations and degrees of current obstruction (33Go).

Homogeneous single-stranded DNA and RNA polymers are highly flexible, so it is more difficult to infer their particular arrangement in the {alpha}HL vestibule than for the more rigid double-stranded DNA constructs. We suggest two distinct ways in which a polymer may enter the vestibule from the cis compartment. One possibility is that entrance is initiated by fluctuation of an end segment into the vestibule, and the other possibility is that an internal segment of the polymer forms a transient loop that fluctuates into the vestibule and initiates polymer entrance. The result of Vercoutere et al. (9Go), where ssDNA constructs with hairpins at both the 3' and 5' ends did not cause resolvable Mid state blockades, argues against loop-initiated entrance. We also note two potential arrangements of the polymer in the vestibule configuration. Either the entire 50-mer is drawn inside the vestibule, or only part of the 50-mer is in the vestibule while the rest of the polymer remains in the cis compartment. The latter possibility is illustrated in Fig. 2. In the duplex DNA experiments, the {alpha}HL vestibule accommodated up to ~24 nucleotides and still demonstrated Mid states (36Go), suggesting that at least half of a 50-nucleotide-long ssDNA or RNA molecule could be drawn into the vestibule. Jung et al. estimated that the accessible volume of the {alpha}HL vestibule, corrected for hydration of the vestibule surface, is ~32,600 Å3 (38Go). Deamer et al. estimated a molecular volume of ~290 Å3 for a single Adenylic nucleotide (39Go). Assuming that each nucleotide is accompanied by four water molecules leads to a rough estimate of ~17,000 Å3 for the volume occupied by a hydrated dA50 molecule. These calculations indicate that the vestibule is, in principle, large enough to hold an entire ssDNA or RNA 50-mer. While it is energetically favorable for the entire polymer to enter the vestibule, the entropic contribution to the free energy is quite complex. As the polymer enters the vestibule there are entropic costs associated with removing each monomer from the bulk cis volume (26Go) and from packing additional monomers into the confining vestibule (40Go). However, there is also entropic benefit associated with moving the polymer from a partially filled state, where a segment is restricted to the entrance, to a completely filled state where the polymer is free to reorient within the confining vestibule (20Go). It is not clear, a priori, which of these effects dominate in our experiments. We will return to the question of how much of the polymer is drawn into the vestibule after presenting our data.

Deep state statistics
Each point in Fig. 3 gives the duration and average current of one Deep state. The panels show representative data from one experiment at 120 mV for each of the six polymers we investigated. We display the data on a log-linear scale to allow visualization of the Deep state duration distributions across many orders of magnitude. Clear grouping behavior is observed in all of the panels. Most of these groups correspond to local maxima in the two-dimensional distributions of Deep state durations and currents. However, the groups of long duration states in the dT50 and rU50 panels do not correspond to local maxima, but appear as groups because of the logarithmic scaling of the duration axis. Previous work strongly suggests that the signals in the maxima-associated groups arise from translocation of polymers through the pore (1Go,14Go,31Go). The general characteristics we observed for these maxima-associated signals are in good agreement with previously reported results (1Go,3Go,14Go,19Go,30Go,31Go). These references present extensive characterization and discussion of the Deep state. In this report, we focus only on the aspects of the Deep state that are important to our characterization and interpretation of blockade structure.


Figure 3
View larger version (18K):
[in this window]
[in a new window]

  		FIGURE 3  Distributions of Deep state durations and currents. Each point corresponds to the duration and normalized average current of one Deep state. Data were acquired with 120 mV applied across the bilayer. We display the data on a semilogarithmic scale to allow visualization of the Deep state distributions across many orders of magnitude in duration. The vertical lines give the locations of minima in the duration distributions (see Fig. 5). We only associate signals longer than these minima with translocation of polymers through the pore. The longest duration groups of signals in the dT50 and rU50 panels do not correspond to local maxima in the distributions. These long signals are well described by exponential distributions with time constants of ~3.7 ms for dT50 (Fig. 4) and ~94 ms for rU50. They appear as groups because of the logarithmic scaling of the duration axis. Previous work found that the two groups in the dA50 and rC50 distributions were related to 3'- or 5'-first translocation.

 

Figure 5
View larger version (12K):
[in this window]
[in a new window]

  		FIGURE 5  rC50 Deep state duration distribution. This distribution shows the same data as the rC50 panel of Fig. 3, but the signals have been binned in duration and are displayed on a linear scale. There is a clear minimum in this duration distribution between our 30-µs experimental cutoff time and the translocation-associated maxima. We observed similar minima in the Deep state duration distributions of other polymers we investigated (Fig. 3, solid vertical lines). We designated signals longer than the minima as translocation-associated Deep states and those shorter than the minima as ambiguous Deep states.

 

Figure 4
View larger version (7K):
[in this window]
[in a new window]

  		FIGURE 4  Duration distribution of the extremely long dT50 Deep states. These signals appear as the rightmost group in the dT50 panel of Fig. 3. As is clear from this plot, the group of long dT50 Deep states does not correspond to an additional maximum in the dT50 Deep state duration distribution. The long signals constitute ~20% of the total number of signals in the dT50 panel of Fig. 3 and are well described by an exponential distribution with a time constant of ~3.7 ms (dashed curve).

 
The durations of the Deep states in the translocation-associated groups for the five homopolymers other than dC50 were appreciably longer than our 30-µs cutoff, suggesting that we were able to detect nearly all translocation events in experiments with these polymers at 120 mV. It is clear from the dC50 panel that many translocation-associated signals were shorter than our 30-µs cutoff time. Events demonstrating these short translocation signals were classified as Other, and we hypothesize that fast translocation accounts for a significant fraction of dC50 Other-type events listed in Table 1. The two groups in the dA50 and rC50 distributions have previously been shown to arise from translocation of polymers in either of two possible orientations, 3'-first or 5'-first (28Go,30Go). We will utilize this orientation information in our subsequent analysis of the Mid state.

The long-duration Deep states in the dT50 and rU50 panels have not been previously reported. For dT50, Deep states with duration longer than 4 ms constitute ~20% of the total number of signals shown in the panel and their distribution is well described by a single exponential with a ~3.7-ms time constant (Fig. 4). Analogously for rU50, Deep states with duration longer than 10 ms constitute ~21% of the total and their distribution is well described by a single exponential with a ~94-ms time constant. Previous analysis of Deep state data found a 10–100-µs exponential time constant for the fall-off of the Deep state duration distribution in the vicinity of the maxima (3Go). However, this analysis did not take into account these extremely long events. The characteristics of the Mid states before and after the long Deep states are similar to those of the translocation-associated signals, suggesting that the long signals also represent translocation. To explain these long Deep states, we hypothesize that during translocation a polymer's forward progress can be significantly inhibited or even stopped entirely for relatively long periods of time. Tangling of the polymer at the neck or binding of one or more bases to the residues lining the β-barrel are two possible causes of such a stoppage. The Deep state would then persist until untangling or unbinding occurred, allowing normal translocation to resume and the subsequent exit of the polymer from the pore into the trans compartment.

A final important feature in the Deep state distributions is that, with the exception of dC50, we observed minima in the duration distributions between our 30 µs experimental cutoff time and the first maxima (Fig. 5). These minima are not readily apparent in Fig. 3 due to the logarithmic scaling, but their locations are given by the solid vertical lines. The rise in the distribution for durations shorter than these minima is not consistent with a simple picture of translocation, wherein one expects a monotonic decrease in probability for durations shorter than the most probable duration (24Go). Similar minima in Deep state duration distributions have been previously reported (1Go,30Go). Signals shorter than the minima were attributed to multiple possible effects such as retraction from the threaded configuration back into the vestibule configuration, very rapid translocation, or translocation of short polynucleotide fragments. Due to the ambiguity in the interpretation of these short Deep states, we only designated Deep states with durations longer than the minima as translocations (Fig. 5). If some of these short signals represent retraction, then we would expect to find events with multiple Deep states. Such events are often observed in experiments with short double-stranded DNA constructs in the {alpha}HL vestibule (36Go,37Go,41Go). We occasionally observed such events in experiments with dT50 and rA50, but their frequency of occurrence was difficult to quantify because they were rare and difficult to identify with our simple data analysis algorithm. Multiple Deep state events in experiments with the other four polymers were almost never observed. The existence of events with multiple Deep states suggests that retraction from the threaded configuration back into the vestibule configuration can occur and be resolved in our experiments. However, this process seems quite rare and it remains unclear to what extent it can account for the short, ambiguous Deep states.

Mid state statistics
Fig. 6 shows the durations and currents of the three types of Mid states, M0, M1, and M2, from experiments at 120 mV with dT50, rA50, and rU50. Each distribution displays a characteristic current, IM, around which the majority of the signals are clustered. We quantified this characteristic current with a Gaussian fit to the largest peak in a weighted distribution of the currents of the Mid states. The duration of each signal was used as the weighting factor, and the longest 5% of signals in each data set were excluded so that the peak location was not dominated by a minority of extremely long signals. IM values for the distributions in Fig. 6 are shown as solid horizontal lines. We focused our analysis on Mid states of similar current by defining a primary group with boundaries at IM ± (w + 2*{sigma}OS/sqrt(nobs)). The term w is an empirical width factor that we set to 7.5 pA. The term 2*{sigma}OS/sqrt(nobs) accounts for the expected increased statistical variation in the estimated average current of shorter signals, where {sigma}OS is the RMS noise observed in the Open state and nobs is the number of samples of the ionic current time series used to estimate the average current of a given Mid state. The dashed curves in Fig. 6 show the boundaries defining the primary groups of signals. We used the average duration of the Mid states in the primary group, tM, to parameterize the Mid state duration distributions. Again, the longest 5% of Mid states within the primary group were excluded from the calculation of tM to prevent the statistic from being dominated by a minority population of very long signals. The tM values for each data set are shown as solid vertical lines in each panel of Fig. 6. M0 and M1 data for dA50, dC50, and rC50 followed the same general patterns as the data shown in Fig. 6, while the occurrence of M2 signals was negligible in experiments with these polymers. Table 2 summarizes the average values of IM and tM at 120 mV for all six of the polymers we investigated.


Figure 6
View larger version (20K):
[in this window]
[in a new window]

  		FIGURE 6  Distributions of Mid state durations and currents. Each point gives the duration and normalized average current of one Mid state. Rows show representative data from one experiment at 120 mV with the indicated polymer. Columns correspond to different types of Mid states (see Fig. 1). The solid horizontal lines show the characteristic Mid state current, IM, and the solid vertical lines show tM, the average duration of the Mid state signals contained within the dashed curves. We associate M0 and M1 states with the presence of a polymer in the {alpha}HL vestibule. Arrows in the M2 panels point out signals in the M2 distributions that have different normalized currents than the M0 or M1 signals. We hypothesize that these differing M2 states are produced when a polymer is held in a configuration where it still occupies part of the β-barrel but is no longer threaded through the neck of the pore (Fig. 2). Table 2 summarizes IM and tM values for the six polymers we investigated.

 

View this table:
[in this window]
[in a new window]

  		TABLE 2  Summary of Mid state characteristics for different polymers with 120 mV applied across the bilayer

 
In comparing the first two columns of Fig. 6, one sees that for a given polymer there is a strong similarity between the M0 and M1 distributions. This similarity is also reflected in the IM and tM values in Table 2 (except for the dT50 tM values). These observations are consistent with the association of both M0 and M1 signals with a polymer in the vestibule configuration. Since we interpret M0 signals as escape from the vestibule configuration back into the cis compartment, we define an escape probability as N(M0)/(N(M0)+N(M1)), where N(M0) and N(M1) are the number of signals in the primary groups of the M0 and M1 distributions, respectively. Average escape probabilities for experiments at 120 mV are listed in Table 2 and range between ~0.5 and ~0.7. This indicates that at 120 mV, a homogeneous, 50-nucleotide long DNA or RNA molecule residing in the vestibule configuration is slightly more likely to escape against the voltage gradient back into the cis compartment than to follow the voltage gradient and translocate through the pore.

The majority of the M2 signals demonstrate current levels that are clearly distinct from M0 and M1 current levels (Fig. 6, arrows). Assuming that the Deep states preceding these M2 signals indicate translocation, we hypothesize that these Mid states occur when a polymer's exit into the trans compartment is hindered while it still occupies part of the β-barrel but is no longer threaded through the inner constriction of the pore (Fig. 2). The mechanism that hinders polymer exit could also occur at a previous stage in the translocation process, when the polymer was still threaded through the inner constriction of the pore and blocking the current at the Deep state. In this case, one might expect that polymers demonstrating longer, more frequent M2 signals would also demonstrate longer Deep states. This correlation is found in our experiments as dT50 and rU50 produce the longest Deep states and the most frequent M2 states. The M2 distributions in Fig. 6 also show a few signals with current levels similar to the characteristic current levels in the M0 and M1 distributions. These M2 signals are consistent with the hypothesis that a polymer can retract from the threaded configuration back into the vestibule configuration. We interpret these signals in the same way we interpret the existence and rarity of events with multiple Deep states—retraction of a threaded polymer back into the vestibule configuration can be observed, but only very rarely.

For the six homopolymers we investigated, the Mid state characteristics are moderately sensitive to polymer composition. M0 and M1 states for dA50, dT50, rC50, and rU50 all show a normalized residual current of IM/IOS {approx} 0.59, whereas dC50 signals show a slightly larger residual current and rA50 signals show a significantly smaller residual current. The characteristic durations of M0 and M1 states for rA50 and dT50 are significantly longer than those of the other four polymers. The largest escape probability was 0.69 for rU50, and the smallest was 0.48 for dT50, with the escape probabilities of the remaining polymers more or less evenly distributed between the two extreme values. These polymer-specific differences in Mid state characteristics could enhance the ability of {alpha}HL to analyze and distinguish polynucleotide samples. However, further work is needed to understand the molecular origins of these differences. For example, the large difference between the Mid state characteristics of dA50 and rA50 suggests that specific interactions between Adenine bases and the residues lining the vestibule are not the primary factors determining the Mid state characteristics of a given polymer.

Voltage dependence of the vestibule configuration
We conducted experiments at a variety of voltages over the range of 100–200 mV. Data obtained at each voltage were analyzed as described in the previous section. Fig. 7 shows the voltage dependencies of IM/IOS, tM, and the escape probability for dA50, rC50, and rU50. The IM/IOS and tM values in Fig. 7 were derived from the M1 distributions as opposed to the M0 distributions because we observed a large number of M1 signals for all polymers at all voltages. As we increased the voltage we saw a moderate but clear decrease in the normalized M1 current level (Fig. 7 a), a substantial increase in the typical M1 duration (Fig. 7 b), and a substantial decrease in the escape probability to near zero at 200 mV (Fig. 7 c). Voltage curves for dC50, dT50, and rA50 are not shown, though their trends are qualitatively the same as those shown in Fig. 7.


Figure 7
View larger version (6K):
[in this window]
[in a new window]

  		FIGURE 7  Voltage dependence of the vestibule configuration for dA50, rC50, and rU50. (a) Normalized M1 current. (b) Mean duration of the primary group of M1 states. (c) Escape probability. The escape probability reflects the likelihood that a polymer which has inserted into the vestibule will escape back into the electrically negative compartment. It was calculated as the ratio of the number of M0 states to the sum of the number of M0 and M1 states. Voltage-dependent data for dC50, dT50, and rA50 shows the same trends as the data in this figure, but we do not display it here.

 
The general voltage-dependent trends shown in Fig. 7 provide further evidence supporting the hypothesis that M0 and M1 signals correspond to a polymer in the vestibule configuration. If a higher applied positive voltage on the trans side pulls more negatively charged bases from the cis side into the vestibule, then the volume exclusion mechanism of current obstruction accounts for the decrease in the normalized Mid state current (Fig. 7 a). This simple mechanism can explain the general trend in the data, but a variety of additional factors may also contribute to the observed Mid state current (42Go). The trends in Fig. 7, b and c, can be explained by a simple kinetic model wherein escape from the vestibule and threading of a polymer end from the vestibule into the inner constriction would be described as activated processes with associated rate-constants β and {gamma}, respectively. If the vestibule configuration obeyed such simple kinetics, its duration would be exponentially distributed with a time constant of {tau}M = (β + {gamma})–1, and the escape probability would be given by β/(β + {gamma}). The trends in tM and the escape probability shown in Fig. 7 could then result if 1), escape is favored at low voltages, while threading is favored at high voltages (β > {gamma} below ~140 mV, {gamma} > β above ~140 mV); 2), the escape rate constant, β, decreases with increasing voltage; and 3), escape has a stronger voltage dependence than threading (over the range 100–200 mV, |{partial}β/{partial}V| > |{partial}{gamma} /{partial}V|).

Polymer orientation and the vestibule configuration
Before making a more detailed comparison between our data and this kinetic model of the vestibule configuration, we must consider the effects of polymer orientation. We previously demonstrated that translocation of an rC50 polymer in a 3'->5' orientation produced stronger obstruction of the ionic current than translocation in a 5'->3' orientation (30Go), and Mathe et al. found an analogous effect for poly-dA translocation (28Go). This orientation-dependent Deep state conductance leads to the two well-resolved groups of signals in the rC50 and dA50 panels of Fig. 3. We utilized this effect to distinguish M1 states that were followed by 3'-first translocation from those that were followed by 5'-first translocation. The resulting M1 duration distributions are shown in Fig. 8. These distributions indicate that 5'-associated M1 states tend to be longer than 3' states, and that this difference is enhanced as the voltage increases. Application of the kinetic model we have proposed requires that the duration distributions of 3' and 5'-associated vestibule configurations be individually well described by single exponentials. Single-exponential fits accurately represent the dA50 data for both orientations and the rC50 data for 5'-associated M1 states. These fits were determined using the maximum likelihood method. M1 states followed by 3'-first translocation for rC50 were rare and short, making it difficult to assess the exponential character of these data. The existence of orientation-dependent differences in the duration of the vestibule configuration for the other four polymers is suggested by the observation that, at voltages above ~140 mV, their overall M1 duration distributions are better described using two exponentials rather than a single exponential.


Figure 8
View larger version (17K):
[in this window]
[in a new window]

  		FIGURE 8  M1 state duration distributions grouped by polymer orientation. Panels a and b show data for rC50 and panels c and d show data for dA50. Data in top row were obtained at 120 mV while the data in the bottom row were obtained at 160 mV. The correlation between Deep state characteristics and polymer orientation shown in Fig. 3 enabled us to group these M1 states by orientation.

 
The correlation between translocation orientation and vestibule configuration lifetime suggests that there are two distinct vestibule configurations, one for each orientation of the polymer in the vestibule. This could occur if only part of the polymer resides in the vestibule, with one end remaining in the cis compartment. Alternatively, if the entire polymer resides in the vestibule then its conformational freedom must be so restricted that it cannot reorient itself within the vestibule. While it remains unclear how much of the polymer resides in the vestibule, it nonetheless seems appropriate to use separate sets of rate constants to model 3'- and 5'-associated vestibule configuration kinetics.

Kinetic model of vestibule configuration
We are able to compare our dA50 data to a kinetic model with the characteristics described above because its 3'- and 5'-associated M1 states can be distinguished and the orientation-dependent M1 duration distributions are individually well described by single exponentials (Fig. 8). However, because we cannot distinguish between 3'- and 5'-associated M0 signals, we are unable to estimate orientation-specific 3'- and 5'-associated escape probabilities and use them to derive estimates for the escape and threading rate constants directly from our data. We can estimate the overall escape probability, which is a weighted average of the orientation-specific escape probabilities. It is then possible to make a semiquantitative analysis of the vestibule configuration kinetics by assuming a particular model for the voltage dependence of the rate constants and then assessing that model's ability to reproduce our voltage-dependent {tau}3, {tau}5, and escape probability data (Fig. 9). We assume a very simple model where the 3'- and 5'-associated escape rate constants, β3 and β5, are identical and decrease exponentially with voltage (Eq. 1), while the threading rate constants, {gamma}3 and {gamma}5, are different but are both independent of voltage (Eqs. 2 and 3):

Formula 1(1)

Formula 2(2)

Formula 3(3)
In Eq. 1, b is the escape rate constant with zero applied voltage, e is the elementary charge, z is the effective charge of a phosphate on the polynucleotide backbone after taking into account the effects of counterions, n is the number of phosphates resident in the vestibule during the transition state for escape, {delta} is the fractional voltage drop in the vestibule, V is the applied voltage, k is the Boltzmann constant, and T is the temperature. The values g3 and g5 are the voltage-independent rate constants for threading of a polymer-end through the neck of the pore. This model is a simplified implementation of the Woodhull formalism (43Go). A similar approach has previously been used to describe the kinetics of polypeptides in the {alpha}HL pore (44Go). To simplify our notation we define a characteristic escape voltage VE {equiv} [ezn{delta}/kT]–1. This leads to the following expressions for the exponential time constants and the overall escape probability:

Formula 4(4)

Formula 5(5)

Formula 6(6)
G3 and G5 are the observed numbers of 3'- and 5'-associated M1 signals, and B is the total number of observed M0 signals. The parameter values g3 = 0.016 µs–1, g5 = 0.006 µs–1, b = 18 µs–1, and VE = 17 mV give reasonable agreement between this model and our data (Fig. 9). These values were found by minimization of the total sum-squared deviation between the data in Fig. 9 and model of Eqs. 4–6GoGo. Weighting factors were set such that the contributions to the total sum-squared deviation were approximately equal for the three curves.


Figure 9
View larger version (8K):
[in this window]
[in a new window]

  		FIGURE 9  Voltage-dependent dA50 data compared to a kinetic model with voltage-dependent rate constants (Eqs. 4–6GoGo). M1 lifetimes at each voltage were derived from exponential fits to the 3' and 5' associated M1 duration distributions. The escape probability data in the bottom panel ({diamond}) is the same as the dA50 escape probability data shown in Fig. 7 c, while the model values (x) were derived from Eq. 6. The model for the escape probability cannot be shown as a continuous curve because it requires experimental data (the observed numbers of 3'- and 5'-associated M1 signals). The model curves and values were generated with the parameters g3 = 0.015 µs–1, g5 = 0.006 µs–1, b = 18 µs–1, and VE =17 mV in Eqs. 4–6GoGo.

 
Introducing voltage dependence into the threading rate constants as {gamma} = gexp(V/VTh) required very weak voltage dependence (VTh > 1000 mV) to give a reasonable description of our data. We found that a variety of other five-parameter models were also able to give reasonable descriptions of our data, but none did appreciably better than the four-parameter model of Eqs. 1–3GoGo. We also explored additional four-parameter models, one where the orientation dependence arose from differing voltage-independent factors in the escape rate constants (β3 = b3exp(–V/VE), β5 = b5exp(–V/VE), and {gamma}3 = {gamma}5 = g), and another where the orientation dependence arose from differing 3'- and 5'-associated escape voltages (β3 = bexp(–V/VE3), β5 = bexp(–V/VE5), and {gamma}3 = {gamma}5 = g). Neither of these four-parameter models gave acceptable descriptions of our data.

Interpretation of vestibule configuration model
The model specified by Eqs. 1–3GoGo strongly supports the assertion that over the range V ~ 100–200 mV, the escape process is much more sensitive to voltage than the threading process. The relatively weak voltage dependence of threading may be related to a small voltage gradient within the vestibule (36Go) and the lack of negatively charged terminal phosphate groups at both the 3' and 5' ends of our synthetic polymers. The arrival of a polymer end at the neck may therefore be a diffusive process which is little biased by the applied voltage. Furthermore, insertion of the uncharged sugar/base portion of the terminal nucleotide far enough into the pore neck to allow the charged phosphate to experience the high electric field could also be a relatively voltage-insensitive step in the threading process. Another aspect of the vestibule configuration suggested by our model is that the orientation dependence of the vestibule configuration lifetime (Fig. 8) can be primarily attributed to different rate constants for 3' vs. 5' threading. Specifically, 3' threading is easier than 5' threading. This finding is consistent with the results of Mathe et al. (28Go), where experiments showed that poly-dA segments passed more rapidly through the {alpha}HL inner constriction when they moved in the 3'->5' direction as compared to the 5'->3' direction. In the same report, molecular dynamics simulations of a polydA molecule threaded through the {alpha}HL pore also found that 3'->5' passage was faster than 5'->3' passage. These simulations revealed a microscopic mechanism for the orientation asymmetry. When the polymer was confined in a tube with diameter <1.5 nm, the DNA bases tilted toward the 5' end of the polymer. Because initial threading of a polymer and passage of a threaded polymer through the neck are very similar processes, we expect that the confinement-induced tilt is what makes 3' threading easier than 5' threading. Because the pore entrance that connects the vestibule to the cis compartment is ~2.5 nm in diameter, we do not expect an asymmetric tilting of the bases to be important to the escape process. Our finding that the voltage-independent factor in the escape rate constant (b in Eq. 1) need not depend on orientation is consistent with this expectation. The orientation independence of escape also suggests that z, n, and {delta} in Eq. 1, the microscopic factors that determine the voltage sensitivity of escape, do not depend strongly on polymer orientation. Finally, we note that assumptions z = 0.3, n = 20, and {delta} = 0.2 lead to the estimate VE = 20 mV, indicating that our crude experimental value of 17 mV is physically plausible. The kinetic model we have presented is quite simple, but it gives a very reasonable description of our data and it suggests a number of interesting properties of the behavior of a polynucleotide molecule when it is confined in the {alpha}HL vestibule.


   CONCLUSIONS
TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 CONCLUSIONS
ACKNOWLEDGEMENTS
 REFERENCES
 
We presented an analysis of the current blockades produced by single-stranded polynucleotide molecules in the {alpha}HL pore that takes into account their substate structure. We modeled the structure of each blockade as a sequence of Mid (current at ~50% of unobstructed level) and Deep (current <~30% of unobstructed level) conductance states. There appear to be two distinct configurations of the polymer in the pore that cause a Mid state. Our data are consistent with the interpretation that the most common Mid state occurs at the beginning of the polymer/pore interaction, when part of a polymer enters the {alpha}HL vestibule from the cis compartment. The ionic current remains in the Mid state until the polymer either 1), escapes back into the cis compartment, causing a return to the Open state current; or 2), threads through the inner constriction of the pore and blocks the current at the Deep state. For some polymers, up to ~80% of the observed current blockade events demonstrated resolvable signals indicative of a polymer in the vestibule configuration. Occasionally we observed Mid states that followed Deep states and demonstrated a different amount of residual current than the vestibule-associated Mid states. We attribute these Mid states to a configuration where the trailing end of the polymer occupies part of the β-barrel region of the pore, but is no longer threaded through the pore's inner constriction. These β-barrel signals were much less frequent than the vestibule-associated signals, occurring in between 3% and 8% of the blockade events for three of the polymers we investigated (rA50, rU50, and dT50), and almost never occurring for the other three polymers (dA50, dC50, and rC50). We observed clear statistical differences in the Mid state characteristics of different polymers, indicating that Mid state analysis could enhance the capability of {alpha}HL as a single-molecule analytical device. Simulation studies of these polymer-specific differences could provide much needed insight into their molecular origins, and accurate prediction of these differences would provide very strong validation of the simulation methodology.

We focused particular attention on the signals indicative of the vestibule configuration because they are a very prominent but little studied component of the current blockades produced by polynucleotide molecules in the {alpha}HL pore. To gain insight into the general physics of the vestibule configuration, we analyzed vestibule configuration data at a variety of voltages between 100 and 200 mV. We observed a clear decrease in the normalized current with increasing voltage. This may be the result of crowding as more bases are pulled into in the vestibule at higher voltages. We also found that the characteristic duration of the vestibule configuration increased with increasing voltage while the probability of polymer-escape from the vestibule back into the electrically negative volume decreased. A simple kinetic picture where the rate constant for escape decreases strongly with voltage while the rate constant for threading has much weaker voltage dependence qualitatively explains these trends. This simple picture was slightly complicated by the observation that the duration of the vestibule configuration depended on which end of the polymer, 3' or 5', subsequently threaded into the inner constriction. We found that a kinetic model with the following characteristics gave a very nice description of our orientation and voltage-dependent dA50 data: 1), the escape rate constant decreases exponentially with voltage but does not depend on orientation; and 2), 3'- and 5'-threading have different rate constants that do not depend on the applied voltage. From this semiquantitative analysis we conclude that, for our unphosphorylated ssDNA and RNA 50-mers, threading has much weaker voltage dependence than escape, threading is sensitive to polymer orientation, and escape is insensitive to polymer orientation. The intricate molecular details revealed by this analysis make an important contribution to the overall understanding of the interaction between single-stranded polynucleotide molecules and the {alpha}HL pore. Further, the vestibule configuration represents an additional aspect of the {alpha}HL system that can be used to study the physical characteristics of nucleic acids. While our data support the conclusion that there are two distinct vestibule configurations, one for each orientation of the polymer, we are presently unable to make any more detailed inferences into the arrangement of the polymer in the vestibule. Experiments with polymers of various lengths will likely give insight into this question, possibly allowing identification of a characteristic polymer length below which the entire polymer is drawn into the vestibule and above which part of the polymer remains in the cis compartment. This information may give new insight into the flexibility and packaging properties of single-stranded polynucleotide molecules on nanometer length scales.


   ACKNOWLEDGEMENTS
TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 CONCLUSIONS
 ACKNOWLEDGEMENTS
REFERENCES
 
We are very grateful to Bertil Hille for loan of the patch-clamp amplifier used in this work, and we thank Bertil Hille and Michael Schurr for careful review of this manuscript.

This work was supported by the National Science Foundation Integrative Graduate Research and Education Traineeship, The University of Washington Royalty Research Fund, The University of Washington Physics Department, and the University of Washington Microbiology Department.


   FOOTNOTES
 
M. Troll's current address is Aaron Thermal Technologies, Seattle, WA 98115.

Editor: David S. Weiss.

Submitted on February 16, 2007; accepted for publication June 25, 2007.


   REFERENCES
TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 CONCLUSIONS
 ACKNOWLEDGEMENTS
 REFERENCES
 
1. Kasianowicz, J. J., E. Brandin, D. Branton, and D. W. Deamer. 1996. Characterization of individual polynucleotide molecules using a membrane channel. Proc. Natl. Acad. Sci. USA. 93:13770–13773.[Abstract/Free Full Text]

2. Li, J. L., M. Gershow, D. Stein, E. Brandin, and J. A. Golovchenko. 2003. DNA molecules and configurations in a solid-state nanopore microscope. Nat. Mater. 2:611–615.[CrossRef][Medline]

3. Meller, A., L. Nivon, E. Brandin, J. Golovchenko, and D. Branton. 2000. Rapid nanopore discrimination between single polynucleotide molecules. Proc. Natl. Acad. Sci. USA. 97:1079–1084.[Abstract/Free Full Text]

4. Wang, H., J. E. Dunning, A. P. H. Huang, J. A. Nyamwanda, and D. Branton. 2004. DNA heterogeneity and phosphorylation unveiled by single-molecule electrophoresis. Proc. Natl. Acad. Sci. USA. 101:13472–13477.[Abstract/Free Full Text]

5. Nakane, J., M. Akeson, and A. Marziali. 2002. Evaluation of nanopores as candidates for electronic analyte detection. Electrophoresis. 23:2592–2601.[CrossRef][Medline]

6. Kasianowicz, J. J., S. E. Henrickson, H. H. Weetall, and B. Robertson. 2001. Simultaneous multianalyte detection with a nanometer-scale pore. Anal. Chem. 73:2268–2272.[Medline]

7. Henrickson, S. E., M. Misakian, B. Robertson, and J. J. Kasianowicz. 2000. Driven DNA transport into an asymmetric nanometer-scale pore. Phys. Rev. Lett. 85:3057–3060.[CrossRef][Medline]

8. Bayley, H., and P. S. Cremer. 2001. Stochastic sensors inspired by biology. Nature. 413:226–230.[CrossRef][Medline]

9. Vercoutere, W., S. Winters-Hilt, H. Olsen, D. Deamer, D. Haussler, and M. Akeson. 2001. Rapid discrimination among individual DNA hairpin molecules at single-nucleotide resolution using an ion channel. Nat. Biotechnol. 19:248–252.[CrossRef][Medline]

10. Howorka, S., S. Cheley, and H. Bayley. 2001. Sequence-specific detection of individual DNA strands using engineered nanopores. Nat. Biotechnol. 19:636–639.[CrossRef][Medline]

11. Nakane, J., M. Wiggin, and A. Marziali. 2004. A nanosensor for transmembrane capture and identification of single nucleic acid molecules. Biophys. J. 87:615–621.[Abstract/Free Full Text]

12. Storm, A. J., J. H. Chen, H. W. Zandbergen, and C. Dekker. 2005. Translocation of double-strand DNA through a silicon oxide nanopore. Phys. Rev. E. 71:051903.[CrossRef]

13. Song, L., M. R. Hobaugh, C. Shustak, S. Cheley, H. Bayley, and J. E. Gouaux. 1996. Structure of staphylococcal alpha-hemolysin, a heptameric transmembrane pore. Science. 274:1859–1866.[Abstract/Free Full Text]

14. Akeson, M., D. Branton, J. J. Kasianowicz, E. Brandin, and D. W. Deamer. 1999. Microsecond timescale discrimination among polycytidylic acid, polyadenylic acid, and polyuridylic acid as homopolymers or as segments within single RNA molecules. Biophys. J. 77:3227–3233.[Abstract/Free Full Text]

15. Howorka, S., L. Movileanu, O. Braha, and H. Bayley. 2001. Kinetics of duplex formation for individual DNA strands within a single protein nanopore. Proc. Natl. Acad. Sci. USA. 98:12996–13001.[Abstract/Free Full Text]

16. DeGuzman, V. S., C. C. Lee, D. W. Deamer, and W. A. Vercoutere. 2006. Sequence-dependent gating of an ion channel by DNA hairpin molecules. Nucleic Acids Res. 34:6425–6437.[Abstract/Free Full Text]

17. Mathe, J., A. Arinstein, Y. Rabin, and A. Meller. 2006. Equilibrium and irreversible unzipping of DNA in a nanopore. Europhys. Lett. 73:128–134.[CrossRef]

18. Sauer-Budge, A. F., J. A. Nyamwanda, D. K. Lubensky, and D. Branton. 2003. Unzipping kinetics of double-stranded DNA in a nanopore. Phys. Rev. Lett. 90:238101.[CrossRef][Medline]

19. Meller, A., and D. Branton. 2002. Single molecule measurements of DNA transport through a nanopore. Electrophoresis. 23:2583–2591.[CrossRef][Medline]

20. Muthukumar, M. 2001. Translocation of a confined polymer through a hole. Phys. Rev. Lett. 86:3188–3191.[CrossRef][Medline]

21. Bonthuis, D. J., J. S. Zhang, B. Hornblower, J. Mathe, B. I. Shklovskii, and A. Meller. 2006. Self-energy-limited ion transport in subnanometer channels. Phys. Rev. Lett. 97:128104.[CrossRef][Medline]

22. Bayley, H., and L. Jayasinghe. 2004. Functional engineered channels and pores (review). Mol. Membr. Biol. 21:209–220.[CrossRef][Medline]

23. Astier, Y., O. Braha, and H. Bayley. 2006. Toward single molecule DNA sequencing: direct identification of ribonucleoside and deoxyribonucleoside 5'-monophosphates by using an engineered protein nanopore equipped with a molecular adapter. J. Am. Chem. Soc. 128:1705–1710.[CrossRef][Medline]

24. Lubensky, D. K., and D. R. Nelson. 1999. Driven polymer translocation through a narrow pore. Biophys. J. 77:1824–1838.[Abstract/Free Full Text]

25. Matysiak, S., A. Montesi, M. Pasquali, A. B. Kolomeisky, and C. Clementi. 2006. Dynamics of polymer translocation through nanopores: theory meets experiment. Phys. Rev. Lett. 96:118103.[CrossRef][Medline]

26. Muthukumar, M. 1999. Polymer translocation through a hole. J. Chem. Phys. 111:10371–10374.[CrossRef]

27. Ambjornsson, T., S. P. Apell, Z. Konkoli, E. A. Di Marzio, and J. J. Kasianowicz. 2002. Charged polymer membrane translocation. J. Chem. Phys. 117:4063–4073.[CrossRef]

28. Mathe, J., A. Aksimentiev, D. R. Nelson, K. Schulten, and A. Meller. 2005. Orientation discrimination of single-stranded DNA inside the alpha-hemolysin membrane channel. Proc. Natl. Acad. Sci. USA. 102:12377–12382.[Abstract/Free Full Text]

29. Slonkina, E., and A. B. Kolomeisky. 2003. Polymer translocation through a long nanopore. J. Chem. Phys. 118:7112–7118.[CrossRef]

30. Butler, T. Z., J. H. Gundlach, and M. A. Troll. 2006. Determination of RNA orientation