Functional significance of the beta hairpin in the insecticidal neurotoxin omega-atracotoxin-Hv1a

omega-Atracotoxin-Hv1a is an insect-specific neurotoxin whose phylogenetic specificity derives from its ability to antagonize insect, but not vertebrate, voltage-gated calcium channels. In order to help understand its mechanism of action and to enhance its utility as a lead compound for insecticide development, we used a combination of protein engineering and site-directed mutagenesis to probe the toxin for key functional regions. First, we constructed a Hairpinless mutant in which the C-terminal beta-hairpin, which is highly conserved in this family of neurotoxins, was excised without affecting the fold of the residual disulfide-rich core of the toxin. The Hairpinless mutant was devoid of insecticidal activity, indicating the functional importance of the hairpin. We subsequently developed a highly efficient system for production of recombinant toxin and then probed the hairpin for key functional residues using alanine-scanning mutagenesis followed by a second round of mutagenesis based on initial "hits" from the alanine scan. This revealed that two spatially proximal residues, Asn(27) and Arg(35), form a contiguous molecular surface that is essential for toxin activity. We propose that this surface of the beta-hairpin is a key site for interaction of the toxin with insect calcium channels.     

Solution structure of agelenin, an insecticidal peptide isolated from the spider Agelena opulenta, and its structural similarities to insect-specific calcium channel inhibitors

Agelenin, isolated from the Agelenidae spider Agelena opulenta, is a peptide composed of 35 amino acids. We determined the three-dimensional structure of agelenin using two-dimensional NMR spectroscopy. The structure is composed of a short antiparallel beta-sheet and four beta-turns, which are stabilized by three disulfide bonds. Agelenin has characteristic residues, Phe9, Ser28 and Arg33, which are arranged similarly to the pharmacophore of the insect channel inhibitor, omega-atracotoxin-Hv1a. These observations suggest that agelenin and omega-atracotoxin-Hv1a bind to insect calcium channels in a similar manner. We also suggest that another mode of action may operate in the channel inhibition by omega-agatoxin-IVA and omega-atracotoxin-Hv2a.     

Synthesis and characterization of Pi4, a scorpion toxin from Pandinus imperator that acts on K+ channels.

Pi4 is a 38-residue toxin cross-linked by four disulfide bridges that has been isolated from the venom of the Chactidae scorpion Pandinus imperator. Together with maurotoxin, Pi1, Pi7 and HsTx1, Pi4 belongs to the alpha KTX6 subfamily of short four-disulfide-bridged scorpion toxins acting on K+ channels. Due to its very low abundance in venom, Pi4 was chemically synthesized in order to better characterize its pharmacology and structural properties. An enzyme-based cleavage of synthetic Pi4 (sPi4) indicated half-cystine pairings between Cys6-Cys27, Cys12-32, Cys16-34 and Cys22-37, which denotes a conventional pattern of scorpion toxin reticulation (Pi1/HsTx1 type). In vivo, sPi4 was lethal after intracerebroventricular injection to mice (LD50 of 0.2 microg per mouse). In vitro, addition of sPi4 onto Xenopus laevis oocytes heterologously expressing various voltage-gated K+ channel subtypes showed potent inhibition of currents from rat Kv1.2 (IC50 of 8 pm) and Shaker B (IC50 of 3 nm) channels, whereas no effect was observed on rat Kv1.1 and Kv1.3 channels. The sPi4 was also found to compete with 125I-labeled apamin for binding to small-conductance Ca(2+)-activated K+ (SK) channels from rat brain synaptosomes (IC50 value of 0.5 microm). sPi4 is a high affinity blocker of the Kv1.2 channel. The toxin was docked (BIGGER program) on the Kv channel using the solution structure of sPi4 and a molecular model of the Kv1.2 channel pore region. The model suggests a key role for residues Arg10, Arg19, Lys26 (dyad), Ile28, Lys30, Lys33 and Tyr35 (dyad) in the interaction and the associated blockage of the Kv1.2 channel.     

Alanine scanning mutagenesis of a type 1 insulin-like growth factor receptor ligand binding site

