Dynamic genetic differentiation drives the widespread structural and functional convergent evolution of snake venom proteinaceous toxins.

BACKGROUND: The explosive radiation and diversification of the advanced snakes (superfamily Colubroidea) was associated with changes in all aspects of the shared venom system. Morphological changes included the partitioning of the mixed ancestral glands into two discrete glands devoted for production of venom or mucous respectively, as well as changes in the location, size and structural elements of the venom-delivering teeth. Evidence also exists for homology among venom gland toxins expressed across the advanced snakes. However, despite the evolutionary novelty of snake venoms, in-depth toxin molecular evolutionary history reconstructions have been mostly limited to those types present in only two front-fanged snake families, Elapidae and Viperidae. To have a broader understanding of toxins shared among extant snakes, here we first sequenced the transcriptomes of eight taxonomically diverse rear-fanged species and four key viperid species and analysed major toxin types shared across the advanced snakes. RESULTS: Transcriptomes were constructed for the following families and species: Colubridae - Helicops leopardinus, Heterodon nasicus, Rhabdophis subminiatus; Homalopsidae - Homalopsis buccata; Lamprophiidae - Malpolon monspessulanus, Psammophis schokari, Psammophis subtaeniatus, Rhamphiophis oxyrhynchus; and Viperidae - Bitis atropos, Pseudocerastes urarachnoides, Tropidolaeumus subannulatus, Vipera transcaucasiana. These sequences were combined with those from available databases of other species in order to facilitate a robust reconstruction of the molecular evolutionary history of the key toxin classes present in the venom of the last common ancestor of the advanced snakes, and thus present across the full diversity of colubroid snake venoms. In addition to differential rates of evolution in toxin classes between the snake lineages, these analyses revealed multiple instances of previously unknown instances of structural and functional convergences. Structural convergences included: the evolution of new cysteines to form heteromeric complexes, such as within kunitz peptides (the beta-bungarotoxin trait evolving on at least two occasions) and within SVMP enzymes (the P-IIId trait evolving on at least three occasions); and the C-terminal tail evolving on two separate occasions within the C-type natriuretic peptides, to create structural and functional analogues of the ANP/BNP tailed condition. Also shown was that the de novo evolution of new post-translationally liberated toxin families within the natriuretic peptide gene propeptide region occurred on at least five occasions, with novel functions ranging from induction of hypotension to post-synaptic neurotoxicity. Functional convergences included the following: multiple occasions of SVMP neofunctionalised in procoagulant venoms into activators of the clotting factors prothrombin and Factor X; multiple instances in procoagulant venoms where kunitz peptides were neofunctionalised into inhibitors of the clot destroying enzyme plasmin, thereby prolonging the half-life of the clots formed by the clotting activating enzymatic toxins; and multiple occasions of kunitz peptides neofunctionalised into neurotoxins acting on presynaptic targets, including twice just within Bungarus venoms. CONCLUSIONS: We found novel convergences in both structural and functional evolution of snake toxins. These results provide a detailed roadmap for future work to elucidate predator-prey evolutionary arms races, ascertain differential clinical pathologies, as well as documenting rich biodiscovery resources for lead compounds in the drug design and discovery pipeline.

Ancient Diversification of Three-Finger Toxins in Micrurus Coral Snakes.

Coral snakes, most notably the genus Micrurus, are the only terrestrial elapid snakes in the Americas. Elapid venoms are generally known for their potent neurotoxicity which is usually caused by Three-Finger Toxin (3FTx) proteins. These toxins can have a wide array of functions that have been characterized from the venom of other elapids. We examined publicly available sequences from Micrurus 3FTx to show that they belong to 8 monophyletic clades that diverged as deep in the 3FTx phylogenetic tree as the other clades with characterized functions. Functional residues from previously characterized clades of 3FTx are not well conserved in most of the Micrurus toxin clades. We also analyzed the patterns of selection on these toxins and find that they have been diversifying at different rates, with some having undergone extreme diversifying selection. This suggests that Micrurus 3FTx may contain a previously underappreciated functional diversity that has implications for the clinical outcomes of bite victims, the evolution and ecology of the genus, as well as the potential for biodiscovery efforts focusing on these toxins.

