Search

Mollusk-hunting (molluscivorous) cone snails belong to a monophyletic group in Conus, a genus of venomous marine snails. The molluscivorous lineage evolved from ancestral worm-hunting (vermivorous) snails ∼18 Ma. To enable the shift to a molluscivorous lifestyle, molluscivorous cone snails must solve biological problems encountered when hunting other gastropods, namely: (i) preventing prey escape and (ii) overcoming the formidable defense of the prey in the form of the molluscan shell, a problem unique to molluscivorous Conus. Here, we show that χ-conotoxins, peptides exclusively expressed in the venoms of molluscivorous Conus, provide solutions to the above problems. Injecting χ-conotoxins into the gastropod mollusk Aplysia californica results in impaired locomotion and uncoordinated hyperactivity. Impaired locomotion impedes escape, and a hyperactive snail will likely emerge from its shell, negating the protection the shell provides. Thus, χ-conotoxins are an evolutionary innovation that accompanied the emergence of molluscivory in Conus and provide solutions to problems posed by hunting other snails., Competing Interests: Conflict of Interest: The authors declare no conflict of interest in the conduct of this research., (© The Author(s) 2024. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

α-Conotoxins are disulfide-rich peptides obtained from the venom of cone snails, which are considered potential molecular probes and drug leads for nAChR-related disorders. However, low specificity towards different nAChR subtypes restricts the further application of many α-conotoxins. In this work, a series of loop1 amino acid-substituted mutants of α-conotoxin RegIIA were synthesized, whose potency and selectivity were evaluated by an electrophysiological approach. The results showed that loop1 alanine scanning mutants [H5A]RegIIA and [P6A]RegIIA blocked rα7 nAChR with IC 50 s of 446 nM and 459 nM, respectively, while their inhibition against rα3β2 and rα3β4 subtypes was negligible, indicating the importance of the fifth and sixth amino acid residues for RegIIA's potency and selectivity. Then, second-generation mutants were designed and synthesized, among which the analogues [H5V]RegIIA and [H5S]RegIIA showed significantly improved selectivity and comparable potency towards rα7 nAChR compared with the native RegIIA. Overall, these findings provide deep insights into the structure-activity relationship of RegIIA, as well as revealing a unique perspective for the further modification and optimization of α-conotoxins and other active peptides.

The two most active disulfide bond isomers of the analgesic αO-conotoxin GeXIVA, namely GeXIVA[1, 2] and GeXIVA[1, 4], were subjected to Asp-scanning mutagenesis to determine the key amino acid residues for activity at the rat α9α10 nicotinic acetylcholine receptor (nAChR). These studies revealed the key role of arginine residues for the activity of GeXIVA isomers towards the α9α10 nAChR. Based on these results, additional analogues with 2-4 mutations were designed and tested. The analogues [T1A,D14A,V28K]GeXIVA[1, 2] and [D14A,I23A,V28K]GeXIVA[1, 4] were developed and showed sub-nanomolar activity for the α9α10 nAChR with IC 50 values of 0.79 and 0.38 nM. The latter analogue had exceptional selectivity for the α9α10 receptor subtype over other nAChR subtypes and can be considered as a drug candidate for further development. Molecular dynamics of receptor-ligand complexes allowed us to make deductions about the possible causes of increases in the affinity of key GeXIVA[1, 4] mutants for the α9α10 nAChR., Competing Interests: Declaration of competing interest The authors declare that they have no competing interests., (Copyright © 2024 Elsevier B.V. All rights reserved.)

Understanding the determinants of α-conotoxin (α-CTX) selectivity for different nicotinic acetylcholine receptor (nAChR) subtypes is a prerequisite for the design of tool compounds to study nAChRs. However, selectivity optimization of these small, disulfide-rich peptides is difficult not only because of an absence of α-CTX/nAChR co-structures but also because it is challenging to predict how a mutation to an α-CTX will alter its potency and selectivity. As a prototypical system to investigate selectivity, we employed the α-CTX LvIA that is 25-fold selective for the α3β2 nAChR over the related α3β4 nAChR subtype, which is a target for nicotine addiction. Using two-electrode voltage clamp electrophysiology, we identified LvIA[D11R] that is 2-fold selective for the α3β4 nAChR, reversing the subtype preference. This effect is specifically due to the change in charge and not shape of LvIA[D11R], as substitution of D11 with citrulline retains selectivity for the α3β2 nAChR. Furthermore, LvIA[D11K] shows a stronger reversal, with 4-fold selectivity for the α3β4 nAChR. Motivated by these findings, using site-directed mutagenesis, we found that β2[K79A] (I79 on β4), but not β2[K78A] (N78 on β4), largely restores the potency of basic mutants at position 11. Finally, to understand the structural basis of this effect, we used AlphaFold2 to generate models of LvIA in complex with both nAChR subtypes. Both models confirm the plausibility of an electrostatic mechanism to explain the data and also reproduce a broad range of potency and selectivity structure-activity relationships for LvIA mutants, as measured using free energy perturbation simulations. Our work highlights how electrostatic interactions can drive α-CTX selectivity and may serve as a strategy for optimizing the selectivity of LvIA and other α-CTXs.

Venomous marine gastropods of the family Conidae are among the most diversified predators in marine realm-in large due to their complex venoms. Besides being a valuable source of bioactive neuropeptides conotoxins, cone-snails venoms are an excellent model for molecular evolution studies, addressing origin of key innovations. However, these studies are handicapped by scarce current knowledge on the tissues involved in venom production, as it is generally assumed the sole prerogative of the venom gland (VG). The role of other secretory glands that are present in all Conus species (salivary gland, SG) or only in some species (accessory salivary gland, ASG) remains poorly understood. Here, for the first time, we carry out a detailed analysis of the VG, SG, and ASG transcriptomes in the vermivorous Conus virgo. We detect multiple transcripts clusters in both the SG and ASG, whose annotations imply venom-related functions. Despite the subsets of transcripts highly-expressed in the VG, SG, and ASG being very distinct, SG expresses an L-, and ASG-Cerm08-, and MEFRR- superfamily conotoxins, all previously considered specific for VG. We corroborate our results with the analysis of published SG and VG transcriptomes from unrelated fish-hunting C. geographus, and C. striatus, possibly fish-hunting C. rolani, and worm-hunting Conus quercinus. In spite of low expression levels of conotoxins, some other specific clusters of putative venom-related peptides are present and may be highly expressed in the SG of these species. Further functional studies are necessary to determine the role that these peptides play in envenomation. In the meantime, our results show importance of routine multi-tissue sampling both for accurate interpretation of tissue-specific venom composition in cone-snails, and for better understanding origin and evolution of venom peptides genes., (© 2023. The Author(s).)

Background: Conus, a highly diverse species of venomous predators, has attracted significant attention in neuroscience and new drug development due to their rich collection of neuroactive peptides called conotoxins. Recent advancements in transcriptome, proteome, and genome analyses have facilitated the identification of conotoxins within Conus' venom glands, providing insights into the genetic features and evolutionary patterns of conotoxin genes. However, the underlying mechanism behind the extraordinary hypervariability of conotoxins remains largely unknown., Results: We analyzed the transcriptomes of 34 Conus species, examining various tissues such as the venom duct, venom bulb, and salivary gland, leading to the identification of conotoxin genes. Genetic variation analysis revealed that a subset of these genes (15.78% of the total) in Conus species underwent positive selection (Ka/Ks > 1, p < 0.01). Additionally, we reassembled and annotated the genome of C. betulinus, uncovering 221 conotoxin-encoding genes. These genes primarily consisted of three exons, with a significant portion showing high transcriptional activity in the venom ducts. Importantly, the flanking regions and adjacent introns of conotoxin genes exhibited a higher prevalence of transposon elements, suggesting their potential contribution to the extensive variability observed in conotoxins. Furthermore, we detected genome duplication in C. betulinus, which likely contributed to the expansion of conotoxin gene numbers. Interestingly, our study also provided evidence of introgression among Conus species, indicating that interspecies hybridization may have played a role in shaping the evolution of diverse conotoxin genes., Conclusions: This study highlights the impact of adaptive evolution and introgressive hybridization on the genetic diversity of conotoxin genes and the evolution of Conus. We also propose a hypothesis suggesting that transposable elements might significantly contribute to the remarkable diversity observed in conotoxins. These findings not only enhance our understanding of peptide genetic diversity but also present a novel approach for peptide bioengineering., (© 2023. BioMed Central Ltd., part of Springer Nature.)

The marine cone snail produces one of the fastest prey strikes in the animal kingdom. It injects highly efficacious venom, often causing prey paralysis and death within seconds. Each snail has hundreds of conotoxins, which serve as a source for discovering and utilizing novel analgesic peptide therapeutics. In this study, we discovered, isolated, and synthesized a novel α3/5-conotoxins derived from the milked venom of Conus obscurus (α-conotoxin OI) and identified the presence of α-conotoxin SI-like sequence previously found in the venom of Conus striatus . Five synthetic analogs of the native α-conotoxin OI were generated. These analogs incorporated single residue or double residue mutations. Three synthetic post-translational modifications (PTMs) were synthetically incorporated into these analogs: N -terminal truncation, proline hydroxylation, and tryptophan bromination. The native α-conotoxin OI demonstrated nanomolar potency in Poecilia reticulata and Homo sapiens muscle-type nicotinic acetylcholine receptor (nAChR) isoforms. Moreover, the synthetic α-[P9K] conotoxin OI displayed enhanced potency in both bioassays, ranging from a 2.85 (LD 50 ) to 18.4 (IC 50 ) fold increase in comparative bioactivity. The successful incorporation of PTMs, with retention of both potency and nAChR isoform selectivity, ultimately pushes new boundaries of peptide bioengineering and the generation of novel α-conotoxin-like sequences.

Venoms of predatory marine cone snails are intensely studied because of the biomedical applications of the neuropeptides that they contain, termed conotoxins. Meanwhile some gastropod lineages have independently acquired secretory glands strikingly similar to the venom gland of cone snails, suggesting that they possess similar venoms. Here we focus on the most diversified of these clades, the genus Vexillum . Based on the analysis of a multi-species proteo-transcriptomic dataset, we show that Vexillum species indeed produce complex venoms dominated by highly diversified short cysteine-rich peptides, vexitoxins. Vexitoxins possess the same precursor organization, display overlapping cysteine frameworks and share several common post-translational modifications with conotoxins. Some vexitoxins show sequence similarity to conotoxins and adopt similar domain conformations, including a pharmacologically relevant inhibitory cysteine knot motif. The Vexillum envenomation gland (gL) is a notably more recent evolutionary novelty than the conoidean venom gland. Thus, we hypothesize lower divergence between vexitoxin genes, and their ancestral 'somatic' counterparts compared to that in conotoxins, and we find support for this hypothesis in the evolution of the vexitoxin cluster V027. We use this example to discuss how future studies on vexitoxins can inform the origin of conotoxins, and how they may help to address outstanding questions in venom evolution.

Somatostatin and its related peptides (SSRPs) form an important family of hormones with diverse physiological roles. The ubiquitous presence of SSRPs in vertebrates and several invertebrate deuterostomes suggests an ancient origin of the SSRP signaling system. However, the existence of SSRP genes outside of deuterostomes has not been established, and the evolutionary history of this signaling system remains poorly understood. Our recent discovery of SSRP-like toxins (consomatins) in venomous marine cone snails (Conus) suggested the presence of a related signaling system in mollusks and potentially other protostomes. Here, we identify the molluscan SSRP-like signaling gene that gave rise to the consomatin family. Following recruitment into venom, consomatin genes experienced strong positive selection and repeated gene duplications resulting in the formation of a hyperdiverse family of venom peptides. Intriguingly, the largest number of consomatins was found in worm-hunting species (>400 sequences), indicating a homologous system in annelids, another large protostome phylum. Consistent with this, comprehensive sequence mining enabled the identification of SSRP-like sequences (and their corresponding orphan receptor) in annelids and several other protostome phyla. These results established the existence of SSRP-like peptides in many major branches of bilaterians and challenge the prevailing hypothesis that deuterostome SSRPs and protostome allatostatin-C are orthologous peptide families. Finally, having a large set of predator-prey SSRP sequences available, we show that although the cone snail's signaling SSRP-like genes are under purifying selection, the venom consomatin genes experience rapid directional selection to target receptors in a changing mix of prey., (© The Author(s) 2022. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

The defensive use of cone snail venom is hypothesised to have first arisen in ancestral worm-hunting snails and later repurposed in a compartmentalised venom duct to facilitate the dietary shift to molluscivory and piscivory. Consistent with its placement in a basal lineage, we demonstrate that the C. distans venom gland lacked distinct compartmentalisation. Transcriptomics revealed C. distans expressed a wide range of structural classes, with inhibitory cysteine knot (ICK)-containing peptides dominating. To better understand the evolution of the venom gland compartmentalisation, we compared C. distans to C. planorbis , the earliest diverging species from which a defence-evoked venom has been obtained, and fish-hunting C. geographus from the Gastridium subgenus that injects distinct defensive and predatory venoms. These comparisons support the hypothesis that venom gland compartmentalisation arose in worm-hunting species and enabled repurposing of venom peptides to facilitate the dietary shift from vermivory to molluscivory and piscivory in more recently diverged cone snail lineages.

Conotoxins are tools used by marine Conus snails to hunt and are a significant repository for marine drug research. Conotoxins highly selectively coordinate different subtypes of various ion channels, and a few have been used in pain management. Although more than 8000 conotoxin genes have been found, the biological activity and function of most have not yet been examined. In this report, we selected the toxin gene QcMNCL-XIII0.1 from our previous investigation and studied it in vitro. First, we successfully prepared active recombinant QcMNCL-XIII0.1 using a TrxA (Thioredoxin A)-assisted folding expression vector based on genetic engineering technology. Animal experiments showed that the recombinant QcMNCL-XIII0.1 exhibited nerve conduction inhibition similar to that of pethidine hydrochloride. With flow cytometry combined fluorescent probe Fluo-4 AM, we found that 10 ng/μL recombinant QcMNCL-XIII0.1 inhibited the fluorescence intensity by 31.07% in the 293T cell model transfected with Cav3.1, implying an interaction between α1G T-type calcium channel protein and recombinant QcMNCL-XIII0.1. This toxin could be an important drug in biomedical research and medicine for pain control.

Venoms are complex mixtures of proteins that have evolved repeatedly in the animal kingdom. Cone snail venoms represent one of the best studied venom systems. In nature, this venom can be dynamically adjusted depending on its final purpose, whether to deter predators or hunt prey. Here, the transcriptome of the venom gland and the proteomes of the predation-evoked and defensive venoms of the molluscivorous cone snail Cylinder ammiralis were catalogued. A total of 242 venom-related transcripts were annotated. The conotoxin superfamilies presenting more different peptides were O1, O2, T, and M, which also showed high expression levels (except T). The three precursors of the J superfamily were also highly expressed. The predation-evoked and defensive venoms showed a markedly distinct profile. A total of 217 different peptides were identified, with half of them being unique to one venom. A total of 59 peptides ascribed to 23 different protein families were found to be exclusive to the predatory venom, including the cono-insulin, which was, for the first time, identified in an injected venom. A total of 43 peptides from 20 protein families were exclusive to the defensive venom. Finally, comparisons of the relative abundance (in terms of number of peptides) of the different conotoxin precursor superfamilies showed that most of them present similar abundance regardless of the diet.

Nicotinic acetylcholine receptor (nAChR) subtypes are key drug targets, but it is challenging to pharmacologically differentiate between them because of their highly similar sequence identities. Furthermore, α-conotoxins (α-CTXs) are naturally selective and competitive antagonists for nAChRs and hold great potential for treating nAChR disorders. Identifying selectivity-enhancing mutations is the chief aim of most α-CTX mutagenesis studies, although doing so with traditional docking methods is difficult due to the lack of α-CTX/nAChR crystal structures. Here, we use homology modeling to predict the structures of α-CTXs bound to two nearly identical nAChR subtypes, α3β2 and α3β4, and use free-energy perturbation (FEP) to re-predict the relative potency and selectivity of α-CTX mutants at these subtypes. First, we use three available crystal structures of the nAChR homologue, acetylcholine-binding protein (AChBP), and re-predict the relative affinities of twenty point mutations made to the α-CTXs LvIA, LsIA, and GIC, with an overall root mean square error (RMSE) of 1.08 ± 0.15 kcal/mol and an R 2 of 0.62, equivalent to experimental uncertainty. We then use AChBP as a template for α3β2 and α3β4 nAChR homology models bound to the α-CTX LvIA and re-predict the potencies of eleven point mutations at both subtypes, with an overall RMSE of 0.85 ± 0.08 kcal/mol and an R 2 of 0.49. This is significantly better than the widely used molecular mechanics-generalized born/surface area (MM-GB/SA) method, which gives an RMSE of 1.96 ± 0.24 kcal/mol and an R 2 of 0.06 on the same test set. Next, we demonstrate that FEP accurately classifies α3β2 nAChR selective LvIA mutants while MM-GB/SA does not. Finally, we use FEP to perform an exhaustive amino acid mutational scan of LvIA and predict fifty-two mutations of LvIA to have greater than 100X selectivity for the α3β2 nAChR. Our results demonstrate the FEP is well-suited to accurately predict potency- and selectivity-enhancing mutations of α-CTXs for nAChRs and to identify alternative strategies for developing selective α-CTXs.

The venom duct origins of predatory and defensive venoms has not been studied for hook-and-line fish hunting cone snails despite the pharmacological importance of their venoms. To better understand the biochemistry and evolution of injected predatory and defensive venoms, we compared distal, central and proximal venom duct sections across three specimens of C. striatus (Pionoconus) using proteomic and transcriptomic approaches. A total of 370 conotoxin precursors were identified from the whole venom duct transcriptome. Milked defensive venom was enriched with a potent cocktail of proximally expressed inhibitory α-, ω- and μ-conotoxins compared to milked predatory venom. In contrast, excitatory κA-conotoxins dominated both the predatory and defensive venoms despite their distal expression, suggesting this class of conotoxin can be selectively expressed from the same duct segment in response to either a predatory or defensive stimuli. Given the high abundance of κA-conotoxins in the Pionoconus clade, we hypothesise that the κA-conotoxins have evolved through adaptive evolution following their repurposing from ancestral inhibitory A superfamily conotoxins to facilitate the dietary shift to fish hunting and species radiation in this clade.

The tetrapeptides Li504 and Li520, differing in the modification of the 4- trans -hydroxylation of proline, are novel conopeptides derived from the venom duct transcriptome of the marine cone snail Conus lividus . These predicted mature peptides are homologous to the active site motif of oxidoreductases that catalyze the oxidation, reduction, and rearrangement of disulfide bonds in peptides and proteins. The estimated reduction potential of the disulfide of Li504 and Li520 is within the range of disulfide reduction potentials of oxidoreductases, indicating that they may catalyze the oxidative folding of conotoxins. Conformational features of Li504 and Li520 include the trans configuration of the Cys1-Pro2/Hyp2 peptide bond with a type 1 turn that is similar to the active site motif of glutaredoxin that regulates the oxidation of cysteine thiols to disulfides. Li504- and Li520-assisted oxidative folding of α-conotoxin ImI confirms that Li520 improves the yield of the natively folded peptide by concomitantly decreasing the yield of the non-native disulfide isomer and thus acts as a miniature disulfide isomerase. The geometry of the Cys1-Hyp2 peptide bond of Li520 shifts between the trans and cis configurations in the disulfide form and thiol/thiolate form, which regulates the deprotonation of the N-terminal cysteine residue. Hydrogen bonding of the hydroxyl group of 4- trans -hydroxyproline with the interpeptide chain unit in the mixed disulfide form may play a vital role in shifting the geometry of the Cys1-Hyp2 peptide bond from cis to trans configuration. The Li520 conopeptide together with similar peptides derived from other species may constitute a new family of "redox-active" conopeptides that are integral components of the oxidative folding machinery of conotoxins.

Marine cone snails are predatory gastropods characterized by a well-developed venom apparatus and highly evolved hunting strategies that utilize toxins to paralyze prey and defend against predators. The venom of each species of cone snail has a large number of pharmacologically active peptides known as conopeptides or conotoxins that are usually unique in each species. Nevertheless, venoms of only very few species have been characterized so far by transcriptomic approaches. In this study, we used transcriptome sequencing technologies and mass spectrometric methods to describe the diversity of venom components expressed by a worm-hunting species, Conus bayani . A total of 82 conotoxin sequences were retrieved from transcriptomic data that contain 54 validated conotoxin sequences clustered into 21 gene superfamilies including divergent gene family, 17 sequences clustered to 6 different conotoxin classes, and 11 conotoxins classified as unassigned gene family. Seven new conotoxin sequences showed unusual cysteine patterns. We were also able to identify 19 peptide sequences using mass spectrometry that completely overlapped with the conotoxin sequences obtained from transcriptome analysis. Importantly, herein we document the presence of 16 proteins that include five post-translational modifying enzymes obtained from transcriptomic data. Our results revealed diverse and novel conopeptides of an unexplored species that could be used extensively in biomedical research due to their therapeutic potentials.

Conotoxins are disulfide-rich peptides found in the venom of cone snails. Due to their exquisite potency and high selectivity for a wide range of voltage and ligand gated ion channels they are attractive drug leads in neuropharmacology. Recently, cone snails were found to have the capability to rapidly switch between venom types with different proteome profiles in response to predatory or defensive stimuli. A novel conotoxin, GXIA (original name G117), belonging to the I 3 -subfamily was identified as the major component of the predatory venom of piscivorous Conus geographus . Using 2D solution NMR spectroscopy techniques, we resolved the 3D structure for GXIA, the first structure reported for the I 3 -subfamily and framework XI family. The 32 amino acid peptide is comprised of eight cysteine residues with the resultant disulfide connectivity forming an ICK+1 motif. With a triple stranded β-sheet, the GXIA backbone shows striking similarity to several tarantula toxins targeting the voltage sensor of voltage gated potassium and sodium channels. Supported by an amphipathic surface, the structural evidence suggests that GXIA is able to embed in the membrane and bind to the voltage sensor domain of a putative ion channel target.

Conotoxins are a group of cysteine-rich, neurotoxic peptides isolated from the venom of marine cone snails. MfVIA is a member of the μO-conotoxin family, and acts as an inhibitor of subtype 1.8 voltage-gated sodium ion channels (Na V 1.8). The unique selectivity of MfVIA as an inhibitor of Na V 1.8 makes it an ideal peptide for elucidation of the physiological functions of this voltage-gated ion channel. Previous experimental studies of point mutations of MfVIA showed that the double mutant [E5K,E8K] exhibited greater activity at Na V 1.8 relative to the wild-type toxin. The present study employs molecular dynamics (MD) simulations to examine the effects of various mutations at these key residues (E5 and E8) on the structure and dynamics of MfVIA. Five double mutants were studied, in which the positions 5 and 8 residues were mutated to amino acids with a range of different physicochemical properties, namely [E5A,E8A], [E5D,E8D], [E5F,E8F], [E5K,E8K], and [E5R,E8R]. Except for [E5D,E8D], all of the mutants tend to show decreased contacts at the N-terminus owing to the loss of the R1-E5 salt bridge relative to that of the wild-type, which subsequently lead to greater exposure and flexibility of the N-terminus for most of the mutant peptides studied, potentially rendering them more able to interact with other species, including Na V 1.8. Molecular docking studies of the peptides to Na V 1.8 via different binding mechanisms suggest that the [E5R, E8R] mutant may be especially worthy of further investigation owing to its predicted binding mode, which differs markedly from those of the other peptides in this study., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2020 Elsevier Inc. All rights reserved.)