The high resolution crystal structure of an N-terminal fragment of the IGF-I receptor, has been reported. While this fragment is itself devoid of ligand binding activity, mutational analysis has indicated that its N terminus (L1, amino acids 1-150) and the C terminus of its cysteine-rich domain (amino acids 190-300) contain ligand binding determinants. Mutational analysis also suggests that amino acids 692-702 from the C terminus of the alpha subunit are critical for ligand binding. A fusion protein, formed from these fragments, binds IGF-I with an affinity similar to that of the whole extracellular domain, suggesting that these are the minimal structural elements of the IGF-I binding site. To further characterize the binding site, we have performed structure directed and alanine-scanning mutagenesis of L1, the cysteine-rich domain and amino acids 692-702. Alanine mutants of residues in these regions were transiently expressed as secreted recombinant receptors and their affinity was determined. In L1 alanine mutants of Asp(8), Asn(11), Tyr(28), His(30), Leu(33), Leu(56), Phe(58), Arg(59), and Trp(79) produced a 2- to 10-fold decrease in affinity and alanine mutation of Phe(90) resulted in a 23-fold decrease in affinity. In the cysteine-rich domain, mutation of Arg(240), Phe(241), Glu(242), and Phe(251) produced a 2- to 10-fold decrease in affinity. In the region between amino acids 692 and 702, alanine mutation of Phe(701) produced a receptor devoid of binding activity and alanine mutations of Phe(693), Glu(693), Asn(694), Leu(696), His(697), Asn(698), and Ile(700) exhibited decreases in affinity ranging from 10- to 30-fold. With the exception of Trp(79), the disruptive mutants in L1 form a discrete epitope on the surface of the receptor. Those in the cysteine-rich domain essential for intact affinity also form a discrete epitope together with Trp(79).     

Brownian dynamics simulations of interaction between scorpion toxin Lq2 and potassium ion channel

The association of the scorpion toxin Lq2 and a potassium ion (K(+)) channel has been studied using the Brownian dynamics (BD) simulation method. All of the 22 available structures of Lq2 in the Brookhaven Protein Data Bank (PDB) determined by NMR were considered during the simulation, which indicated that the conformation of Lq2 affects the binding between the two proteins significantly. Among the 22 structures of Lq2, only 4 structures dock in the binding site of the K(+) channel with a high probability and favorable electrostatic interactions. From the 4 candidates of the Lq2-K(+) channel binding models, we identified a good three-dimensional model of Lq2-K(+) channel complex through triplet contact analysis, electrostatic interaction energy estimation by BD simulation and structural refinement by molecular mechanics. Lq2 locates around the extracellular mouth of the K(+) channel and contacts the K(+) channel using its beta-sheet rather than its alpha-helix. Lys27, a conserved amino acid in the scorpion toxins, plugs the pore of the K(+) channel and forms three hydrogen bonds with the conserved residues Tyr78(A-C) and two hydrophobic contacts with Gly79 of the K(+) channel. In addition, eight hydrogen-bonds are formed between residues Arg25, Cys28, Lys31, Arg34 and Tyr36 of Lq2 and residues Pro55, Tyr78, Gly79, Asp80, and Tyr82 of K(+) channel. Many of them are formed by side chains of residues of Lq2 and backbone atoms of the K(+) channel. Thirteen hydrophobic contacts exist between residues Met29, Asn30, Lys31 and Tyr36 of Lq2 and residues Pro55, Ala58, Gly79, Asp80 and Tyr82 of the K(+) channel. These favorable interactions stabilize the association between the two proteins. These observations are in good agreement with the experimental results and can explain the binding phenomena between scorpion toxins and K(+) channels at the level of molecular structure. The consistency between the BD simulation and the experimental data indicates that our three-dimensional model of Lq2-K(+) channel complex is reasonable and can be used in further biological studies such as rational design of blocking agents of K(+) channels and mutagenesis in both toxins and K(+) channels.     

Mapping function to structure in a channel-blocking peptide: electrostatic mutants of charybdotoxin.

Electrostatic interactions between charybdotoxin (CTX), a specific peptide pore blocker of K+ channels, and a Ca(2+)-activated K+ channel were investigated with a genetically manipulable recombinant CTX. Point mutations at certain charged residues showed only small effects on the binding affinity of the toxin molecule: Lys11, Glu12, Arg19, His21, Lys31, and Lys32. Replacement by Gln at Arg25, Lys27, or Lys34 strongly decreased the affinity of the toxin. These affinity changes were mainly due to large increases of toxin dissociation rates without much effect on association rates, as if close-range interactions between the toxin and its receptor site of the channel were disrupted. We also found that the neutralization of Lys27 to Gln removed the toxin's characteristic voltage dependence in dissociation rate. Mutation and functional mapping of charged residues revealed a molecular surface of CTX which makes direct contact with the extracellular mouth of the K+ channel.     

