Novel interactions identified between micro -Conotoxin and the Na+ channel domain I P-loop: implications for toxin-pore binding geometry

Historical Perspective of the Characterization of Conotoxins Targeting Voltage-Gated Sodium Channels

Marine toxins have potent actions on diverse sodium ion channels regulated by transmembrane voltage (voltage-gated ion channels) or by neurotransmitters (nicotinic acetylcholine receptor channels). Studies of these toxins have focused on varied aspects of venom peptides ranging from evolutionary relationships of predator and prey, biological actions on excitable tissues, potential application as pharmacological intervention in disease therapy, and as part of multiple experimental approaches towards an understanding of the atomistic characterization of ion channel structure. This review examines the historical perspective of the study of conotoxin peptides active on sodium channels gated by transmembrane voltage, which has led to recent advances in ion channel research made possible with the exploitation of the diversity of these marine toxins.

Structural and Functional Analyses of Cone Snail Toxins

Cone snails are marine gastropod mollusks with one of the most powerful venoms in nature. The toxins, named conotoxins, must act quickly on the cone snails´ prey due to the fact that snails are extremely slow, reducing their hunting capability. Therefore, the characteristics of conotoxins have become the object of investigation, and as a result medicines have been developed or are in the trialing process. Conotoxins interact with transmembrane proteins, showing specificity and potency. They target ion channels and ionotropic receptors with greater regularity, and when interaction occurs, there is immediate physiological decompensation. In this review we aimed to evaluate the structural features of conotoxins and the relationship with their target types.

Predicting Structural Details of the Sodium Channel Pore Basing on Animal Toxin Studies

Eukaryotic voltage-gated sodium channels play key roles in physiology and are targets for many toxins and medically important drugs. Physiology, pharmacology, and general architecture of the channels has long been the subject of intensive research in academia and industry. In particular, animal toxins such as tetrodotoxin, saxitoxin, and conotoxins have been used as molecular probes of the channel structure. More recently, X-ray structures of potassium and prokaryotic sodium channels allowed elaborating models of the toxin-channel complexes that integrated data from biophysical, electrophysiological, and mutational studies. Atomic level cryo-EM structures of eukaryotic sodium channels, which became available in 2017, show that the selectivity filter structure and other important features of the pore domain have been correctly predicted. This validates further employments of toxins and other small molecules as sensitive probes of fine structural details of ion channels.

The Role of Individual Disulfide Bonds of μ-Conotoxin GIIIA in the Inhibition of NaV1.4

μ-Conotoxin GIIIA, a peptide toxin isolated from Conus geographus, preferentially blocks the skeletal muscle sodium channel Na_V1.4. GIIIA folds compactly to a pyramidal structure stabilized by three disulfide bonds. To assess the contributions of individual disulfide bonds of GIIIA to the blockade of Na_V1.4, seven disulfide-deficient analogues were prepared and characterized, each with one, two, or three pairs of disulfide-bonded Cys residues replaced with Ala. The inhibitory potency of the analogues against Na_V1.4 was assayed by whole cell patch-clamp on rNa_V1.4, heterologously expressed in HEK293 cells. The corresponding IC₅₀ values were 0.069 ± 0.005 μM for GIIIA, 2.1 ± 0.3 μM for GIIIA-1, 3.3 ± 0.2 μM for GIIIA-2, and 15.8 ± 0.8 μM for GIIIA-3 (-1, -2 and -3 represent the removal of disulfide bridges Cys3–Cys15, Cys4–Cys20 and Cys10–Cys21, respectively). Other analogues were not active enough for IC₅₀ measurement. Our results indicate that all three disulfide bonds of GIIIA are required to produce effective inhibition of Na_V1.4, and the removal of any one significantly lowers its sodium channel binding affinity. Cys10–Cys21 is the most important for the Na_V1.4 potency.

Computational Study of Binding of μ-Conotoxin GIIIA to Bacterial Sodium Channels NaVAb and NaVRh

