Email updates

Keep up to date with the latest news and content from BMC Structural Biology and BioMed Central.

Open Access Research article

Differential molecular information of maurotoxin peptide recognizing IKCa and Kv1.2 channels explored by computational simulation

Hong Yi12, Su Qiu1, Yingliang Wu1, Wenxin Li1 and Baoshan Wang2*

Author Affiliations

1 State Key Laboratory of Virology, College of Life Sciences, Wuhan University, Wuhan, 430072, PR China

2 College of Chemistry and Molecular Sciences, Wuhan University, Wuhan 430072, PR China

For all author emails, please log on.

BMC Structural Biology 2011, 11:3  doi:10.1186/1472-6807-11-3


The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1472-6807/11/3


Received:16 September 2010
Accepted:25 January 2011
Published:25 January 2011

© 2011 Yi et al; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

Scorpion toxins are invaluable tools for ion channel research and are potential drugs for human channelopathies. However, it is still an open task to determine the molecular basis underlying the diverse interactions between toxin peptides and ion channels. The inhibitory peptide Maurotoxin (MTX) recognized the distantly related IKCa and Kv1.2 channel with approximately the same potency and using the same functional residues, their differential binding mechanism remain elusive. In this study, we applied computational methods to explore the differential binding modes of MTX to Kv1.2 and IKCa channels, which would help to understand the diversity of channel-toxin interactions and accelerate the toxin-based drug design.

Results

A reasonably stable MTX-IKCa complex was obtained by combining various computational methods and by in-depth comparison with the previous model of the MTX-Kv1.2 complex. Similarly, MTX adopted the β-sheet structure as the interacting surface for binding both channels, with Lys23 occluding the pore. In contrast, the other critical residues Lys27, Lys30, and Tyr32 of MTX adopted distinct interactions when associating with the IKCa channel. In addition, the residues Gln229, Ala230, Ala233, and Thr234 on the IKCa channel turret formed polar and non-polar interactions with MTX, whereas the turret of Kv1.2 was almost not involved in recognizing MTX. In all, the pairs of interacting residues on MTX and the IKCa channel of the bound complex indicated that electrostatic and Van der Waal interactions contributed equally to the formation of a stable MTX-IKCa complex, in contrast to the MTX-Kv1.2 binding that is dominantly mediated by electrostatic forces.

Conclusions

Despite sharing similar pharmacological profiles toward both IKCa and Kv1.2 channels, MTX adopted totally diverging modes in the two association processes. All the molecular information unveiled here could not only offer a better understanding about the structural differences between the IKCa and Kv1.2 channels, but also provide novel structural clews that will help in the designing of more selective molecular probes to discriminate between these two channels.

Background

Scorpion venoms produce a large variety of peptide toxins that target ion channels [1-5]. Especially, the widespread use of scorpion-venom peptides acting on K+-channels as neuroscience tools and excellent ligand models has tremendously increased our knowledge in many fields, including exploration of the 3-dimensional structures and elucidation of the pharmacological characteristics of K+ channels [4,6-8]. In addition, peptide toxins are increasingly recognized as valuable sources of new drugs for channelopathies [9,10]. Although natural toxins often lack sufficient efficacy and specificity toward an individual channel type, most peptide toxins adopt a cysteine-stabilized α/β scaffold; thus, they could serve as good candidates for further structure-based drug design [4,10]. However, crystal structures for many medically important potassium channels have not been determined, which makes the rational designing of K+-channel modulators difficult. Therefore, applying computational methods to model reasonably stable structures of channel-peptide toxin complexes could be a good alternative, which would greatly help to highlight the diversity of channel-toxin interactions and provide structural information for toxin-based drug design.

The intermediate-conductance calcium-activated potassium channels (IKCa) act as positive modulators of cell proliferation by hyperpolarizing the cell membrane in T and B cells, fibroblasts, and vascular smooth muscle cells [11-13]. Furthermore, blocking of IKCa channels has been shown to be a potential therapeutic strategy against autoimmune disorders involving these tissues [13-15]. However, almost all the peptidic and small molecular IKCa blockers could not discriminate well between the IKCa channel and other related Kv-family channels and, thus, lack the specificity needed for further drug development [12,13,15,16].

Maurotoxin (MTX), a peptide derived from the venom of the scorpion Scorpio maurus palmatus, is the most potent peptidic blocker of the IKCa channel [17,18]. In addition, MTX could distinguish the IKCa channel from the other calcium-activated channels and the Kv1-family channels, except for the voltage-gated Kv1.2 channel [17-21]. Interestingly, although the IKCa channel is entirely different from the Kv1.2 channel in tissue contribution and physiological function [11,12,15,16], MTX shows very similar pharmacological profiles in recognizing these two channels with approximately the same potency and using the same functional residues [17,18,21]. In this study, we aimed to interpret the differential binding mechanisms of MTX with reference to the IKCa and Kv1.2 channels, which would provide a deep insight into the topological differences of these two channels and offer important clues for designing inhibitors that are more selective toward the therapeutic IKCa channel.

