Denmotoxin, a three-finger toxin from the colubrid snake Boiga dendrophila (Mangrove Catsnake) with bird-specific activity

Preliminary Biochemical and Venomic Characterization of the Venom of Phalotris lemniscatus (Serpentes, Colubridae)

BACKGROUND

For many decades, research on snake venom toxinology focused mainly on the venoms of Viperidae and Elapidae species, which were traditionally the only ones considered as venomous. However, much less interest has been given to the venom produced by opisthoglyphous colubrid snakes, since they were typically considered of no clinical relevance.

OBJECTIVE

The aim of this work is to perform a preliminary biochemical and venomic characterization of the venom of the colubrid snake Phalotris lemniscatus, a species that has been responsible for two relevant cases of envenomation in Uruguay.

METHODS

We extracted venom from collected specimens and performed different biochemical and proteomic assays to understand its toxin composition.

RESULTS

We found that the venom of P. lemniscatus is composed of protein families typically present in snake venoms, such as metallo and serine preoteases, L-amino acid oxidases, phospholipases A2s, Ctype lectines-like, Kunitz-type proteins and three-finger toxins. Activity assays demonstrated a highly active gelatinolytic component as well as a potent capability to induce blood coagulation.

CONCLUSION

The results indicate that the venom of P. lemniscatus contains hemotoxic activities and components that resemble those found in Viperidae (Bothrops) snakes and that can induce a clinically relevant accident. Further studies are needed to better understand the venom composition of this colubrid snake and its most active compounds.

Neurobiology and therapeutic utility of neurotoxins targeting postsynaptic mechanisms of neuromuscular transmission

The neuromuscular junction (NMJ) is the principal site for the translation of motor neurochemical signals to muscle activity. Therefore, the release and sensing machinery of acetylcholine (ACh) along with muscle contraction are two of the main targets of natural toxins and pathogens, causing paralysis. Given pharmacology and medical advances, the active ingredients of toxins that target postsynaptic mechanisms have become of major interest, showing promise as drug leads. Herein, we review key facets of prevalent toxins modulating the mechanisms of ACh sensing and generation of the postsynaptic response, with muscle contraction. We consider the correlation between their outstanding selectivity and potency plus effects on motor function, and discuss emerging data advocating their usage for the development of therapies alleviating neuromuscular dysfunction.

The Urgent Need to Develop Novel Strategies for the Diagnosis and Treatment of Snakebites

Snakebite envenoming (SBE) is a priority neglected tropical disease, which kills in excess of 100,000 people per year. Additionally, many millions of survivors also suffer through disabilities and long-term health consequences. The only treatment for SBE, antivenom, has a number of major associated problems, not least, adverse reactions and limited availability. This emphasises the necessity for urgent improvements to the management of this disease. Administration of antivenom is too frequently based on symptomatology, which results in wasting crucial time. The majority of SBE-affected regions rely on broad-spectrum polyvalent antivenoms that have a low content of case-specific efficacious immunoglobulins. Research into small molecular therapeutics such as varespladib/methyl-varespladib (PLA2 inhibitors) and batimastat/marimastat (metalloprotease inhibitors) suggest that such adjunctive treatments could be hugely beneficial to victims. Progress into toxin-specific monoclonal antibodies as well as alternative binding scaffolds such as aptamers hold much promise for future treatment strategies. SBE is not implicit during snakebite, due to venom metering. Thus, the delay between bite and symptom presentation is critical and when symptoms appear it may often already be too late to effectively treat SBE. The development of reliable diagnostical tools could therefore initiate a paradigm shift in the treatment of SBE. While the complete eradication of SBE is an impossibility, mitigation is in the pipeline, with new treatments and diagnostics rapidly emerging. Here we critically review the urgent necessity for the development of diagnostic tools and improved therapeutics to mitigate the deaths and disabilities caused by SBE.

Evolution and Medical Significance of LU Domain-Containing Proteins

Proteins containing Ly6/uPAR (LU) domains exhibit very diverse biological functions and have broad taxonomic distributions in eukaryotes. In general, they adopt a characteristic three-fingered folding topology with three long loops projecting from a disulfide-rich globular core. The majority of the members of this protein domain family contain only a single LU domain, which can be secreted, glycolipid anchored, or constitute the extracellular ligand binding domain of type-I membrane proteins. Nonetheless, a few proteins contain multiple LU domains, for example, the urokinase receptor uPAR, C4.4A, and Haldisin. In the current review, we will discuss evolutionary aspects of this protein domain family with special emphasis on variations in their consensus disulfide bond patterns. Furthermore, we will present selected cases where missense mutations in LU domain-containing proteins leads to dysfunctional proteins that are causally linked to genesis of human disease.

Toxins for decoding interface selectivity in nicotinic acetylcholine receptors

Nicotinic acetylcholine receptors (nAChRs) are pentameric ligand-gated ion channels that play crucial roles in neurotransmission and regulate complex processes in brain functions, including anxiety, learning and memory, food intake, drug addiction, cognition and nociception. To perform these and other functions, a diverse array of nAChR subtypes are generated by homomeric or heteromeric assembly of 17 homologous nAChR subunits. Agonists, acetylcholine and nicotine, bind to the interface formed between two α subunits and between α and non-α subunits to activate the nAChR and allow cation influx. The diversity of subunit interfaces determines the channel properties, the responses to different agonists/antagonists, desensitization and downstream signaling and thus, define specialized properties and functions. Over the last several decades, snake venom neurotoxins have contributed to the purification, localization and characterization of molecular details of various nAChRs. Utkin et al. have described the purification and characterization of αδ-bungarotoxins, a novel class of neurotoxins in a recent paper published in the Biochemical Journal [Biochem. J. (2019) 476, 1285-1302]. These toxins from Bungarus candidus venom preferably bind to α-δ site with two orders of magnitude higher affinity compared with α-γ or α-ε sites. The subtle changes in the structure of αδ-bungarotoxins led to variation in interface selectivity. Such new classes of antagonists will offer us great opportunity to delineate the pharmacophores and design new highly selective antagonists. Thus, their findings provide new impetus to re-evaluate molecular details of pharmacological properties of α-neurotoxins with careful consideration towards subtype-, interface- and species-selectivity.

Identification of a α-helical molten globule intermediate and structural characterization of β-cardiotoxin, an all β-sheet protein isolated from the venom of Ophiophagus hannah (king cobra)

β-Cardiotoxin is a novel member of the snake venom three-finger toxin (3FTX) family. This is the first exogenous protein to antagonize β-adrenergic receptors and thereby causing reduction in heart rates (bradycardia) when administered into animals, unlike the conventional cardiotoxins as reported earlier. 3FTXs are stable all β-sheet peptides with 60-80 amino acid residues. Here, we describe the three-dimensional crystal structure of β-cardiotoxin together with the identification of a molten globule intermediate in the unfolding pathway of this protein. In spite of the overall structural similarity of this protein with conventional cardiotoxins, there are notable differences observed at the loop region and in the charge distribution on the surface, which are known to be critical for cytolytic activity of cardiotoxins. The molten globule intermediate state present in the thermal unfolding pathway of β-cardiotoxin was however not observed during the chemical denaturation of the protein. Interestingly, circular dichroism (CD) and NMR studies revealed the presence of α-helical secondary structure in the molten globule intermediate. These results point to substantial conformational plasticity of β-cardiotoxin, which might aid the protein in responding to the sometimes conflicting demands of structure, stability, and function during its biological lifetime.

Novel long-chain neurotoxins from Bungarus candidus distinguish the two binding sites in muscle-type nicotinic acetylcholine receptors

αδ-Bungarotoxins, a novel group of long-chain α-neurotoxins, manifest different affinity to two agonist/competitive antagonist binding sites of muscle-type nicotinic acetylcholine receptors (nAChRs), being more active at the interface of α-δ subunits. Three isoforms (αδ-BgTx-1-3) were identified in Malayan Krait (Bungarus candidus) from Thailand by genomic DNA analysis; two of them (αδ-BgTx-1 and 2) were isolated from its venom. The toxins comprise 73 amino acid residues and 5 disulfide bridges, being homologous to α-bungarotoxin (α-BgTx), a classical blocker of muscle-type and neuronal α7, α8, and α9α10 nAChRs. The toxicity of αδ-BgTx-1 (LD₅₀ = 0.17-0.28 µg/g mouse, i.p. injection) is essentially as high as that of α-BgTx. In the chick biventer cervicis nerve-muscle preparation, αδ-BgTx-1 completely abolished acetylcholine response, but in contrast with the block by α-BgTx, acetylcholine response was fully reversible by washing. αδ-BgTxs, similar to α-BgTx, bind with high affinity to α7 and muscle-type nAChRs. However, the major difference of αδ-BgTxs from α-BgTx and other naturally occurring α-neurotoxins is that αδ-BgTxs discriminate the two binding sites in the Torpedo californica and mouse muscle nAChRs showing up to two orders of magnitude higher affinity for the α-δ site as compared with α-ε or α-γ binding site interfaces. Molecular modeling and analysis of the literature provided possible explanations for these differences in binding mode; one of the probable reasons being the lower content of positively charged residues in αδ-BgTxs. Thus, αδ-BgTxs are new tools for studies on nAChRs.

