A new level of conotoxin diversity, a non-native disulfide bond connectivity in alpha-conotoxin AuIB reduces structural definition but increases biological activity

Conotoxin Prediction: New Features to Increase Prediction Accuracy

Conotoxins are toxic, disulfide-bond-rich peptides from cone snail venom that target a wide range of receptors and ion channels with multiple pathophysiological effects. Conotoxins have extraordinary potential for medical therapeutics that include cancer, microbial infections, epilepsy, autoimmune diseases, neurological conditions, and cardiovascular disorders. Despite the potential for these compounds in novel therapeutic treatment development, the process of identifying and characterizing the toxicities of conotoxins is difficult, costly, and time-consuming. This challenge requires a series of diverse, complex, and labor-intensive biological, toxicological, and analytical techniques for effective characterization. While recent attempts, using machine learning based solely on primary amino acid sequences to predict biological toxins (e.g., conotoxins and animal venoms), have improved toxin identification, these methods are limited due to peptide conformational flexibility and the high frequency of cysteines present in toxin sequences. This results in an enumerable set of disulfide-bridged foldamers with different conformations of the same primary amino acid sequence that affect function and toxicity levels. Consequently, a given peptide may be toxic when its cysteine residues form a particular disulfide-bond pattern, while alternative bonding patterns (isoforms) or its reduced form (free cysteines with no disulfide bridges) may have little or no toxicological effects. Similarly, the same disulfide-bond pattern may be possible for other peptide sequences and result in different conformations that all exhibit varying toxicities to the same receptor or to different receptors. We present here new features, when combined with primary sequence features to train machine learning algorithms to predict conotoxins, that significantly increase prediction accuracy.

Usage of Cell-Free Protein Synthesis in Post-Translational Modification of μ-Conopeptide PIIIA

The post-translational modifications of conopeptides are the most complicated modifications to date and are well-known and closely related to the activity of conopeptides. The hydroxylation of proline in conopeptides affects folding, structure, and biological activity, and prolyl 4 hydroxylase has been characterized in Conus literatus. However, the hydroxylation machinery of proline in conopeptides is still unclear. In order to address the hydroxylation mechanism of proline in μ-PIIIA, three recombinant plasmids encoding different hybrid precursors of μ-PIIIA were constructed and crossly combined with protein disulfide isomerase, prolyl 4 hydroxylase, and glutaminyl cyclase in a continuous exchange cell-free protein system. The findings showed that prolyl 4 hydroxylase might recognize the propeptide of μ-PIIIA to achieve the hydroxylation of proline, while the cyclization of glutamate was also formed. Additionally, in Escherichia coli, the co-expression plasmid encoding prolyl 4 hydroxylase and the precursor of μ-PIIIA containing pro and mature regions were used to validate the continuous exchange cell-free protein system. Surprisingly, in addition to the two hydroxyproline residues and one pyroglutamyl residue, three disulfide bridges were formed using Trx as a fusion tag, and the yield of the fusion peptide was approximately 20 mg/L. The results of electrophysiology analysis indicated that the recombinant μ-PIIIA without C-terminal amidate inhibited the current of hNaV1.4 with a 939 nM IC50. Our work solved the issue that it was challenging to quickly generate post-translationally modified conopeptides in vitro. This is the first study to demonstrate that prolyl 4 hydroxylase catalyzes the proline hydroxylation through recognition in the propeptide of μ-PIIIA, and it will provide a new way for synthesizing multi-modified conopeptides with pharmacological activity.

Characterization of the Native Disulfide Isomers of the Novel χ-Conotoxin PnID: Implications for Further Increasing Conotoxin Diversity

χ-Conotoxins are known for their ability to selectively inhibit norepinephrine transporters, an ability that makes them potential leads for treating various neurological disorders, including neuropathic pain. PnID, a peptide isolated from the venom of Conus pennaceus, shares high sequence homology with previously characterized χ-conotoxins. Whereas previously reported χ-conotoxins seem to only have a single native disulfide bonding pattern, PnID has three native isomers due to the formation of different disulfide bond patterns during its maturation in the venom duct. In this study, the disulfide connectivity and three-dimensional structure of these disulfide isomers were explored using regioselective synthesis, chromatographic coelution, and solution-state nuclear magnetic resonance spectroscopy. Of the native isomers, only the isomer with a ribbon disulfide configuration showed pharmacological activity similar to other χ-conotoxins. This isomer inhibited the rat norepinephrine transporter (IC50 = 10 ± 2 µM) and has the most structural similarity to previously characterized χ-conotoxins. In contrast, the globular isoform of PnID showed more than ten times less activity against this transporter and the beaded isoform did not display any measurable biological activity. This study is the first report of the pharmacological and structural characterization of an χ-conotoxin from a species other than Conus marmoreus and is the first report of the existence of natively-formed conotoxin isomers.

A Single Amino Acid Replacement Boosts the Analgesic Activity of α-Conotoxin AuIB through the Inhibition of the GABABR-Coupled N-Type Calcium Channel

α-conotoxin AuIB is the only one of the 4/6 type α-conotoxins (α-CTxs) that inhibits the γ-aminobutyric acid receptor B (GABA_BR)-coupled N-type calcium channel (Ca_V2.2). To improve its inhibitory activity, a series of variants were synthesized and evaluated according to the structure-activity relationships of 4/7 type α-CTxs targeting GABA_BR-coupled Ca_V2.2. Surprisingly, only the substitution of Pro7 with Arg results in a 2-3-fold increase in the inhibition of GABA_BR-coupled Ca_V2.2 (IC₅₀ is 0.74 nM); substitutions of position 9-12 with basic or hydrophobic amino acid and the addition of hydrophobic amino acid Leu or Ile at the second loop to mimic 4/7 type α-CTxs all failed to improve the inhibitory activity of AuIB against GABA_BR-coupled Ca_V2.2. Interestingly, the most potent form of AuIB[P7R] has disulfide bridges of "1-4, 2-3" (ribbon), which differs from the "1-3, 2-4" (globular) in the isoforms of wildtype AuIB. In addition, AuIB[P7R](globular) displays potent analgesic activity in the acetic acid writhing model and the partial sciatic nerve injury (PNL) model. Our study demonstrated that 4/6 type α-CTxs, with the disulfide bridge connectivity "1-4, 2-3," are also potent inhibitors for GABA_BR-coupled Ca_V2.2, exhibiting potent analgesic activity.