Combined computational methods were used to investigate the details of the interactions between MTX and the IKCa channel; the structural details were further compared with the previous model of the MTX-Kv1.2 complex [22]. A stable structure of the MTX-IKCa complex was obtained by using ligand docking, clustering analysis, and molecular dynamics simulation (MDS) methods. The validity of the final MTX-IKCa complex was supported by good accordance between the computational alanine-scanning results and the experimental data. On comparison, although with similar pharmacological profiles, MTX adopted very different modes for associating with the IKCa and the Kv1.2 channels. In both MTX-IKCa and MTX-Kv1.2 complexes, MTX adopted the β-sheet domain as the interaction surface with the Lys23 occluding the pore. However, the other key residue, Tyr32, was positioned quite differently in these two complexes. Meanwhile, the turret region of the IKCa channel played an important role in binding with MTX, which is different from the noninvolvement of the Kv1.2 turret during interaction with MTX. In addition, due to the different physicochemical profiles of the two channels, electrostatic and van der Waals (vdW) interactions made different contributions to the free energies of binding in MTX-Kv1.2 and MTX-IKCa complexes. All these structural and energetic discrepancies constitute the key determinants responsible for the binding specificity of MTX to the IKCa and Kv1.2 channels, which could help design MTX derivatives that would discriminate between these two channels.

Results and Discussion

Different MTX binding modes towards IKCa and Kv1.2 channels

Given the similar blocking activities that MTX showed toward IKCa and Kv1.2 channels [17,18], we first tested whether MTX recognized the IKCa channel in the same binding mode as for the Kv1.2 channel. An MTX-IKCa complex was modeled on the basis of the MTX-Kv1.2 complex structure [22] using the distance-restraint homologous modeling method and then subjected to unrestrained MDS to test its stability. Previously, MTX adopted the β-sheet domain as the interacting surface, with Lys23 occluding the pore of the Kv1.2 channel (Figure 1A, left panel). And the most important residue, Tyr32, of MTX kept sticking into the pocket formed by residues on the turret and pore region of the Kv1.2 channel during the 8-ns production run [22] (Figure 1A). Thus, the interaction between MTX and the Kv1.2 channel was beneficial for the stability of the turret region. On the contrary, in the MTX-IKCa complex, in spite of starting from a similar conformation of sticking into the pocket formed by the turret and the pore region (Figure 1B, left panel), the Tyr32 was unwelcome in this position and gradually bent away from the "pocket" after 4-ns MDS (Figure 1B, middle panel). Therefore, the experimentally important residues, Lys23 and Tyr32, failed to form strong interactions with the channel, and MTX deviated from the central position (Figure 1B, right panel). As a result, the turret on the other side of the IKCa channel failed to maintain stability and bent outward (see Figure 1C the RMSD of IKCa channel). Obviously, this complex structure could not explain the experiment results at all. Therefore, although MTX showed similar pharmacological profiles toward both Kv1.2 and IKCa channels, it did not use the same binding mode in associating with these two channels.

thumbnailFigure 1. The stability of an MTX-IKCa candidate that resembles the MTX-Kv1.2 complex. (A) The structural comparison of the MTX-Kv1.2 complex before and after 8-ns unrestrained MDS. (B) The structural comparison of the MTX-Kv1.2-like candidate for the MTX-IKCa complex before and after 4-ns unrestrained MDS. Complexes before and after MDS are distinguished by colors. (C) The RMSd fluctuation of MTX (left), IKCa (middle) and complex (right) of the MTX-Kv1.2-like candidate for the MTX-IKCa complex in 4-ns MDS.

MTX-IKCa complex from docking and MDS

We next applied a routine molecular docking and clustering analysis to screen plausible MTX-IKCa complexes [22-25]. All the 35 nuclear magnetic resonance conformations of MTX with different side-chain positions were used in the ZDOCK program, and 35, 000 complexes were generated in total.

The 35, 000 MTX-IKCa complexes fell into four main binding modes, according to the orientation of the MTX β-sheet domain (Figure 2, top panel). As indicated in our previous study, the starting position of the residue Tyr32 of MTX is essential for obtaining a final stable complex [22]. Thus, we mapped all possible orientations of the aromatic ring of Tyr32 in the four groups of complexes (Figure 2, middle panel), including the following: (1) lying slantways over the linker connecting S6 and the selectivity filter of the IKCa channel; (2) hanging between two subunits of the channel; (3) plugging the pore of the channel along with Lys23; (4) hanging upright over the linker connecting S6 and the selectivity filter of the IKCa channel. Next, each of the candidates in the four binding modes was subjected to energy minimization, followed by 500-ps unrestrained MDS, to introduce more flexibility and investigate the potential structural fluctuations. For better discrimination among the four binding modes, the alanine-scanning method in the MM-GBSA of Amber-8 was applied for comparison of the computational and experimental data of functional residues (Figure 2, bottom panel). Finally, an MTX-IKCa model that best fit the experimental data was selected from the binding mode II. To make it more equilibrated and stable, the selected MTX-IKCa model was subjected to a further 10-ns unrestrained MDS (Figure 3A). According to the major transitions during the MDs, we also evaluated the free energy convergence by sampling along the whole simulation (Figure 3B). The results showed that both the peptide and the channel were induced to fit in their nanosecond-scale diffusional encounter, and then guided to form a specific stable complex.

thumbnailFigure 2. Structural view of the four selected MTX-IKCa channel-binding modes along with the comparison of the calculated and experimental effects. Each row of figures is as follows: (top) the different spatial orientations of MTX in its complex with the IKCa channel from four main binding modes; (middle) the representative positions of the Lys23 and Tyr32 side chains in four main binding modes distinguished by different colors; (bottom) the comparison of the calculated and experimental effects for the five alanine mutations of MTX on the binding affinity toward the IKCa channel after 500-ps unrestrained MDS. The calculated results are normalized values of ΔΔGbinding, whereas the experimental results are obtained as kbT ln [IC50(mutant)/IC50(wt)].

