Charybdotoxin unbinding from the mKv1.3 potassium channel: a combined computational and experimental study

Inhibitory Effects of 2-Aminoethoxydiphenyl Borate (2-APB) on Three KV1 Channel Currents

2-Aminoethoxydiphenyl borate (2-APB), a boron-containing compound, is a multitarget compound with potential as a drug precursor and exerts various effects in systems of the human body. Ion channels are among the reported targets of 2-APB. The effects of 2-APB on voltage-gated potassium channels (K_V) have been reported, but the types of K_V channels that 2-APB inhibits and the inhibitory mechanism remain unknown. In this paper, we discovered that 2-APB acted as an inhibitor of three representative human K_V1 channels. 2-APB significantly blocked A-type Kv channel K_V1.4 in a concentration-dependent manner, with an IC₅₀ of 67.3 μM, while it inhibited the delayed outward rectifier channels K_V1.2 and K_V1.3, with IC₅₀s of 310.4 μM and 454.9 μM, respectively. Further studies on K_V1.4 showed that V549, T551, A553, and L554 at the cavity region and N-terminal played significant roles in 2-APB's effects on the K_V1.4 channel. The results also indicated the importance of fast inactivation gating in determining the different effects of 2-APB on three channels. Interestingly, a current facilitation phenomenon by a short prepulse after 2-APB application was discovered for the first time. The docked modeling revealed that 2-APB could form hydrogen bonds with different sites in the cavity region of three channels, and the inhibition constants showed a similar trend to the experimental results. These findings revealed new molecular targets of 2-APB and demonstrated that 2-APB's effects on K_V1 channels might be part of the reason for the diverse bioactivities of 2-APB in the human body and in animal models of human disease.

Insight on the interaction between the scorpion toxin blocker Discrepin on potassium voltage-gated channel Kv4.3 by molecular dynamics simulations

Discrepin is a 38-residue α-toxin extracted from the venom of the Venezuelan scorpion Tityus discrepans, which inhibits ionic transit in the voltage-dependent potassium channels (Kv) of A-type current. The effect of specific residues on the IC₅₀ between Discrepine and Kv4.3, the main component of A-type currents, is known; however, the molecular details of the toxin-channel interaction are not known. In this work, we present interaction models between Discrepin (wt) and two peptide variants (V6K/D20K and K13A) on the pore-forming domain of the Kv4.3 channel obtained from homology, docking, and molecular dynamics modeling techniques. The free energy calculations in these models correspond to the order of the experimentally determined IC₅₀ values. Our studies shed light on the role of the K13 residue as responsible for occluding the Kv4.3 selectivity filter and the importance of the V6K mutation in the approach and stabilization of toxin-channel complex interactions.Communicated by Ramaswamy H. Sarma.

Olfactory bulb-targeted quantum dot (QD) bioconjugate and Kv1.3 blocking peptide improve metabolic health in obese male mice

The olfactory system is a driver of feeding behavior, whereby olfactory acuity is modulated by the metabolic state of the individual. The excitability of the major output neurons of the olfactory bulb (OB) can be modulated through targeting a voltage-dependent potassium channel, Kv1.3, which responds to changes in metabolic factors such as insulin, glucose, and glucagon-like peptide-1. Because gene-targeted deletion or inhibition of Kv1.3 in the periphery has been found to increase energy metabolism and decrease body weight, we hypothesized that inhibition of Kv1.3 selectively in the OB could enhance excitability of the output neurons to evoke changes in energy homeostasis. We thereby employed metal-histidine coordination to self-assemble the Kv1.3 inhibitor margatoxin (MgTx) to fluorescent quantum dots (QDMgTx) as a means to label cells in vivo and test changes in neuronal excitability and metabolism when delivered to the OB. Using patch-clamp electrophysiology to measure Kv1.3 properties in heterologously expressed cells and native mitral cells in OB slices, we found that QDMgTx had a fast rate of inhibition, but with a reduced IC_50, and increased action potential firing frequency. QDMgTx was capable of labeling cloned Kv1.3 channels but was not visible when delivered to native Kv1.3 in the OB. Diet-induced obese mice were observed to reduce body weight and clear glucose more quickly following osmotic mini-pump delivery of QDMgTx/MgTx to the OB, and following MgTx delivery, they increased the use of fats as fuels (reduced respiratory exchange ratio). These results suggest that enhanced excitability of bulbar output neurons can drive metabolic responses.

