Prey Shifts Drive Venom Evolution in Cone Snails.

Venom systems are complex traits that have independently emerged multiple times in diverse plant and animal phyla. Within each venomous lineage there typically exists interspecific variation in venom composition where several factors have been proposed as drivers of variation, including phylogeny and diet. Understanding these factors is of broad biological interest and has implications for the development of anti-venom therapies and venom-based drug discovery. Because of their high species richness and the presence of several major evolutionary prey shifts, venomous marine cone snails (genus Conus) provide an ideal system to investigate drivers of interspecific venom variation. Here, by analyzing the venom gland expression profiles of ∼3,000 toxin genes from 42 species of cone snail, we elucidate the role of prey-specific selection pressures in shaping venom variation. By analyzing overall venom composition and individual toxin structures, we demonstrate that the shifts from vermivory to piscivory in Conus are complemented by distinct changes in venom composition independent of phylogeny. In vivo injections of venom from piscivorous cone snails in fish further showed a higher potency compared to venom of non-piscivores demonstrating a selective advantage. Together, our findings provide compelling evidence for the role of prey shifts in directing the venom composition of cone snails and expand our understanding of the mechanisms of venom variation and diversification.

Analogs of α-conotoxin PnIC selectively inhibit α7ß2- over α7-only subtype nicotinic acetylcholine receptors via a novel allosteric mechanism.

This study was undertaken to identify and characterize the first ligands capable of selectively identifying nicotinic acetylcholine receptors containing α7 and ß2 subunits (α7ß2-nAChR subtype). Basal forebrain cholinergic neurons express α7ß2-nAChR. Here, they appear to mediate neuronal dysfunction induced by the elevated levels of oligomeric amyloid-ß associated with early Alzheimer's disease. Additional work indicates that α7ß2-nAChR are expressed across several further critically important cholinergic and GABAergic neuronal circuits within the central nervous system. Further studies, however, are significantly hindered by the inability of currently available ligands to distinguish heteromeric α7ß2-nAChR from the closely related and more widespread homomeric α7-only-nAChR subtype. Functional screening using two-electrode voltage-clamp electrophysiology identified a family of α7ß2-nAChR-selective analogs of α-conotoxin PnIC (α-CtxPnIC). A combined electrophysiology, functional kinetics, site-directed mutagenesis, and molecular dynamics approach was used to further characterize the α7ß2-nAChR selectivity and site of action of these α-CtxPnIC analogs. We determined that α7ß2-nAChR selectivity of α-CtxPnIC analogs arises from interactions at a site distinct from the orthosteric agonist-binding site shared between α7ß2- and α7-only-nAChR. As numerous previously identified α-Ctx ligands are competitive antagonists of orthosteric agonist-binding sites, this study profoundly expands the scope of use of α-Ctx ligands (which have already provided important nAChR research and translational breakthroughs). More immediately, analogs of α-CtxPnIC promise to enable, for the first time, both comprehensive mapping of the distribution of α7ß2-nAChR and detailed investigations of their physiological roles.

A Serendipitous Path to Pharmacology.

My path to research in neuropharmacology has been a coalescing of my training as a molecular biologist and my intense interest in an esoteric group of animals, the fish-hunting cone snails. Attempting to bridge these two disparate worlds has led me to an idiosyncratic career as a pharmacologist.

Identification of sodium channel toxins from marine cone snails of the subgenera Textilia and Afonsoconus.

