Mechanism of μ-conotoxin PIIIA binding to the voltage-gated Na+ channel NaV1.4

Functional effects of drugs and toxins interacting with NaV1.4

NaV1.4 is a voltage-gated sodium channel subtype that is predominantly expressed in skeletal muscle cells. It is essential for producing action potentials and stimulating muscle contraction, and mutations in NaV1.4 can cause various muscle disorders. The discovery of the cryo-EM structure of NaV1.4 in complex with β1 has opened new possibilities for designing drugs and toxins that target NaV1.4. In this review, we summarize the current understanding of channelopathies, the binding sites and functions of chemicals including medicine and toxins that interact with NaV1.4. These substances could be considered novel candidate compounds or tools to develop more potent and selective drugs targeting NaV1.4. Therefore, studying NaV1.4 pharmacology is both theoretically and practically meaningful.

Voltage-Gated Sodium Channel Inhibition by µ-Conotoxins

µ-Conotoxins are small, potent pore-blocker inhibitors of voltage-gated sodium (NaV) channels, which have been identified as pharmacological probes and putative leads for analgesic development. A limiting factor in their therapeutic development has been their promiscuity for different NaV channel subtypes, which can lead to undesirable side-effects. This review will focus on four areas of µ-conotoxin research: (1) mapping the interactions of µ-conotoxins with different NaV channel subtypes, (2) µ-conotoxin structure-activity relationship studies, (3) observed species selectivity of µ-conotoxins and (4) the effects of µ-conotoxin disulfide connectivity on activity. Our aim is to provide a clear overview of the current status of µ-conotoxin research.

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.

Preparation of Conotoxin-Encapsulated Chitosan Nanoparticles and Evaluation of Their Skin Permeability

μ-Conotoxin CnIIIC (conotoxin, CTX)-loaded chitosan nanoparticles (CTX-NPs) were prepared using the ionic cross-linking method. The CTX-NPs were spherical and well with a polydispersity index of 0.292 ± 0.039, drug loading efficiency of 25.9 ± 1.2%, and encapsulation efficiency of 95.6 ± 1.3%. In vitro release studies showed that the release behavior of CTX-NPs in a pH 5.0 acetate buffer followed zero-order kinetics. In vitro transdermal experiments using Franz diffusion cells mounted with mouse abdominal skin demonstrated that the cumulative intradermal deposition amount of CTX per unit area in 8 h (D₈) and permeability coefficient (P_f) of CTX loaded on CTX-NPs were 2.30- and 7.71-times that of the CTX solution. In vivo transdermal experiments in mice showed that the amount of CTX deposited in the skin after 8 h of CTX saline administration was significantly lower than that of CTX deposited in the skin after administration of CTX-NPs. In vitro fluorescence labeling transdermal studies through Franz diffusion cells mounted with mouse abdominal skin indicated that CTX-NPs aggregated at hair follicles. Skin irritation tests in mice indicated that the irritation due to CTX-NPs was negligible. The cytotoxicity experiment showed that the viability of Balb/c 3T3 cells with CTX-NPs containing 230 μg/mL (0.08 μM) CTX was greater than 75%. CTX-NPs increase intradermal deposition of CTX by accumulating in hair follicles, which has positive implications for transdermal penetration of CTX.

Brevetoxin and Conotoxin Interactions with Single-Domain Voltage-Gated Sodium Channels from a Diatom and Coccolithophore

The recently characterized single-domain voltage-gated ion channels from eukaryotic protists (EukCats) provide an array of novel channel proteins upon which to test the pharmacology of both clinically and environmentally relevant marine toxins. Here, we examined the effects of the hydrophilic µ-CTx PIIIA and the lipophilic brevetoxins PbTx-2 and PbTx-3 on heterologously expressed EukCat ion channels from a marine diatom and coccolithophore. Surprisingly, none of the toxins inhibited the peak currents evoked by the two EukCats tested. The lack of homology in the outer pore elements of the channel may disrupt the binding of µ-CTx PIIIA, while major structural differences between mammalian sodium channels and the C-terminal domains of the EukCats may diminish interactions with the brevetoxins. However, all three toxins produced significant negative shifts in the voltage dependence of activation and steady state inactivation, suggesting alternative and state-dependent binding conformations that potentially lead to changes in the excitability of the phytoplankton themselves.