Proteomic Investigation of Cape Cobra (Naja nivea) Venom Reveals First Evidence of Quaternary Protein Structures.

Naja nivea (N. nivea) is classed as a category one snake by the World Health Organization since its envenomation causes high levels of mortality and disability annually. Despite this, there has been little research into the venom composition of N. nivea, with only one full venom proteome published to date. Our current study separated N. nivea venom using size exclusion chromatography before utilizing a traditional bottom-up proteomics approach to unravel the composition of the venom proteome. As expected by its clinical presentation, N. nivea venom was found to consist mainly of neurotoxins, with three-finger toxins (3FTx), making up 76.01% of the total venom proteome. Additionally, cysteine-rich secretory proteins (CRISPs), vespryns (VESPs), cobra venom factors (CVFs), 5'-nucleotidases (5'NUCs), nerve growth factors (NGFs), phospholipase A2s (PLA2), acetylcholinesterases (AChEs), Kunitz-type serine protease inhibitor (KUN), phosphodiesterases (PDEs), L-amino acid oxidases (LAAOs), hydrolases (HYDs), snake venom metalloproteinases (SVMPs), and snake venom serine protease (SVSP) toxins were also identified in decreasing order of abundance. Interestingly, contrary to previous reports, we find PLA2 toxins in N. nivea venom. This highlights the importance of repeatedly profiling the venom of the same species to account for intra-species variation. Additionally, we report the first evidence of covalent protein complexes in N. nivea venom, which likely contribute to the potency of this venom.

Anti-platelet activity of a three-finger toxin (3FTx) from Indian monocled cobra (Naja kaouthia) venom.

A low molecular weight anti-platelet peptide (6.9 kDa) has been purified from Naja kaouthia venom and was named KT-6.9. MALDI-TOF/TOF mass spectrometry analysis revealed the homology of KT-6.9 peptide sequence with many three finger toxin family members. KT-6.9 inhibited human platelet aggregation process in a dose dependent manner. It has inhibited ADP, thrombin and arachidonic acid induced platelet aggregation process in dose dependent manner, but did not inhibit collagen and ristocetin induced platelet aggregation. Strong inhibition (70%) of the ADP induced platelet aggregation by KT-6.9 suggests competition with ADP for its receptors on platelet surface. Anti-platelet activity of KT-6.9 was found to be 25 times stronger than that of anti-platelet drug clopidogrel. Binding of KT-6.9 to platelet surface was confirmed by surface plasma resonance analysis using BIAcore X100. Binding was also observed by a modified sandwich ELISA method using anti-KT-6.9 antibodies. KT-6.9 is probably the first 3 FTx from Indian monocled cobra venom reported as a platelet aggregation inhibitor.

Biochemical and computational approaches to understand venom toxin-toxin interaction.

Snake venoms are mainly composed of proteins and peptides (venom toxins). The venom transcriptomes and proteomes have been extensively investigated; however, venom toxin-toxin interactions remain poorly characterized. We detected the interaction of venom Asp49-PLA2 and 3FTx using biochemical and computational approaches. A stable structure of Asp49-PLA2-3FTx was identified, and the interface of Asp49-PLA2 and 3FTx was analyzed. The approaches will shed light on understanding the venom complexity and deciphering the synergistic effects of venom toxins.

Venomics of the Duvernoy's gland secretion of the false coral snake Rhinobothryum bovallii (Andersson, 1916) and assessment of venom lethality towards synapsid and diapsid animal models.