Trans-toxin ion-sensitivity of charybdotoxin-blocked potassium-channels reveals unbinding transitional states

In silico and in vitro studies have made progress in understanding protein–protein complex formation; however, the molecular mechanisms for their dissociation are unclear. Protein–protein complexes, lasting from microseconds to years, often involve induced-fit, challenging computational or kinetic analysis. Charybdotoxin (CTX), a peptide from the Leiurus scorpion venom, blocks voltage-gated K⁺-channels in a unique example of binding/unbinding simplicity. CTX plugs the external mouth of K⁺-channels pore, stopping K⁺-ion conduction, without inducing conformational changes. Conflicting with a tight binding, we show that external permeant ions enhance CTX-dissociation, implying a path connecting the pore, in the toxin-bound channel, with the external solution. This sensitivity is explained if CTX wobbles between several bound conformations, producing transient events that restore the electrical and ionic trans-pore gradients. Wobbling may originate from a network of contacts in the interaction interface that are in dynamic stochastic equilibria. These partially-bound intermediates could lead to distinct, and potentially manipulable, dissociation pathways.

Computational Protein-Protein Docking Reveals the Therapeutic Potential of Kunitz-type Venom against hKv1.2 Binding Sites

BACKGROUND & OBJECTIVE

Kunitz-type venoms are bioactive proteins isolated from a wide variety of venomous animals. These venoms are involved in protease inhibitory activity or potassium channel blocking activity. Therefore, they are reported as an important source for lead drug candidates towards protease or channel associated diseases like neurological, metabolic and cardiovascular disorders.

METHODS

This study aimed to check the inhibitory action of Kunitz-type venoms against potassium channels using computational tools.

RESULTS

Among potassium channels, Human Voltage-Gated Potassium Channel 1.2 (hKv1.2) was used as a receptor whereas Kunitz-type peptides from the venoms of various species were selected as ligand dataset.

CONCLUSION

This study helped in finding the binding interface between the receptor and ligand dataset for their potential therapeutic use in treating potassium channelopathies.

Molecular Simulations of Disulfide-Rich Venom Peptides with Ion Channels and Membranes

Disulfide-rich peptides isolated from the venom of arthropods and marine animals are a rich source of potent and selective modulators of ion channels. This makes these peptides valuable lead molecules for the development of new drugs to treat neurological disorders. Consequently, much effort goes into understanding their mechanism of action. This paper presents an overview of how molecular simulations have been used to study the interactions of disulfide-rich venom peptides with ion channels and membranes. The review is focused on the use of docking, molecular dynamics simulations, and free energy calculations to (i) predict the structure of peptide-channel complexes; (ii) calculate binding free energies including the effect of peptide modifications; and (iii) study the membrane-binding properties of disulfide-rich venom peptides. The review concludes with a summary and outlook.

Computational Studies of Venom Peptides Targeting Potassium Channels

Small peptides isolated from the venom of animals are potential scaffolds for ion channel drug discovery. This review article mainly focuses on the computational studies that have advanced our understanding of how various toxins interfere with the function of K⁺ channels. We introduce the computational tools available for the study of toxin-channel interactions. We then discuss how these computational tools have been fruitfully applied to elucidate the mechanisms of action of a wide range of venom peptides from scorpions, spiders, and sea anemone.

Scorpion toxins prefer salt solutions