Did evolution create a flexible ligand-binding cavity in the urokinase receptor through deletion of a plesiotypic disulfide bond?

The urokinase receptor (uPAR) is a founding member of a small protein family with multiple Ly6/uPAR (LU) domains. The motif defining these LU domains contains five plesiotypic disulfide bonds stabilizing its prototypical three-fingered fold having three protruding loops. Notwithstanding the detailed knowledge on structure-function relationships in uPAR, one puzzling enigma remains unexplored. Why does the first LU domain in uPAR (DI) lack one of its consensus disulfide bonds, when the absence of this particular disulfide bond impairs the correct folding of other single LU domain-containing proteins? Here, using a variety of contemporary biophysical methods, we found that reintroducing the two missing half-cystines in uPAR DI caused the spontaneous formation of the corresponding consensus 7-8 LU domain disulfide bond. Importantly, constraints due to this cross-link impaired (i) the binding of uPAR to its primary ligand urokinase and (ii) the flexible interdomain assembly of the three LU domains in uPAR. We conclude that the evolutionary deletion of this particular disulfide bond in uPAR DI may have enabled the assembly of a high-affinity urokinase-binding cavity involving all three LU domains in uPAR. Of note, an analogous neofunctionalization occurred in snake venom α-neurotoxins upon loss of another pair of the plesiotypic LU domain half-cystines. In summary, elimination of the 7-8 consensus disulfide bond in the first LU domain of uPAR did have significant functional and structural consequences.

Habu coagulotoxicity: Clinical implications of the functional diversification of Protobothrops snake venoms upon blood clotting factors

Venom can affect any part of the body reachable via the bloodstream. Toxins which specifically act upon the coagulation cascade do so either by anticoagulant or procoagulant mechanisms. Here we investigated the coagulotoxic effects of six species within the medically important pit viper genus Protobothrops (Habu) from the Chinese mainland and Japanese islands, a genus known to produce hemorrhagic shock in envenomed patients. Differential coagulotoxicity was revealed: P. jerdonii and P. mangshanensis produced an overall net anticoagulant effect through the pseudo-procoagulant clotting of fibrinogen; P. flavoviridis and P. tokarensis exhibit a strong anticoagulant activity through the destructive cleavage of fibrinogen; and while P. elegans and P. mucrosquamatus both cleaved the A-alpha and B-beta chains of fibrinogen they did not exhibit strong anticoagulant activity. These variations in coagulant properties were congruent with phylogeny, with the closest relatives exhibiting similar venom effects in their action upon fibrinogen. Ancestral state reconstruction indicated that anticoagulation mediated by pseudo-procoagulant cleavage of fibrinogen is the basal state, while anticoagulation produced by destructive cleavage of fibrinogen is the derived state within this genus. This is the first in depth study of its kind highlighting extreme enzymatic variability, functional diversification and clotting diversification within one genus surrounding one target site, governed by variability in co-factor dependency. The documentation that the same net overall function, anticoagulation, is mediated by differential underlying mechanics suggests limited antivenom cross-reactivity, although this must be tested in future work. These results add to the body of knowledge necessary to inform clinical management of the envenomed patient.

Coagulotoxicity of Bothrops (Lancehead Pit-Vipers) Venoms from Brazil: Differential Biochemistry and Antivenom Efficacy Resulting from Prey-Driven Venom Variation

Lancehead pit-vipers (Bothrops genus) are an extremely diverse and medically important group responsible for the greatest number of snakebite envenomations and deaths in South America. Bothrops atrox (common lancehead), responsible for majority of snakebites and related deaths within the Brazilian Amazon, is a highly adaptable and widely distributed species, whose venom variability has been related to several factors, including geographical distribution and habitat type. This study examined venoms from four B. atrox populations (Belterra and Santarém, PA; Pres. Figueiredo, AM and São Bento, MA), and two additional Bothrops species (B. jararaca and B. neuwiedi) from Southeastern region for their coagulotoxic effects upon different plasmas (human, amphibian, and avian). The results revealed inter– and intraspecific variations in coagulotoxicity, including distinct activities between the three plasmas, with variations in the latter two linked to ecological niche occupied by the snakes. Also examined were the correlated biochemical mechanisms of venom action. Significant variation in the relative reliance upon the cofactors calcium and phospholipid were revealed, and the relative dependency did not significantly correlate with potency. Relative levels of Factor X or prothrombin activating toxins correlated with prey type and prey escape potential. The antivenom was shown to perform better in neutralising prothrombin activation activity than neutralising Factor X activation activity. Thus, the data reveal new information regarding the evolutionary selection pressures shaping snake venom evolution, while also having significant implications for the treatment of the envenomed patient. These results are, therefore, an intersection between evolutionary biology and clinical medicine.

Three-Finger Toxin Diversification in the Venoms of Cat-Eye Snakes (Colubridae: Boiga)

The Asian genus Boiga (Colubridae) is among the better studied non-front-fanged snake lineages, because their bites have minor, but noticeable, effects on humans. Furthermore, B. irregularis has gained worldwide notoriety for successfully invading Guam and other nearby islands with drastic impacts on the local bird populations. One of the factors thought to allow B. irregularis to become such a noxious pest is irditoxin, a dimeric neurotoxin composed of two three-finger toxins (3FTx) joined by a covalent bond between two newly evolved cysteines. Irditoxin is highly toxic to diapsid (birds and reptiles) prey, but roughly 1000 × less potent to synapsids (mammals). Venom plays an important role in the ecology of all species of Boiga, but it remains unknown if any species besides B. irregularis produce irditoxin-like dimeric toxins. In this study, we use transcriptomic analyses of venom glands from five species [B. cynodon, B. dendrophila dendrophila, B. d. gemmicincta, B. irregularis (Brisbane population), B. irregularis (Sulawesi population), B. nigriceps, B. trigonata] and proteomic analyses of B. d. dendrophila and a representative of the sister genus Toxicodryas blandingii to investigate the evolutionary history of 3FTx within Boiga and its close relative. We found that 92.5% of Boiga 3FTx belong to a single clade which we refer to as denmotoxin-like because of the close relation between these toxins and the monomeric denmotoxin according to phylogenetic, sequence clustering, and protein similarity network analyses. We show for the first time that species beyond B. irregularis secrete 3FTx with additional cysteines in the same position as both the A and B subunits of irditoxin. Transcripts with the characteristic mutations are found in B. d. dendrophila, B. d. gemmicincta, B. irregularis (Brisbane population), B. irregularis (Sulawesi population), and B. nigriceps. These results are confirmed by proteomic analyses that show direct evidence of dimerization within the venom of B. d. dendrophila, but not T. blandingii. Our results also suggest the possibility of novel dimeric toxins in other genera such as Telescopus and Trimorphodon. All together, this suggests that the origin of these peculiar 3FTx is far earlier than was appreciated and their evolutionary history has been complex.

Nanobodies as novel therapeutic agents in envenomation

BACKGROUND

An effective therapy against envenoming should be a priority in view of the high number scorpion stings and snakebites. Serum therapy is still widely applied to treat the envenomation victims; however this approach suffers from several shortcomings. The employment of monoclonal antibodies might be an outcome as these molecules are at the core of a variety of applications from protein structure determination to cancer treatment. The progress of activities in the twilight zone between genetic and antibody engineering have led to the development of a unique class of antibody fragments. These molecules possess several benefits and lack many possible disadvantages over classical antibodies. Within recombinant antibody formats, nanobodies or single domain antigen binding fragments derived from heavy chain only antibodies in camelids occupy a privileged position.

SCOPE OF REVIEW

In this paper we will briefly review the common methods of envenomation treatment and focus on details of various in vivo research activities that investigate the performance of recombinant, monoclonal nanobodies in venom neutralization.

MAJOR CONCLUSIONS