Single-Disulfide Conopeptide Czon1107, an Allosteric Antagonist of the Human α3β4 Nicotinic Acetylcholine Receptor

Conopeptides are peptides in the venom of marine cone snails that are used for capturing prey or as a defense against predators. A new cysteine-poor conopeptide, Czon1107, has exhibited non-competitive inhibition with an undefined allosteric mechanism in the human (h) α3β4 nicotinic acetylcholine receptors (nAChRs). In this study, the binding mode of Czon1107 to hα3β4 nAChR was investigated using molecular dynamics simulations coupled with mutagenesis studies of the peptide and electrophysiology studies on heterologous hα3β4 nAChRs. Overall, this study clarifies the structure–activity relationship of Czon1107 and hα3β4 nAChR and provides an important experimental and theoretical basis for the development of new peptide drugs.

Structural and functional insights into the inhibition of human voltage-gated sodium channels by μ-conotoxin KIIIA disulfide isomers

μ-Conotoxins are components of cone snail venom, well-known for their analgesic activity through potent inhibition of voltage-gated sodium channel (Na_V) subtypes, including Na_V1.7. These small, disulfide-rich peptides are typically stabilized by three disulfide bonds arranged in a ‘native’ CysI-CysIV, CysII-CysV, CysIII-CysVI pattern of disulfide connectivity. However, μ-conotoxin KIIIA, the smallest and most studied μ-conotoxin with inhibitory activity at Na_V1.7, forms two distinct disulfide bond isomers during thermodynamic oxidative folding, including Isomer 1 (CysI-CysV, CysII-CysIV, CysIII-CysVI) and Isomer 2 (CysI-CysVI, CysII-CysIV, CysIII-CysV), but not the native μ-conotoxin arrangement. To date, there has been no study on the structure and activity of KIIIA comprising the native μ-conotoxin disulfide bond arrangement. Here, we evaluated the synthesis, potency, sodium channel subtype selectivity, and 3D structure of the three isomers of KIIIA. Using a regioselective disulfide bond-forming strategy, we synthetically produced the three μ-conotoxin KIIIA isomers displaying distinct bioactivity and Na_V subtype selectivity across human Na_V channel subtypes 1.2, 1.4, and 1.7. We show that Isomer 1 inhibits Na_V subtypes with a rank order of potency of Na_V1.4 > 1.2 > 1.7 and Isomer 2 in the order of Na_V1.4≈1.2 > 1.7, while the native isomer inhibited Na_V1.4 > 1.7≈1.2. The three KIIIA isomers were further evaluated by NMR solution structure analysis and molecular docking with hNa_V1.2. Our study highlights the importance of investigating alternate disulfide isomers, as disulfide connectivity affects not only the overall structure of the peptides but also the potency and subtype selectivity of μ-conotoxins targeting therapeutically relevant Na_V subtypes.

Toxic Peptide From Palythoa caribaeorum Acting on the TRPV1 Channel Prevents Pentylenetetrazol-Induced Epilepsy in Zebrafish Larvae

PcActx peptide, identified from the transcriptome of zoantharian Palythoa caribaeorum, was clustered into the phylogeny of analgesic polypeptides from sea anemone Heteractis crispa (known as APHC peptides). APHC peptides were considered as inhibitors of transient receptor potential cation channel subfamily V member 1 (TRPV1). TRPV1 is a calcium-permeable channel expressed in epileptic brain areas, serving as a potential target for preventing epileptic seizures. Through in silico and in vitro analysis, PcActx peptide was shown to be a potential TRPV1 channel blocker. In vivo studies showed that the linear and oxidized PcActx peptides caused concentration-dependent increases in mortality of zebrafish larvae. However, monotreatment with PcActx peptides below the maximum tolerated doses (MTD) did not affect locomotor behavior. Moreover, PcActx peptides (both linear and oxidized forms) could effectively reverse pentylenetetrazol (PTZ)-induced seizure-related behavior in zebrafish larvae and prevent overexpression of c-fos and npas4a at the mRNA level. The excessive production of ROS induced by PTZ was markedly attenuated by both linear and oxidized PcActx peptides. It was also verified that the oxidized PcActx peptide was more effective than the linear one. In particular, oxidized PcActx peptide notably modulated the mRNA expression of genes involved in calcium signaling and γ-aminobutyric acid (GABA)ergic-glutamatergic signaling, including calb1, calb2, gabra1, grm1, gria1b, grin2b, gat1, slc1a2b, gad1b, and glsa. Taken together, PcActx peptide, as a novel neuroactive peptide, exhibits prominent anti-epileptic activity, probably through modulating calcium signaling and GABAergic-glutamatergic signaling, and is a promising candidate for epilepsy management.

Rigidity of loop 1 contributes to equipotency of globular and ribbon isomers of α-conotoxin AusIA

α-Conotoxins are small disulfide-rich peptides targeting nicotinic acetylcholine receptors (nAChRs) characterised by a C^IC^II-X_m-C^III-X_n-C^IV framework that invariably adopt the native globular conformations which is typically most potent. α-Conotoxins are divided into several structural subgroups based on the number of residues within the two loops braced by the disulfide bonds (m/n), with the 4/7 and 4/3 subgroups dominating. AusIA is a relatively rare α5/5-conotoxin isolated from the venom of Conus australis. Surprisingly, the ribbon isomer displayed equipotency to the wild-type globular AusIA at human α7-containing nAChR. To understand the molecular basis for equipotency, we determined the co-crystal structures of both isomers at Lymnea stagnalis acetylcholine binding protein. The additional residue in the first loop of AusIA was found to be a critical determinant of equipotency, with 11-fold and 86-fold shifts in potency in favour of globular AusIA over ribbon AusIA observed following deletion of Ala4 or Arg5, respectively. This divergence in the potency between globular AusIA and ribbon AusIA was further enhanced upon truncation of the non-conserved Val at the C-termini. Conversely, equipotency could be replicated in LsIA and TxIA [A10L] following insertion of an Ala in the first loop. These findings provide a new understanding of the role the first loop in ribbon and globular α-conotoxins can play in directing α-conotoxin nAChR pharmacology.