There is a wide variety of ion channel types with various types of blockers, making research in this field very complicated. To reduce this complexity, it is essential to study ion channels and their blockers independently. Scorpion toxins, a major class of blockers, are charged short peptides with high affinities for potassium channels. Their high selectivity and inhibitory properties make them an important pharmacological tool for treating autoimmune or nervous system disorders. Scorpion toxins typically have highly charged surfaces and-like other proteins-an intrinsic ability to bind ions (Friedman J Phys Chem B 115(29):9213-9223, 1996; Baldwin Biophys J 71(4):2056-2063, 1996; Vrbka et al. Proc Natl Acad Sci USA 103(42):15440-15444, 2006a; Vrbka et al. J Phys Chem B 110(13):7036-43, 2006b). Thus, their effects on potassium channels are usually investigated in various ionic solutions. In this work, computer simulations of protein structures were performed to analyze the structural properties of the key residues (i.e., those that are presumably involved in contact with the surfaces of the ion channels) of 12 scorpion toxins. The presence of the two most physiologically abundant cations, Na(+) and K(+), was considered. The results indicated that the ion-binding properties of the toxin residues vary. Overall, all of the investigated toxins had more stable structures in ionic solutions than in water. We found that both the number and length of elements in the secondary structure varied depending on the ionic solution used (i.e., in the presence of NaCl or KCl). This study revealed that the ionic solution should be chosen carefully before performing experiments on these toxins. Similarly, the influence of these ions should be taken into consideration in the design of toxin-based pharmaceuticals.

Computational approaches for designing potent and selective analogs of peptide toxins as novel therapeutics

Peptide toxins provide valuable therapeutic leads for many diseases. As they bind to their targets with high affinity, potency is usually ensured. However, toxins also bind to off-target receptors, causing potential side effects. Thus, a major challenge in generating drugs from peptide toxins is ensuring their specificity for their intended targets. Computational methods can play an important role in solving such design problems through construction of accurate models of receptor-toxin complexes and calculation of binding free energies. Here we review the computational methods used for this purpose and their application to toxins targeting ion channels. We describe ShK and HsTX1 toxins, high-affinity blockers of the voltage-gated potassium channel Kv1.3, which could be developed as therapeutic agents for autoimmune diseases.

Binding modes of two scorpion toxins to the voltage-gated potassium channel kv1.3 revealed from molecular dynamics

Molecular dynamics (MD) simulations are used to examine the binding modes of two scorpion toxins, margatoxin (MgTx) and hongotoxin (HgTx), to the voltage gated K+ channel, Kv1.3. Using steered MD simulations, we insert either Lys28 or Lys35 of the toxins into the selectivity filter of the channel. The MgTx-Kv1.3 complex is stable when the side chain of Lys35 from the toxin occludes the channel filter, suggesting that Lys35 is the pore-blocking residue for Kv1.3. In this complex, Lys28 of the toxin forms one additional salt bridge with Asp449 just outside the filter of the channel. On the other hand, HgTx forms a stable complex with Kv1.3 when the side chain of Lys28 but not Lys35 protrudes into the filter of the channel. A survey of all the possible favorable binding modes of HgTx-Kv1.3 is carried out by rotating the toxin at 3° intervals around the channel axis while the position of HgTx-Lys28 relative to the filter is maintained. We identify two possible favorable binding modes: HgTx-Arg24 can interact with either Asp433 or Glu420 on the vestibular wall of the channel. The dissociation constants calculated from the two binding modes of HgTx-Kv1.3 differ by approximately 20 fold, suggesting that the two modes are of similar energetics.

Molecular dynamics simulations of scorpion toxin recognition by the Ca(2+)-activated potassium channel KCa3.1

The Ca(2+)-activated channel of intermediate-conductance (KCa3.1) is a target for antisickling and immunosuppressant agents. Many small peptides isolated from animal venoms inhibit KCa3.1 with nanomolar affinities and are promising drug scaffolds. Although the inhibitory effect of peptide toxins on KCa3.1 has been examined extensively, the structural basis of toxin-channel recognition has not been understood in detail. Here, the binding modes of two selected scorpion toxins, charybdotoxin (ChTx) and OSK1, to human KCa3.1 are examined in atomic detail using molecular dynamics (MD) simulations. Employing a homology model of KCa3.1, we first determine conduction properties of the channel using Brownian dynamics and ascertain that the simulated results are in accord with experiment. The model structures of ChTx-KCa3.1 and OSK1-KCa3.1 complexes are then constructed using MD simulations biased with distance restraints. The ChTx-KCa3.1 complex predicted from biased MD is consistent with the crystal structure of ChTx bound to a voltage-gated K(+) channel. The dissociation constants (Kd) for the binding of both ChTx and OSK1 to KCa3.1 determined experimentally are reproduced within fivefold using potential of mean force calculations. Making use of the knowledge we gained by studying the ChTx-KCa3.1 complex, we attempt to enhance the binding affinity of the toxin by carrying out a theoretical mutagenesis. A mutant toxin, in which the positions of two amino acid residues are interchanged, exhibits a 35-fold lower Kd value for KCa3.1 than that of the wild-type. This study provides insight into the key molecular determinants for the high-affinity binding of peptide toxins to KCa3.1, and demonstrates the power of computational methods in the design of novel toxins.

