Unfolding Kinetics of the Human Telomere i-Motif Under a 10 pN Force Imposed by the α-Hemolysin Nanopore Identify Transient Folded-State Lifetimes at Physiological pH.

Cytosine (C)-rich DNA can adopt i-motif folds under acidic conditions, with the human telomere i-motif providing a well-studied example. The dimensions of this i-motif are appropriate for capture in the nanocavity of the α-hemolysin (α-HL) protein pore under an electrophoretic force. Interrogation of the current vs time (i-t) traces when the i-motif interacts with α-HL identified characteristic signals that were pH dependent. These features were evaluated from pH 5.0 to 7.2, a region surrounding the transition pH of the i-motif (6.1). When the i-motif without polynucleotide tails was studied at pH 5.0, the folded structure entered the nanocavity of α-HL from either the top or bottom face to yield characteristic current patterns. Addition of a 5' 25-mer poly-2'-deoxyadensosine tail allowed capture of the i-motif from the unfolded terminus, and this was used to analyze the pH dependency of unfolding. At pH values below the transition point, only folded strands were observed, and when the pH was increased above the transition pH, the number of folded events decreased, while the unfolded events increased. At pH 6.8 and 7.2 4% and 2% of the strands were still folded, respectively. The lifetimes for the folded states at pH 6.8 and 7.2 were 21 and 9 ms, respectively, at 160 mV electrophoretic force. These lifetimes are sufficiently long to affect enzymes operating on DNA. Furthermore, these transient lifetimes are readily obtained using the α-HL nanopore, a feature that is not easily achievable by other methods.

Introduction

The instructional information contained in DNA and RNA polymers can direct a wealth of secondary structures via hydrogen bonding, π stacking, and coordination with metals. Secondary structural variability observed for these polymers includes many different forms of hairpins, duplexes, cruciform folds, triple helices, G-quadruplexes, triplexes, and i-motifs.[1] The intercalation of cytosine (C)-rich strands at low pH yields i-motif folds that are finding broad utility in chemistry and biology.[2,3] Specifically, hemiprotonation of C-rich strands at low pH allows C+•C base pairs to form that initiate folding of i-motifs (Figure A). The pH dependency in formation of the i-motif topology causes the folding process to be highly cooperative.[2−4] These folds have immense potential in biosensor development and other analytical applications, due to the ease with which these folds can be manipulated by pH.[5] Biological studies have identified i-motifs in human promoter regions that are proposed to aid in the control of gene transcription.[3] For example, the Hurley laboratory has identified the BCL2 promoter i-motif as a cellular target for drugs to bind and control transcription of this gene.[6,7] The ability to regulate gene transcription via i-motif interacting drugs has enormous potential to combat a number of cancers.[8] Realization of the full potential of i-motif scaffolds in analytical applications and understanding their role in biological processes requires an extensive set of tools. These tools allow us to better understand the sequence dependency and how changes in physical conditions affect i-motif structures, their unfolding kinetics and thermodynamics as well as how small molecules interact with them.

Figure 1

Structures of the human telomere i-motif and the α-hemolysin nanopore. (A) The C+•C base pair scheme, the human telomere i-motif line structure, and cartoon drawing based on pdb 1ELN.[9] (B) Space filling model of the α-hemolysin nanopore (pdb 7AHL)[10] and the human telomere i-motif along with their dimensions for comparison.

Typical analytical tools for studying i-motif folds include spectroscopic methods (e.g., UV–vis, CD, and fluorescence), NMR, gel electrophoresis, thermal denaturation, and mass spectrometry.[11−13] These tools have identified numerous sequence-dependent folds for i-motifs, their optimum pH ranges for folding, and alterations of the physical parameters that affect i-motif stability.[13,14] Generally, hemiprotonation of C occurs at pH < 6, leading to favorable i-motif folding conditions, although molecular crowding agents and selection of certain sequences have found i-motif folds that are stable at higher pH values.[13−15] Recently, studies of individual i-motif folds by optical tweezers evaluated their pH-dependent tensile strengths.[16,17] Another single-molecule approach for studying DNA and RNA secondary structures is the alpha-hemolysin (α-HL) protein nanopore; however, this tool has not been employed for studying the i-motif scaffold.