Extremely Potent Block of Bacterial Voltage-Gated Sodium Channels by µ-Conotoxin PIIIA

µ-Conotoxin PIIIA, in the sub-picomolar, range inhibits the archetypal bacterial sodium channel NaChBac (NavBh) in a voltage- and use-dependent manner. Peptide µ-conotoxins were first recognized as potent components of the venoms of fish-hunting cone snails that selectively inhibit voltage-gated skeletal muscle sodium channels, thus preventing muscle contraction. Intriguingly, computer simulations predicted that PIIIA binds to prokaryotic channel NavAb with much higher affinity than to fish (and other vertebrates) skeletal muscle sodium channel (Nav 1.4). Here, using whole-cell voltage clamp, we demonstrate that PIIIA inhibits NavBac mediated currents even more potently than predicted. From concentration-response data, with [PIIIA] varying more than 6 orders of magnitude (10⁻¹² to 10⁻⁵ M), we estimated an IC₅₀ = ~5 pM, maximal block of 0.95 and a Hill coefficient of 0.81 for the inhibition of peak currents. Inhibition was stronger at depolarized holding potentials and was modulated by the frequency and duration of the stimulation pulses. An important feature of the PIIIA action was acceleration of macroscopic inactivation. Docking of PIIIA in a NaChBac (NavBh) model revealed two interconvertible binding modes. In one mode, PIIIA sterically and electrostatically blocks the permeation pathway. In a second mode, apparent stabilization of the inactivated state was achieved by PIIIA binding between P2 helices and trans-membrane S5s from adjacent channel subunits, partially occluding the outer pore. Together, our experimental and computational results suggest that, besides blocking the channel-mediated currents by directly occluding the conducting pathway, PIIIA may also change the relative populations of conducting (activated) and non-conducting (inactivated) states.

Effect of Conformational Diversity on the Bioactivity of µ-Conotoxin PIIIA Disulfide Isomers

Cyclic µ-conotoxin PIIIA, a potent blocker of skeletal muscle voltage-gated sodium channel NaV1.4, is a 22mer peptide stabilized by three disulfide bonds. Combining electrophysiological measurements with molecular docking and dynamic simulations based on NMR solution structures, we investigated the 15 possible 3-disulfide-bonded isomers of µ-PIIIA to relate their blocking activity at NaV1.4 to their disulfide connectivity. In addition, three µ-PIIIA mutants derived from the native disulfide isomer, in which one of the disulfide bonds was omitted (C4-16, C5-C21, C11-C22), were generated using a targeted protecting group strategy and tested using the aforementioned methods. The 3-disulfide-bonded isomers had a range of different conformational stabilities, with highly unstructured, flexible conformations with low or no channel-blocking activity, while more constrained molecules preserved 30% to 50% of the native isomer's activity. This emphasizes the importance and direct link between correct fold and function. The elimination of one disulfide bond resulted in a significant loss of blocking activity at NaV1.4, highlighting the importance of the 3-disulfide-bonded architecture for µ-PIIIA. µ-PIIIA bioactivity is governed by a subtle interplay between an optimally folded structure resulting from a specific disulfide connectivity and the electrostatic potential of the conformational ensemble.

In Silico Analysis of the Subtype Selective Blockage of KCNA Ion Channels through the µ-Conotoxins PIIIA, SIIIA, and GIIIA

Understanding subtype specific ion channel pore blockage by natural peptide-based toxins is crucial for developing such compounds into promising drug candidates. Herein, docking and molecular dynamics simulations were employed in order to understand the dynamics and binding states of the µ-conotoxins, PIIIA, SIIIA, and GIIIA, at the voltage-gated potassium channels of the KV1 family, and they were correlated with their experimental activities recently reported by Leipold et al. Their different activities can only adequately be understood when dynamic information about the toxin-channel systems is available. For all of the channel-bound toxins investigated herein, a certain conformational flexibility was observed during the molecular dynamic simulations, which corresponds to their bioactivity. Our data suggest a similar binding mode of µ-PIIIA at KV1.6 and KV1.1, in which a plethora of hydrogen bonds are formed by the Arg and Lys residues within the α-helical core region of µ-PIIIA, with the central pore residues of the channel. Furthermore, the contribution of the K+ channel’s outer and inner pore loops with respect to the toxin binding. and how the subtype specificity is induced, were proposed.