The α9α10 nicotinic acetylcholine receptor (nAChR) has been characterized as an effective anti-pain target that functions through a non-opioid mechanism. However, as a pentameric ion channel comprised of two different subunits, the specific targeting of α9α10 nAChRs has proven challenging. Previously the 13-amino-acid peptide, RgIA, was shown to block α9α10 nAChRs with high potency and specificity. This peptide, characterized from the venom of the carnivorous marine snail, Conus regius, produced analgesia in several rodent models of chronic pain. Despite promising pre-clinical data in behavioral assays, the number of specific α9α10 nAChR antagonists remains small and the physiological mechanisms of analgesia remain cryptic. In this study, we implement amino-acid substitutions to definitively characterize the chemical properties of RgIA that contribute to its activity against α9α10 nAChRs. Using this mutational approach, we determined the vital role of biochemical side-chain properties and amino acids in the second loop that are amenable to substitutions to further engineer next-generation analogs for the blockade of α9α10 nAChRs., (Copyright © 2020 Elsevier Inc. All rights reserved.)

The venom of various Conus species is composed of a rich variety of unique bioactive peptides, commonly referred to as conotoxins (conopeptides). Most conopeptides have specific receptors or ion channels as physiologically relevant targets. In this paper, high-throughput transcriptome sequencing was performed to analyze putative conotoxin transcripts from the venom duct of a vermivorous cone snail species, Conus litteratus native to the South China Sea. A total of 128 putative conotoxins were identified, most of them belonging to 22 known superfamilies, with 43 conotoxins being regarded as belonging to new superfamilies. Notably, the M superfamily was the most abundant in conotoxins among the known superfamilies. A total of 15 known cysteine frameworks were also described. The largest proportion of cysteine frameworks were VI/VII (C-C-CC-C-C), IX (C-C-C-C-C-C) and XIV (C-C-C-C). In addition, five novel cysteine patterns were also discovered. Simple sequence repeat detection results showed that di-nucleotide was the major type of repetition, and the codon usage bias results indicated that the codon usage bias of the conotoxin genes was weak, but the M, O1, O2 superfamilies differed in codon preference. Gene cloning indicated that there was no intron in conotoxins of the B1- or J superfamily, one intron with 1273-1339 bp existed in a mature region of the F superfamily, which is different from the previously reported gene structure of conotoxins from other superfamilies. This study will enhance our understanding of conotoxin diversity, and the new conotoxins discovered in this paper will provide more potential candidates for the development of pharmacological probes and marine peptide drugs.

