Design of bioactive peptides from naturally occurring μ-conotoxin structures

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.

Synthesis and insecticidal activity of cysteine-free conopeptides from Conus betulinus

The cone snail Conus betulinus is a vermivorous species that is widely distributed in the South China Sea. Its crude venom contains various peptides used to prey on marine worms. In previous studies, a systematic analysis of the peptide toxin sequences from C. betulinus was carried out using a multiomics technique. In this study, 10 cysteine-free peptides that may possess insecticidal activity were selected from a previously constructed conopeptide library of C. betulinus using the CPY-Fe conopeptide as a template. These conopeptides were prepared by solid-phase peptide synthesis (SPPS), then characterized by the reverse-phase high performance liquid chromatography (HPLC) and mass spectrometry. Insect cytotoxicity and injection experiments revealed that these cysteine-free peptides exerted favorable insecticidal effects, and two of them (Bt010 and Bt016) exhibited high insecticidal efficacy with LD50 of 9.07 nM and 10.93 nM, respectively. In addition, the 3D structures of these peptides were predicted by homology modeling, and a phylogenetic tree was constructed based on the nucleotide data of conopeptides to analyze the relationships among structures, functions, and evolution. A preliminary mechanism for the insecticidal activity of the cysteine-free conopeptides was predicted by molecular docking. To the best of our knowledge, this is the first study to report the insecticidal activity of cysteine-free conopeptides derived from Conus betulinus, signaling that they could potentially be developed into bioinsecticides with desirable properties such as easy preparation, low cost, and high potency.

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.

µ-Conotoxins Targeting the Human Voltage-Gated Sodium Channel Subtype NaV1.7

µ-Conotoxins are small, potent, peptide voltage-gated sodium (Na_V) channel inhibitors characterised by a conserved cysteine framework. Despite promising in vivo studies indicating analgesic potential of these compounds, selectivity towards the therapeutically relevant subtype Na_V1.7 has so far been limited. We recently identified a novel µ-conotoxin, SxIIIC, which potently inhibits human Na_V1.7 (hNa_V1.7). SxIIIC has high sequence homology with other µ-conotoxins, including SmIIIA and KIIIA, yet shows different Na_V channel selectivity for mammalian subtypes. Here, we evaluated and compared the inhibitory potency of µ-conotoxins SxIIIC, SmIIIA and KIIIA at hNa_V channels by whole-cell patch-clamp electrophysiology and discovered that these three closely related µ-conotoxins display unique selectivity profiles with significant variations in inhibitory potency at hNa_V1.7. Analysis of other µ-conotoxins at hNa_V1.7 shows that only a limited number are capable of inhibition at this subtype and that differences between the number of residues in loop 3 appear to influence the ability of µ-conotoxins to inhibit hNa_V1.7. Through mutagenesis studies, we confirmed that charged residues in this region also affect the selectivity for hNa_V1.4. Comparison of µ-conotoxin NMR solution structures identified differences that may contribute to the variance in hNa_V1.7 inhibition and validated the role of the loop 1 extension in SxIIIC for improving potency at hNa_V1.7, when compared to KIIIA. This work could assist in designing µ-conotoxin derivatives specific for hNa_V1.7.

Small cyclic sodium channel inhibitors

Voltage-gated sodium (Na_V) channels play crucial roles in a range of (patho)physiological processes. Much interest has arisen within the pharmaceutical industry to pursue these channels as analgesic targets following overwhelming evidence that Na_V channel subtypes Na_V1.7-Na_V1.9 are involved in nociception. More recently, Na_V1.1, Na_V1.3 and Na_V1.6 have also been identified to be involved in pain pathways. Venom-derived disulfide-rich peptide toxins, isolated from spiders and cone snails, have been used extensively as probes to investigate these channels and have attracted much interest as drug leads. However, few peptide-based leads have made it as drugs due to unfavourable physiochemical attributes including poor in vivo pharmacokinetics and limited oral bioavailability. The present work aims to bridge the gap in the development pipeline between drug leads and drug candidates by downsizing these larger venom-derived Na_V inhibitors into smaller, more "drug-like" molecules. Here, we use molecular engineering of small cyclic peptides to aid in the determination of what drives subtype selectivity and molecular interactions of these downsized inhibitors across Na_V subtypes. We designed a series of small, stable and novel Na_V probes displaying Na_V subtype selectivity and potency in vitro coupled with potent in vivo analgesic activity, involving yet to be elucidated analgesic pathways in addition to Na_V subtype modulation.

Aspartic Acid Isomerization Characterized by High Definition Mass Spectrometry Significantly Alters the Bioactivity of a Novel Toxin from Poecilotheria