Structures of several voltage-gated sodium (NaV) channels from bacteria have been determined recently, but the same feat might not be achieved for the mammalian counterparts in the near future. Thus, at present, computational studies of the mammalian NaV channels have to be performed using homology models based on the bacterial crystal structures. A successful homology model for the mammalian NaV1.4 channel was recently constructed using the extensive mutation data for binding of μ-conotoxin GIIIA to NaV1.4, which was further validated through studies of binding of other μ-conotoxins and ion permeation. Understanding the similarities and differences between the bacterial and mammalian NaV channels is an important issue, and the NaV1.4-GIIIA system provides a good opportunity for such a comparison. To this end, we study the binding of GIIIA to the bacterial channels NaVAb and NaVRh. The complex structures are obtained using docking and molecular dynamics simulations, and the dissociation of GIIIA is studied through umbrella sampling simulations. The results are compared to those obtained from the NaV1.4-GIIIA system, and the differences in the binding modes arising from the changes in the selectivity filters are highlighted.

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.

Folding similarity of the outer pore region in prokaryotic and eukaryotic sodium channels revealed by docking of conotoxins GIIIA, PIIIA, and KIIIA in a NavAb-based model of Nav1.4

Analyses of toxin binding to a homology model of Nav1.4 indicate similar folding of the outer pore region in eukaryotic and prokaryotic sodium channels.

Voltage-gated sodium channels are targets for many drugs and toxins. However, the rational design of medically relevant channel modulators is hampered by the lack of x-ray structures of eukaryotic channels. Here, we used a homology model based on the x-ray structure of the NavAb prokaryotic sodium channel together with published experimental data to analyze interactions of the μ-conotoxins GIIIA, PIIIA, and KIIIA with the Nav1.4 eukaryotic channel. Using Monte Carlo energy minimizations and published experimentally defined pairwise contacts as distance constraints, we developed a model in which specific contacts between GIIIA and Nav1.4 were readily reproduced without deformation of the channel or toxin backbones. Computed energies of specific interactions between individual residues of GIIIA and the channel correlated with experimental estimates. The predicted complexes of PIIIA and KIIIA with Nav1.4 are consistent with a large body of experimental data. In particular, a model of Nav1.4 interactions with KIIIA and tetrodotoxin (TTX) indicated that TTX can pass between Nav1.4 and channel-bound KIIIA to reach its binding site at the selectivity filter. Our models also allowed us to explain experimental data that currently lack structural interpretations. For instance, consistent with the incomplete block observed with KIIIA and some GIIIA and PIIIA mutants, our computations predict an uninterrupted pathway for sodium ions between the extracellular space and the selectivity filter if at least one of the four outer carboxylates is not bound to the toxin. We found a good correlation between computational and experimental data on complete and incomplete channel block by native and mutant toxins. Thus, our study suggests similar folding of the outer pore region in eukaryotic and prokaryotic sodium channels.

Systematic study of binding of μ-conotoxins to the sodium channel NaV1.4

Voltage-gated sodium channels (NaV) are fundamental components of the nervous system. Their dysfunction is implicated in a number of neurological disorders, such as chronic pain, making them potential targets for the treatment of such disorders. The prominence of the NaV channels in the nervous system has been exploited by venomous animals for preying purposes, which have developed toxins that can block the NaV channels, thereby disabling their function. Because of their potency, such toxins could provide drug leads for the treatment of neurological disorders associated with NaV channels. However, most toxins lack selectivity for a given target NaV channel, and improving their selectivity profile among the NaV1 isoforms is essential for their development as drug leads. Computational methods will be very useful in the solution of such design problems, provided accurate models of the protein-ligand complex can be constructed. Using docking and molecular dynamics simulations, we have recently constructed a model for the NaV1.4-μ-conotoxin-GIIIA complex and validated it with the ample mutational data available for this complex. Here, we use the validated NaV1.4 model in a systematic study of binding other μ-conotoxins (PIIIA, KIIIA and BuIIIB) to NaV1.4. The binding mode obtained for each complex is shown to be consistent with the available mutation data and binding constants. We compare the binding modes of PIIIA, KIIIA and BuIIIB to that of GIIIA and point out the similarities and differences among them. The detailed information about NaV1.4-μ-conotoxin interactions provided here will be useful in the design of new NaV channel blocking peptides.