Electrophysiological characterization of Ts6 and Ts7, K⁺ channel toxins isolated through an improved Tityus serrulatus venom purification procedure

In Brazil, Tityus serrulatus (Ts) is the species responsible for most of the scorpion related accidents. Among the Ts toxins, the neurotoxins with action on potassium channels (α-KTx) present high interest, due to their effect in the envenoming process and the ion channel specificity they display. The α-KTx toxins family is the most relevant because its toxins can be used as therapeutic tools for specific target cells. The improved isolation method provided toxins with high resolution, obtaining pure Ts6 and Ts7 in two chromatographic steps. The effects of Ts6 and Ts7 toxins were evaluated in 14 different types of potassium channels using the voltage-clamp technique with two-microelectrodes. Ts6 toxin shows high affinity for Kv1.2, Kv1.3 and Shaker IR, blocking these channels in low concentrations. Moreover, Ts6 blocks the Kv1.3 channel in picomolar concentrations with an IC₅₀ of 0.55 nM and therefore could be of valuable assistance to further designing immunosuppressive therapeutics. Ts7 toxin blocks multiple subtypes channels, showing low selectivity among the channels analyzed. This work also stands out in its attempt to elucidate the residues important for interacting with each channel and, in the near future, to model a desired drug.

Complex structures between the N-type calcium channel (CaV2.2) and ω-conotoxin GVIA predicted via molecular dynamics

The N-type voltage-gated Ca(2+) channel CaV2.2 is one of the important targets for pain management. ω-Conotoxins isolated from venoms of cone snails, which specifically inhibit CaV2.2, are promising scaffolds for novel analgesics. The inhibitory action of ω-conotoxins on CaV2.2 has been examined experimentally, but the modes of binding of the toxins to this and other related subfamilies of Ca(2+) channels are not understood in detail. Here molecular dynamics simulations are used to construct models of ω-conotoxin GVIA in complex with a homology model of the pore domain of CaV2.2. Three different binding modes in which the side chain of Lys2, Arg17, or Lys24 from the toxin protrudes into the selectivity filter of CaV2.2 are considered. In all the modes, the toxin forms a salt bridge with an aspartate residue of subunit II just above the EEEE ring of the selectivity filter. Using the umbrella sampling technique and potential of mean force calculations, the half-maximal inhibitory concentration (IC50) values are calculated to be 1.5 and 0.7 nM for the modes in which Lys2 and Arg17 occlude the ion conduction pathway, respectively. Both IC50 values compare favorably with the values of 0.04-1.0 nM determined experimentally. The similar IC50 values calculated for the different binding modes demonstrate that GVIA can inhibit CaV2.2 with alternative binding modes. Such a multiple-binding mode mechanism may be common for ω-conotoxins.

Computational methods of studying the binding of toxins from venomous animals to biological ion channels: theory and applications

The discovery of new drugs that selectively block or modulate ion channels has great potential to provide new treatments for a host of conditions. One promising avenue revolves around modifying or mimicking certain naturally occurring ion channel modulator toxins. This strategy appears to offer the prospect of designing drugs that are both potent and specific. The use of computational modeling is crucial to this endeavor, as it has the potential to provide lower cost alternatives for exploring the effects of new compounds on ion channels. In addition, computational modeling can provide structural information and theoretical understanding that is not easily derivable from experimental results. In this review, we look at the theory and computational methods that are applicable to the study of ion channel modulators. The first section provides an introduction to various theoretical concepts, including force-fields and the statistical mechanics of binding. We then look at various computational techniques available to the researcher, including molecular dynamics, brownian dynamics, and molecular docking systems. The latter section of the review explores applications of these techniques, concentrating on pore blocker and gating modifier toxins of potassium and sodium channels. After first discussing the structural features of these channels, and their modes of block, we provide an in-depth review of past computational work that has been carried out. Finally, we discuss prospects for future developments in the field.