Research in toxinology has created a pharmacological paradox. With an estimated 220,000 venomous animals worldwide, the study of peptidyl toxins provides a vast number of effector molecules. However, due to the complexity of the protein-protein interactions, there are fewer than ten venom-derived molecules on the market. Structural characterization and identification of post-translational modifications are essential to develop biological lead structures into pharmaceuticals. Utilizing advancements in mass spectrometry, we have created a high definition approach that fuses conventional high-resolution MS-MS with ion mobility spectrometry (HDMS^E) to elucidate these primary structure characteristics. We investigated venom from ten species of “tiger” spider (Genus: Poecilotheria) and discovered they contain isobaric conformers originating from non-enzymatic Asp isomerization. One conformer pair conserved in five of ten species examined, denominated PcaTX-1a and PcaTX-1b, was found to be a 36-residue peptide with a cysteine knot, an amidated C-terminus, and isoAsp33Asp substitution. Although the isomerization of Asp has been implicated in many pathologies, this is the first characterization of Asp isomerization in a toxin and demonstrates the isomerized product’s diminished physiological effects. This study establishes the value of a HDMS^E approach to toxin screening and characterization.

Where cone snails and spiders meet: design of small cyclic sodium-channel inhibitors

A 13 aa residue voltage-gated sodium (Na_V) channel inhibitor peptide, Pn, containing 2 disulfide bridges was designed by using a chimeric approach. This approach was based on a common pharmacophore deduced from sequence and secondary structural homology of 2 Na_V inhibitors: Conus kinoshitai toxin IIIA, a 14 residue cone snail peptide with 3 disulfide bonds, and Phoneutria nigriventer toxin 1, a 78 residue spider toxin with 7 disulfide bonds. As with the parent peptides, this novel Na_V channel inhibitor was active on Na_V1.2. Through the generation of 3 series of peptide mutants, we investigated the role of key residues and cyclization and their influence on Na_V inhibition and subtype selectivity. Cyclic PnCS1, a 10 residue peptide cyclized via a disulfide bond, exhibited increased inhibitory activity toward therapeutically relevant Na_V channel subtypes, including Na_V1.7 and Na_V1.9, while displaying remarkable serum stability. These peptides represent the first and the smallest cyclic peptide Na_V modulators to date and are promising templates for the development of toxin-based therapeutic agents.-Peigneur, S., Cheneval, O., Maiti, M., Leipold, E., Heinemann, S. H., Lescrinier, E., Herdewijn, P., De Lima, M. E., Craik, D. J., Schroeder, C. I., Tytgat, J. Where cone snails and spiders meet: design of small cyclic sodium-channel inhibitors.

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.

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.

Theoretical and simulation studies on voltage-gated sodium channels

Voltage-gated sodium (Na_v) channels are indispensable membrane elements for the generation and propagation of electric signals in excitable cells. The successes in the crystallographic studies on prokaryotic Na_v channels in recent years greatly promote the mechanistic investigation of these proteins and their eukaryotic counterparts. In this paper, we mainly review the progress in computational studies, especially the simulation studies, on these proteins in the past years.

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.

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

Several subtypes of voltage-gated Na+ (NaV) channels are important targets for pain management. μ-Conotoxins isolated from venoms of cone snails are potent and specific blockers of different NaV channel isoforms. The inhibitory effect of μ-conotoxins on NaV channels has been examined extensively, but the mechanism of toxin specificity has not been understood in detail. Here the known structure of μ-conotoxin PIIIA and a model of the skeletal muscle channel NaV1.4 are used to elucidate elements that contribute to the structural basis of μ-conotoxin binding and specificity. The model of NaV1.4 is constructed based on the crystal structure of the bacterial NaV channel, NaVAb. Six different binding modes, in which the side chain of each of the basic residues carried by the toxin protrudes into the selectivity filter of NaV1.4, are examined in atomic detail using molecular dynamics simulations with explicit solvent. The dissociation constants (Kd) computed for two selected binding modes in which Lys9 or Arg14 from the toxin protrudes into the filter of the channel are within 2 fold; both values in close proximity to those determined from dose response data for the block of NaV currents. To explore the mechanism of PIIIA specificity, a double mutant of NaV1.4 mimicking NaV channels resistant to μ-conotoxins and tetrodotoxin is constructed and the binding of PIIIA to this mutant channel examined. The double mutation causes the affinity of PIIIA to reduce by two orders of magnitude.

Ala-7, His-10 and Arg-12 are crucial amino acids for activity of a synthetically engineered μ-conotoxin