Conformational dynamics in peptide toxins: Implications for receptor interactions and molecular design

Peptide toxins are potent and often exquisitely selective probes of the structure and function of ion channels and receptors, and are therefore of significant interest to the pharmaceutical and biotech industries as both pharmacological tools and therapeutic leads. The three-dimensional structures of peptide toxins are essential as a basis for understanding their structure-activity relationships and their binding to target receptors, as well as in guiding the design of analogues with modified potency and/or selectivity for key targets. NMR spectroscopy has played a key role in elucidating the structures of peptide toxins and probing their structure-function relationships. In this article, we highlight the additional important contribution of NMR to characterising the dynamics of peptide toxins. We also compare the information available from NMR measurements with that afforded by molecular dynamics simulations. We describe several examples of the importance of dynamics measurements over a range of timescales for understanding the structure-function relationships of peptide toxins and their receptor engagement. Peptide toxins that inhibit the voltage-gated potassium channel K_V1.3 with pM affinities display different degrees of conformational flexibility, even though they contain multiple disulfide bonds, and this flexibility can affect the relative orientation of residues that have been shown to be critical for channel binding. Information on the dynamic properties of peptide toxins is important in the design of analogues or mimetics where receptor-bound structures are not available.

Globular and ribbon isomers of Conus geographus α-conotoxins antagonize human nicotinic acetylcholine receptors

The short disulfide-rich α-conotoxins derived from the venom of Conus snails comprise a conserved C^IC^II(m)C^III(n)C^IV cysteine framework (m and n, number of amino acids) and the majority antagonize nicotinic acetylcholine receptors (nAChRs). Depending on disulfide connectivity, α-conotoxins can exist as either globular (C^I-C^III, C^II-C^IV), ribbon (C^I-C^IV, C^II-C^III) or bead (C^I-C^II, C^III-C^IV) isomers. In the present study, C. geographus α-conotoxins GI, GIB, G1.5 and G1.9 were chemically synthesized as globular and ribbon isomers and their activity investigated at human nAChRs expressed in Xenopus oocytes using the two-electrode voltage clamp recording technique. Both the globular and ribbon isomers of the 3/5 (m/n) α-conotoxins GI and GIB selectively inhibit heterologous human muscle-type α1β1δε nAChRs, whereas G1.5, a 4/7 α-conotoxin, selectively antagonizes neuronal (non-muscle) nAChR subtypes particularly human α3β2, α7 and α9α10 nAChRs. In contrast, globular and ribbon isomers of G1.9, a novel C-terminal elongated 4/8 α-conotoxin exhibited no activity at the human nAChR subtypes studied. This study reinforces earlier observations that 3/5 α-conotoxins selectively target the muscle nAChR subtypes, although interestingly, GIB is also active at α7 and α9 α10 nAChRs. The 4/7 α-conotoxins target human neuronal nAChR subtypes whereas the pharmacology of the 4/8 α-conotoxin remains unknown.

Linking molecular evolution to molecular grafting

Molecular grafting is a strategy for the engineering of molecular scaffolds into new functional agents, such as next-generation therapeutics. Despite its wide use, studies so far have focused almost exclusively on demonstrating its utility rather than understanding the factors that lead to either poor or successful grafting outcomes. Here, we examine protein evolution and identify parallels between the natural process of protein functional diversification and the artificial process of molecular grafting. We discuss features of natural proteins that are correlated to innovability-the capacity to acquire new functions-and describe their implications to molecular grafting scaffolds. Disulfide-rich peptides are used as exemplars because they are particularly promising scaffolds onto which new functions can be grafted. This article provides a perspective on why some scaffolds are more suitable for grafting than others, identifying opportunities on how molecular grafting might be improved.

LC-Trapped Ion Mobility Spectrometry-TOF MS Differentiation of 2- and 3-Disulfide-Bonded Isomers of the μ-Conotoxin PIIIA

Disulfide bonds within cysteine-rich peptides are important for their stability and biological function. In this respect, the correct disulfide connectivity plays a decisive role. The differentiation of individual disulfide-bonded isomers by traditional high-performance liquid chromatography (HPLC) and mass spectrometry (MS) is limited due to the similarity in physicochemical properties of the isomers sharing the same amino acid sequence. By using trapped ion mobility spectrometry-mass spectrometry (TIMS-MS), several 2- and 3-disulfide-bonded isomers of the μ-conotoxin PIIIA were investigated for their distinguishability by collision cross section (CCS) values and their characteristic mobilogram traces. The isomers could be differentiated by TIMS-MS and also identified in mixing experiments. Thus, TIMS-MS provides a highly valuable and enriching addition to standard HPLC and MS analysis of conformational isomers of disulfide-rich peptides and proteins.

α-Conotoxin Peptidomimetics: Probing the Minimal Binding Motif for Effective Analgesia

Several analgesic α-conotoxins have been isolated from marine cone snails. Structural modification of native peptides has provided potent and selective analogues for two of its known biological targets-nicotinic acetylcholine and γ-aminobutyric acid (GABA) G protein-coupled (GABAB) receptors. Both of these molecular targets are implicated in pain pathways. Despite their small size, an incomplete understanding of the structure-activity relationship of α-conotoxins at each of these targets has hampered the development of therapeutic leads. This review scrutinises the N-terminal domain of the α-conotoxin family of peptides, a region defined by an invariant disulfide bridge, a turn-inducing proline residue and multiple polar sidechain residues, and focusses on structural features that provide analgesia through inhibition of high-voltage-activated Ca2+ channels. Elucidating the bioactive conformation of this region of these peptides may hold the key to discovering potent drugs for the unmet management of debilitating chronic pain associated with a wide range of medical conditions.

Graph-Directed Approach for Downselecting Toxins for Experimental Structure Determination