Computational studies of marine toxins targeting ion channels

Toxins from marine animals offer novel drug leads for treatment of diseases involving ion channels. Computational methods could be very helpful in this endeavour in several ways, e.g., (i) constructing accurate models of the channel-toxin complexes using docking and molecular dynamics (MD) simulations; (ii) determining the binding free energies of toxins from umbrella sampling MD simulations; (iii) predicting the effect of mutations from free energy MD simulations. Using these methods, one can design new analogs of toxins with improved affinity and selectivity properties. Here we present a review of the computational methods and discuss their applications to marine toxins targeting potassium and sodium channels. Detailed examples from the potassium channel toxins—ShK from sea anemone and κ-conotoxin PVIIA—are provided to demonstrate capabilities of the computational methods to give accurate descriptions of the channel-toxin complexes and the energetics of their binding. An example is also given from sodium channel toxins (µ-conotoxin GIIIA) to illustrate the differences between the toxin binding modes in potassium and sodium channels.

Why the Drosophila Shaker K+ channel is not a good model for ligand binding to voltage-gated Kv1 channels

The Drosophila Shaker K(+) channel is the first cloned voltage-gated potassium channel and has, therefore, played an important role in structural and functional studies of those channels. While such a role is well justified for ion permeation, it is not clear whether this also extends to ligand binding. Despite the high degree of homology among Shaker and Kv1 channels, κ-conotoxin PVIIA (κ-PVIIA) binds to Shaker with high affinity but not to Kv1 channels. Here we address this issue by studying binding of κ-PVIIA to Shaker and Kv1 channels using molecular dynamics (MD) simulations. The structures of the channel-toxin complexes are constructed via docking and refinement with MD. The binding mode of each complex is characterized and compared to available mutagenesis data to validate the complex models. The potential of mean force for dissociation of the Shaker-κ-PVIIA complex is calculated from umbrella sampling MD simulations, and the corresponding binding free energy is determined, which provides further validation of the complex structure. Comparison of the Shaker and Kv1 complex models shows that a few mutations in the turret and extended regions are sufficient to abolish the observed sensitivity of Shaker to κ-PVIIA. This study demonstrates that Shaker is not always a good model for Kv1 channels for ligand binding. It also provides insights into the binding of the toxin to potassium channels that will be useful for improving affinity and selectivity properties of Kv1 channels.

Molecular approaches for structural characterization of a new potassium channel blocker from Tityus stigmurus venom: cDNA cloning, homology modeling, dynamic simulations and docking

Potassium channels are involved in the maintenance of resting membrane potential, control of cardiac and neuronal excitability, neurotransmitters release, muscle contractility and hormone secretion. The Tityus stigmurus scorpion is widely distributed in Northeastern Brazil and known to cause severe human envenomations, inducing pain, hypoesthesia, edema, erythema, paresthesia, headaches and vomiting. Most potassium channel blocking peptides that have been purified from scorpion venoms contain 30-40 amino acids with three or four disulfide bridges. These peptides belong to α-KTx subfamily. On the other hand, the β-KTx subfamily is poorly characterized, though it is very representative in some scorpion venoms. A transcriptomic approach of T.stigmurus scorpions developed by our group revealed the repertoire of possible molecules present in the venom, including many toxins of the β-KTx subfamily. One of the ESTs found, named TSTI0003C has a cDNA sequence of 538 bp codifying a mature protein with 47 amino acid residues, corresponding to 5299 Da. This β-KTx peptide is a new member of the BmTXKβ-related toxins, and was here named TstKMK. The three-dimensional structure of this potassium channel toxin of the T. stigmurus scorpion was obtained by computational modeling and refined by molecular dynamic simulations. Furthermore, we have made docking simulations using a Shaker kV-1.2 potassium channel from rats as receptor model and proposed which amino acid residues and interactions could be involved in its blockade.

Structural basis of the selective block of Kv1.2 by maurotoxin from computer simulations