The α-HL nanopore is a mushroom-shaped protein channel that, once embedded in a lipid bilayer, forms a well-defined, nanometer-sized hole through which ions can flow in an electrical field (Figure B).[10] On the cis side of α-HL is the opening (d ∼ 3.0 nm) to the nanocavity or vestibule that has a large volume (39.5 nm3)[18] for capturing and analyzing folded DNA structures. At the base of the nanocavity is a ∼1.4 nm central constriction leading to a β-barrel that exits the trans side of the protein.[18] Some protein nanopores allow single-stranded DNA to translocate readily, a feature being utilized as a next-generation DNA sequencing platform.[19−22] In contrast, folded DNA structures cannot pass the central constriction causing them to be trapped. When molecules are trapped and interact with the protein channel, they impede the ion flow.[23−25] Thus, the change in ion flow and the blockage duration time when molecules, such as DNA, interact with the nanopore yield valuable information about the molecular interactions. Specifically, the α-HL nanopore has successfully identified features of DNA secondary structures that have been challenging to observe by other bulk measurements. For example, the human telomere G-quadruplex can fold to two nearly identical hybrid structures in KCl-containing solutions, and the α-HL nanopore can readily identify the populations of these two topologies.[26] The thrombin-binding aptamer G-quadruplex studied with α-HL identified ion-regulated unfolding kinetic values that are beneficial for the development of fine-tuned biosensors.[27]

A second advantage to studies of DNA secondary structures with the α-HL nanopore is achieved by the size-selective properties of this nanochannel.[18] For instance, if molecules are too big (d > 3.0 nm) they do not enter the mouth of the nanocavity. Attachment of a long, single-stranded DNA tail allows folded DNA structures to be drawn into the pore. Large folds interact with the pore on the outside, while small structures can be trapped via the tail tethered through the β-barrel that allows the initial orientation to be known. This feature has verified that hairpins with two tails (internal hairpins) and propeller folded G-quadruplexes are too big to enter the protein channel when electrophoretically driven toward the channel.[28−30] Smaller structures, such as blunt-ended hairpins,[31,32] DNA duplexes,[33,34] and small G-quadruplexes enter the nanocavity of α-HL.[27,28,35,36] Encapsulation of these structures inside the nanocavity provides current levels and residence times that are diagnostic of specific structure- and sequence-dependent interactions. Other protein and solid-state nanopores have been utilized for interrogating DNA, aptamers, and protein structures when they interact with these small aperture pores under an electrical bias.[37−45] Single-molecule techniques using nanopores have the further advantage of providing information about the equilibrium or distribution of molecules in solution.[26,46,47] These studies provide background knowledge that can be leveraged for evaluating data collected from uncharted DNA structures, such as the i-motif scaffold, when they were allowed to interact with α-HL under an electrophoretic force.

In the current studies, the α-HL nanopore was utilized to interrogate the i-motif that folds from the human telomere repeat sequence (5′-CCC TAA-3′). Studies were conducted to understand whether and how polynucleotide tails attached to this structure modulate the DNA–protein interactions. Because i-motifs cooperatively fold in a pH-dependent fashion, titration studies were conducted to monitor this folding and unfolding process in situ. Additionally, the residence time that an i-motif was trapped in the nanocavity before unraveling provides a measure of the structure’s stability, allowing us to probe the unraveling kinetics of this structure as a function of pH. These studies identify unique and characteristic properties of i-motifs and their interactions with α-HL. This knowledge will be beneficial for stochastic biosensor development and for understanding the kinetic properties of these structures that may regulate genetic processes.

Results and Discussion