α7 nicotinic acetylcholine receptors (nAChR) is an important nicotinic acetylcholine receptors subtype and closely associated with cognitive disorders, such as Alzheimer's and schizophrenia disease. The mutant ArIB (V11L, V16A) of α-conotoxin ArIB with 17-amino acid residues specifically targets α7 nAChR with no obvious effect on other nAChR subtypes. In the study, the synthetic gene encoding mature peptide of ArIB and mutant ArIB (V11L, V16A) carried a fusion protein Trx and 6 × His-tag was separately inserted in pET-32a (+) vector and transformed into Escherichia coli strain BL21(DE3) pLysS for expression. The expressions of Trx-ArIB-His 6 and Trx-ArIB (V11L, V16A)-His 6 were soluble in Escherichia coli , which were purified by Ni-NTA affinity chromatography column and cleaved by enterokinase to release rArIB and rArIB (V11L, V16A). Then, rArIB and rArIB (V11L, V16A) were purified by high-performance liquid chromatography (HPLC) and identified by matrix-assisted laser desorption ionization-time of flight mass spectrometry (MALDI-TOF MS). Bioactivity of rArIB and rArIB (V11L, V16A) was assessed by two-electrode voltage-clamp electrophysiology in Xenopus laevis oocytes expressing human nAChR subtypes. The results indicated that the yield of the fusion proteins was approximately 50 mg/L and rArIB (V11L, V16A) antagonized the α7 nAChR subtype selectively with 8-nM IC 50 . In summary, this study provides an efficient method to biosynthesize α-conotoxin ArIB and rArIB (V11L, V16A) in Escherichia coli , which could be economical to obtain massively bioactive disulfide-rich polypeptides at fast speed.