Snails In Silico: A Review of Computational Studies on the Conopeptides

Marine cone snails are carnivorous gastropods that use peptide toxins called conopeptides both as a defense mechanism and as a means to immobilize and kill their prey. These peptide toxins exhibit a large chemical diversity that enables exquisite specificity and potency for target receptor proteins. This diversity arises in terms of variations both in amino acid sequence and length, and in posttranslational modifications, particularly the formation of multiple disulfide linkages. Most of the functionally characterized conopeptides target ion channels of animal nervous systems, which has led to research on their therapeutic applications. Many facets of the underlying molecular mechanisms responsible for the specificity and virulence of conopeptides, however, remain poorly understood. In this review, we will explore the chemical diversity of conopeptides from a computational perspective. First, we discuss current approaches used for classifying conopeptides. Next, we review different computational strategies that have been applied to understanding and predicting their structure and function, from machine learning techniques for predictive classification to docking studies and molecular dynamics simulations for molecular-level understanding. We then review recent novel computational approaches for rapid high-throughput screening and chemical design of conopeptides for particular applications. We close with an assessment of the state of the field, emphasizing important questions for future lines of inquiry.

µ-Conotoxins Modulating Sodium Currents in Pain Perception and Transmission: A Therapeutic Potential

The Conus genus includes around 500 species of marine mollusks with a peculiar production of venomous peptides known as conotoxins (CTX). Each species is able to produce up to 200 different biological active peptides. Common structure of CTX is the low number of amino acids stabilized by disulfide bridges and post-translational modifications that give rise to different isoforms. µ and µO-CTX are two isoforms that specifically target voltage-gated sodium channels. These, by inducing the entrance of sodium ions in the cell, modulate the neuronal excitability by depolarizing plasma membrane and propagating the action potential. Hyperexcitability and mutations of sodium channels are responsible for perception and transmission of inflammatory and neuropathic pain states. In this review, we describe the current knowledge of µ-CTX interacting with the different sodium channels subtypes, the mechanism of action and their potential therapeutic use as analgesic compounds in the clinical management of pain conditions.

Functionalized Fullerene Targeting Human Voltage-Gated Sodium Channel, hNav1.7

Mutations of hNa_v1.7 that cause its activities to be enhanced contribute to severe neuropathic pain. Only a small number of hNa_v1.7 specific inhibitors have been identified, most of which interact with the voltage-sensing domain of the voltage-activated sodium ion channel. In our previous computational study, we demonstrated that a [Lys₆]-C₈₄ fullerene binds tightly (affinity of 46 nM) to Na_vAb, the voltage-gated sodium channel from the bacterium Arcobacter butzleri. Here, we extend this work and, using molecular dynamics simulations, demonstrate that the same [Lys₆]-C₈₄ fullerene binds strongly (2.7 nM) to the pore of a modeled human sodium ion channel hNa_v1.7. In contrast, the fullerene binds only weakly to a mutated model of hNa_v1.7 (I1399D) (14.5 mM) and a model of the skeletal muscle hNa_v1.4 (3.7 mM). Comparison of one representative sequence from each of the nine human sodium channel isoforms shows that only hNa_v1.7 possesses residues that are critical for binding the fullerene derivative and blocking the channel pore.

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.

Lysine and the Na+/K+ Selectivity in Mammalian Voltage-Gated Sodium Channels