Molecular dynamics study of binding of µ-conotoxin GIIIA to the voltage-gated sodium channel Na(v)1.4

Homology models of mammalian voltage-gated sodium (NaV) channels based on the crystal structures of the bacterial counterparts are needed to interpret the functional data on sodium channels and understand how they operate. Such models would also be invaluable in structure-based design of therapeutics for diseases involving sodium channels such as chronic pain and heart diseases. Here we construct a homology model for the pore domain of the NaV1.4 channel and use the functional data for the binding of µ-conotoxin GIIIA to NaV1.4 to validate the model. The initial poses for the NaV1.4-GIIIA complex are obtained using the HADDOCK protein docking program, which are then refined in molecular dynamics simulations. The binding mode for the final complex is shown to be in broad agreement with the available mutagenesis data. The standard binding free energy, determined from the potential of mean force calculations, is also in good agreement with the experimental value. Because the pore domains of NaV1 channels are highly homologous, the model constructed for NaV1.4 will provide an excellent template for other NaV1 channels.

Interactions of disulfide-deficient selenocysteine analogs of μ-conotoxin BuIIIB with the α-subunit of the voltage-gated sodium channel subtype 1.3

Inhibitors of the α-subunit of the voltage-gated sodium channel subtype 1.3 (NaV 1.3) are of interest as pharmacological tools for the study of neuropathic pain associated with spinal cord injury and have potential therapeutic applications. The recently described μ-conotoxin BuIIIB (μ-BuIIIB) from Conus bullatus was shown to block NaV 1.3 with submicromolar potency (Kd = 0.2 μm), making it one of the most potent peptidic inhibitors of this subtype described to date. However, oxidative folding of μ-BuIIIB results in numerous folding isoforms, making it difficult to obtain sufficient quantities of the active form of the peptide for detailed structure-activity studies. In the present study, we report the synthesis and characterization of μ-BuIIIB analogs incorporating a disulfide-deficient, diselenide-containing scaffold designed to simplify synthesis and facilitate structure-activity studies directed at identifying amino acid residues involved in NaV 1.3 blockade. Our results indicate that, similar to other μ-conotoxins, the C-terminal residues (Trp16, Arg18 and His20) are most crucial for NaV 1 blockade. At the N-terminus, replacement of Glu3 by Ala resulted in an analog with an increased potency for NaV 1.3 (Kd = 0.07 μm), implicating this position as a potential site for modification for increased potency and/or selectivity. Further examination of this position showed that increased negative charge, through γ-carboxyglutamate replacement, decreased potency (Kd = 0.33 μm), whereas replacement with positively-charged 2,4-diamonobutyric acid increased potency (Kd = 0.036 μm). These results provide a foundation for the design and synthesis of μ-BuIIIB-based analogs with increased potency against NaV 1.3.

Probing gating mechanisms of sodium channels using pore blockers

Several classes of small molecules and peptides bind at the central pore of voltage-gated sodium channels either from the extracellular or intracellular side of the membrane and block ion conduction through the pore. Biophysical studies that shed light on the chemical nature, accessibility, and kinetics of binding of these naturally occurring and synthetic compounds reveal a wealth of information about how these channels gate. Here, we discuss insights into the structural underpinnings of gating of the channel pore and its coupling to the voltage sensors obtained from pore blockers including site 1 neurotoxins and local anesthetics.

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.

Developing a comparative docking protocol for the prediction of peptide selectivity profiles: investigation of potassium channel toxins