The Duvernoy's gland secretory proteome of the false coral snake Rhinobothryum bovallii (Costa Rica), unveiled applying bottom-up venomics, comprises a handful of toxins belonging to only three protein families, three-finger toxin (3FTx), cysteine-rich secretory protein (CRISP), and snake venom metalloprotease (PIII-SVMP). Except for small differences in the relative abundance of the PIII-SVMPs, which may be due to individual variability, no evidence of geographic variability or ontogenetic changes was found among the venom proteomes of the juvenile and adult R. bovallii snakes sampled. Major monomeric (86.5%mol) and minor dimeric (2.8%mol) 3FTxs dominate the toxin arsenal of the Costa Rican false coral snake. The remaining 10.7% of the venom proteome comprises CRISP (8.2%) and PIII-SVMP (2.4%) molecules. In vivo lethality assays showed that R. bovallii produces venom that is non-toxic towards mammalian prey, and which exerts a different toxic effect on domestic chicken chicks and baby green iguana. Toxicovenomic analysis of R. bovallii venom in the iguana model identified two 3FTx RP-HPLC fractions that faithfully mimicked the irreversible immobilizing effect of the whole venom. BIOLOGICAL SIGNIFICANCE: With more than 2200 species in family Colubridae (sensu lato), rear-fanged snakes comprise approximately two-thirds of the extant species of advanced snakes. Snakebites from venomous snakes that are of medical concern are predominantly from front-fanged snakes of families Viperidae and Elapidae. On the other hand, rear-fanged snakes have been conventionally considered non-venomous, and thus their venoms have remained a largely untapped area of venomics. However, increasing documentation of life-threadening, even fatal, envenomings from rear-fanged snakes has sparked interest in their venoms. Appying bottom-up venomics we have revealed that the Duvernoy's gland secretory proteome of R. bovallii comprises a handful of toxins belonging to only three protein families, with slow-acting three-finger toxins (3FTx) that are non-toxic towards mammalian prey and show preference towards diapsid taxa representing the dominant structural and functional proteins. Our work documents for the first time 3FTxs exerting different effect in an avian model than in a reptile model. Besides, the 3FTx fractions that faithfully mimicked the irreversible iguana-immobilizing effect of the whole venom were identified through toxicovenomic analysis of R. bovallii venom on Iguana iguana. Our work underscores the importance of using biologically-relevant animal toxicity models for investigating the biological roles of venoms in an evolutionary-ecological context.

An Investigation of Three-Finger Toxin-nAChR Interactions through Rosetta Protein Docking.

Three-finger toxins (3FTX) are a group of peptides that affect multiple receptor types. One group of proteins affected by 3FTX are nicotinic acetylcholine receptors (nAChR). Structural information on how neurotoxins interact with nAChR is limited and is confined to a small group of neurotoxins. Therefore, in silico methods are valuable in understanding the interactions between 3FTX and different nAChR subtypes, but there are no established protocols to model 3FTX-nAChR interactions. We followed a homology modeling and protein docking protocol to address this issue and tested its success on three different systems. First, neurotoxin peptides co-crystallized with acetylcholine binding protein (AChBP) were re-docked to assess whether Rosetta protein-protein docking can reproduce the native poses. Second, experimental data on peptide binding to AChBP was used to test whether the docking protocol can qualitatively distinguish AChBP-binders from non-binders. Finally, we docked eight peptides with known α7 and muscle-type nAChR binding properties to test whether the protocol can explain the differential activities of the peptides at the two receptor subtypes. Overall, the docking protocol predicted the qualitative and some specific aspects of 3FTX binding to nAChR with reasonable success and shed light on unknown aspects of 3FTX binding to different receptor subtypes.

A Decoy-Receptor Approach Using Nicotinic Acetylcholine Receptor Mimics Reveals Their Potential as Novel Therapeutics Against Neurotoxic Snakebite.

