Impact Factor 4.225 | CiteScore 5.0
More on impact ›

Brief Research Report ARTICLE

Front. Pharmacol., 07 July 2020 | https://doi.org/10.3389/fphar.2020.01010

Tuning Scorpion Toxin Selectivity: Switching From KV1.1 to KV1.3

  • 1Shemyakin-Ovchinnikov Institute of Bioorganic Chemistry, Russian Academy of Sciences, Moscow, Russia
  • 2Toxicology and Pharmacology, KU Leuven, Leuven, Belgium
  • 3School of Biomedicine, Far Eastern Federal University, Vladivostok, Russia
  • 4Department of Applied Mathematics, National Research University Higher School of Economics, Moscow, Russia
  • 5Moscow Institute of Physics and Technology (State University), Dolgoprudny, Russia

Voltage-gated potassium channels (KVs) perform vital physiological functions and are targets in different disorders ranging from ataxia and arrhythmia to autoimmune diseases. An important issue is the search for and production of selective ligands of these channels. Peptide toxins found in scorpion venom named KTx excel in both potency and selectivity with respect to some potassium channel isoforms, which may present only minute differences in their structure. Despite several decades of research the molecular determinants of KTx selectivity are still poorly understood. Here we analyze MeKTx13-3 (Kalium ID: α-KTx 3.19) from the lesser Asian scorpion Mesobuthus eupeus, a high-affinity KV1.1 blocker (IC50 ~2 nM); it also affects KV1.2 (IC50 ~100 nM), 1.3 (~10 nM) and 1.6 (~60 nM). By constructing computer models of its complex with KV1.1–1.3 channels we identify specific contacts between the toxin and the three isoforms. We then perform mutagenesis to disturb the identified contacts with KV1.1 and 1.2 and produce recombinant MeKTx13-3_AAAR, which differs by four amino acid residues from the parent toxin. As predicted by the modeling, this derivative shows decreased activity on KV1.1 (IC50 ~550 nM) and 1.2 (~200 nM). It also has diminished activity on KV1.6 (~1500 nM) but preserves KV1.3 affinity as measured using the voltage-clamp technique on mammalian channels expressed in Xenopus oocytes. In effect, we convert a selective KV1.1 ligand into a new specific KV1.3 ligand. MeKTx13-3 and its derivatives are attractive tools to study the structure-function relationship in potassium channel blockers.

Introduction

It is believed that potassium (K+) channels arose near the time of life origin on the earth. K+ channels are key membrane proteins of all living organisms, and about 80 genes encoding the main α-subunits are found in mammalian genomes (Alexander et al., 2019). The most prevalent family of K+ channels in mammals is voltage-gated potassium channels (KVs) that includes 40 isoforms (Attali et al., 2019). These proteins control neuronal excitability, heart rate, muscle contraction, hormonal secretion, cell proliferation, etc. It is not surprising that modulation of KVs provokes changes in the physiology of a cell or even of the whole organism (Hille, 2001).

KV1.3 is one of the most studied and pharmacologically important isoforms of K+ channels. At least two major directions of biomedical research are associated with this type of channel. First, KV1.3 in T lymphocytes is a validated target for diverse autoimmune diseases, such as multiple sclerosis, rheumatoid arthritis, and type 1 diabetes (Chandy et al., 2004; Beeton et al., 2006; Feske et al., 2012). Second, this protein is a crucial participant in a number of cancers because it is necessary for cell proliferation, malignant angiogenesis, and metastasis (Pardo and Stühmer, 2014; Chandy and Norton, 2016; Teisseyre et al., 2019). For both of these directions selective and effective inhibitors of KV1.3 are desirable. Novel ligands and their derivatives are considered as promising molecular instruments in KV1.3 research and are exploited as templates in drug design (Wulff and Zhorov, 2008; Chandy and Norton, 2017; Prosdocimi et al., 2019). Active compounds affecting KV1.3 can be obtained from different natural sources, such as plant extracts and animal venoms (King, 2011; Norton and Chandy, 2017), as well as synthesized de novo (Schmitz et al., 2005; Hendrickx et al., 2020).

