Prolonged Residence Time of a Noncovalent Molecular Adapter, ß-Cyclodextrin, within the Lumen of Mutant {alpha}-Hemolysin Pores

Li-Qun Gua, Stephen Cheleya, and Hagan Bayleya,b
a Department of Medical Biochemistry and Genetics, The Texas A&M University System Health Science Center, College Station, TX 77843
b Department of Chemistry, Texas A&M University, College Station, TX 77843

Correspondence to: Hagan Bayley, Department of Medical Biochemistry and Genetics, The Texas A&M University System Health Science Center, 440 Reynolds Medical Building, College Station, TX 77843-1114. Tel.: (979) 845-7047; Fax:(979) 847-9481 E-mail:bayley{at}tamu.edu.


  Abstract
Top
Abstract
INTRODUCTION
MATERIALS AND METHODS
RESULTS
DISCUSSION
Acknowledgements
References

Noncovalent molecular adapters, such as cyclodextrins, act as binding sites for channel blockers when lodged in the lumen of the {alpha}-hemolysin ({alpha}HL) pore, thereby offering a basis for the detection of a variety of organic molecules with {alpha}HL as a sensor element. ß-Cyclodextrin (ßCD) resides in the wild-type {alpha}HL pore for several hundred microseconds. The residence time can be extended to several milliseconds by the manipulation of pH and transmembrane potential. Here, we describe mutant homoheptameric {alpha}HL pores that are capable of accommodating ßCD for tens of seconds. The mutants were obtained by site-directed mutagenesis at position 113, which is a residue that lies near a constriction in the lumen of the transmembrane ß barrel, and fall into two classes. Members of the tight-binding class, M113D, M113N, M113V, M113H, M113F and M113Y, bind ßCD ~104-fold more avidly than the remaining {alpha}HL pores, including WT-{alpha}HL. The lower Kd values of these mutants are dominated by reduced values of koff. The major effect of the mutations is most likely a remodeling of the binding site for ßCD in the vicinity of position 113. In addition, there is a smaller voltage-sensitive component of the binding, which is also affected by the residue at 113 and may result from transport of the neutral ßCD molecule by electroosmotic flow. The mutant pores for which the dwell time of ßCD is prolonged can serve as improved components for stochastic sensors.

Key Words: {alpha}-toxin, molecular adapter, mutagenesis, pore, stochastic sensing


  INTRODUCTION
Top
Abstract
INTRODUCTION
MATERIALS AND METHODS
RESULTS
DISCUSSION
Acknowledgements
References

Considerable effort has been spent on the engineering of soluble proteins, such as enzymes and antibodies, both by natural and unnatural amino acid substitution, and by targeted chemical modification. By contrast, the engineering of transmembrane channels and pores is a relatively unexplored area (Bayley 1999 ) that is worth examination because membrane proteins differ in environment and structure compared with soluble proteins. Further, just as in the case of soluble proteins, applications for engineered channels and pores are emerging in several areas of biotechnology (Bayley 1999 ).

To further efforts in both the protein engineering and the biotechnology of membrane proteins, we have been working with the bacterial pore-forming toxin, staphylococcal {alpha}-hemolysin ({alpha}HL).* The {alpha}HL pore is a heptamer made up of identical subunits of 293 amino acids (see Fig 1 A). Roughly globular molecules with molecular masses of up to ~2,000 D (Fussle et al. 1981 ), or larger elongated polymers such as single-stranded nucleic acids (Kasianowicz et al. 1996 ), can pass through a channel centered on the molecular sevenfold axis of the pore (see Fig 1 A). Before the crystal structure of the {alpha}HL pore was solved (Song et al. 1996 ), we used biochemical knowledge to engineer triggers and switches into the protein through which the assembly pathway could be controlled by a variety of agents including metal ions, enzymes, and light (Bayley 1995 , Bayley 1997 ). Potential applications of these molecules include the use of metal-regulated pores for reversible cell permeabilization (Russo et al. 1997 ; Otto-Bruc et al. 1998 ; Eroglu et al. 2000 ) and the development of cytotoxic agents activated by cell-surface proteases (Walker and Bayley 1994 ; Panchal et al. 1996 ).



View larger version (48K):
[in this window]
[in a new window]
 
Figure 1. Representations of staphylococcal {alpha}-hemolysin ({alpha}HL) and ß-cyclodextrin (ßCD). (A) Sagittal section through the WT-{alpha}HL pore showing the location of Met-113 (dark green), Glu-111 (red), Lys-147 (blue), Asn-139 (yellow), and Leu-135 (green). (B) Structure of ßCD. (C) Schematic of the WT-{alpha}HL pore showing ßCD lodged in the lumen of the channel. The location is based on mutagenesis data (Gu et al. 1999 , Gu et al. 2001 ), including that described here. (D) Sequence of the transmembrane ß barrel of {alpha}HL-RL2. WT residues are shown in parentheses.

The availability of the structure of the {alpha}HL pore (Song et al. 1996 ) allowed us to focus on the properties of fully assembled pores altered by direct mutagenesis (Braha et al. 1997 ) or by targeted covalent modification (Movileanu et al. 2000 ; Howorka et al. 2001a , Howorka et al. 2001b ). From the viewpoint of biotechnology, the primary goal has been to produce pores that respond to various analytes as components for stochastic sensors operating at the single molecule level (Bayley et al. 2000 ; Bayley and Cremer 2001 ). For example, direct mutagenesis has permitted the detection of divalent metal ions (Braha et al. 1997 , Braha et al. 2000 ), whereas the covalent attachment of ligands has allowed the detection of proteins (Movileanu et al. 2000 ) and DNA (Howorka et al. 2001a ).