Scanning mutagenesis of a Janus-faced atracotoxin reveals a bipartite surface patch that is essential for neurotoxic function

The Janus-faced atracotoxins (J-ACTXs) are a family of insect-specific excitatory neurotoxins isolated from the venom of Australian funnel web spiders. In addition to a strikingly asymmetric distribution of charged residues, from which their name is derived, these toxins contain an extremely rare vicinal disulfide bond. To shed light on the mechanism of action of these toxins and to enhance their utility as lead compounds for insecticide development, we developed a recombinant expression system for the prototypic family member, J-ACTX-Hv1c, and mapped the key functional residues using site-directed mutagenesis. An alanine scan using a panel of 24 mutants provided the first complete map of the bioactive surface of a spider toxin and revealed that the entire J-ACTX-Hv1c pharmacophore is restricted to seven residues that form a bipartite surface patch on one face of the toxin. However, the primary pharmacophore, or hot spot, is formed by just five residues (Arg(8), Pro(9), Tyr(31), and the Cys(13)-Cys(14) vicinal disulfide). The Arg(8)-Tyr(31) diad in J-ACTX-Hv1c superimposes closely on the Lys-(Tyr/Phe) diad that is spatially conserved across a range of structurally dissimilar K(+) channel blockers, which leads us to speculate that the J-ACTXs might target an invertebrate K(+) channel.     

Structural and functional consequences of the presence of a fourth disulfide bridge in the scorpion short toxins: solution structure of the potassium channel inhibitor HsTX1

We have determined the three-dimensional structure of the potassium channel inhibitor HsTX1, using nuclear magnetic resonance and molecular modeling. This protein belongs to the scorpion short toxin family, which essentially contains potassium channel blockers of 29 to 39 amino acids and three disulfide bridges. It is highly active on voltage-gated Kv1.3 potassium channels. Furthermore, it has the particularity to possess a fourth disulfide bridge. We show that HsTX1 has a fold similar to that of the three-disulfide-bridged toxins and conserves the hydrophobic core found in the scorpion short toxins. Thus, the fourth bridge has no influence on the global conformation of HsTX1. Most residues spatially analogous to those interacting with voltage-gated potassium channels in the three-disulfide-bridged toxins are conserved in HsTX1. Thus, we propose that Tyr21, Lys23, Met25, and Asn26 are involved in the biological activity of HsTX1. As an additional positively charged residue is always spatially close to the aromatic residue in toxins blocking the voltage-gated potassium channels, and as previous mutagenesis experiments have shown the critical role played by the C-terminus in HsTX1, we suggest that Arg33 is also important for the activity of the four disulfide-bridged toxin. Docking calculations confirm that, if Lys23 and Met25 interact with the GYGDMH motif of Kv1.3, Arg33 can contact Asp386 and, thus, play the role of the additional positively charged residue of the toxin functional site. This original configuration of the binding site of HsTX1 for Kv1.3, if confirmed experimentally, offers new structural possibilities for the construction of a molecule blocking the voltage-gated potassium channels.     

Action of derivatives of mu-conotoxin GIIIA on sodium channels. Single amino acid substitutions in the toxin separately affect association and dissociation rates.

We have studied binding and block of sodium channels by 12 derivatives of the 22-residue peptide mu-conotoxin GIIIA (mu-CTX) in which single amino acids were substituted as follows: Arg or Lys by Gln, Gln-18 by Lys, Asp by Asn, and HO-Pro by Pro. Derivatives were synthesized as described by Becker et al. [(1989) Eur. J. Biochem. 185, 79]. Binding was measured by displacement of labeled saxitoxin from eel electroplax membranes (100 mM choline chloride, 10 mM HEPES-NaOH, pH 7.4). Blocking kinetics were evaluated from steady-state, single-channel recordings from rat skeletal muscle sodium channels incorporated into planar, neutral phospholipid/decane bilayers (200 mM NaCl, 10 mM HEPES-NaOH, pH 7.0). Blocking events generally appeared as periods of seconds to minutes in which current through the single channel was completely eliminated. A notable exception was seen for the substitution Arg-13-Gln for which the "blocked" events showed measurable conductances of about 20-40% of the open state. The substitution of Arg-13 reduced binding to electroplax membranes to undetectable levels and increased the apparent dissociation constant determined for skeletal muscle channels by greater than 80-fold compared with the native peptide. Other substitutions caused smaller decreases in affinity. The decreased potency of the toxin derivatives resulted both from increases in the rates of dissociation from the channel, and from decreases in association rates. Our data support the suggestion by Sato et al. [(1991) J. Biol. Chem. 265, 16989] that Arg-13 associates intimately with the binding site on the channel. In addition, our results suggest that certain residues affect almost exclusively the approach and docking of the toxin with its binding site, others appear to be important only to the strength of the association once binding has taken place, and yet others affect both.     