Voltage-gated sodium (NaV) channels are transmembrane proteins that play a critical role in electrical signaling in the nervous system and other excitable tissues. µ-Conotoxins are peptide toxins from the venoms of marine cone snails (genus Conus) that block NaV channels with nanomolar potency. Most species of the subgenera Textilia and Afonsoconus are difficult to acquire; therefore, their venoms have yet to be comprehensively interrogated for µ-conotoxins. The goal of this study was to find new µ-conotoxins from species of the subgenera Textilia and Afonsoconus and investigate their selectivity at human NaV channels. Using RNA-seq of the venom gland of Conus (Textilia) bullatus, we identified 12 µ-conotoxin (or µ-conotoxin-like) sequences. Based on these sequences we designed primers which we used to identify additional µ-conotoxin sequences from DNA extracted from historical specimens of species from Textilia and Afonsoconus. We synthesized six of these µ-conotoxins and tested their activity on human NaV1.1-NaV1.8. Five of the six synthetic peptides were potent blockers of human NaV channels. Of these, two peptides (BuIIIB and BuIIIE) were potent blockers of hNaV1.3. Three of the peptides (BuIIIB, BuIIIE and AdIIIA) had submicromolar activity at hNaV1.7. This study serves as an example of the identification of new peptide toxins from historical DNA and provides new insights into structure-activity relationships of µ-conotoxins with activity at hNaV1.3 and hNaV1.7.

Using Constellation Pharmacology to Characterize a Novel α-Conotoxin from Conus ateralbus.

The venom of cone snails has been proven to be a rich source of bioactive peptides that target a variety of ion channels and receptors. α-Conotoxins (αCtx) interact with nicotinic acetylcholine receptors (nAChRs) and are powerful tools for investigating the structure and function of the various nAChR subtypes. By studying how conotoxins interact with nAChRs, we can improve our understanding of these receptors, leading to new insights into neurological diseases associated with nAChRs. Here, we describe the discovery and characterization of a novel conotoxin from Conus ateralbus, αCtx-AtIA, which has an amino acid sequence homologous to the well-described αCtx-PeIA, but with a different selectivity profile towards nAChRs. We tested the synthetic αCtx-AtIA using the calcium imaging-based Constellation Pharmacology assay on mouse DRG neurons and found that αCtx-AtIA significantly inhibited ACh-induced calcium influx in the presence of an α7 positive allosteric modulator, PNU-120596 (PNU). However, αCtx-AtIA did not display any activity in the absence of PNU. These findings were further validated using two-electrode voltage clamp electrophysiology performed on oocytes overexpressing mouse α3ß4, α6/α3ß4 and α7 nAChRs subtypes. We observed that αCtx-AtIA displayed no or low potency in blocking α3ß4 and α6/α3ß4 receptors, respectively, but improved potency and selectivity to block α7 nAChRs when compared with αCtx-PeIA. Through the synthesis of two additional analogs of αCtx-AtIA and subsequent characterization using Constellation Pharmacology, we were able to identify residue Trp18 as a major contributor to the activity of the peptide.

Reconstructing the Origins of the Somatostatin and Allatostatin-C Signaling Systems Using the Accelerated Evolution of Biodiverse Cone Snail Toxins.

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.

Chronicling changes in the somatosensory neurons after peripheral nerve injury.

Current drug discovery efforts focus on identifying lead compounds acting on a molecular target associated with an established pathological state. Concerted molecular changes that occur in specific cell types during disease progression have generally not been identified. Here, we used constellation pharmacology to investigate rat dorsal root ganglion neurons using two models of peripheral nerve injury: chronic constriction injury (CCI) and spinal nerve ligation (SNL). In these well-established models of neuropathic pain, we show that the onset of chronic pain is accompanied by a dramatic, previously unreported increase in the number of bradykinin-responsive neurons, with larger increases observed after SNL relative to CCI. To define the neurons with altered expression, we charted the temporal course of molecular changes following 1, 3, 6, and 14 d after SNL injury and demonstrated that specific molecular changes have different time courses during the progression to a pain state. In particular, ATP receptors up-regulated on day 1 postinjury, whereas the increase in bradykinin receptors was gradual after day 3 postinjury. We specifically tracked changes in two subsets of neurons: peptidergic and nonpeptidergic nociceptors. Significant increases occurred in ATP responses in nAChR-expressing isolectin B4+ nonpeptidergic neurons 1 d postinjury, whereas peptidergic neurons did not display any significant change. We propose that remodeling of ion channels and receptors occurs in a concerted and cell-specific manner, resulting in the appearance of bradykinin-responsive neuronal subclasses that are relevant to chronic pain.