The stochastic sensing of small organic molecules was a substantial challenge that was solved by the introduction of noncovalent molecular adapters (Fig 1B and Fig C; Gu et al. 1999 ). Adapter molecules, such as cyclodextrins, can become lodged in the lumen of the {alpha}HL pore, where they remain available for the host–guest interactions that are well-known to occur in solution (D'Souza and Lipkowitz 1998 ). Hence, binding of a cyclodextrin within the lumen of the {alpha}HL pore causes a reduction in current flowing through the pore. Additional transient reductions in the current are brought about by the subsequent binding of organic analyte molecules to the cyclodextrin, permitting the analyte to be identified and quantified (Gu et al. 1999 ). Other molecules, notably cyclic peptides, can act as adapters (Sanchez-Quesada et al. 2000 ). The adapters can bring about further changes in the properties of the {alpha}HL pore; for example, they can alter ion selectivity (Gu et al. 2000 ; Sanchez-Quesada et al. 2000 ).

Because of the utility of noncovalent molecular adapters, it is important to maximize their dwell times within the lumen of the pore. The dwell time of ß-cyclodextrin (ßCD) within the pore can be extended to several milliseconds at low pH and high negative transmembrane potentials or at high pH and high positive potentials (Gu and Bayley 2000 ). However, it would be advantageous to lengthen the residency to seconds or longer. We have already shown that the mutant M113N binds ßCD with high affinity (Gu et al. 1999 ) and, recently, this finding was extended when a nanocavity was built between two different cyclodextrins lodged in the transmembrane ß barrel (Gu et al. 2001 ). Met-113 is located near the constriction at the internal end of the ß barrel (Fig 1 A); here, we report in detail the results of mutagenesis of Met-113 to each of the remaining nineteen natural amino acids.


  MATERIALS AND METHODS
Top
Abstract
INTRODUCTION
MATERIALS AND METHODS
RESULTS
DISCUSSION
Acknowledgements
References

Mutagenesis of {alpha}HL at Position 113
All Met-113 mutants were made by cutting pT7-{alpha}HL-RL2 (Cheley et al. 1999 ) with SacII and HpaI and replacing the excised internal fragment with duplex DNA formed from 5'-GGAATTCGATTGATACAAAAGAGTATxyzAGTACGTT-3' (sense) and 5'-AACGTACTz'y'x'ATACTCTTTTGTATCAATCGAATTCCGC-3' (antisense), where xyz and z'y'x' represent, respectively, the codon and anticodon replacements for the following: Ala (xyz = GCA/ x'y'z' = TGC), Cys (TGC/GCA), Asp (GAT/ATC), Glu (GAG/CTC), Phe (TTT/AAA), Gly (GGG/CCC), His (CAT/ATG), Ile (ATT/AAT), Lys (AAA/TTT), Leu (CTA/TAG), Asn (AAT/ATT), Pro (CCA/TGG), Gln (CAG/CTG), Arg (AGA/TCT), Ser (AGC/GCT), Thr (ACT/AGT), Val (GTT/AAC), Trp (TGG/CCA), and Tyr (TAT/ATA). Because these changes were made in the RL2 background (Cheley et al. 1999 ), each mutant contains the following additional replacements with respect to WT-{alpha}HL: Lys-8 -> Ala, Val-124 -> Leu, Gly-130 -> Ser, Asn-139 -> Gln, and Ile-142 -> Leu.

Additional {alpha}HL Mutants
Replacements of charged residues near the constriction were also made by constructing the mutants K147N, E111N/K147N, and E111N/M113N/K147N. First, mutants {alpha}HL-E111N and {alpha}HL-K147N were constructed by cutting pT7-{alpha}HL-RL2 with the enzyme pairs SacII-HpaI and AflII-XhoI, respectively, and replacing the excised internal fragments with, respectively, the duplexes formed from 5'-GGAATTCGATTGATACAAAAAATTATATGAGTACGTT-3' (sense) and 5'-AACGTACTCATATAATTTTTTGTATCAATCGAATTCCGC-3' (antisense) and from 5'-TTAATTATGTTCAACCTGATTTCAAAACAATTC-3' (sense) and 5'-TCGAGAATTGTTTTGAAATCAGGTTGAACATAA-3' (antisense). To construct the double mutant E111N/K147N, pT7-{alpha}HL-K147N was digested with HpaI and HindIII, and the resulting internal fragment was inserted into pT7-{alpha}HL-E111N that had been cut with the same enzymes. The triple mutant E111N/M113N/K147N was prepared by cutting pT7-{alpha}HL-K147N with SacII and HpaI and replacing the small fragment with duplex DNA prepared from 5'-GGAATTCGATTGATACAAAAAATTATAATAGTACGTT-3' (sense) and 5'-AACGTACTATTATAATTTTTTGTATCAATCGAATTCCGC-3' (antisense). Because these changes were made in the RL2 background (Cheley et al. 1999 ), each mutant contained the additional replacements listed above.

Mutants M113N(WT), N139Q(WT), and M113N/N139Q (WT) were constructed as described previously (Gu et al. 2001 ). {alpha}HL-L135N was constructed by cutting pT7-{alpha}HL-RL1 (Cheley et al. 1997 ) with SpeI and ApaI and replacing the fragment with duplex DNA prepared from 5' CTAGTAAAATTGGAGGCAATATTGGGGCCCAGG-3' (sense) and 5' GGCCCCAATATTGCCTCCAATTTTA-3' (antisense). The StuI site that is present in the {alpha}HL-RL1 gene, but absent in {alpha}HL-L135N, was used to screen for cassette replacement before DNA sequencing. RL1 encodes the same amino sequence as RL2, but contains different restriction sites. Therefore, in the text, L135N is treated as having the RL2 background. {alpha}HL-M113N/L135N was constructed in a similar fashion. pT7-{alpha}HL-M113N (RL2 background; see above) was cut with SpeI and ApaI, and the fragment was replaced with duplex DNA prepared from the same oligonucleotides used to construct {alpha}HL-L135N (see above). As before, plasmids were screened by digestion with StuI.

DNA Sequencing
The genes of all {alpha}HL mutants used in this work were sequenced entirely. Sequencing was performed with SC001 (upstream forward primer) 5'-CACTATAGGGAGACCACAACGG-3' and SC003 (internal forward primer) 5'-CAGGGTTTTCACCAGACTTCGC-3'. No changes were found, except for those intended.

{alpha}HL Homoheptamer Formation and Purification
Heptameric WT-{alpha}HL was formed by treating monomeric {alpha}HL, purified from Staphylococcus aureus, with deoxycholate (Bhakdi et al. 1981 ; Walker et al. 1992 ) and isolated from SDS–polyacrylamide gels as described previously (Braha et al. 1997 ). The remaining {alpha}HL polypeptides were synthesized in vitro by coupled transcription and translation (IVTT) and assembled into homoheptamers by the inclusion of rabbit red cell membranes during synthesis, as described previously (Cheley et al. 1999 ). The heptamers were purified by SDS-PAGE and stored in 50-µl aliquots at -80°C (Cheley et al. 1999 ). All mutant polypeptides were synthesized and purified twice.

Planar Bilayer Recordings
Planar bilayer recordings were made as described at 22 ± 2°C (Gu and Bayley 2000 ; Gu et al. 2000 ). Buffers for bilayer recording contained 1 M NaCl and 10 mM dibasic sodium phosphate (Sigma-Aldrich), in deionized water (Millipore Corp.), and were titrated to pH 7.5 with aqueous HCl (EMScience). The bilayer was formed from 1,2-diphytanoyl-sn-glycero-phosphocholine (Avanti Polar Lipids) over an orifice 50–100 µm in diameter (Montal and Mueller 1972 ). The orifice had been pretreated with hexadecane in pentane, and the lipid was transferred to the chambers in pentane. Protein was added to the cis chamber, which was at ground. A positive potential indicates a higher potential in the trans chamber of the apparatus, and a positive current is one in which cations flow from the trans to the cis side. Experiments were initiated by the addition of heptameric {alpha}HL to the cis chamber, to a final concentration of 3–30 ng ml-l, with stirring until a single channel appeared. ßCD (Sigma-Aldrich) was added to the trans chamber at 40 µM, unless otherwise specified. The amplifier's internal low-pass Bessel filter was set at 5 kHz. Data were acquired at a sampling rate of 20 kHz. The current recordings were analyzed essentially as described previously (Gu and Bayley 2000 ; Gu et al. 2000 ). To determine each set of kinetic constants, three or more experiments had been performed and, in each case, data were analyzed that were acquired for at least 2 min for weak binding mutants (class 2) and at least 1 h for tight binding mutants (class 1). {tau}on and {tau}off for ßCD for each mutant, for data obtained at either -40 mV or +40 mV and 40 µM ßCD, were obtained from dwell-time histograms fitted to single exponentials by the Levenberg-Marquardt procedure. In all cases, the coefficient of determination of the fits was R >= 0.85. The data were replotted in semilogarithmic form for display in the paper. Separate segments of the data yielded similar {tau} values, suggesting that stationary kinetics prevailed. Kinetic constants were calculated by using koff = 1/{tau}off, kon = 1/{tau}on[ßCD], and Kd = koff / kon, where [ßCD] is the concentration of ßCD. Values for unitary conductance, kon and koff, and Kd are quoted as the mean ± SD. Ion selectivity (PK+/PCl-) was determined from reversal potentials measured with the following solution: cis 1,000 mM KCl, 10 mM potassium phosphate, pH 7.5 (dibasic salt titrated with HCl), trans 200 mM KCl, and 10 mM potassium phosphate, pH 7.5 (Gu et al. 2000 ).

Online Supplemental Material
Supplemental figures are available at http://www.jgp.org/cgi/content/full/118/5/481/DC1. They show examples of dwell-time histograms used to obtain values of {tau}on for M113E, M113K, M113N, and M113W. The figures also display histograms obtained from independent segments of the recordings to show that the channel kinetics were stationary.


  RESULTS
Top
Abstract
INTRODUCTION
MATERIALS AND METHODS
RESULTS
DISCUSSION
Acknowledgements
References

The WT-{alpha}HL Pore and {alpha}HL-RL2, with the Background Used for Mutagenesis, Have Similar Properties
The Met-113 replacements examined in this work were made in {alpha}HL-RL2. RL2 is the product of a semisynthetic gene that was originally devised to permit ready cassette mutagenesis of the sequence encoding the transmembrane ß barrel. RL2, and the mutants described here, contain the following additional mutations over WT-{alpha}HL: Lys-8 -> Ala, Val-124 -> Leu, Gly-130 -> Ser, Asn-139 -> Gln, and Ile-142 -> Leu (Fig 1 D). Therefore, the properties of the WT and RL2 {alpha}HL pores were first compared. They were found to be similar in terms of unitary conductance values and rectification ratios (g+40mV/g-40mV; Table 1), but their properties do differ in detail. For example, the extent of channel block by ßCD at +40 mV is lower for RL2 (30% residual current) than for WT-{alpha}HL (35% residual current), and whereas WT-{alpha}HL•ßCD is incompletely blocked by a variety of organic molecules (e.g., adamantane-1-carboxylic acid; Gu et al. 1999 ), RL2•ßCD is almost completely blocked by the same molecules (unpublished data).


 
View this table:
[in this window]
[in a new window]
 
Table 1. Conductance Values and Kinetic Parameters for the Interaction of ßCD with WT-{alpha}HL and Met-113 Mutantsa

ßCD Exhibits a Prolonged Dwell Time within the Lumen of Certain Met-113 Replacement Mutants
The interaction of ßCD was tested with pores formed from 20 {alpha}HL polypeptides, representing Met and all 19 natural amino acid substitutions at position 113. As documented above, the mutations were made in the "RL2" background. In single-channel recordings, in all 20 cases, the addition of 40 µM ßCD to the trans compartment produced reversible partial blockades of the ionic current (Fig 2). The dwell times for ßCD in the lumen of the pore ({tau}off, level 2) were strikingly prolonged for certain mutants (Fig 2 and Fig 3 A), whereas the dwell times for the unoccupied state of each mutant (inter-event interval, {tau}on, level 1) did not vary as much and are similar to {tau}on values for ßCD and WT-{alpha}HL (Fig 3 A).



View larger version (47K):
[in this window]
[in a new window]
 
Figure 2. Representative current traces from single {alpha}HL pores showing the blockade of Met-113 mutants by ßCD. All traces were recorded under symmetrical conditions in buffer containing 1 M NaCl, 10 mM sodium phosphate, pH 7.5; 40 µM ßCD was added to the trans chamber. (left) Traces recorded at -40 mV; (right) traces recorded at +40 mV. The mutants shown are all derived from RL2 (see the first paragraph of RESULTS). The broken line indicates zero current. The mutants are ordered (top to bottom) according to increasing affinity for ßCD.




View larger version (39K):
[in this window]
[in a new window]
 
Figure 3. Dwell times for the interaction between ßCD and the Met-113 mutants. (A) Values of {tau}off, the dwell time of ßCD in the pore, and {tau}on, the inter-event interval at 40 µM trans ßCD. (closed gray circle) -40 mV; ({circ}) +40 mV. The mutants are ordered (left to right) according to increasing affinity for ßCD. (B) Dependence of {tau}off and {tau}on on the concentration of ßCD (trans) for WT-{alpha}HL and selected mutants.

Interestingly, most of the pores tested could be grouped into two classes according to their {tau}off values. ßCD binds to mutant pores of class 1, namely M113D, M113N, M113V, M113H, M113F, and M113Y, with extended dwell times ({tau}off) compared with WT-{alpha}HL (Fig 3 A). For example at -40 mV, the dwell time of M113N•ßCD ({tau}off = 27 s) is over 104-fold longer than that of WT-{alpha}HL•ßCD ({tau}off = 0.76 ms) or RL2•ßCD ({tau}off = 1.0 ms). By contrast, ßCD binds to the mutants in class 2 (all but two of the mutants not in class 1) with a dwell time similar to that of WT-{alpha}HL or RL2. For example at -40 mV, the dwell time for M113K•ßCD is {tau}off = 1.4 ms. M113W and M113P are borderline cases and were not included in class 1 or 2.

In contrast to {tau}off, the {tau}on values of the mutants do not vary greatly and are similar to WT-{alpha}HL. For example at -40 mV, the inter-event interval for M113N at 40 µM ßCD ({tau}on = 91 ms) is comparable to that of WT-{alpha}HL ({tau}on = 68 ms) or RL2 ({tau}on = 93 ms). Therefore, the mutants in class 1, with prolonged {tau}off values, show high affinity for ßCD, whereas mutants in class 2, with short {tau}off values, show a similar affinity for ßCD to WT-{alpha}HL (Fig 2 and Table 1). For example at -40 mV, M113N in class 1 (Kd = 1.3 ± 0.6 x 10-7 M) binds ßCD over 104-fold more tightly than WT-{alpha}HL (Kd = 3.4 ± 0.4 x 10-3 M) and RL2 (Kd = 3.3 ± 0.2 x 10-3 M), whereas M113K in class 2 (Kd = 2.1 ± 0.10-3 M) is similar to WT-{alpha}HL in affinity.

ßCD Binds to {alpha}HL Pores in a Simple Bimolecular Interaction
Level 2 is the only major current blockade level seen in the amplitude histograms of all 20 mutants, which suggests that there is only one major binding site for ßCD within the lumen of each {alpha}HL pore (Gu and Bayley 2000 ). The binding kinetics are also in keeping with this interpretation. As in the case of WT-{alpha}HL (Gu and Bayley 2000 ), {tau}on and {tau}off for ßCD for each mutant could be fitted by single-exponential distributions for data obtained at either -40 or +40 mV (1 M NaCl and 10 mM sodium phosphate, pH 7.5) and 40 µM ßCD in the trans chamber (Supplemental Material). Therefore, the kinetics of the interaction between ßCD and the mutant {alpha}HL pores most likely obey the simple Scheme 1, with koff = 1/{tau}off, kon = 1/{tau}on[ßCD], and Kd = koff/kon, where [ßCD] is the concentration of ßCD. The kinetic constants were calculated accordingly (Table 1).


Scheme 1.

In selected cases, WT-{alpha}HL, RL2, M113E, M113K, M113N and M113W, the concentration dependence of 1/{tau}on was examined and found to be proportional to [ßCD] (Fig 3 B), in further support of a simple bimolecular interaction between the pores and ßCD.

For mutants in class 1, short additional blockades could be seen while ßCD was bound at either -40 or +40 mV. Because their frequency of occurrence was independent of ßCD concentration (tested for M113N and M113D; unpublished data), these events are probably not due to the binding of a second ßCD. Instead, they may correspond to a second conformation of the occupied state, {alpha}HL•ßCD. Because, these events occupied <3% of the total ßCD binding time, they were merged with level 2 for the kinetic analysis.

ßCD Affinity and Side Chain Properties at Position 113
The properties of the side chains at position 113 in the tight-binding class 1 mutants show no obvious relationship with one other. For example, when 1/Kd, kon and koff are plotted versus side-chain van der Waals volume no pattern is seen, although this is a convenient way to display the data (Fig 4). Patterns are revealed when the data are examined in detail. Small side chains of <85 Å3 fall into class 2. Whereas M113N and M113D fall into class 1, mutants with homologous substitution, M113Q and M113E, are in class 2. Similarly, of mutants with nonaromatic hydrophobic amino acid substituents, the one with the lowest volume side chain, M113V, is in class 1, whereas the others M113L, M113I, and RL2 itself (Met-113) are in class 2. Mutants with aromatic substitutions fall into class 1, with the exception of the mutant with the bulkiest substituent, M113W, which has a Kd value intermediate between those characteristic of the two classes.



View larger version (39K):
[in this window]
[in a new window]
 
Figure 4. Plots of kinetic constants for the Met-113 mutants versus the van der Waals volume (Creighton 1993 ) of the residue at position 113. (A) 1/Kd; (B) kon; and (C) koff. (closed gray circle) -40 mV; ({circ}) +40 mV. Some of the points are obscured, but the values can be found in Table 1.

Single-channel conductance values also show no clear correlation with side-chain van der Waals volume (Fig 5A and Fig B, and Table 1). Although the presence of seven copies of the bulkiest residue, Trp, at position 113 does yield the pore with the lowest conductance, the effect is not dramatic. The rectification ratios (g+40 mV/g-40 mV) also appear to be independent of the van der Waals volume (Fig 5 C). The extent of block by ßCD (Fig 5 A, B, and D, and Table 1) and the rectification ratios with ßCD bound (Fig 5 C) are also uncorrelated with van der Waals volume. There is also no correlation between the examined electrical properties and whether a mutant is a member of class 1 or class 2.



View larger version (28K):
[in this window]
[in a new window]
 
Figure 5. Plots of single-channel conductance values for the Met-113 mutants with and without ßCD bound versus the van der Waals volume of the residue at position 113. (A) Conductance values at -40 mV: (closed gray circle), {alpha}HL conductance; ({circ}) {alpha}HL•ßCD conductance. (B) Conductance values at +40 mV. (C) Ratio of conductance at +40 mV to conductance at -40 mV (g+40/g-40). (D) Conductance of the pore with ßCD bound divided by the conductance of the pore itself. (closed gray square), -40 mV; ({square}) +40 mV.

The Interaction of ßCD with Pores Mutated at Position 113 Is Independent of Various Other Mutations in the ß Barrel
To check the effect of additional mutations in the ß barrel on the affinity of ßCD, five additional mutants were tested: M113N(WT), N139Q(WT), M113N/N139Q(WT), L135N(RL2), and M113N/L135N(RL2) (Fig 6). As noted earlier, mutants with the RL2 background contain the following additional mutations over WT-{alpha}HL: Lys-8 -> Ala, Val-124 -> Leu, Gly-130 -> Ser, Asn-139 -> Gln, and Ile-142 -> Leu.



View larger version (18K):
[in this window]
[in a new window]
 
Figure 6. Representative current traces from single {alpha}HL pores showing blockades by ßCD. All traces were recorded under symmetrical conditions in buffer containing 1 M NaCl, 10 mM sodium phosphate, pH 7.5. ßCD (40 µM) was present on the trans side of the membrane. (A) WT-{alpha}HL, (B) N139Q(WT), (C) L135N(RL2), (D) M113N(WT), (E) M113N/N139Q(WT), and (F) M113N/L135N(RL2). The broken line indicates zero current.

WT-{alpha}HL, N139Q(WT), and L135N(RL2) share the residue Met-113, and all three mutants bind ßCD weakly: WT-{alpha}HL, Kd = 3.4 x 10-3 M; N139Q(WT), Kd = 2.9 x 10-3 M; and L135N(RL2), Kd = 3.9 x 10-3 M. By contrast, M113N(WT), M113N/N139Q(WT) (also named PNQ; Gu et al. 2001 ) and M113N/L135N(RL2) share Asn at position 113, and all three bind ßCD strongly: M113N(WT), Kd = 2.1 x 10-7 M; M113N/N139Q(WT), Kd = 3.1 x 10-7 M; and M113N/L135N(RL2), Kd = 9.3 x 10-7 M. Therefore, the affinity of the {alpha}HL pore for ßCD is dependent on the residue at position 113, but independent of other sites, notably positions 135 and 139, at least for the substitutions examined here. Residues 135 and 139 project into the lumen of the pore (Fig 1) and we have previously shown that the mutation N139Q forms a binding site for a different cyclodextrin, hepta-6-sulfato-ßCD, but not for ßCD (Gu et al. 2001 ).

A similar pattern was observed when the extent of current block by ßCD was examined for the same mutants. ßCD reduced the conductance of WT-{alpha}HL from 651 to 240 pA (37% residual current), N139Q(WT) from 635 to 253 pA (40%), and L135N(RL2) from 640 to 243 pA (38%). In contrast, the residual conductance values for the pores containing the M113N mutation were as follows: M113N(WT) from 623 to 283 pA with ßCD (46%); M113N/N139Q(WT) from 668 to 285 pA (43%); and M113N/L135N(RL2) from 645 to 283 pA (45%). Therefore, the extent of current block by ßCD is more dependent on the amino acid at position 113 than those at positions 135 and 139.

The Voltage Dependence of ßCD Binding Is Correlated with the Charge Selectivity of the Mutant Pores
The affinity of each mutant {alpha}HL pore for ßCD is voltage dependent. ßCD binds to some mutants more weakly at positive transmembrane potentials than at negative potentials, but binds to others with the opposite dependence on potential (Fig 2 and Table 1). For example, ßCD binds to M113E more strongly at +40 mV (Kd = 3.3 ± 0.5 x 10-3 M) than at -40 mV (Kd = 1.2 ± 0.2 x 10-2 M), whereas ßCD binds to M113R more weakly at +40 mV (Kd = 1.3 ± 0.1 x 10-2 M) than at -40 mV (Kd = 1.6 ± 0.1 x 10-3 M).

These data were quantified by using ß1/K = log (Kd-40/Kd+40). ß1/K > 0 reflects a stronger affinity for ßCD at +40 mV than at -40 mV, and ß1/K < 0, reflects the opposite. Similarly, {alpha} = log (PK+/PCl-) was used as a measure of the charge selectivity of each mutant. Where {alpha} > 0, a pore is cation selective, and where {alpha} < 0, a pore is anion selective (Table 2). When {alpha} and ß1/K were displayed on a scatter plot, they were seen to be correlated (Fig 7 A). This means that ßCD binds to a cation-selective pore more strongly at positive potentials than at negative potentials, whereas ßCD binds to an anion-selective pore more strongly at negative potentials than at positive. Although the effect of mutagenesis on the affinity for ßCD was largely reflected in koff, the smaller effect of voltage was manifested in both kon and koff (Fig 7b and Fig c).





View larger version (53K):
[in this window]
[in a new window]
 
Figure 7. Relationships of the voltage dependence of kinetic constants for the interaction of {alpha}HL and ßCD, and the charge selectivity of the pore. (A) Plot of ß1/K versus {alpha}, where ß1/K = log[(1/Kd+40)/(1/Kd-40)] = log(Kd-40/Kd+40) for the mutants in Table 2. ß1/K > 0 reflects a stronger affinity at +40 mV than -40 mV, and ß1/K < 0, the opposite. {alpha} = log(PK+/PCl-), a measure of the charge selectivity of each mutant. Where {alpha} > 0, a pore is cation selective, and where {alpha} < 0, a pore is anion selective. A line was fitted to the data by linear regression. The correlation coefficient is R2 = 0.80. (B) Plot of ßkon versus {alpha}, where ßkon = log(kon+40/kon-40). R2 = 0.79. (C) Plot of ßkoff versus {alpha}, where ßkoff = log(koff+40/koff-40). R2 = 0.77.


 
View this table:
[in this window]
[in a new window]
 
Table 2. Charge Selectivity of Pores and Voltage Dependence of Affinity for ßCD


  DISCUSSION
Top
Abstract
INTRODUCTION
MATERIALS AND METHODS
RESULTS
DISCUSSION
Acknowledgements
References

ßCD Binds at a Single Site Near Residue 113 in the {alpha}HL Pore
Previous work demonstrated that ßCD binds within the lumen of the {alpha}HL pore (Gu et al. 1999 ), where it reduces the single-channel conductance (Gu et al. 1999 ; Gu and Bayley 2000 ), alters ion selectivity (Gu et al. 2000 ) and acts as a blocker site for various small organic molecules (Gu et al. 1999 , Gu et al. 2001 ). Mutagenesis experiments implied that ßCD binds in the vicinity of residue 113. The seven Met-113 side chains in the WT pore project into the lumen of the transmembrane ß barrel near its cis end (Fig 1 A). When the Met-113 residues are replaced with Asn, in M113N, ßCD binds 4.4 x 104 times more tightly (Gu et al. 1999 , Gu et al. 2001 ).

The present work, in which all possible natural amino acid substitutions at position 113 were examined and six were found to bind ßCD tightly (Table 1 and Fig 4), supports the idea that ßCD binds at or near residue 113. The variation in the affinity of ßCD (Kd) for the mutants at 113 spans a range of about five orders of magnitude, far greater than the roughly two orders observed for ßCD and WT-{alpha}HL when pH and transmembrane potential were varied (Gu and Bayley 2000 ). Further, mutations at positions in the ß barrel other than 113, notably 135 and 139, have little effect on the interaction with ßCD. WT-{alpha}HL (Leu-135, Asn-139), N139Q(WT) (Leu-135, Gln-139), and L135N(RL2) (Asn-135, Gln-139) all share Met-113 and all bind ßCD weakly. M113N(WT), M113N/N139Q(WT), and M113N/L135N(RL2) all share Asn-113 and all bind ßCD strongly. Similarly, the three pores with Met-113, show conductance blockades by ßCD in the range of 37–40%, whereas the three pores with Asn-113 have blockades in the range of 43–46%.

We cannot rule out the possibility that ßCD binds at a site removed from position 113 and that mutations at 113 induce conformational changes that affect the binding site. Indeed, the existence of two distinct classes of binding mutants (see next section) is suggestive of low and high affinity states. However, because a second cyclodextrin-binding site can be engineered nearer the trans entrance of the ß barrel (Gu et al. 2001 ) and because ßCD cannot bind from the cis side, there are few options for the location of the site under consideration here other than near, but not necessarily in contact with, residue 113. Obviously, additional support for the location of the ßCD site through structural studies would be most welcome.

The exact location aside, several arguments suggest that ßCD binds at a single site, rather than multiple sites, within the lumen of each {alpha}HL mutant. First, for each mutant, there is only one major conductance state that can be assigned to {alpha}HL•ßCD. Second, there is no additional noise associated with the {alpha}HL•ßCD state, suggesting that the rather rigid ßCD molecule is firmly held at the binding site (Fig 2). Third, for each mutant, whether binding is strong or weak, the dwell time histograms for {tau}off and {tau}on can be fitted to single exponentials (Supplemental Material), which is consistent with simple bimolecular kinetics for the interaction of ßCD with {alpha}HL.

The Kd Values for ßCD and Met-113 Replacement Mutants Fall into Two Major Classes
The Kd values for ßCD and the Met-113 mutants of {alpha}HL fall into two classes: class 1, tight binding mutants, mean Kd = 6.8 x 10-7 M; and class 2, weak binding mutants, mean Kd = 4.1 x 10-3 M. The major determinant of whether a mutant falls into class 1 or class 2 is the dissociation rate constant (koff), values of which fall into the same two classes (Fig 4). The association rate constant (kon) is hardly changed by mutagenesis (Fig 4), suggesting that it may reflect transfer of ßCD through the entrance to the pore, which would be expected to be little affected by mutations at position 113. The nature of the mutations that cause tight binding do not fall into an easily recognizable group. However, patterns do appear when the data are examined in detail. For example, whereas M113N and M113D fall into class 1, mutants with homologous substitutions, M113Q and M113E, are in class 2. The one mutant in class 1 with a nonaromatic hydrophobic side chain, M113V, is also the smallest of its kind. Mutants with aromatic substitutions, except the bulky tryptophan, are tight binding. As mentioned above, the existence of two distinct classes of binding mutants is suggestive of low and high affinity states of {alpha}HL for ßCD. However, this would require a conformational change that is not readily detected in measurements of primary electrical properties, such as single-channel conductance (Fig 5).

The Interaction of ßCD with {alpha}HL Is Voltage Dependent and Correlated with the Charge Selectivity of the Mutant Pores
The interaction of ßCD with the WT-{alpha}HL pore is voltage- and pH-dependent (Gu and Bayley 2000 ). At low pH values, ßCD (trans) binds more tightly at negative potentials; at high pH values, ßCD binds more tightly at positive potentials. Because ßCD is a neutral molecule, Woodhull's mechanism for the voltage-dependent binding of a charged blocker was ruled out. Further, because the dissociation rate constant of ßCD from its binding site varies continuously with voltage, a mechanism involving the voltage-dependent interconversion of two different states was also discounted. Instead, a continuous change in the free energy of {alpha}HL (and/or {alpha}HL•ßCD) as a function of the membrane potential was postulated (Gu and Bayley 2000 ). Voltage-dependent block by neutral molecules has been observed in other systems, but it has not been investigated in detail (Bezrukov et al. 2000 ).

Like the manipulation of pH, mutagenesis provides another way to change the charge distribution in a protein. In this work, we found that negative substituents at position 113 favor the binding of ßCD at positive potentials, whereas positive substituents favor binding at negative potentials (Table 2). It is possible that the charge status of residue 113 can determine the affinity of the pore for ßCD in the same way that unidentified charged groups do in WT-{alpha}HL (Gu and Bayley 2000 ). Therefore, an effect of membrane potential on the structure of the pore, as proposed for the pH-dependent properties of WT-{alpha}HL (Gu and Bayley 2000 ), is also a reasonable explanation for the effects of mutagenesis.

However, the mutagenesis experiments prompt a second possible explanation. The affinity of ßCD is correlated with the charge selectivity of the pore, which in {alpha}HL is modulated in a predictable manner by the charge at position 113 (Table 2). For example, M113K with seven more positive charges than WT-{alpha}HL, is more anion selective than WT-{alpha}HL, whereas M113E with seven additional negative charges is cation selective. The anion-selective mutants bind ßCD, applied from the trans chamber, more tightly at positive applied potentials, whereas the cation-selective mutants bind ßCD more tightly at negative potentials (Fig 7). Just as in the case of pH (Gu and Bayley 2000 ), both kon and koff are affected (Fig 7). Consistent with these data, pH also alters the charge selectivity of the {alpha}HL pore (Table II; Krasilnikov et al. 1997 ). For example, {alpha}HL is weakly anion selective at pH 7.5 (PK+/PCl- = 0.77), more anion selective at pH 5.0 (PK+/PCl- = 0.34), and cation selective at pH 11.0 (PK+/PCl- = 2.2). Again, the voltage dependence of the affinity for ßCD is correlated with charge selectivity (Table 2).

These results show that ßCD binding is favored when association occurs in the direction of the net movement of ions and dissociation occurs against the net ion flow. For example, in the case of M113E, a cation-selective channel, net ion flow is from trans to cis in a positive-applied potential, and ßCD binding from the trans side of the membrane is enhanced under these conditions. Thus, the second possible explanation is that the ßCD molecules move into the pore carried by water flow induced by ion movement—an electroosmotic effect (Katchalsky and Curran 1965 ). Because the main barrier to reaching the binding site is entry into the ß barrel, it is likely that the water flow would cause accumulation or depletion of the ßCD at the trans entrance compared with the bulk concentration and, hence, promote or impede binding.

Water flow caused by electroosmosis is considerable. Applying various simplifications, water flow in a charge-selective pore is given by (Katchalsky and Curran 1965 ):

where D is the diffusion coefficient of the mobile ion; r is the radius of pore; {eta} is the viscosity of water; I is the current flow through pore; and R, T, and F have their usual meanings. At a current of 20 pA, J ~109 s-1. Because 20 pA in a charge-selective pore corresponds to the transport of 1.25 x 108 ions s-1, approximately eight water molecules are transported for each ion that moves through the pore. This is in keeping with experimental findings. For example, in the case of a single cation-selective nafion-filled pore in mica of ~50 µm in radius, ~10 water molecules flow per Na+ ion transported (Bath et al. 1998 , Bath et al. 2000 ). In the same experiments, neutral organic molecules were shown to move with the water. For example, in the case of 0.2 M hydroquinone, about one molecule was transported for every 1,000 water molecules. For a 40-µM solute, one molecule would move with every 5 x 106 waters. Therefore, if the water flow were ~109 s-1 (see above), this would result in an appreciable concentration of ßCD at the trans entrance.

To further examine the likelihood of an electroosmotic effect, the mutant E111N/K147N was studied. E111N/K147N is a weakly cation-selective pore (PK+/PCl- = 1.20). As expected, ßCD applied from the trans chamber binds to E111N/K147N more strongly at positive than at negative potentials (ß1/K = 0.10). In most cases, ßCD binds to {alpha}HL pores only from the trans side of the bilayer. However, in the case of E111N/K147N, ßCD can also bind from the cis side and in this case, the voltage dependency was reversed (ß1/K = -0.15). Although this result is available only for a weakly ion-selective pore, it is in keeping with an electroosmotic mechanism, as opposed to an effect of voltage on the protein. The effects of voltage are also independent of whether a mutant is in class 1 or class 2, and this also favors an electroosmotic mechanism. Otherwise, the effects of voltage would have to be similar on both the conformation of class 1 and that of class 2.

In conclusion, we have shown that long residence times for the noncovalent molecular adapter ßCD within the lumen of the {alpha}HL pore can be achieved at neutral pH by mutagenesis at position 113. For example, the mean residence time ({tau}off) for ßCD bound to M113N at -40 mV is 27 s. In addition to the greatly enhanced affinity of the class 1 mutants, which most likely arises from alterations of the binding site for ßCD in the vicinity of residue 113, the binding can be further enhanced by a voltage-dependent mechanism that may originate in an electroosmotic effect. The latter explanation requires verification through more detailed study. The improved residence times for ßCD will contribute to our ability to build nanostructures from the {alpha}HL pore (Gu et al. 2001 ) and to use it as a component of stochastic sensors (Bayley and Cremer 2001 ).


  Footnotes

The online version of this article contains supplemental material.
* Abbreviations used in this paper: {alpha}HL, staphylococcal {alpha}-hemolysin; ßCD, ß-cyclodextrin.


  Acknowledgements
Top
Abstract
INTRODUCTION
MATERIALS AND METHODS
RESULTS
DISCUSSION
Acknowledgements
References

We thank Orit Braha, Alan Finkelstein, Chuck Martin, Liviu Movileanu, and the reviewers for their remarks.

This paper was supported by the US Department of Energy, the National Institutes of Health, the Office of Naval Research (MURI 1999), and the Texas Advanced Technology Program.

Submitted: 6 July 2001
Revised: 7 September 2001
Accepted: 11 September 2001


  References
Top
Abstract
INTRODUCTION
MATERIALS AND METHODS
RESULTS
DISCUSSION
Acknowledgements
References

    Bath, B.R., Lee, R.D., White, H.S., and Scott, E.R. 1998. Imaging molecular transport in porous membranes. Observation and analysis of electroosmotic flow in individual pores using the scanning electrochemical microscope. Anal. Chem. 70:1047-1058.

    Bath, B.R., White, H.S., and Scott, E.R. 2000. Electrically facilitated molecular transport. Analysis of the relative contributions of diffusion, migration, and electroosmosis to solute transport in an ion-exchange membrane. Anal. Chem. 72:433-442[Medline].

    Bayley, H. 1995. Pore-forming proteins with built-in triggers and switches. Bioorg. Chem. 23:340-345.

    Bayley, H. 1997. Building a door into cells. Sci. Am. 277:62-67[Medline].

    Bayley, H. 1999. Designed membrane channels and pores. Curr. Opin. Biotechnol. 10:94-103[Medline].

    Bayley, H., Braha, O., and Gu, L.-Q. 2000. Stochastic sensing with protein pores. Adv. Mater. 12:139-142.

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

    Bezrukov, S.M., Kullman, L., and Winterhalter, M. 2000. Probing sugar translocation through maltoporin at the single channel level. FEBS Lett. 476:224-228[Medline].

    Bhakdi, S., Füssle, R., and Tranum-Jensen, J. 1981. Staphylococcal {alpha}-toxin: oligomerization of hydrophilic monomers to form amphiphilic hexamers induced through contact with deoxycholate micelles. Proc. Natl. Acad. Sci. USA. 78:5475-5479[Medline].

    Braha, O., Walker, B., Cheley, S., Kasianowicz, J.J., Song, L., Gouaux, J.E., and Bayley, H. 1997. Designed protein pores as components for biosensors. Chem. Biol. 4:497-505[Medline].

    Braha, O., Gu, L.-Q., Zhou, L., Lu, X., Cheley, S., and Bayley, H. 2000. Simultaneous stochastic sensing of divalent metal ions. Nat. Biotechnol. 17:1005-1007.

    Cheley, S., Malghani, M.S., Song, L., Hobaugh, M., Gouaux, J.E., Yang, J., and Bayley, H. 1997. Spontaneous oligomerization of a staphylococcal a-hemolysin conformationally constrained by removal of residues that form the transmembrane b barrel. Protein Eng. 10:1433-1443[Abstract].

    Cheley, S., Braha, O., Lu, X., Conlan, S., and Bayley, H. 1999. A functional protein pore with a "retro" transmembrane domain. Protein Sci. 8:1257-1267[Abstract].

    Creighton, T.E. 1993. Proteins: Structures and Molecular Properties. New York, Freeman, , pp. 507 p.

    D'Souza, V.T., and Lipkowitz, K.B. 1998. Cyclodextrins. Chem. Rev. 98:1741-2076.

    Eroglu, A., Russo, M.J., Bieganski, R., Fowler, A., Cheley, S., Bayley, H., and Toner, M. 2000. Intracellular trehalose improves the survival of cryopreserved mammalian cells. Nat. Biotechnol. 18:163-167[Medline].

    Füssle, R., Bhakdi, S., Sziegoleit, A., Tranum-Jensen, J., Kranz, T., and Wellensiek, H.-J. 1981. On the mechanism of membrane damage by Staphylococcus aureus {alpha}-toxin. J. Cell Biol. 91:83-94[Abstract].

    Gu, L.-Q., and Bayley, H. 2000. Interaction of the non-covalent molecular adapter, ß-cyclodextrin, with the staphylococcal {alpha}-hemolysin pore. Biophys. J. 79:1967-1975[Abstract/Full Text].

    Gu, L.-Q., Braha, O., Conlan, S., Cheley, S., and Bayley, H. 1999. Stochastic sensing of organic analytes by a pore-forming protein containing a molecular adapter. Nature. 398:686-690[Medline].

    Gu, L.-Q., Dalla Serra, M., Vincent, J.B., Vigh, G., Cheley, S., Braha, O., and Bayley, H. 2000. Reversal of charge selectivity in transmembrane protein pores by using non-covalent molecular adapters. Proc. Natl. Acad. Sci. USA. 97:3959-3964[Abstract/Full Text].

    Gu, L.-Q., Cheley, S., and Bayley, H. 2001. Capture of a single molecule in a nanocavity. Science. 291:636-640[Abstract/Full Text].

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

    Howorka, S., Movileanu, L., Braha, O., and Bayley, H. 2001b. Kinetics of duplex formation for individual DNA strands within a single protein nanopore. Proc. Natl. Acad. Sci. USA. In press.

    Kasianowicz, J.J., Brandin, E., Branton, D., and DeAm, D.W. 1996. Characterization of individual polynucleotide molecules using a membrane channel. Proc. Natl. Acad. Sci. USA. 93:13770-13773[Abstract/Full Text].

    Katchalsky, A., and Curran, P.F. 1965. Nonequilibrium Thermodynamics in Biophysics. Cambridge, MA, Harvard University Press, , pp. 248 p.

    Krasilnikov, O.V., Capistrano, M.-F.P., Yuldasheva, L.N., and Nogueira, R.A. 1997. Influence of Cys-130 S. aureus alpha-toxin on planar lipid bilayer and erythrocyte membranes. J. Membr. Biol. 156:157-172[Medline].

    Montal, M., and Mueller, P. 1972. Formation of bimolecular membranes from lipid monolayers and study of their electrical properties. Proc. Natl. Acad. Sci. USA. 69:3561-3566[Medline].

    Movileanu, L., Howorka, S., Braha, O., and Bayley, H. 2000. Detecting protein analytes that modulate transmembrane movement of a polymer chain within a single protein pore. Nat. Biotechnol. 18:1091-1095[Medline].

    Otto-Bruc, A.E., Fariss, R.N., Van Hooser, J.P., and Palczewski, K. 1998. Phosphorylation of photolyzed rhodopsin is calcium-insensitive in retina permeabilized by {alpha}-toxin. Proc. Natl. Acad. Sci. USA. 95:15014-15019[Abstract/Full Text].

    Panchal, R.G., Cusack, E., Cheley, S., and Bayley, H. 1996. Tumor protease-activated, pore-forming toxins from a combinatorial library. Nat. Biotechnol. 14:852-856[Medline].

    Russo, M.J., Bayley, H., and Toner, M. 1997. Reversible permeabilization of plasma membranes with an engineered switchable pore. Nat. Biotechnol. 15:278-282[Medline].

    Sanchez-Quesada, J., Ghadiri, M.R., Bayley, H., and Braha, O. 2000. Cyclic peptides as molecular adapters for a pore-forming protein. J. Am. Chem. Soc. 122:11758-11766.

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

    Walker, B.J., and Bayley, H. 1994. A pore-forming protein with a protease-activated trigger. Prot. Eng. 7:91-97[Abstract].

    Walker, B.J., Krishnasastry, M., Zorn, L., Kasianowicz, J.J., and Bayley, H. 1992. Functional expression of the {alpha}-hemolysin of Staphylococcus aureus in intact Escherichia coli and in cell lysates. J. Biol. Chem. 267:10902-10909[Abstract].