Initial studies were conducted with the minimal human telomere i-motif sequence 5′-CCC TAA CCC TAA CCC TAA CCC-3′ that does not have polynucleotide tails protruding from the core of the fold (Figure A). To obtain satisfactory current levels for interpreting interactions of the i-motif with α-HL under an electrophoretic force, the ionic strength was held at 1.0 M KCl and was appropriately buffered at each pH studied. Before conducting α-HL nanopore measurements, a set of pH-dependent CD studies were performed with this sequence under the high ionic strength conditions that favor DNA capture in the nanopore to determine the transition pH for folding this i-motif. Measurement of the CD spectra from 220–320 nm as the pH changed by 0.25 increments from 5.0 to 7.5 provided signals that changed with the pH (Figure ). Plots of the molar ellipticity at 286 nm vs pH demonstrate a transition with an inflection point at 6.1 (pHT). This value is similar to literature sources for this sequence under lower ionic strength conditions.[9] Knowledge of this transition pH allowed α-HL experiments to be designed that ensured monitoring of the fully folded and fully unfolded states. Lastly, to better understand how the thermodynamics of the i-motif are affected by the high ionic strength, thermal melting (Tm) measurements were conducted. At pH 5.0 and physiologically relevant ionic strength (140 mM KCl and 12 mM NaCl), the Tm was 43.6 °C, and at the high ionic strength conditions (1.0 M KCl) used for the nanopore studies, the Tm was 46.7 °C (Figure B). The high salt concentration increased the Tm by 3.4 °C, indicating that the high ionic strength of the nanopore conditions has a minor effect on increasing the thermodynamic stability of the folded state.

Figure 2

pH-Dependent CD spectra and ionic strength-dependent Tm values for the human telomere i-motif. (A) The CD spectra were recorded on a 10 μM solution of the sequence 5′-CCC TAA CCC TAA CCC TAA CCC-3′ in a 25 mM three buffer system (acetate, phosphate, and borate) with 1.0 M KCl at 20 °C. The pH values were adjusted from 5.0 to 7.5 in 0.25 pH unit increments. The inset represents the molar ellipticity at 286 nm vs pH that determined the transition pH (pHT) to be 6.1. (B) Table of measured Tm values for the i-motif sequence at physiological ionic strength (140 mM KCl and 12 mM NaCl) vs the high ionic strength for the nanopore analysis (1.0 M) at pH 5.0. Values were determined by monitoring the temperature-dependent unfolding process by UV–vis at 295 nm.

Analysis of the fully folded i-motif by α-HL was conducted at pH 5.0. When the i-motif was added to the cis side of the nanopore and electrophoretically driven into the nanocavity under a −120 mV bias (cis vs trans), only very long blockages to the current were observed. These blockages persisted for >2 min (Figure S1). Studies with greater driving force (−160 mV) gave similar results. The observation that the i-motif folds were trapped in the nanocavity is consistent with the exterior dimensions of the i-motif (d ∼ 2.0 nm, Figure B) being greater than the central constriction of α-HL (d ∼ 1.4 nm, Figure B). Furthermore, these initial observations at pH 5.0 identify the i-motif to be very stable and the kinetics of unraveling to be too slow to be monitored in a reasonable amount of time. All nanopore data reported were determined on at least three α-HL nanopores, and the number of events was dependent on the experiment as described below.

To obtain populations of blockage currents required for conducting statistical analysis of the results, the i-motifs were trapped in the nanocavity and held there for 2 s (Figure A). After the 2 s hold, the voltage bias was reversed to force the structure to exit the cis side of the channel, followed by reversing the current back again to allow trapping and studying the next DNA molecule. Analysis of >200 of these trapping events found two types of current vs time (i–t) patterns to emerge (Figure B and Figures S2 & S3). Close inspection of the i–t patterns revealed the characteristic patterns of these events that were termed Types 1 and 2 (Figures B). Type 1 event profiles were interpreted as entry of the i-motif into the channel causing the open channel current (Io) to initially go to an intermediate blockage state (IM) that lasted >100 μs and then began to transition to a deep blockage current (ID) that also lasted >100 μs. An alteration in the two blocking current levels was maintained throughout the 2 s duration for which each Type 1 event monitored. The current values reported were normalized by the Io to give percent residual currents (i.e., %I/Io). Accordingly, the %IM/Io was 36 ± 1%, and the %ID/Io was 10 ± 1% for the Type 1 events. In contrast, the Type 2 events initiated with an intermediate current blockage level of %IM/Io = 31 ± 1%, and this level was the dominate current recorded during the event (Figure B). Occasionally, a deeper blockage to the current was recorded that had a value of %ID/Io = 10 ± 1%, and these blockages lasted for <100 μs (Figure B). For the Type 2 events, this current pattern persisted through the entire 2 s during which they were recorded. The distribution of Type 1 and 2 events was 5.6:1, respectively (Figure C). These were the only two event types observed for the i-motif at pH 5.0, leading us to hypothesize the interpretation below.