thumbnailFigure 3. The stability and validity of the final MTX-IKCa complex. (A) Root-mean-square deviations (Å) of the α-carbons in the final MTX-IKCa complex from the starting complex during the 10-ns unrestrained MD simulations. (B) The free energy convergence during the 10-ns unrestrained MD simulations, trajectories were sampled every ns along the simulation time. (C) Calculated and experimental effects for 5 single-alanine mutations on the blocking activity of MTX toward the IKCa channel after 10-ns MDS. The calculated results are normalized values of ΔΔGbinding, whereas the experimental results are obtained as kbT ln [IC50(mutant)/IC50(wt)]. (D) An overview of the final model of the MTX-IKCa complex, with the labels of all the experimentally determined residues.

To check the confidence of our MTX-IKCa model, the ΔΔGbinding of eight single mutations of MTX were calculated and compared with the experimental data [17,18]. An overall high degree of correlation was found between the calculations and the experiments involving mutational effects (Figure 3C). Replacing the Lys23 residue of MTX with alanine caused the most noticeable decrease of 4.4 kcal/mol in the calculated binding energy, which is well in accordance with the experimental data of 4.32 kcal/mol. Substitution of another important residue, Tyr32, with alanine significantly reduced the MTX affinity by over 1000 fold [18], and the calculated ΔΔGbinding value of 4.31 kcal/mol corresponded well with the experimental data of 4.21 kcal/mol. However, MTX affinity for IKCa was decreased by less than 10 fold by the S2A, T4A, S6A, K7A, Y10A, and G33A mutants [18]. This is strongly supportive for the little change in binding energy when these residues were mutated to alanine; this is because Ser2, Thr4, Ser6, and Lys7 were located at the N-terminal of MTX, whereas Tyr10 was in the middle of the α-helix of MTX, all outside the interface of MTX.

Consistent with the findings in previous docking experiments of MTX onto IKCa channel [26], in the final MTX-IKCa complex, the peptide used its β-sheet as the interacting surface, with Lys23 as a structurally conserved pore-blocking residue (Figure 3D). This phenomenon was also observed when MTX associated with the Kv1.2 channel in both our study [22] and the previous docking results by Visan [18]. All these studies underline the key role of MTX β-sheet region in IKCa and Kv1.2 channel recognition. Interestingly, such importance is strongly supported by the previous experimental data that when substituting the β-sheet region of MTX with that of another toxin, HsTX1, its activity toward the two channels almost disappeared [27].

Despite that MTX used the β-sheet to interact with both IKCa and Kv1.2 channels, when analyzing the conformation of other bioactive residues, the molecular information for the recognition of IKCa by MTX showed several distinct features, compared to those for the Kv1.2 channel.

Differential molecular information contained by Tyr32 and Lys7 in MTX

The mutant-cycle experiment showed that two major functional residues, namely, Lys23 and Tyr32, in MTX block the IKCa channel [17,18]. The key role of these residues can be shown by their structural conformation in the final MTX-IKCa complex. Although the pore-blocking characteristic of Lys23 was common when MTX recognized both the IKCa and the Kv1.2 channels, the conformation of Tyr32 of MTX differed greatly in each of these contexts (Figure 4A, B). As indicated in the second section, the residue Tyr32 is essential for MTX-associating potassium channels, and its position was proposed to be the main factor determining the stability of the complex structure after a long-term unrestrained production run. In the final MTX-IKCa complex, the favorable position for the residue Tyr32 was lying on the linker connecting the selectivity filter and the S6 helix of the channel, forming strong polar and nonpolar interactions with Gly254 (D chain), Asp255 (C and D chains), Val256 (D chain), and Val257 (D chain) on the pore region of the channel (Figure 4A). In comparison, the Tyr32 in the MTX-Kv1.2 complex, however, differed significantly from that in the MTX-IKCa complex by sticking into the pocket formed by Arg354 and F358 on the turret and Asp379, Met380, Val381, and Thr383 on the pore region of the Kv1.2 channel [22] (Figure 4B). Thus, the differential chemical environment for Tyr32 contributed greatly to the different binding modes of MTX when associating with the IKCa and Kv1.2 channels.

thumbnailFigure 4. The distinct recognition mechanisms of MTX toward the IKCa and Kv.1.2 channels. (A), (C) Interaction details of Tyr32 and Lys7 of MTX in the MTX-IKCa complex within a contact distance of 5Å, respectively. (B), (D) Interaction details of Tyr32 and Lys7 of MTX in the MTX-Kv1.2 complex within a contact distance of 5Å, respectively.