Snakebite is a neglected tropical disease that causes 138,000 deaths each year. Neurotoxic snake venoms contain small neurotoxins, including three-finger toxins (3FTxs), which can cause rapid paralysis in snakebite victims by blocking postsynaptic transmission via nicotinic acetylcholine receptors (nAChRs). These toxins are typically weakly immunogenic and thus are often not effectively targeted by current polyclonal antivenom therapies. We investigated whether nAChR mimics, also known as acetylcholine binding proteins (AChBPs), could effectively capture 3FTxs and therefore be developed as a novel class of snake-generic therapeutics for combatting neurotoxic envenoming. First, we identified the binding specificities of 3FTx from various medically important elapid snake venoms to nAChR using two recombinant nAChR mimics: the AChBP from Lymnaea stagnalis and a humanized neuronal α7 version (α7-AChBP). We next characterized these AChBP-bound and unbound fractions using SDS-PAGE and mass spectrometry. Interestingly, both mimics effectively captured long-chain 3FTxs from multiple snake species but largely failed to capture the highly related short-chain 3FTxs, suggesting a high level of binding specificity. We next investigated whether nAChR mimics could be used as snakebite therapeutics. We showed that while α7-AChBP alone did not protect against Naja haje (Egyptian cobra) venom lethality in vivo, it significantly prolonged survival times when coadministered with a nonprotective dose of antivenom. Thus, nAChR mimics are capable of neutralizing specific venom toxins and may be useful adjunct therapeutics for improving the safety and affordability of existing snakebite treatments by reducing therapeutic doses. Our findings justify exploring the future development of AChBPs as potential snakebite treatments.

New insights into the phylogeographic distribution of the 3FTx/PLA2 venom dichotomy across genus Micrurus in South America.

Micrurus is a monophyletic genus of venomous coral snakes of the family Elapidae. The ~80 recognized species within this genus are endemic to the Americas, and are distributed from southeastern United States to northern Argentina. Although relatively few bites are recorded due to their reclusive nature, semi-fossorial habits, and their occurrence in sparsely populated areas, coral snakes possess powerful venoms that target the cholinergic system and, if early treatment is missed, can cause neuromuscular paralysis, respiratory failure, and death by asphyxiation within hours of envenoming. The to-date proteomically characterized 18 micrurine venoms exhibit a puzzling phenotypic dichotomy, characterized by the toxin arsenal being dominated either by pre-synaptically acting PLA2s or post-synaptic 3FTxs, and a general, but imperfect, distributional pattern of these venom phenotypes along the North-South axis of the American continent. The lack of perfect phylogenetic clustering suggests that phylogeny may not be the sole factor driving the evolution of the divergent venom phenotypes across Micrurus venoms. To shed new light on the origin and expression pattern of the 3FTx/PLA2 venom dichotomy, we have conducted a comparative proteomics analysis of venoms from the Brazilian ribbon coral snake, Micrurus lemniscatus carvalhoi, sourced from different localities in the Brazilian states of São Paulo; the Caatinga coral snake, M. ibiboboca, from central Bahia state (Brazil); two Micrurus specimens of uncertain taxonomy collected in the Brazilian states of Alagoas and Rio de Janeiro; and the Western ribbon coral snake, M. l. helleri, from Leticia, the southernmost town of the Colombian Department of Amazonas. Venoms from São Paulo and Rio de Janeiro showed 3FTx-predominant phenotypes, while in venoms from Leticia, Alagoas and Bahia PLA2s represented the major toxin family. Comparative venom proteomics suggests that both Micrurus venom phenotypes exhibit a high degree of toxin evolvability. Mapping the 3FTx/PLA2 dichotomy across the Americas points to a phylogeographic pattern for venom phenotypes consistent with, but more complex than, the North-South distribution hypothesis anticipated in previous investigations. BIOLOGICAL SIGNIFICANCE: New World coral snakes (Micrurus: Elapidae) produce potent venoms that target pre- and post-synaptically cholinergic nerve terminals resulting in neuromuscular paralysis, and in severe envenomings, may lead to death from asphyxiation by respiratory arrest. Presynaptic ß-neurotoxins of group IA PLA2 protein subfamily and postsynaptic α-neurotoxins with 3FTx fold are the major components (>80%) of coral snake venoms. Micrurine venoms exhibit a puzzling phenotypic venom dichotomy, characterized by the dominant expression of either α- or ß-neurotoxins. The distribution of these alternative compositional profiles has been fragmentarily studied both across Micrurus phylogeny and along the North-South axis of the genus radiation in the American continent, from southern United States to Northern Argentina. The unpredictability of the neurotoxin profile across the distribution range of the coral snakes represents a difficulty for applying the most appropriate treatment upon a coral snakebite. A deep knowledge of the phylogeographic distribution and the evolution of dichotomic Micrurus venoms would be useful for tracing the evolutionary path to their present day phenotypes, rationalizing the patchy cross-reactivity of current Micrurus antivenoms, and improving the efficacy of antivenoms to neutralize coral snake envenomings.