Figure 3

Analysis of the human telomere i-motif with the α-HL nanopore at pH 5.0. (A) A 20 s continuous i–t trace showing the raw data. (B) Expanded views for the Type 1 and 2 events to show their unique i–t patterns. (C) Current vs IM counts to identify the Type 1 and 2 events based on their unique differences. (D) Proposed model that describes the process leading to these two event types. Data were recorded at −120 mV (cis vs trans) at 22 °C in 1.0 M KCl electrolyte solution buffered with 10 mM phosphate and 5 mM citrate at pH 5.0.

Inspection of the space filling model of the human telomere i-motif (Figure B) identifies it to be nearly cylindrically shaped with a height of 3.1 nm and a diameter of 2.0 nm similar to that of duplex B-form DNA.[9] Furthermore, the top and bottom of this cylinder have different features: on one side, the 5′ and 3′ termini exist along with one edgewise loop, while the other side presents two edgewise loops (Figure A,B). Therefore, we hypothesize that these two event types result from entry of the i-motif into the nanocavity from the side with the 5′ and 3′ termini or from the side with the two edgewise loops to yield the two different i–t patterns (Figure D). The different sides should interact with the narrow central constriction differently, resulting in the two different event types. Specifically, the side with one edgewise loop is proposed to be more flexible leading to more frequent and longer-lived interactions with the central constriction, such as those observed in the Type 1 events. In contrast, the side with two edgewise loops is conformationally restricted due to a loop–loop base pair that caps the structure[48] leading to less frequent interactions with the central constriction, such as those observed in the Type 2 events.

Support for this orientation-specific hypothesis was achieved by studying an i-motif sequence that had two additional nucleotides on both the 5′ and 3′ ends (sequence = 5′- AA CCC TAA CCC TAA CCC TAA CCC TA-3′, the underlined nucleotides being the ones added). These additional nucleotides do not change the pH-dependent folding properties of the sequence (Figure S4). These additional nucleotides should cause a change in the event type that enters the pore from the side with the 5′ and 3′ tails; thus, identifying the event type resulting from this face interacting with the central constriction. Examination of the results when the new sequence was studied with α-HL found that Type 2 events remained the same and that Type 1 events were changed (Figure S5). New features of the Type 1 events with the additional two-nucleotide tails include similar IM and ID currents, while the duration times dramatically changed (Figure ). Due to the highly stochastic nature of these new events, a thorough analysis was not conducted. These new events could not be unraveled at increased voltage (−160 mV) up to 2 min (Figure S6), suggesting the stability of the i-motif was not affected by the tails. Because the height of the i-motif is 3.1 nm and the mouth of the nanocavity is 3.0 nm, we suspect entry from the side would not be a favorable mode for entry of the i-motif into the nanocavity. There were no additional event types that would support this type of process.

Figure 4

Comparison of Type 1 events for i-motifs without tails and with two nucleotide tails on the 5′ and 3′ ends. The i–t traces represent 20 ms of two representative events. Data were recorded at −120 mV (cis vs trans) at 22 °C in 1.0 M KCl electrolyte solution buffered with 10 mM phosphate and 5 mM citrate at pH 5.0.