An integrative approach to the facile functional classification of dorsal root ganglion neuronal subclasses.

Somatosensory neurons have historically been classified by a variety of approaches, including structural, anatomical, and genetic markers; electrophysiological properties; pharmacological sensitivities; and more recently, transcriptional profile differentiation. These methodologies, used separately, have yielded inconsistent classification schemes. Here, we describe phenotypic differences in response to pharmacological agents as measured by changes in cytosolic calcium concentration for the rapid classification of neurons in vitro; further analysis with genetic markers, whole-cell recordings, and single-cell transcriptomics validated these findings in a functional context. Using this general approach, which we refer to as tripartite constellation analysis (TCA), we focused on large-diameter dorsal-root ganglion (L-DRG) neurons with myelinated axons. Divergent responses to the K-channel antagonist, κM-conopeptide RIIIJ (RIIIJ), reliably identified six discrete functional cell classes. In two neuronal subclasses (L1 and L2), block with RIIIJ led to an increase in [Ca] i Simultaneous electrophysiology and calcium imaging showed that the RIIIJ-elicited increase in [Ca] i corresponded to different patterns of action potentials (APs), a train of APs in L1 neurons, and sporadic firing in L2 neurons. Genetically labeled mice established that L1 neurons are proprioceptors. The single-cell transcriptomes of L1 and L2 neurons showed that L2 neurons are Aδ-low-threshold mechanoreceptors. RIIIJ effects were replicated by application of the Kv1.1 selective antagonist, Dendrotoxin-K, in several L-DRG subclasses (L1, L2, L3, and L5), suggesting the presence of functional Kv1.1/Kv1.2 heteromeric channels. Using this approach on other neuronal subclasses should ultimately accelerate the comprehensive classification and characterization of individual somatosensory neuronal subclasses within a mixed population.

Conotoxin κM-RIIIJ, a tool targeting asymmetric heteromeric Kv1 channels.

The vast complexity of native heteromeric K+ channels is largely unexplored. Defining the composition and subunit arrangement of individual subunits in native heteromeric K+ channels and establishing their physiological roles is experimentally challenging. Here we systematically explored this "zone of ignorance" in molecular neuroscience. Venom components, such as peptide toxins, appear to have evolved to modulate physiologically relevant targets by discriminating among closely related native ion channel complexes. We provide proof-of-principle for this assertion by demonstrating that κM-conotoxin RIIIJ (κM-RIIIJ) from Conus radiatus precisely targets "asymmetric" Kv channels composed of three Kv1.2 subunits and one Kv1.1 or Kv1.6 subunit with 100-fold higher apparent affinity compared with homomeric Kv1.2 channels. Our study shows that dorsal root ganglion (DRG) neurons contain at least two major functional Kv1.2 channel complexes: a heteromer, for which κM-RIIIJ has high affinity, and a putative Kv1.2 homomer, toward which κM-RIIIJ is less potent. This conclusion was reached by (i) covalent linkage of members of the mammalian Shaker-related Kv1 family to Kv1.2 and systematic assessment of the potency of κM-RIIIJ block of heteromeric K+ channel-mediated currents in heterologous expression systems; (ii) molecular dynamics simulations of asymmetric Kv1 channels providing insights into the molecular basis of κM-RIIIJ selectivity and potency toward its targets; and (iii) evaluation of calcium responses of a defined population of DRG neurons to κM-RIIIJ. Our study demonstrates that bioactive molecules present in venoms provide essential pharmacological tools that systematically target specific heteromeric Kv channel complexes that operate in native tissues.