Nanobodies bind to their cognate target with high specificity and affinity, they can be produced in large quantities from microbial expression systems and are very robust even when challenged with harsh environmental conditions. Upon administering, they rapidly distribute throughout the body and seem to be well tolerated in humans posing low immunogenicity.

GENERAL SIGNIFICANCE

Scorpion and snake envenomation is a major issue in developing countries and nanobodies as a venom-neutralizing agent can be considered as a valuable and promising candidate in envenomation therapy.

Transcriptome-facilitated proteomic characterization of rear-fanged snake venoms reveal abundant metalloproteinases with enhanced activity

High-throughput technologies were used to identify venom gland toxin expression and to characterize the venom proteomes of two rear-fanged snakes, Ahaetulla prasina (Asian Green Vine Snake) and Borikenophis portoricensis (Puerto Rican Racer). Sixty-nine complete toxin-coding transcripts from 12 venom protein superfamilies (A. prasina) and 50 complete coding transcripts from 11 venom protein superfamilies (B. portoricensis) were identified in the venom glands. However, only 18% (A. prasina) and 32% (B. portoricensis) of the translated protein isoforms were detected in the proteome of these venoms. Both venom gland transcriptomes and venom proteomes were dominated by P-III metalloproteinases. Three-finger toxins, cysteine-rich secretory proteins, and C-type lectins were present in moderate amounts, but other protein superfamilies showed very low abundances. Venoms contained metalloproteinase activity comparable to viperid snake venom levels, but other common venom enzymes were absent or present at negligible levels. Western blot analysis showed metalloproteinase and cysteine-rich secretory protein epitopes shared with the highly venomous Boomslang (Dispholidus typus). The abundance of metalloproteinases emphasizes the important trophic role of these toxins. Comprehensive, transcriptome-informed definition of proteomes and functional characterization of venom proteins in rear-fanged snake families help to elucidate toxin evolution and provide models for protein structure-function analyses.

Adaptive evolution of distinct prey-specific toxin genes in rear-fanged snake venom

Venom proteins evolve rapidly, and as a trophic adaptation are excellent models for predator-prey evolutionary studies. The key to a deeper understanding of venom evolution is an integrated approach, combining prey assays with analysis of venom gene expression and venom phenotype. Here, we use such an approach to study venom evolution in the Amazon puffing snake, Spilotes sulphureus, a generalist feeder. We identify two novel three-finger toxins: sulditoxin and sulmotoxin 1. These new toxins are not only two of the most abundant venom proteins, but are also functionally intriguing, displaying distinct prey-specific toxicities. Sulditoxin is highly toxic towards lizard prey, but is non-toxic towards mammalian prey, even at greater than 22-fold higher dosage. By contrast, sulmotoxin 1 exhibits the reverse trend. Furthermore, evolutionary analysis and structural modelling show highest sequence variability in the central loop of these proteins, probably driving taxon-specific toxicity. This is, to our knowledge, the first case in which a bimodal and contrasting pattern of toxicity has been shown for proteins in the venom of a single snake in relation to diet. Our study is an example of how toxin gene neofunctionalization can result in a venom system dominated by one protein superfamily and still exhibit flexibility in prey capture efficacy.

Orphan Three-Finger Toxins Bind at Tissue Factor-Factor VIIa Interface to Inhibit Factor X Activation: Identification of Functional Site by Docking

Three-finger toxins (3FTxs) contribute to toxicity of venomous snakes belonging to the family Elapidae. Currently, functions of a considerable proportion of 3FTxs are still unknown. Here, we describe the function of orphan group I 3FTxs consisting of four members. We also identified a new member of this group by sequencing a transcript isolated from Naja naja venom. This transcript, named najalexin, is identical to that previously described 3FTx from Naja atra venom gland, and shared high sequence identity with ringhalexin from Hemachatus haemachatus and a hypothetical protein from Ophiophagus hannah (here named as ophiolexin). The three-dimensional structure, as predicted by molecular modeling, showed that najalexin and ophiolexin share the same conserved structural organization as ringhalexin and other 3FTxs. Since ringhalexin inhibits the activation of factor X by the tissue factor-factor VIIa complex (TF-FVIIa), we evaluated the interaction of this group of 3FTxs with all components using in silico protein-protein docking studies. The binding of orphan group I 3FTxs to TF-FVIIa complex appears to be driven by their interaction with TF. They bind to fibronectin domain closer to the 170-loop of the FVIIa heavy chain to inhibit factor X activation. The docking studies reveal that functional site residues Tyr7, Lys9, Glu12, Lys26, Arg34, Leu35, Arg40, Val55, Asp56, Cys57, Cys58, and Arg65 on these 3FTxs are crucial for interaction. In silico replacement of these residues by Ala resulted in significant effects in the binding energies. Furthermore, these functional residues are not found in other groups of 3FTxs, which exhibit distinct pharmacological properties.

Neurotoxicity fingerprinting of venoms using on-line microfluidic AChBP profiling

Venoms from snakes are rich sources of highly active proteins with potent affinity towards a variety of enzymes and receptors. Of the many distinct toxicities caused by envenomation, neurotoxicity plays an important role in the paralysis of prey by snakes as well as by venomous sea snails and insects. In order to improve the analytical discovery component of venom toxicity profiling, this paper describes the implementation of microfluidic high-resolution screening (HRS) to obtain neurotoxicity fingerprints from venoms that facilitates identification of the neurotoxic components of envenomation. To demonstrate this workflow, 47 snake venoms were profiled using the acetylcholine binding protein (AChBP) to mimic the target of neurotoxic proteins, in particular nicotinic acetylcholine receptors (nAChRs). In the microfluidic HRS system, nanoliquid chromatographic (nanoLC) separations were on-line connected to both AChBP profiling and parallel mass spectrometry (MS). For virtually all neurotoxic elapid snake venoms tested, we obtained bioactivity fingerprints showing major and minor bioactive zones containing masses consistent with three-finger toxins (3FTxs), whereas, viperid and colubrid venoms showed little or no detectable bioactivity. Our findings demonstrate that venom interactions with AChBP correlate with the severity of neurotoxicity observed following human envenoming by different snake species. We further, as proof of principle, characterized bioactive venom peptides from a viperid (Daboia russelli) and an elapid (Aspidelaps scutatus scutatus) snake by nanoLC-MS/MS, revealing that different toxin classes interact with the AChBP, and that this binding correlates with the inhibition of α7-nAChR in calcium-flux cell-based assays. The on-line post-column binding assay and subsequent toxin characterization methodologies described here provide a new in vitro analytic platform for rapidly investigating neurotoxic snake venom proteins.

Synergistic strategies of predominant toxins in snake venoms

Synergism is a significant phenomenon present in snake venoms that may be an evolving strategy to potentiate toxicities. Synergism exists between different toxins or toxin complexes in various snake venoms, with phospholipaseA₂s (PLA₂s) (toxins or subunits) the main enablers. The predominant toxins, snake venom PLA₂s, metalloproteases (SVMPs), serine proteases (SVSPs) and three-finger toxins (3FTxs), play essential roles in synergistic processes. The hypothetical mechanisms of synergistic effect can be generalized under the effects of amplification and chaperoning. The Toxicity Score is among the few quantitative methods to assess synergism. Selection of toxins involved in synergistically enhanced toxicity as the targets are important for development of novel antivenoms or inhibitors.

Rattling the border wall: Pathophysiological implications of functional and proteomic venom variation between Mexican and US subspecies of the desert rattlesnake Crotalus scutulatus

While some US populations of the Mohave rattlesnake (Crotalus scutulatus scutulatus) are infamous for being potently neurotoxic, the Mexican subspecies C. s. salvini (Huamantlan rattlesnake) has been largely unstudied beyond crude lethality testing upon mice. In this study we show that at least some populations of this snake are as potently neurotoxic as its northern cousin. Testing of the Mexican antivenom Antivipmyn showed a complete lack of neutralisation for the neurotoxic effects of C. s. salvini venom, while the neurotoxic effects of the US subspecies C. s. scutulatus were time-delayed but ultimately not eliminated. These results document unrecognised potent neurological effects of a Mexican snake and highlight the medical importance of this subspecies, a finding augmented by the ineffectiveness of the Antivipmyn antivenom. These results also influence our understanding of the venom evolution of Crotalus scutulatus, suggesting that neurotoxicity is the ancestral feature of this species, with the US populations which lack neurotoxicity being derived states.

Transcriptomics-guided bottom-up and top-down venomics of neonate and adult specimens of the arboreal rear-fanged Brown Treesnake, Boiga irregularis, from Guam