In the next study, the tailless i-motif was initially analyzed with a protein channel at −120 mV and pH 5.0, and then a predetermined aliquot of NaOH was added to the cis side of the channel to raise the pH to 6.8 (Figure ). At pH 6.8, this i-motif should be nearly completely unfolded based on the CD studies (Figure A). Upon addition of NaOH and within a few seconds (mixing time in the analysis chamber), new events were observed characteristic of unfolded, single-stranded DNA translocation, and the long dwell-time blockages were reduced to <5% of the events. The new events were very fast and gave a wide distribution of blocking levels that was expected because of the short length of this strand (21 nucleotides, see Figure ). For example, at −100 mV the ID distribution was broad and centered at %ID/Io = 11 ± 1%. The mean translocation time for a population of events (>2000 events) was a Gaussian distribution with a time constant of 0.03 ± 0.01 ms (Figure S7). Support for translocation of these unfolded strands through the channel was achieved by a voltage-dependent study that found as the voltage was increased, the translocation time decreased (Figure S7). The results of these studies at pH 5.0 and 6.8 verify that the nanopore can easily distinguish the folded and unfolded states of the i-motif strand.

Figure 5

Comparison of the human telomere i-motif studied at pH 5.0 and 6.8 with the α-HL nanopore. Data were recorded at −120 mV (cis vs trans) at 22 °C in in 1.0 M KCl electrolyte solution buffered with 10 mM phosphate and 5 mM citrate at pH 5.0 or 6.8 on the sequence 5′-CCC TAA CCC TAA CCC TAA CCC-3′.

Previous studies from our laboratory and others have found that the addition of a long (>25 nucleotide), single-stranded polynucleotide tail to bulky secondary structures such as hairpins, duplexes, and G-quadruplexes forces the structures to preferentially enter from the single-stranded terminus.[28,29,34,41,49−51] Also, the tail allows trapping the folds for very long times, and it allows the entry orientation to be known.[28,34] Accordingly, a 25-mer poly-2′-deoxyadenosine tail (dA25) was added to the 5′-end of the i-motif. Selection of the dA25 tail was advantageous because the structure of poly dA strands is not affected by low pH. First, we determined the addition of the 25-mer tail to the i-motif had no effect on the transition pH between the single-stranded and folded states based on CD analysis (Figure S8). Next, the i-motif bearing the tail was studied with α-HL at pH 5.0 under a −120 mV bias (Figure A). In these studies, two event types were observed that were incapable of translocating through the channel at −120 mV (Figure B); these studies determined the events to persist for >2 min (Figure S9). The two event types were termed Types 1′ and 2′ (Figure B). The Type 1′ events were characterized as having a static deep blockage to the current with %ID/Io = 10 ± 2%, while the Type 2′ events fluctuated between an intermediate current level with %IM/Io = 31 ± 1% and a deep blockage current level with %ID/Io = 10 ± 1% (Figure B). The intermediate current level was the dominate current recorded with only short pulses (<100 μs) of deep blockage current levels. The Type 2′ events behaved nearly identically to the Type 2 events recorded for the i-motif without a long poly dA tail (Figures B and 6B). On the basis of a comparison to our previous results,[28,30] the Type 1′ events are consistent with entry of the poly dA tail into the protein channel first, in which the tail occupies the entire β-barrel, and the i-motif resides in the nanocavity (Figure B,C). These studies with the 25-mer tail further support the model proposed in which the Type 1 and 2 events result from side-specific entry of the i-motif into the nanocavity (Figure D). The 5′-poly dA tail only affects the Type 1′ events, the side of the i-motif from which the tail protrudes from the i-motif core; whereas, Type 2′ events are nearly identical to the tailless i-motif Type 2 events.

Figure 6

Analysis of the i-motif with a 5′ 25-mer poly dA tail at pH 5.0 with the α-HL nanopore. (A) A 30 s continuous i–t trace. (B) Expanded views of the Type 1′ and 2′ events and graphical models showing the i-motif interactions with α-HL that yield these event types. (C) Comparison of the IM histograms for the i-motif with a 5′-tail (top) and tailless i-motif (bottom) to show that the IM values for the Type 2 and 2′ events are nearly identical. Data were recorded at −120 mV (cis vs trans) at 22 °C in 1.0 M KCl electrolyte solution buffered with 10 mM phosphate and 5 mM citrate at pH 5.0.

To obtain a sufficient number of events for statistical analysis, the event profiles were recorded for 2 s followed by ejection of the sample from the nanocavity by reversing the electrical bias. Once the channel was cleared, the bias was returned to its original state (−120 mV, cis vs trans) to capture another i-motif. Analysis of this population of data identified the Type 1′ and 2′ events to occur in an 11:1 ratio, respectively (Figure C). This ratio identifies the poly dA tail to enter the pore more readily than the i-motif end of the strand. Now that the event types have been identified, pH-dependent studies could next be conducted.

