Scorpion venomics: a 2019 overview.

Introduction: A few scorpions are dangerous to humans. Their medical relevance was the initial driving force for venom research. By classical biochemistry and molecular cloning, several venom peptides and their coding transcripts were characterized, mainly those related to toxins. The discovery of other components with novel activities and potential applications has revitalized the interest in the field in the last decade and a half. Nontoxic scorpion species have also attracted major interest.Areas covered: Advances in the identification of scorpion venom components via high-throughput venomics (genomics, transcriptomics and proteomics) up to 2019 are summarized. A classification system for venom-related transcripts and proteins, together with an intuitive systematic nomenclature for RNAseq-generated transcripts are proposed. Venom components classified as Na+, K+, Ca2+, Cl- and TRP channel toxins, enzymes, protease inhibitors, host defense peptides and other peptidic molecules are briefly reviewed, giving a comprehensive picture of the venom.Expert opinion: Modern high-throughput technologies applied to scorpion venom studies have resulted in a dramatic increase in both, the number and diversity of available sequences, leading to a deeper understanding of the composition of scorpion venoms. Still, many newly-discovered venom constituents remain to be characterized, to complete the puzzle of scorpion venoms.

Profiling the resting venom gland of the scorpion Tityus stigmurus through a transcriptomic survey.

BACKGROUND: The scorpion Tityus stigmurus is widely distributed in Northeastern Brazil and known to cause severe human envenoming, inducing pain, hyposthesia, edema, erythema, paresthesia, headaches and vomiting. The present study uses a transcriptomic approach to characterize the gene expression profile from the non-stimulated venom gland of Tityus stigmurus scorpion. RESULTS: A cDNA library was constructed and 540 clones were sequenced and grouped into 153 clusters, with one or more ESTs (expressed sequence tags). Forty-one percent of ESTs belong to recognized toxin-coding sequences, with transcripts encoding antimicrobial toxins (AMP-like) being the most abundant, followed by alfa KTx- like, beta KTx-like, beta NaTx-like and alfa NaTx-like. Our analysis indicated that 34% of the transcripts encode "other possible venom molecules", which correspond to anionic peptides, hypothetical secreted peptides, metalloproteinases, cystein-rich peptides and lectins. Fifteen percent of ESTs are similar to cellular transcripts. Sequences without good matches corresponded to 11%. CONCLUSIONS: This investigation provides the first global view of gene expression of the venom gland from Tityus stigmurus under resting conditions. This approach enables characterization of a large number of venom gland component molecules, which belong either to known or non yet described types of venom peptides and proteins from the Buthidae family.

Solution structure of BmKalphaTx11, a toxin from the venom of the Chinese scorpion Buthus martensii Karsch.

The solution structure of BmKalphaTx11 presented by this paper is distinctive from any other structures of wide-type scorpion alpha-toxins reported so far, for its trans-9,10 peptide bond conformation is accompanied by 'protruding' topology of the 'NC-domain'. The orientation of the C-tail of BmKalphaTx11 is obviously different from that of classical alpha-toxins (e.g., AaH2, BmK-M8), despite the fact that they share common trans conformation of peptide bond between residues 9 and 10. Accordingly, there must be other structural factors dominating the orientation of the C-tail except the conformation of peptide bond 9-10. Our study reveals that residues at position 58 play an important role in it, and different type of residues at this position (e.g., Lys, Arg, Met, Ile) result in different spatial relationship between the C-terminus and the 'five-residue-turn' and then different topology of the 'NC-domain', therefore residues at position 58 are believed to function as structure and bioactivity switch for specificity of scorpion alpha-toxins. The mechanism for stabilizing the geometry of the 'NC-domain' in wide-type scorpion alpha-toxins is also discussed.

Bioinformatic characterizations and prediction of K+ and Na+ ion channels effector toxins.