The Brown Treesnake (Boiga irregularis) is an arboreal, nocturnal, rear-fanged venomous snake native to northern and eastern regions of Australia, Papua New Guinea and the Solomon Islands. It was inadvertently introduced onto the island of Guam during the late 1940's to early 1950's, and it has caused massive declines and extirpations of the native bird, lizard, and mammal populations. In the current study, we report the characterization of the venom proteome of an adult and a neonate B. irregularis specimens from Guam by a combination of venom gland transcriptomic and venomic analyses. Venom gland transcriptomic analysis of an adult individual identified toxins belonging to 18 protein families, with three-finger toxin isoforms being the most abundantly expressed transcripts, comprising 94% of all venom protein transcript reads. Transcripts for PIII-metalloproteinases, C-type lectins, cysteine-rich secretory proteins, acetylcholinesterases, natriuretic peptides, ficolins, phospholipase A₂ (PLA₂) inhibitors, PLA₂s, vascular endothelial growth factors, Kunitz-type protease inhibitors, cystatins, phospholipase Bs, cobra venom factors, waprins, SVMP inhibitors, matrix metalloproteinases, and hyaluronidases were also identified, albeit, at very low abundances ranging from 0.05% to 1.7% of the transcriptome. The venom proteomes of neonate and adult B. irregularis were also both overwhelmingly (78 and 84%, respectively) dominated by monomeric and dimeric 3FTxs, followed by moderately abundant (21% (N) and 13% (A)) CRISPs, low abundance (1% (N), 3% (A)) PIII-SVMPs, and very low abundance (<0.01%) PLA₂ and SVMP inhibitors. The differences in relative toxin abundances identified between neonate and adult snakes likely correlates to shifts in prey preference between the two age classes, from nearly-exclusively lizards to lizards, birds and small mammals. Immunoaffinity antivenomics with experimentally designed rabbit anti-Brown Treesnake (anti-BTS) venom IgGs against homologous venom from adult snakes demonstrated that CRISPs, PIII-SVMPs, and 60-70% of 3FTxs were effectively immunocaptured. Western blot analysis showed that all venom proteins were recognized by anti-BTS IgGs, and cross-reactivity with other rear-fanged snake venoms was also observed. Incubation of anti-BTS venom IgGs with crude B. irregularis venom resulted in a significant decrease in proteolytic (SVMP) activity against azocasein. These results provide the first comparative venomic and anti-venomic analysis of neonate and adult B. irregularis from Guam, further highlighting evolutionary trends in venom composition among rear-fanged venomous snakes.

SIGNIFICANCE PARAGRAPH

The Brown Treesnake (Boiga irregularis) has caused extensive ecological and economic damage to the island of Guam where it has become a classic example of the negative impacts of invasive species. In the current study, we report the first combined transcriptomic and proteomic analysis of B. irregularis venom of Guam origin. The transcriptome of an adult snake contained toxin sequences belonging to 18 protein families, with three-finger toxin (3FTx) isoforms being the most abundant and representing 94% of all venom protein transcript reads. Our bottom-up and top-down venomic analyses confirmed that 3FTxs are the major components of B. irregularis venom, and a comparative analysis of neonate and adult venoms demonstrate a clear ontogenetic shift in toxin abundance, likely driven by dietary variation between the two age classes. Second-generation antivenomics and Western blot analysis using purified anti-Brown Treesnake rabbit serum IgGs (anti-BTS IgGs) showed strong immunoreactivity toward B. irregularis venom. Interestingly, our anti-BTS IgGs did not cross-react with 3FTxs found in several other rear-fanged snake venoms, or against 3FTxs in the venom of the elapid Ophiophagus hannah, indicating that epitopes in these 3FTx molecules are quite distinct.

Identification and characterisation of novel inhibitors on extrinsic tenase complex from Bungarus fasciatus (banded krait) venom

Snake venoms are excellent sources of pharmacologically active proteins and peptides, and hence are potential sources of leads for drug developments. It has been previously established that krait (Bungarus genus) venoms contain mainly neurotoxins. A screening for anticoagulants showed that Bungarus fasciatus venom exhibits potent anticoagulant effect in standard clotting assays. Through sequential fractionation of the venom by size exclusion and high performance liquid chromatographies, coupled with functional screening for anticoagulant activities, we have isolated and purified two anticoagulant proteins, termed BF-AC1 ( Bungarus fasciatus anticoagulant 1) and BFAC2. They have potent inhibitory activities (IC50 of 10 nM) on the extrinsic tenase complex. Structurally, these proteins each has two subunits covalently held together by disulfide bond(s). The N-terminal sequences of the individual subunits of BF-AC1 and BF-AC2 showed that the larger subunit is homologous to phospholipase A2, while the smaller subunit is homologous to Kunitz type serine proteinase inhibitor. Functionally, in addition to their anticoagulant activity, these proteins showed presynaptic neurotoxic effects in both in vivo and ex vivo experiments. Thus, BF-AC1 and BF-AC2 are structurally and functionally similar to β-bungarotoxins, a class of neurotoxins. The enzymatic activity of phospholipase A2 subunit plays a significant role in the anticoagulant activities. This is the first report on the anticoagulant activity of β-bungarotoxins and these results expand on the existing catalogue of haemostatically active snake venom proteins.

Animal toxins for channelopathy treatment

Ion channels are transmembrane proteins that allow passive flow of ions inside and/or outside of cells or cell organelles. Except mutations lead to nonfunctional protein production or abolished receptor entrance on the membrane surface an altered channel may have two principal conditions that can be corrected. The channel may conduct fewer ions through (loss-of-function mutations) or too many ions (gain-of-function mutations) compared to a normal channel. Toxins from animal venoms are specialised molecules that are generally oriented toward interactions with ion channels. This is a result of long coevolution between predators and their prey. On the molecular level, toxins activate or inhibit ion channels, so they are ideal molecules for restoring conductance in mutated channels. Another aspect of this long coevolution is that a broad variety of toxins have been fine tuned to recognize the channels of different species, keeping many amino acids substitution among sequences. Many peptide ligands with high selectivity to specific receptor subtypes have been isolated from animal venoms, some of which are absolutely non-toxic to humans and mammalians. It is expected that molecules that are selective to each known receptor can be found in animal venoms, but the pool of toxins currently does not override all receptors described as being involved in channelopathies. Modern investigating methods have enhanced the search process for selective ligands. One prominent method is a site-directed mutagenesis of existing toxins to change the selectivity or/and affinity to the selected receptor, which has shown positive results. This article is part of the Special Issue entitled 'Channelopathies.'

Correlation between ontogenetic dietary shifts and venom variation in Australian brown snakes (Pseudonaja)

Venom is a key evolutionary trait, as evidenced by its widespread convergent evolution across the animal kingdom. In an escalating prey-predator arms race, venoms evolve rapidly to guarantee predatory or defensive success. Variation in venom composition is ubiquitous among snakes. Here, we tested variation in venom activity on substrates relevant to blood coagulation among Pseudonaja (brown snake) species, Australian elapids responsible for the majority of medically important human envenomations in Australia. A functional approach was employed to elucidate interspecific variation in venom activity in all nine currently recognised species of Pseudonaja. Fluorometric enzymatic activity assays were performed to test variation in whole venom procoagulant activity among species. Analyses confirmed the previously documented ontogenetic shift from non-coagulopathic venom in juveniles to coagulopathic venom as adults, except for the case of P. modesta, which retains non-coagulopathic venom as an adult. These shifts in venom activity correlate with documented ontogenetic shifts in diet among brown snakes from specialisation on reptilian prey as juveniles (and throughout the life cycle of P. modesta), to a more generalised diet in adults that includes mammals. The results of this study bring to light findings relevant to both clinical and evolutionary toxinology.

Nicotinic acetylcholine receptor inhibitors derived from snake and snail venoms

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