Discovery of a Potent Conorfamide from Conus episcopatus Using a Novel Zebrafish Larvae Assay.

Natural products such as conotoxins have tremendous potential as tools for biomedical research and for the treatment of different human diseases. Conotoxins are peptides present in the venoms of predatory cone snails that have a rich diversity of pharmacological functions. One of the major bottlenecks in natural products research is the rapid identification and evaluation of bioactive molecules. To overcome this limitation, we designed a set of light-induced behavioral assays in zebrafish larvae to screen for bioactive conotoxins. We used this screening approach to test several unique conotoxins derived from different cone snail clades and discovered that a conorfamide from Conus episcopatus, CNF-Ep1, had the most dramatic alterations in the locomotor behavior of zebrafish larvae. Interestingly, CNF-Ep1 is also bioactive in several mouse assay systems when tested in vitro and in vivo. Our novel screening platform can thus accelerate the identification of bioactive marine natural products, and the first compound discovered using this assay has intriguing properties that may uncover novel neuronal circuitry.

Integration of the pSLT Plasmid into the Salmonella Chromosome Results in a Temperature-Sensitive Growth Defect Due to Aberrant DNA Replication.

A mutant of Salmonella enterica serovar Typhimurium was isolated that simultaneously affected two metabolic pathways as follows: NAD metabolism and DNA repair. The mutant was isolated as resistant to a nicotinamide analog and as temperature-sensitive for growth on minimal glucose medium. In this mutant, Salmonella's 94-kb virulence plasmid pSLT had recombined into the chromosome upstream of the NAD salvage pathway gene pncA This insertion blocked most transcription of pncA, which reduced uptake of the nicotinamide analog. The pSLT insertion mutant also exhibited phenotypes associated with induction of the SOS DNA repair system, including an increase in filamentous cells, higher exonuclease III and catalase activities, and derepression of SOS gene expression. Genome sequencing revealed increased read coverage extending out from the site of pSLT insertion. The two pSLT replication origins are likely initiating replication of the chromosome near the normal replication terminus. Too much replication initiation at the wrong site is probably causing the observed growth defects. Accordingly, deletion of both pSLT replication origins restored growth at higher temperatures.IMPORTANCE In studies that insert a second replication origin into the chromosome, both origins are typically active at the same time. In contrast, the integrated pSLT plasmid initiated replication in stationary phase after normal chromosomal replication had finished. The gradient in read coverage extending out from a single site could be a simple but powerful tool for studying replication and detecting chromosomal rearrangements. This technique may be of particular value when a genome has been sequenced for the first time to verify correct assembly.

The three-dimensional structure of an H-superfamily conotoxin reveals a granulin fold arising from a common ICK cysteine framework.

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 CysI-CysIV/CysII-CysV/CysIII-CysVI 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.

Boholamide A, an APD-Class, Hypoxia-Selective Cyclodepsipeptide.

Calcium homeostasis is implicated in some cancers, leading to the possibility that selective control of calcium might lead to new cancer drugs. On the basis of this idea, we designed an assay using a glioblastoma cell line and screened a collection of 1000 unique bacterial extracts. Isolation of the active compound from a hit extract led to the identification of boholamide A (1), a 4-amido-2,4-pentadieneoate (APD)-class peptide. Boholamide A (1) applied in the nanomolar range induces an immediate influx of Ca2+ in glioblastoma and neuronal cells. APD-class natural products are hypoxia-selective cytotoxins that primarily target mitochondria. Like other APD-containing compounds, 1 is hypoxia selective. Since APD natural products have received significant interest as potential chemotherapeutic agents, 1 provides a novel APD scaffold for the development of new anticancer compounds.

trans-Anethole of Fennel Oil is a Selective and Nonelectrophilic Agonist of the TRPA1 Ion Channel.