Conotoxins are short, cysteine-rich peptides of great interest as novel therapeutic leads and of great concern as lethal biological agents due to their high affinity and specificity for various receptors involved in neuromuscular transmission. Currently, of the approximately 6000 known conotoxin sequences, only about 3% have associated structural characterization, which leads to a bottleneck in rapid high-throughput screening (HTS) for identification of potential leads or threats. In this work, we combine a graph-based approach with homology modeling to expand the library of conotoxin structures and to identify those conotoxin sequences that are of the greatest value for experimental structural characterization. The latter would allow for the rapid expansion of the known structural space for generating high quality template-based models. Our approach generalizes to other evolutionarily-related, short, cysteine-rich venoms of interest. Overall, we present and validate an approach for venom structure modeling and experimental guidance and employ it to produce a 290%-larger library of approximate conotoxin structures for HTS. We also provide a set of ranked conotoxin sequences for experimental structure determination to further expand this library.

Discovery of Cyclotides from Australasian Plants

This article is part of a special issue celebrating the contributions of Professor Paul Alewood to peptide science. We begin by providing a summary of collaborative projects between the Alewood and Craik groups at The University of Queensland and highlighting the impacts of some of these studies. In particular, studies on the discovery, synthesis, structures, and bioactivities of disulfide-rich toxins from animal venoms have led to a greater understanding of the biology of ion channels and to applications of these bioactive peptides in drug design. The second part of the article focuses on plant-derived disulfide-rich cyclic peptides, known as cyclotides, and includes an analysis of the geographical distribution of Australasian plant species that contain cyclotides as well as an analysis of the diversity of cyclotide sequences found in Australasian plants. This should provide a useful resource for researchers to access native cyclotides and explore their chemistry and biology.

Ribbon α-Conotoxin KTM Exhibits Potent Inhibition of Nicotinic Acetylcholine Receptors

KTM is a 16 amino acid peptide with the sequence WCCSYPGCYWSSSKWC. Here, we present the nuclear magnetic resonance (NMR) structure and bioactivity of this rationally designed α-conotoxin (α-CTx) that demonstrates potent inhibition of rat α3β2-nicotinic acetylcholine receptors (rα3β2-nAChRs). Two bioassays were used to test the efficacy of KTM. First, a qualitative PC12 cell-based assay confirmed that KTM acts as a nAChR antagonist. Second, bioactivity evaluation by two-electrode voltage clamp electrophysiology was used to measure the inhibition of rα3β2-nAChRs by KTM (IC50 = 0.19 ± 0.02 nM), and inhibition of the same nAChR isoform by α-CTx MII (IC50 = 0.35 ± 0.8 nM). The three-dimensional structure of KTM was determined by NMR spectroscopy, and the final set of 20 structures derived from 32 distance restraints, four dihedral angle constraints, and two disulfide bond constraints overlapped with a mean global backbone root-mean-square deviation (RMSD) of 1.7 ± 0.5 Å. The structure of KTM did not adopt the disulfide fold of α-CTx MII for which it was designed, but instead adopted a flexible ribbon backbone and disulfide connectivity of C2-C16 and C3-C8 with an estimated 12.5% α-helical content. In contrast, α-CTx MII, which has a native fold of C2-C8 and C3-C16, has an estimated 38.1% α-helical secondary structure. KTM is the first reported instance of a Framework I (CC-C-C) α-CTx with ribbon connectivity to display sub-nanomolar inhibitory potency of rα3β2-nAChR subtypes.

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.

α-Conotoxin VnIB from Conus ventricosus is a potent and selective antagonist of α6β4* nicotinic acetylcholine receptors

α6-containing (α6*) nicotinic acetylcholine receptors (nAChRs) are expressed throughout the periphery and the central nervous system and constitute putative therapeutic targets in pain, addiction and movement disorders. The α6β2* nAChRs are relatively well studied, in part due to the availability of target specific α-conotoxins (α-Ctxs). In contrast, all native α-Ctxs identified that potently block α6β4 nAChRs exhibit higher potencies for the closely related α6β2β3 and/or α3β4 subtypes. In this study, we have identified a novel peptide from Conus ventricosus with pronounced selectivity for the α6β4 nAChR. The peptide-encoding gene was cloned from genomic DNA and the predicted mature peptide, α-Ctx VnIB, was synthesized. The functional properties of VnIB were characterized at rat and human nAChRs expressed in Xenopus oocytes by two-electrode voltage clamp electrophysiology. VnIB potently inhibited ACh-evoked currents at rα6β4 and rα6/α3β4 nAChRs, displayed ∼20-fold and ∼250-fold lower potencies at rα3β4 and rα6/α3β2β3 receptors, respectively, and exhibited negligible effects at eight other nAChR subtypes. Interestingly, even higher degrees of selectivity were observed for hα6/α3β4 over hα6/α3β2β3 and hα3β4 receptors. Finally, VnIB displayed fast binding kinetics at rα6/α3β4 (on-rate t_½ = 0.87 min^-1, off-rate t_½ = 2.7 min^-1). The overall preference of VnIB for β4* over β2* nAChRs is similar to the selectivity profiles of other 4/6 α-Ctxs. However, in contrast to previously identified native α-Ctxs targeting α6* nAChRs, VnIB displays pronounced selectivity for α6β4 nAChRs over both α3β4 and α6β2β3 receptors. VnIB thus represents a novel molecular probe for elucidating the physiological role and therapeutic properties of α6β4* nAChRs.

Periplasmic Expression of 4/7 α-Conotoxin TxIA Analogs in E. coli Favors Ribbon Isomer Formation - Suggestion of a Binding Mode at the α7 nAChR