A survey on some biochemical and pharmacological activities of venom from two Colombian colubrid snakes, Erythrolamprus bizona (Double-banded coral snake mimic) and Pseudoboa neuwiedii (Neuwied's false boa)

Colombian colubrid snake venoms have been poorly studied. They represent a great resource of biological, ecological, toxinological and pharmacological research. We assessed some enzymatic properties and neuromuscular effects of Erythrolamprus bizona and Pseudoboa neuwiedii venoms from Colombia. Proteolytic, amidolytic and phospholipase A₂ (PLA₂) activities were analyzed using colorimetric assays and the neuromuscular activity was analyzed in chick biventer cervicis (BC) preparations. The venom of both species showed very low PLA₂ and amidolytic activities; however, both exhibited high proteolytic activity, which in E. bizona venom surpassed that of P. neuwiedii venom. E. bizona and P. neuwiedii venoms provoked partial neuromuscular blockade, which was more prominent in P. neuwiedii venom. E. bizona venom (30 μg/ml) induced a significant potentiation of the contracture response to exogenous ACh (110 μM), which was not accompanied by twitch height alteration, whereas the highest venom concentration (100 μg/ml) inhibited contracture responses to both ACh and KCl (40 mM). In contrast, P. neuwiedii venom (30 and 100 μg/ml) caused significant reduction in the contracture responses to exogenous ACh and KCl. The morphological analyses showed high myotoxic effects in the muscle fibers of BC incubated with either venoms; however, they are more prominent in the P. neuwiedii venom. Our results suggest that the myotoxicity of the venom of the two Colombian species can be ascribed to their high proteolytic activity. An interesting data was the potentiation of the ACh-induced contracture, but not the twitch height, caused by E. bizona venom, at a concentration that is harmless to muscle fibers integrity. This phenomenon remains to be further elucidated, and suggest that a possible involvement of post-synaptic receptors cannot be discarded. This work is a contribution to expand the knowledge on colubrid venoms; it allows envisaging that the two venoms offer the potential to go further in the identification of their components and biological targets.

Comparative venom gland transcriptomics of Naja kaouthia (monocled cobra) from Malaysia and Thailand: elucidating geographical venom variation and insights into sequence novelty

Background

The monocled cobra (Naja kaouthia) is a medically important venomous snake in Southeast Asia. Its venom has been shown to vary geographically in relation to venom composition and neurotoxic activity, indicating vast diversity of the toxin genes within the species. To investigate the polygenic trait of the venom and its locale-specific variation, we profiled and compared the venom gland transcriptomes of N. kaouthia from Malaysia (NK-M) and Thailand (NK-T) applying next-generation sequencing (NGS) technology.

Methods

The transcriptomes were sequenced on the Illumina HiSeq platform, assembled and followed by transcript clustering and annotations for gene expression and function. Pairwise or multiple sequence alignments were conducted on the toxin genes expressed. Substitution rates were studied for the major toxins co-expressed in NK-M and NK-T.

Results and discussion

The toxin transcripts showed high redundancy (41–82% of the total mRNA expression) and comprised 23 gene families expressed in NK-M and NK-T, respectively (22 gene families were co-expressed). Among the venom genes, three-finger toxins (3FTxs) predominated in the expression, with multiple sequences noted. Comparative analysis and selection study revealed that 3FTxs are genetically conserved between the geographical specimens whilst demonstrating distinct differential expression patterns, implying gene up-regulation for selected principal toxins, or alternatively, enhanced transcript degradation or lack of transcription of certain traits. One of the striking features that elucidates the inter-geographical venom variation is the up-regulation of α-neurotoxins (constitutes ∼80.0% of toxin’s fragments per kilobase of exon model per million mapped reads (FPKM)), particularly the long-chain α-elapitoxin-Nk2a (48.3%) in NK-T but only 1.7% was noted in NK-M. Instead, short neurotoxin isoforms were up-regulated in NK-M (46.4%). Another distinct transcriptional pattern observed is the exclusively and abundantly expressed cytotoxin CTX-3 in NK-T. The findings suggested correlation with the geographical variation in proteome and toxicity of the venom, and support the call for optimising antivenom production and use in the region. Besides, the current study uncovered full and partial sequences of numerous toxin genes from N. kaouthia which have not been reported hitherto; these include N. kaouthia-specific l-amino acid oxidase (LAAO), snake venom serine protease (SVSP), cystatin, acetylcholinesterase (AChE), hyaluronidase (HYA), waprin, phospholipase B (PLB), aminopeptidase (AP), neprilysin, etc. Taken together, the findings further enrich the snake toxin database and provide deeper insights into the genetic diversity of cobra venom toxins.

The three-finger toxin fold: a multifunctional structural scaffold able to modulate cholinergic functions

Three-finger fold toxins are miniproteins frequently found in Elapidae snake venoms. This fold is characterized by three distinct loops rich in β-strands and emerging from a dense, globular core reticulated by four highly conserved disulfide bridges. The number and diversity of receptors, channels, and enzymes identified as targets of three-finger fold toxins is increasing continuously. Such manifold diversity highlights the specific adaptability of this fold for generating pleiotropic functions. Although this toxin superfamily disturbs many biological functions by interacting with a large diversity of molecular targets, the most significant target is the cholinergic system. By blocking the activity of the nicotinic and muscarinic acetylcholine receptors or by inhibiting the enzyme acetylcholinesterase, three-finger fold toxins interfere most drastically with neuromuscular junction functioning. Several of these toxins have become powerful pharmacological tools for studying the function and structure of their molecular targets. Most importantly, since dysfunction of these receptors/enzyme is involved in many diseases, exploiting the three-finger scaffold to create novel, highly specific therapeutic agents may represent a major future endeavor. This is an article for the special issue XVth International Symposium on Cholinergic Mechanisms.

Venomics: integrative venom proteomics and beyond*

Venoms are integrated phenotypes that evolved independently in, and are used for predatory and defensive purposes by, a wide phylogenetic range of organisms. The same principles that contribute to the evolutionary success of venoms, contribute to making the study of venoms of great interest in such diverse fields as evolutionary ecology and biotechnology. Evolution is profoundly contingent, and nature also reinvents itself continuosly. Changes in a complex phenotypic trait, such as venom, reflect the influences of prior evolutionary history, chance events, and selection. Reconstructing the natural history of venoms, particularly those of snakes, which will be dealt with in more detail in this review, requires the integration of different levels of knowledge into a meaningful and comprehensive evolutionary framework for separating stochastic changes from adaptive evolution. The application of omics technologies and other disciplines have contributed to a qualitative and quantitative advance in the road map towards this goal. In this review we will make a foray into the world of animal venoms, discuss synergies and complementarities of the different approaches used in their study, and identify current bottlenecks that prevent inferring the evolutionary mechanisms and ecological constraints that molded snake venoms to their present-day variability landscape.

What killed Karl Patterson Schmidt? Combined venom gland transcriptomic, venomic and antivenomic analysis of the South African green tree snake (the boomslang), Dispholidus typus

Background

Non-front-fanged colubroid snakes comprise about two-thirds of extant ophidian species. The medical significance of the majority of these snakes is unknown, but at least five species have caused life-threatening or fatal human envenomings. However, the venoms of only a small number of species have been explored.

Methods

A combined venomic and venom gland transcriptomic approach was employed to characterise of venom of Dispholidus typus (boomslang), the snake that caused the tragic death of Professor Karl Patterson Schmidt. The ability of CroFab™ antivenom to immunocapture boomslang venom proteins was investigated using antivenomics.

Results

Transcriptomic-assisted proteomic analysis identified venom proteins belonging to seven protein families: three-finger toxin (3FTx); phospholipase A₂ (PLA₂); cysteine-rich secretory proteins (CRISP); snake venom (SV) serine proteinase (SP); C-type lectin-like (CTL); SV metalloproteinases (SVMPs); and disintegrin-like/cysteine-rich (DC) proteolytic fragments. CroFab™ antivenom efficiently immunodepleted some boomslang SVMPs.

Conclusions

The present work is the first to address the overall proteomic profile of D. typus venom. This study allowed us to correlate the toxin composition with the toxic activities of the venom. The antivenomic analysis suggested that the antivenom available at the time of the unfortunate accident could have exhibited at least some immunoreactivity against the boomslang SVMPs responsible for the disseminated intravascular coagulation syndrome that caused K.P. Schmidt's fatal outcome.

General significance

This study may stimulate further research on other non-front-fanged colubroid snake venoms capable of causing life-threatening envenomings to humans, which in turn should contribute to prevent fatal human accidents, such as that unfortunately suffered by K.P. Schmidt.

•

The venom proteome of Dispholidus typus (boomslang) is reported.

•

Transcriptomic-assisted proteomic analysis identified venom proteins belonging to seven protein families.

•

Boomslang venom proteome is dominated (75%) by snake venom PIII-metalloproteinases (PIII-SVMPs).

•

CroFab™ antivenom efficiently immunodepleted some boomslang PIII-SVMPs.

Small Packages, Big Returns: Uncovering the Venom Diversity of Small Invertebrate Conoidean Snails

Venomous organisms used in research were historically chosen based on size and availability. This opportunity-driven strategy created a species bias in which snakes, scorpions, and spiders became the primary subjects of venom research. Increasing technological advancements have enabled interdisciplinary studies using genomics, transcriptomics, and proteomics to expand venom investigation to animals that produce small amounts of venom or lack traditional venom producing organs. One group of non-traditional venomous organisms that have benefitted from the rise of -omic technologies is the Conoideans. The Conoidean superfamily of venomous marine snails includes, the Terebridae, Turridae (s.l), and Conidae. Conoidea venom is used for both predation and defense, and therefore under strong selection pressures. The need for conoidean venom peptides to be potent and specific to their molecular targets has made them important tools for investigating cellular physiology and bioactive compounds that are beneficial to improving human health. A convincing case for the potential of Conoidean venom is made with the first commercially available conoidean venom peptide drug Ziconotide (Prialt®), an analgesic derived from Conus magus venom that is used to treat chronic pain in HIV and cancer patients. Investigation of conoidean venom using -omics technology provides significant insights into predator-driven diversification in biodiversity and identifies novel compounds for manipulating cellular communication, especially as it pertains to disease and disorders.

The Snake with the Scorpion's Sting: Novel Three-Finger Toxin Sodium Channel Activators from the Venom of the Long-Glanded Blue Coral Snake (Calliophis bivirgatus)

Millions of years of evolution have fine-tuned the ability of venom peptides to rapidly incapacitate both prey and potential predators. Toxicofera reptiles are characterized by serous-secreting mandibular or maxillary glands with heightened levels of protein expression. These glands are the core anatomical components of the toxicoferan venom system, which exists in myriad points along an evolutionary continuum. Neofunctionalisation of toxins is facilitated by positive selection at functional hotspots on the ancestral protein and venom proteins have undergone dynamic diversification in helodermatid and varanid lizards as well as advanced snakes. A spectacular point on the venom system continuum is the long-glanded blue coral snake (Calliophis bivirgatus), a specialist feeder that preys on fast moving, venomous snakes which have both a high likelihood of prey escape but also represent significant danger to the predator itself. The maxillary venom glands of C. bivirgatus extend one quarter of the snake's body length and nestle within the rib cavity. Despite the snake's notoriety its venom has remained largely unstudied. Here we show that the venom uniquely produces spastic paralysis, in contrast to the flaccid paralysis typically produced by neurotoxic snake venoms. The toxin responsible, which we have called calliotoxin (δ-elapitoxin-Cb1a), is a three-finger toxin (3FTx). Calliotoxin shifts the voltage-dependence of NaV1.4 activation to more hyperpolarised potentials, inhibits inactivation, and produces large ramp currents, consistent with its profound effects on contractile force in an isolated skeletal muscle preparation. Voltage-gated sodium channels (NaV) are a particularly attractive pharmacological target as they are involved in almost all physiological processes including action potential generation and conduction. Accordingly, venom peptides that interfere with NaV function provide a key defensive and predatory advantage to a range of invertebrate venomous species including cone snails, scorpions, spiders, and anemones. Enhanced activation or delayed inactivation of sodium channels by toxins is associated with the extremely rapid onset of tetanic/excitatory paralysis in envenomed prey animals. A strong selection pressure exists for the evolution of such toxins where there is a high chance of prey escape. However, despite their prevalence in other venomous species, toxins causing delay of sodium channel inhibition have never previously been described in vertebrate venoms. Here we show that NaV modulators, convergent with those of invertebrates, have evolved in the venom of the long-glanded coral snake. Calliotoxin represents a functionally novel class of 3FTx and a structurally novel class of NaV toxins that will provide significant insights into the pharmacology and physiology of NaV. The toxin represents a remarkable case of functional convergence between invertebrate and vertebrate venom systems in response to similar selection pressures. These results underscore the dynamic evolution of the Toxicofera reptile system and reinforces the value of using evolution as a roadmap for biodiscovery.

Understanding Biological Roles of Venoms Among the Caenophidia: The Importance of Rear-Fanged Snakes

Snake venoms represent an adaptive trophic response to the challenges confronting a limbless predator for overcoming combative prey, and this chemical means of subduing prey shows several dominant phenotypes. Many front-fanged snakes, particularly vipers, feed on various vertebrate and invertebrate prey species, and some of their venom components (e.g., metalloproteinases, cobratoxin) appear to have been selected for "broad-brush" incapacitation of different prey taxa. Using proteomic and genomic techniques, the compositional diversity of front-fanged snakes is becoming well characterized; however, this is not the case for most rear-fanged colubroid snakes. Because these species consume a high diversity of prey, and because venoms are primarily a trophic adaptation, important clues for understanding specific selective pressures favoring venom component composition will be found among rear-fanged snake venoms. Rear-fanged snakes typically (but not always) produce venoms with lower complexity than front-fanged snakes, and there are even fewer dominant (and, arguably, biologically most relevant) venom protein families. We have demonstrated taxon-specific toxic effects, where lizards and birds show high susceptibility while mammals are largely unaffected, for both Old World and New World rear-fanged snakes, strongly indicating a causal link between toxin evolution and prey preference. New data are presented on myotoxin a, showing that the extremely rapid paralysis induced by this rattlesnake toxin is specific for rodents, and that myotoxin a is ineffectual against lizards. Relatively few rear-fanged snake venoms have been characterized, and basic natural history data are largely lacking, but directed sampling of specialized species indicates that novel compounds are likely among these specialists, particularly among those species feeding on invertebrate prey such as scorpions and centipedes. Because many of the more than 2200 species of colubroid snakes are rear-fanged, and many possess a Duvernoy's venom gland, understanding the nature of their venoms is foundational to understanding venom evolution in advanced snakes.

Colubrid Venom Composition: An -Omics Perspective

Snake venoms have been subjected to increasingly sensitive analyses for well over 100 years, but most research has been restricted to front-fanged snakes, which actually represent a relatively small proportion of extant species of advanced snakes. Because rear-fanged snakes are a diverse and distinct radiation of the advanced snakes, understanding venom composition among "colubrids" is critical to understanding the evolution of venom among snakes. Here we review the state of knowledge concerning rear-fanged snake venom composition, emphasizing those toxins for which protein or transcript sequences are available. We have also added new transcriptome-based data on venoms of three species of rear-fanged snakes. Based on this compilation, it is apparent that several components, including cysteine-rich secretory proteins (CRiSPs), C-type lectins (CTLs), CTLs-like proteins and snake venom metalloproteinases (SVMPs), are broadly distributed among "colubrid" venoms, while others, notably three-finger toxins (3FTxs), appear nearly restricted to the Colubridae (sensu stricto). Some putative new toxins, such as snake venom matrix metalloproteinases, are in fact present in several colubrid venoms, while others are only transcribed, at lower levels. This work provides insights into the evolution of these toxin classes, but because only a small number of species have been explored, generalizations are still rather limited. It is likely that new venom protein families await discovery, particularly among those species with highly specialized diets.

Canopy Venom: Proteomic Comparison among New World Arboreal Pit-Viper Venoms

Central and South American pitvipers, belonging to the genera Bothrops and Bothriechis, have independently evolved arboreal tendencies. Little is known regarding the composition and activity of their venoms. In order to close this knowledge gap, venom proteomics and toxin activity of species of Bothriechis, and Bothrops (including Bothriopsis) were investigated through established analytical methods. A combination of proteomics and bioactivity techniques was used to demonstrate a similar diversification of venom composition between large and small species within Bothriechis and Bothriopsis. Increasing our understanding of the evolution of complex venom cocktails may facilitate future biodiscoveries.

Full-Length Venom Protein cDNA Sequences from Venom-Derived mRNA: Exploring Compositional Variation and Adaptive Multigene Evolution

Envenomation of humans by snakes is a complex and continuously evolving medical emergency, and treatment is made that much more difficult by the diverse biochemical composition of many venoms. Venomous snakes and their venoms also provide models for the study of molecular evolutionary processes leading to adaptation and genotype-phenotype relationships. To compare venom complexity and protein sequences, venom gland transcriptomes are assembled, which usually requires the sacrifice of snakes for tissue. However, toxin transcripts are also present in venoms, offering the possibility of obtaining cDNA sequences directly from venom. This study provides evidence that unknown full-length venom protein transcripts can be obtained from the venoms of multiple species from all major venomous snake families. These unknown venom protein cDNAs are obtained by the use of primers designed from conserved signal peptide sequences within each venom protein superfamily. This technique was used to assemble a partial venom gland transcriptome for the Middle American Rattlesnake (Crotalus simus tzabcan) by amplifying sequences for phospholipases A₂, serine proteases, C-lectins, and metalloproteinases from within venom. Phospholipase A₂ sequences were also recovered from the venoms of several rattlesnakes and an elapid snake (Pseudechis porphyriacus), and three-finger toxin sequences were recovered from multiple rear-fanged snake species, demonstrating that the three major clades of advanced snakes (Elapidae, Viperidae, Colubridae) have stable mRNA present in their venoms. These cDNA sequences from venom were then used to explore potential activities derived from protein sequence similarities and evolutionary histories within these large multigene superfamilies. Venom-derived sequences can also be used to aid in characterizing venoms that lack proteomic profiles and identify sequence characteristics indicating specific envenomation profiles. This approach, requiring only venom, provides access to cDNA sequences in the absence of living specimens, even from commercial venom sources, to evaluate important regional differences in venom composition and to study snake venom protein evolution.

This work demonstrates that full-length venom protein messenger RNAs are present in secreted venoms and can be used to acquire full-length protein sequences of toxins from both front-fanged (Elapidae, Viperidae) and rear-fanged (Colubridae) snake venoms, eliminating the need to use venom glands. Full-length transcripts were obtained from venom samples that were fresh, newly lyophilized, old, field desiccated or commercially prepared, representing a significant advance over previous attempts which produced only partial sequence transcripts. Transcripts for all major venom protein families (metalloproteinases, serine proteases, C-type lectins, phospholipases A₂ and three-finger toxins) responsible for clinically significant snakebite symptoms were obtained from venoms. These sequences aid in the identification and characterization of venom proteome profiles, allowing for the identification of peptide sequences, specific isoforms, and novel venom proteins. The application of this technique will help to provide venom protein sequences for many snake species, including understudied rear-fanged snakes. Venom protein transcripts offer important insights into potential snakebite envenomation profiles and the molecular evolution of venom protein multigene families. By requiring only venom to obtain venom protein cDNAs, the approach detailed here will provide access to cDNA-based protein sequences from commercial and other venom sources, facilitating study of snake venom protein composition and evolution.

Ringhalexin from Hemachatus haemachatus: A novel inhibitor of extrinsic tenase complex

Anticoagulant therapy is used for the prevention and treatment of thromboembolic disorders. Blood coagulation is initiated by the interaction of factor VIIa (FVIIa) with membrane-bound tissue factor (TF) to form the extrinsic tenase complex which activates FX to FXa. Thus, it is an important target for the development of novel anticoagulants. Here, we report the isolation and characterization of a novel anticoagulant ringhalexin from the venom of Hemachatus haemachatus (African Ringhals Cobra). Amino acid sequence of the protein indicates that it belongs to the three-finger toxin family and exhibits 94% identity to an uncharacterized Neurotoxin-like protein NTL2 from Naja atra. Ringhalexin inhibited FX activation by extrinsic tenase complex with an IC50 of 123.8 ± 9.54 nM. It is a mixed-type inhibitor with the kinetic constants, Ki and Ki' of 84.25 ± 3.53 nM and 152.5 ± 11.32 nM, respectively. Ringhalexin also exhibits a weak, irreversible neurotoxicity on chick biventer cervicis muscle preparations. Subsequently, the three-dimensional structure of ringhalexin was determined at 2.95 Å resolution. This study for the first time reports the structure of an anticoagulant three-finger toxin. Thus, ringhalexin is a potent inhibitor of the FX activation by extrinsic tenase complex and a weak, irreversible neurotoxin.

An analysis of venom ontogeny and prey-specific toxicity in the Monocled Cobra (Naja kaouthia)

Venoms of snakes of the family Elapidae (cobras, kraits, mambas, and relatives) are predominantly composed of numerous phospholipases A2 (PLA2s) and three-finger toxins (3FTxs), some of which are lethal while others are not significantly toxic. Currently, the only identified prey-specific toxins are several nonconventional 3FTxs, and given the large diversity of 3FTxs within Monocled Cobra (Naja kaouthia) venom, it was hypothesized that several 3FTxs, previously found to be non-toxic or weakly toxic 3FTxs in murine models, could potentially be toxic towards non-murine prey. Additionally, it was hypothesized that ontogenetic dietary shifts will be correlated with observable changes in specific 3FTx isoform abundance. Adult and juvenile N. kaouthia venom composition was investigated using ion-exchange FPLC, 1D and 2D SDS-PAGE, mass spectrometry, and various enzymatic and LD50 assays. Alpha-cobratoxin (α-elapitoxin) was the only significantly toxic (LD50 < 1 μg/g) 3FTx found in N. kaouthia venom and was equally toxic toward both lizard and mouse models. The abundance and diversity of 3FTxs and most enzyme activities did not vary between adult and juvenile cobra venoms; however, total venom PLA2 activity and specific PLA2 isoforms did vary, with juveniles lacking several of the least acidic PLA2s, and these differences could have both biological (related to predation) and clinical (antivenom efficacy) implications. Nevertheless, the ubiquitous presence of α-cobratoxin in both adult and juvenile cobra venoms, with high toxicity toward both reptiles and mammals, represents a venom compositional strategy wherein a single potent toxin effectively immobilizes a variety of prey types encountered across life history stages.

Proteomic characterization and comparison of venoms from two elapid snakes (Bungarus multicinctus and Naja atra) from China

Bungarus multicinctus (many-banded krait) and Naja atra (Chinese cobra) are widely distributed and medically important venomous snakes in China; however, their venom proteomic profiles have not been fully compared. Here, we fractionated crude venoms and analyzed them using a combination of proteomic techniques. Three-finger toxins (3-FTx) and phospholipase A2 (PLA2) were most abundant in both species, respectively accounting for 32.6% and 66.4% of total B. multicinctus venom, and 84.3% and 12.2% of total N. atra venom. Venoms from these two species contained one common protein family and six less abundant species-specific protein families. The proteomic profiles of B. multicinctus and N. atra venoms and analysis of toxicological activity in mice suggested that 3-FTx and PLA2 are the major contributors to clinical symptoms caused by envenomation. The venoms differed in enzymatic activity, likely the result of inter-specific variation in the amount of related venom components. Antivenomics assessment revealed that a small number of venom components (3-FTxs and PLA2s in B. multicinctus, and 3-FTxs in N. atra) could not be immunocaptured completely, suggesting that we should pay attention to enhancing the immune response of these components in designing commercial antivenoms for B. multicinctus and N. atra.

BIOLOGICAL SIGNIFICANCE

The proteomic profiles of venoms from two medically important snake species - B. multicinctus and N. atra - have been explored. Quantitative and qualitative differences are evident in both venoms when proteomic profiles and transcriptomic results are compared; this is a reminder that combined approaches are needed to explore the precise composition of snake venom. Two protein families (3-FTx and PLA2) of high abundance in these snake venoms are major players in the biochemical and pharmacological effects of envenomation. Elucidation of the proteomic profiles of these snake venoms is helpful in understanding composition-function relationships and will facilitate the clinical application of antivenoms.

Identification and molecular characterization of five putative toxins from the venom gland of the snake Philodryas chamissonis (Serpentes: Dipsadidae)

Philodryas chamissonis is a rear-fanged snake endemic to Chile. Its bite produces mild to moderate symptoms with proteolytic and anti-coagulant effects. Presently, the composition of the venom, as well as, the biochemical and structural characteristics of its toxins, remains unknown. In this study, we cloned and reported the first full-length sequences of five toxin-encoding genes from the venom gland of this species: Type III snake venom metalloprotease (SVMP), snake venom serine protease (SVSP), Cysteine-rich secretory protein (CRISP), α and β subunits of C-type lectin-like protein (CLP) and C-type natriuretic peptide (NP). These genes are highly expressed in the venom gland and their sequences exhibited a putative signal peptide, suggesting that these are components of the venom. These putative toxins had different evolutionary relationships with those reported for some front-fanged snakes, being SVMP, SVSP and CRISP of P. chamissonis closely related to the toxins present in Elapidae species, while NP was more related to those of Viperidae species. In addition, analyses suggest that the α and β subunits of CLP of P. chamissonis might have a α-subunit scaffold in common with Viperidae species, whose highly variable C-terminal region might have allowed the diversification in α and β subunits. Our results provide the first molecular description of the toxins possibly implicated in the envenomation of prey and humans by the bite of P. chamissonis.

Discovery of toxin-encoding genes from the false viper Macropisthodon rudis, a rear-fanged snake, by transcriptome analysis of venom gland

Although rear-fanged snakes are often considered as non-threatening to humans, some species are lethal or medically hazardous. The toxin components and bioactivities of front-fanged snakes have been extensively studied; however, only limited research has explored the venoms of rear-fanged snakes. The false viper, Macropisthodon rudis, is widespread in southern China, but little is known about the toxins that this snake produces. Here, we analyzed the transcriptome of the venom gland of M. rudis using high-throughput sequencing with an illumina HiSeq 2000. The raw data were assembled and annotated using public databases. The Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways and gene ontology (GO) were analyzed. Using sequence comparisons, snake venom metalloproteinases (SVMPs) and a phosphodiesterase (PDE) were discovered in the venom gland of M. rudis.

Venom-gland transcriptome and venom proteome of the Malaysian king cobra (Ophiophagus hannah)

Background

The king cobra (Ophiophagus hannah) is widely distributed throughout many parts of Asia. This study aims to investigate the complexity of Malaysian Ophiophagus hannah (MOh) venom for a better understanding of king cobra venom variation and its envenoming pathophysiology. The venom gland transcriptome was investigated using the Illumina HiSeq™ platform, while the venom proteome was profiled by 1D-SDS-PAGE-nano-ESI-LCMS/MS.

Results

Transcriptomic results reveal high redundancy of toxin transcripts (3357.36 FPKM/transcript) despite small cluster numbers, implying gene duplication and diversification within restricted protein families. Among the 23 toxin families identified, three-finger toxins (3FTxs) and snake-venom metalloproteases (SVMPs) have the most diverse isoforms. These 2 toxin families are also the most abundantly transcribed, followed in descending order by phospholipases A₂ (PLA₂s), cysteine-rich secretory proteins (CRISPs), Kunitz-type inhibitors (KUNs), and L-amino acid oxidases (LAAOs). Seventeen toxin families exhibited low mRNA expression, including hyaluronidase, DPP-IV and 5’-nucleotidase that were not previously reported in the venom-gland transcriptome of a Balinese O. hannah. On the other hand, the MOh proteome includes 3FTxs, the most abundantly expressed proteins in the venom (43 % toxin sbundance). Within this toxin family, there are 6 long-chain, 5 short-chain and 2 non-conventional 3FTx. Neurotoxins comprise the major 3FTxs in the MOh venom, consistent with rapid neuromuscular paralysis reported in systemic envenoming. The presence of toxic enzymes such as LAAOs, SVMPs and PLA₂ would explain tissue inflammation and necrotising destruction in local envenoming. Dissimilarities in the subtypes and sequences between the neurotoxins of MOh and Naja kaouthia (monocled cobra) are in agreement with the poor cross-neutralization activity of N. kaouthia antivenom used against MOh venom. Besides, the presence of cobra venom factor, nerve growth factors, phosphodiesterase, 5’-nucleotidase, and DPP-IV in the venom proteome suggests its probable hypotensive action in subduing prey.

Conclusion

This study reports the diversity and abundance of toxins in the venom of the Malaysian king cobra (MOh). The results correlate with the pathophysiological actions of MOh venom, and dispute the use of Naja cobra antivenoms to treat MOh envenomation. The findings also provide a deeper insight into venom variations due to geography, which is crucial for the development of a useful pan-regional antivenom.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1828-2) contains supplementary material, which is available to authorized users.