In addition, rather than forming strong electrostatic interactions with the aspartic acid residues in the pore region of the channel, the side chain of the Lys7 of MTX pointed to an opposite orientation from the IKCa channel, contacting no residue of the channel within a 5-Å distance (Figure 4C). This is different from the result of the mutant-cycle analysis that Lys7 is situated near the Asp239 of the IKCa channel [17,18]. However, as the K7A mutation only affected the blocking activity of MTX by less than 10 fold [17,18], it is possible that the Lys7 just faced its alternative partner Asp239 in the interface reorganization process, but does not contact Asp239 directly in the final conformation. Such a position of the Lys7 while associating with the IKCa channel differed significantly from that in the MTX-Kv1.2 complex, in which the Lys7 formed strong polar interactions with the Asp373 at the pore region of the channel [22] (Figure 4D). This distinctness was in consistence with the experimental data that the blocking activity of the MTX-K7A mutant decreased by about 100 fold in the case of the Kv1.2 channel, but decreased by only less than 10 fold in the case of the IKCa channel [18].

Differential molecular information of other residues in MTX

Considering that the β-sheets of MTX constitute its main channel-interacting surfaces, two other residues Lys27 and Lys30 were also found to play outstanding roles in the recognition process (Figure 5A, B). In the model of the MTX-IKCa complex, the Lys30 mainly interacted with Gly254 (D chain) and Val257 (A chain) at the pore region of the IKCa channel (Figure 5A); in addition, it also formed some electrostatic interactions with the backbone of Asp255 (D chain). Another basic residue Lys27, however, did not electrostatically interact with any acidic residues of the channel, but rather stuck into a pocket formed mostly by nonpolar residues, including Thr234 and Trp221 on the turret; and Val256, Val257, and Pro258 on the linker connecting the S6 helix and the selectivity filter (Figure 5A). Such finding is completely in line with the previous docking result of MTX onto IKCa channel [26], which highlighted the important contribution of Lys27 and Lys30 in the interaction, due to their proximity to Asp255, Val257 and Asp239 of the channel.

thumbnailFigure 5. The comparison for the pairs of MTX-IKCa and MTX-Kv1.2 interacting residues in the bound complexes. Lys27 and Lys30 of MTX are surrounded by residues from the IKCa channel (A) and Kv1.2 channel (B) within a distance of 5Å, respectively.

However, these interaction modes mediated through Lys27 and Lys30 of the MTX in recognizing the IKCa channel differed obviously from those involving MTX and the Kv1.2 channel, in which both Lys27 and Lys30 formed strong polar interactions with the channel [22]. Within a contact distance of 5Å, the Lys27 and Lys30 of MTX, respectively, contacted closely with the conserved acidic residue, Asp379, in the pore region of the Kv1.2 channel and formed strong electrostatic interactions (Figure 5B).

The different functional roles for the Lys27 and Lys30 of MTX when recognizing the IKCa and the Kv1.2 channels could be further illustrated by calculating their mutation effects. As indicated in Table 1, either changing Lys27 or Lys30 of MTX into alanine resulted in a big change in the MTX-IKCa and MTX-Kv1.2 interactive energies. However, in the model of the MTX-IKCa complex, the changes in the vdW energies and electrostatic energies almost synergistically contributed to the large values of the altered interactive energies when replacing Lys27 with alanine (Table 1), in accordance with the observation that the Lys27 of MTX mainly interacted with the nonpolar residues in IKCa channel. In contrast, in the MTX-Kv1.2 complex, the increase in the interactive energies caused by mutating Lys27 and Lys30 were both dominantly constituted by the significant change in the electrostatic energies (Table 1), resulting from the interacting residue pairs Lys27-Asp373 and Lys30-Asp379 (Figure 5B).

Table 1. Effects of K27A and K30A mutants on interactive energies (Kcal/mol)

Distinct channel vestibules constitute different recognition modes toward maurotoxin

The α-KTx family of K+-channel blockers has been proved to function as informative molecular probes for the structure-function analysis of K+ channels. Although the IKCa and Kv1.2 channels have distinct tissue distributions and biophysical features [11-16], both can be blocked by MTX with a similar pharmacology profile. Thus, identifying the differential determinants that are responsible for the MTX binding of the IKCa and Kv1.2 channels could help discover the different topologies of a mechanistically interesting part of these two channels: the outer vestibule of the ion-conduction pore.

Sequence comparison showed that the IKCa channel has a longer turret region than the Kv1.2 channel, and the sequence identity is rather low (Figure 6A). These features were further shown by the structural analysis of unbound IKCa and Kv1.2 channel vestibules (Figure 6B), whereas the differential conformation of the channel turret suggested a different functional role for the channel vestibule during the toxin-recognition process. This hypothesis was supported by comparing our model of the IKCa-MTX complex with the previous Kv.1.2-MTX model [22] (Figure 6C, D). We mapped all the toxin-interacting residues in the IKCa channel turret (see Figure 6C and additional file 1: Figure S1) and found that the channel turret played an important role in stabilizing the final IKCa-MTX complex. Four residues, Gln229, Ala230, Ala233, and Thr234 in the IKCa turret formed polar and non-polar interactions with the corresponding residues of MTX. Interactions have been found between Gln229 of IKCa and Gly18, Pro20 of MTX, Ala230 of IKCa and Pro20 of MTX, Ala233 of IKCa and Val1, Thr17 of MTX, Thr234 of IKCa and Lys27, Tyr32, Gly33 of MTX, respectively. Thus, the important role of the IKCa channel's turret in recognizing MTX is in sharp contrast with the almost nil involvement of the Kv1.2 channel turret, because only the Tyr32 of MTX was able to contact the Kv1.2 turret within the distance of 5Å (Figure 6D) and several mutations on the Kv1.2 channel turret hardly affected the binding affinity of MTX [18,22]. In addition, it is noticeable that the total net charge of the turret of the IKCa channel is almost neutral, due to the average distribution of positive and negative residues in the turret (Figure 6A, E). This is very different compared with the Kv1.x family, all of which contained four extremely negative-charged turrets (Figure 6A, F). Actually, the IKCa channel turret, with an important functional role and neutral charge, could well explain the equal contribution of electrostatic and vdW interactions in mediating the recognition process between the IKCa channel and MTX (Table 2), and the highly negative-charged Kv1.2 turrets are responsible for the dominance of electrostatic recognition in the binding of MTX to the Kv1.2 channel (Table 2). All these data suggest that the diverse vestibules of the IKCa and Kv1.2 channels, including the sequence length, sequence identity, vestibule conformation, and property of molecular surface jointly determine their different interaction modes.

