Mapping the chemical and sequence space of the ShKT superfamily

Diversity analysis of sea anemone peptide toxins in different tissues of Heteractis crispa based on transcriptomics

Peptide toxins found in sea anemones venom have diverse properties that make them important research subjects in the fields of pharmacology, neuroscience and biotechnology. This study used high-throughput sequencing technology to systematically analyze the venom components of the tentacles, column, and mesenterial filaments of sea anemone Heteractis crispa, revealing the diversity and complexity of sea anemone toxins in different tissues. A total of 1049 transcripts were identified and categorized into 60 families, of which 91.0% were proteins and 9.0% were peptides. Of those 1049 transcripts, 416, 291, and 307 putative proteins and peptide precursors were identified from tentacles, column, and mesenterial filaments respectively, while 428 were identified when the datasets were combined. Of these putative toxin sequences, 42 were detected in all three tissues, including 33 proteins and 9 peptides, with the majority of peptides being ShKT domain, β-defensin, and Kunitz-type. In addition, this study applied bioinformatics approaches to predict the family classification, 3D structures, and functional annotation of these representative peptides, as well as the evolutionary relationships between peptides, laying the foundation for the next step of peptide pharmacological activity research.

Structure-function relationships in ShKT domain peptides: ShKT-Ts1 from the sea anemone Telmatactis stephensoni

Diverse structural scaffolds have been described in peptides from sea anemones, with the ShKT domain being a common scaffold first identified in ShK toxin from Stichodactyla helianthus. ShK is a potent blocker of voltage-gated potassium channels (KV 1.x), and an analog, ShK-186 (dalazatide), has completed Phase 1 clinical trials in plaque psoriasis. The ShKT domain has been found in numerous other species, but only a tiny fraction of ShKT domains has been characterized functionally. Despite adopting the canonical ShK fold, some ShKT peptides from sea anemones inhibit KV 1.x, while others do not. Mutagenesis studies have shown that a Lys-Tyr (KY) dyad plays a key role in KV 1.x blockade, although a cationic residue followed by a hydrophobic residue may also suffice. Nevertheless, ShKT peptides displaying an ShK-like fold and containing a KY dyad do not necessarily block potassium channels, so additional criteria are needed to determine whether new ShKT peptides might show activity against potassium channels. In this study, we used a combination of NMR and molecular dynamics (MD) simulations to assess the potential activity of a new ShKT peptide. We determined the structure of ShKT-Ts1, from the sea anemone Telmatactis stephensoni, examined its tissue localization, and investigated its activity against a range of ion channels. As ShKT-Ts1 showed no activity against KV 1.x channels, we used MD simulations to investigate whether solvent exposure of the dyad residues may be informative in rationalizing and potentially predicting the ability of ShKT peptides to block KV 1.x channels. We show that either a buried dyad that does not become exposed during MD simulations, or a partially exposed dyad that becomes buried during MD simulations, correlates with weak or absent activity against KV 1.x channels. Therefore, structure determination coupled with MD simulations, may be used to predict whether new sequences belonging to the ShKT family may act as potassium channel blockers.

Structure of the voltage-gated potassium channel KV1.3: Insights into the inactivated conformation and binding to therapeutic leads

The voltage-gated potassium channel KV1.3 is an important therapeutic target for the treatment of autoimmune and neuroinflammatory diseases. The recent structures of KV1.3, Shaker-IR (wild-type and inactivating W434F mutant) and an inactivating mutant of rat KV1.2-KV2.1 paddle chimera (KVChim-W362F+S367T+V377T) reveal that the transition of voltage-gated potassium channels from the open-conducting conformation into the non-conducting inactivated conformation involves the rupture of a key intra-subunit hydrogen bond that tethers the selectivity filter to the pore helix. Breakage of this bond allows the side chains of residues at the external end of the selectivity filter (Tyr447 and Asp449 in KV1.3) to rotate outwards, dilating the outer pore and disrupting ion permeation. Binding of the peptide dalazatide (ShK-186) and an antibody-ShK fusion to the external vestibule of KV1.3 narrows and stabilizes the selectivity filter in the open-conducting conformation, although K+ efflux is blocked by the peptide occluding the pore through the interaction of ShK-Lys22 with the backbone carbonyl of KV1.3-Tyr447 in the selectivity filter. Electrophysiological studies on ShK and the closely-related peptide HmK show that ShK blocks KV1.3 with significantly higher potency, even though molecular dynamics simulations show that ShK is more flexible than HmK. Binding of the anti-KV1.3 nanobody A0194009G09 to the turret and residues in the external loops of the voltage-sensing domain enhances the dilation of the outer selectivity filter in an exaggerated inactivated conformation. These studies lay the foundation to further define the mechanism of slow inactivation in KV channels and can help guide the development of future KV1.3-targeted immuno-therapeutics.

Genomic, functional and structural analyses elucidate evolutionary innovation within the sea anemone 8 toxin family

BACKGROUND