Voltage-gated sodium (Na_v) channels are critical in the generation and transmission of neuronal signals in mammals. The crystal structures of several prokaryotic Na_v channels determined in recent years inspire the mechanistic studies on their selection upon the permeable cations (especially between Na⁺ and K⁺ ions), a property that is proposed to be mainly determined by residues in the selectivity filter. However, the mechanism of cation selection in mammalian Na_v channels lacks direct explanation at atomic level due to the difference in amino acid sequences between mammalian and prokaryotic Na_v homologues, especially at the constriction site where the DEKA motif has been identified to determine the Na⁺/K⁺ selectivity in mammalian Na_v channels but is completely absent in the prokaryotic counterparts. Among the DEKA residues, Lys is of the most importance since its mutation to Arg abolishes the Na⁺/K⁺ selectivity. In this work, we modeled the pore domain of mammalian Na_v channels by mutating the four residues at the constriction site of a prokaryotic Na_v channel (Na_vRh) to DEKA, and then mechanistically investigated the contribution of Lys in cation selection using molecular dynamics simulations. The DERA mutant was generated as a comparison to understand the loss of ion selectivity caused by the K-to-R mutation. Simulations and free energy calculations on the mutants indicate that Lys facilitates Na⁺/K⁺ selection by electrostatically repelling the cation to a highly Na⁺-selective location sandwiched by the carboxylate groups of Asp and Glu at the constriction site. In contrast, the electrostatic repulsion is substantially weakened when Lys is mutated to Arg, because of two intrinsic properties of the Arg side chain: the planar geometric design and the sparse charge distribution of the guanidine group.

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.

Investigation of α-conotoxin unbinding using umbrella sampling

α-Conotoxins, a class of short and disulfide rich peptide toxins, specifically and potently block nicotinic acetylcholine receptors (nAChRs). In this study umbrella sampling was performed to study the unbinding pathways and potential of mean force (PMF) of α-conotoxin ImI and PNIA(A10L,D14K). Our results suggest that (i) the unbinding pathways of ImI and PNIA(A10L,D14K) are similar despite of their different disulfide framework and structure, and (ii) α-conotoxin unbinding requires large conformation perturbation of the C-loop and the backbone flexibility of the C-loop can affect the binding or unbinding kinetics of the α-conotoxins. In addition, (iii) umbrella sampling gave correct ranking of the binding affinities of ImI and PNIA(A10L,D14K) indicating its efficacy on prediction of the binding affinities of α-conotoxins and implicating its potential application in design of more potent α-conotoxin analogs.

Mechanism of Ion Permeation in Mammalian Voltage-Gated Sodium Channels

Recent determination of the crystal structures of bacterial voltage-gated sodium (NaV) channels have raised hopes that modeling of the mammalian counterparts could soon be achieved. However, there are substantial differences between the pore domains of the bacterial and mammalian NaV channels, which necessitates careful validation of mammalian homology models constructed from the bacterial NaV structures. Such a validated homology model for the NaV1.4 channel was constructed recently using the extensive mutagenesis data available for binding of μ-conotoxins. Here we use this NaV1.4 model to study the ion permeation mechanism in mammalian NaV channels. Linking of the DEKA residues in the selectivity filter with residues in the neighboring domains is found to be important for keeping the permeation pathway open. Molecular dynamics simulations and potential of mean force calculations reveal that there is a binding site for a Na+ ion just inside the DEKA locus, and 1-2 Na+ ions can occupy the vestibule near the EEDD ring. These sites are separated by a low free energy barrier, suggesting that inward conduction occurs when a Na+ ion in the vestibule goes over the free energy barrier and pushes the Na+ ion in the filter to the intracellular cavity, consistent with the classical knock-on mechanism. The NaV1.4 model also provides a good description of the observed Na+/K+ selectivity.

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.

Voltage-Gated Sodium Channels: Structure, Function, Pharmacology, and Clinical Indications

The tremendous therapeutic potential of voltage-gated sodium channels (Na(v)s) has been the subject of many studies in the past and is of intense interest today. Na(v)1.7 channels in particular have received much attention recently because of strong genetic validation of their involvement in nociception. Here we summarize the current status of research in the Na(v) field and present the most relevant recent developments with respect to the molecular structure, general physiology, and pharmacology of distinct Na(v) channel subtypes. We discuss Na(v) channel ligands such as small molecules, toxins isolated from animal venoms, and the recently identified Na(v)1.7-selective antibody. Furthermore, we review eight characterized ligand binding sites on the Na(v) channel α subunit. Finally, we examine possible therapeutic applications of Na(v) ligands and provide an update on current clinical studies.

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.

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.