Purification and Characterization of Nk-3FTx: A Three Finger Toxin from the Venom of North East Indian Monocled Cobra

Snake venom three finger toxins (3FTxs) are a non-enzymatic family of venom proteins abundantly found in elapids. We have purified a 7579.5 ± 0.591 Da 3FTx named as Nk-3FTx from the venom of Naja kaouthia of North East India origin. The primary structure was determined by a combination of N-terminal sequencing and electrospray ionization liquid chromatography-mass spectrometry/mass spectrometry. Biochemical and biological characterization reveal that it is nontoxic to human cell lines and exhibit mild anticoagulant activity when tested on citrated human plasma. Nk-3FTx was found to affect the compound action potential (CAP) and nerve conduction velocity of isolated toad sciatic nerve. This is the first report of a non-conventional 3FTx from Naja kaouthia venom that reduces CAP for its neurotoxic effect. Further studies can be carried out to understand the mechanism of action and to explore its potential therapeutic application.

Cabinet of Curiosities: Venom Systems and Their Ecological Function in Mammals, with a Focus on Primates

Venom delivery systems (VDS) are common in the animal kingdom, but rare amongst mammals. New definitions of venom allow us to reconsider its diversity amongst mammals by reviewing the VDS of Chiroptera, Eulipotyphla, Monotremata, and Primates. All orders use modified anterior dentition as the venom delivery apparatus, except Monotremata, which possesses a crural system. The venom gland in most taxa is a modified submaxillary salivary gland. In Primates, the saliva is activated when combined with brachial gland exudate. In Monotremata, the crural spur contains the venom duct. Venom functions include feeding, intraspecific competition, anti-predator defense and parasite defense. Including mammals in discussion of venom evolution could prove vital in our understanding protein functioning in mammals and provide a new avenue for biomedical and therapeutic applications and drug discovery.