The ShK toxin from Stichodactyla helianthus has established the therapeutic potential of sea anemone venom peptides, but many lineage-specific toxin families in Actiniarians remain uncharacterised. One such peptide family, sea anemone 8 (SA8), is present in all five sea anemone superfamilies. We explored the genomic arrangement and evolution of the SA8 gene family in Actinia tenebrosa and Telmatactis stephensoni, characterised the expression patterns of SA8 sequences, and examined the structure and function of SA8 from the venom of T. stephensoni.

RESULTS

We identified ten SA8-family genes in two clusters and six SA8-family genes in five clusters for T. stephensoni and A. tenebrosa, respectively. Nine SA8 T. stephensoni genes were found in a single cluster, and an SA8 peptide encoded by an inverted SA8 gene from this cluster was recruited to venom. We show that SA8 genes in both species are expressed in a tissue-specific manner and the inverted SA8 gene has a unique tissue distribution. While the functional activity of the SA8 putative toxin encoded by the inverted gene was inconclusive, its tissue localisation is similar to toxins used for predator deterrence. We demonstrate that, although mature SA8 putative toxins have similar cysteine spacing to ShK, SA8 peptides are distinct from ShK peptides based on structure and disulfide connectivity.

CONCLUSIONS

Our results provide the first demonstration that SA8 is a unique gene family in Actiniarians, evolving through a variety of structural changes including tandem and proximal gene duplication and an inversion event that together allowed SA8 to be recruited into the venom of T. stephensoni.

Toxin-like neuropeptides in the sea anemone Nematostella unravel recruitment from the nervous system to venom

The sea anemone Nematostella vectensis (Anthozoa, Cnidaria) is a powerful model for characterizing the evolution of genes functioning in venom and nervous systems. Although venom has evolved independently numerous times in animals, the evolutionary origin of many toxins remains unknown. In this work, we pinpoint an ancestral gene giving rise to a new toxin and functionally characterize both genes in the same species. Thus, we report a case of protein recruitment from the cnidarian nervous to venom system. The ShK-like1 peptide has a ShKT cysteine motif, is lethal for fish larvae and packaged into nematocysts, the cnidarian venom-producing stinging capsules. Thus, ShK-like1 is a toxic venom component. Its paralog, ShK-like2, is a neuropeptide localized to neurons and is involved in development. Both peptides exhibit similarities in their functional activities: They provoke contraction in Nematostella polyps and are toxic to fish. Because ShK-like2 but not ShK-like1 is conserved throughout sea anemone phylogeny, we conclude that the two paralogs originated due to a Nematostella-specific duplication of a ShK-like2 ancestor, a neuropeptide-encoding gene, followed by diversification and partial functional specialization. ShK-like2 is represented by two gene isoforms controlled by alternative promoters conferring regulatory flexibility throughout development. Additionally, we characterized the expression patterns of four other peptides with structural similarities to studied venom components and revealed their unexpected neuronal localization. Thus, we employed genomics, transcriptomics, and functional approaches to reveal one venom component, five neuropeptides with two different cysteine motifs, and an evolutionary pathway from nervous to venom system in Cnidaria.

Histidine-Rich Defensins from the Solanaceae and Brasicaceae Are Antifungal and Metal Binding Proteins

Plant defensins are best known for their antifungal activity and contribution to the plant immune system. The defining feature of plant defensins is their three-dimensional structure known as the cysteine stabilized alpha-beta motif. This protein fold is remarkably tolerant to sequence variation with only the eight cysteines that contribute to the stabilizing disulfide bonds absolutely conserved across the family. Mature defensins are typically 46–50 amino acids in length and are enriched in lysine and/or arginine residues. Examination of a database of approximately 1200 defensin sequences revealed a subset of defensin sequences that were extended in length and were enriched in histidine residues leading to their classification as histidine-rich defensins (HRDs). Using these initial HRD sequences as a query, a search of the available sequence databases identified over 750 HRDs in solanaceous plants and 20 in brassicas. Histidine residues are known to contribute to metal binding functions in proteins leading to the hypothesis that HRDs would have metal binding properties. A selection of the HRD sequences were recombinantly expressed and purified and their antifungal and metal binding activity was characterized. Of the four HRDs that were successfully expressed all displayed some level of metal binding and two of four had antifungal activity. Structural characterization of the other HRDs identified a novel pattern of disulfide linkages in one of the HRDs that is predicted to also occur in HRDs with similar cysteine spacing. Metal binding by HRDs represents a specialization of the plant defensin fold outside of antifungal activity.

Transcriptomic Analysis of Four Cerianthid (Cnidaria, Ceriantharia) Venoms

Tube anemones, or cerianthids, are a phylogenetically informative group of cnidarians with complex life histories, including a pelagic larval stage and tube-dwelling adult stage, both known to utilize venom in stinging-cell rich tentacles. Cnidarians are an entirely venomous group that utilize their proteinaceous-dominated toxins to capture prey and defend against predators, in addition to several other ecological functions, including intraspecific interactions. At present there are no studies describing the venom for any species within cerianthids. Given their unique development, ecology, and distinct phylogenetic-placement within Cnidaria, our objective is to evaluate the venom-like gene diversity of four species of cerianthids from newly collected transcriptomic data. We identified 525 venom-like genes between all four species. The venom-gene profile for each species was dominated by enzymatic protein and peptide families, which is consistent with previous findings in other cnidarian venoms. However, we found few toxins that are typical of sea anemones and corals, and furthermore, three of the four species express toxin-like genes closely related to potent pore-forming toxins in box jellyfish. Our study is the first to provide a survey of the putative venom composition of cerianthids and contributes to our general understanding of the diversity of cnidarian toxins.