Motor neuron degeneration and spinal cord demyelination are hallmark pathological events in Amyotrophic Lateral Sclerosis (ALS). Endogenous retrovirus-K (ERVK) expression has an established association with ALS neuropathology, with murine modeling pointing to a role for the ERVK envelope ( env ) gene in disease processes. Here, we describe a novel viral protein cryptically encoded within the ERVK env transcript, which resembles two distinct cysteine-rich neurotoxic proteins: conotoxin proteins found in marine snails and the Human Immunodeficiency Virus (HIV) Tat protein. Consistent with Nuclear factor-kappa B (NF-κB)-induced retrotransposon expression, the ERVK conotoxin-like protein (CTXLP) is induced by inflammatory signaling. CTXLP is found in the nucleus, impacting innate immune gene expression and NF-κB p65 activity. Using human autopsy specimens from patients with ALS, we further showcase CTXLP expression in degenerating motor cortex and spinal cord tissues, concomitant with inflammation linked pathways, including enhancement of necroptosis marker mixed lineage kinase domain-like (MLKL) protein and oligodendrocyte maturation/myelination inhibitor Nogo-A. These findings identify CTXLP as a novel ERVK protein product, which may act as an effector in ALS neuropathology.

The transcriptomes of the venom glands of 13 closely related species of vermivorous cones endemic to West Africa from genera Africonus and Varioconus were sequenced and venom repertoires compared within a phylogenetic framework using one Kalloconus species as outgroup. The total number of conotoxin precursors per species varied between 108 and 221. Individuals of the same species shared about one-fourth of the total conotoxin precursors. The number of common sequences was drastically reduced in the pairwise comparisons between closely related species, and the phylogenetical signal was totally eroded at the inter-generic level (no sequence was identified as shared derived), due to the intrinsic high variability of these secreted peptides. A common set of four conotoxin precursor superfamilies (T, O1, O2 and M) was expanded in all studied cone species, and thus, they are considered the basic venom toolkit for hunting and defense in the West African vermivorous cone snails. Maximum-likelihood ancestral character reconstructions inferred shared conotoxin precursors preferentially at internal nodes close to the tips of the phylogeny (between individuals and between closely related species) as well as in the common ancestor of Varioconus . Besides the common toolkit, the two genera showed significantly distinct catalogues of conotoxin precursors in terms of type of superfamilies present and the abundance of members per superfamily, but had similar relative expression levels indicating functional convergence. Differential expression comparisons between vermivorous and piscivorous cones highlighted the importance of the A and S superfamilies for fish hunting and defense.

ConoMode is a database for complex three-dimensional (3D) structures of conopeptides binding with their target proteins. Conopeptides, a large family of peptides from the venom of marine snails of the Conus genus, have exceptionally diverse sequences, and their high specificity to block ion channels makes them crucial as drug leads and tools for physiological studies. ConoMode is a specialized archive for the collection of 3D coordinate data for the conopeptides and their binding target proteins from published literature and the Protein Data Bank. These 3D structures can be determined using experimental methods such as X-ray crystallography and electron microscopy and computational methods including docking, homology modeling and molecular dynamics simulations. The binding modes for the conopeptides determined using computational modeling must be validated based on experimental data. The 3D coordinate data from ConoMode can be searched, visualized, downloaded and uploaded. Currently, ConoMode manages 19 conopeptide sequences (from 10 Conus species), 15 protein sequences and 37 3D structures. ConoMode utilizes a modern technical framework to provide a good user experience on mobile devices with touch interaction features. Furthermore, the database is fully optimized for unstructured data and flexible data models. Database URL: http://conomode.qnlm.ac/conomode/conomode/index., (© The Author(s) 2020. Published by Oxford University Press.)