Natural compounds interacting with nicotinic acetylcholine receptors: from low-molecular weight ones to peptides and proteins

Nicotinic acetylcholine receptors (nAChRs) fulfill a variety of functions making identification and analysis of nAChR subtypes a challenging task. Traditional instruments for nAChR research are d-tubocurarine, snake venom protein α-bungarotoxin (α-Bgt), and α-conotoxins, neurotoxic peptides from Conus snails. Various new compounds of different structural classes also interacting with nAChRs have been recently identified. Among the low-molecular weight compounds are alkaloids pibocin, varacin and makaluvamines C and G. 6-Bromohypaphorine from the mollusk Hermissenda crassicornis does not bind to Torpedo nAChR but behaves as an agonist on human α7 nAChR. To get more selective α-conotoxins, computer modeling of their complexes with acetylcholine-binding proteins and distinct nAChRs was used. Several novel three-finger neurotoxins targeting nAChRs were described and α-Bgt inhibition of GABA-A receptors was discovered. Information on the mechanisms of nAChR interactions with the three-finger proteins of the Ly6 family was found. Snake venom phospholipases A₂ were recently found to inhibit different nAChR subtypes. Blocking of nAChRs in Lymnaea stagnalis neurons was shown for venom C-type lectin-like proteins, appearing to be the largest molecules capable to interact with the receptor. A huge nAChR molecule sensible to conformational rearrangements accommodates diverse binding sites recognizable by structurally very different compounds.