thumbnailFigure 6. The comparison for the sequence, conformation, and functional role of the outer vestibules of the IKCa and Kv1.2 channels. (A) Sequence alignment for the outer vestibule of the IKCa and Kv1.2 channels. The conserved residues are grey-shaded, acidic residues red-shaded, basic residues blue-shaded. (B) Left and right are, respectively, the side and top views of the ribbon structures of the outer vestibules of the IKCa and the Kv1.2 channels, respectively. The IKCa channel is pink-colored, and the Kv1.2 channel is green-colored. For clarity, only two chains of the channels are shown in the side view. (C), (D) Functional roles for the turrets of the IKCa and Kv1.2 channels, respectively, in associating with MTX, respectively. All the toxin-interacting residues within a contact distance of 5Å are labeled on the turret of the IKCa and Kv1.2 channels. (E), (F) Stereoview displaying the molecular surfaces of the IKCa and Kv1.2 channels, respectively. The molecular surfaces are colored according to their electrostatic potentials, with basic residues in blue, acidic residues in red, and neutral residues in white.

Additional file 1. The interaction details between the D chain of IKCa channel and MTX. (A) Gln229 on the channel turret interacts with Gly18 and Pro20 of MTX. (B) Ala230 on the channel turret interacts with Pro20 of MTX. (C) Ala233 on the channel turret interacts with Thr17 of MTX. (D) Thr234 on the channel turret interacts with Tyr32 and Gly33 of MTX.

Format: JPEG Size: 2.3MB Download fileOpen Data

Table 2. Different polar and nonpolar contributions to the binding free energies in MTX-IKCa and MTX-Kv1.2 complexes

Conclusions

Through combined computational methods, including ZDOCK, clustering analysis, and MDS, a reasonably stable MTX-IKCa complex structure was obtained. Further study of this structure showed that in spite of sharing similar pharmacological profiles toward both IKCa and Kv1.2 channels, MTX associated with the IKCa channel in a quite different mode compared to that of MTX interacting with the Kv1.2 channel. In the bound complex, MTX assumed the β-sheet domain as the interaction surface with the Lys23 occluding the pore of the IKCa channel in a manner similar to its interaction with the Kv1.2 channel. However, the conformation of another key residue Tyr32, which was the key to the stability of the complex structure, differed greatly when MTX recognized the IKCa channel, compared to the process with the Kv1.2 channel. It continued lying on the linker connecting the selectivity filter and the S6 helix of the IKCa channel, forming strong polar and nonpolar interactions with residues on the pore region of the channel. In addition, the Lys7 of MTX is possibly involved in the toxin-channel interface reorganization process; however, it does not contact any residues of the IKCa channel directly in the final conformation. This is in contrast with the fact that the Lys7 of MTX formed strong polar interactions with the Asp373 at the pore region of the Kv1.2 channel. In addition, electrostatic and vdW interactions contributed equally to the binding of MTX with IKCa, whereas the MTX-Kv1.2 association featured dominant electrostatic contribution. Such conformational and energetic differences in recognition could be well explained by the different functional roles of the channel vestibules. The longer, neutral-charged IKCa channel turret played an important role in stabilizing the final IKCa-MTX complex, with four residues--Gln229, Ala230, Ala233, and Thr234 forming polar and non-polar interactions with MTX. On the contrary, the shorter Kv1.2-channel turret is highly negatively charged and is barely involved in recognizing MTX. In all, the differences in the binding mechanisms of MTX toward the IKCa and Kv1.2 channels unveiled in this study could offer a better understanding of the physicochemical properties and conformational distinctness of the two channels and thus give a hint for designing MTX-derived inhibitors to discriminate between these two channels.

Methods

Atomic Coordinates and Molecular Docking

The atomic coordinates of MTX (PDB code: 1TXM) were downloaded from the PDB [28]. The previous segment-assembly homology model was applied to obtain the structure of the pore region of the IKCa channel [24]. This model was then subjected to 5-nanosecond (ns) MDS for equilibration.

To improve the docking performance, all 35 conformations of MTX were used to dock with the equilibrated IKCa structure through the ZDOCK program [29], a fast Fourier transform (FFT)-based, initial-stage rigid-body molecular-docking algorithm. Each docking produced 1000 candidate complexes, thus 35000 candidate MTX-IKCa complexes were obtained and used for the clustering analysis. According to the orientation of the MTX β-sheet domain, the 35000 complexes were then divided into four main binding modes. Clustering analysis and experimental data-based screening [17,18] were then carried out on all the complexes to select the possible hits from all modes. Candidates from each binding mode were then subjected to a 500-step energy minimization using the Sander module of the Amber-8 suit of programs [30]. By calculating the ligand-receptor binding energies with the ANAL program of Amber-8, appropriate candidate complexes were identified for further MDS study.