To conduct the pH-dependent studies, the i-motif with a 25-mer poly dA tail was placed on the cis side of the channel initially at pH 5.0, and the bias was set to −120 mV (Figure A). Once a sufficient population of data was recorded (>200 events), a predetermined aliquot of NaOH was added to the cis side of the channel to raise the pH to 5.7 (Figure A). Again, a sufficient population of data was recorded (>200 events). The only events observed were the folded i-motifs (Figure A). Comparison of the data at pH 5.0 and 5.7 determined them to be virtually identical (Figure A), which was consistent with the pH-dependent CD studies (Figure S8) that identify these two pH conditions to behave the same (Figures S12 and S13). Next, the pH was adjusted to 6.3, a value near the transition pH of 6.1, and three event types were observed. Two long event types identical to Types 1′ and 2′ were detected (28 ± 4%) as well as fast events identical to unfolded, single-stranded DNA translocation (72 ± 14%, Figures B and S16). This observation is consistent with the event profile expected during the transition phase of an i-motif from folded to the unfolded state. In the final pH analysis, NaOH was added to achieve pH 6.8 (Figures C and S17). Under these conditions, the i-motif gave predominantly unfolded states (96 ± 2%) and a small population of folded states (4 ± 1%, Figure C). Analysis of the population of times for single-stranded DNA translocation at pH 6.3 and 6.8 found them to be Gaussian distributed and to have nearly identical time constants (Figure S16). These studies provided knowledge of the current patterns expected as the i-motif goes through a transition from the folded to unfolded states.

Figure 7

The pH-dependent analysis of the i-motif with a poly dA tail via the α-HL nanopore at pH 5.0, 5.7, 6.3, and 6.8. (A) Typical i–t profiles observed for 40 s at pH 5.0 and 5.7. The data for both pH studies were very similar, and therefore, only data for pH 5.0 are shown. (B) Typical i–t profiles observed for 10 s at pH 6.3. (C) Typical i–t profiles observed for 6 s at pH 6.8. All i–t profiles were recorded at −120 mV (cis vs trans) at 22 °C in 1.0 M KCl electrolyte solution buffered with 10 mM phosphate and 5 mM citrate at the designated pH. Errors for the event frequencies are ∼20% based on analysis of events from three different protein channels.

In the last set of experiments, the i-motif with a 25-mer poly dA tail was analyzed at −160 mV (cis vs trans) bias while repeating the pH-dependent study. Interestingly, with the 25-mer tail, the i-motif at pH 5.0 could be unfolded at −160 mV when the tail enters first (Type 1′ events), a feat not achievable in 3 min at −120 mV (Figure S9). A voltage-dependent study (−120, −160, and −180 mV) confirmed that the events represented DNA strands unzipping and translocating through the pore, because as the voltage was increased, the time of the Type 1′ events decreased (Figures S9–S11). The less frequent Type 2′ events were readily identifiable and ejected from the channel and omitted from the analysis that follows. Additionally, when the pH was high enough to see single-stranded DNA translocations, these fast events were also not part of the analysis that follows. In the first case, the pH was initially at 5.0. The Type 1′ events represent DNA captured in the pore and held until they unzipped to obtain 150 events that ranged in translocation time from ∼5 s up to 150 s (Figure A,B). The median translocation time at pH 5.0 was determined to be 25.0 s (Figure A,B, n = 150 events). The broad distribution of the events is exponentially distributed, and to best describe the data, it is shown graphically in box plots in Figure A. The boxes represent the lifetimes for 75% of the data, and the bars show the 95% confidence interval. Inspection of these plots illustrates that as the pH was increased the lifetimes decreased. Further, we only report the median lifetimes in the text and in Figure B, while providing the error analysis graphically in Figure A. Next, the pH was increased to 5.7, and a population of 200 Type 1′ events was analyzed. Type 1′ events represented 100% of the population. From this population, the distribution of times ranged from ∼1 s up to 130 s, and the median time for translocation was determined to be 14.3 s (Figure A,B). Studies at pH 6.0 found a population of 250 Type 1′ events representing 70% of the population with a median lifetime of 8.1 s distributed from ∼1 to 40 s. When the pH was adjusted to 6.3, 300 Type 1′ events (28% of the population) were collected and analyzed. From this population of Type 1′ events, the spread in translocation time was <1 s up to 20 s with a median translocation time of 2.7 s (Figure A,B). Analysis of the i-motif unzipping at pH 6.8 showed that 4% of the population was Type 1′, and analysis of 300 of these events gave a median lifetime of 21 ms. The last pH studied was 7.2, in which Type 1′ events represented <2% of the population. The median lifetime of this population was 9.0 ms. These translocation times show a clear pH dependence that demonstrates as the pH was increased, the event time decreased. These low numbers of events are challenging to conduct statistical analysis; however, to a first approximation, the event populations appear to be exponentially distributed.