BACKGROUND: K+ and Na+ channel toxins constitute a large set of polypeptides, which interact with their ion channel targets. These polypeptides are classified in two different structural groups. Recently a new structural group called birtoxin-like appeared to contain both types of toxins has been described. We hypothesized that peptides of this group may contain two conserved structural motifs in K+ and/or Na+ channels scorpion toxins, allowing these birtoxin-like peptides to be active on K+ and/or Na+ channels. RESULTS: Four multilevel motifs, overrepresented and specific to each group of K+ and/or Na+ ion channel toxins have been identified, using GIBBS and MEME and based on a training dataset of 79 sequences judged as representative of K+ and Na+ toxins.Unexpectedly birtoxin-like peptides appeared to present a new structural motif distinct from those present in K+ and Na+ channels Toxins. This result, supported by previous experimental data, suggests that birtoxin-like peptides may exert their activity on different sites than those targeted by classic K+ or Na+ toxins.Searching, the nr database with these newly identified motifs using MAST, retrieved several sequences (116 with e-value < 1) from various scorpion species (test dataset). The filtering process left 30 new and highly likely ion channel effectors.Phylogenetic analysis was used to classify the newly found sequences. Alternatively, classification tree analysis, using CART algorithm adjusted with the training dataset, using the motifs and their 2D structure as explanatory variables, provided a model for prediction of the activity of the new sequences. CONCLUSION: The phylogenetic results were in perfect agreement with those obtained by the CART algorithm.Our results may be used as criteria for a new classification of scorpion toxins based on functional motifs.

Scorpion toxins interact with nicotinic acetylcholine receptors.

Neurotoxins are among the main components of scorpion and snake venoms. Scorpion neurotoxins affect voltage-gated ion channels, while most snake neurotoxins target ligand-gated ion channels, mainly nicotinic acetylcholine receptors (nAChRs). We report that scorpion venoms inhibit α-bungarotoxin binding to both muscle-type nAChR from Torpedo californica and neuronal human α7 nAChR. Toxins inhibiting nAChRs were identified as OSK-1 (α-KTx family) from Orthochirus scrobiculosus and HelaTx1 (κ-KTx family) from Heterometrus laoticus, both being blockers of voltage-gated potassium channels. With an IC50 of 1.6 µm, OSK1 inhibits acetylcholine-induced current through mouse muscle-type nAChR heterologously expressed in Xenopus oocytes. Other well-characterized scorpion toxins from these families also bind to Torpedo nAChR with micromolar affinities. Our results indicate that scorpion neurotoxins present target promiscuity.

Chinese-scorpion (Buthus martensi Karsch) toxin BmK alphaIV, a novel modulator of sodium channels: from genomic organization to functional analysis.

In the present study, BmK alphaIV, a novel modulator of sodium channels, was cloned from venomous glands of the Chinese scorpion (Buthus martensi Karsch) and expressed successfully in Escherichia coli. The BmK alphaIV gene is composed of two exons separated by a 503 bp intron. The mature polypeptide contains 66 amino acids. BmK alphaIV has potent toxicity in mice and cockroaches. Surface-plasmon-resonance analysis found that BmK alphaIV could bind to both rat cerebrocortical synaptosomes and cockroach neuronal membranes, and shared similar binding sites on sodium channels with classical AaH II (alpha-mammal neurotoxin from the scorpion Androctonus australis Hector), BmK AS (beta-like neurotoxin), BmK IT2 (the depressant insect-selective neurotoxin) and BmK abT (transitional neurotoxin), but not with BmK I (alpha-like neurotoxin). Two-electrode voltage clamp recordings on rNav1.2 channels expressed in Xenopus laevis oocytes revealed that BmK alphaIV increased the peak amplitude and prolonged the inactivation phase of Na+ currents. The structural and pharmacological properties compared with those of other scorpion alpha-toxins suggests that BmK alphaIV represents a novel subgroup or functional hybrid of alpha-toxins and might be an evolutionary intermediate neurotoxin for alpha-toxins.