During the development of selective peptides against highly homologous targets, a reliable tool is sought that can predict information on both mechanisms of binding and relative affinities. These tools must first be tested on known profiles before application on novel therapeutic candidates. We therefore present a comparative docking protocol in HADDOCK using critical motifs, and use it to “predict” the various selectivity profiles of several major αKTX scorpion toxin families versus K_v1.1, K_v1.2 and K_v1.3. By correlating results across toxins of similar profiles, a comprehensive set of functional residues can be identified. Reasonable models of channel-toxin interactions can be then drawn that are consistent with known affinity and mutagenesis. Without biological information on the interaction, HADDOCK reproduces mechanisms underlying the universal binding of αKTX-2 toxins, and K_v1.3 selectivity of αKTX-3 toxins. The addition of constraints encouraging the critical lysine insertion confirms these findings, and gives analogous explanations for other families, including models of partial pore-block in αKTX-6. While qualitatively informative, the HADDOCK scoring function is not yet sufficient for accurate affinity-ranking. False minima in low-affinity complexes often resemble true binding in high-affinity complexes, despite steric/conformational penalties apparent from visual inspection. This contamination significantly complicates energetic analysis, although it is usually possible to obtain correct ranking via careful interpretation of binding-well characteristics and elimination of false positives. Aside from adaptations to the broader potassium channel family, we suggest that this strategy of comparative docking can be extended to other channels of interest with known structure, especially in cases where a critical motif exists to improve docking effectiveness.

Importance of position 8 in μ-conotoxin KIIIA for voltage-gated sodium channel selectivity

μ-Conotoxin KIIIA from Conus kinoshitai is a 16-residue peptide that acts as a potent pore blocker of several voltage-gated sodium channels (Na(v)). In order to obtain more selective blockers and to investigate the role of Trp at position 8, we substituted this residue with Arg, Gln and Glu. KIIIA and analogues were tested on a range of Na(v) expressed in Xenopus laevis oocytes. The rank order of potency for KIIIA was: rNa(v)1.4 ≥ rNa(v)1.2 > mNa(v)1.6 > rNa(v)1.3, with IC(50) values of 48 ± 6 nm, 61 ± 5 nm, 183 ± 31 nm and 3.6 ± 0.3 μm, respectively, whereas no effect was seen on hNa(v)1.5 and hNa(v)1.8 at a concentration of 10 μm. Replacement of Trp8 resulted in more selective blockers with a preference for neuronal sodium channels over the skeletal sodium channel. The activity on rNa(v)1.4 was reduced about 40-, 70- and 200-fold for [W8R]KIIIA, [W8Q]KIIIA and [W8E]KIIIA, respectively. All analogues showed a completely reversible block of rNa(v)1.2, as opposed to the partial reversibility of KIIIA. At saturating concentrations, complete block of rNa(v)1.2 was never achieved. The residual current was lower than 10%, except for [W8E]KIIIA. KIIIA had no effect on the voltage dependence of activation of rNa(v)1.2, whereas all analogues caused a depolarizing shift. Overall, this study shows that Trp8 is a key residue in the pharmacophore. Replacement of Trp8 enables more selective blockers to be obtained for neuronal sodium channels. Trp is a key determinant for the reversibility of block of rNa(v)1.2.

Cooccupancy of the outer vestibule of voltage-gated sodium channels by micro-conotoxin KIIIA and saxitoxin or tetrodotoxin

The guanidinium alkaloids tetrodotoxin (TTX) and saxitoxin (STX) are classic ligands of voltage-gated sodium channels (VGSCs). Like TTX and STX, micro-conotoxin peptides are pore blockers but with greater VGSC subtype selectivity. micro-Conotoxin KIIIA blocks the neuronal subtype Na(V)1.2 with nanomolar affinity and we recently discovered that KIIIA and its mutant with one fewer positive charge, KIIIA[K7A], could act synergistically with TTX in a ternary peptide x TTX x Na(V) complex. In the complex, the peptide appeared to trap TTX in its normal binding site such that TTX could not readily dissociate from the channel until the peptide had done so; in turn, the presence of TTX accelerated the rate at which peptide dissociated from the channel. In the present study we examined the inhibition of Na(V)1.2, exogenously expressed in Xenopus oocytes, by STX (a divalent cation) and its sulfated congener GTX2/3 (with a net +1 charge). Each could form a ternary complex with KIIIA and Na(V)1.2, as previously found with TTX (a monovalent cation), but only when STX or GTX2/3 was added before KIIIA. The KIIIA x alkaloid x Na(V) complex was considerably less stable with STX than with either GTX2/3 or TTX. In contrast, ternary KIIIA[K7A] x alkaloid x Na(V) complexes could be formed with either order of ligand addition and were about equally stable with STX, GTX2/3, or TTX. The most parsimonious interpretation of the overall results is that the alkaloid and peptide are closely apposed in the ternary complex. The demonstration that two interacting ligands ("syntoxins") occupy adjacent sites raises the possibility of evolving a much more sophisticated neuropharmacology of VGSCs.