Peptides derived from animal venoms provide important research tools for biochemical and pharmacological characterization of receptors, ion channels, and transporters. Some venom peptides have been developed into drugs (such as the synthetic ω-conotoxin MVIIA, ziconotide) and several are currently undergoing clinical trials for various clinical indications. Challenges in the development of peptides include their usually limited supply from natural sources, cost-intensive chemical synthesis, and potentially complicated stereoselective disulfide-bond formation in the case of disulfide-rich peptides. In particular, if extended structure–function analysis is performed or incorporation of stable isotopes for NMR studies is required, the comparatively low yields and high costs of synthesized peptides might constitute a limiting factor. Here we investigated the expression of the 4/7 α-conotoxin TxIA, a potent blocker at α3β2 and α7 nicotinic acetylcholine receptors (nAChRs), and three analogs in the form of maltose binding protein fusion proteins in Escherichia coli. Upon purification via nickel affinity chromatography and release of the toxins by protease cleavage, HPLC analysis revealed one major peak with the correct mass for all peptides. The final yield was 1–2 mg of recombinant peptide per liter of bacterial culture. Two-electrode voltage clamp analysis on oocyte-expressed nAChR subtypes demonstrated the functionality of these peptides but also revealed a 30 to 100-fold potency decrease of expressed TxIA compared to chemically synthesized TxIA. NMR spectroscopy analysis of TxIA and two of its analogs confirmed that the decreased activity was due to an alternative disulfide linkage rather than the missing C-terminal amidation, a post-translational modification that is common in α-conotoxins. All peptides preferentially formed in the ribbon conformation rather than the native globular conformation. Interestingly, in the case of the α7 nAChR, but not the α3β2 subtype, the loss of potency could be rescued by an R5D substitution. In conclusion, we demonstrate efficient expression of functional but alternatively folded ribbon TxIA variants in E. coli and provide the first structure–function analysis for a ribbon 4/7-α-conotoxin at α7 and α3β2 nAChRs. Computational analysis based on these data provide evidence for a ribbon α-conotoxin binding mode that might be exploited to design ligands with optimized selectivity.

α-Conotoxin GI triazole-peptidomimetics: potent and stable blockers of a human acetylcholine receptor

A conotoxin peptidomimetic was developed as a potential muscle relaxant that is highly potent and blood plasma stable.

Stoichiometry dependent inhibition of rat α3β4 nicotinic acetylcholine receptor by the ribbon isomer of α-conotoxin AuIB

The ribbon isomer of α-conotoxin AuIB has 10-fold greater potency than the wild-type globular isomer at inhibiting nicotinic acetylcholine receptors (nAChRs) in rat parasympathetic neurons, and unlike its globular isoform, ribbon AuIB only targets a specific stoichiometry of the α3β4 nAChR subtype. Previous electrophysiological recordings of AuIB indicated that ribbon AuIB binds to the α3(+)α3(-) interface within the nAChR extracellular domain, which is displayed by the (α3)₃(β4)₂ stoichiometry but not by (α3)₂(β4)₃. This specificity for a particular stoichiometry is remarkable and suggests that ribbon isoforms of α-conotoxins might have great potential in drug design. In this study, we investigated the binding mode and structure-activity relationships of ribbon AuIB using a combination of molecular modeling and electrophysiology recording to determine the features that underpin its selectivity. An alanine scan showed that positions 4 and 9 of ribbon AuIB are the main determinants of the interaction with (α3)₃(β4)₂ nAChR. Our computational models indicate that the first loop of ribbon AuIB binds in the "aromatic box" of the acetylcholine orthosteric binding site, similar to that of globular AuIB. In contrast, the second loop and the termini of the ribbon isomer have different orientations and interactions in the binding sites to those of the globular isomer. The structure-activity relationships reported herein should be useful to design peptides displaying a ribbon α-conotoxin scaffold for inhibition of nAChR subtypes that have hitherto been difficult to selectively target.

Neuronal Nicotinic Acetylcholine Receptor Modulators from Cone Snails

Marine cone snails are a large family of gastropods that have evolved highly potent venoms for predation and defense. The cone snail venom has exceptional molecular diversity in neuropharmacologically active compounds, targeting a range of receptors, ion channels, and transporters. These conotoxins have helped to dissect the structure and function of many of these therapeutically significant targets in the central and peripheral nervous systems, as well as unravelling the complex cellular mechanisms modulated by these receptors and ion channels. This review provides an overview of α-conotoxins targeting neuronal nicotinic acetylcholine receptors. The structure and activity of both classical and non-classical α-conotoxins are discussed, along with their contributions towards understanding nicotinic acetylcholine receptor (nAChR) structure and function.

Aqueous ionic liquids influence the disulfide bond isoform equilibrium in conotoxin AuIB: a consequence of the Hofmeister effect?

The appearance of several disulfide bond isoforms in multiple cysteine containing venom peptides poses a significant challenge in their synthesis and purification under laboratory conditions. Recent experiments suggest that careful tuning of solvent and temperature conditions can propel the disulfide bond isoform equilibrium in favor of the most potent, native form. Certain aqueous ionic liquids (ILs) have proven significantly useful as solvents for this purpose, while exceptions have also been noted. To elucidate the molecular level origin behind such a preference, we report a detailed explicit solvent replica exchange molecular dynamics study of a conotoxin, AuIB, in pure water and four different aqueous IL solutions (~45-60% v/v). The ILs studied here are comprised of cations like 1-ethyl-3-methyl-imidazolium (Im₂₁⁺) or 1-butyl-3-methyl-imidazolium (Im₄₁⁺) coupled with either acetate (OAc^-) or chloride (Cl^-) as the counter anion. Our simulations unfold interesting features of the conformational spaces sampled by the peptide and its solvation in pure water and aqueous IL solutions. Detailed investigation into populations of the globular disulfide bond isoform of AuIB in aqueous IL solutions reveal distinct trends which might be related to the Hofmeister effect of the cation and anion of the IL and of specific interactions of the aqueous IL solutions with the peptide. In accordance with experimental observations, the aqueous [Im₂₁][OAc] solution is found to promote the highest globular isoform population in AuIB.

α-Conotoxins to explore the molecular, physiological and pathophysiological functions of neuronal nicotinic acetylcholine receptors

The vast diversity of neuronal nicotinic acetylcholine subunits expressed in the central and peripheral nervous systems, as well as in non-neuronal tissues, constitutes a formidable challenge for researchers and clinicians to decipher the role of particular subtypes, including complex subunit associations, in physiological and pathophysiological functions. Many natural products target the nAChRs, but there is no richer source of nicotinic ligands than the venom of predatory gastropods known as cone snails. Indeed, every single species of cone snail was shown to produce at least one type of such α-conotoxins. These tiny peptides (10-25 amino acids), constrained by disulfide bridges, proved to be unvaluable tools to investigate the structure and function of nAChRs, some of them having also therapeutic potential. In this review, we provide a recent update on the pharmacology and subtype specificity of several major α-conotoxins.