Profundiconus is the most divergent cone snail genus and its unique phylogenetic position, sister to the rest of the family Conidae, makes it a key taxon for examining venom evolution and diversity. Venom gland and foot transcriptomes of Profundiconus cf. vaubani and Profundiconus neocaledonicus were de novo assembled, annotated, and analyzed for differential expression. One hundred and thirty-seven venom components were identified from P. cf. vaubani and 82 from P. neocaledonicus , with only four shared by both species. The majority of the transcript diversity was composed of putative peptides, including conotoxins, profunditoxins, turripeptides, insulin, and prohormone-4. However, there were also a significant percentage of other putative venom components such as chymotrypsin and L-rhamnose-binding lectin. The large majority of conotoxins appeared to be from new gene superfamilies, three of which are highly different from previously reported venom peptide toxins. Their low conotoxin diversity and the type of insulin found suggested that these species, for which no ecological information are available, have a worm or molluscan diet associated with a narrow dietary breadth. Our results indicate that Profundiconus venom is highly distinct from that of other cone snails, and therefore important for examining venom evolution in the Conidae family.

α-Conotoxins (α-CTxs) are small disulfide-rich peptides from venom of Conus species that target nicotinic acetylcholine receptors (nAChRs). The muscle-type nAChRs have been recognized as a potential target for several diseases, such as myogenic disorders, muscle dystrophies, and myasthenia gravis. EI, an α4/7-CTx, mainly blocks α1β1δε nAChRs and has an extra N-terminal extension of three amino acids. In this study, the alanine scanning (Ala-scan) mutagenesis was applied in order to identify key residues of EI for binding with mouse α1β1δε nAChR. The Ala-substituted analogues were tested for their abilities of modulating muscle and neuronal nAChRs in Xenopus laevis oocytes using two-electrode voltage clamp (TEVC) recordings. Electrophysiological results indicated that the vital residues for functional activity of EI were His-7, Pro-8, Met-12, and Pro-15. These changes exhibited a significant decrease in potency of EI against mouse α1β1δε nAChR. Interestingly, replacing the critical serine (Ser) at position 13 with an alanine (Ala) residue resulted in a 2-fold increase in potency at the α1β1δε nAChR, and showed loss of activity on α3β2 and α3β4 nAChRs. Selectivity and potency of [S13A] EI was improved compared with wild-type EI (WT EI). In addition, the structure-activity relationship (SAR) of EI revealed that the "Arg1-Asn2-Hyp3" residues at the N-terminus conferred potency at the muscle-type nAChRs, and the deletion analogue △1-3 EI caused a total loss of activity at the α1β1δε nAChR. Circular dichroism (CD) spectroscopy studies demonstrated that activity loss of truncated analogue △1-3 EI for α1β1δε nAChR is attributed to disturbance of the secondary structure. In this report, an Ala-scan mutagenesis strategy is presented to identify crucial residues that are significantly affecting potency of E1 for mouse α1β1δε nAChR. It may also be important in remodeling of some novel ligands for inhibiting muscle-type nAChRs.

The transcriptomes of the venom glands of two individuals of the magician's cone, Pionoconus magus , from Okinawa (Japan) were sequenced, assembled, and annotated. In addition, RNA-seq raw reads available at the SRA database from one additional specimen of P. magus from the Philippines were also assembled and annotated. The total numbers of identified conotoxin precursors and hormones per specimen were 118, 112, and 93. The three individuals shared only five identical sequences whereas the two specimens from Okinawa had 30 sequences in common. The total number of distinct conotoxin precursors and hormones for P. magus was 275, and were assigned to 53 conotoxin precursor and hormone superfamilies, two of which were new based on their divergent signal region. The superfamilies that had the highest number of precursors were M (42), O1 (34), T (27), A (18), O2 (17), and F (13), accounting for 55% of the total diversity. The D superfamily, previously thought to be exclusive of vermivorous cones was found in P. magus and contained a highly divergent mature region. Similarly, the A superfamily alpha 4/3 was found in P. magus despite the fact that it was previously postulated to be almost exclusive of the genus Rhombiconus . Differential expression analyses of P. magus compared to Chelyconus ermineus , the only fish-hunting cone from the Atlantic Ocean revealed that M and A2 superfamilies appeared to be more expressed in the former whereas the O2 superfamily was more expressed in the latter.

The a-Conotoxins are peptide toxins that are found in the venom of marine cone snails and they are potent antagonists of various subtypes of nicotinic acetylcholine receptors (nAChRs). Because nAChRs have an important role in regulating transmitter release, cell excitability, and neuronal integration, nAChR dysfunctions have been implicated in a variety of severe pathologies. We describe the isolation and characterization of α-conotoxin MilIA, the first conopeptide from the venom of Conus milneedwardsi. The peptide was characterized by electrophysiological screening against several types of cloned nAChRs that were expressed in Xenopus laevis oocytes. MilIA, which is a member of the α3/5 family, is an antagonist of muscle type nAChRs with a high selectivity for muscle versus neuronal subtype nAChRs. Several analogues were designed and investigated for their activity in order to determine the key epitopes of MilIA. Native MilIA and analogues both showed activity at the fetal muscle type nAChR. Two single mutations (Met9 and Asn10) allowed for MilIA to strongly discriminate between the two types of muscle nAChRs. Moreover, one analogue, MilIA [∆1,M2R, M9G, N10K, H11K], displayed a remarkable enhanced potency when compared to native peptide. The key residues that are responsible for switching between muscle and neuronal nAChRs preference were elucidated. Interestingly, the same analogue showed a preference for α9α10 nAChRs among the neuronal types.

Understanding why some groups of organisms are more diverse than others is a central goal in macroevolution. Evolvability, or the intrinsic capacity of lineages for evolutionary change, is thought to influence disparities in species diversity across taxa. Over macroevolutionary time scales, clades that exhibit high evolvability are expected to have higher speciation rates. Cone snails (family: Conidae, $>$900 spp.) provide a unique opportunity to test this prediction because their toxin genes can be used to characterize differences in evolvability between clades. Cone snails are carnivorous, use prey-specific venom (conotoxins) to capture prey, and the genes that encode venom are known and diversify through gene duplication. Theory predicts that higher gene diversity confers a greater potential to generate novel phenotypes for specialization and adaptation. Therefore, if conotoxin gene diversity gives rise to varying levels of evolvability, conotoxin gene diversity should be coupled with macroevolutionary speciation rates. We applied exon capture techniques to recover phylogenetic markers and conotoxin loci across 314 species, the largest venom discovery effort in a single study. We paired a reconstructed timetree using 12 fossil calibrations with species-specific estimates of conotoxin gene diversity and used trait-dependent diversification methods to test the impact of evolvability on diversification patterns. Surprisingly, we did not detect any signal for the relationship between conotoxin gene diversity and speciation rates, suggesting that venom evolution may not be the rate-limiting factor controlling diversification dynamics in Conidae. Comparative analyses showed some signal for the impact of diet and larval dispersal strategy on diversification patterns, though detection of a signal depended on the dataset and the method. If our results remain true with increased taxonomic sampling in future studies, they suggest that the rapid evolution of conid venom may cause other factors to become more critical to diversification, such as ecological opportunity or traits that promote isolation among lineages., (© The Author(s) 2019. Published by Oxford University Press, on behalf of the Society of Systematic Biologists. All rights reserved. For permissions, please email: journals.permissions@oup.com.)