Dissecting the energetics of an antibody-antigen interface by alanine shaving and molecular grafting.

Alanine-scanning mutagenesis on human growth hormone (hGH) identified 5 primary determinants (Arg 8, Asn 12, Arg 16, Asp 112, and Asp 116) for binding to a monoclonal antibody (MAb 3) (Jin L, Fendly BM, Wells JA, 1992, J Mol Biol 226:851-865). To further analyze the energetic importance of residues surrounding these five, we mutated all neighboring residues to alanine in groups of 7-16 (a procedure we call alanine shaving). Even the most extremely mutated variant, with 16 alanine substitutions, caused less than a 10-fold reduction in binding affinity to MAb3. By comparison, mutating any 1 of the 5 primary determinants to alanine caused a 6- to > 500-fold reduction in affinity. Replacing any of the 4 charged residues (Arg 8, Arg 16, Asp 112, and Asp 116) with a homologous residue (i.e., Arg to Lys or Asp to Glu) caused nearly as large a reduction in affinity as the corresponding alanine replacement. It was possible to graft the 5 primary binding determinants onto a nonbinding homologue of hGH, human placental lactogen (hPL), which has 86% sequence identity to hGH. The grafted hPL mutant bound 10-fold less tightly than hGH to MAb3 but bound as well as hGH when 2 additional framework mutations were introduced. Attempts to recover binding affinity by grafting the MAb3 epitope onto more distantly related scaffolds having a similar 4-helix bundle motif, such as human prolactin (23% sequence identity) or granulocyte colony-stimulating factor, were unsuccessful.(ABSTRACT TRUNCATED AT 250 WORDS)     

Adenosine inhibits calcium channel currents via A1 receptors on salamander retinal ganglion cells in a mini-slice preparation

The effects of adenosine on high-voltage-activated calcium channel currents in tiger salamander retinal ganglion cells were investigated in a mini-slice preparation. Adenosine produced a concentration-dependent decrease in the amplitude of calcium channel current with a maximum inhibition of 26%. The effects of adenosine on calcium channel current were both time- and voltage-dependent. In cells dialyzed with GTP-gamma-s, adenosine caused a sustained and irreversible inhibition of calcium channel current, suggesting involvement of a GTP-binding protein. The inhibitory effect of adenosine on calcium channel current was blocked by the A1 antagonist 8-cyclopentyltheophylline (DPCPX, 1-10 microm), but not by the A2 antagonist 3-7-dimethyl-1-propargylxanthine (DMPX, 10 microm), and was mimicked by the A1 agonist N6-cyclohexyladenosine (CHA, 1 microm) but not by the A2 agonist 5'-(N-cyclopropyl) carbox-amidoadenosine (CPCA, 1 microm). Adenosine's inhibition of calcium channel current was not affected by the L-type calcium channel blocker nifedipine (5 microm). However, adenosine's inhibition of calcium channel current was reduced to approximately 10% after application of omega-conotoxin GVIA (1 microm), suggesting that adenosine inhibits N-type calcium channels. These results show that adenosine acts on an A1 adenosine receptor subtype via a G protein-coupled pathway to inhibit the component of calcium channel current carried in N-type calcium channels.     

BmBKTx1, a novel Ca2+-activated K+ channel blocker purified from the Asian scorpion Buthus martensi Karsch

BmBKTx1 is a novel short chain toxin purified from the venom of the Asian scorpion Buthus martensi Karsch. It is composed of 31 residues and is structurally related to SK toxins. However, when tested on the cloned rat SK2 channel, it only partially inhibited rSK2 currents, even at a concentration of 1 microm. To screen for other possible targets, BmBKTx1 was then tested on isolated metathoracic dorsal unpaired median neurons of Locusta migratoria, in which a wide variety of ion channels are expressed. The results suggested that BmBKTx1 could specifically block voltage-gated Ca(2+)-activated K(+) currents (BK-type). This was confirmed by testing the BmBKTx1 effect on the alpha subunits of BK channels of the cockroach (pSlo), fruit fly (dSlo), and human (hSlo), heterologously expressed in HEK293 cells. The IC(50) for channel blocking by BmBKTx1 was 82 nm for pSlo and 194 nm for dSlo. Interestingly, BmBKTx1 hardly affected hSlo currents, even at concentrations as high as 10 microm, suggesting that the toxin might be insect specific. In contrast to most other scorpion BK blockers that also act on the Kv1.3 channel, BmBKTx1 did not affect this channel as well as other Kv channels. These results show that BmBKTx1 is a novel kind of blocker of BK-type Ca(2+)-activated K(+) channels. As the first reported toxin active on the Drosophila Slo channel dSlo, it will also greatly facilitate studying the physiological role of BK channels in this model organism.     

Low Voltage Activation of KCa1.1 Current by Cav3-KCa1.1 Complexes

Calcium-activated potassium channels of the KCa1.1 class are known to regulate repolarization of action potential discharge through a molecular association with high voltage-activated calcium channels. The current study examined the potential for low voltage-activated Cav3 (T-type) calcium channels to interact with KCa1.1 when expressed in tsA-201 cells and in rat medial vestibular neurons (MVN) in vitro. Expression of the channel α-subunits alone in tsA-201 cells was sufficient to enable Cav3 activation of KCa1.1 current. Cav3 calcium influx induced a 50 mV negative shift in KCa1.1 voltage for activation, an interaction that was blocked by Cav3 or KCa1.1 channel blockers, or high internal EGTA. Cav3 and KCa1.1 channels coimmunoprecipitated from lysates of either tsA-201 cells or rat brain, with Cav3 channels associating with the transmembrane S0 segment of the KCa1.1 N-terminus. KCa1.1 channel activation was closely aligned with Cav3 calcium conductance in that KCa1.1 current shared the same low voltage dependence of Cav3 activation, and was blocked by voltage-dependent inactivation of Cav3 channels or by coexpressing a non calcium-conducting Cav3 channel pore mutant. The Cav3-KCa1.1 interaction was found to function highly effectively in a subset of MVN neurons by activating near –50 mV to contribute to spike repolarization and gain of firing. Modelling data indicate that multiple neighboring Cav3-KCa1.1 complexes must act cooperatively to raise calcium to sufficiently high levels to permit KCa1.1 activation. Together the results identify a novel Cav3-KCa1.1 signaling complex where Cav3-mediated calcium entry enables KCa1.1 activation over a wide range of membrane potentials according to the unique voltage profile of Cav3 calcium channels, greatly extending the roles for KCa1.1 potassium channels in controlling membrane excitability.     

Custom distinctions in the interaction of G-protein beta subunits with N-type (CaV2.2) versus P/Q-type (CaV2.1) calcium channels

Inhibition of N- (Cav2.2) and P/Q-type (Cav2.1) calcium channels by G-proteins contribute importantly to presynaptic inhibition as well as to the effects of opiates and cannabinoids. Accordingly, elucidating the molecular mechanisms underlying G-protein inhibition of voltage-gated calcium channels has been a major research focus. So far, inhibition is thought to result from the interaction of multiple proposed sites with the Gbetagamma complex (Gbetagamma). Far less is known about the important interaction sites on Gbetagamma itself. Here, we developed a novel electrophysiological paradigm, "compound-state willing-reluctant analysis," to describe Gbetagamma interaction with N- and P/Q-type channels, and to provide a sensitive and efficient screen for changes in modulatory behavior over a broad range of potentials. The analysis confirmed that the apparent (un)binding kinetics of Gbetagamma with N-type are twofold slower than with P/Q-type at the voltage extremes, and emphasized that the kinetic discrepancy increases up to ten-fold in the mid-voltage range. To further investigate apparent differences in modulatory behavior, we screened both channels for the effects of single point alanine mutations within four regions of Gbeta1, at residues known to interact with Galpha. These residues might thereby be expected to interact with channel effectors. Of eight mutations studied, six affected G-protein modulation of both N- and P/Q-type channels to varying degrees, and one had no appreciable effect on either channel. The remaining mutation was remarkable for selective attenuation of effects on P/Q-, but not N-type channels. Surprisingly, this mutation decreased the (un)binding rates without affecting its overall affinity. The latter mutation suggests that the binding surface on Gbetagamma for N- and P/Q-type channels are different. Also, the manner in which this last mutation affected P/Q-type channels suggests that some residues may be important for "steering" or guiding the protein into the binding pocket, whereas others are important for simply binding to the channel.     

RNA recognition mechanism of the minimal active domain of the human immunodeficiency virus type-2 nucleocapsid protein

NCp8 of HIV-2 contains two CCHC-type zinc fingers connected by a linker, and is involved in many critical steps of the virus life cycle. It was previously shown that the first zinc finger flanked by the linker is the minimal active domain for specific binding to viral RNA. In our previous study, we determined the three-dimensional structure of NCp8-f1, including the minimal active domain, and found that a hydrogen bond between Asn(11) N(delta)H and Arg(27) O stabilized the conformation of the linker in the vicinity of the zinc finger [Kodera et al. (1998) Biochemistry 37, 17704-17713]. In this study, RNA binding activities of NCp8-f1 and three types of its mutant peptides were analysed by native PAGE assay. The activity and three-dimensional structure of NCp8-f1/N11A, in which alanine is substituted for Asn(11) thereby affecting the conformation of the linker, was analyzed and compared with those of NCp8-f1. We demonstrated that the existence of Arg(4) and/or Lys(5) and Arg(26) and/or Arg(27) were necessary for binding RNA. Furthermore, the linker's flexible orientation, which is controlled by the hydrogen bond between Asn(11) N(delta)H and Arg(27) O, appears to be a structural basis for NCp8 existing as a multi-functional protein.     

Early response of cultured lepidopteran cells to exposure to delta-endotoxin from Bacillus thuringiensis: involvement of calcium and anionic channels.

The role of ion channels in the initial steps following exposure of SF-9 lepidopteran insect cells in culture to the delta-endotoxin CryIC from the insecticidal bacterium Bacillus thuringiensis was investigated using single ionic channel measurements and microspectrofluorescence of the calcium-sensitive probe fura-2. It was found that: (1) the toxin triggers an immediate rise in intracellular calcium; (2) the surge is due to calcium entering the cells via calcium channels; (3) the toxin recruits or introduces anionic channels in the cell's plasma membrane in a time-dependent manner. These channels, not seen in the absence of the toxin, are induced by toxin exposure to either side of the cell membrane. They have a conductance of 26 picosiemens (pS) and are mainly permeable to chloride. This study provides the first evidence of the primary role of calcium and chloride ions in the action of delta-endotoxin on cultured insect cells.     

Brownian dynamics simulations of the recognition of the scorpion toxin maurotoxin with the voltage-gated potassium ion channels

The recognition of the scorpion toxin maurotoxin (MTX) by the voltage-gated potassium (Kv1) channels, Kv1.1, Kv1.2, and Kv1.3, has been studied by means of Brownian dynamics (BD) simulations. All of the 35 available structures of MTX in the Protein Data Bank (http://www.rcsb.org/pdb) determined by nuclear magnetic resonance were considered during the simulations, which indicated that the conformation of MTX significantly affected both the recognition and the binding between MTX and the Kv1 channels. Comparing the top five highest-frequency structures of MTX binding to the Kv1 channels, we found that the Kv1.2 channel, with the highest docking frequencies and the lowest electrostatic interaction energies, was the most favorable for MTX binding, whereas Kv1.1 was intermediate, and Kv1.3 was the least favorable one. Among the 35 structures of MTX, the 10th structure docked into the binding site of the Kv1.2 channel with the highest probability and the most favorable electrostatic interactions. From the MTX-Kv1.2 binding model, we identified the critical residues for the recognition of these two proteins through triplet contact analyses. MTX locates around the extracellular mouth of the Kv1 channels, making contacts with its beta-sheets. Lys23, a conserved amino acid in the scorpion toxins, protrudes into the pore of the Kv1.2 channel and forms two hydrogen bonds with the conserved residues Gly401(D) and Tyr400(C) and one hydrophobic contact with Gly401(C) of the Kv1.2 channel. The critical triplet contacts for recognition between MTX and the Kv1.2 channel are Lys23(MTX)-Asp402(C)(Kv1), Lys27(MTX)-Asp378(D)(Kv1), and Lys30(MTX)-Asp402(A)(Kv1). In addition, six hydrogen-bonding interactions are formed between residues Lys23, Lys27, Lys30, and Tyr32 of MTX and residues Gly401, Tyr400, Asp402, Asp378, and Thr406 of Kv1.2. Many of them are formed by side chains of residues of MTX and backbone atoms of the Kv1.2 channel. Five hydrophobic contacts exist between residues Pro20, Lys23, Lys30 and Tyr32 of MTX and residues Asp402, Val404, Gly401, and Arg377 of the Kv1.2 channel. The simulation results are in agreement with the previous molecular biology experiments and explain the binding phenomena between MTX and Kv1 channels at the molecular level. The consistency between the results of the BD simulations and the experimental data indicated that our three-dimensional model of the MTX-Kv1.2 channel complex is reasonable and can be used in additional biological studies, such as rational design of novel therapeutic agents blocking the voltage-gated channels and in mutagenesis studies in both the toxins and the Kv1 channels. In particular, both the BD simulations and the molecular mechanics refinements indicate that residue Asp378 of the Kv1.2 channel is critical for its recognition and binding functionality toward MTX. This phenomenon has not been appreciated in the previous mutagenesis experiments, indicating this might be a new clue for additional functional study of Kv1 channels.     

Unique Mechanism of the Interaction between Honey Bee Toxin TPNQ and rKir1.1 Potassium Channel Explored by Computational Simulations: Insights into the Relative Insensitivity of Channel towards Animal Toxins

Background

The 21-residue compact tertiapin-Q (TPN_Q) toxin, a derivative of honey bee toxin tertiapin (TPN), is a potent blocker of inward-rectifier K⁺ channel subtype, rat Kir1.1 (rKir1.1) channel, and their interaction mechanism remains unclear.

Principal Findings

Based on the flexible feature of potassium channel turrets, a good starting rKir1.1 channel structure was modeled for the accessibility of rKir1.1 channel turrets to TPN_Q toxin. In combination with experimental alanine scanning mutagenesis data, computational approaches were further used to obtain a reasonable TPN_Q toxin-rKir1.1 channel complex structure, which was completely different from the known binding modes between animal toxins and potassium channels. TPN_Q toxin mainly adopted its helical domain as the channel-interacting surface together with His12 as the pore-blocking residue. The important Gln13 residue mainly contacted channel residues near the selectivity filter, and Lys20 residue was surrounded by a polar “groove” formed by Arg118, Thr119, Glu123, and Asn124 in the channel turret. On the other hand, four turrets of rKir1.1 channel gathered to form a narrow pore entryway for TPN_Q toxin recognition. The Phe146 and Phe148 residues in the channel pore region formed strong hydrophobic protrusions, and produced dominant nonpolar interactions with toxin residues. These specific structure features of rKir1.1 channel vestibule well matched the binding of potent TPN_Q toxin, and likely restricted the binding of the classical animal toxins.

Conclusions/Significance

The TPN_Q toxin-rKir1.1 channel complex structure not only revealed their unique interaction mechanism, but also would highlight the diverse animal toxin-potassium channel interactions, and elucidate the relative insensitivity of rKir1.1 channel towards animal toxins.     

Study of the binding residues between ANEPII and insect sodium channel receptor

The present study aimed at determining the functional characteristics of anti-neuroexcitation peptide II (ANEPII). The depressant insect toxin ANEPII from the Chinese scorpion Buthus martensii Karsch had an effect on insect sodium channels. Previous studies showed that scorpion depressant toxins induce insect flaccid paralysis upon binding to receptor site-4, so we tried to predict the functional residues involved using computational techniques. In this study, three-dimensional structure modeling of ANEPII and site-4 of the insect sodium channel were carried out by homology modeling, and these models were used as the starting point for nanosecond-duration molecular dynamics simulations. Docking studies of ANEPII in the sodium channel homology model were conducted, and likely ANEPII binding loci were investigated. Based on these analyses, the residues Tyr34, Trp36, Gly39, Leu40, Trp53, Asn58, Gly61 and Gly62 were predicted to interact with sodium channel receptor and to act as functional residues.