Comparative venomics of Brazilian coral snakes: Micrurus frontalis, Micrurus spixii spixii, and Micrurus surinamensis.

A comparative venom proteomic analysis of the Brazilian southern coral snake, M. frontalis, the Amazon coral snake M. spixii spixii, and the aquatic coral snake M. surinamensis is reported. Venoms from M. frontalis and M. s. spixii were composed mainly (>90% of the total venom proteome) by 3FTxs and PLA2s in different proportions, and minor proteins from 2 to 5 protein families. Conversely, the aquatic coral snake expressed a streamlined (95%) 3FTx venom with low abundance (4.2%) of PLA2 molecules. A compositional-lethal activity for natural prey correlation analysis suggests that M. surinamensis venom may has evolved under strong pressure to quickly immobilize aquatic prey. On the other hand, venoms from M. frontalis and M. s. spixii, whose diet consist mainly of amphisbaenians and colubrid snakes, may have been shaped through balancing selection. Our work provides strong evidence for the occurrence in M. frontalis venom, but not in those from M. s. spixi and M. surinamensis, of a KUN-PLA2 complex homologue to heterodimeric venom toxins from some long-tailed monadal coral snakes that target acid-sensing receptors ASIC1a/2 evoking pain. The M. frontalis protein would represent the first example of a KUN-PLA2 heterodimer in a South American short-tailed triadal coral snake venom.

Coralsnake Venomics: Analyses of Venom Gland Transcriptomes and Proteomes of Six Brazilian Taxa.

Venom gland transcriptomes and proteomes of six Micrurus taxa (M. corallinus, M. lemniscatus carvalhoi, M. lemniscatus lemniscatus, M. paraensis, M. spixii spixii, and M. surinamensis) were investigated, providing the most comprehensive, quantitative data on Micrurus venom composition to date, and more than tripling the number of Micrurus venom protein sequences previously available. The six venomes differ dramatically. All are dominated by 2-6 toxin classes that account for 91-99% of the toxin transcripts. The M. s. spixii venome is compositionally the simplest. In it, three-finger toxins (3FTxs) and phospholipases A2 (PLA2s) comprise >99% of the toxin transcripts, which include only four additional toxin families at levels ≥0.1%. Micrurus l. lemniscatus venom is the most complex, with at least 17 toxin families. However, in each venome, multiple structural subclasses of 3FTXs and PLA2s are present. These almost certainly differ in pharmacology as well. All venoms also contain phospholipase B and vascular endothelial growth factors. Minor components (0.1-2.0%) are found in all venoms except that of M. s. spixii. Other toxin families are present in all six venoms at trace levels (<0.005%). Minor and trace venom components differ in each venom. Numerous novel toxin chemistries include 3FTxs with previously unknown 8- and 10-cysteine arrangements, resulting in new 3D structures and target specificities. 9-cysteine toxins raise the possibility of covalent, homodimeric 3FTxs or heterodimeric toxins with unknown pharmacologies. Probable muscarinic sequences may be reptile-specific homologs that promote hypotension via vascular mAChRs. The first complete sequences are presented for 3FTxs putatively responsible for liberating glutamate from rat brain synaptosomes. Micrurus C-type lectin-like proteins may have 6-9 cysteine residues and may be monomers, or homo- or heterodimers of unknown pharmacology. Novel KSPIs, 3× longer than any seen previously, appear to have arisen in three species by gene duplication and fusion. Four species have transcripts homologous to the nociceptive toxin, (MitTx) α-subunit, but all six species had homologs to the ß-subunit. The first non-neurotoxic, non-catalytic elapid phospholipase A2s are reported. All are probably myonecrotic. Phylogenetic analysis indicates that the six taxa diverged 15-35 million years ago and that they split from their last common ancestor with Old World elapines nearly 55 million years ago. Given their early diversification, many cryptic micrurine taxa are anticipated.