Expanding the pharmacological profile of κ-hefutoxin 1 and analogues: A focus on the inhibitory effect on the oncogenic channel Kv10.1

Peptide toxins, such as scorpion peptides, are interesting lead compounds in the search for novel drugs. In this paper, the focus is on the scorpion peptide κ-hefutoxin 1. This peptide displays a cysteine-stabilized helix-loop-helix fold (CSα/α) and is known to be a weak K_v1.x inhibitor. Due to the low affinity of κ-hefutoxin 1 for these channels, it is assumed that the main target(s) of κ-hefutoxin 1 remain(s) unknown. In order to identify novel targets, electrophysiological measurements and antifungal assays were performed. The effect of κ-hefutoxin 1 was previously evaluated on a panel of 11 different voltage-gated potassium channels. Here, we extended this target screening with the oncogenic potassium channel K_v10.1. κ-Hefutoxin 1 was able to inhibit this channel in a dose-dependent manner (IC₅₀∼26μM). Although the affinity is rather low, this is the first peptide toxin ever described to be a K_v10.1 inhibitor. The structure-activity relationship of κ-hefutoxin 1 on K_v10.1 was investigated by testing eight κ-hefutoxin 1 variants using the two-electrode voltage clamp technique. Several important amino acid residues were identified; the functional dyad residues (Tyr⁵ and Lys¹⁹), N-terminal residues (Gly¹ and His²) and the amidated C-terminal residue (Cys²²). Since the CSα/α fold is also found in a class of antifungal plant peptides, the α-hairpinines, we investigated the antifungal activity of κ-hefutoxin 1. κ-Hefutoxin 1 showed low activity against the plant pathogen Fusarium culmorum and no activity against three other yeast and fungal species, even at high concentrations (∼100μM).

αO-Conotoxin GeXIVA disulfide bond isomers exhibit differential sensitivity for various nicotinic acetylcholine receptors but retain potency and selectivity for the human α9α10 subtype

Nicotinic acetylcholine receptor (nAChR) subtypes exhibit distinct neuropharmacological properties that are involved in a range of neuropathological conditions, including pain, addiction, epilepsy, autism, schizophrenia, Tourette's syndrome, Alzheimer's and Parkinson's diseases, as well as many types of cancer. The α9α10 nAChR is a potential target in chronic pain, wound healing, the pathophysiology of the auditory system, and breast and lung cancers. αO-conotoxin GeXIVA is a potent antagonist of rat α9α10 nAChRs, with the 'bead' disulfide bond isomer displaying the lowest IC50 of the three possible isomers. In the rat chronic constriction injury model of neuropathic pain, this isomer reduced mechanical hyperalgesia without affecting motor performance. Here, we report the effects of the three disulfide bond isomers of GeXIVA on human α9α10 nAChRs, other human nAChR subtypes, various rat nAChR subtypes, and 10 rat α9α10 nAChR mutants. The three isomers displayed only ∼5-fold difference in potency on the human vs rat α9α10 receptors and had similar affinities at wild-type rat α9α10 nAChRs and all 10 α9α10 receptor mutants. From these findings, the binding site and mechanism of action of GeXIVA on rat and human α9α10 nAChR was deduced to be different from that of other conotoxins targeting this nAChR subtype. GeXIVA is therefore a unique ligand that might prove useful for further probing of binding sites on the α9α10 nAChR. This article is part of the Special Issue entitled 'Venom-derived Peptides as Pharmacological Tools.'

Role of CysI-CysIII Disulfide Bond on the Structure and Activity of α-Conotoxins at Human Neuronal Nicotinic Acetylcholine Receptors

α-Conotoxins preferentially antagonize muscle and neuronal nicotinic acetylcholine receptors (nAChRs). Native α-conotoxins have two disulfide links, CI-CIII and CII-CIV, and owing to the inherent properties of disulfide bonds, α-conotoxins have been systematically engineered to improve their chemical and biological properties. In this study, we explored the possibility of simplifying the disulfide framework of α-conotoxins Vc1.1, BuIA, ImI, and AuIB, by introducing [C2H,C8F] modification to the CI-CIII bond. We therefore explored the possibility of using hydrophobic packing of standard amino acid side chains to replace disulfide bonds as an alternative strategy to nonnatural amino acid cross-links. The impact of CI-CIII disulfide bond replacement on the conformation of the α-conotoxins was investigated using molecular dynamics (MD) simulations and nuclear magnetic resonance chemical shift index study. Two-electrode voltage clamp techniques and MD simulations were used to study the impact of disulfide bond deletion on the activities of the peptides at human neuronal nAChRs. All disulfide-deleted variants except ImI[C2H,C8F] had reduced potency for inhibiting nAChRs. Our results suggest that the CI-CIII disulfide bond is important to stabilize the secondary structure of α-conotoxins as well as their interaction with neuronal nAChR targets. Results from this study enrich our understanding of the function of the CI-CIII disulfide bond and are useful in guiding future structural engineering of the α-conotoxins.

Nicotinic acetylcholine receptor inhibitors derived from snake and snail venoms

The nicotinic acetylcholine receptor (nAChR) represents the prototype of ligand-gated ion channels. It is vital for neuromuscular transmission and an important regulator of neurotransmission. A variety of toxic compounds derived from diverse species target this receptor and have been of elemental importance in basic and applied research. They enabled milestone discoveries in pharmacology and biochemistry ranging from the original formulation of the receptor concept, the first isolation and structural analysis of a receptor protein (the nAChR) to the identification, localization, and differentiation of its diverse subtypes and their validation as a target for therapeutic intervention. Among the venom-derived compounds, α-neurotoxins and α-conotoxins provide the largest families and still represent indispensable pharmacological tools. Application of modified α-neurotoxins provided substantial structural and functional details of the nAChR long before high resolution structures were available. α-bungarotoxin represents not only a standard pharmacological tool and label in nAChR research but also for unrelated proteins tagged with a minimal α-bungarotoxin binding motif. A major advantage of α-conotoxins is their smaller size, as well as superior selectivity for diverse nAChR subtypes that allows their development into ligands with optimized pharmacological and chemical properties and potentially novel drugs. In the following, these two groups of nAChR antagonists will be described focusing on their respective roles in the structural and functional characterization of nAChRs and their development into research tools. In addition, we provide a comparative overview of the diverse α-conotoxin selectivities that can serve as a practical guide for both structure activity studies and subtype classification. This article is part of the Special Issue entitled 'Venom-derived Peptides as Pharmacological Tools.'