The marine gastropod Hemifusus tuba is served as a luxury food in Asian countries and used in traditional Chinese medicine to treat lumbago and deafness. The lack of genomic data on H. tuba is a barrier to aquaculture development and functional characteristics of potential bioactive molecules are poorly understood. In the present study, we used high-throughput sequencing technologies to generate the first transcriptomic database of H. tuba . A total of 41 unique conopeptides were retrieved from 44 unigenes, containing 6-cysteine frameworks belonging to four superfamilies. Duplication of mature regions and alternative splicing were also found in some of the conopeptides, and the de novo assembly identified a total of 76,306 transcripts with an average length of 824.6 nt, of which including 75,620 (99.1%) were annotated. In addition, simple sequence repeats (SSRs) detection identified 14,000 unigenes containing 20,735 SSRs, among which, 23 polymorphic SSRs were screened. Thirteen of these markers could be amplified in Hemifusus ternatanus and seven in Rapana venosa . This study provides reports of conopeptide genes in Buccinidae for the first time as well as genomic resources for further drug development, gene discovery and population resource studies of this species.

Venomous marine cone snails produce peptide toxins (conotoxins) that bind ion channels and receptors with high specificity and therefore are important pharmacological tools. Conotoxins contain conserved cysteine residues that form disulfide bonds that stabilize their structures. To gain structural insight into the large, yet poorly characterized conotoxin H-superfamily, we used NMR and CD spectroscopy along with MS-based analyses to investigate H-Vc7.2 from Conus victoriae , a peptide with a VI/VII cysteine framework. This framework has Cys I -Cys IV /Cys II -Cys V /Cys III -Cys VI connectivities, which have invariably been associated with the inhibitor cystine knot (ICK) fold. However, the solution structure of recombinantly expressed and purified H-Vc7.2 revealed that although it displays the expected cysteine connectivities, H-Vc7.2 adopts a different fold consisting of two stacked β-hairpins with opposing β-strands connected by two parallel disulfide bonds, a structure homologous to the N-terminal region of the human granulin protein. Using structural comparisons, we subsequently identified several toxins and nontoxin proteins with this "mini-granulin" fold. These findings raise fundamental questions concerning sequence-structure relationships within peptides and proteins and the key determinants that specify a given fold., (© 2019 Nielsen et al.)

Species interactions are fundamental ecological forces that can have significant impacts on the evolutionary trajectories of species. Nonetheless, the contribution of predator-prey interactions to genetic and phenotypic divergence remains largely unknown. Predatory marine snails of the family Conidae exhibit specializations for different prey items and intraspecific variation in prey utilization patterns at geographic scales. Because cone snails utilize venom to capture prey and venom peptides are direct gene products, it is feasible to examine the evolution of genes associated with changes in resource utilization. Here, we compared feeding ecologies and venom duct transcriptomes of individuals from three populations of Conus miliaris , a species that exhibits geographic variation in prey utilization and dietary breadth, in order to determine the extent to which dietary differences are correlated with differences in venom composition, and if expanded niche breadth is associated with increased variation in venom composition. While populations showed little to no overlap in resource utilization, taxonomic richness of prey was greatest at Easter Island. Changes in dietary breadth were associated with differences in expression patterns and increased genetic differentiation of toxin-related genes. The Easter Island population also exhibited greater diversity of toxin-related transcripts, but did not show increased variance in expression of these transcripts. These results imply that differences in dietary breadth contribute more to the structural and regulatory differentiation of venoms than differences in diet.

α-Conotoxins selectively bind to nicotinic acetylcholine receptors (nAChRs), which are therapeutic targets due to their important role in signaling transmission in excitable cells. A previous experimental study has demonstrated that carboxylation of the C-terminal of α-conotoxin LsIA reduces its potency to inhibit human α7 nAChR relative to naturally amidated LsIA. However, little is known about the contribution of conformational changes in the receptor and interactions, induced by C-terminal amidation/carboxylation of conotoxins, to selective binding to nAChRs, since most conotoxins and some disulfide-rich peptides from other conotoxin subfamilies possess a naturally amidated C-terminal. In this study, we employ homology modeling and molecular dynamics (MD) simulations to propose the determinants for differential interactions between amidated and carboxylated LsIAs with α7 nAChR. Our findings indicate an overall increased number of contacts favored by binding of amidated LsIA versus its carboxylated counterpart. Toxin-receptor pairwise interactions, which may play a role in enhancing the potency of the former, include ARG10-TRP77, LEU141 and CYS17-GLN79 via persistent hydrogen bonds and cation-π interactions, which are weakened in the carboxylated form due to a strong intramolecular salt-bridge formed by ARG10 and carboxylated C-terminus. The binding of amidated LsIA also induces enhanced movements in loop C and the juxtamembrane Cys-loop that are closely associated with receptor function. Additionally, the impacts of binding of LsIA on the overall structure and inter-subunit contacts were examined using inter-residue network analysis, suggesting a clockwise tilting of the α7 C and F loops upon binding to carboxylated LsIA, which is absent for amidated LsIA binding. The predicted molecular mechanism of LsIA binding to the α7 receptor may provide new insights into the important role of the C-terminal in the binding potency of conotoxins at neuronal nAChRs for pharmacological purposes., Competing Interests: The authors declare no conflict of interest.

The venom of each Conus species consists of a diverse array of neurophysiologically active peptides, which are mostly unique to the examined species. In this study, we performed high-throughput transcriptome sequencing to extract and analyze putative conotoxin transcripts from the venom ducts of 3 vermivorous cone snails ( C. caracteristicus , C. generalis , and C. quercinus ), which are resident in offshore waters of the South China Sea. In total, 118, 61, and 48 putative conotoxins (across 22 superfamilies) were identified from the 3 Conus species, respectively; most of them are novel, and some possess new cysteine patterns. Interestingly, a series of 45 unassigned conotoxins presented with a new framework of C-C-C-C-C-C, and their mature regions were sufficiently distinct from any other known conotoxins, most likely representing a new superfamily. O- and M-superfamily conotoxins were the most abundant in transcript number and transcription level, suggesting their critical roles in the venom functions of these vermivorous cone snails. In addition, we identified numerous functional proteins with potential involvement in the biosynthesis, modification, and delivery process of conotoxins, which may shed light on the fundamental mechanisms for the generation of these important conotoxins within the venom duct of cone snails.

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.

α-Conotoxin RgIA is a selective and potent competitive antagonist of rat α9α10 nicotinic acetylcholine receptors (nAChR), but it is much less potent towards human α9α10 nAChR. Furthermore, RgIA is susceptible to proteolytic degradation due to containing four arginine residues. These disadvantages greatly limit its use for clinical applications. The purpose of this research was to identify critical stereocenters of RgIA and discover more stable analogues, enhancing its bioavailability by using the d-amino acid scan method. The activity of each variant was investigated against rat and human α9α10 nAChRs, which were expressed in Xenopus oocytes . Experimental assays showed that 14 out of 15 analogues had a substantial reduction in potency towards rat α9α10 nAChR. Noticeably, analogue 13 retained full biological activity compared with RgIA. Meanwhile, two other analogues, 14 and 15 , of which l-Args were substituted with d-Args, exhibited a significantly increased potency towards human α9α10 nAChR, although these analogues showed decreased activities against rat α9α10 nAChR. Additionally, these three analogues exhibited a high resistance against enzymatic degradation in human serum and simulated intestinal fluid (SIF). Collectively, our findings suggest that a d-amino acid scan is a useful strategy for investigating how the side-chain chirality of amino acids affects the structure and function of peptides and may facilitate the development of more stable analogues to increase therapeutic potential.