Identification and characterization of B-cell epitopes of 3FTx and PLA(2) toxins from Micrurus corallinus snake venom.

The main goal of this work was to develop a strategy to identify B-cell epitopes on four different three finger toxins (3FTX) and one phospholipase A2 (PLA2) from Micrurus corallinus snake venom. 3FTx and PLA2 are highly abundant components in Elapidic venoms and are the major responsibles for the toxicity observed in envenomation by coral snakes. Overlapping peptides from the sequence of each toxin were prepared by SPOT method and three different anti-elapidic sera were used to map the epitopes. After immunogenicity analysis of the spot-reactive peptides by EPITOPIA, a computational method, nine sequences from the five toxins were chemically synthesized and antigenically and immunogenically characterized. All the peptides were used together as immunogens in rabbits, delivered with Freund's adjuvant for a first cycle of immunization and Montanide in the second. A good antibody response against individual synthetic peptides and M. corallinus venom was achieved. Anti-peptide IgGs were also cross-reactive against Micrurus frontalis and Micrurus lemniscatus crude venoms. In addition, anti-peptide IgGs inhibits the lethal and phospholipasic activities of M. corallinus crude venom. Our results provide a rational basis to the identification of neutralizing epitopes on coral snake toxins and show that their corresponding synthetic peptides could improve the generation of immuno-therapeutics. The use of synthetic peptide for immunization is a reasonable approach, since it enables poly-specificity, low risk of toxic effects and large scale production.

Functional proteomic approach to discover geographic variations of king cobra venoms from Southeast Asia and China.

This study deciphers the geographic variations of king cobra (Ophiophagus hannah) venom using functional proteomics. Pooled samples of king cobra venom (abbreviated as Ohv) were obtained from Indonesia, Malaysia, Thailand, and two provinces of China, namely Guangxi and Hainan. Using two animal models to test and compare the lethal effects, we found that the Chinese Ohvs were more fatal to mice, while the Southeast Asian Ohvs were more fatal to lizards (Eutropis multifasciata). Various phospholipases A2 (PLA2s), three-finger toxins (3FTxs) and Kunitz-type inhibitors were purified from these Ohvs and compared. Besides the two Chinese Ohv PLA2s with known sequences, eight novel PLA2s were identified from the five Ohv samples and their antiplatelet activities were compared. While two 3FTxs (namely oh-55 and oh-27) were common in all the Ohvs, different sets of 3FTx markers were present in the Chinese and Southeast Asian Ohvs. All the Ohvs contain the Kunitz inhibitor, OH-TCI, while only the Chinese Ohvs contain the inhibitor variant, Oh11-1. Relative to the Chinese Ohvs which contained more phospholipases, the Southeast Asian Ohvs had higher metalloproteinase, acetylcholine esterase, and alkaline phosphatase activities. BIOLOGICAL SIGNIFICANCE: Remarkable variations in five king cobra geographic samples reveal fast evolution and dynamic translational regulation of the venom which probably adapted to different prey ecology as testified by the lethal tests on mice and lizards. Our results predict possible variations of the king cobra envenoming to human and the importance of using local antivenin for snakebite treatment.