Docking of mu-conotoxin GIIIA in the sodium channel outer vestibule

mu-Conotoxin GIIIA (mu-CTX) is a high-affinity ligand for the outer vestibule of selected isoforms of the voltage-gated Na(+) channel. The detailed bases for the toxin's high affinity binding and isoform selectivity are unclear. The outer vestibule is lined by four pore-forming (P) loops, each with an acidic residue near the mouth of the vestibule. mu-CTX has seven positively charged residues that may interact with these acidic P-loop residues. Using pair-wise alanine replacement of charged toxin and channel residues, in conjunction with double mutant cycle analysis, we determined coupling energies for specific interactions between each P-loop acidic residue and selected toxin residues to systematically establish quantitative restraints on the toxin orientation in the outer vestibule. Xenopus oocytes were injected with the mutant or native Na(+) channel mRNA, and currents measured by two-electrode voltage clamp. Mutant cycle analysis revealed novel, strong, toxin-channel interactions between K9/E403, K11/D1241, K11/D1532, and R19/D1532. Experimentally determined coupling energies for interacting residue pairs provided restraints for molecular dynamics simulations of mu-CTX docking. Our simulations suggest a refined orientation of the toxin in the pore, with toxin basic side-chains playing key roles in high-affinity binding. This modeling also provides a set of testable predictions for toxin-channel interactions, hitherto not described, that may contribute to high-affinity binding and channel isoform selectivity.

Saxitoxin, a toxic marine natural product that targets a multitude of receptors

Saxitoxin (STX) was discovered early last century and can contaminate seafood and drinking water, and over time has become an invaluable research tool and an internationally regulated chemical weapon. Among natural products, toxins obtain a unique reputation from their high affinity and selectivity for their target pharmacological receptor, which for STX has long been considered to only be the voltage gated sodium channel. In recent times however, STX has been discovered to also bind to calcium and potassium channels, neuronal nitric oxide synthase, STX metabolizing enzymes and two circulatory fluid proteins, namely a transferrin-like family of proteins and a unique protein found in the blood of pufferfish.

Voltage-gated sodium and calcium channels in nerve, muscle, and heart

Voltage-gated ion channels are membrane proteins which underlie rapid electrical signals among neurons and the spread of excitation in skeletal muscle and heart. We outline some recent advances in the study of voltage-sensitive sodium and calcium channels. Investigations are providing insight into the changes in molecular conformation associated with open-closed gating of the channels, the mechanisms by which they allow only specific ion species to pass through and carry an electric current, and the pathological consequences of small perturbations in channel structure which result from genetic mutations. Determination of three-dimensional structures, coupled with molecular manipulations by site-directed mutagenesis, and parallel electrophysiological analyses of currents through the ion channels, are providing an understanding of the roles and function of these channels at an unprecedented level of molecular detail. Crucial to these advances are studies of bacterial homologues of ion channels from man and other eukaryotes, and the use of naturally occurring peptide toxins which target different ion channel types with exquisite specificity.

Using the deadly mu-conotoxins as probes of voltage-gated sodium channels

Mu-Conotoxins (mu-CTX) are potent Na channel inhibitory peptides isolated from the venom of the predatory marine snail Conus geographus. Mu-CTXs exert their biological action by physically occluding the ion-conducting pore of voltage-gated Na (Na(v)) channels with a 1:1 stoichiometry in an all-or-none fashion. This article reviews our current knowledge of the mechanism of mu-CTX and the associated structural and functional insights into its molecular target--Na(v) channels.