Transient receptor potential (TRP) cation channels are molecular targets of various natural products. TRPA1, a member of TRP channel family, is specifically activated by natural products such as allyl isothiocyanate (mustard oil), cinnamaldehyde (cinnamon), and allicin (garlic). In this study, we demonstrated that TRPA1 is also a target of trans-anethole in fennel oil (FO) and fennel seed extract. Similar to FO, trans-anethole selectively elicited calcium influx in TRPA1-expressing mouse sensory neurons of the dorsal root and trigeminal ganglia. These FO- and anethole-induced calcium responses were blocked by a selective TRPA1 channel antagonist, HC-030031. Moreover, both FO and trans-anethole induced calcium influx and transmembrane currents in HEK293 cells stably overexpressing human TRPA1 channels, but not in regular HEK293 cells. Mutation of the amino acids S873 and T874 binding site of human TRPA1 significantly attenuated channel activation by trans-anethole, whereas pretreating with glutathione, a nucleophile, did not. Conversely, activation of TRPA1 by the electrophile allyl isothiocyanate was abolished by glutathione, but was ostensibly unaffected by mutation of the ST binding site. Finally, it was found that trans-anethole was capable of desensitizing TRPA1, and unlike allyl isothiocyanate, it failed to induce nocifensive behaviors in mice. We conclude that trans-anethole is a selective, nonelectrophilic, and seemingly less-irritating agonist of TRPA1.

Characterization of the First Conotoxin from Conus ateralbus, a Vermivorous Cone Snail from the Cabo Verde Archipelago.

Conus ateralbus is a cone snail endemic to the west side of the island of Sal, in the Cabo Verde Archipelago off West Africa. We describe the isolation and characterization of the first bioactive peptide from the venom of this species. This 30AA venom peptide is named conotoxin AtVIA (δ-conotoxin-like). An excitatory activity was manifested by the peptide on a majority of mouse lumbar dorsal root ganglion neurons. An analog of AtVIA with conservative changes on three amino acid residues at the C-terminal region was synthesized and this analog produced an identical effect on the mouse neurons. AtVIA has homology with δ-conotoxins from other worm-hunters, which include conserved sequence elements that are shared with δ-conotoxins from fish-hunting Conus. In contrast, there is no comparable sequence similarity with δ-conotoxins from the venoms of molluscivorous Conus species. A rationale for the potential presence of δ-conotoxins, that are potent in vertebrate systems in two different lineages of worm-hunting cone snails, is discussed.

Rapid expansion of the protein disulfide isomerase gene family facilitates the folding of venom peptides.

Formation of correct disulfide bonds in the endoplasmic reticulum is a crucial step for folding proteins destined for secretion. Protein disulfide isomerases (PDIs) play a central role in this process. We report a previously unidentified, hypervariable family of PDIs that represents the most diverse gene family of oxidoreductases described in a single genus to date. These enzymes are highly expressed specifically in the venom glands of predatory cone snails, animals that synthesize a remarkably diverse set of cysteine-rich peptide toxins (conotoxins). Enzymes in this PDI family, termed conotoxin-specific PDIs, significantly and differentially accelerate the kinetics of disulfide-bond formation of several conotoxins. Our results are consistent with a unique biological scenario associated with protein folding: The diversification of a family of foldases can be correlated with the rapid evolution of an unprecedented diversity of disulfide-rich structural domains expressed by venomous marine snails in the superfamily Conoidea.

Metabolic model for diversity-generating biosynthesis.

A conventional metabolic pathway leads to a specific product. In stark contrast, there are diversity-generating metabolic pathways that naturally produce different chemicals, sometimes of great diversity. We demonstrate that for one such pathway, tru, each ensuing metabolic step is slower, in parallel with the increasing potential chemical divergence generated as the pathway proceeds. Intermediates are long lived and accumulate progressively, in contrast with conventional metabolic pathways, in which the first step is rate-limiting and metabolic intermediates are short-lived. Understanding these fundamental differences enables several different practical applications, such as combinatorial biosynthesis, some of which we demonstrate here. We propose that these principles may provide a unifying framework underlying diversity-generating metabolism in many different biosynthetic pathways.