MDS study

All the simulations in this work were carried out using the Amber-8 program [30] on a 64-CPU Dawning TC4000L cluster (Beijing, China). The generalized Born model [31], which has been successfully used to study other toxin-channel interactions [22,24,25,32,33], was applied in this study.

All the candidate complexes selected by the screening process went through 400-picosecond (ps) equilibration and 500-ps unrestrained simulations to introduce more flexibility. The equilibration steps were taken by gradually reducing the force constant--from 5.0 (kcal/mol)/Å2 for restraining all the heavy atoms, to 0.02 (kcal/mol)/Å2 for heavy atoms of the backbone only. The temperature was set at 300 K, with a cutoff distance of 12 Å. For the most reasonably stable complex selected after a 500-ps unrestrained simulation, an additional 10-ns unrestrained simulation was conducted to introduce enough flexibility and to probe into the interaction details. Throughout all the energy minimization and simulation processes, the ff99 force field (Parm 99) [34] was applied.

During the simulation, the membrane around the channel has not been taken into account. It is because that the scorpion peptide binds to the extracellular part of the channel according to mutagenesis studies and solid-state NMR results [6,35-37], where the interaction is hardly affected by the membrane and the transmenbrane segment of channel. Other study groups have also used the same membrane-ignoring measures in molecular simulation studies of toxin-channel interactions [22-25,32,33,38-40]. However, the importance of the membrane in the functioning of channels has been increasingly recognized. A transmembrane protein system could be more reliable if the role of the membrane were taken into account.

Calculation of Free energy of Binding by the Molecular mechanics--Generalized Born Surface Area method

In the molecular mechanics--generalized born surface area (MM-GBSA) method of AMBER-8 [30], the free energy of binding of the reaction A + B → AB is calculated using the following thermodynamic cycle:

<a onClick="popup('http://www.biomedcentral.com/1472-6807/11/3/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1472-6807/11/3/mathml/M2">View MathML</a>

(1)

<a onClick="popup('http://www.biomedcentral.com/1472-6807/11/3/mathml/M3','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1472-6807/11/3/mathml/M3">View MathML</a>

(2)

<a onClick="popup('http://www.biomedcentral.com/1472-6807/11/3/mathml/M4','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1472-6807/11/3/mathml/M4">View MathML</a>

(3)

<a onClick="popup('http://www.biomedcentral.com/1472-6807/11/3/mathml/M5','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1472-6807/11/3/mathml/M5">View MathML</a>

(4)

where T is the temperature, S is the solute entropy, ΔGgas is the interaction energy between A and B in the gaseous phase, and <a onClick="popup('http://www.biomedcentral.com/1472-6807/11/3/mathml/M6','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1472-6807/11/3/mathml/M6">View MathML</a>, <a onClick="popup('http://www.biomedcentral.com/1472-6807/11/3/mathml/M7','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1472-6807/11/3/mathml/M7">View MathML</a>, and <a onClick="popup('http://www.biomedcentral.com/1472-6807/11/3/mathml/M8','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1472-6807/11/3/mathml/M8">View MathML</a> are the solvation free energies of A, B, and AB, which are estimated using the GBSA method [30]. That is, <a onClick="popup('http://www.biomedcentral.com/1472-6807/11/3/mathml/M9','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1472-6807/11/3/mathml/M9">View MathML</a>, and so forth. ΔGGB and ΔGSA are the electrostatic and nonpolar terms, respectively. ΔEbond, ΔEangle, and ΔEtorsion are contributions to the intramolecular energy ΔEintra of the complex. EvdW is vdW interaction energy. Because of the constant contribution of -TΔS for each docked complex, we quote ΔG*binding for ΔGbinding + TΔS in the discussion. To verify the quality and validity of the resulting MTX-IKCa complexes, the relative free energy of binding, ΔG*binding, was calculated using MM-GBSA for postprocessing-collected snapshots from the MD trajectories. In this work, 30 snapshots from the last 30-ps MDS were used for analysis of the free energy of binding.

Authors' contributions

HY conceived of the study, analyzed the data and drafted the manuscript. SQ performed the docking analysis, carried out the molecular simulation studies and analyzed the data. YLW and WXL participated in its design and coordination, and helped to draft the manuscript. BSW supervised the study design, coordination and edited the manuscript. All authors read and approved the final manuscript.

Acknowledgements

This work was supported by grants from the National Natural Sciences Foundation of China (number 30900265, 30770519 and 30973636), the China Postdoctoral Science Foundation (number 20090451075) and the National Basic Research Program of China (2010CB529800).