Temperature-dependent molecular dynamics study reveals an ionic liquid induced 310 - to α-helical switch in a neurotoxin

Thermal melting and recooling of AuIB, a neurotoxic conopeptide and a highly potent nonaddictive pain reliever is investigated thoroughly in water and an ionic liquid (IL) 1-butyl-3-methylimidazolium Chloride, [Im₄₁ ][Cl] by classical molecular dynamics simulations. Structural evolution of AuIB in water and the IL is observed at different temperatures between 305 and 400 K, to explore how highly viscous ionic solvents affect the peptide structure as compared to conventional solvent water. At 305 K, unlike water, the coercive effect of IL frustrates AuIB secondary structural motifs significantly. As the temperature is raised, a very interesting IL induced conformational transition from 3₁₀ - to α-helix is noticed in the peptide, presumably triggered by a significant restructuring of the peptide H-bond network. The backbone length distributions of the peptide indicate that the IL induced conformational switching is accompanied by a reduction of the axial rise of the helical region, encompassing the residues Pro-6 to Ala-10. Further, we estimated the void space available to the peptide for its structural relaxation within the first solvation shell of ∼5 Å in water as well as in IL. A temperature increase by 100 K, opens up an estimated void volume of ∼70 ų , equivalent to the volume of approximately six water molecules, around the peptide in IL. Cooling simulations of AuIB point to the crucial interplay between thermodynamically favored AuIB conformers and their kinetic control. This study provides a comprehensive understanding of the ionic solvation of biomolecules reinforcing previous experimental findings.

Clustering of disulfide-rich peptides provides scaffolds for hit discovery by phage display: application to interleukin-23

BACKGROUND

Disulfide-rich peptides (DRPs) are found throughout nature. They are suitable scaffolds for drug development due to their small cores, whose disulfide bonds impart extraordinary chemical and biological stability. A challenge in developing a DRP therapeutic is to engineer binding to a specific target. This challenge can be overcome by (i) sampling the large sequence space of a given scaffold through a phage display library and by (ii) panning multiple libraries encoding structurally distinct scaffolds. Here, we implement a protocol for defining these diverse scaffolds, based on clustering structurally defined DRPs according to their conformational similarity.

RESULTS

We developed and applied a hierarchical clustering protocol based on DRP structural similarity, followed by two post-processing steps, to classify 806 unique DRP structures into 81 clusters. The 20 most populated clusters comprised 85% of all DRPs. Representative scaffolds were selected from each of these clusters; the representatives were structurally distinct from one another, but similar to other DRPs in their respective clusters. To demonstrate the utility of the clusters, phage libraries were constructed for three of the representative scaffolds and panned against interleukin-23. One library produced a peptide that bound to this target with an IC₅₀ of 3.3 μM.

CONCLUSIONS

Most DRP clusters contained members that were diverse in sequence, host organism, and interacting proteins, indicating that cluster members were functionally diverse despite having similar structure. Only 20 peptide scaffolds accounted for most of the natural DRP structural diversity, providing suitable starting points for seeding phage display experiments. Through selection of the scaffold surface to vary in phage display, libraries can be designed that present sequence diversity in architecturally distinct, biologically relevant combinations of secondary structures. We supported this hypothesis with a proof-of-concept experiment in which three phage libraries were constructed and panned against the IL-23 target, resulting in a single-digit μM hit and suggesting that a collection of libraries based on the full set of 20 scaffolds increases the potential to identify efficiently peptide binders to a protein target in a drug discovery program.

Role of individual disulfide bridges in the conformation and activity of spinoxin (α-KTx6.13), a potassium channel toxin from Heterometrus spinifer scorpion venom

Spinoxin (SPX; α-KTx6.13), isolated from venom of the scorpion Heterometrus spinifer, is a K⁺ channel-specific peptide toxin (KTx), which adopts a cysteine-stabilized α/β scaffold that is cross-linked by four disulfide bridges (Cys1-Cys5, Cys2-Cys6, Cys3-Cys7, and Cys4-Cys8). To investigate the role of the individual disulfide bonds in the structure-activity relationship of SPX, we synthesized four SPX analogs in which each pair of cysteine residues was replaced by alanine residues. The analysis of circular dichroism spectra and inhibitory activity against Kv1.3 channels showed that the SPX analogs lacking any of three specific disulfide bonds (Cys1-Cys5, Cys2-Cys6, and Cys3-Cys7) were unable to form the native secondary structure and completely lost inhibitory activities. Thus, we conclude that Cys1-Cys5, Cys2-Cys6, and Cys3-Cys7 are required for the inhibition of the Kv1.3 channel by SPX. In contrast, the analog lacking Cys4-Cys8 retained both native secondary structure and inhibitory activity. Interestingly, one of the isomers of the analog lacking Cys1-Cys5 also showed inhibitory activities, although its inhibition was ∼18-fold weaker than native SPX. This isomer had an atypical disulfide bond pairing (Cys3-Cys4 and Cys7-Cys8) that corresponds to that of maurotoxin (MTX), another α-KTx6 family member. These results indicate that the Cys1-Cys5 and Cys2-Cys6 bonds are important for restricting the toxin from forming an atypical (MTX-type) disulfide bond pairing among the remaining four cysteine residues (Cys3, Cys4, Cys7, and Cys8) in native SPX.

Im10A, a short conopeptide isolated from Conus imperialis and possesses two highly concentrated disulfide bridges and analgesic activity