Adaptive evolution of insect selective excitatory ß-type sodium channel neurotoxins from scorpion venom.

Insect selective excitatory ß-type sodium channel neurotoxins from scorpion venom (ß-NaScTxs) are composed of about 70-76 amino acid residues and share a common scaffold stabilized by four unique disulfide bonds. The phylogenetic analysis of these toxins was hindered by limited sequence data. In our recent study, two new insect selective excitatory ß-NaScTxs, LmIT and ImIT, were isolated from Lychas mucronatus and Isometrus maculatus, respectively. With the sequences previously reported, we examined the adaptive molecular evolution of insect selective excitatory ß-NaScTxs by estimating the nonsynonymous-to-synonymous rate ratio (ω=dN/dS). The results revealed 12 positively selected sites in the genes of insect selective excitatory ß-NaScTxs. Moreover, these positively selected sites match well with the sites important for interacting with sodium channels, as demonstrated in previous mutagenesis study. These results reveal that adaptive evolution after gene duplication is one of the most important genetic mechanisms of scorpion neurotoxin diversification.

CysView: protein classification based on cysteine pairing patterns.

CysView is a web-based application tool that identifies and classifies proteins according to their disulfide connectivity patterns. It accepts a dataset of annotated protein sequences in various formats and returns a graphical representation of cysteine pairing patterns. CysView displays cysteine patterns for those records in the data with disulfide annotations. It allows the viewing of records grouped by connectivity patterns. CysView's utility as an analysis tool was demonstrated by the rapid and correct classification of scorpion toxin entries from GenPept on the basis of their disulfide pairing patterns. It has proved useful for rapid detection of irrelevant and partial records, or those with incomplete annotations. CysView can be used to support distant homology between proteins. CysView is publicly available at http://research.i2r.a-star.edu.sg/CysView/.

Categorization of venoms according to bonding properties: An immunological overview.

In this report, we present a study on the antigenic cross-reactivity of various venoms from the most dangerous Egyptian snakes and scorpions belonging to families Elapidae, Viperidae and Buthidae. The study was carried out with special reference to bonding properties between venoms and antivenoms and their involvement in the formation of specific and/or cross-reactive interactions. The homologous polyclonal antivenoms showed high reactivity to the respective venoms and cross-reacted with varying degrees to other non-homologous venoms. Assorting the antivenoms according to their susceptibility to dissociation by different concentrations of NH4SCN revealed that most of the antibodies involved in homologous venom-antivenom interactions were highly avid; building up strong venom-antivenom bonding. Whereas cross-reactions due to heterologous interactions were mediated by less avid antibodies that ultimately led to the formation of venom-antivenom bonding of different power strengths depending on the antigenic similarity and hence on the phylogenetic relationship of the tested venom. A new parameter evaluating high and low avid interactions, designated as H/L value, for each antigen-antibody bonding was initiated and used as an indicator of bonding strength between different venom-antivenom partners. H/L values were many folds higher than 1 for homologous and closely related venoms, 1 or around 1 for cross-reactive venoms, whereas venoms from unrelated remote sources recorded H/L values far less than 1. Using well defined polyclonal antivenoms, H/L value was successfully used to assign eight unknown venoms to their animal families and the results were confirmed by species-specific ELISA and immunoblotting assays.

Structure-function relationships of peptides forming the calcin family of ryanodine receptor ligands.