Californiconus californicus , previously named Conus californicus , has always been considered a unique species within cone snails, because of its molecular, toxicological and morphological singularities; including the wide range of its diet, since it is capable of preying indifferently on fish, snails, octopus, shrimps, and worms. We report here a new cysteine pattern conotoxin assigned to the O1-superfamily capable of inhibiting the growth of Mycobacterium tuberculosis (Mtb). The conotoxin was tested on a pathogen reference strain (H37Rv) and multidrug-resistant strains, having an inhibition effect on growth with a minimal inhibitory concentration (MIC) range of 3.52⁻0.22 μM, similar concentrations to drugs used in clinics. The peptide was purified from the venom using reverse phase high-performance liquid chromatography (RP-HPLC), a partial sequence was constructed by Edman degradation, completed by RACE and confirmed with venom gland transcriptome. The 32-mer peptide containing eight cysteine residues was named O1&#95;cal29b, according to the current nomenclature for this type of molecule. Moreover, transcriptomic analysis of O-superfamily toxins present in the venom gland of the snail allowed us to assign several signal peptides to O2 and O3 superfamilies not described before in C. californicus , with new conotoxins frameworks.

Nicotinic acetylcholine receptors (nAChRs) are found throughout the mammalian body and have been studied extensively because of their implication in a myriad of diseases. α-Conotoxins (α-CTxs) are peptide neurotoxins found in the venom of marine snails of genus Conus . α-CTxs are potent and selective antagonists for a variety of nAChR isoforms. Over the past 40 years, α-CTxs have proven to be valuable molecular probes capable of differentiating between closely related nAChR subtypes and have contributed greatly to understanding the physiological role of nAChRs in the mammalian nervous system. Here, we review the amino acid composition and structure of several α-CTxs that selectively target nAChR isoforms and explore strategies and outcomes for introducing mutations in native α-CTxs to direct selectivity and enhance binding affinity for specific nAChRs. This review will focus on structure-activity relationship studies involving native α-CTxs that have been rationally mutated and molecular interactions that underlie binding between ligand and nAChR isoform.

The piscivorous cone snail Conus tulipa has evolved a net-hunting strategy, akin to the deadly Conus geographus , and is considered the second most dangerous cone snail to humans. Here, we present the first venomics study of C. tulipa venom using integrated transcriptomic and proteomic approaches. Parallel transcriptomic analysis of two C. tulipa specimens revealed striking differences in conopeptide expression levels (2.5-fold) between individuals, identifying 522 and 328 conotoxin precursors from 18 known gene superfamilies. Despite broad overlap at the superfamily level, only 86 precursors (11%) were common to both specimens. Conantokins (NMDA antagonists) from the superfamily B1 dominated the transcriptome and proteome of C. tulipa venom, along with superfamilies B2, A, O1, O3, con-ikot-ikot and conopressins, plus novel putative conotoxins precursors T1.3, T6.2, T6.3, T6.4 and T8.1. Thus, C. tulipa venom comprised both paralytic (putative ion channel modulating α-, ω-, μ-, δ-) and non-paralytic (conantokins, con-ikot-ikots, conopressins) conotoxins. This venomic study confirms the potential for non-paralytic conotoxins to contribute to the net-hunting strategy of C. tulipa.

The primary objective of this study was to realize the large-scale discovery of conotoxin sequences from different organs (including the venom duct, venom bulb and salivary gland) of the vermivorous Oak cone snail, Conus quercinus . Using high-throughput transcriptome sequencing, we identified 133 putative conotoxins that belong to 34 known superfamilies, of which nine were previously reported while the remaining 124 were novel conotoxins, with 17 in new and unassigned conotoxin groups. A-, O₁-, M-, and I₂- superfamilies were the most abundant, and the cysteine frameworks XIII and VIII were observed for the first time in the A- and I₂-superfamilies. The transcriptome data from the venom duct, venom bulb and salivary gland showed considerable inter-organizational variations. Each organ had many exclusive conotoxins, and only seven of all the inferred mature peptides were common in the three organs. As expected, most of the identified conotoxins were synthesized in the venom duct at relatively high levels; however, a number of conotoxins were also identified in the venom bulb and the salivary gland with very low transcription levels. Therefore, various organs have different conotoxins with high diversity, suggesting greater contributions from several organs to the high-throughput discovery of new conotoxins for future drug development.

Cone snails (genus Conus ) are venomous marine snails that inject prey with a lethal cocktail of conotoxins, small, secreted, and cysteine-rich peptides. Given the diversity and often high affinity for their molecular targets, consisting of ion channels, receptors or transporters, many conotoxins have become invaluable pharmacological probes, drug leads, and therapeutics. Transcriptome sequencing of Conus venom glands followed by de novo assembly and homology-based toxin identification and annotation is currently the state-of-the-art for discovery of new conotoxins. However, homology-based search techniques, by definition, can only detect novel toxins that are homologous to previously reported conotoxins. To overcome these obstacles for discovery, we have created Conus Pipe, a machine learning tool that utilizes prominent chemical characters of conotoxins to predict whether a certain transcript in a Conus transcriptome, which has no otherwise detectable homologs in current reference databases, is a putative conotoxin. By using Conus Pipe on RNASeq data of 10 species, we report 5148 new putative conotoxin transcripts that have no homologues in current reference databases. 896 of these were identified by at least three out of four models used. These data significantly expand current publicly available conotoxin datasets and our approach provides a new computational avenue for the discovery of novel toxin families.

The α3β4 nicotinic acetylcholine receptor (nAChR) is an important target implicated in various disease states. α-Conotoxin TxID (1) is the most potent antagonist of α3β4 nAChR, but it also exhibits inhibition of α6/α3β4 nAChR. The results of alanine scanning of 1 suggested a vital role for Ser9 in the selectivity of the peptide. In this study, Ser9 was substituted with a series of 14 amino acids, including some non-natural amino acids, displaying different physicochemical characteristics to further improve the selectivity of 1 toward α3β4 nAChR. The pharmacological activities of the mutants were evaluated using an electrophysiological approach. The best selectivity was obtained with [S9K]TxID, 12, which inhibited α3β4 nAChR with an IC 50 of 6.9 nM and had no effects on other nAChRs. Molecular modeling suggested a possible explanation for the high selectivity of 12 toward α3β4 nAChR, providing deeper insight into the interaction between α-conotoxins and nAChRs as well as potential treatments for nAChR-related diseases.

Complete mitochondrial genome (mitogenome) sequence of a worm-hunting cone snail, Conus quercinus, was reported in this study. Its mitogenome, the longest one (16,460 bp) among reported Conus specie, is composed of 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, two ribosomal RNA (rRNA) genes and one D-loop region. The mitochondrial gene arrangement is highly-conserved and identical to other reported. However, the D-loop region of C. quercinus is the longest (943 bp) with the higher A+T content (71.3%) and a long AT tandem repeat stretch (68 bp). Subsequent phylogenetic analysis demonstrated that three different dietary types (vermivorous, molluscivorous and piscivorous) of cone snails are clustered separately, suggesting that the phylogenetics of cone snails is related to their dietary types. In conclusion, our current work improves our understanding of the mitogenomic structure and evolutionary status of the vermivorous C. quercinus, which support the putative hypothesis that the Conus ancestor was vermivorous., Competing Interests: The authors have declared that no competing interests exist.

&alpha;-Conotoxin TxID was discovered from Conus textile by gene cloning, which has 4/6 inter-cysteine loop spacing and selectively inhibits &alpha;3&beta;4 nicotinic acetylcholine receptor (nAChR) subtype. However, TxID is susceptible to modification due to it containing a methionine (Met) residue that easily forms methionine sulfoxide (MetO) in oxidative environment. In this study, we investigated how Met-11 and its derivatives affect the activity of TxID using a combination of electrophysiological recordings and molecular modelling. The results showed most TxID analogues had substantially decreased activities on &alpha;3&beta;4 nAChR with more than 10-fold potency loss and 5 of them demonstrated no inhibition on &alpha;3&beta;4 nAChR. However, one mutant, [M11I]TxID, displayed potent inhibition at &alpha;3&beta;4 nAChR with an IC 50 of 69 nM, which only exhibited 3.8-fold less compared with TxID. Molecular dynamics simulations were performed to expound the decrease in the affinity for &alpha;3&beta;4 nAChR. The results indicate replacement of Met with a hydrophobic moderate-sized Ile in TxID is an alternative strategy to reduce the impact of Met oxidation, which may help to redesign conotoxins containing methionine residue.