In the present study, we isolated, synthesized and NMR structurally characterized a novel conopeptide Im10A consisting of 11 amino acids (NTICCEGCMCY-NH2) from Conus imperialis. Unlike other conopeptides with four cysteine residues, Im10A had only two residues in loop 1 and one residue in loop 2 (CC-loop1-C-loop2-C), which formed a stable disulfide connectivity "I-IV, II- III" (framework X) with a type I β-turn. Interestingly, Im10A exhibited 50.7% analgesic activity on rat partial sciatic nerve ligation (PNL) at 2h after Im10A administration. However, 10μM Im10A exhibited no apparent effect on neuronal nicotinic acetylcholine receptor, and it did not target DRG voltage-dependent sodium, potassium and calcium ion channels and opioid receptor. To our knowledge, Im10A had the most concentrated disulfide bridges among conopeptides with four cysteine residues. This finding provided a new motif for the future development of biomimetic compounds.

The role of defensive ecological interactions in the evolution of conotoxins

Venoms comprise of complex mixtures of peptides evolved for predation and defensive purposes. Remarkably, some carnivorous cone snails can inject two distinct venoms in response to predatory or defensive stimuli, providing a unique opportunity to study separately how different ecological pressures contribute to toxin diversification. Here, we report the extraordinary defensive strategy of the Rhizoconus subgenus of cone snails. The defensive venom from this worm-hunting subgenus is unusually simple, almost exclusively composed of αD-conotoxins instead of the ubiquitous αA-conotoxins found in the more complex defensive venom of mollusc- and fish-hunting cone snails. A similarly compartmentalized venom gland as those observed in the other dietary groups facilitates the deployment of this defensive venom. Transcriptomic analysis of a Conus vexillum venom gland revealed the αD-conotoxins as the major transcripts, with lower amounts of 15 known and four new conotoxin superfamilies also detected with likely roles in prey capture. Our phylogenetic and molecular evolution analysis of the αD-conotoxins from five subgenera of cone snails suggests they evolved episodically as part of a defensive strategy in the Rhizoconus subgenus. Thus, our results demonstrate an important role for defence in the evolution of conotoxins.

Cloning, synthesis, and characterization of αO-conotoxin GeXIVA, a potent α9α10 nicotinic acetylcholine receptor antagonist

We identified a previously unidentified conotoxin gene from Conus generalis whose precursor signal sequence has high similarity to the O1-gene conotoxin superfamily. The predicted mature peptide, αO-conotoxin GeXIVA (GeXIVA), has four Cys residues, and its three disulfide isomers were synthesized. Previously pharmacologically characterized O1-superfamily peptides, exemplified by the US Food and Drug Administration-approved pain medication, ziconotide, contain six Cys residues and are calcium, sodium, or potassium channel antagonists. However, GeXIVA did not inhibit calcium channels but antagonized nicotinic AChRs (nAChRs), most potently on the α9α10 nAChR subtype (IC50 = 4.6 nM). Toxin blockade was voltage-dependent, and kinetic analysis of toxin dissociation indicated that the binding site of GeXIVA does not overlap with the binding site of the competitive antagonist α-conotoxin RgIA. Surprisingly, the most active disulfide isomer of GeXIVA is the bead isomer, comprising, according to NMR analysis, two well-resolved but uncoupled disulfide-restrained loops. The ribbon isomer is almost as potent but has a more rigid structure built around a short 310-helix. In contrast to most α-conotoxins, the globular isomer is the least potent and has a flexible, multiconformational nature. GeXIVA reduced mechanical hyperalgesia in the rat chronic constriction injury model of neuropathic pain but had no effect on motor performance, warranting its further investigation as a possible therapeutic agent.

Structure and function of μ-conotoxins, peptide-based sodium channel blockers with analgesic activity

μ-Conotoxins block voltage-gated sodium channels (VGSCs) and compete with tetrodotoxin for binding to the sodium conductance pore. Early efforts identified µ-conotoxins that preferentially blocked the skeletal muscle subtype (NaV1.4). However, the last decade witnessed a significant increase in the number of µ-conotoxins and the range of VGSC subtypes inhibited (NaV1.2, NaV1.3 or NaV1.7). Twenty µ-conotoxin sequences have been identified to date and structure-activity relationship studies of several of these identified key residues responsible for interactions with VGSC subtypes. Efforts to engineer-in subtype specificity are driven by in vivo analgesic and neuromuscular blocking activities. This review summarizes structural and pharmacological studies of µ-conotoxins, which show promise for development of selective blockers of NaV1.2, and perhaps also NaV1.1,1.3 or 1.7.

Therapeutic conotoxins: a US patent literature survey

INTRODUCTION

Conotoxins are a large family of bioactive peptides derived from cone snail venom. They target specific classes of ion channels and other membrane proteins and may have therapeutic value, primarily in the management of pain.

AREAS COVERED

The authors surveyed the US patent literature covering conotoxins, and their potential therapeutic applications. They describe the various subclasses of conotoxins that are the subject of current patent applications and their therapeutic indications. Limitations that may preclude broader application of these molecules are discussed and strategies for overcoming these limitations are presented.

EXPERT OPINION

Despite more than 25 years of intense global conotoxin research, only one molecule has successfully reached the market. Several other conotoxin-derived candidates failed in clinical trials, indicating that 'from the bench into the clinic' translation has been more difficult than originally anticipated. Nevertheless, we are optimistic that the potent activities of these molecules and the potential for improving their biopharmaceutical properties may lead to next-generation drug candidates with favorable pharmacological properties.

Bioinformatics-Aided Venomics

Venomics is a modern approach that combines transcriptomics and proteomics to explore the toxin content of venoms. This review will give an overview of computational approaches that have been created to classify and consolidate venomics data, as well as algorithms that have helped discovery and analysis of toxin nucleic acid and protein sequences, toxin three-dimensional structures and toxin functions. Bioinformatics is used to tackle specific challenges associated with the identification and annotations of toxins. Recognizing toxin transcript sequences among second generation sequencing data cannot rely only on basic sequence similarity because toxins are highly divergent. Mass spectrometry sequencing of mature toxins is challenging because toxins can display a large number of post-translational modifications. Identifying the mature toxin region in toxin precursor sequences requires the prediction of the cleavage sites of proprotein convertases, most of which are unknown or not well characterized. Tracing the evolutionary relationships between toxins should consider specific mechanisms of rapid evolution as well as interactions between predatory animals and prey. Rapidly determining the activity of toxins is the main bottleneck in venomics discovery, but some recent bioinformatics and molecular modeling approaches give hope that accurate predictions of toxin specificity could be made in the near future.