Cone snail toxins or conotoxins are often small cysteine-rich peptides which have shown to be highly selective ligands for a wide range of ion channels such as voltage-gated sodium channels (Na(V)s). Na(V)s participate in a wide range of electrophysiological processes. Consequently, their malfunction has been associated with numerous diseases. The development of subtype-selective modulators of Na(V)s remains highly important in the treatment of such disorders. In order to expand our knowledge in the search for novel therapeutics to treat Na(V)-related diseases, we explored the field of peptide engineering. In the current study, the impact of well considered point mutations into a bioactive peptide that was found to be a very potent and selective inhibitor of Na(V)s (i.e. Midi R2) was examined. We designed two peptides, named Midi R2[A7G] and Midi R2[H10A, R12A] which have mutations at position 7, and both 10 and 12, respectively. Electrophysiological recordings indicated that an Ala to Gly mutation at position 7 increased IC50-values from the nanomolar range to the micromolar range. For Midi R2[H10A, R12A] at a concentration of 10 μM, activity is even reduced to 0-10% for all of the tested Na(V)-channels. Circular dichroism measurements proved that overall structural conformations did not change. These findings suggest that the minimal space between the second and the third intercysteine loop of Midi R2 is the sequence RRWARDHSR and that His at position 10 and Arg at position 12 are crucial amino acids for the potency and specificity of Midi R2. In this way, new insights into the structure-activity relationships of μ-conotoxins were found.

NMR investigation of disulfide containing peptides and proteins

Peptides and proteins with disulfide bonds are abundant in all kingdoms and play essential role in many biological events. Because small disulfide-rich peptides (proteins) are usually difficult to crystallize, nuclear magnetic resonance (NMR) is by far one of the most powerful techniques for the determination of their solution structure. Besides the “static” three-dimensional structure, NMR has unique opportunities to acquire additional information about molecular dynamics and folding at atomic resolution. Nowadays it is becoming increasingly evident, that “excited”, “disordered” or “fuzzy” protein states may exhibit biological function and disulfide proteins are also promising targets for such studies. In this short two-three years overview those disulfide peptides and proteins were cited from the literature that were studied by NMR. Though we may have missed some, their structural diversity and complexity as well as their wide repertoire of biological functions is impressive. We emphasised especially antimicrobial peptides and peptide based toxins in addition to some biologically important other structures. Besides the general NMR methods we reviewed some contemporary techniques suitable for disclosing the peculiar properties of disulfide bonds. Interesting dynamics and folding studies of disulfide proteins were also mentioned. It is important to disclose the essential structure, dynamics, function aspects of disulfide proteins since this aids the design of new compounds with improved activity and reduced toxicity. Undoubtedly, NMR has the potential to accelerate the development of new disulfide peptides/proteins with pharmacological activity.

Strategies for the development of conotoxins as new therapeutic leads

Peptide toxins typically bind to their target ion channels or receptors with high potency and selectivity, making them attractive leads for therapeutic development. In some cases the native peptide as it is found in the venom from which it originates can be used directly, but in many instances it is desirable to truncate and/or stabilize the peptide to improve its therapeutic properties. A complementary strategy is to display the key residues that make up the pharmacophore of the peptide toxin on a non-peptidic scaffold, thereby creating a peptidomimetic. This review exemplifies these approaches with peptide toxins from marine organisms, with a particular focus on conotoxins.

Expanding chemical diversity of conotoxins: peptoid-peptide chimeras of the sodium channel blocker μ-KIIIA and its selenopeptide analogues

The μ-conotoxin KIIIA is a three disulfide-bridged blocker of voltage-gated sodium channels (VGSCs). The Lys(7) residue in KIIIA is an attractive target for manipulating the selectivity and efficacy of this peptide. Here, we report the design and chemical synthesis of μ-conopeptoid analogues (peptomers) in which we replaced Lys(7) with peptoid monomers of increasing side-chain size: N-methylglycine, N-butylglycine and N-octylglycine. In the first series of analogues, the peptide core contained all three disulfide bridges; whereas in the second series, a disulfide-depleted selenoconopeptide core was used to simplify oxidative folding. The analogues were tested for functional activity in blocking the Nav1.2 subtype of mammalian VGSCs exogenously expressed in Xenopus oocytes. All six analogues were active, with the N-methylglycine analogue, [Sar(7)]KIIIA, the most potent in blocking the channels while favouring lower efficacy. Our findings demonstrate that the use of N-substituted Gly residues in conotoxins show promise as a tool to optimize their pharmacological properties as potential analgesic drug leads.

Mammalian neuronal sodium channel blocker μ-conotoxin BuIIIB has a structured N-terminus that influences potency