Privileged frameworks from snake venom

Venom as a form of chemical prey capture is a key innovation that has underpinned the explosive radiation of the advanced snakes (Caenophidia). Small venom proteins are often rich in disulfide bonds thus facilitating stable molecular scaffolds that present key functional residues on the protein surface. New toxin types are initially developed through the venom gland over-expression of normal body proteins, their subsequent gene duplication and diversification that leads to neofunctionalisation as random mutations modify their structure and function. This process has led to preferentially selected (privileged) cysteine-rich scaffolds that enable the snake to build arrays of toxins many of which may lead to therapeutic products and research tools. This review focuses on cysteine-rich small proteins and peptides found in snake venoms spanning natriuretic peptides to phospholipase enzymes, while highlighting their three-dimensional structures and biological functions as well as their potential as therapeutic agents or research tools.

6-bromohypaphorine from marine nudibranch mollusk Hermissenda crassicornis is an agonist of human α7 nicotinic acetylcholine receptor

6-Bromohypaphorine (6-BHP) has been isolated from the marine sponges Pachymatisma johnstoni, Aplysina sp., and the tunicate Aplidium conicum, but data on its biological activity were not available. For the nudibranch mollusk Hermissenda crassicornis no endogenous compounds were known, and here we describe the isolation of 6-BHP from this mollusk and its effects on different nicotinic acetylcholine receptors (nAChR). Two-electrode voltage-clamp experiments on the chimeric α7 nAChR (built of chicken α7 ligand-binding and glycine receptor transmembrane domains) or on rat α4β2 nAChR expressed in Xenopus oocytes revealed no action of 6-BHP. However, in radioligand analysis, 6-BHP competed with radioiodinated α-bungarotoxin for binding to human α7 nAChR expressed in GH4C1 cells (IC50 23 ± 1 μM), but showed no competition on muscle-type nAChR from Torpedo californica. In Ca2+-imaging experiments on the human α7 nAChR expressed in the Neuro2a cells, 6-BHP in the presence of PNU120596 behaved as an agonist (EC50 ~80 μM). To the best of our knowledge, 6-BHP is the first low-molecular weight compound from marine source which is an agonist of the nAChR subtype. This may have physiological importance because H. crassicornis, with its simple and tractable nervous system, is a convenient model system for studying the learning and memory processes.