To expand our capacity to discover venom sequences from the genomes of venomous organisms, we applied targeted sequencing techniques to selectively recover venom gene superfamilies and nontoxin loci from the genomes of 32 cone snail species (family, Conidae), a diverse group of marine gastropods that capture their prey using a cocktail of neurotoxic peptides (conotoxins). We were able to successfully recover conotoxin gene superfamilies across all species with high confidence (> 100× coverage) and used these data to provide new insights into conotoxin evolution. First, we found that conotoxin gene superfamilies are composed of one to six exons and are typically short in length (mean = ∼85 bp). Second, we expanded our understanding of the following genetic features of conotoxin evolution: 1) positive selection, where exons coding the mature toxin region were often three times more divergent than their adjacent noncoding regions, 2) expression regulation, with comparisons to transcriptome data showing that cone snails only express a fraction of the genes available in their genome (24-63%), and 3) extensive gene turnover, where Conidae species varied from 120 to 859 conotoxin gene copies. Finally, using comparative phylogenetic methods, we found that while diet specificity did not predict patterns of conotoxin evolution, dietary breadth was positively correlated with total conotoxin gene diversity. Overall, the targeted sequencing technique demonstrated here has the potential to radically increase the pace at which venom gene families are sequenced and studied, reshaping our ability to understand the impact of genetic changes on ecologically relevant phenotypes and subsequent diversification.

Nicotinic acetylcholine receptors (nAChRs) are ligand-gated ion channels, which are widely distributed in the central and peripheral nervous system. The α6β2* nAChR is an important subtype, which is closely associated with nicotine addiction and movement disorders etc. α-conotoxin TxIB with 16-amino acid residues specifically targets α6β2* nAChR with no obvious effect on other nAChR subtypes. However, chemical synthesis of TxIB is expensive, and the quantity of native TxIB extracted from cone snail is limited. In the present study, we attempted to obtain TxIB using biological method based on the recombinant expression in Escherichia coli (E. coli). The synthetic gene encoding mature peptide of TxIB was inserted in pET-31b(+) vector and transformed into E. coli strain BLR(DE3)pLysS for expression. The recombinant fusion protein KSI-TxIB-His 6 (KSI, ketosteroid isomerase) was expressed successfully as inclusion body in E. coli, which was purified by Ni-NTA affinity chromatography column and cleaved by cyanogen bromide (CNBr) to release recombinant α-conotoxin TxIB (rTxIB). Then, rTxIB was purified by reverse-phase high-performance liquid chromatography (RP-HPLC) and was identified by electrospray ionization mass spectrometry (ESI-MS). Pharmacological activity of rTxIB was assessed by electrophysiological approaches. The results indicated that it preserved about 50% of potency, but, was even more important, had the same selectivity as the natural conotoxin which may provide an alternative method for quantity production of small peptides with low cost on the premise of not changing their potency., (© 2017 John Wiley & Sons Ltd.)

Cone snail venoms have separately evolved for predation and defense. Despite remarkable inter- and intra-species variability, defined sets of synergistic venom peptides (cabals) are considered essential for prey capture by cone snails. To better understand the role of predatory cabals in cone snails, we used a high-throughput proteomic data mining and visualisation approach. Using this approach, the relationship between the predatory venom peptides from nine C. purpurascens was systematically analysed. Surprisingly, potentially synergistic levels of κ-PVIIA and δ-PVIA were only identified in five of nine specimens. In contrast, the remaining four specimens lacked significant levels of these known excitotoxins and instead contained high levels of the muscle nAChR blockers ψ-PIIIE and αA-PIVA. Interestingly, one of nine specimens expressed both cabals, suggesting that these sub-groups might represent inter-breeding sub-species of C. purpurascens. High throughput cluster analysis also revealed these two cabals clustered with distinct groups of venom peptides that are presently uncharacterised. This is the first report showing that the cone snails of the same species can deploy two separate and distinct predatory cabals for prey capture and shows that the cabals deployed by this species can be more complex than presently realized. Our semi-automated proteomic analysis facilitates the deconvolution of complex venoms to identify co-evolved families of peptides and help unravel their evolutionary relationships in complex venoms.

Marine drugs have developed rapidly in recent decades. Cone snails, a group of more than 700 species, have always been one of the focuses for new drug discovery. These venomous snails capture prey using a diverse array of unique bioactive neurotoxins, usually named as conotoxins or conopeptides. These conotoxins have proven to be valuable pharmacological probes and potential drugs due to their high specificity and affinity to ion channels, receptors, and transporters in the nervous systems of target prey and humans. Several research groups, including ours, have examined the venom gland of cone snails using a combination of transcriptomic and proteomic sequencing, and revealed the existence of hundreds of conotoxin transcripts and thousands of conopeptides in each Conus species. Over 2000 nucleotide and 8000 peptide sequences of conotoxins have been published, and the number is still increasing quickly. However, more than 98% of these sequences still lack 3D structural and functional information. With the rapid development of genomics and bioinformatics in recent years, functional predictions and investigations on conotoxins are making great progress in promoting the discovery of novel drugs. For example, ω-MVIIA was approved by the U.S. Food and Drug Administration in 2004 to treat chronic pain, and nine more conotoxins are at various stages of preclinical or clinical evaluation. In short, the genus Conus , the big family of cone snails, has become an important genetic resource for conotoxin identification and drug development.

The genus Conus comprises approximately 700 species of venomous marine cone snails that are highly efficient predators of worms, snails, and fish. In evolutionary terms, cone snails are relatively young with the earliest fossil records occurring in the Lower Eocene, 55 Ma. The rapid radiation of cone snail species has been accompanied by remarkably high rates of toxin diversification. To shed light on the molecular mechanisms that accompany speciation, we investigated the toxin repertoire of two sister species, Conus andremenezi and Conus praecellens, that were until recently considered a single variable species. A total of 196 and 250 toxin sequences were identified in the venom gland transcriptomes of C. andremenezi and C. praecellens belonging to 25 and 29 putative toxin gene superfamilies, respectively. Comparative analysis with closely (Conus tribblei and Conus lenavati) and more distantly related species (Conus geographus) suggests that speciation is associated with significant diversification of individual toxin genes (exogenes) whereas the expression pattern of toxin gene superfamilies within lineages remains largely conserved. Thus, changes within individual toxin sequences can serve as a sensitive indicator for recent speciation whereas changes in the expression pattern of gene superfamilies are likely to reflect more dramatic differences in a species' interaction with its prey, predators, and competitors., (© The Author 2017. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

As a class of peptides with 10 cysteine residues (-C-CC-C-CC-C-C-C-C-), D-superfamily conotoxins (D-conotoxins) can specifically act on nicotinic acetylcholine receptors (nAChRs). According to the conserved signal peptides of D-conotoxins, seven D-conotoxin precursor sequences with a previously undescribed Cys arrangement (-C-C-C-CC-C-C-C-C-C-) were identified by PCR-RACE methodology in the present study. The alignment of sequences revealed that signal peptide regions were same as D-VxXXA from Conus vexillum, and their mature peptides were almost different from the D-conotoxins. Analyses of the evolutionary tree demonstrated that they had low homology to those reported conotoxins with 10 cysteine residues (less than 35%) and lied in a separate branch in the evolutionary tree. Furthermore, a previously undescribed D-superfamily conotoxin Lt28.1 was further expressed in Pichia pastoris and then functionally characterized. The results showed that the recombinant Lt28.1 targeted α9α10 nAChRs but not other nAChRs subtypes. These findings defined a new branch of D-superfamily and expanded our knowledge of targets and potential application of D-conotoxins., (Copyright © 2017 Elsevier Inc. All rights reserved.)