Calcins are a novel family of scorpion peptides that bind with high affinity to ryanodine receptors (RyRs) and increase their activity by inducing subconductance states. Here, we provide a comprehensive analysis of the structure-function relationships of the eight calcins known to date, based on their primary sequence, three-dimensional modeling, and functional effects on skeletal RyRs (RyR1). Primary sequence alignment and evolutionary analysis show high similarity among all calcins (≥78.8% identity). Other common characteristics include an inhibitor cysteine knot (ICK) motif stabilized by three pairs of disulfide bridges and a dipole moment (DM) formed by positively charged residues clustering on one side of the molecule and neutral and negatively charged residues segregating on the opposite side. [(3)H]Ryanodine binding assays, used as an index of the open probability of RyRs, reveal that all eight calcins activate RyR1 dose-dependently with Kd values spanning approximately three orders of magnitude and in the following rank order: opicalcin1 > opicalcin2 > vejocalcin > hemicalcin > imperacalcin > hadrucalcin > maurocalcin >> urocalcin. All calcins significantly augment the bell-shaped [Ca(2+)]-[(3)H]ryanodine binding curve with variable effects on the affinity constants for Ca(2+) activation and inactivation. In single channel recordings, calcins induce the appearance of a subconductance state in RyR1 that has a unique fractional value (∼20% to ∼60% of the full conductance state) but bears no relationship to binding affinity, DM, or capacity to stimulate Ca(2+) release. Except for urocalcin, all calcins at 100 nM concentration stimulate Ca(2+) release and deplete Ca(2+) load from skeletal sarcoplasmic reticulum. The natural variation within the calcin family of peptides offers a diversified set of high-affinity ligands with the capacity to modulate RyRs with high dynamic range and potency.

Accurate prediction of scorpion toxin functional properties from primary structures.

Scorpion toxins are common experimental tools for studies of biochemical and pharmacological properties of ion channels. The number of functionally annotated scorpion toxins is steadily growing, but the number of identified toxin sequences is increasing at much faster pace. With an estimated 100,000 different variants, bioinformatic analysis of scorpion toxins is becoming a necessary tool for their systematic functional analysis. Here, we report a bioinformatics-driven system involving scorpion toxin structural classification, functional annotation, database technology, sequence comparison, nearest neighbour analysis, and decision rules which produces highly accurate predictions of scorpion toxin functional properties.

Overview of the Knottin scorpion toxin-like peptides in scorpion venoms: Insights on their classification and evolution.

Scorpion venoms include several compounds with different pharmacological activities. Within these compounds, toxins affecting ion channels are among the most studied. They are all peptides that have been classified based on their 3D structure, chain size and function. Usually, they show a spatial arrangement characterized by the presence of a cysteine-stabilized alpha beta motif; most of them affect Na(+) and K(+) ion-channels. These features have been revised in several occasions before, but a complete phylogenetic analysis of the disulfide containing peptides is not been done. In the present contribution, two databases (Pfam and InterPro) including more than 800 toxins from different scorpions were analyzed. Pfam database included toxins from several organisms other than scorpions such as insects and plants, while InterPro included only scorpion toxins. Our results suggest that Na(+) toxins have evolved independently from those of K(+) toxins no matter the length of the peptidic chains. These preliminary results suggest that current classification needs a more detailed revision, in order to have better characterized toxin families, so the new peptides obtained from transcriptomic analyses would be properly classified.

A large number of novel Ergtoxin-like genes and ERG K+-channels blocking peptides from scorpions of the genus Centruroides.

Twenty-three novel sequences similar to Ergtoxin (ErgTx) were obtained by direct sequencing of peptides or deduced from gene cloned using cDNAs of venomous glands of Centruroides (C.) elegans, C. exilicauda, C. gracilis, C. limpidus limpidus, C. noxius and C. sculpturatus. These peptides have from 42 to 47 amino acid residues cross-linked by four disulfide bridges. They share sequence similarities (60-98% compared with ErgTx1) and were shown to block ERG K(+)-channels of F-11 clone (N18TG-2xrat DRG) cultured cells. An unrooted phylogenetic tree analysis of these peptides showed that they conform at least five different subfamilies, of which three are novel subfamilies.

Tc1, from Tityus cambridgei, is the first member of a new subfamily of scorpion toxin that blocks K(+)-channels.

A new peptide, Tc1, containing only 23 amino acids closely packed by three disulfide bridges was isolated from the Amazonian scorpion Tityus cambridgei. It blocks reversibly the Shaker B K(+)-channels with a K(d) of 65 nM and displaces binding of noxiustoxin to mouse brain synaptosome membranes. It is the shortest known peptide from scorpion venom that recognizes K(+)-channels and constitutes a new structural subfamily of toxin, classified as alphaKTx 13.1.