The 34-residue polypeptide maurotoxin (MTx) isolated from scorpion venoms selectively inhibits the current of the voltage-gated potassium channel Kv1.2 by occluding the ion conduction pathway. Here using molecular dynamics simulation as a docking method, the binding modes of MTx to three closely related channels (Kv1.1, Kv1.2 and Kv1.3) are examined. We show that MTx forms more favorable electrostatic interactions with the outer vestibule of Kv1.2 compared to Kv1.1 and Kv1.3, consistent with the selectivity of MTx for Kv1.2 over Kv1.1 and Kv1.3 observed experimentally. One salt bridge in the bound complex of MTx-Kv1.2 forms and breaks in a simulation period of 20ns, suggesting the dynamic nature of toxin-channel interactions. The toxin selectivity likely arises from the differences in the shape of the channel outer vestibule, giving rise to distinct orientations of MTx on block. Potential of mean force calculations show that MTx blocks Kv1.1, Kv1.2 and Kv1.3 with an IC₅₀ value of 6 µM, 0.6nM and 18 µM, respectively.

Structure, function, and chemical synthesis of Vaejovis mexicanus peptide 24: a novel potent blocker of Kv1.3 potassium channels of human T lymphocytes

Animal venoms are rich sources of ligands for studying ion channels and other pharmacological targets. Proteomic analyses of the soluble venom from the Mexican scorpion Vaejovis mexicanus smithi showed that it contains more than 200 different components. Among them, a 36-residue peptide with a molecular mass of 3864 Da (named Vm24) was shown to be a potent blocker of Kv1.3 of human lymphocytes (K(d) ∼ 3 pM). The three-dimensional solution structure of Vm24 was determined by nuclear magnetic resonance, showing the peptide folds into a distorted cystine-stabilized α/β motif consisting of a single-turn α-helix and a three-stranded antiparallel β-sheet, stabilized by four disulfide bridges. The disulfide pairs are formed between Cys6 and Cys26, Cys12 and Cys31, Cys16 and Cys33, and Cys21 and Cys36. Sequence analyses identified Vm24 as the first example of a new subfamily of α-type K(+) channel blockers (systematic number α-KTx 23.1). Comparison with other Kv1.3 blockers isolated from scorpions suggests a number of structural features that could explain the remarkable affinity and specificity of Vm24 toward Kv1.3 channels of lymphocytes.

Charybdotoxin and margatoxin acting on the human voltage-gated potassium channel hKv1.3 and its H399N mutant: an experimental and computational comparison

The effect of the pore-blocking peptides charybdotoxin and margatoxin, both scorpion toxins, on currents through human voltage-gated hK(v)1.3 wild-type and hK(v)1.3_H399N mutant potassium channels was characterized by the whole-cell patch clamp technique. In the mutant channels, both toxins hardly blocked current through the channels, although they did prevent C-type inactivation by slowing down the current decay during depolarization. Molecular dynamics simulations suggested that the fast current decay in the mutant channel was a consequence of amino acid reorientations behind the selectivity filter and indicated that the rigidity-flexibility in that region played a key role in its interactions with scorpion toxins. A channel with a slightly more flexible selectivity filter region exhibits distinct interactions with scorpion toxins. Our studies suggest that the toxin-channel interactions might partially restore rigidity in the selectivity filter and thereby prevent the structural rearrangements associated with C-type inactivation.

Affinity and selectivity of ShK toxin for the Kv1 potassium channels from free energy simulations

The voltage-gated potassium channel Kv1.3 is an attractive target for treatment of autoimmune diseases. ShK toxin from sea anemone is one of the most potent blockers of Kv1.3, and therefore ShK and its analogues have been proposed as therapeutic leads for such diseases. Increasing the selectivity of the proposed leads for Kv1.3 over other Kv1 channels is a major issue in this endeavor. Here we study binding of ShK toxin to Kv1 channels using free energy simulation methods. Homology models for Kv1.1 and Kv1.3 channels are constructed using the crystal structure of Kv1.2. The initial poses for the Kv1.x-ShK complexes are obtained using HADDOCK, which are then refined via molecular dynamics simulations. The binding mode in each complex is characterized by identifying the strongly interacting residues, which compare well with available mutagenesis studies. For each complex, the potential of mean force is calculated from umbrella sampling simulations, and the corresponding absolute binding free energy is determined. The computed binding free energies are in good agreement with the experimental data, which increases the confidence on the model complexes. The insights gained on Kv1.x-ShK binding modes will be valuable in the development of new ShK analogues with better selectivity properties.