The voltage-gated potassium channel KV1.3 as a therapeutic target for venom-derived peptides

The voltage-gated potassium channel K_V1.3 is a well-established therapeutic target for a range of autoimmune diseases, in addition to being the site of action of many venom-derived peptides. Numerous studies have documented the efficacy of venom peptides that target K_V1.3, in particular from sea anemones and scorpions, in animal models of autoimmune diseases such as rheumatoid arthritis, psoriasis and multiple sclerosis. Moreover, an analogue of the sea anemone peptide ShK (known as dalazatide) has successfully completed Phase 1 clinical trials in mild-to-moderate plaque psoriasis. In this article we consider other potential therapeutic applications of inhibitors of K_V1.3, including in inflammatory bowel disease and neuroinflammatory conditions such as Alzheimer's and Parkinson's diseases, as well as fibrotic diseases. We also summarise strategies for facilitating the entry of peptides to the central nervous system, given that this will be a pre-requisite for the treatment of most neuroinflammatory diseases. Venom-derived peptides that have been reported recently to target K_V1.3 are also described. The increasing number of autoimmune and other conditions in which K_V1.3 is upregulated and is therefore a potential therapeutic target, combined with the fact that many venom-derived peptides are potent inhibitors of K_V1.3, suggests that venoms are likely to continue to serve as a rich source of new pharmacological tools and therapeutic leads targeting this channel.

Evolution of Sequence-Diverse Disordered Regions in a Protein Family: Order within the Chaos

Approaches for studying the evolution of globular proteins are now well established yet are unsuitable for disordered sequences. Our understanding of the evolution of proteins containing disordered regions therefore lags that of globular proteins, limiting our capacity to estimate their evolutionary history, classify paralogs, and identify potential sequence–function relationships. Here, we overcome these limitations by using new analytical approaches that project representations of sequence space to dissect the evolution of proteins with both ordered and disordered regions, and the correlated changes between these. We use the fasciclin-like arabinogalactan proteins (FLAs) as a model family, since they contain a variable number of globular fasciclin domains as well as several distinct types of disordered regions: proline (Pro)-rich arabinogalactan (AG) regions and longer Pro-depleted regions.

Sequence space projections of fasciclin domains from 2019 FLAs from 78 species identified distinct clusters corresponding to different types of fasciclin domains. Clusters can be similarly identified in the seemingly random Pro-rich AG and Pro-depleted disordered regions. Sequence features of the globular and disordered regions clearly correlate with one another, implying coevolution of these distinct regions, as well as with the N-linked and O-linked glycosylation motifs. We reconstruct the overall evolutionary history of the FLAs, annotated with the changing domain architectures, glycosylation motifs, number and length of AG regions, and disordered region sequence features. Mapping these features onto the functionally characterized FLAs therefore enables their sequence–function relationships to be interrogated. These findings will inform research on the abundant disordered regions in protein families from all kingdoms of life.

Tentacle Transcriptomes of the Speckled Anemone (Actiniaria: Actiniidae: Oulactis sp.): Venom-Related Components and Their Domain Structure

Cnidarians are one of the oldest known animal lineages (ca. 700 million years), with a unique envenomation apparatus to deliver a potent mixture of peptides and proteins. Some peptide toxins from cnidarian venom have proven therapeutic potential. Here, we use a transcriptomic/proteomic strategy to identify sequences with similarity to known venom protein families in the tentacles of the endemic Australian 'speckled anemone' (Oulactis sp.). Illumina RNASeq data were assembled de novo. Annotated sequences in the library were verified by cross-referencing individuals' transcriptomes or protein expression evidence from LC-MS/MS data. Sequences include pore-forming toxins, phospholipases, peptidases, neurotoxins (sodium and potassium channel modulators), cysteine-rich secretory proteins and defensins (antimicrobial peptides). Fewer than 4% of the sequences in the library occurred across the three individuals examined, demonstrating high sequence variability of an individual's arsenal. We searched for actinoporins in Oulactis sp. to assess sequence similarity to the only described toxins (OR-A and -G) for this genus and examined the domain architecture of venom-related peptides and proteins. The novel putative actinoporin of Oulactis sp. has a greater similarity to other species in the Actiniidae family than to O. orientalis. Venom-related sequences have an architecture that occurs in single, repeat or multi-domain combinations of venom-related (e.g. ShK-like) and non-venom (e.g. whey acid protein) domains. This study has produced the first transcriptomes for an endemic Australian sea anemone species and the genus Oulactis, while identifying nearly 400 novel venom-related peptides and proteins for future structural and functional analyses and venom evolution studies.