Definition of the alpha-KTx15 subfamily.

Three novel scorpion toxins, Aa1 from Androctonus australis, BmTX3 from Buthus martensi and AmmTX3 from Androctonus mauretanicus were shown able to selectively block A-type K+ currents in cerebellum granular cells or cultured striatum neurons from rat brain. In electrophysiology experiments, the transient A-current completely disappeared when 1 microM of the toxins was applied to the external solution whereas the sustained K+ current was unaffected. The three toxins shared high sequence homologies (more than 94%) and constituted a new 'short-chain' scorpion toxin subfamily: alpha-KTx15. Monoiododerivative of 125I-sBmTX3 specifically bound to rat brain synaptosomes. Under equilibrium binding conditions, maximum binding was 14 fmol/mg of protein and the dissociation constant (Kd) was 0.21 nM. This Kd value was confirmed by kinetic experiments (kon = 6.0 x 10(6) M(-1) s(-1) and koff = 6.0 x 10(-4) s(-1)). Competitions with AmmTX3 and Aa1 with 125I-sBmTX3 bound to its receptor on rat brain synaptosomes showed that they fully inhibited the 125I-sBmTX3 binding (Ki values of 20 and 44 pM, respectively), demonstrating unambiguously that the three molecules shared the same target in rat brain. A panel of toxins described as specific ligands for different K+, Na+ and Ca2+ channels were not able to displace 125I-sBmTX3 from its binding site. Thus, 125I-sBmTX3 is a new ligand for a still unidentified target in rat brain. In autoradiography, the distribution of 125I-sBmTX3 binding sites in the adult rat brain indicated a high density of 125I-sBmTX3 receptors in the striatum, hippocampus, superior colliculus, and cerebellum.

Immunological cross reactivity & paraspecificity of the scorpion Heterometrus bengalensis antivenom.

Immunological cross-reactivity and paraspecificity of scorpion H. bengalensis antivenom were studied to find out the intergeneric therapeutic relationship between the venom of other scorpions in West Bengal Buthus tamulus, Lychas laevifrons and Heterometrus swammerdami. Of these scorpions, Buthus tamulus and Lychas laevifrons failed to show any cross reactivity. However, H. swammerdami venom showed cross-reactivity with H. bengalensis antiserum as revealed from immunogeldiffusion and immunoelectrophoresis. This antiserum protected H. swammerdami venom-induced lethality in mice, blocked the contractile response in smooth muscles and antagonised the venom-induced neuromuscular blockade in rat phrenic nerve diaphragm and chick biventer cervicis.

Molecular cloning and nucleotide sequence analysis of encoded anti-insect toxin BotIT2 from the scorpion Buthus occitanus tunetanus venom.

Numerous toxins from scorpion venoms are much more toxic to insects than to other animal classes, and possess high affinity to Na+ channels. Many of them active on insects were purified from the venom of Buthus occitanus tunetanus. Using amino acid sequences of BotIT2 and RACE-PCR amplification (Rapid amplification of cDNA ends) technique, we isolated, identified and sequenced the nucleotide sequence from the venom glands of the scorpion Buthus occitanus tunetanus. The cDNA encodes a precursor of an insect toxin of 60 amino acid residues. The deduced nucleotide sequence toxin was identical to the determined amino acid sequence of BotIT2. BotIT2 is more similar to the excitatory toxins in its mode of action and to the depressant toxins in its primary structure.

Scorpion venom components that affect ion-channels function.