Among the μ-conotoxins that block vertebrate voltage-gated sodium channels (VGSCs), some have been shown to be potent analgesics following systemic administration in mice. We have determined the solution structure of a new representative of this family, μ-BuIIIB, and established its disulfide connectivities by direct mass spectrometric collision induced dissociation fragmentation of the peptide with disulfides intact. The major oxidative folding product adopts a 1-4/2-5/3-6 pattern with the following disulfide bridges: Cys5-Cys17, Cys6-Cys23, and Cys13-Cys24. The solution structure reveals that the unique N-terminal extension in μ-BuIIIB, which is also present in μ-BuIIIA and μ-BuIIIC but absent in other μ-conotoxins, forms part of a short α-helix encompassing Glu3 to Asn8. This helix is packed against the rest of the toxin and stabilized by the Cys5-Cys17 and Cys6-Cys23 disulfide bonds. As such, the side chain of Val1 is located close to the aromatic rings of Trp16 and His20, which are located on the canonical helix that displays several residues found to be essential for VGSC blockade in related μ-conotoxins. Mutations of residues 2 and 3 in the N-terminal extension enhanced the potency of μ-BuIIIB for NaV1.3. One analogue, [d-Ala2]BuIIIB, showed a 40-fold increase, making it the most potent peptide blocker of this channel characterized to date and thus a useful new tool with which to characterize this channel. On the basis of previous results for related μ-conotoxins, the dramatic effects of mutations at the N-terminus were unanticipated and suggest that further gains in potency might be achieved by additional modifications of this region.

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.

Distinct disulfide isomers of μ-conotoxins KIIIA and KIIIB block voltage-gated sodium channels

In the preparation of synthetic conotoxins containing multiple disulfide bonds, oxidative folding can produce numerous permutations of disulfide bond connectivities. Establishing the native disulfide connectivities thus presents a significant challenge when the venom-derived peptide is not available, as is increasingly the case when conotoxins are identified from cDNA sequences. Here, we investigate the disulfide connectivity of μ-conotoxin KIIIA, which was predicted originally to have a [C1-C9,C2-C15,C4-C16] disulfide pattern based on homology with closely related μ-conotoxins. The two major isomers of synthetic μ-KIIIA formed during oxidative folding were purified and their disulfide connectivities mapped by direct mass spectrometric collision-induced dissociation fragmentation of the disulfide-bonded polypeptides. Our results show that the major oxidative folding product adopts a [C1-C15,C2-C9,C4-C16] disulfide connectivity, while the minor product adopts a [C1-C16,C2-C9,C4-C15] connectivity. Both of these peptides were potent blockers of Na(V)1.2 (K(d) values of 5 and 230 nM, respectively). The solution structure for μ-KIIIA based on nuclear magnetic resonance data was recalculated with the [C1-C15,C2-C9,C4-C16] disulfide pattern; its structure was very similar to the μ-KIIIA structure calculated with the incorrect [C1-C9,C2-C15,C4-C16] disulfide pattern, with an α-helix spanning residues 7-12. In addition, the major folding isomers of μ-KIIIB, an N-terminally extended isoform of μ-KIIIA identified from its cDNA sequence, were isolated. These folding products had the same disulfide connectivities as μ-KIIIA, and both blocked Na(V)1.2 (K(d) values of 470 and 26 nM, respectively). Our results establish that the preferred disulfide pattern of synthetic μ-KIIIA and μ-KIIIB folded in vitro is 1-5/2-4/3-6 but that other disulfide isomers are also potent sodium channel blockers. These findings raise questions about the disulfide pattern(s) of μ-KIIIA in the venom of Conus kinoshitai; indeed, the presence of multiple disulfide isomers in the venom could provide a means of further expanding the snail's repertoire of active peptides.

Conotoxins targeting neuronal voltage-gated sodium channel subtypes: potential analgesics?

Voltage-gated sodium channels (VGSC) are the primary mediators of electrical signal amplification and propagation in excitable cells. VGSC subtypes are diverse, with different biophysical and pharmacological properties, and varied tissue distribution. Altered VGSC expression and/or increased VGSC activity in sensory neurons is characteristic of inflammatory and neuropathic pain states. Therefore, VGSC modulators could be used in prospective analgesic compounds. VGSCs have specific binding sites for four conotoxin families: μ-, μO-, δ- and ί-conotoxins. Various studies have identified that the binding site of these peptide toxins is restricted to well-defined areas or domains. To date, only the μ- and μO-family exhibit analgesic properties in animal pain models. This review will focus on conotoxins from the μ- and μO-families that act on neuronal VGSCs. Examples of how these conotoxins target various pharmacologically important neuronal ion channels, as well as potential problems with the development of drugs from conotoxins, will be discussed.