Constellation pharmacology: a new paradigm for drug discovery.

Constellation pharmacology is a cell-based high-content phenotypic-screening platform that utilizes subtype-selective pharmacological agents to elucidate the cell-specific combinations (constellations) of key signaling proteins that define specific cell types. Heterogeneous populations of native cells, in which the different individual cell types have been identified and characterized, are the foundation for this screening platform. Constellation pharmacology is useful for screening small molecules or for deconvoluting complex mixtures of biologically active natural products. This platform has been used to purify natural products and discover their molecular mechanisms. In the ongoing development of constellation pharmacology, there is a positive feedback loop between the pharmacological characterization of cell types and screening for new drug candidates. As constellation pharmacology is used to discover compounds with novel targeting-selectivity profiles, those new compounds then further help to elucidate the constellations of specific cell types, thereby increasing the content of this high-content platform.

Insights into the origins of fish hunting in venomous cone snails from studies of Conus tessulatus.

Prey shifts in carnivorous predators are events that can initiate the accelerated generation of new biodiversity. However, it is seldom possible to reconstruct how the change in prey preference occurred. Here we describe an evolutionary "smoking gun" that illuminates the transition from worm hunting to fish hunting among marine cone snails, resulting in the adaptive radiation of fish-hunting lineages comprising ∼100 piscivorous Conus species. This smoking gun is δ-conotoxin TsVIA, a peptide from the venom of Conus tessulatus that delays inactivation of vertebrate voltage-gated sodium channels. C. tessulatus is a species in a worm-hunting clade, which is phylogenetically closely related to the fish-hunting cone snail specialists. The discovery of a δ-conotoxin that potently acts on vertebrate sodium channels in the venom of a worm-hunting cone snail suggests that a closely related ancestral toxin enabled the transition from worm hunting to fish hunting, as δ-conotoxins are highly conserved among fish hunters and critical to their mechanism of prey capture; this peptide, δ-conotoxin TsVIA, has striking sequence similarity to these δ-conotoxins from piscivorous cone snail venoms. Calcium-imaging studies on dissociated dorsal root ganglion (DRG) neurons revealed the peptide's putative molecular target (voltage-gated sodium channels) and mechanism of action (inhibition of channel inactivation). The results were confirmed by electrophysiology. This work demonstrates how elucidating the specific interactions between toxins and receptors from phylogenetically well-defined lineages can uncover molecular mechanisms that underlie significant evolutionary transitions.

Specialized insulin is used for chemical warfare by fish-hunting cone snails.

More than 100 species of venomous cone snails (genus Conus) are highly effective predators of fish. The vast majority of venom components identified and functionally characterized to date are neurotoxins specifically targeted to receptors, ion channels, and transporters in the nervous system of prey, predators, or competitors. Here we describe a venom component targeting energy metabolism, a radically different mechanism. Two fish-hunting cone snails, Conus geographus and Conus tulipa, have evolved specialized insulins that are expressed as major components of their venoms. These insulins are distinctive in having much greater similarity to fish insulins than to the molluscan hormone and are unique in that posttranslational modifications characteristic of conotoxins (hydroxyproline, γ-carboxyglutamate) are present. When injected into fish, the venom insulin elicits hypoglycemic shock, a condition characterized by dangerously low blood glucose. Our evidence suggests that insulin is specifically used as a weapon for prey capture by a subset of fish-hunting cone snails that use a net strategy to capture prey. Insulin appears to be a component of the nirvana cabal, a toxin combination in these venoms that is released into the water to disorient schools of small fish, making them easier to engulf with the snail's distended false mouth, which functions as a net. If an entire school of fish simultaneously experiences hypoglycemic shock, this should directly facilitate capture by the predatory snail.