The number and types of venom components that affect ion-channel function are reviewed. These are the most important venom components responsible for human intoxication, deserving medical attention, often requiring the use of specific anti-venoms. Special emphasis is given to peptides that recognize Na(+)-, K(+)- and Ca(++)-channels of excitable cells. Knowledge generated by direct isolation of peptides from venom and components deduced from cloned genes, whose amino acid sequences are deposited into databanks are nowadays in the order of 1.5 thousands, out of an estimate biodiversity closed to 300,000. Here the diversity of components is briefly reviewed with mention to specific references. Structural characteristic are discussed with examples taken from published work. The principal mechanisms of action of the three different types of peptides are also reviewed. Na(+)-channel specific venom components usually are modifier of the open and closing kinetic mechanisms of the ion-channels, whereas peptides affecting K(+)-channels are normally pore blocking agents. The Ryanodine Ca(++)-channel specific peptides are known for causing sub-conducting stages of the channels conductance and some were shown to be able to internalize penetrating inside the muscle cells.

Characterization of the first K⁺ channel blockers from the venom of the Moroccan scorpion Buthus occitanus Paris.

The availability of a large variety of specific blockers, which inhibit different K(+) currents, would help to elucidate their differences in physiological function. Short peptide toxins isolated from scorpion venoms are able to block voltage-dependent or Ca(2+)-activated K(+) channels. Here, we have studied the venom of the Moroccan scorpion Buthus occitanus Paris (BoP) in order to find new peptides, which could enlarge our structure-function relationship knowledge on the Kv1.3 blocker Kaliotoxin (KTX) that belongs to the α-KTx3.1 family. Indeed and since more a decade, KTX is widely used by international investigators because it exhibits a quite sharp specificity and a high-affinity for the Kv1.3 channel, which is not only a neuronal channel but also a therapeutic target for diverse autoimmune diseases such as multiple sclerosis, type 1 diabetes, and rheumatoid arthritis. The BoP venom was first investigated using HPLC and MALDI-TOF/MS. Further, the HPLC fractions were screened by ELISA with antibodies raised against KTX. These antibodies recognized at least three components toxic in mice by intracerebroventricular injection. They were further pharmacologically characterized by competition using (125)I-KTX bound to its specific binding sites on rat brain synaptosomes. A single component (4161 Da) inhibited totally the (125)I-KTX binding and with high-affinity (IC50 = 0.1 nM), while the two other components poorly competed with (IC50 > 100 nM). These toxins were sequenced in full by Edman's degradation. The high-affinity ligand (BoPKTX) shares 86% sequence identity with KTX and was classified as toxin α-KTx3.17. The two others peptides (BoP1 and BoP2, 4093 Da and 4121 Da, respectively) only differ by a Lys/Arg mutation. Their amino acid sequences were related to Martentoxin, which has been characterized from the Chinese scorpion Buthus martenzi Karch and described as both a BKCa and Kv1.3 blocker. Accordingly, they belong to the α-KTx16 family.

Accelerated evolution and functional divergence of scorpion short-chain K+ channel toxins after speciation.

The α-KTx14 subfamily of scorpion toxins is a group of short-chain polypeptides affecting K(+) channels, including five known members which are restrictedly distributed in Mesobuthus martensii. Here, we describe seven new α-KTx14 peptides from M. martensii and its sibling species Mesobuthus eupeus, two of which (termed MarKTX-3 and MeuKTX-1) were chemically synthesized and refolded for structural and functional studies. Electrophysiological recordings of effects of these two peptides on an array of voltage-gated potassium channels revealed that MarKTX-3 was capable of inhibiting five mammalian K(v)1 isoforms (rK(v)1.1-rK(v)1.5) and the Drosophila Shaker channel with low potency whereas MeuKTX-1 lacks such activity. Circular dichroism spectroscopy analysis combined with homology modeling demonstrates that MarKTX-3 and MeuKTX-1 both adopt a similar cysteine-stabilized α-helical and ß-sheet fold. Evolutionary analysis indicates accelerated amino acid substitutions in the mature-peptide-encoding regions of orthologous α-KTx14 peptides after speciation, thereby providing evidences for adaptive evolution and functional divergence of this subfamily.