Figure 8

The pH-dependent times measured for the unraveling of the i-motif with a 25-mer poly dA tail at high bias. (A) Plots of unraveling time vs pH for individual molecules of the human telomere i-motif with a 5′-tail. (B) Table of median unraveling times vs pH. All measurements were recorded at −160 mV (cis vs trans) at 22 °C in a solution buffered with 10 mM phosphate and 5 mM citrate and 1.0 M KCl electrolyte.

The pH-dependent results provide insight about the lifetime of the i-motif folds when an electrophoretic force was applied. Comparison of the median i-motif lifetimes between pH 5.0 and 5.7 identifies that at the higher pH, the median lifetime decreases by nearly 2-fold (25.0 s vs 14.3 s, Figure B). The observation that the i-motif lifetime under a force decreased as the pH was increased was due to the increased likelihood of deprotonating the C+•C base pairs at the higher pH. The high degree of cooperativity[5] in the i-motif folding process supports the observation of the shorter lifetime at pH 5.7 relative to pH 5.0 (Figure B). The trend in i-motif lifetime continued to decrease as the pH was increased (Figure B). When the pH was 6.8, the measured lifetime was 21 ms, and at pH 7.2 the folded lifetime was 9.0 ms. Measurement of i-motif lifetimes as a function of pH is typically achieved by NMR measurements[52] that require large sample sizes or SPR analysis,[53] both of which require sophisticated, expensive instruments. Further, NMR and SPR results provide average lifetime values without an added stimulus or force. The millisecond folded lifetimes at pH 6.8 and 7.2 would be very challenging to measure by NMR, AFM, or SPR, if at all possible. Lastly, we evaluated the folded-state lifetime at pH 6.8 under physiological ionic strength (150 mM) and observed the median folded-state lifetime to be similar to that observed at 1.0 M ionic strength (Figure S21). The α-HL nanopore results presented here profile molecules individually and identify the fact that the i-motif lifetime population was broad under a force load and yields a nearly exponential distribution of times. This distribution was expected, due to the first-order processes that regulate unzipping of the strands in the α-HL nanopore.[29,49]

The results establish the α-HL nanopore as a convenient method for examining i-motif lifetimes as a function of pH under a force. The present approach for generating the force is analogous to the forces generated by RNA polymerases and DNA helicases on genomic DNA.[54,55] These enzymes apply forces on DNA from one direction, similar to the way the electrophoretic force does in the α-HL nanopore.[29] Further, the force applied by the nanopore at −160 mV is ∼10 pN,[56] a value that is similar to that imposed by RNA polymerase II (5–20 pN)[54] and DNA helicases (6–16 pN).[55] Utilization of the nanopore method allows measurement of the i-motif folded lifetime distribution with an applied force similar to these DNA processing enzymes. Understanding the lifetime an i-motif will maintain its folded state at physiological pH values and will help gauge the ability of the fold to interfere with biological processes. For example, RNA polymerase II synthesizes mRNA at rate of ∼15 ms per nucleotide.[57] Any DNA structure that can maintain its fold under the force of RNA polymerase II will cause it to stall, impacting the rate of mRNA synthesis. In this study, we determined that at pH values similar to those found in the nucleus (pH 6.8–7.0),[58] the human telomere i-motif has a folded-state lifetime of ∼20 ms (Figure B). This is long enough to cause RNA polymerase II to stall. This knowledge is readily obtained by nanopore measurements, yet quite challenging to determine by NMR, SPR, or optical tweezer methods. The human telomere i-motif does not regulate transcription; however, future α-HL nanopore studies on the BCL2(6,7) or c-MYC(59) i-motifs will allow one to understand the lifetimes of these folds under controlled physical conditions, such as those presented in these studies. Determination of these values will provide a window of times in which i-motif folds can regulate transcription, for example, or the length of time these folds may exist and slow helicase activity. Knowledge of these lifetimes will aid in our understanding of these biologically relevant structures that are targets for drug design to treat a number of cancers.[3]