References

  1. Darbon H, Blanc E, Sabatier JM: Three-dimensional structure of scorpion toxins: Towards a new model of interaction with potassium channels.

    Perspect Drug Discovery Des 1999, 16:41-60. Publisher Full Text OpenURL

  2. Possani LD, Selisko B, Gurrola GB: Structure and function of scorpion toxins affecting K+-channels.

    Perspect Drug Discovery Des 1999, 16:15-40. Publisher Full Text OpenURL

  3. Possani LD, Merino E, Corona M, Bolivar F, Becerril B: Peptides and genes coding for scorpion toxins that affect ion-channels.

    Biochimie 2000, 82:861-868. PubMed Abstract | Publisher Full Text OpenURL

  4. Rodriguez de la Vega RC, Merino E, Becerril B, Possani LD: Novel interactions between K+ channels and scorpion toxins.

    Trends Pharmacol Sci 2003, 24:222-227. PubMed Abstract | Publisher Full Text OpenURL

  5. Giangiacomo KM, Ceralde Y, Mullmann TJ: Molecular basis of alpha-KTx specificity.

    Toxicon 2004, 43:877-886. PubMed Abstract | Publisher Full Text OpenURL

  6. Hidalgo P, MacKinnon R: Revealing the architecture of a K+ channel pore through mutant cycles with a peptide inhibitor.

    Science 1995, 268:307-310. PubMed Abstract | Publisher Full Text OpenURL

  7. MacKinnon R, Cohen SL, Kuo A, Lee A, Chait BT: Structural conservation in prokaryotic and eukaryotic potassium channels.

    Science 1998, 280:106-109. PubMed Abstract | Publisher Full Text OpenURL

  8. Ruta V, Jiang Y, Lee A, Chen J, MacKinnon R: Functional analysis of an archaebacterial voltage-dependent K+ channel.

    Nature 2003, 422:180-185. PubMed Abstract | Publisher Full Text OpenURL

  9. Garber K: Peptide leads new class of chronic pain drugs.

    Nat Biotechnol 2005, 23:399. PubMed Abstract | Publisher Full Text OpenURL

  10. Lewis RJ, Garcia ML: Therapeutic potential of venom peptides.

    Nat Rev Drug Discov 2003, 2:790-802. PubMed Abstract | Publisher Full Text OpenURL

  11. Wulff H, Knaus HG, Pennington M, Chandy KG: K+ channel expression during B cell differentiation: implications for immunomodulation and autoimmunity.

    J Immunol 2004, 173:776-786. PubMed Abstract | Publisher Full Text OpenURL

  12. Wei AD, Gutman GA, Aldrich R, Chandy KG, Grissmer S, Wulff H: International Union of Pharmacology. LII. Nomenclature and molecular relationships of calcium-activated potassium channels.

    Pharmacol Rev 2005, 57:463-472. PubMed Abstract | Publisher Full Text OpenURL

  13. Wulff H, Kolski-Andreaco A, Sankaranarayanan A, Sabatier JM, Shakkottai V: Modulators of small- and intermediate-conductance calcium-activated potassium channels and their therapeutic indications.

    Curr Med Chem 2007, 14:1437-1457. PubMed Abstract | Publisher Full Text OpenURL

  14. Wulff H, Miller MJ, Hansel W, Grissmer S, Cahalan MD, Chandy KG: Design of a potent and selective inhibitor of the intermediate-conductance Ca2+-activated K+ channel, IKCa1: a potential immunosuppressant.

    Proc Natl Acad Sci USA 2000, 97:8151-8156. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  15. Wulff H, Zhorov BS: K+ channel modulators for the treatment of neurological disorders and autoimmune diseases.

    Chem Rev 2008, 108:1744-1773. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  16. Gutman GA, Chandy KG, Grissmer S, Lazdunski M, McKinnon D, Pardo LA, Robertson GA, Rudy B, Sanguinetti MC, Stuhmer W, Wang X: International Union of Pharmacology. LIII. Nomenclature and molecular relationships of voltage-gated potassium channels.

    Pharmacol Rev 2005, 57:473-508. PubMed Abstract | Publisher Full Text OpenURL

  17. Castle NA, London DO, Creech C, Fajloun Z, Stocker JW, Sabatier JM: Maurotoxin: a potent inhibitor of intermediate conductance Ca2+-activated potassium channels.

    Mol Pharmacol 2003, 63:409-418. PubMed Abstract | Publisher Full Text OpenURL

  18. Visan V, Fajloun Z, Sabatier JM, Grissmer S: Mapping of maurotoxin binding sites on hKv1.2, hKv1.3, and hIKCa1 channels.

    Mol Pharmacol 2004, 66:1103-1112. PubMed Abstract | Publisher Full Text OpenURL

  19. Kharrat R, Mabrouk K, Crest M, Darbon H, Oughideni R, Martin-Eauclaire MF, Jacquet G, el Ayeb M, Van Rietschoten J, Rochat H, Sabatier JM: Chemical synthesis and characterization of maurotoxin, a short scorpion toxin with four disulfide bridges that acts on K+ channels.

    Eur J Biochem 1996, 242:491-498. PubMed Abstract | Publisher Full Text OpenURL

  20. Kharrat R, Mansuelle P, Sampieri F, Crest M, Oughideni R, Van Rietschoten J, Martin-Eauclaire MF, Rochat H, El Ayeb M: Maurotoxin, a four disulfide bridge toxin from Scorpio maurus venom: purification, structure and action on potassium channels.

    FEBS Lett 1997, 406:284-290. PubMed Abstract | Publisher Full Text OpenURL

  21. Avdonin V, Nolan B, Sabatier JM, De Waard M, Hoshi T: Mechanisms of maurotoxin action on Shaker potassium channels.

    Biophys J 2000, 79:776-787. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  22. Yi H, Qiu S, Cao Z, Wu Y, Li W: Molecular basis of inhibitory peptide maurotoxin recognizing Kv1.2 channel explored by ZDOCK and molecular dynamic simulations.

    Proteins 2008, 70:844-854. PubMed Abstract | Publisher Full Text OpenURL

  23. Wu Y, Cao Z, Yi H, Jiang D, Mao X, Liu H, Li W: Simulation of the interaction between ScyTx and small conductance calcium-activated potassium channel by docking and MM-PBSA.

    Biophys J 2004, 87:105-112. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  24. Yi H, Cao Z, Yin S, Dai C, Wu Y, Li W: Interaction simulation of hERG K+ channel with its specific BeKm-1 peptide: insights into the selectivity of molecular recognition.

    J Proteome Res 2007, 6:611-620. PubMed Abstract | Publisher Full Text OpenURL

  25. Qiu S, Yi H, Liu H, Cao Z, Wu Y, Li W: Molecular Information of charybdotoxin blockade in the large conductance calcium-activated potassium channel.

    J Chem Inf Model 2009, 49:1831-1838. PubMed Abstract | Publisher Full Text OpenURL

  26. Visan V, Sabatier JM, Grissmer S: Block of maurotoxin and charybdotoxin on human intermediate-conductance calcium-activated potassium channels (hIKCa1).

    Toxicon 2004, 43:973-980. PubMed Abstract | Publisher Full Text OpenURL

  27. Regaya I, Beeton C, Ferrat G, Andreotti N, Darbon H, De Waard M, Sabatier JM: Evidence for domain-specific recognition of SK and Kv channels by MTX and HsTx1 scorpion toxins.

    J Biol Chem 2004, 279:55690-55696. PubMed Abstract | Publisher Full Text OpenURL

  28. Berman HM, Westbrook J, Feng Z, Gilliland G, Bhat TN, Weissig H, Shindyalov IN, Bourne PE: The Protein Data Bank.

    Nucleic Acids Res 2000, 28:235-242. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  29. Chen R, Li L, Weng Z: ZDOCK: an initial-stage protein-docking algorithm.

    Proteins 2003, 52:80-87. PubMed Abstract | Publisher Full Text OpenURL

  30. Case DA, Darden TA, Cheatham TEI, Simmerling CL, Wang J, Duke RE, Luo R, Merz KM, Wang B, Pearlman DA, et al.: Amber 8. San Francisco: University of California; 2004. OpenURL

  31. Darden TA, Case DA: Molecular Dynamics Simulations of Nucleic Acids with a Generalized Born solvation model.

    J Am Chem Soc 2000, 122:2489-2498. Publisher Full Text OpenURL

  32. Han S, Yi H, Yin SJ, Chen ZY, Liu H, Cao ZJ, Wu YL, Li WX: Structural basis of a potent peptide inhibitor designed for Kv1.3 channel, a therapeutic target of autoimmune disease.

    J Biol Chem 2008, 283(27):19058-65. PubMed Abstract | Publisher Full Text OpenURL

  33. Yin SJ, Jiang L, Yi H, Han S, Yang DW, Liu ML, Liu H, Cao ZJ, Wu YL, Li WX: Different residues in channel turret determining the selectivity of ADWX-1 inhibitor peptide between Kv1.1 and Kv1.3 channels.

    J Proteome Res 2008, 7:4890-4897. PubMed Abstract | Publisher Full Text OpenURL

  34. Wang JM, Cieplak P, Kollman PA: How well does a restrained electrostatic potential (RESP) model perform in calculating conformational energies of organic and biological molecules?

    J Comput Chem 2000, 21:1049-1074. Publisher Full Text OpenURL

  35. Aiyar J, Withka JM, Rizzi JP, Singleton DH, Andrews GC, Lin W, Boyd J, Hanson DC, Simon M, Dethlefs B, et al.: Topology of the pore-region of a K+ channel revealed by the NMR-derived structures of scorpion toxins.

    Neuron 1995, 15:1169-1181. PubMed Abstract | Publisher Full Text OpenURL

  36. Gross A, MacKinnon R: Agitoxin footprinting the shaker potassium channel pore.

    Neuron 1996, 16:399-406. PubMed Abstract | Publisher Full Text OpenURL

  37. Lange A, Giller K, Hornig S, Martin-Eauclaire MF, Pongs O, Becker S, Baldus M: Toxin-induced conformational changes in a potassium channel revealed by solid-state NMR.

    Nature 2006, 440:959-962. PubMed Abstract | Publisher Full Text OpenURL

  38. Cui M, Shen JH, Briggs JM, Luo XM, Tan XJ, Jiang HL, Chen KX, Ji RY: Brownian dynamics simulations of interaction between scorpion toxin Lq2 and potassium ion channel.

    Biophys J 2001, 80:1659-1669. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  39. Cui M, Shen JH, Briggs JM, Fu W, Wu JL, Zhang YM, Luo XM, Chi ZW, Ji RY, Jiang HL, Chen KX: Brownian dynamics simulations of the recognition of the scorpion toxin P05 with the small-conductance calcium-activated potassium channels.

    J Mol Bio 2002, 318:417-428. Publisher Full Text OpenURL

  40. Eriksson MA, Roux B: Modeling the structure of agitoxin in complex with the Shaker K+ channel: a computational approach based on experimental distance restraints extracted from thermodynamic mutant cycles.

    Biophys J 2002, 83:2595-2609. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL