![]() |
|
|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Institut de Pharmacologie Moléculaire et Cellulaire Centre National de la Recherche Scientifique Unité Mixte de Recherche 6097, Valbonne, France (L.R., S.D., M.L., P.E.); and Laboratory of Toxicology, University of Leuven, Leuven, Belgium (F.B., S.Z., J.T.)
Received June 21, 2005; accepted November 1, 2005
| Abstract |
|---|
|
|
|---|
1, Nav1.2/
1, Nav1.3/
1, Nav1.4/
1, Nav1.5/
1, and Nav1.8/
1). These novel toxins modulate voltage-gated sodium channels with properties similar to those of typical gating-modifier toxins, both by causing a depolarizing shift in gating kinetics and by blocking the inward component of the sodium current. PaurTx3 is one of the most potent peptide modulators of voltage-gated sodium channels described thus far from spider venom, modulating Nav1.2 with an IC50 value of 0.6 ± 0.1 nM. CcoTx1 and CcoTx2, differing by only one amino acid, are potent modulators of different voltage-gated sodium channel subtypes from the central nervous system, except for Nav1.3, which is only affected by CcoTx2. The potency of CcoTx3 is lower, although this toxin seems to be more selective for the tetrodotoxin-resistant channel subtype Nav1.5/
1 (IC50 = 447 ± 32 nM). In addition to these results, molecular modeling indicates that subtle differences in toxin surfaces may relate to their different pharmacological profiles. Furthermore, an evolutionary trace analysis of these toxins and other structurally related three-disulfide spider toxins provides clues for the exploration of toxin-channel interaction and future structure-function research.
subunit associated with up to four known different
subunits. The
subunits are classified according to sequence homology as Nav1.1 to Nav1.9 and a less defined subunit is called Nax (Yu and Catterall, 2003
-subunits are TTX-sensitive (TTX-s). Subtype localization varies: Nav1.1, Nav1.2, and Nav1.3 are principally found in the central nervous system, whereas Nav1.6, Nav1.7, Nav1.8, and Nav1.9 are mostly distributed in the peripheral nervous system. Nav1.4 can be found in skeletal muscle, and Nav1.5 is present predominantly in cardiac muscle (French and Terlau, 2004
As crucial components of the development of action potentials, VGSCs are one of the foremost targets of animal venoms or plant neurotoxins (Possani et al., 1999
; Bosmans et al., 2002
; Lewis and Garcia, 2003
; French and Terlau, 2004
). Toxins from various organisms have been used to describe eight different receptor sites on the
subunit of VGSCs, all of which are linked to specific effects on channel function (Wang and Wang, 2003
). However, toxin characterization is often limited to whole-cell sodium currents and binding studies on neuronal membranes. For most of these toxins, the precise pattern of subtype selectivity is either unknown or at best fragmentary.
In the past 2 years, a number of novel spider toxins have been demonstrated to modulate VGSCs (Escoubas and Rash, 2004
). Although the overall structure of these toxins, based on the inhibitory cystine knot (ICK) fold, is fairly similar, their electrophysiological properties vary greatly (Escoubas et al., 2000b
; Rash and Hodgson, 2002
; Sollod et al., 2005
). Hainantoxins I, II, IV, and V as well as huwentoxin IV do not alter activation or inactivation kinetics (Li et al., 2003
, 2004
; Peng et al., 2002
). They reduce the current amplitude and have therefore been hypothesized to act by occluding the ion conduction pathway, similarly to toxins acting on site 1. Conversely, the protoxins and jingzhaotoxin-III shift the activation voltage toward more positive values, without affecting inactivation (Middleton et al., 2002
; Xiao et al., 2004
). The protoxins inhibit several VGSC subtypes but can also hamper the function of selected voltage-gated calcium or potassium channels. At this moment, the exact site of action of these toxins on the VGSC
subunit remains the object of speculation. It is thought that they bind to DII-S3/DII-S4 near site 4.
To find more selective modulators of different VGSCs subtypes, we have screened a large number of tarantula venoms on six cloned VGSC subtypes and investigated the activity of the four most potent toxins purified from two of those venoms. To examine potential promiscuous action of the toxins, we also tested them on mouse DRG neurons as a source of voltage-gated calcium channels and on Kv1.3 as a representative voltage-gated potassium channel. Mouse central injections allowed us to examine their overall neurotoxicity pattern. To complement our pharmacological analysis, homology modeling and an evolutionary trace (ET) analysis of these new peptides and other structurally related spider toxins were performed. The results reported here offer new insights into the mode of action of this family of VGSC modulators and the associated structure-function relationships.
| Materials and Methods |
|---|
|
|
|---|
20% of volume) were redissolved in distilled water to 10 times the initial venom volume (1:10 dilution), centrifuged (14,000 rpm, 20 min), filtered on 0.45-µm microfilters (Millipore Corporation, Billerica, MA) and stored at 20°C. Toxin Purification. A total of 110 µl of crude C. cornuatus venom and 10 µl of crude P. auratus venom was fractionated by C8 reversed-phase semipreparative HPLC (10 x 250 mm; Nacalai Tesque, Kyoto, Japan) using a linear gradient of water (A)/acetonitrile (B) in constant 0.1% trifluoroacetic acid (015% B in 15 min, 1550% in 70 min, 5060% in 10 min, 6090% in 10 min, 2 ml/min). Fractions were hand-collected by monitoring the effluent signal (215 nm) and dried in a vacuum centrifuge. Aliquots of each fraction were assayed for activity against cloned VGSCs expressed in Xenopus laevis oocytes (see below). Active fractions were further separated by cation-exchange chromatography on an SP5PW column (4.6 x 70 mm) (TOSOH, Tokyo, Japan), with a linear gradient of ammonium acetate in water from 20 mM (A) to 2 M (B) (0% B for 10 min, 060% in 60 min, 6095% B in 10 min, 1 ml/min, detection at 280 nm). A third purification step was conducted for each active peak using a C4 reversed-phase column (Develosil 300C4, 4.6 x 250 mm, Nomura Chemical, Seto, Japan) and a linear gradient of acetonitrile in water and 0.1% TFA (0% B for 5 min, 040% in 80 min, 4090% in 10 min, 1 ml/min). All solvents used were of HPLC grade. Separations were conducted on a HP1100 system (Hewlett Packard, Palo Alto, CA) coupled to a diode-array detector.
Toxin Characterization. For N-terminal sequence determination, peptides were reduced (55 mM dithiothreitol, 60°C, 1 h) and alkylated with iodoacetamide (RT, 30 min), desalted by reversed-phase high performance liquid chromatography (RP-HPLC) on a Purospher STAR column (C18, 55 x 4 mm, 3 µm; Merck, Whitehouse Station, NJ) (water/acetonitrile/0.1% TFA, 040% B in 20 min, 1 ml/min) and submitted to automated N-terminal sequencing on a gas-phase sequencer (model 477A; Applied Biosystems, Foster City, CA).
For confirmation of the sequences and sequencing of the C terminus, reduced and alkylated peptides were digested with trypsin (500 pmol of peptide, trypsin ratio of 1:50 (w/w) in 25 mM NH4HCO3 buffer, pH 8.3, 12 h at 37°C). After quenching the reaction with 0.1% TFA in water, an aliquot was analyzed by MALDI-TOF mass spectrometry, and the mixture of proteolytic peptides was separated by reversed-phase HPLC on the same analytical column (water/acetonitrile/0.1% TFA, 060% B in 60 min, 1 ml/min). Each fraction was then analyzed by MALDI-TOF mass spectrometry. Peptides corresponding to the unknown part of the sequence according to mass calculations done in GPMAW (http://welcome.to/GPMAW) were submitted to automated Edman sequencing.
Sequence similarities were determined by a search of nonredundant protein databases, via the BLAST server (http://www.ncbi.nlm.nih.gov/) and by comparison with sequences obtained from literature data and compiled in an in-house spider peptide toxin database. Sequence alignments were calculated using ClustalX 1.8.
Mass Spectrometry. Peptides were analyzed by MALDI-TOF mass spectrometry on an Applied Biosystems Voyager DE-PRO system, in reflector mode using recrystallized
-cyano-4-hydroxycinnamic acid matrix. Mass spectra were calibrated with internal peptide standards and analyzed in the Data Explorer software. Masses were also measured by electrospray ionization on a Thermo Finnigan LCQ DecaXP electrospray ion trap mass spectrometer (Thermo Electron Corporation, Waltham, MA).
Mouse Intracerebroventricular Injections. Activity against vertebrates was evaluated by i.c.v. injection in mice. C57/Bl6 mice were briefly anesthetized with diethyl ether, injected in the left cerebral ventricle with 10 µl of sample (500 pmol of toxin in 10 µl of double-distilled H2O) or control (10 µl double-distilled H2O) and placed in glass jars for observation. Mice were continuously monitored for symptoms of neurotoxicity during the first hour after injection or until death. All the experiments involving the use of animals complied with University Ethics Regulations.
DRG Preparation and Patch-Clamp Recording of Calcium Currents. Dorsal root ganglia (DRG) were dissected from adult mice (710 weeks) and enzymatically dissociated with collagenase (0.2% type II; Worthington Biochemicals, Freehold, NJ) and trypsin (2.5 mg/ml; Seromed, Berlin, Germany). Cells were plated on 35-mm Petri dishes coated with poly-(L-lysine) (Sigma, St. Louis, MO) and maintained in culture at 37°C (95% air/5% CO2) in Ham's F12 medium (Invitrogen, Carlsbad, CA) containing 10% fetal calf serum (MP Biomedicals, Irvine, CA) and 1% penicillin/streptomycin (Invitrogen). Electrophysiological experiments were carried out 1 or 2 days after plating.
Currents were sampled at 3.3 kHz for whole-cell recordings using pClamp8 software (Molecular Devices, Sunnyvale, CA). The pipette solution contained 150 mM CsCl, 1 mM MgCl2, 5 mM EGTA, 3 mM ATP, and 10 mM HEPES-CsOH, pH 7.35, and the bath solution contained: 140 mM TEACl, 5 mM BaCl2, 10 mM glucose, and 10 mM HEPES-TEAOH, pH 7.35. Solutions were applied near the neuron using a rapid perfusion system. The holding potential was 80 mV and high voltage-activated calcium currents were elicited by depolarizing steps between 20 and +50 mV. Averaged data are presented as mean ± S.E.M.
Sodium Channel Expression. For expression in X. laevis oocytes, the Nav1.5 and
1 genes were previously subcloned into pSP64T. For in vitro transcription, Nav1.5/pSP64T and Nav1.8/pSP64T were first linearized with XbaI and
1/pSP64T with EcoRI. Capped cRNAs were synthesized from the linearized plasmid using the SP6 mMESSAGE mMACHINE transcription kit (Ambion, Austin, TX). The Nav1.1/pLCT1, Nav1.2/pLCT1, Nav1.3/pNa3T, and Nav1.4/pUI-2 vectors were linearized with NotI and transcribed with the T7 mMESSAGE mMACHINE kit (Ambion) (Kayano et al., 1988
; Smith and Goldin, 1998
; Bosmans et al., 2002
).
Potassium Channel Expression. The vector pCI.neo containing the gene for Kv1.3 was linearized with NotI and transcribed using the large-scale T7 mMESSAGE mMACHINE transcription kit (Ambion) (Swanson et al., 1990
).
Electrophysiological Studies on Cloned Channels. The harvesting of stage VVI oocytes from the ovarian lobes of anesthetized female X. laevis frogs was done as described previously (Bosmans et al., 2002
). Oocytes were injected with 50 nl of cRNA at a concentration of 1 ng/nl using a microinjector from Drummond Scientific (Broomall, PA). The solution used for incubating the oocytes contained 96 mM NaCl, 2 mM KCl, 1.8 mM CaCl2, 2 mM MgCl2, and 5 mM HEPES, pH 7.4, supplemented with 50 mg/l gentamycin sulfate and 180 mg/l theophylline.
Two-electrode voltage-clamp (TEVC) recordings were performed at room temperature (1822°C) using a GeneClamp 500 amplifier (Molecular Devices) controlled by a pClamp data acquisition system (Molecular Devices). Whole-cell currents from oocytes were recorded 2 to 4 days after injection. Voltage and current electrodes were filled with 3 M KCl. Resistances of both electrodes were kept as low as possible (<0.5 M
). Bath solution composition was 96 mM NaCl, 2 mM KCl, 1.8 mM CaCl2, 2 mM MgCl2, and 5 mM HEPES, pH 7.4. Currents were filtered at 1 kHz with a four-pole, low-pass Bessel filter and sampled at 5 kHz. Leak subtraction was performed using a P/4 protocol. Currents were evoked in oocytes expressing the cloned VGSCs by depolarizations between 70 and 40 mV, using 10-mV increments from a holding potential of 90 mV. To avoid overestimation of a potential toxin-induced shift in the current-voltage relationship as a result of inadequate voltage control when measuring large sodium currents in oocytes, only results from cells with currents lower than 1.5 µA were considered in Table 1. To obtain IC50 values on VGSCs, the percentage of toxin-induced block was plotted against the concentration of toxin used and a fit with the Hill equation yielded the IC50 values. The percentage of toxin-induced block was always measured at the same voltage, which also gave the maximum current under control conditions (no toxin: *). This method was preferred over choosing a fixed voltage (e.g., 20 mV) for all VGSCs to avoid a distorted calculation of toxin-induced block of channels that activate at more positive voltages (e.g., Nav1.8). Furthermore, over the tested voltage range (between 20 and +30 mV), the percentage of toxin-induced block, compared among all the VGSC isoforms, was found to be invariable [data derived from Fig. 4 (not shown)]. Kv1.3 currents were evoked by depolarizations to 0 mV from a holding potential of 90 mV. All toxins were tested on at least three oocytes (n
3). Data manipulation was performed in pClamp8 (Molecular Devices) and Origin software (OriginLab Corp., Northampton, MA). Averaged data are presented as mean ± S.E.M.
|
|
|
Evolutionary Trace Analysis. The sequences of 33 inhibitory spider venom peptides conforming to the ICK fold were used in the evolutionary trace analysis (Craik et al., 2001
). These sequences are supplied online as a Supplemental Table. The multiple sequence alignment was performed using the program ClustalX and was manually refined. The ET analysis was carried out based on a phylogenetic tree created with the unweighted pair group method with arithmetic mean. To summarize briefly, the ET method divides all residues of aligned sequences into three classes: neutral, conserved, and class-specific, based on the comparison of consensus sequences for groups of proteins that originate from a common node defined by the evolutionary time cut-off in a phylogenetic tree (for details, see Zhu et al., 2004
). Class-specific trace residues identified were mapped onto the structure of HnTX-IV (PDB code 1NIY
[PDB]
) in the PyMOL molecular modeling program (http://pymol.sourceforge.net).
|
| Results |
|---|
|
|
|---|
20 tarantula venoms, both C. cornuatus and P. auratus venoms consistently displayed high and reproducible inhibitory activity against a variety of cloned VGSCs. The bioassay-guided fractionation of C. cornuatus venom resulted in the identification of three active fractions. Two peptides were isolated from fraction 24 and fraction 25/26 yielded another one, which were named, respectively, CcoTx1, CcoTx2, and CcoTx3) (see Fig. 1A). Bioassay-guided fractionation of P. auratus venom yielded a single active peptide named PaurTx3 (see Fig. 1B). Purity of the toxins (>99%) was assessed in the last two steps of purification using orthogonal methods (ion-exchange followed by reversed-phase), which yielded single, symmetrical peaks, and by MALDI-TOF mass spectrometry, which did not reveal the presence of contaminants. Automated sequencing of all reduced and alkylated peptides yielded almost complete sequences and sequencing of several overlapping tryptic peptides completed the C-terminal part of the sequences. The four toxins are 32 to 39 amino acids long, and they each possess six cysteine residues forming three disulfide bridges (based on molecular mass data). CcoTx1, CcoTx2 and PaurTx3 are basic peptides (respective calculated pI 10.48, 10.07, 10.18) and CcoTx3 is an acidic peptide (pI 6.05). The molecular masses measured by MALDI-TOF mass spectrometry (PaurTx3, 4055.88 Da; CcoTx1, 4041.79 Da; CcoTx2, 4089.61 Da; and CcoTx3, 4321.84 Da) are in perfect accordance with those calculated from the sequence data. They indicate amidation at the C terminus of CcoTx1 (calc. 4041.74 Da), CcoTx2 (calc. 4089.78 Da), and CcoTx3 (calc. 4321.85 Da) (1 Da difference with the molecular mass calculated with a carboxylic C terminus). Likewise, mass spectrometry indicates a free carboxyl-terminal form for PaurTx3 (calc. 4055.87 Da).
CcoTx1, CcoTx2, and PaurTx3 are highly similar toxins; CcoTx1 and CcoTx2 are isoforms that are different only at position 32, with an aspartic acid residue in CcoTx1 and a tyrosine in CcoTx2. CcoTx3 shows a completely different primary sequence. The pairing of cysteine residues was not determined experimentally in the present study, because the very high similarity of the ceratotoxins and phrixotoxin 3 with toxins described previously strongly suggests that they conform to the canonical motif CICIV, CIICV, CIIICVI found in all members of this structural group (Escoubas and Rash, 2004
). Our sequences were analyzed with CysView, a web-based tool that identifies and classifies proteins according to their disulfide connectivity patterns (http://research.i2r.a-star.edu.sg/CysView). The results returned by the program suggested that CcoTx1, CcoTx2, CcoTx3, and PaurTx3 indeed adhere to the aforementioned disulfide bridge canonical motif.
Sequence Similarity. CcoTx1, CcoTx2, and PaurTx3 share a high similarity with other peptides previously isolated from other tarantula venoms (Fig. 1C). The most similar toxins are Huwentoxin IV (74%, 74%, and 77%) and Hainantoxin IV (77%, 75%, and 82%), VGSC inhibitors from the Chinese tarantulas Ornithoctonus huwena and Selenocosmia hainana, respectively (Peng et al., 2002
; Li et al., 2004
), and GsMTX4 (74% for PaurTx3), an inhibitor of mechanosensitive channels from the venom of the Chilean tarantula Grammostola spatulata (Suchyna et al., 2000
).
CcoTx3 shares a high similarity with toxin SNX482 (76%), an R-type voltage-gated calcium channel blocker from the African tarantula Hysterocrates gigas (Newcomb et al., 1998
) and with the Kv2.1 channel inhibitor hanatoxin 1 (60%) from the venom of G. spatulata (Swartz and MacKinnon, 1995
) (see Fig. 1C).
Mouse Neurotoxicity. Mice injected with 500 pmol of CcoTx1, CcoTx2, and PaurTx3 displayed exactly the same symptoms of neurotoxicity: injection of the toxins was immediately followed by general ataxia, lack of response to stimuli, and semiparalysis. After a few minutes, the mice were unable to stand, and breathing was reduced in rhythm and intensity. Symptoms gradually increased with progressive slowing of breathing and flaccid paralysis; death occurred within 10 to 20 min of injection. Animals remained totally flaccid, and during the course of intoxication, no symptoms of excitatory neurotoxicity (spasms, convulsions, rapid movements) were observed. Conversely, CcoTx3 injection resulted in markedly different symptoms characterized by an erect, elevated tail, initial partial ataxia, followed by recovery over approximately 1 h after injection and the progressive development of shaking. Although paralysis subsided, the body tremors never ceased and persisted until the end of the experiment.
Activity on Calcium Currents. Using the whole cell configuration of the patch-clamp technique on sensory neurons from adult mice, high-voltage activated (HVA) calcium currents were activated between 20 and +50 mV from a holding potential of 80 mV to obtain current-voltage relationships (I-V curves). Perfusion of CcoTx1, CcoTx2, CcoTx3, and PaurTx3 at 100 nM concentration did not inhibit HVA calcium currents evoked at 10 or 0 mV in sensory neurons (2.4 ± 3, 4 ± 4, 10 ± 7, and 2.5 ± 5% inhibition, respectively, n = 4 each) (data not shown).
Activity on Sodium Currents. Using the TEVC technique on Xenopus laevis oocytes, the effects of CcoTx1, CcoTx2, CcoTx3, and PaurTx3 were compared on six different cloned VGSCs coexpressed with the
1 subunit (Nav1.1/
1, Nav1.2/
1, Nav1.3/
1, Nav1.4/
1, Nav1.5/
1, and Nav1.8/
1) (see Figs. 2 and 3, Table 1). Clones for Nav1.6 and Nav1.7 were not available at the time of this study and the Nav1.9 channel currently fails to express in heterologous systems. When measured at the voltage of maximum sodium influx, CcoTx1 decreases Nav1.1, Nav1.2, Nav1.4, and Nav1.5 currents with variable potency. At the maximum concentration tested (2 µM), there is no observable inhibition of Nav1.3 and only
55% inhibition of Nav1.8. CcoTx1 seems to be very highly selective for the neuronal channel Nav1.2, with an IC50 value of 3 ± 1 nM, whereas IC50 values for the other subtypes are at least 2 orders of magnitude higher. The I-V curves of all TTX-s VGSCs are shifted to more positive potentials when CcoTx1 is applied, and voltage of half-maximal activation (V1/2) increases substantially (see Table 1). Nav1.1 is influenced to a lesser degree (V1/2, +1.7 mV), and the I-V curve of Nav1.8 displays no shift at all.
|
40% at a concentration of 2 µM. CcoTx2 is similar to CcoTx1 in that it displays the strongest inhibition on Nav1.2. The most striking result concerning CcoTx1 and CcoTx2 is the dramatic effect of the single mutation of Asp32 to Tyr32 on Nav1.3 activity. Whereas CcoTx1 elicits no observable inhibition of this particular subtype, Nav1.3 current is inhibited by CcoTx2 with an IC50 of 88 ± 25 nM, pointing at the crucial role of this C-terminal residue in channel subtype selectivity.
PaurTx3, the most potent toxin in our study, affects all studied VGSCs with different IC50 values (see Table 1). Here again, Nav1.8 is the least sensitive to the toxin, with a current reduction of
65% at a concentration of 2 µM. PaurTx3 is particularly potent against Nav1.2 (IC50 = 0.6 ± 0.1 nM) and is also by far the most potent toxin for Nav1.5 in our study (IC50 = 72 ± 10 nM). As with CcoTx1 and CcoTx2, PaurTx3 is much less potent against Nav1.1 (IC50 = 610 ± 63 nM) and Nav1.4 (IC50 = 288 ± 58 nM). All I-V curves except the one for Nav1.8 are also shifted to more positive potentials for this toxin (Fig. 4).
No obvious alteration in inactivation rate is observed in all these experiments, and there was also no noticeable shift in reversal potential (Erev), suggesting that toxin interaction with the VGSCs does not cause a change in ion selectivity.
Additional experiments were performed using CcoTx2 to clarify whether these toxins act as pore blockers or gating modifiers. One of the characteristics of gating modifier toxins is that whereas they mediate complete inhibition of all inward currents, outward currents can still be observed despite the presence of the toxin (McDonough et al., 1997
; Bourinet et al., 2001
). CcoTx2 (50 nM) apparently changes the voltage dependence of gating of Nav1.2/
1, so channels do not open in response to the moderate depolarizations that evoke inward current through unblocked channels (Fig. 5A). Nevertheless, channels can still be opened by sufficiently large depolarizations to +120 mV. However, a major limitation for studying sodium currents in oocytes is the presence of endogenous outward currents at very positive voltages. We tried to minimize this effect as shown in Fig. 5B by plotting a trace at +80 mV and subtracting the endogenous currents measured on uninjected oocytes at this voltage (n = 5). As shown in Fig. 5B, sodium current is still present at +80 mV that inactivates over time, and the currents elicited at voltages more positive than the Erev are superimposed in both control and toxin application experiments. To test whether CcoTx2 (50 nM) was still bound to the channel (Nav1.2/
1) after large depolarizations, a test pulse to 10 mV was given 25 ms after a 50-ms depolarization to +120 mV (Fig. 5C). The current at 10 mV was still inhibited by CcoTx2 after the large depolarization, which could be interpreted to mean that CcoTx2 was still bound to the channel. We also observed that the rate of development of inhibition upon CcoTx2 exposure occurred with a time constant of >20 s for this cell, which is slow compared with the 25-ms interval between the pulse to +120 mV and the start of the second pulse to 10 mV. Thus, it is far more likely that the inhibition during the second pulse to 10 mV reflects continuous presence of CcoTx2 on or nearby the channel throughout the pulse sequence rather than toxin unbinding during the pulse to +120 mV followed by rebinding before the following pulse to 10 mV. The trace shown here at +120 mV has not been corrected for the presence of endogenous currents at very positive voltages, and therefore an outward current is seen. In the inset, the deactivation of these endogenous currents is displayed.
|
The application of 2 µM CcoTx3 does not result in current inhibition for Nav1.1, Nav1.2, Nav1.3, or Nav1.4. As with CcoTx1 and CcoTx2, CcoTx3 does reduce Nav1.8 currents by
45%. The only VGSC that is somewhat sensitive to CcoTx3 is the TTX-r Nav1.5, with an IC50 value of 447 ± 32 nM. In addition, CcoTx3 barely shifts the I-V curve of Nav1.5. In contrast, the I-V curve of Nav1.4 was substantially shifted but only at concentrations of 2 µM or more. The rather low potency of CcoTx3 on the tested VGSCs can probably be attributed to its very different primary sequence and suggests either a preferential affinity for other pharmacological targets such as calcium channels or potassium channels or selectivity for VGSCs of a different phylogenetic origin.
Activity on Potassium Currents. Using the TEVC technique on X. laevis oocytes, the effects of CcoTx1, CcoTx2, CcoTx3, and PaurTx3 were studied on Kv1.3, a member of the voltage-gated potassium channel family. Currents were evoked by depolarizations to 0 mV from a holding potential of 90 mV and then clamped back to 50 mV. Neither CcoTx1, CcoTx2, CcoTx3, nor PaurTx3 had any substantial effect on Kv1.3 currents, causing less than 20% block at a concentration of 2 µM (n = 3) (data not shown). Extensive testing on other members of the Kv family could not be undertaken because of insufficient material.
Structural Analysis. In accordance with their primary sequence similarity, examination of the models of the four toxins shows that they indeed adopt an ICK fold motif, with backbones closely matching those of other spider toxins studied by NMR (Craik et al., 2001
; Escoubas and Rash, 2004
). Further examination of the toxin surfaces points out the similarity of the three basic toxins (CcoTx1, CcoTx2, and PaurTx3) with other toxins acting on voltage-gated channels, such as SGTx1, a tarantula toxin inhibiting activation of Kv2.1 channels (Lee et al., 2004
) (see Fig. 6). These toxins seem to be organized into two distinct faces, with the face interacting with the channel bearing a hydrophobic patch surrounded by a crown of charged residues. The importance of this structural organization in gating-modifier ICK toxins was first proposed by Takahashi et al. (2000
) to explain the interaction between hanatoxin and drk1 (Kv2.1). The model was validated in recent mutagenesis studies that singled out the essential role of the hydrophobic protrusion on the surface of SGTx1 (Wang et al., 2004
). This protrusion consists of an "aromatic sandwich" encompassing both hydrophobic and two or three aromatic residues. Sequence alignment shows that similar residues can be found in equivalent positions in loop 4 of the ceratotoxins and PaurTx3 between cysteines IV and V. The three-dimensional models show that they form a similar hydrophobic patch on the surface of the molecule. This area is borne by the first N-terminal
-turn and the first two
-sheets of the toxin, stabilized by the disulfide bridges. Hence, the central hydrophobic patch is essentially formed by residues in the second
-sheet (after Cys IV).
|
In addition, the important role of the single amino acid mutation between CcoTx1 and CcoTx2 demonstrates the crucial importance of the C terminus of the toxin in subtype selectivity. This residue is located on the face of the toxin opposite the hydrophobic patch and therefore highlights the complexity of toxin-channel interaction that does not involve solely the face bearing the hydrophobic patch.
Both the primary sequence and the three-dimensional model of CcoTx3 clearly distinguish it as a different type of ICK toxin. The surface features and in particular the charge distribution of CcoTx3 are very different from those of CcoTx1, CcoTx2, and PaurTx3, clearly reflecting its acidic nature and probable different channel selectivity. In addition, the hydrophobic surface homologous to that of the other ceratotoxins is discontinuous and extends to the C-terminal loop. Its high sequence similarity with SNX482, a blocker of R-type voltage-gated calcium channels, and the results of our electrophysiological study on VGSCs, may point at possible activity against voltage-gated calcium channels. However, no significant inhibition of the HVA calcium currents in DRG neurons was observed. The ultimate target of the peptide might thus be a subtype that is not expressed in this cell model or is responsible only for a minor fraction of the total current. No significant effect was observable on Kv1.3 either.
Evolutionary Trace Analysis. To complete this study, a bioinformatics-based investigation of toxin residues that are potentially important in toxin/channel interaction and functional diversification was undertaken. It was based on the evolutionary trace method (ET), which identifies sequence positions at which variations among related proteins correlate with evolutionary divergence. After ET analysis of 33 sequences with a common ICK motif, the amino acids that were calculated as being of significant importance are highlighted in red in Fig. 7. These residues are mapped onto the structure of HnTX-IV (Li et al., 2004
). Examination of the three-dimensional model makes it clear that almost all evolutionarily important residues, either structurally or pharmacologically, are located on face A of the toxin. Only a leucine at position 3 is located on face B, which is defined by a 180° rotation of face A around the y-axis of the toxin. A remarkable feature is the strong basic node that is present in face A and consists of an arginine (Arg26) and two lysines (Lys27 and Lys32), among which Lys32 is the only significantly important residue located in a
-sheet. Furthermore, it seems that the final loop between the two antiparallel
-sheets is an important functional surface, because four of the six residues in this loop (Ser25, Arg26, Lys27, and Trp30) are highlighted by our ET analysis.
|
| Discussion |
|---|
|
|
|---|
Toxins as Tools to Study VGSCs. Both CcoTx1 and CcoTx2 are able to inhibit Nav1.2 inward currents in the low nanomolar range, but display only moderate activity against Nav1.1, Nav1.4, Nav1.5, and low activity against Nav1.8 (Figs. 2, 3, and 5 and Table 1). Their pharmacological profile indicates that these peptides, when used at concentrations near their IC50, could be used to selectively characterize Nav1.2 currents in neurons in which Nav1.2 is coexpressed with Nav1.1 and Nav1.3. While CcoTx2 is an effective modulator of both Nav1.2 (IC50 = 8 ± 1 nM) and Nav1.3 (IC50 = 88 ± 25 nM), CcoTx1, which differs by only one C-terminal residue, does not inhibit Nav1.3 currents at concentrations up to 2 µM. Not only will the combined application of both toxins permit observation of Nav1.1-related currents only, their separate application should also allow discrimination of Nav1.2 from Nav1.3.
With even higher affinity toward Nav1.2, PaurTx3 is a complementary tool for the study of this channel (Figs. 2, 3, and 4 and Table 1). The overall pharmacological profile of PaurTx3 is similar to that of CcoTx2. It inhibits Nav1.2 (IC50 = 0.6 ± 0.1 nM), and it is also the most potent inhibitor of Nav1.5 in our study.
Our results demonstrate the crucial role of the C-terminal part of these peptides in their interaction with the VGSC and especially in recognizing Nav1.3. For the toxins reported here, the presence of either an acidic group (Asp32 in CcoTx1) or an aromatic one (Tyr32 in CcoTx2) dramatically changes the selectivity through a minor change in its overall structure and surface charge (Fig. 6). Similarly to CcoTx2, the C-terminal acidic residue of CcoTx1 (Asp32) is mutated to a polar uncharged residue (Gln32) in PaurTx3. As in CcoTx2, this modification results in affinity toward Nav1.3 and seems to be critical in the recognition of this particular channel, perhaps by the formation of a hydrogen bond network or the presence of the
-electron cloud of the Tyr32 ring.
Although discovered through the same process of bioassay-guided purification, CcoTx3 has the least activity against the tested VGSCs (Figs. 2 and 3 and Table 1). Examination of its primary structure reveals a sequence totally different from that of the other ceratotoxins but closely related to SNX482 (Fig. 1C). It could be hypothesized that CcoTx3 has a different, high-affinity primary target. The bias in activity-guided purification can be explained by the high abundance of the toxin in the venom and possible "target promiscuity" as shown previously for other ion channel spider toxins. It has been demonstrated previously that at elevated concentrations (>500 nM), ion channel selectivity can be modified, with recognition of conserved motifs such as the voltage-sensing domain in other channel families (Li-Smerin and Swartz, 1998
). Nevertheless, under certain conditions, CcoTx3 might still be of use in the study of the VGSC family as a selective ligand of moderate affinity for Nav1.5.
Modulation of Channel Function. Application of CcoTx1, CcoTx2, and PaurTx3 results in a shift of the current-voltage relationship of the studied VGSCs (Figs. 4 and 5 and Table 1), a characteristic that has also been reported for related spider toxins such as protoxin I, protoxin II, and Jingzhaotoxin III (Middleton et al., 2002
; Xiao et al., 2004
; Smith et al., 2005
). In those studies, it was assumed that these toxins bind to the extracellular linker between DII-S3 and DII-S4 of the VGSC. A similar modification of ion channel activation was also observed in voltage-gated calcium channels with the spider toxins
-agatoxin-IVA,
-grammotoxin-SIA, and the scorpion toxin kurtoxin (Li-Smerin and Swartz, 1998
; Winterfield and Swartz, 2000
). These toxins have thus been classified as gating modifiers. However, other related ICK toxins isolated from tarantula venoms, such as hainantoxins I, III, IV, V, and huwentoxin IV have been reported to entirely block sodium currents without affecting either activation or inactivation. They have been classified as pore blockers, acting by occlusion of the channel, in a manner similar to that of site 1 toxins (Li et al., 2003
, 2004
; Xiao and Liang, 2003
). These studies raise the intriguing possibility of different binding sites and modes of action on the VGSC for these ICK spider toxins.
Our experiments on Nav1.2 using CcoTx2 (Fig. 5) reveal a complete inhibition of all inward currents, but outward currents can still be observed despite the presence of the toxin. Furthermore, as seen in Fig. 5C it seems that CcoTx2 stays continuously bound to or is near the channel even at very positive voltages (+120 mV). In addition, CcoTx1, CcoTx2, and PaurTx3 affect neither channel inactivation nor reversal potential. Compared with certain voltage-gated calcium channel toxins, which have been characterized as gating modifiers, our toxins are quite similar in their mode of action. Both
-grammotoxin-SIA (McDonough et al., 1997
) and SNX482 (Bourinet et al., 2001
) completely block inward calcium currents. However, they do not inhibit outward currents. Considering these results, we suggest that the toxins described in the present report can also be generally classified as gating modifiers.
On the other hand, we show that the ceratotoxins and PaurTx3 do not cause a shift of the voltage-current relationship in Nav1.8 (Fig. 4) but rather seem to act in a manner resembling that of pore blockers. One simple hypothesis for this divergent behavior could be that Nav1.8, which already has a much more positive activation potential by itself, does not possess a high sequence similarity with other VGSCs, and therefore lacks a conserved voltage-sensing domain (Li-Smerin and Swartz, 1998
; Winterfield and Swartz, 2000
). A concomitant supporting fact is that Nav1.8 is also insensitive to the "classic" toxins from scorpions, spider (F. Bosmans and J. Tytgat, unpublished data), and sea anemones (Bosmans et al., 2002
).
Because of fragmented data in the literature, it remains difficult to properly compare the "activation blockers" of VGSCs. Both the protoxins and Jingzhaotoxin III have only been reported to block channel activation (Middleton et al., 2002
; Xiao et al., 2004
). The decrease of inward currents or the binding site of these toxins has not been thoroughly described. A way of classifying our toxins other than as gating modifiers could therefore be as "voltage-dependent blockers", which block only inward sodium currents in a voltage-dependent manner. All the above actually suggests a slightly different mode of action for these peptides and certainly warrants further investigation.
A recent study focused on the molecular surface of tarantula toxins interacting with the voltage sensor of voltage-gated potassium channels (Wang et al., 2004
). The authors stated that SGTx1 interacts with the lipid membrane via its hydrophobic protrusion with a crown of surrounding residues containing important residues (Arg3, Arg22) for contacting the voltage sensor. A similar mechanism of partitioning in the lipid membrane was also recently described for ProTx-II (Smith et al., 2005
). CcoTx1, CcoTx2, and PaurTx3 possess a similar hydrophobic protrusion. This could mean that these toxins interact with the lipid membrane too as they modulate the VGSC, because effects appear within 20 to 50 s of application. This mechanism of action would in turn suggest the possibility that CcoTx2 does dissociate from the channel at high voltages (Fig. 5C) but stays close to the VGSC and remains bound to the membrane so that it can bind again thereafter.
ET Analysis. The ET analysis of related ICK toxin sequences outlines the potential importance of basic residues and particularly of amino acids in loops 1 and 4 (Fig. 7). When linking this theoretical model with the result of the mutant cycle analysis of SGTx1, an excellent correlation is seen (Wang et al., 2004
). Mutagenesis of SGTx1 outlined the crucial importance of residues in loops 1 and 4, in particular Arg3, Arg22, and Asp24. These residues are all singled out by our ET analysis. Mutagenesis, however, points to the importance of aromatic residues in loops 1 and 4 in defining the hydrophobic contact surface of the toxin, whereas the ET analysis seems to be biased toward the more basic residues. This is exemplified by the strong basic node that is present in face A of the toxin. Consistent with the hydrophobic patch contact hypothesis, our ET analysis indicates that face A of the toxins encompasses all important residues, whereas face B includes a leucine at position 3 that can make a pharmacological or structural contribution. Clearly, an ET analysis based solely on a combination of primary sequences of peptides is not sufficient to completely predict the crucial determinants of toxin function. However, ET analysis can be useful in orienting mutagenesis studies in toxins from which little is known regarding their interaction surface.
| Acknowledgements |
|---|
1 subunit, and E. Cuypers for helpful discussions. | Footnotes |
|---|
Accession numbers for Swiss-Prot database are p84507, p84508, p84509, and p84510 for CcoTx1, CcoTx2, CcoTx3, and PaurTx3, respectively.
Article, publication date, and citation information can be found at http://molpharm.aspetjournals.org.
ABBREVIATIONS: VGSC, voltage-gated sodium channel; TTX, tetrodotoxin; TTX-r, tetrodotoxin-resistant; TTX-s, tetrodotoxin-sensitive; ICK, inhibitory cystine knot; DRG, dorsal root ganglion; RP, reversed phase; HPLC, high-performance liquid chromatography; TFA, trifluoroacetic acid; MALDI-TOF, matrix-assisted laser desorption ionization/time of flight; TEVC, two-electrode voltage-clamp; PDB, Protein Data Bank; HnTx, hainantoxin; PaurTx, phrixotoxin; CcoTx, ceratotoxin; I-V, current-voltage; SNX482, Ca2+ channel-blocking toxin from Hysterocrates gigas spider venom.
The online version of this article (available at http://molpharm.aspetjournals.org) contains supplemental material. ![]()
1 Current affiliation: Institute of Zoology, Chinese Academy of Sciences, Beijing, P.R. of China. ![]()
Address correspondence to: Jan Tytgat, University of Leuven, Laboratory of Toxicology, Van Evenstraat 4, Leuven 3000, Belgium. E-mail: jan.tytgat{at}pharm.kuleuven.be
| References |
|---|
|
|
|---|
Bourinet E, Stotz SC, Spaetgens RL, Dayanithi G, Lemos J, Nargeot J, and Zamponi GW (2001) Interaction of SNX482 with domains III and IV inhibits activation gating of alpha(1E) (CaV2.3) calcium channels. Biophys J 81: 7988.
Craik DJ, Daly NL, and Waine C (2001) The cystine knot motif in toxins and implications for drug design. Toxicon 39: 4360.[Medline]
Craner MJ, Newcombe J, Black JA, Hartle C, Cuzner ML, and Waxman SG (2004) Molecular changes in neurons in multiple sclerosis: altered axonal expression of Nav1.2 and Nav1.6 sodium channels and Na+/Ca2+ exchanger. Proc Natl Acad Sci USA 101: 81688173.
Escoubas P, De Weille JR, Lecoq A, Diochot S, Waldmann R, Champigny G, Moinier D, Menez A, and Lazdunski M (2000a) Isolation of a tarantula toxin specific for a class of proton-gated Na+ channels. J Biol Chem 275: 2511625121.
Escoubas P, Diochot S, Celerier ML, Nakajima T, and Lazdunski M (2002) Novel tarantula toxins for subtypes of voltage-dependent potassium channels in the Kv2 and Kv4 subfamilies. Mol Pharmacol 62: 4857.
Escoubas P, Diochot S, and Corzo G (2000b) Structure and pharmacology of spider venom neurotoxins. Biochimie 82: 893907.[Medline]
Escoubas P and Rash L (2004) Tarantulas: eight-legged pharmacists and combinatorial chemists. Toxicon 43: 555574.[Medline]
French RJ and Terlau H (2004) Sodium channel toxinsreceptor targeting and therapeutic potential. Curr Med Chem 11: 30533064.[Medline]
Head C and Gardiner M (2003) Paroxysms of excitement: sodium channel dysfunction in heart and brain. Bioessays 25: 981993.[CrossRef][Medline]
Julius D and Basbaum AI (2001) Molecular mechanisms of nociception. Nature (Lond) 413: 203210.[CrossRef][Medline]
Kayano T, Noda M, Flockerzi V, Takahashi H, and Numa S (1988) Primary structure of rat brain sodium channel III deduced from the cDNA sequence. FEBS Lett 228: 187194.[CrossRef][Medline]
Lee CW, Kim S, Roh SH, Endoh H, Kodera Y, Maeda T, Kohno T, Wang JM, Swartz KJ, and Kim JI (2004) Solution structure and functional characterization of SGTx1, a modifier of Kv2.1 channel gating. Biochemistry 43: 890897.[Medline]
Lewis RJ and Garcia ML (2003) Therapeutic potential of venom peptides. Nat Rev Drug Discov 2: 790802.[CrossRef][Medline]
Li D, Xiao Y, Hu W, Xie J, Bosmans F, Tytgat J, and Liang S (2003) Function and solution structure of hainantoxin-I, a novel insect sodium channel inhibitor from the Chinese bird spider Selenocosmia hainana. FEBS Lett 555: 616622.[CrossRef][Medline]
Li D, Xiao Y, Xu X, Xiong X, Lu S, Liu Z, Zhu Q, Wang M, Gu X, and Liang S (2004) Structure-activity relationships of hainantoxin-IV and structure determination of active and inactive sodium channel blockers. J Biol Chem 279: 3773437740.
Li-Smerin Y and Swartz KJ (1998) Gating modifier toxins reveal a conserved structural motif in voltage-gated Ca2+ and K+ channels. Proc Natl Acad Sci USA 95: 85858589.
McDonough SI, Lampe RA, Keith RA, and Bean BP (1997) Voltage-dependent inhibition of N- and P-type calcium channels by the peptide toxin
-grammotoxin-SIA. Mol Pharmacol 52: 10951104.
Middleton RE, Warren VA, Kraus RL, Hwang JC, Liu CJ, Dai G, Brochu RM, Kohler MG, Gao YD, Garsky VM, et al. (2002) Two tarantula peptides inhibit activation of multiple sodium channels. Biochemistry 41: 1473414747.[CrossRef][Medline]
Nassar MA, Stirling LC, Forlani G, Baker MD, Matthews EA, Dickenson AH, and Wood JN (2004) Nociceptor-specific gene deletion reveals a major role for Nav1.7 (PN1) in acute and inflammatory pain. Proc Natl Acad Sci USA 101: 1270612711.
Newcomb R, Szoke B, Palma A, Wang G, Chen X, Hopkins W, Cong R, Miller J, Urge L, Tarczy-Hornoch K, et al. (1998) Selective peptide antagonist of the class E calcium channel from the venom of the tarantula Hysterocrates gigas. Biochemistry 37: 1535315362.[CrossRef][Medline]
Peng K, Shu Q, Liu Z, and Liang S (2002) Function and solution structure of huwentoxin-IV, a potent neuronal tetrodotoxin (TTX)-sensitive sodium channel antagonist from Chinese bird spider Selenocosmia huwena. J Biol Chem 277: 4756447571.
Possani LD, Becerril B, Delepierre M, and Tytgat J (1999) Scorpion toxins specific for Na+-channels. Eur J Biochem 264: 287300.[Medline]
Rash LD and Hodgson WC (2002) Pharmacology and biochemistry of spider venoms. Toxicon 40: 225254.[Medline]
Smith JJ, Alphy S, Seibert AL, and Blumenthal KM (2005) Differential phospholipid binding by site 3 and site 4 toxins: Implications for structural variability between voltage-sensitive sodium channel domains. J Biol Chem. 25: 1112711133.
Smith RD and Goldin AL (1998) Functional analysis of the rat I sodium channel in Xenopus oocytes. J Neurosci 18: 811820.
Sollod BL, Wilson D, Zhaxybayeva O, Gogarten JP, Drinkwater R, and King GF (2005) Were arachnids the first to use combinatorial peptide libraries? Peptides 26: 131139.[CrossRef][Medline]
Suchyna TM, Johnson JH, Hamer K, Leykam JF, Gage DA, Clemo HF, Baumgarten CM, and Sachs F (2000) Identification of a peptide toxin from Grammostola spatulata spider venom that blocks cation-selective stretch-activated channels. J Gen Physiol 115: 583598.
Swanson R, Marshall J, Smith JS, Williams JB, Boyle MB, Folander K, Luneau CJ, Antanavage J, Oliva C, Buhrow SA and et al. (1990) Cloning and expression of cDNA and genomic clones encoding three delayed rectifier potassium channels in rat brain. Neuron 4: 929939.[CrossRef][Medline]
Swartz KJ and MacKinnon R (1995) An inhibitor of the Kv2.1 potassium channel isolated from the venom of a Chilean tarantula. Neuron 15: 941949.[CrossRef][Medline]
Takahashi H, Kim JI, Min HJ, Sato K, Swartz KJ, and Shimada I (2000) Solution structure of hanatoxin1, a gating modifier of voltage-dependent K+ channels: common surface features of gating modifier toxins. J Mol Biol 297: 771780.[CrossRef][Medline]
Wang JM, Roh SH, Kim S, Lee CW, Kim JI, and Swartz KJ (2004) Molecular surface of tarantula toxins interacting with voltage sensors in K(v) channels. J Gen Physiol 123: 455467.
Wang SY and Wang GK (2003) Voltage-gated sodium channels as primary targets of diverse lipid-soluble neurotoxins. Cell Signal 15: 151159.[CrossRef][Medline]
Winterfield JR and Swartz KJ (2000) A hot spot for the interaction of gating modifier toxins with voltage-dependent ion channels. J Gen Physiol 116: 637644.
Wood JN, Boorman JP, Okuse K, and Baker MD (2004) Voltage-gated sodium channels and pain pathways. J Neurobiol 61: 5571.[CrossRef][Medline]
Xiao Y and Liang S (2003) Inhibition of neuronal tetrodotoxin-sensitive Na+ channels by two spider toxins: hainantoxin-III and hainantoxin-IV. Eur J Pharmacol 477: 17.[CrossRef][Medline]
Xiao Y, Tang J, Yang Y, Wang M, Hu W, Xie J, Zeng X, and Liang S (2004) Jingzhaotoxin-III, a novel spider toxin inhibiting activation of voltage-gated sodium channel in rat cardiac myocytes. J Biol Chem 279: 2622026226.
Yu FH and Catterall WA (2003) Overview of the voltage-gated sodium channel family. Genome Biol 4: 207.[CrossRef]