Scorpion venom serves as an abundant source of toxins acting as K+ channel ligands (KTx), which have evolved and been selected for a highly efficient interaction with their molecular targets, including KVs (Kuzmenkov et al., 2015a). According to Kalium database (https://kaliumdb.org/), these compounds are polypeptides containing 23 to 78 amino acid residues and cross-linked by two to four intramolecular disulfide bonds (Kuzmenkov et al., 2016a; Tabakmakher et al., 2019). A dominating number of known KTx adopts the CSα/β (cysteine-stabilized α-helix and β-sheet) fold, but some of them present other types of fold (Mouhat et al., 2004; Kuzmenkov et al., 2015a). A number of structural and pharmacological findings pinpointed important determinants in the interfaces of KVs and KTx contact that contribute to toxin selectivity (Aiyar et al., 1995; Hidalgo and MacKinnon, 1995; Gross and MacKinnon, 1996; Giangiacomo et al., 2004). Perhaps most importantly, obtaining the crystal structure of the KV1.2/KV2.1 paddle chimera in complex with charybdotoxin (ChTx) (Banerjee et al., 2013) highlighted the key amino acids involved in the interaction and opened new opportunities for scaffold engineering of more selective KTx (Han et al., 2008; Kuzmenkov et al., 2018).

Earlier, we identified and purified MeKTx13-3 toxin (Kalium ID: a-KTx 3.19, UniProt ID: C0HJQ6, 37 residues, three disulfide bonds) from the venom of the lesser Asian scorpion Mesobuthus eupeus (Kuzmenkov et al., 2015b). We performed pharmacological profiling of this KTx on several isoforms of KVs and found that it is active on KV1.1–1.3 and 1.6 with half-maximal inhibitory concentration (IC50) values of ~2, 100, 10, and 60 nM, respectively. The toxin preferably blocked KV1.1, however, cross-reactivity with KV1.3 was also observed (Kuzmenkov et al., 2019). Since a large number of KTx inhibits both KV1.1 and 1.3 in a similar manner (Mouhat et al., 2005; Takacs et al., 2009; Gao et al., 2010), the goal of our present work is to identify molecular determinants responsible for the interaction with different channel isoforms and switch the selectivity of MeKTx13-3 from KV1.1 to KV1.3.

Materials and Methods

Ethics Statement

This study strictly complied with the World Health Organization’s International Guiding Principles for Biomedical Research Involving Animals. The research was carried out in AAALAC accredited organization according to the standards of the Guide for Care and Use of Laboratory Animals (8th edition, Institute for Laboratory Research of Animals). All experiments were approved by the Institutional Policy on the Use of Laboratory Animals of the Shemyakin-Ovchinnikov Institute of Bioorganic Chemistry Russian Academy of Sciences (Protocol Number 267/2018; date of approval: 28 February 2019).

Homology Modeling of Toxins and Their Complexes With KVs

Since the amino acid sequence of MeKTx13-3 is identical to that of BmKTX (Romi-Lebrun et al., 1997), the known 3D structure of the latter (PDB ID: 1BKT) (Renisio et al., 2000) was used in our work. KV1.1 model was generated in MODELLER 9.19 (Webb and Sali, 2016) using the KV1.2 structure (PDB ID: 3LUT) (Chen et al., 2010) as a template. KV1.3 model has been generated previously (Kuzmenkov et al., 2017; Kuzmenkov et al., 2018; Berkut et al., 2019) using an analogous procedure.

Complexes of MeKTx13-3 with KVs were modeled considering that the toxin interacts with the channels similarly to ChTx, one of the most thoroughly studied KTx (Goldstein et al., 1994). The model of the complex of MeKTx13-3 with KV1.2 was built on the basis of the KV1.2/2.1–ChTx complex crystal structure (Banerjee et al., 2013): the structure of MeKTx13-3 was spatially aligned with the structure of channel-bound ChTx, which was subsequently replaced by the aligned toxin. Complexes with KV1.1 and 1.3 were generated similarly, but the first step was spatial alignment of the channel models with the KV1.2/2.1 chimera (Kuzmenkov et al., 2017; Kuzmenkov et al., 2018; Berkut et al., 2019).

Molecular Dynamics Simulations

The resulting complexes of MeKTx13-3 with KVs were placed inside a lipid bilayer mimicking a neuronal membrane. We used a pre-equilibrated fragment of bilayer (7.0 × 7.0 × 13.5 nm3; 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine/1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine/cholesterol, POPC : POPE : Chl; 100:50:50 molecules, respectively, solvated with 14172 water molecules) that has been described in detail in our previous works (Berkut et al., 2019); some phospholipid and Chl molecules were removed to provide room for the protein. The TIP3P water model (Jorgensen et al., 1983) and the required number of Na+ ions (to maintain electroneutrality) were used for resolvation. All systems were equilibrated (heated up to 37°C) during 100 ps of molecular dynamics (MD) simulation. Positions of the channel Cα atoms of residues not involved in the channel pore vestibule, as well as the Nϵ atom of Lys26 in MeKTx13-3 were restrained during the equilibration to prevent destabilization of the initial complex. Systems were then subjected to 500 ns of MD. All simulations were performed with the GROMACS software (Abraham et al., 2015) (version 2018) using the AMBER99SB-ILDN parameters set (Klepeis et al., 2010). Simulations were carried out with a time step of 2 fs, imposing 3D periodic boundary conditions, in the isothermal-isobaric (NPT) ensemble with a semi-isotropic pressure of 1 bar using the Berendsen pressure coupling algorithm (Berendsen et al., 1984), and at a constant temperature of 37°C using the V-rescale thermostat (Bussi et al., 2007). Van der Waals interactions were truncated using a 1.5-nm spherical cut-off function. Electrostatic interactions were treated with the PME algorithm. During the simulation, the position of the Nϵ atom of Lys26 in each complex was restrained inside the channel pore.

Determination of Interaction Energy and Intermolecular Contacts

We determined the intermolecular contacts during MD and estimated residual contributions to intermolecular interaction energy based on MD trajectory using our in-house software package IMPULSE (Krylov et al., in preparation) analogously to the procedures described in our previous study (Berkut et al., 2019). Briefly, H-bonds were assigned using the parameters set from the hbond utility of GROMACS software (Abraham et al., 2015) (the distance D—A ≤ 0.35 nm and the angle D—H—A ≥ 150° for the hydrogen bond D—H···A, where D and A are the hydrogen bond donor and acceptor, respectively); salt bridges, cation-π, stacking, and hydrophobic contacts were calculated using algorithms described in our previous works (Pyrkov and Efremov, 2007; Pyrkov et al., 2009). The AMBER99SB-ILDN parameters set (Klepeis et al., 2010) and 1.5 nm cutoff distance for Lennard-Jones and electrostatic interactions were used during the intermolecular short-range non-bonded interaction energy estimation, the latter being the sum of the Lennard-Jones and electrostatic terms. All drawings of 3D structures were prepared with the PyMOL Molecular Graphics System, version 1.8 (Schrödinger, LLC). Graphical representation of interaction energy profiles was performed using Python built-in libraries and the NumPy package.

Toxin Isolation From Scorpion Venom

Natural MeKTx13-3 (α-KTx 1.19) was isolated from the same stock of M. eupeus venom and following the same procedure as described previously (Kuzmenkov et al., 2015b; Kuzmenkov et al., 2019).

Recombinant Peptide Production

Recombinant MeKTx13-3 and its derivative were produced using an approach elaborated previously (Pluzhnikov et al., 2007). Briefly, the peptides were produced in a bacterial expression system as fusions with the carrier protein thioredoxin (Trx) (McCoy and LaVallie, 1994) and recombinant human enteropeptidase light chain (Gasparian et al., 2011) was used to cleave the fusions.

DNA sequences encoding MeKTx13-3 and its derivative were constructed from synthetic oligonucleotides by PCR in two steps (see Supplemental Figure S1). On the first step the target PCR fragments were amplified in 5 cycles using two forward primers and two reverse primers (F1, F2, R1, and R2). The four primers altogether constitute a full gene sequence. For the second step, PCR mixtures from the first step were diluted 100 times, and 1 µl of the dilution was used as a matrix; only the terminal primers (F1 and R1) were used for the amplification (Supplemental Table S1). The resulting PCR fragments were cloned into the expression vector pET-32b (Novagen) using KpnI and BamHI restriction enzymes to produce the vectors coding for the target polypeptides.

Escherichia coli SHuffle T7 Express cells (New England Biolabs) were transformed using the corresponding expression vectors and cultured at 30°C in LB medium to the mid-log phase. Expression was then induced by 0.2 mM Isopropyl β-d-1-thiogalactopyranoside. Cells were cultured at room temperature (24°C) overnight (16 h) and harvested by centrifugation. The cell pellet was resuspended in 300 mM NaCl, 50 mM Tris-HCl buffer (pH 8.0) and ultrasonicated. The lysate was applied to a HisPur Cobalt Resin (ThermoFisher Scientific); and the Trx-fusion proteins were purified according to the manufacturer’s protocol.

Fusion proteins were dissolved in 50 mM Tris-HCl (pH 8.0) to a concentration of 1 mg/ml. Protein cleavage with human enteropeptidase light chain (1 U of enzyme per 1 mg of substrate) was performed overnight (16 h) at 37°C. Recombinant peptides were purified by reversed-phase HPLC on a Jupiter C5 column (4.6 × 250 mm; Phenomenex) in a linear gradient of acetonitrile concentration (0–60% in 60 min) in the presence of 0.1% trifluoroacetic acid. The purity of the target peptides was checked by MALDI MS and analytical chromatography on a Vydac C18 column (4.6 × 250 mm; Separations Group) in the same acetonitrile gradient.

Mass Spectrometry

Molecular mass measurements for natural and recombinant peptides were performed using MALDI on an Ultraflex TOF-TOF (Bruker Daltonik) spectrometer as described earlier (Kuzmenkov et al., 2016b). 2,5-Dihydroxybenzoic acid (Sigma-Aldrich) was used as a matrix. Measurements were carried out in both linear and reflector modes. Mass spectra were analyzed with the Data Analysis 4.3 and Data Analysis Viewer 4.3 software (Bruker).

Ion Channel Expression in Xenopus Oocytes

All procedures were performed in agreement with the guidelines of ARRIVE (Animal Research: Reporting of In Vivo Experiments) and the “European convention for the protection of vertebrate animals used for experimental and other scientific purposes” (Strasbourg, 18.III.1986).

The major pipeline of ion channel expression in oocytes was described previously (Peigneur et al., 2011). Briefly, for the expression of KV genes (rat (r)KV1.1, rKV1.2, human (h)KV1.3, rKV1.4, rKV1.5, and rKV1.6) in Xenopus laevis oocytes, linearized plasmids containing the respective gene sequences were transcribed using the T7 mMESSAGE-mMACHINE transcription kit (Ambion). 50 nl of cRNA solution (1 ng/nl) were injected into oocytes using a micro-injector (Drummond Scientific). The oocytes were incubated in ND96 solution: 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.

Electrophysiological Recordings

Two-electrode voltage-clamp recordings were performed at room temperature (18–22°C) using a Geneclamp 500 amplifier (Molecular Devices) controlled by a pClamp data acquisition system (Axon Instruments) as described (Peigneur et al., 2011). Bath solution composition was ND96. KV currents were evoked by 250-ms depolarization to 0 mV from a holding potential of −90 mV, followed by 250-ms pulses to −50 mV. For current–voltage relationship studies, currents were evoked by 10-mV depolarization steps. Concentration–response curves were constructed, in which the percentage of current inhibition was plotted as a function of toxin concentration. Data were fitted with the Hill equation: y = 100/[1 + (IC50/[toxin])h], where y is the amplitude of the toxin-induced effect, [toxin] is toxin concentration, and h is the Hill coefficient. Comparison of two sample means was performed using a paired Student’s t-test (p-value of 0.05 was used as a threshold of significance). All data were obtained in at least three independent experiments (n ≥ 3) and are presented as mean ± standard error of the mean.

Results

Computational Study

Amino acid sequence of MeKTx13-3 is identical to BmKTX that was isolated from Mesobuthus martensii, a close relative of M. eupeus (Romi-Lebrun et al., 1997). Since the 3D structure of BmKTX is known (PDB ID: 1BKT) (Renisio et al., 2000), we used it to generate models of MeKTx13-3 in complex with KV1.1–1.3. The models were then subjected to MD simulations (Figures 1A–C).

FIGURE 1
www.frontiersin.org

Figure 1 (A–C) Modeled structure of MeKTx13-3 in complex with KV1.1–1.3. (A) Overall structure of the KV1.3–MeKTx13-3 complex after 100-ns MD simulation inside a hydrated lipid bilayer membrane. Four channel α-subunits with identical sequences are color-coded. The pore domain helices of the channel subunit in the foreground and voltage-sensing domain (VSD) of the adjacent subunit, as well as extended extracellular loops of the VSDs are omitted for clarity. Lipids are shown in a semi-transparent space-filling representation; atoms are colored: oxygen, red; phosphorus, orange; nitrogen, blue; hydrogen of amino group, white; carbon of POPC, light-yellow; carbon of POPE, yellow; and carbon of cholesterol, beige. Some lipids are omitted for clarity. MeKTx13-3 is presented in pink; residue Lys26 (plugs the channel pore) is shown as sticks. (B, C) Close-up view on the channel pore vestibule area in complexes KV1.1–MeKTx13-3 and KV1.2–MeKTx13-3, respectively. Channels are shown in a semi-transparent representation. Lys26 and residues involved in the intermolecular contacts not present in the KV1.3–MeKTx13-3 complex are shown as sticks. Hydrogen bonds and salt bridges are shown as dashed yellow lines. Lipids are omitted for clarity. (D) Amino acid sequence alignment of the extracellular pore region of KV1.1–1.3 channels. Residue numbering is above each sequence; different residues are shaded gray; sequences of S5-P loops containing channel-specific residues are underlined.

To shed light on the molecular differences in MeKTx13-3 interaction with KV isoforms, we analyzed intermolecular contacts and residual contributions to interaction energy during the MD simulations using our in-house software package IMPULSE (Krylov et al., in preparation). It was observed that in complex with KV1.2 the toxin does not form any cation-π or stacking contacts, unlike in complex with KV1.1 or 1.3 (see Supplemental Table S2). This observation is in good agreement with electrophysiological measurements that displayed the preferred activity of MeKTx13-3 against KV1.1 and 1.3 (Kuzmenkov et al., 2019). We noted the following specific contacts between the toxin and particular channel isoforms.

(1) MeKTx13-3 residue Lys15 forms an H-bond and a salt bridge with KV1.1 residue Glu353 (Figure 1B). Analogous contact with Asp351 (Figure 1D) is not observed in complex with KV1.2, apparently due to the fact that the neighboring residue Arg350 repulses Lys15 and prevents this contact formation. No contact with Thr375 was found in complex with KV1.3 either, presumably due to the small size of the threonine side chain and the lack of electrostatic attraction to reach Lys15.

(2) MeKTx13-3 residue Lys18 forms an H-bond with KV1.1 residue Ser354 and a cation-π contact with His355; in complex with KV1.2 it forms an H-bond with Gln353 (Figures 1B, C). In the complex of MeKTx13-3 with KV1.3 the conformation of the loop containing channel-specific residues Thr375 and Gly377 (Figure 1D) during MD is such that it does not reach Lys18, so no specific contacts are observed.

(3) MeKTx13-3 residue Gln12 forms an H-bond with KV1.2 residue Asp351 (Figure 1C). Analogous contact with Thr375 (Figure 1D) is not observed in complex with KV1.3 because the short threonine side chain does not reach Gln12. No contact with Glu353 is found in complex with KV1.1 because the conformation of the loop containing channel-specific residues during MD prevents reaching Gln12.

In addition to the observed channel-specific contacts, MeKTx13-3 residue Asp33 makes a significant positive contribution to the binding energy (negatively affects the affinity) (see Supplemental Figure S2) due to electrostatic repulsion with a conserved negatively charged residue in the channel vestibule (Asp377/375/399 in KV1.1/1.2/1.3, see Figure 1D).

Pharmacological profiling and detailed complex structure analysis allow us to propose several point substitutions in MeKTx13-3 for switching toxin selectivity. Since Gln12, Lys15, and Lys18 form hydrogen bonds, salt bridges, and cation-π interactions in complexes with KV1.1 and 1.2, the general idea of the modifications was to abolish the formation of these polar contacts. Moreover, since these toxin residues do not form specific contacts with KV1.3, it is reasonable to assume that such modifications will not affect the affinity to this channel isoform significantly. Therefore, we suggested a derivative of MeKTx13-3 in which Gln12, Lys15, and Lys18 are replaced by Ala to reduce toxin affinity to KV1.1 and 1.2, and Asp33 is replaced by Arg to increase its affinity to KVs (Figure 2A).

FIGURE 2
www.frontiersin.org

Figure 2 Production of MeKTx13-3 and its derivative. (A) Amino acid sequence alignment of MeKTx13-3 and MeKTx13-3_AAAR. Gray shading indicates the positions where replacements were introduced. Cysteine residues are in bold, and lines above the sequences indicate disulfide bonds. Take a note that recombinant analogue of MeKTx13-3 does not bear the C-terminal amidation of the natural toxin. (B, C) Reversed-phase HPLC separation of recombinant MeKTx13-3 and MeKTx13-3_AAAR after digestion by enteropeptidase. For target peptides the measured molecular masses are indicated.

Recombinant Toxin Production

The natural toxin was purified from crude venom as we described previously (Kuzmenkov et al., 2015b; Kuzmenkov et al., 2019). Recombinant MeKTx13-3 and its derivative MeKTx13-3_AAAR (MeKTx13-3 with the following replacements: Gln12Ala, Lys15Ala, Lys18Ala, and Asp33Arg) were obtained according to our common protocol (Pluzhnikov et al., 2007) using E. coli SHuffle B strain as an expression system. Synthetic genes encoding the peptides were cloned into the pET-32b expression vector, and Trx was used as a fusion partner to ensure a high yield of the disulfide-containing peptides with native conformation. The target peptides were produced as a result of fusion protein cleavage by enteropeptidase followed by separation using reversed-phase HPLC and identification by MALDI MS (Figures 2B, C). The measured molecular masses of the purified peptides were equal to the calculated values. The final yield of the peptides was ~5 mg per 1 l of bacterial culture. Note that the natural toxin is C-terminally amidated (Kuzmenkov et al., 2015b), whereas this modification is missing from the recombinant peptides.

Electrophysiology

We first compared the activity of the natural and recombinant MeKTx13-3 (Supplemental Figure S3 and Supplemental Table S3) on KV1.1. Recombinant peptide was less active than the native toxin (IC50 values of 6.7 ± 2.7 and 1.9 ± 0.2 nM, respectively), which is probably due to the lack of C-terminal amidation in the former.

We then estimated the activities of the obtained recombinant peptides MeKTx13-3 and MeKTx13-3_AAAR at a concentration of 10 nM against six isoforms of KVs (KV1.1–1.6) expressed in X. laevis oocytes (Figure 3A). Recombinant MeKTx13-3 as well as the natural toxin inhibited almost completely (by >95%) potassium currents through KV1.1; ~25%, ~50%, and ~15% of current through KV1.2, 1.3, and 1.6 was blocked. At the same concentration of 10 nM MeKTx13-3_AAAR blocked KV1.3 by ~50%, whereas only ~15%, ~5% and ~10% was blocked in KV1.1, 1.2 and 1.6. Neither KV1.4 nor KV1.5 were affected by any of the peptides.

FIGURE 3
www.frontiersin.org

Figure 3 Electrophysiological profiling of MeKTx13-3 and MeKTx13-3_AAAR activities. (A) Representative traces of currents through KV1.1–1.6 in control (black) and after application of 10 nM toxin (blue). (B) Concentration–response curves of MeKTx13-3 (left) and MeKTx13-3_AAAR (right) on KV1.1–1.3 and 1.6 obtained by electrophysiological measurements. IC50 values are listed in Table 1.

Finally, we constructed dose–response curves for the susceptible channels (Figure 3B; see Table 1 for IC50 values and Hill coefficients). MeKTx13-3_AAAR demonstrated a comparable activity with native MeKTx13-3 on KV1.3 (IC50 = 8.9 ± 0.9 nM for the natural toxin and 9.1 ± 0.4 nM for the mutant), whereas its affinity to KV1.1 decreased dramatically (IC50 = 541.5 ± 48.6 nM instead of 1.9 ± 0.2 nM for natural MeKTx13-3). MeKTx13-3_AAAR also showed reduced activity on KV1.2 (IC50 = 208.2 ± 15.2 nM compared to 105.9 ± 14.6 nM for MeKTx13-3) and KV1.6 (IC50 = 1522.3 ± 183.4 nM instead of 63.4 ± 4.5 nM).

TABLE 1
www.frontiersin.org

Table 1 IC50 Values (in nM) Calculated for MeKTx13-3 and Its Derivative against KV1.1–1.6.

Discussion

Animal venom serves a rich source of promising compounds affecting ion channels, which can be utilized as potential drug hits (Wulff et al., 2019). Detailed studies based on either toxin or channel mutagenesis are essential for (i) the understanding of fine molecular contacts between the toxins and very close channel isoforms, and (ii) design and production of more selective ligands. Hence prediction of critical amino acids involved in toxin–channel complex formation in silico is a convenient and powerful approach for following mutagenesis studies (Kuyucak and Norton, 2014). For instance, such computationally guided assay helped to design highly selective peptide drug hits or leads, such as ShK-192 (Pennington et al., 2009) and HsTX1[R14A] (Rashid et al., 2015).

Here, we have designed and produced a derivative of scorpion toxin MeKTx13-3 with its selectivity switched from KV1.1 to KV1.3. MeKTx13-3 is one of a limited number of known animal toxins that possess selectivity to homotetrameric KV1.1 (Kuzmenkov et al., 2019). We introduced several substitutions according to computer modeling experiments. To convert the selectivity of MeKTx13-3 we replaced three amino acids by alanine (Gln12Ala, Lys15Ala, and Lys18Ala) to abolish the formation of H-bonds, salt bridges, or cation-π interactions in the complexes with KV1.1 and 1.2. Moreover, to prevent the electrostatic repulsion between the negatively charged Asp33 in MeKTx13-3 and the conserved aspartic acid residue in the P-S6 loop of KVs we introduced the Asp33Arg replacement.

We produced not only MeKTx13-3_AAAR, but also recombinant MeKTx13-3 to confirm that (i) the peptide folding is correct, and (ii) lack of the C-terminal amidation does not impact toxin activity dramatically. These points are of importance, because misfolded toxins usually lose activity, whereas C-terminal amidation can boost ligand potency (Lebrun et al., 1997; Coelho et al., 2014). Voltage-clamp recordings in X. laevis oocytes showed that recombinant MeKTx13-3 is ~3.5 times less potent than the natural toxin. We attribute this decrease in activity to the amidation. It is well known that this post-translational modification can affect the activity of peptides. The effects may vary from dramatic to negligible, with most apparent cases found in hormones (Merkler, 1994). As for potassium channel blockers, the C-terminal amidation of ShK from the sea anemone Stichodactyla helianthus resulted in ~4 times decreased potency against KV1.3 (Pennington et al., 2012). Conversely, the amidation of HsTX1 (α-KTx 6.3) from the scorpion Heterometrus spinifer increases the activity against KV1.3 five-fold (Lebrun et al., 1997). The effects in our case are subtle, and we decided to use the recombinantly produced peptide MeKTx13-3 lacking the amidation in further studies.

The designed derivative MeKTx13-3_AAAR was also tested on six KV isoforms and as we expected the selectivity of this analogue shifted towards KV1.3. A graphic approach to estimate toxin specificity to a particular channel isoform (KV1.3 in our case) is using the selectivity factor, i.e. the ratio of IC50 (or Kd) values for two channels. MeKTx13-3 displays KV1.1/KV1.3 selectivity factor of approximately 0.2, while for MeKTx13-3_AAAR this parameter has changed to 60 (Supplemental Figure S4). There is a number of more KV1.3-specific toxins and their derivatives, for instance, HsTX1, Vm24, or moka1, demonstrating selectivity factors (KV1.1/KV1.3) of 500 or even 1000 (Romi-Lebrun et al., 1997; Takacs et al., 2009; Varga et al., 2012).

We aligned the sequences of several dozen toxins from KTx subfamilies α-KTx1–4, 11, 12, 15–18, and 21–24 with comparable length and cysteine pattern to MeKTx13-3 and found that two of the residues mutated in our study (Gln12, the first amino acid before the second Cys residue; and Lys18, the first amino acid after the third Cys) are quite conserved. In numerous toxins Gln/Glu and Lys/Arg are located in these positions, respectively. On the other hand, Lys15 between second and third Cys and Asp33 between the fifth and sixth Cys are variable. Within α-KTx 3 subfamily (to which MeKTx13-3 belongs) these positions display a similar pattern (see Supplemental Figure S4). All four residues seem to contribute to the bioactive surface of the toxins, and at present, there is no apparent explanation as to why two of them are more conserved than the other two. Moreover, we cannot infer any correlation between these residues and toxin potency or selectivity.

One may argue that the same substitutions as we introduced to MeKTx13-3 might bring about similar changes of selectivity in other α-KTx 3 toxins. However, and quite unfortunately, our current understanding of KV–KTx interactions does not allow to predict the specificity of toxins from primary structure. It appears that in each case molecular modeling experiments and a detailed analysis of the contacts are required. This is because one change in the sequence may actually affect how other residues interact—due to sterical hindrances or packing effects, electrostatic attraction or repulsion, H-bond reshuffling, or local folding rearrangements, all of which are not easily discernible from primary structures. For instance, in case of MeKTx13-3 Lys15 seems to make a salt bridge with Glu353 in KV1.1 (Figure 1B). Simple consideration of the primary structures would predict the same salt bridge in KV1.2 since it has Asp351 in the same position of the alignment (Figure 1D). This contact is not established in molecular modeling however, due to the neighboring Arg350, which repulses Lys15. Similarly, MeKTx13-3 residues Lys18 and Gln12 make channel-specific contacts due to the different folding of the channel extracellular loops, not just the amino acid substitutions. In conclusion, we hope to have demonstrated here the possibility of switching toxin specificity between two very close channel isoforms based on careful in silico design.

Data Availability Statement

The datasets generated for this study can be found in the article/Supplementary Material.

Ethics Statement

This study strictly complied with the World Health Organization’s International Guiding Principles for Biomedical Research Involving Animals. The research was carried out in AAALAC accredited organization according to the standards of the Guidefor Care and Use of Laboratory Animals (8th edition, Institute forLaboratory Research of Animals). The use of the frogs was in accordance with the license number LA1210239 of the Laboratory of Toxicology & Pharmacology, University of Leuven. The use of Xenopus laevis was approved by the Ethical Committee for animal experiments of the University of Leuven (P186/2019). All animal care and experimental procedures agreed with the guidelines of ‘European convention for the protection of vertebrate animals used for experimental and other scientific purposes’ (Strasbourg, 18.III.1986).

Author Contributions

AK and AV designed research. AG, AK, EP-J, SP, VT, and AV analyzed data. AJK, EP-J, SP, VT, AG, and AC performed research. AK, VT, and AV wrote the paper. RE and JT supervised molecular modeling and electrophysiology, respectively.

Funding

This work was supported by the Russian Science Foundation (grant no. 18-74-00125). JT was funded by grants G0E7120N, GOC2319N and GOA4919N from the F.W.O Vlaanderen. SP was supported by KU Leuven funding (PDM/19/164). Fundação de Amparo à Pesquisa do Estado de São Paulo (São Paulo Research Foundation, Brazil, scholarship to EP-J n. 2016/04761-4), and Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (Coordination for the Improvement of Higher Education Personnel, Brazil; scholarship to EP-J n. 88881.186830/2018-01).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

The molecular dynamics simulations were carried out using the computational facilities of the Supercomputer Center “Polytechnical” at the St. Petersburg Polytechnic University and IACP FEB RAS Shared Resource Center “Far Eastern Computing Resource” equipment (https://cc.dvo.ru). Supercomputer calculations were sponsored in the framework of the Basic Research Program at the National Research University Higher School of Economics and Russian Academic Excellence Project “5-100”. We thank O. Pongs for sharing the rKv1.2, rKv1.4, rKv1.5, and rKv1.6 cDNA and we are grateful to M.L. Garcia for sharing the hKv1.3 clone.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphar.2020.01010/full#supplementary-material

References

Abraham, M. J., Murtola, T., Schulz, R., Páll, S., Smith, J. C., Hess, B., et al. (2015). GROMACS: High performance molecular simulations through multi-level parallelism from laptops to supercomputers. SoftwareX 1, 19–25. doi: 10.1016/j.softx.2015.06.001

CrossRef Full Text | Google Scholar

Aiyar, J., Withka, J. M., Rizzi, J. P., Singleton, D. H., Andrews, G. C., Lin, W., et al. (1995). Topology of the pore-region of a K+ channel revealed by the NMR-derived structures of scorpion toxins. Neuron 15, 1169–1181. doi: 10.1016/0896-6273(95)90104-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Alexander, S. P. H., Mathie, A., Peters, J. A., Veale, E. L., Striessnig, J., Kelly, E., et al. (2019). THE CONCISE GUIDE TO PHARMACOLOGY 2019/20: Ion channels. Br. J. Pharmacol. 176, S142–S228. doi: 10.1111/bph.14749

PubMed Abstract | CrossRef Full Text | Google Scholar

Attali, B., Chandy, K. G., Giese, M. H., Grissmer, S., Gutman, G. A., Jan, L. Y., et al. (2019). “Voltage-gated potassium channels (version 2019.4) in the IUPHAR/BPS Guide to Pharmacology Database,” Available at: https://wwws.ed.ac.uk/gtopdb-cite/article/view/3233. doi: 10.2218/gtopdb/F81/2019.4

CrossRef Full Text | Google Scholar

Banerjee, A., Lee, A., Campbell, E., Mackinnon, R. (2013). Structure of a pore-blocking toxin in complex with a eukaryotic voltage-dependent K(+) channel. Elife 2, e00594. doi: 10.7554/eLife.00594

PubMed Abstract | CrossRef Full Text | Google Scholar

Beeton, C., Wulff, H., Standifer, N. E., Azam, P., Mullen, K. M., Pennington, M. W., et al. (2006). Kv1.3 channels are a therapeutic target for T cell-mediated autoimmune diseases. Proc. Natl. Acad. Sci. U. S. A. 103, 17414–17419. doi: 10.1073/pnas.0605136103

PubMed Abstract | CrossRef Full Text | Google Scholar

Berendsen, H. J. C., Postma, J. P. M., van Gunsteren, W. F., DiNola, A., Haak, J. R. (1984). Molecular dynamics with coupling to an external bath. J. Chem. Phys. 81, 3684–3690. doi: 10.1063/1.448118

CrossRef Full Text | Google Scholar

Berkut, A. A., Chugunov, A. O., Mineev, K. S., Peigneur, S., Tabakmakher, V. M., Krylov, N. A., et al. (2019). Protein surface topography as a tool to enhance the selective activity of a potassium channel blocker. J. Biol. Chem. 294, 18349–18359. doi: 10.1074/jbc.RA119.010494

PubMed Abstract | CrossRef Full Text | Google Scholar

Bussi, G., Donadio, D., Parrinello, M. (2007). Canonical sampling through velocity rescaling. J. Chem. Phys. 126, 014101. doi: 10.1063/1.2408420

PubMed Abstract | CrossRef Full Text | Google Scholar

Chandy, K. G., Norton, R. S. (2016). Channelling potassium to fight cancer. Nature 537, 497–499. doi: 10.1038/nature19467

PubMed Abstract | CrossRef Full Text | Google Scholar

Chandy, K. G., Norton, R. S. (2017). Peptide blockers of Kv1.3 channels in T cells as therapeutics for autoimmune disease. Curr. Opin. Chem. Biol. 38, 97–107. doi: 10.1016/j.cbpa.2017.02.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Chandy, G., Wulff, H., Beeton, C., Pennington, M., Gutman, G. A., Cahalan, M. D. (2004). K+ channels as targets for specific immunomodulation. Trends Pharmacol. Sci. 25, 280–289. doi: 10.1016/j.tips.2004.03.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, X., Wang, Q., Ni, F., Ma, J. (2010). Structure of the full-length Shaker potassium channel Kv1.2 by normal-mode-based X-ray crystallographic refinement. Proc. Natl. Acad. Sci. U. S. A. 107, 11352–11357. doi: 10.1073/pnas.1000142107

PubMed Abstract | CrossRef Full Text | Google Scholar

Coelho, V. A., Cremonez, C. M., Anjolette, F. A. P., Aguiar, J. F., Varanda, W. A., Arantes, E. C. (2014). Functional and structural study comparing the C-terminal amidated β-neurotoxin Ts1 with its isoform Ts1-G isolated from Tityus serrulatus venom. Toxicon 83, 15–21. doi: 10.1016/j.toxicon.2014.02.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Feske, S., Skolnik, E. Y., Prakriya, M. (2012). Ion channels and transporters in lymphocyte function and immunity. Nat. Rev. Immunol. 12, 532–547. doi: 10.1038/nri3233

PubMed Abstract | CrossRef Full Text | Google Scholar

Gao, B., Peigneur, S., Tytgat, J., Zhu, S. (2010). A potent potassium channel blocker from Mesobuthus eupeus scorpion venom. Biochimie 92, 1847–1853. doi: 10.1016/j.biochi.2010.08.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Gasparian, M. E., Bychkov, M. L., Dolgikh, D. A., Kirpichnikov, M. P. (2011). Strategy for improvement of enteropeptidase efficiency in tag removal processes. Protein Expr. Purif. 79, 191–196. doi: 10.1016/j.pep.2011.04.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Giangiacomo, K. M., Ceralde, Y., Mullmann, T. J. (2004). Molecular basis of α-KTx specificity. Toxicon 43, 877–886. doi: 10.1016/j.toxicon.2003.11.029

PubMed Abstract | CrossRef Full Text | Google Scholar

Goldstein, S. A., Pheasant, D. J., Miller, C. (1994). The charybdotoxin receptor of a Shaker K+ channel: peptide and channel residues mediating molecular recognition. Neuron 12, 1377–1388. doi: 10.1016/0896-6273(94)90452-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Gross, A., MacKinnon, R. (1996). Agitoxin Footprinting the Shaker Potassium Channel Pore. Neuron 16, 399–406. doi: 10.1016/S0896-6273(00)80057-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Han, S., Yi, H., Yin, S.-J., Chen, Z.-Y., Liu, H., Cao, Z.-J., et al. (2008). Structural basis of a potent peptide inhibitor designed for Kv1.3 channel, a therapeutic target of autoimmune disease. J. Biol. Chem. 283, 19058–19065. doi: 10.1074/jbc.M802054200

PubMed Abstract | CrossRef Full Text | Google Scholar

Hendrickx, L. A., Dobričić, V., Toplak, Ž., Peigneur, S., Mašič, L. P., Tomašič, T., et al. (2020). Design and characterization of a novel structural class of Kv1.3 inhibitors. Bioorg. Chem. 98, 103746. doi: 10.1016/j.bioorg.2020.103746

PubMed Abstract | CrossRef Full Text | Google Scholar

Hidalgo, P., MacKinnon, R. (1995). Revealing the architecture of a K+ channel pore through mutant cycles with a peptide inhibitor. Science 268, 307–310. doi: 10.1126/science.7716527

PubMed Abstract | CrossRef Full Text | Google Scholar

Hille, B. (2001). Ion Channels of Excitable Membranes. Sinauer Associates, Inc.

Google Scholar

Jorgensen, W. L., Chandrasekhar, J., Madura, J. D., Impey, R. W., Klein, M. L. (1983). Comparison of simple potential functions for simulating liquid water. J. Chem. Phys. 79, 926–935. doi: 10.1063/1.445869

CrossRef Full Text | Google Scholar

King, G. F. (2011). Venoms as a platform for human drugs: translating toxins into therapeutics. Expert Opin. Biol. Ther. 11, 1469–1484. doi: 10.1517/14712598.2011.621940

PubMed Abstract | CrossRef Full Text | Google Scholar

Klepeis, J. L., Lindorff-Larsen, K., Shaw, D. E., Palmo, K., Dror, R. O., Maragakis, P., et al. (2010). Improved side-chain torsion potentials for the Amber ff99SB protein force field. Proteins Struct. Funct. Bioinforma. 78, 1950–1958. doi: 10.1002/prot.22711

CrossRef Full Text | Google Scholar

Kuyucak, S., Norton, R. S. (2014). Computational approaches for designing potent and selective analogs of peptide toxins as novel therapeutics. Future Med. Chem. 6, 1645–1658. doi: 10.4155/FMC.14.98

PubMed Abstract | CrossRef Full Text | Google Scholar

Kuzmenkov, A., II, Grishin, E. V., Vassilevski, A. A. (2015a). Diversity of Potassium Channel Ligands: Focus on Scorpion Toxins. Biochem (Mosc.) 80, 1764–1799. doi: 10.1134/S0006297915130118

CrossRef Full Text | Google Scholar

Kuzmenkov, A., II, Vassilevski, A. A., Kudryashova, K. S., Nekrasova, O. V., Peigneur, S., Tytgat, J., et al. (2015b). Variability of potassium channel blockers in Mesobuthus eupeus scorpion venom with focus on Kv1.1: An integrated transcriptomic and proteomic study. J. Biol. Chem. 290, 12195–12209. doi: 10.1074/jbc.M115.637611

PubMed Abstract | CrossRef Full Text | Google Scholar

Kuzmenkov, A., II, Krylov, N. A., Chugunov, A. O., Grishin, E. V., Vassilevski, A. A. (2016a). Kalium: a database of potassium channel toxins from scorpion venom. Database (Oxford) 2016, pii: baw056. doi: 10.1093/database/baw056

CrossRef Full Text | Google Scholar

Kuzmenkov, A., II, Sachkova, M. Y., Kovalchuk, S., II, Grishin, E. V., Vassilevski, A. A. (2016b). Lachesana tarabaevi, an expert in membrane-active toxins. Biochem. J. 473, 2495–2506. doi: 10.1042/BCJ20160436

PubMed Abstract | CrossRef Full Text | Google Scholar

Kuzmenkov, A., II, Peigneur, S., Chugunov, A. O., Tabakmakher, V. M., Efremov, R. G., Tytgat, J., et al. (2017). C-Terminal residues in small potassium channel blockers OdK1 and OSK3 from scorpion venom fine-tune the selectivity. Biochim. Biophys. Acta Proteins Proteomics 1865, 465–472. doi: 10.1016/j.bbapap.2017.02.001

CrossRef Full Text | Google Scholar

Kuzmenkov, A., II, Nekrasova, O. V., Peigneur, S., Tabakmakher, V. M., Gigolaev, A. M., Fradkov, A. F., et al. (2018). KV1.2 channel-specific blocker from Mesobuthus eupeus scorpion venom: Structural basis of selectivity. Neuropharmacology 143, 228–238. doi: 10.1016/j.neuropharm.2018.09.030

PubMed Abstract | CrossRef Full Text | Google Scholar

Kuzmenkov, A., II, Peigneur, S., Tytgat, J., Vassilevski, A. A. (2019). Pharmacological Characterisation of MeKTx13-2 and MeKTx13-3, Peptide Ligands of Potassium Channels From the Scorpion Mesobuthus eupeus Venom. Russ. J. Physiol. 105, 1452–1462. doi: 10.1134/S0869813919110074

CrossRef Full Text | Google Scholar

Lebrun, B., Romi-Lebrun, R., Martin-Eauclaire, M. F., Yasuda, A., Ishiguro, M., Oyama, Y., et al. (1997). A four-disulphide-bridged toxin, with high affinity towards voltage-gated K+ channels, isolated from Heterometrus spinnifer (Scorpionidae) venom. Biochem. J. 328, 321–327. doi: 10.1042/bj3280321

PubMed Abstract | CrossRef Full Text | Google Scholar

McCoy, J., LaVallie, E. (1994). Expression and Purification of Thioredoxin Fusion Proteins. Curr. Protoc. Mol. Biol. 28, Unit16.8. doi: 10.1002/0471142727.mb1608s28

CrossRef Full Text | Google Scholar

Merkler, D. (1994). C-terminal Amidated Peptides: Production by the in Vitro Enzymatic Amidation of Glycine-Extended Peptides and the Importance of the Amide to Bioactivity. Enzyme Microb. Technol. 16, 450–456. doi: 10.1016/0141-0229(94)90014-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Mouhat, S., Jouirou, B., Mosbah, A., De Waard, M., Sabatier, J.-M. (2004). Diversity of folds in animal toxins acting on ion channels. Biochem. J. 378, 717–726. doi: 10.1042/bj20031860

PubMed Abstract | CrossRef Full Text | Google Scholar

Mouhat, S., Visan, V., Ananthakrishnan, S., Wulff, H., Andreotti, N., Grissmer, S., et al. (2005). K+ channel types targeted by synthetic OSK1, a toxin from Orthochirus scrobiculosus scorpion venom. Biochem. J. 385, 95–104. doi: 10.1042/BJ20041379

PubMed Abstract | CrossRef Full Text | Google Scholar

Norton, R. S., Chandy, K. G. (2017). Venom-derived peptide inhibitors of voltage-gated potassium channels. Neuropharmacology 127, 124–138. doi: 10.1016/j.neuropharm.2017.07.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Pardo, L. A., Stühmer, W. (2014). The roles of K(+) channels in cancer. Nat. Rev. Cancer 14, 39–48. doi: 10.1038/nrc3635

PubMed Abstract | CrossRef Full Text | Google Scholar

Peigneur, S., Billen, B., Derua, R., Waelkens, E., Debaveye, S., Béress, L., et al. (2011). A bifunctional sea anemone peptide with Kunitz type protease and potassium channel inhibiting properties. Biochem. Pharmacol. 82, 81–90. doi: 10.1016/j.bcp.2011.03.023

PubMed Abstract | CrossRef Full Text | Google Scholar

Pennington, M. W., Beeton, C., Galea, C. A., Smith, B. J., Chi, V., Monaghan, K. P., et al. (2009). Engineering a stable and selective peptide blocker of the Kv1.3 channel in T lymphocytes. Mol. Pharmacol. 75, 762–773. doi: 10.1124/mol.108.052704

PubMed Abstract | CrossRef Full Text | Google Scholar

Pennington, M., Harunur Rashid, M., Tajhya, R., Beeton, C., Kuyucak, S., Norton, R. (2012). A C-terminally Amidated Analogue of ShK Is a Potent and Selective Blocker of the Voltage-Gated Potassium Channel Kv1.3. FEBS Lett. 586, 3996–4001. doi: 10.1016/J.FEBSLET.2012.09.038

PubMed Abstract | CrossRef Full Text | Google Scholar

Pluzhnikov, K., Vassilevski, A., Korolkova, Y., Fisyunov, A., Iegorova, O., Krishtal, O., et al. (2007). ω-Lsp-IA, a novel modulator of P-type Ca2+ channels. Toxicon 50, 993–1004. doi: 10.1016/j.toxicon.2007.07.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Prosdocimi, E., Checchetto, V., Leanza, L. (2019). Targeting the Mitochondrial Potassium Channel Kv1.3 to Kill Cancer Cells: Drugs, Strategies, and New Perspectives. SLAS Discovery Adv. Life Sci. R D 24, 882–892. doi: 10.1177/2472555219864894

CrossRef Full Text | Google Scholar

Pyrkov, T., Efremov, R. (2007). A Fragment-Based Scoring Function to Re-rank ATP Docking Results. Int. J. Mol. Sci. 8, 1083–1094. doi: 10.3390/i8111083

CrossRef Full Text | Google Scholar

Pyrkov, T. V., Chugunov, A. O., Krylov, N. A., Nolde, D. E., Efremov, R. G. (2009). PLATINUM: a web tool for analysis of hydrophobic/hydrophilic organization of biomolecular complexes. Bioinformatics 25, 1201–1202. doi: 10.1093/bioinformatics/btp111

PubMed Abstract | CrossRef Full Text | Google Scholar

Rashid, M. H., Huq, R., Tanner, M. R., Chhabra, S., Khoo, K. K., Estrada, R., et al. (2015). A potent and Kv1.3-selective analogue of the scorpion toxin HsTX1 as a potential therapeutic for autoimmune diseases. Sci. Rep. 4, 4509. doi: 10.1038/srep04509

CrossRef Full Text | Google Scholar

Renisio, J. G., Romi-Lebrun, R., Blanc, E., Bornet, O., Nakajima, T., Darbon, H. (2000). Solution structure of BmKTX, a K+ blocker toxin from the Chinese scorpion Buthus Martensi. Proteins 38, 70–78. doi: 10.1002/(SICI)1097-0134(20000101)38:1<70::AID-PROT8>3.0.CO;2-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Romi-Lebrun, R., Lebrun, B., Martin-Eauclaire, M. F., Ishiguro, M., Escoubas, P., Wu, F. Q., et al. (1997). Purification, characterization, and synthesis of three novel toxins from the Chinese scorpion Buthus martensi, which act on K+ channels. Biochemistry 36, 13473–13482. doi: 10.1021/bi971044w

PubMed Abstract | CrossRef Full Text | Google Scholar

Schmitz, A., Sankaranarayanan, A., Azam, P., Schmidt-Lassen, K., Homerick, D., Hänsel, W., et al. (2005). Design of PAP-1, a Selective Small Molecule Kv1.3 Blocker, for the Suppression of Effector Memory T Cells in Autoimmune Diseases. Mol. Pharmacol. 68, 1254–1270. doi: 10.1124/mol.105.015669

PubMed Abstract | CrossRef Full Text | Google Scholar

Tabakmakher, V. M., Krylov, N. A., Kuzmenkov, A., II, Efremov, R. G., Vassilevski, A. A. (2019). Kalium 2.0, a comprehensive database of polypeptide ligands of potassium channels. Sci. Data 6, 73. doi: 10.1038/s41597-019-0074-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Takacs, Z., Toups, M., Kollewe, A., Johnson, E., Cuello, L. G., Driessens, G., et al. (2009). A designer ligand specific for Kv1.3 channels from a scorpion neurotoxin-based library. Proc. Natl. Acad. Sci. U. S. A. 106, 22211–22216. doi: 10.1073/pnas.0910123106

PubMed Abstract | CrossRef Full Text | Google Scholar

Teisseyre, A., Palko-Labuz, A., Sroda-Pomianek, K., Michalak, K. (2019). Voltage-Gated Potassium Channel Kv1.3 as a Target in Therapy of Cancer. Front. Oncol. 9, 933. doi: 10.3389/fonc.2019.00933

PubMed Abstract | CrossRef Full Text | Google Scholar

Varga, Z., Gurrola-Briones, G., Papp, F., Rodríguez de la Vega, R. C., Pedraza-Alva, G., Tajhya, R. B., et al. (2012). Vm24, a natural immunosuppressive peptide, potently and selectively blocks Kv1.3 potassium channels of human T cells. Mol. Pharmacol. 82, 372–382. doi: 10.1124/mol.112.078006

PubMed Abstract | CrossRef Full Text | Google Scholar

Webb, B., Sali, A. (2016). Comparative Protein Structure Modeling Using MODELLER. Curr. Protoc. Bioinf. 54, 5.6.1–5.6.37. doi: 10.1002/cpbi.3

CrossRef Full Text | Google Scholar

Wulff, H., Zhorov, B. S. (2008). K+ channel modulators for the treatment of neurological disorders and autoimmune diseases. Chem. Rev. 108, 1744–1773. doi: 10.1021/cr078234p

PubMed Abstract | CrossRef Full Text | Google Scholar

Wulff, H., Christophersen, P., Colussi, P., Chandy, K. G., Yarov-Yarovoy, V. (2019). Antibodies and venom peptides: new modalities for ion channels. Nat. Rev. Drug Discovery 18, 339–357. doi: 10.1038/s41573-019-0013-8

CrossRef Full Text | Google Scholar

Keywords: scorpion venom, neurotoxin, voltage-gated potassium channel, potassium channel blocker, molecular modeling, molecular dynamics

Citation: Gigolaev AM, Kuzmenkov AI, Peigneur S, Tabakmakher VM, Pinheiro-Junior EL, Chugunov AO, Efremov RG, Tytgat J and Vassilevski AA (2020) Tuning Scorpion Toxin Selectivity: Switching From KV1.1 to KV1.3. Front. Pharmacol. 11:1010. doi: 10.3389/fphar.2020.01010

Received: 14 April 2020; Accepted: 23 June 2020;
Published: 07 July 2020.

Edited by:

Heike Wulff, University of California, Davis, United States

Reviewed by:

Yonghua Ji, Shanghai University, China
Craig Doupnik, USF Health, United States

Copyright © 2020 Gigolaev, Kuzmenkov, Peigneur, Tabakmakher, Pinheiro-Junior, Chugunov, Efremov, Tytgat and Vassilevski. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Alexander A. Vassilevski, avas@ibch.ru

These authors have contributed equally to this work