Structure of a pore-blocking toxin in complex with a eukaryotic voltage-dependent K+ channel
Abstract
Pore-blocking toxins inhibit voltage-dependent K+ channels (Kv channels) by plugging the ion-conduction pathway. We have solved the crystal structure of paddle chimera, a Kv channel in complex with charybdotoxin (CTX), a pore-blocking toxin. The toxin binds to the extracellular pore entryway without producing discernable alteration of the selectivity filter structure and is oriented to project its Lys27 into the pore. The most extracellular K+ binding site (S1) is devoid of K+ electron-density when wild-type CTX is bound, but K+ density is present to some extent in a Lys27Met mutant. In crystals with Cs+ replacing K+, S1 electron-density is present even in the presence of Lys27, a finding compatible with the differential effects of Cs+ vs K+ on CTX affinity for the channel. Together, these results show that CTX binds to a K+ channel in a lock and key manner and interacts directly with conducting ions inside the selectivity filter.
https://doi.org/10.7554/eLife.00594.001eLife digest
The deadly toxins produced by many creatures, including spiders, snakes, and scorpions, work by blocking the ion channels that are essential for the normal operation of many different types of cells. Ion channels are proteins and, as their name suggests, they allow ions—usually sodium, potassium, or calcium ions—to move in and out of cells. They are especially important for cells that generate or respond to electrical signals, such as neurons and the cells in heart muscle.
Ion channels are located in the lipid membranes that surround all cells, and the ions enter or leave the cell via a pore that runs through the channel protein. They can be opened and closed (or ‘gated’) in different ways: some ion channels open and close in response to voltages, whereas others are gated by biomolecules, such as neurotransmitters, that bind to them.
Now, Banerjee et al. have used x-ray crystallography to study the structure of the complex that is formed when charybdotoxin (CTX), a toxin that is found in scorpion venom, blocks a voltage-gated potassium channel. Previous studies have shown that CTX binds to the channel on the extracellular side of the pore. Banerjee et al. show that the toxin fits into the entrance to the channel like a key into a lock, which means the toxin is preformed to fit the shape of the channel.
The potassium ion channel is made up of four subunits, and the pore contains four ion-binding sites that form a ‘selectivity filter’: it is this filter that ensures that only potassium ions can pass through the channel when it is open. When CTX binds to the channel, a lysine residue poised at a critical position on the toxin is so close to the outermost ion-binding site that it prevents potassium ions binding to the site. The structure determined by Banerjee et al. explains many previous findings, including the fact that ions entering the pore from inside the cell can disrupt the binding between the toxin and the ion channel protein. It remains to be seen if the toxins that target the pore of other types of ion channels work in the same way.
https://doi.org/10.7554/eLife.00594.002Introduction
Poisonous animals such as tarantula spiders, green mamba snakes, or the deathstalker scorpion rely on their venom for efficient defense and precapture strategies. These venoms are stored in dedicated glands and are rapidly delivered through specialized apparatus via subcutaneous, intramuscular, or intravenous routes. What is the underlying cause of toxicity of these venoms? The principal toxic components in all of them are peptidic in nature. These venoms usually contain libraries of hundreds of peptide-based toxins that together encompass a high degree of stereochemical diversity (Han et al., 2008; Liang, 2008; Rodriguez de la Vega et al., 2010). Only a small fraction of these molecules, however, have been pharmacologically characterized thus far. The targets of these toxins are typically a variety of ion channels—voltage-gated Na+(Nav), K+(Kv), and Ca2+(Cav) channels, and cell-surface ‘receptor’ ion channels, such as the nicotinic acetylcholine (Ach) receptor (Billen et al., 2008; King et al., 2008; Mouhat et al., 2008; Kasheverov et al., 2009). The remarkable molecular diversity of these toxins is borne out by the fact that multiple different toxins can target different components of the same ion channel/receptor. However, the end result is alteration of the normal physiology of the ion channel/receptor, thereby eliciting the desired reaction of the venom.
Potassium channels (K+ channels), a large and diverse class of ion channels, are targets of a large number of toxins that have been characterized up to now (Carbone et al., 1982; Miller et al., 1985; Galvez et al., 1990; Garcia et al., 1994; Swartz and MacKinnon, 1995). Most K+ channels are tetrameric in architecture—four pore domains together form an ion-conduction pathway through the membrane (Figure 1A; MacKinnon, 1991; Doyle et al., 1998). In addition, in the voltage-gated family of K+ channels (Kv channels), each channel monomer has attached onto the N-terminal end of the pore domain a transmembrane voltage sensor domain that senses the transmembrane voltage difference (Figure 1A; Papazian et al., 1987; Jiang et al., 2003). Kv channels are targeted by toxins primarily at two distinct sites—pore-blocking toxins that bind at the extracellular mouth of pore domains and gating-modifier toxins that bind to voltage sensor domains (Figure 1A; MacKinnon et al., 1990; Goldstein et al., 1994; Swartz and MacKinnon, 1997).
Scorpion venom specifically has been an abundant source of pore-blocking toxins for K+ channels. These are small peptides, typically ranging from 30 to 40 residues in length, held together by three or four disulfide bonds in a rigid architecture (Figure 1B; Bontems et al., 1991; Johnson and Sugg, 1992; Fernandez et al., 1994). Pore-blocking toxins have profoundly impacted research in the K+ channel field primarily in two ways. First, they have enabled purification of specific novel K+ channels such as the BK channel, a Ca2+ and voltage-gated K+ channel (Garcia et al., 1997). Second, they provided our first knowledge about channel subunit stoichiometry and the shape of the extracellular K+ pore entryway at a time when no three-dimensional structure was available for any ion channel (MacKinnon, 1991; Goldstein et al., 1994; Gross et al., 1994; Stampe et al., 1994; Hidalgo and MacKinnon, 1995; Naranjo and Miller, 1996; Ranganathan et al., 1996).
Charybdotoxin (CTX; Figure 1B), a pore-blocking toxin for K+ channels, is a 37-residue peptide isolated from the venom of the scorpion Leiurus quinquestriatus (Miller et al., 1985). Early experiments with CTX inhibition of the BK channel revealed that CTX binds to the extracellular surface of the channel with a 1:1 channel:toxin stoichiometry, that both the open and closed states of the channel are competent for toxin binding, and that electrostatic interactions play an important role in enhancing the toxin’s affinity (Anderson et al., 1988). Furthermore, CTX affinity was found to be voltage dependent, a property later shown to result from the destabilization of the toxin-channel complex by permeant ions entering from the intracellular side (an effect called ‘trans-enhanced dissociation’; MacKinnon and Miller, 1988). Ions that were unable to traverse the ion conduction pathway did not elicit trans-enhanced dissociation. These observations led to a hypothesis that CTX physically occludes the ion-conduction pathway, and in doing so brings a positive charge on CTX close to a K+ ion-binding site near the extracellular side (MacKinnon and Miller, 1988). The positive charge was later identified as Lys27, a residue that is conserved in all members of the CTX-like toxin family (Figure 2A; Park and Miller, 1992; Goldstein and Miller, 1993). Studies with other members of the CTX toxin family, most extensively, Agitoxin2 (AgTx2), supported the conclusion that they function in a manner similar to CTX (Garcia et al., 1994; Krezel et al., 1995; Hidalgo and MacKinnon, 1995; Ranganathan et al., 1996). Most notably, conservation of the toxin shape and the functionally important lysine suggested that they all bind with a similar orientation on the K+ channel and inhibit through a common mechanism, whereby a lysine amino group functions as a K+ ion mimic to block the pore (Miller, 1995; Figure 1B).
Double-mutant cycle studies between toxins (CTX and AgTx2) and the Shaker K+ channel provided numerous pairwise restraints for mapping the extracellular-facing pore surface (Goldstein et al., 1994; Gross et al., 1994; Stampe et al., 1994; Hidalgo and MacKinnon, 1995; Naranjo and Miller, 1996; Ranganathan et al., 1996). NMR-derived models using the KcsA K+ channel also provided valuable structural data (Takeuchi et al., 2003; Yu et al., 2005). However, models derived from the double-mutant cycle and NMR data were largely silent as to the influence of toxin on the conducting ions. Here, we have used x-ray crystallography to determine the structure of a complex between CTX and the paddle chimera, a mutant of the Kv1.2 K+ channel from rat brain, with particular focus on the influence of toxin on the selectivity filter structure and distribution of ions in the pore (Figure 2B; Alabi et al., 2007; Long et al., 2007).
Results
Overall architecture of the toxin-channel complex
Electrophysiological studies of paddle chimera in planar lipid bilayers had revealed that CTX inhibits paddle chimera with high affinity (∼20 nM Kd; Tao and MacKinnon, 2008). We crystallized the complex of paddle chimera with CTX by mixing together separately purified preparations of the channel and the toxin, and setting up cocrystallization trials. The highest resolution data were obtained from the complex of paddle chimera with the selenomethionine derivative of CTX. We used this dataset to solve the structure of the toxin complex of paddle chimera to 2.5 Å resolution. The architecture of the paddle chimera channel typifies the family of eukaryotic Kv channels such as Shaker, with four pore domains together forming the ion-conduction pathway through the membrane and four voltage sensor domains surrounding the pore (Figure 1A; Long et al., 2007). The voltage-sensors are linked to the cytoplasmic T1 domains that form a cytosolic tetrameric interface. Each channel-forming α-subunit is associated with an accessory β-subunit on the cytoplasmic side. The toxin channel complex shares the same overall architecture, with the fourfold symmetry axis of the channel tetramer coinciding with the fourfold crystallographic symmetry axis. There are two molecules of the α- and β-heterodimeric complex in the asymmetric unit (Figure 3). We refer to these as molecule A (top) and molecule B (bottom). Consequently, the symmetry operations generate two distinct tetramers in the lattice. For the toxin-channel complex, an initial omit map without any toxin in the model clearly shows electron density corresponding to the toxin at the pore entryway (Figure 4A). However, the electron density for the toxin bound to molecule A, henceforth referred to as toxin A, was clearer and so we used this map for building a model of the toxin-channel complex.
Building a model for the toxin
The consequences of an asymmetric toxin, binding to a fourfold symmetric channel, are that the toxin can bind to the tetrameric channel in four statistically distinguishable but structurally and energetically equivalent orientations. In this case, these orientations are related by the fourfold symmetry axis. Thus, the observed electron density map is a superposition of the electron densities for four such individual orientations of the toxin (Figure 4A). Moreover, since in the absence of any external constraints, each orientation is populated with one-fourth occupancy, the electron density for each orientation is inherently weak. Certain secondary structural features in the map were discernible and allowed approximate placement of a CTX molecule, whose structure was determined previously using NMR (Figures 1B, 4A; Bontems et al., 1991).
It was apparent from the outset that this initial model would require additional data to achieve a reasonable level of accuracy. To obtain the additional data, we incorporated electron-dense marker atoms individually, at several sites on the toxin and purified each heavy atom-modified toxin either by peptide synthesis followed by refolding and purification or by overexpression in Escherichia coli and following literature procedures (Park et al., 1991). We used three different heavy atom markers—replacement of a disulfide by a diselenide, 4-iodophenylalanine, and selenomethionine (Figure 4B). We tested each of these derivatives in a planar bilayer system, and they efficiently blocked the paddle chimera channel (Figure 4—figure supplement 1). We then crystallized each toxin derivative with paddle chimera and collected single-wavelength anomalous diffraction data at an appropriate wavelength for each derivative. For each derivative, an anomalous difference electron density map showed four heavy-atom peaks, corresponding to the four orientations of the toxin (Figure 5A). We collected datasets from crystals with each derivative at a distinct site on the toxin, a total of three datasets (Table 1). We next used our highest resolution dataset (of all the different toxin derivative complexes) and roughly placed the toxin in the omit electron density map (Figure 4A). We used this approximate initial placement of the toxin to determine which one of the four symmetry-equivalent peaks in the anomalous difference map (for each marker) corresponded to which orientation of our initial toxin placement. This provided a set of three experimental constraints corresponding to the three individual heavy atom peak positions (for three different markers; Figure 5B). The coordinates of the corresponding three sites from the known NMR structure of the toxin (Bontems et al., 1991) were used as three reference constraints. We then used RMSD-based superposition to minimize the sum of distances between the peaks in the map and the predicted positions on the toxin (Figure 5C; final RMSD 1.6 Å). This procedure yielded a constrained placement of the toxin, which was subsequently refined by rigid-body, coordinate-based B-factor refinement in CNS with manual adjustments, where appropriate (Figure 5D).
As noted previously, the density for toxin A in the omit map is better defined than for toxin B, and thus we chose toxin A density for initial placement and for building and refining the model of the toxin. However, when we superimposed the toxin-bound channel A onto channel B using the pore domains of channels A and B for the superposition, we observed a very reasonable model for toxin B, which agrees well with the omit electron density map for toxin B and places the side-chain of Met29 very near to one of the four experimentally observed heavy atom peaks for toxin B in the selenomethionine dataset (Figure 5—figure supplement 1). Since channels A and B are independent molecules in the asymmetric unit and only the heavy atom peak positions for toxin A were used for the initial placement of the toxin, this provided further validation of our model and boosted our confidence in the placement of the toxin. We did not include a model for toxin B in our final model since the overall electron density is not as well defined as for toxin A and building toxin B caused a small increase in Rfree.
Structural features of the toxin-channel complex
Superposition of the channel in the toxin complex onto the channel in the toxin-free structure (Long et al., 2007) shows that the channel undergoes no discernible structural changes (Figure 6A,F); RMSD 0.33 Å, residues 321–414, main chain atoms, molecule A; RMSD 0.16 Å, residues 321–414, main chain atoms, molecule B. This is consistent with the idea that the toxin fits into the mouth of the channel in a lock and key manner. The oblate-shaped toxin binds asymmetrically to the mouth of the pore domain of the channel such that the wider end is closer to the symmetry axis of the channel than the tapered end (Figure 6B). The helical part of the toxin molecule faces away from the channel and the edge of the toxin formed by the residues 25–29 on a β-strand faces toward the channel. The inherent architecture of this class of toxins is such that the three disulfide bonds that hold the folded toxin together are the main buried components in the structure and most of the side chains are displayed on the surface of the toxin. These side chains are in a position to engage into a number of different kinds of interactions with the channel molecule (Figures 7A,B). Closer to the fourfold symmetry axis of the channel, the aromatic ring of Tyr36 is positioned to pack simultaneously against Asp375 and Val377 of one subunit and Met29 is able to pack against Asp375 of an adjacent subunit (Figure 7A,B). Closer to the periphery, the side chain of Arg25 is within close proximity of Gln353, and the peptide backbone near Thr8-Thr9 is held against Gln353 of another subunit (Figure 7A,B). There are also residues that should be involved in long-range electrostatic interactions, that is, Arg25 is within 5.5 Å of Asp359. The guanidium headgroup of Arg25 could also, in principle, engage in electrostatic interactions with an ordered lipid that is present in the structure. Although the density of the lipid headgroup is not clear in this case, it is likely that the headgroup will be placed close enough for the Arg25 to make electrostatic contact with it. Arg34 is another residue that approaches the channel closely enough to make electrostatic interactions, that is, with the carbonyl oxygen of Asp375. Arg34 is also within H-bonding distance of Gln353 and within long-range electrostatic contact of Asp375. Asn30 is another residue that appears to approach close to Asp375 to enable a weak H-bonding interaction.
The most striking aspect of the toxin-channel complex concerns the distribution of K+ ions in the selectivity filter of the channel. In the structure of the paddle chimera, as well as other high-resolution K+ channel structures, the selectivity filter contains four distinct ion-binding sites or positions, S1 through S4, S1 being the most extracellular (Figure 6C; Zhou et al., 2001; Long et al., 2007; Nishida et al., 2007). From analyses of high-resolution diffraction data on KcsA, the prototypical K+ channel pore, it was inferred that during conduction, these ion-binding sites are occupied alternately in what are referred to as 1,3 and 2,4 configurations (Morais-Cabral et al., 2001). In accordance, the occupancy of each site was experimentally determined to be roughly 0.5 (Zhou and MacKinnon, 2003). In the toxin-channel complex, only sites S2 through S4 have discernible electron density (Figure 6D). The top ion-binding site appears empty. It is important to note here that this observation holds strictly true for both channel molecules in the asymmetric unit. We have several datasets for all the different derivative toxin complexes, and this observation holds true for all of them as well, in both channel molecules A and B, in the asymmetric unit.
Why is the distribution of ions dramatically different in the toxin-bound complex? The toxin binds at the mouth of the pore positioned in such a manner to project the side chain of Lys27 straight into the pore, allowing the amino group to approach the top of site S1 (Figure 6E and Figure 6—figure supplement 1). Thus, Lys27 is within the range to make hydrogen-bonding interactions with all four carbonyl oxygen atoms that would otherwise constitute the top-half layer of coordinating ligands for a K+ at site S1 (Figure 6E; Zhou et al., 2001). We suspect therefore two reasons why K+ is disfavored at site S1. First, there is electrostatic repulsion from the closely placed positively charged amino group, and second, carbonyl oxygen atoms that would otherwise constitute half of the coordination are not fully available for coordination with a K+ at site S1. The structure thus offers a simple rationale for the altered distribution of ions in the selectivity filter of the toxin-bound channel.
Structure of the toxin-channel complex in Cs+ and structure with a Lys27 mutant of CTX
The altered ion distribution informs us that the toxin interacts with ions inside the pore. This interaction is compatible with the electrophysiological observation that intracellular ions destabilize extracellular toxin—the trans-enhanced dissociation effect (MacKinnon and Miller, 1988; Park and Miller, 1992). To further correlate the ion distribution in the crystal with trans-enhanced dissociation in electrophysiology experiments, we crystallized the channel with a mutant toxin, Lys27Met. Miller and coworkers had shown that Lys27 mutants reduce toxin affinity and abolish trans-enhanced dissociation (Park and Miller, 1992; Goldstein et al., 1994). Lys27Met CTX indeed inhibits the paddle chimera channel with reduced affinity (∼630 nM) but still forms a complex at the high concentrations of channel (∼20 μM) and toxin (60–80 μM) present in crystallization trials (Figure 8A,B). Electron density is present at S1 in an ion omit map, compatible with the presence of an ion, although the density at S1 is weaker relative to the ion densities at S2–S4, as if the occupancy at S1 is reduced (Figure 8C). Thus, it appears that the ability of toxin to interact with ions in the pore (and thus alter the ion distribution) is directly connected to the ability of ions from the intracellular solution to destabilize toxin on the extracellular side through interaction with Lys27. A simple mechanistic explanation could be that K+ and toxin—via Lys27—compete for stabilizing interactions at S1 in the selectivity filter.
The above mechanistic proposal provides motive to wonder what happens when K+ is replaced with Cs+, because Cs+ in crystal structures of K+ channels binds at only three sites and with unusually high occupancy at S1 (Zhou and MacKinnon, 2003). A CTX complex with the paddle chimera channel in the presence of Cs+ shows that toxin is bound but that Cs+ adopts its expected distribution (i.e., similar to its distribution in KcsA in the absence of toxin) with an ion at S1 (Figure 8D,E). It is difficult to tell where in the electron density map the amino group of Lys27 resides, but it is clear that Cs+ competes effectively for site S1, despite the presence of CTX. The mechanism of competition put forth above predicts that CTX should not bind with high affinity in the presence of Cs+. As shown, this prediction holds: CTX inhibits with a nearly 10-fold reduced affinity in the presence of Cs+ compared to K+ (Figure 8A).
Discussion
We report here the first x-ray structure of a K+ channel bound to a toxin. Many of the original studies on pore-blocking toxins for K+ channels were carried out with eukaryotic voltage-gated K+ channels, such as Shaker and Kv1.3, closely related in sequence to the mutant version of the eukaryotic Kv1.2 channel that we employed in our structural studies (MacKinnon and Miller, 1989; MacKinnon, 1991; Goldstein and Miller, 1992; Stampe et al., 1992, 1994; Goldstein et al., 1994; Gross et al., 1994; Aiyar et al., 1995, 1996; Hidalgo and MacKinnon, 1995; Gross and MacKinnon, 1996; Naini and Miller, 1996; Naranjo and Miller, 1996; Ranganathan et al., 1996; MacKinnon et al., 1998). Our first major finding is that we do not see discernible changes in the structure of the channel between the toxin-bound and the toxin-free paddle chimera structures (Figure 6A,F). This, we note, is in contrast to the solid-state NMR ‘structure’ of a KcsA mutant with kaliotoxin, where such rearrangements in the channel were proposed (Lange et al., 2006). In this NMR study, chemical shift changes induced by toxin were interpreted as resulting from dihedral angle changes (i.e., structural); however, such chemical shift changes could have other origins (i.e., electrostatic). Moreover, no channel-toxin distance restraints were included in the determination of this solid-state NMR ‘structure’ (Lange et al., 2006). In the crystal structure presented here, the good precomplex complementarity between the shape of the channel entryway and the shape of the toxin (i.e., the absence of channel structural change upon toxin binding) suggests a possible explanation for two prominent features of pore-blocking toxins. First, toxins can bind with relatively high affinity to their target channels because binding free energy is not ‘spent’ bringing about a protein conformational change. And second, single mutations in the ‘toxin receptor’ region of the channel can drastically alter the affinity for the channel by disrupting the good fit (Goldstein et al., 1994; Garcia et al., 1997).
The relatively static architecture of the K+ channel pore entryway undoubtedly reflects the requirement of a well-ordered selectivity filter structure to select K+ ions. In hindsight, this static toxin receptor on Kv channels lends credence to the idea that was originally proposed for using the pore-blocking toxins as ‘molecular slide calipers’ for gauging, at the resolution of mutagenesis, the shape of the protein surface on the extracellular part of the pore domain (Goldstein et al., 1994; Stampe et al., 1994; Hidalgo and MacKinnon, 1995; Ranganathan et al., 1996). A large body of data has emerged from these studies on mutagenesis-based electrophysiological measurements of toxin-channel interactions in the Kv channel family, mainly using the toxin-channel pairs Shaker–CTX, Shaker–AgTx2, and to a lesser extent, Kv1.3–CTX (Goldstein et al., 1994; Stocker and Miller, 1994; Hidalgo and MacKinnon, 1995; Aiyar et al., 1995; Naini and Miller, 1996; Naranjo and Miller, 1996; Ranganathan et al., 1996; Rauer et al., 2000). We have mapped these data onto our structure. Among these measurements, the data obtained with mutant cycle analyses are more reliable in gauging residue proximities on either side of the toxin-channel interface. The mapping has been done in the following manner: for the channel, we have used a sequence alignment (Figure 2B) to match the corresponding residue in Shaker or Kv1.3 onto paddle chimera. Data from the studies using CTX are shown in Figure 9A and those using AgTx2 in Figure 9B. For AgTx2, we used the known structure of AgTx2 and the guidelines in Krezel et al. (1995) to superimpose AgTx2 onto CTX in the structure of the paddle chimera–CTX complex. A few data points were derived from lysine scanning, which alters the length of the wild-type residue appreciably. For the lysine mutations, for the sake of representation, the corresponding residue in paddle chimera was mutated in silico (in Coot) to lysine, and the rotamer of lysine with the least distance between the corresponding toxin residue was chosen. Our structure is overall in excellent agreement with not only the CTX data but also with the AgTx2 data.
We wish to discuss two cases of proximity deduced by the mutant cycle analyses and how the structure reveals atomic insights into them. From analyzing mutants of AgTx2 and Shaker, MacKinnon and Ranganathan derived a coupling energy of >3 kT for the Gly10Val(AgTx2)–Phe425Gly(Shaker) pair and ∼1.5 kT for the Gly10Val(AgTx2)–Thr449Cys(Shaker) pair (Ranganathan et al., 1996). Phe425 and Thr449 in Shaker map onto Gln353 and Val377 in paddle chimera (Figure 2B). An inspection of the structure reveals clearly that Gly10 is close to Gln353 (Figure 10A), which is consistent with the high coupling energy. Intriguingly, Val377 is not within the first layer of surrounding residues contacted by Gly10. However, Val377 is within close proximity to Gln353 such that a mutation of Gly10 to valine would incur a clash of Gln353 with Val377 (Figure 10A). It is worthwhile noting that pairwise mutant cycle analysis per se does not distinguish between direct interactions and such interactions mediated by a third residue. However, the relative magnitude of the coupling energies, in hindsight, is consistent with such a mode of interaction.
Another case worth highlighting is the coupling energy, ∼1.5 kT, between the Lys27Met(AgTx2)–Tyr445Phe(Shaker) pair. Tyr445 in Shaker maps to Tyr373 in paddle chimera (Figure 2B; Ranganathan et al., 1996). Since mutation of Tyr to Phe incurs only a loss of the hydroxyl group, one likely interpretation of this would have been that Lys27 contacts the hydroxyl group of Tyr445. Intriguingly, Lys27 only contacts the backbone carbonyl oxygen of Tyr445, a point farthest from the hydroxyl group (Figure 10B). However, since this is a tightly packed part of the protein structure with the hydroxyl group being a part of an interaction network, mutation of the hydroxyl is felt at the backbone carbonyl by the toxin. MacKinnon and Ranganathan also found that the interaction between these two residues is dependent on the K+ concentration. This is also consistent with Lys27 inserting into the selectivity filter and displacing K+ from site S1.
The structure also shows the chemical rationale behind an experimental observation that is decades old concerning the role of a highly conserved lysine in voltage-dependent block by CTX and other pore-blocking toxins in this family. In studies with the BK channel, MacKinnon and Miller (1988) first observed that permeant ions coming through the channel from the intracellular side enhance the dissociation of a toxin bound to its receptor on the extracellular side. This trans-enhanced dissociation effect of permeant ions on CTX block of K+ channels was subsequently confirmed for Shaker as well (Goldstein and Miller, 1993). In order to explain this phenomenon, MacKinnon and Miller put forward the hypothesis that the toxin binds at a site close to the ion permeation pathway and places a positive charge close to one of the ion-binding sites in the channel. Additionally, it was also observed that this voltage-dependent and permeant ion–dependent block was nearly completely abolished when Lys27 was mutated to a neutral Asn or Gln residue (Park and Miller, 1992). Mutation of no other residue on the toxin had the same effect. This implied that the electrostatic interaction between the K+ in the ion-binding site and the toxin was wholly mediated by this single Lys residue, and the structure of the toxin-paddle chimera complex shows exactly why that is so. In the structure of the WT CTX-K+ complex, the top ion-binding site S1 has no density corresponding to K+ ions (Figure 6C) and the positively charged amino group is positioned to form contacts with S1 instead (i.e., making hydrogen-bonding interactions with the carbonyl oxygen atoms comprising the top layer of coordinating ligands at S1; Figure 6D). From the structure, it is clear that only the side chain of Lys27 approaches close enough to the S1 site to exert such an effect (Figure 6—figure supplement 1). It is worth noting here that the long linear alkyl chain and the tetrahedral disposition of hydrogens around the terminal nitrogen of Lysine make it uniquely and chemically suited for making maximal contacts with the top layer of coordinating carbonyl oxygen atoms at S1. Because this region is buried within the protein, this is likely to be a highly stabilizing interaction. Even a conservative mutation of Lys27 to Arginine causes ∼1000-fold destabilization of the toxin-channel complex in Shaker (Goldstein et al., 1994).
The CTX family of K+ channel toxins has remarkable diversity in sequence and binding affinities to specific channel subtypes. It is interesting to note, in light of the present structure, the basic mechanistic principles underlying the mode of action of these toxins. In a strikingly simple but effective strategy, the toxins target a functional aspect common to all the K+ channels—ion conduction. At the structural level, this is executed by presenting the amino group of Lys27, a highly conserved residue (Figure 2A and Figure 6—figure supplement 1), to the top ion-binding site in the selectivity filter. This acts as a tethered surrogate cation and effectively plugs the ion-conduction pathway. How is this surrogate cation brought to this site in the first place? The structures of these miniproteins are held together by rigid disulfide bonds, that are highly conserved in this and other families of small peptide-based toxins. This provides a relatively rigid scaffold. Through evolutionary sampling of intervening less-conserved residues, individual toxins have gained the ability to engage in specific interactions by long-range electrostatic interactions and a few subtype-specific close contacts, thereby ensuring efficient channel block. This is a remarkable example of combinatorial diversity in nature with the constraints of a rigid scaffold and a conserved mechanism to target one of the most important classes of sensory molecules in biology. Scorpion toxins have evolved to fit like a lock and key into the pore entryway of potassium channels, and disrupt ion conduction through presentation of a lysine amino group that competes with potassium in the selectivity filter.
Materials and methods
Molecular biology
Mutagenesis
All the mutants were generated using the QuikChange site-directed mutagenesis kit (Agilent, Santa clara, USA) and incorporation of the mutation(s) was verified by sequencing.
Protein expression and purification
Expression and purification of Kv2.1 paddle–Kv1.2 chimera channel
Request a detailed protocolKv2.1 paddle–Kv1.2 chimera channel was expressed and purified as described with slight modifications (Long et al., 2007). In brief, the channel was coexpressed with the rat β2-core gene in Pichia pastoris, extracted with DDM, and purified on a cobalt affinity column followed by gel filtration on a Superdex-200 column using a 1:1 mixture of Cymal6 and Cymal7 detergents. The paddle chimera-β2 complex thus purified was subsequently concentrated to 8–9 mg/ml (Centricon-100; Millipore, Billerica, MA).
Expression and purification of wild-type CTX, Ly27Met CTX, and selenomethionine derivative of CTX
Request a detailed protocolWild-type and Lys27Met CTX were expressed in BL21(DE3) E. coli as a T7 gene 9 fusion protein and purified as described (Park et al., 1991). Trypsin was used to cleave the inhibitor from the carrier at the Factor Xa site. The selenomethionine derivative was made using the same method except for modifications that were made in the cell growth conditions to prepare the selenomethionine derivative using the methionine biosynthesis inhibition method as described (van Duyne et al., 1993).
Synthesis and refolding of 4-iodophenylalanine and diselenide mutants of CTX
Materials
Fmoc protected amino acids, Wang resins, and 2-(1H-benzotriazol-1-yl)-1,1,3,3-tetramethyluronium hexafluorophosphate (HBTU) were purchased from Novabiochem (Switzerland). Fmoc-4-iodophenylalanine, Fmoc-(Se-p-methoxybenzyl)–selenocysteine, and Pyroglutamic acid were purchased from Anaspec (Freemont, CA). HPLC grade MeCN, dichloromethane (DCM), dimethylformamide (DMF), MeOH, and N-methylpyrrolidinone (NMP) were purchased from Fischer Scientific (Pittsburgh, PA). Trifluoroacetic acid (TFA) was purchased from Halocarbon (River Edge, NJ). Biotech grade piperidine was from Sigma-Aldrich (St Louis, MO). N,N-Diisopropylethylamine (DIPEA) was from Applied Biosystems (Foster City, CA). All other reagents were purchased from Sigma-Aldrich at the highest available purity.
4-Iodophenylalanine mutant
Request a detailed protocolThe linear precursor was synthesized using Fmoc-chemistry on Wang resin at a 0.25-mmol scale on a CEM liberty microwave-assisted peptide synthesizer (CEM, Matthews, NC) using standard instrument protocols except a milder coupling cycle was used for the coupling steps involving 4-iodophenylalanine. Following synthesis, the peptide was cleaved from the resin with reagent K for 4 hr at room temperature. The free peptide was then filtered to remove spent resin beads, precipitated with ice-cold diethyl ether, washed twice with ice-cold ether, and the remaining ether was removed under vacuum. The crude peptide was subsequently dissolved in water to a concentration of 0.5 mg/ml, and the pH of the solution was adjusted to 8.1 with N-methylmorpholine. Oxidative folding was performed by the addition of 0.5 mM GSH, and the peptide was allowed to oxidize in the presence of air. The oxidative folding was monitored using analytical scale RP-HPLC and ESI-MS. The cyclized peptide was purified by reverse-phase HPLC with a Vydac C-18 reverse-phase preparative column (Vydac 218TP1022) using a linear gradient of 14% solvent B to 40% solvent B over a 60-min period (solvent A—0.1% TFA in H2O; solvent B—0.1% TFA in 90% MeCN/H2O). The cyclized peak eluted as the major product and was further characterized by analytical reverse-phase HPLC using C18 reverse-phase analytical column (Vydac 218TP5415) and ESI-MS. From 100 mg of crude peptide, the typical yield of purified toxin was approximately 0.6–0.7 mg.
Diselenide mutant
Request a detailed protocolThe linear precursor of the diselenide mutant was synthesized on Wang resin at a 0.25-mmol scale on a CEM liberty microwave-assisted peptide synthesizer (CEM) except the two steps where SeCys was coupled onto the solid phase were performed manually to conserve reagent. The coupling cycles used on the CEM synthesizer all through the synthesis were milder than standard synthesizer cycles. Following a published protocol (Walewska et al., 2009), the peptide was cleaved from the resin using reagent K and DTNP, followed by thiolation in presence of DTT. The crude peptide was purified using reverse-phase HPLC with a Vydac C-18 reverse phase preparative column (Vydac 218TP1022) using a linear gradient of 14% solvent B to 40% solvent B over a 60-min period (solvent A—0.1% TFA in H2O; solvent B—0.1% TFA in 90% MeCN/H2O). It was subjected to cyclization as described above for the 4-iodophenylalanine mutant, and the cyclized product was again purified by reverse-phase HPLC with a Vydac C-18 reverse-phase preparative column (Vydac 218TP1022) using a linear gradient of 14% solvent B to 40% solvent B over a 60-min period (solvent A—0.1% TFA in H2O; solvent B—0.1% TFA in 90% MeCN/H2O). The product was further characterized by analytical reverse-phase HPLC using C18 reverse-phase analytical column (Vydac 218TP5415) and ESI-MS. From ∼120 mg of crude peptide, the typical yield of purified toxin was approximately 0.6–0.8 mg.
Crystallization of the toxin–paddle chimera complexes
Request a detailed protocolConcentrated toxin (typically approximately 2–3 mM in water) was supplemented with detergents and lipids so that the final concentration of detergents and lipids in the toxin were the same as in the protein. The concentrated paddle chimera-β2 complex at 8–9 mg/ml was mixed with toxin such that there was a threefold to fourfold molar excess of toxin over the tetrameric channel. Crystallization trials were set up by mixing the toxin-channel complex with the crystallization solution in a 2:1 ratio and supplemented with 10% vol/vol of 40 mM CHAPS. The complexes were crystallized using the hanging drop vapor diffusion method over reservoirs containing 0.1 ml crystallization solution at 20°C. The crystallization solution contained 28–32% PEG400 and 50 mM Tris–HCl, pH 8.8–9.1. Crystals appeared typically within 1–3 days and were directly frozen in liquid nitrogen after overnight equilibration against a reservoir solution containing 33% PEG400 and 50 mM Tris–HCl, pH 8.5.
Structure determination and building a model for CTX
Request a detailed protocolAll diffraction data including selenium and iodine anomalous diffraction data were collected at beamline X29 (Brookhaven NSLS), and images were processed with HKL2000 (Otwinowski and Minor, 1997). Data were further processed using the CCP4 suite (Dodson et al., 1997). The crystals were isomorphous to the paddle chimera channel, belonging to the P4212 space group. The paddle chimera structural model (PDB ID 2 R9R; Long et al., 2007) without ions was used as a starting model. Rigid body refinement in CNS (Brunger et al., 1998; Brunger, 2007; keeping the same test set as wild type) generated an acceptable model with R-free below 35%. Further rounds of iterative refinement in CNS and manual rebuilding using COOT (Emsley and Cowtan, 2004) generated a more complete model with residual R-free below ∼28%. Crystals with the selenomethionine CTX derivative diffracted to the highest resolution of all three heavy atom derivatives and was chosen for subsequent refinement and for building a model for the CTX.
At this stage, a weighted toxin-omit 2Fo − Fc electron density map clearly revealed the presence of the toxin at the mouth of the pore. However, the observed electron density map is a superposition of the electron densities for four individual orientations in which the toxin can bind to the tetrameric channel. Moreover, each orientation is only fractionally occupied, thus making the electron density for each orientation inherently weak. Certain secondary structural features in the map were discernible (Figure 4A) and allowed approximate placement of a CTX molecule using the available NMR structure of CTX (Bontems et al., 1991). However, it was obvious from the outset that this initial model would require additional data to achieve a reasonable level of accuracy. We used anomalous electron density maps from three different heavy atom derivatives of CTX to obtain further constraints that helped us in improving our initial placement of the toxin (see section “Building a model for the toxin” in main text for details). This improved placement was used as a starting point for crystallographic refinement, which was subsequently refined by rigid-body, coordinate-based B-factor refinement in CNS with manual adjustments, where appropriate. During refinement, the interactions between symmetry-related toxin molecules were prevented by using ‘igroup’ statements in CNS refinement scripts.
In particular, the asymmetric unit contains two α- and β-heterodimers (molecules A and B, see Figure 3). The crystallographic symmetry operation thus creates two independent tetrameric channels, each containing four α-monomers and four β-monomers. The model for the toxin was built for the toxin bound to channel A since the electron density for this toxin was much better defined. The asymmetric unit contained one orientation of the toxin modeled at 0.25 occupancy. Consequently, the fourfold symmetry operation generated the other three possible orientations for the toxin molecule to bind to the channel. The independence of molecules A and B in the asymmetric unit and anomalous electron density from the selenomethionine derivative of CTX provided additional validation for our model of the toxin (see section “Building a model for the toxin” in main text for details). Crystallographic data and refinement statistics are shown in data collection and refinement statistics table. Figures were made using PYMOL (www.pymol.org) and COOT (Emsley and Cowtan, 2004).
Structure determination and model building for the Lys27Met CTX–paddle chimera complex and CTX–paddle chimera complex in CsCl
Request a detailed protocolDiffraction data were collected at beamline X29 (Brookhaven NSLS), the data were processed, and the structure was solved as above. An initial omit map calculated using a model without any ions or toxin clearly revealed the presence of the toxin at the mouth of the pore in each case. However, in order to build a model for the toxin in each case, we started from the model of the toxin as built above and carried out subsequent rigid-body, coordinate-based B factor refinements in CNS (Brunger et al., 1998; Brunger, 2007) with slight manual adjustments of the model in COOT (Emsley and Cowtan, 2004) wherever deemed appropriate from prominent features in the electron density map.
Electrophysiological studies of toxin-channel interactions
Request a detailed protocolElectrophysiology of paddle chimera channels was performed essentially as described (Tao and MacKinnon, 2008), except in the experiments for determining Kd’s of paddle chimera with wild-type CTX, Lys27Met CTX, and wild-type CTX in the presence of CsCl, the planar bilayer membranes were painted with 20 mg/ml 1,2-diphytanoyl-sn-glycero-3-phosphocoline (DPhPC; Figure 6A). For experiments in CsCl, the KCl in the cis- and the trans-chambers was replaced by CsCl.
The membranes were held at −110 mV and repeatedly pulsed to +110 mV test voltage. CTX inhibition data measured with different concentrations of CTX or Lys27Met mutant of CTX were fit to the equation:
The modified form of the equation was used because even at the highest concentration of the toxin, we observed residual current representing ∼10% of the initial current that is due largely to the channels facing in the opposite direction. Statistical fit and determination of Kd values were done using the IgorPro (Wavemetrics, Portland, USA) software.
For experiments showing paddle chimera block by a single concentration of heavy atom derivatives of CTX (Figure 4—figure supplement 1B–D), the planar lipid bilayers were formed of 20 mg/ml 3:1 (wt:wt) 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine (POPE): 1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-(1′-rac-glycerol; POPG). For the experiment showing paddle chimera block by single concentration of wild-type CTX in Figure 4—figure supplement 1A, the planar lipid bilayer was formed of 20 mg/ml 1,2-diphytanoyl-sn-glycero-3-phosphocoline (DPhPC).
Data availability
-
Crystal structure of Kv1.2-2.1 paddle chimera channel in complex with CharybdotoxinID 4JTA. Publicly available at the RCSB Protein Data Bank (http://www.rcsb.org/pdb/).
-
Crystal structure of Kv1.2-2.1 paddle chimera channel in complex with Lys27Met mutant of CharybdotoxinID 4JTD. Publicly available at the RCSB Protein Data Bank (http://www.rcsb.org/pdb/).
-
Crystal structure of Kv1.2-2.1 paddle chimera channel in complex with Charybdotoxin in Cs+ID 4JTC. Publicly available at the RCSB Protein Data Bank (http://www.rcsb.org/pdb/).
-
Shaker family voltage dependent potassium channel (kv1.2-kv2.1 paddle chimera channel) in association with beta subunitID 2R9R. Publicly available at the RCSB Protein Data Bank (http://www.rcsb.org/pdb/).
-
Analysis of side-chain organization on a refined model of charybdotoxin: structural and functional implicationsID 2CRD. Publicly available at the RCSB Protein Data Bank (http://www.rcsb.org/pdb/).
References
-
Animal peptides targeting voltage-activated sodium channelsCurr Pharm Des 14:2492–2502.https://doi.org/10.2174/138161208785777423
-
Version 1.2 of the crystallography and NMR systemNat Protoc 2:2728–2733.https://doi.org/10.1038/nprot.2007.406
-
Crystallography & NMR system: a new software suite for macromolecular structure determinationActa Crystallogr D Biol Crystallogr 54:905–921.https://doi.org/10.1107/S0907444998003254
-
Coot: model-building tools for molecular graphicsActa Crystallogr D Biol Crystallogr 60:2126–2132.https://doi.org/10.1107/S0907444904019158
-
Purification and characterization of a unique, potent, peptidyl probe for the high conductance calcium-activated potassium channel from venom of the scorpion Buthus tamulusJ Biol Chem 265:11083–11090.
-
Pharmacology of potassium channelsAdv Pharmacol 39:425–471.https://doi.org/10.1016/S1054-3589(08)60078-2
-
Conus venoms—a rich source of peptide-based therapeuticsCurr Pharm Des 14:2462–2479.https://doi.org/10.2174/138161208785777469
-
Proteome and peptidome profiling of spider venomsExpert Rev Proteomics 5:731–746.https://doi.org/10.1586/14789450.5.5.731
-
Animal toxins acting on voltage-gated potassium channelsCurr Pharm Des 14:2503–2518.https://doi.org/10.2174/138161208785777441
-
Processing of X-ray diffraction data collected in oscillation modeMacromol Crystallogr Pt A 276:307–326.
-
Mining on scorpion venom biodiversityToxicon 56:1155–1161.https://doi.org/10.1016/j.toxicon.2009.11.010
-
Electrostatic distance geometry in a K+ channel vestibuleProc Natl Acad Sci USA 91:9509–9513.https://doi.org/10.1073/pnas.91.20.9509
Article and author information
Author details
Funding
Howard Hughes Medical Institute
- Roderick MacKinnon
National Institutes of Health (GM 43949)
- Roderick MacKinnon
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank Tom Muir (Princeton University) for access to CEM peptide synthesizer and Miquel Vila-Perello in the Muir lab for help and advice on peptide synthesis. We thank H Robinson at beamline X29 (National Synchrotron Light Source, Brookhaven National Laboratory) for assistance at the synchrotron; J Butterwick for comments on the manuscript. RM is an Investigator in the Howard Hughes Medical Institute. Supported by NIH GM43949 awarded to RM. The atomic coordinates and structure factors have been deposited to the Protein Data Bank with the accession numbers 4JTA (selenomethionine CTX with paddle chimera), 4JTD (Lys27Met CTX with paddle chimera), and 4JTC (CTX with paddle chimera in Cesium Chloride).
Copyright
© 2013, Banerjee et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 4,596
- views
-
- 189
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Structural Biology and Molecular Biophysics
Proteins occurring in significantly high concentrations in cellular environments (over 100 mg/ml) and functioning in crowded cytoplasm, often face the prodigious challenges of aggregation which are the pathological hallmark of aging and are critically responsible for a wide spectrum of rising human diseases. Here, we combine a joint-venture of complementary wet-lab experiment and molecular simulation to discern the potential ability of adenosine triphosphate (ATP) as solubilizer of protein aggregates. We show that ATP prevents both condensation of aggregation-prone intrinsically disordered protein Aβ40 and promotes dissolution of preformed aggregates. Computer simulation links ATP’s solubilizing role to its ability to modulate protein’s structural plasticity by unwinding protein conformation. We show that ATP is positioned as a superior biological solubilizer of protein aggregates over traditional chemical hydrotropes, potentially holding promises in therapeutic interventions in protein-aggregation-related diseases. Going beyond its conventional activity as energy currency, the amphiphilic nature of ATP enables its protein-specific interaction that would enhance ATP’s efficiency in cellular processes.
-
- Structural Biology and Molecular Biophysics
How the sequences of intrinsically disordered proteins (IDPs) code for functions is still an enigma. Dynamics, in particular residue-specific dynamics, holds crucial clues. Enormous efforts have been spent to characterize residue-specific dynamics of IDPs, mainly through NMR spin relaxation experiments. Here, we present a sequence-based method, SeqDYN, for predicting residue-specific backbone dynamics of IDPs. SeqDYN employs a mathematical model with 21 parameters: one is a correlation length and 20 are the contributions of the amino acids to slow dynamics. Training on a set of 45 IDPs reveals aromatic, Arg, and long-branched aliphatic amino acids as the most active in slow dynamics whereas Gly and short polar amino acids as the least active. SeqDYN predictions not only provide an accurate and insightful characterization of sequence-dependent IDP dynamics but may also serve as indicators in a host of biophysical processes, including the propensities of IDP sequences to undergo phase separation.