Lastly, the α-HL nanopore is a stochastic sensor that can identify how solutes in a solution change when the equilibrium is perturbed.[23,60] In these studies, the perturbation was pH, and the change was the folding state of the i-motif. The ability to monitor this change in situ provides sensing capabilities of pH. Sensors built around the pH-dependent folding of i-motifs have been utilized to monitor a number of enzymatic reactions and quantify analytes in solution.[2,61−63] The present work provides an alternative method to determine how the equilibrium of i-motif folding is affected by external stimuli for reporting purposes.

Conclusions

The present work utilized the nanocavity of the α-HL nanopore (Figure B) to monitor the pH-dependent folded structures for the human telomere i-motif (Figure A). When the C-rich sequence without a tail was folded (pH 5.0), the structure entered the nanocavity from two different directions that gave characteristic i–t signatures (Figure C). These folded states were found to be long-lived and required a reversal of the electrical bias to eject them from the pore. Upon increasing the pH to 6.8, this sequence unfolded to the single-stranded state, and strands quickly translocated through the pore (Figure ). Upon addition of a 5′ 25-mer poly dA tail, similar i–t patterns were observed at −120 mV (Figure B). A titration of the pH was conducted to monitor the transition from the folded to unfolded state (Figure ). This experiment lead to increasing the voltage to −160 mV, at which point the applied force was capable of unraveling the i-motif folds, allowing determination of the folded-state lifetime under an electrophoretic force. These studies identified that even though the i-motif was folded in a range from pH 5.0 to 7.2, the lifetime of the fold under a force was pH dependent (Figure B). The present results establish a method for analyzing lifetimes of i-motif sequences found in gene promoter regions that are important for biological regulation.[3] Furthermore, the stochastic sensing capability of the α-HL nanopore[23] provides a convenient way to interrogate the folded state of an i-motif in solution when the folding equilibrium is dependent on the physical properties of the solution (i.e., pH) or an analyte that specifically binds these sequences.[2,61−63] Future studies with α-HL and i-motifs are anticipated to be applicable in analytical sensing applications as well as biophysical studies that aim to understand the biology of nucleic acids.

Methods

Ion channel recordings, DNA preparation procedures, and ion channel measurements were conducted as previously reported.[29,64,65] Details about DNA purification, nanopore construction, and CD spectroscopy procedures can be found in the Supporting Information. The ion channel measurements with different pH values were initiated at pH 5.0 (10 mM KPi and 5 mM citrate) with 1.0 M KCl electrolyte at 22 °C. After collection of enough data at pH 5.0 with one protein channel, 5 μL of 0.2 M NaOH was added to the nanopore chamber to reach pH 5.7. After mixing the solution in the nanopore chamber with a pipet tip, data were collected at the new pH. This process was repeated for pH 6.3 and 6.8. All pH-dependent measurements were collected with the same α-HL nanopore. The pH values for each titration step were determined by placing all components of the analysis chamber in an Eppendorf tube, and then a pH meter was used to determine how much NaOH was needed at each step to obtain the desired pH. The data were collected with a 100 kHz filter with an acquisition rate of 500 kHz. For presentation purposes the data were refiltered to 2 or 10 kHz. All data were processed QuB 1.5.0.31 and plotted using OriginPro 9.1.