Three-fingered RAVERs: Rapid Accumulation of Variations in Exposed Residues of snake venom toxins

Selection Across the Three-Dimensional Structure of Venom Proteins from North American Scolopendromorph Centipedes

Gene duplication followed by nucleotide differentiation is one of the simplest mechanisms to develop new functions for genes. However, the evolutionary processes underlying the divergence of multigene families remain controversial. We used multigene families found within the diversity of toxic proteins in centipede venom to test two hypotheses related to venom evolution: the two-speed mode of venom evolution and the rapid accumulation of variation in exposed residues (RAVER) model. The two-speed mode of venom evolution proposes that different types of selection impact ancient and younger venomous lineages with negative selection being the predominant form in ancient lineages and positive selection being the dominant form in younger lineages. The RAVER hypothesis proposes that, instead of different types of selection acting on different ages of venomous lineages, the different types of selection will selectively contribute to amino acid variation based on whether the residue is exposed to the solvent where it can potentially interact directly with toxin targets. This hypothesis parallels the longstanding understanding of protein evolution that suggests that residues found within the structural or active regions of the protein will be under negative or purifying selection, and residues that do not form part of these areas will be more prone to positive selection. To test these two hypotheses, we compared the venom of 26 centipedes from the order Scolopendromorpha from six currently recognized species from across North America using both transcriptomics and proteomics. We first estimated their phylogenetic relationships and uncovered paraphyly among the genus Scolopendra and evidence for cryptic diversity among currently recognized species. Using our phylogeny, we then characterized the diverse venom components from across the identified clades using a combination of transcriptomics and proteomics. We conducted selection-based analyses in the context of predicted three-dimensional properties of the venom proteins and found support for both hypotheses. Consistent with the two-speed hypothesis, we found a prevalence of negative selection across all proteins. Consistent with the RAVER hypothesis, we found evidence of positive selection on solvent-exposed residues, with structural and less-exposed residues showing stronger signal for negative selection. Through the use of phylogenetics, transcriptomics, proteomics, and selection-based analyses, we were able to describe the evolution of venom from an ancient venomous lineage and support principles of protein evolution that directly relate to multigene family evolution.

Using organ-on-a-chip technology to study haemorrhagic activities of snake venoms on endothelial tubules

Snakebite envenomation is a major public health issue which causes severe morbidity and mortality, affecting millions of people annually. Of a diverse range of clinical manifestations, local and systemic haemorrhage are of particular relevance, as this may result in ischemia, organ failure and even cardiovascular shock. Thus far, in vitro studies have failed to recapitulate the haemorrhagic effects observed in vivo. Here, we present an organ-on-a-chip approach to investigate the effects of four different snake venoms on a perfused microfluidic blood vessel model. We assess the effect of the venoms of four snake species on epithelial barrier function, cell viability, and contraction/delamination. Our findings reveal two different mechanisms by which the microvasculature is being affected, either by disruption of the endothelial cell membrane or by delamination of the endothelial cell monolayer from its matrix. The use of our blood vessel model may shed light on the key mechanisms by which tissue-damaging venoms exert their effects on the capillary vessels, which could be helpful for the development of effective treatments against snakebites.

Specific Amino Acid Residues in the Three Loops of Snake Cytotoxins Determine Their Membrane Activity and Provide a Rationale for a New Classification of These Toxins

Cytotoxins (CTs) are three-finger membrane-active toxins present mainly in cobra venom. Our analysis of the available CT amino acid sequences, literature data on their membrane activity, and conformational equilibria in aqueous solution and detergent micelles allowed us to identify specific amino acid residues which interfere with CT incorporation into membranes. They include Pro9, Ser28, and Asn/Asp45 within the N-terminal, central, and C-terminal loops, respectively. There is a hierarchy in the effect of these residues on membrane activity: Pro9 > Ser28 > Asn/Asp45. Taking into account all the possible combinations of special residues, we propose to divide CTs into eight groups. Group 1 includes toxins containing all of the above residues. Their representatives demonstrated the lowest membrane activity. Group 8 combines CTs that lack these residues. For the toxins from this group, the greatest membrane activity was observed. We predict that when solely membrane activity determines the cytotoxic effects, the activity of CTs from a group with a higher number should exceed that of CTs from a group with a lower number. This classification is supported by the available data on the cytotoxicity and membranotropic properties of CTs. We hypothesize that the special amino acid residues within the loops of the CT molecule may indicate their involvement in the interaction with non-lipid targets.

Development of a membrane-disruption assay using phospholipid vesicles as a proxy for the detection of cellular membrane degradation

Snakebite envenoming is a global health issue that affects millions of people worldwide, and that causes morbidity rates surpassing 450,000 individuals annually. Patients suffering from snakebite morbidities may experience permanent disabilities such as pain, blindness and amputations. The (local) tissue damage that causes these life-long morbidities is the result of cell- and tissue-damaging toxins present in the venoms. These compounds belong to a variety of toxin classes and may affect cells in various ways, for example, by affecting the cell membrane. In this study, we have developed a high-throughput in vitro assay that can be used to study membrane disruption caused by snake venoms using phospholipid vesicles from egg yolk as a substrate. Resuspended chicken egg yolk was used to form these vesicles, which were fluorescently stained to allow monitoring of the degradation of egg yolk vesicles on a plate reader. The assay proved to be suitable for studying phospholipid vesicle degradation of crude venoms and was also tested for its applicability for neutralisation studies of varespladib, which is a PLA2 inhibitor. We additionally made an effort to identify the responsible toxins using liquid chromatography, followed by post-column bioassaying and protein identification using high-throughput venomics. We successfully identified various toxins in the venoms of C. rhodostoma and N. mossambica, which are likely to be involved in the observed vesicle-degrading effect. This indicates that the assay can be used for screening the membrane degrading activity of both crude and fractionated venoms as well as for neutralisation studies.

Assessing Target Specificity of the Small Molecule Inhibitor MARIMASTAT to Snake Venom Toxins: A Novel Application of Thermal Proteome Profiling

New treatments that circumvent the pitfalls of traditional antivenom therapies are critical to address the problem of snakebite globally. Numerous snake venom toxin inhibitors have shown promising cross-species neutralization of medically significant venom toxins in vivo and in vitro. The development of high-throughput approaches for the screening of such inhibitors could accelerate their identification, testing, and implementation and thus holds exciting potential for improving the treatments and outcomes of snakebite envenomation worldwide. Energetics-based proteomic approaches, including thermal proteome profiling and proteome integral solubility alteration (PISA) assays, represent "deep proteomics" methods for high throughput, proteome-wide identification of drug targets and ligands. In the following study, we apply thermal proteome profiling and PISA methods to characterize the interactions between venom toxin proteoforms in Crotalus atrox (Western Diamondback Rattlesnake) and the snake venom metalloprotease (SVMP) inhibitor marimastat. We investigate its venom proteome-wide effects and characterize its interactions with specific SVMP proteoforms, as well as its potential targeting of non-SVMP venom toxin families. We also compare the performance of PISA thermal window and soluble supernatant with insoluble precipitate using two inhibitor concentrations, providing the first demonstration of the utility of a sensitive high-throughput PISA-based approach to assess the direct targets of small molecule inhibitors for snake venom.

Synthetic development of a broadly neutralizing antibody against snake venom long-chain α-neurotoxins

Snakebite envenoming is a major global public health concern for which improved therapies are urgently needed. The antigenic diversity present in snake venom toxins from various species presents a considerable challenge to the development of a universal antivenom. Here, we used a synthetic human antibody library to find and develop an antibody that neutralizes long-chain three-finger α-neurotoxins produced by numerous medically relevant snakes. Our antibody bound diverse toxin variants with high affinity, blocked toxin binding to the nicotinic acetylcholine receptor in vitro, and protected mice from lethal venom challenge. Structural analysis of the antibody-toxin complex revealed a binding mode that mimics the receptor-toxin interaction. The overall workflow presented is generalizable for the development of antibodies that target conserved epitopes among antigenically diverse targets, and it offers a promising framework for the creation of a monoclonal antibody-based universal antivenom to treat snakebite envenoming.

A Review of Rattlesnake Venoms

Venom components are invaluable in biomedical research owing to their specificity and potency. Many of these components exist in two genera of rattlesnakes, Crotalus and Sistrurus, with high toxicity and proteolytic activity variation. This review focuses on venom components within rattlesnakes, and offers a comparison and itemized list of factors dictating venom composition, as well as presenting their known characteristics, activities, and significant applications in biosciences. There are 64 families and subfamilies of proteins present in Crotalus and Sistrurus venom. Snake venom serine proteases (SVSP), snake venom metalloproteases (SVMP), and phospholipases A2 (PLA2) are the standard components in Crotalus and Sistrurus venom. Through this review, we highlight gaps in the knowledge of rattlesnake venom; there needs to be more information on the venom composition of three Crotalus species and one Sistrurus subspecies. We discuss the activity and importance of both major and minor components in biomedical research and drug development.

Evolution of Three-Finger Toxin Genes in Neotropical Colubrine Snakes (Colubridae)

Snake venom research has historically focused on front-fanged species (Viperidae and Elapidae), limiting our knowledge of venom evolution in rear-fanged snakes across their ecologically diverse phylogeny. Three-finger toxins (3FTxs) are a known neurotoxic component in the venoms of some rear-fanged snakes (Colubridae: Colubrinae), but it is unclear how prevalent 3FTxs are both in expression within venom glands and more broadly among colubrine species. Here, we used a transcriptomic approach to characterize the venom expression profiles of four species of colubrine snakes from the Neotropics that were dominated by 3FTx expression (in the genera Chironius, Oxybelis, Rhinobothryum, and Spilotes). By reconstructing the gene trees of 3FTxs, we found evidence of putative novel heterodimers in the sequences of Chironius multiventris and Oxybelis aeneus, revealing an instance of parallel evolution of this structural change in 3FTxs among rear-fanged colubrine snakes. We also found positive selection at sites within structural loops or "fingers" of 3FTxs, indicating these areas may be key binding sites that interact with prey target molecules. Overall, our results highlight the importance of exploring the venoms of understudied species in reconstructing the full evolutionary history of toxins across the tree of life.

Domain loss enabled evolution of novel functions in the snake three-finger toxin gene superfamily

Three-finger toxins (3FTXs) are a functionally diverse family of toxins, apparently unique to venoms of caenophidian snakes. Although the ancestral function of 3FTXs is antagonism of nicotinic acetylcholine receptors, redundancy conferred by the accumulation of duplicate genes has facilitated extensive neofunctionalization, such that derived members of the family interact with a range of targets. 3FTXs are members of the LY6/UPAR family, but their non-toxin ancestor remains unknown. Combining traditional phylogenetic approaches, manual synteny analysis, and machine learning techniques (including AlphaFold2 and ProtT5), we have reconstructed a detailed evolutionary history of 3FTXs. We identify their immediate ancestor as a non-secretory LY6, unique to squamate reptiles, and propose that changes in molecular ecology resulting from loss of a membrane-anchoring domain and changes in gene expression, paved the way for the evolution of one of the most important families of snake toxins.

Development of a high-throughput in vitro screening method for the assessment of cell-damaging activities of snake venoms

Snakebite envenoming is a globally important public health issue that has devastating consequences on human health and well-being, with annual mortality rates between 81,000 and 138,000. Snake venoms may cause different pathological effects by altering normal physiological processes such as nervous transfer and blood coagulation. In addition, snake venoms can cause severe (local) tissue damage that may result in life-long morbidities, with current estimates pointing towards an additional 450,000 individuals that suffer from permanent disabilities such as amputations, contractions and blindness. Despite such high morbidity rates, research to date has been mainly focusing on neurotoxic and haemotoxic effects of snake venoms and considerably less on venom-induced tissue damage. The molecular mechanisms underlaying this pathology include membrane disruption and extracellular matrix degradation. This research describes methods used to study the (molecular) mechanisms underlaying venom-induced cell- and tissue damage. A selection of cellular bioassays and fluorescent microscopy were used to study cell-damaging activities of snake venoms in multi-well plates, using both crude and fractionated venoms. A panel of 10 representative medically relevant snake species was used, which cover a large part of the geographical regions most heavily affected by snakebite. The study comprises both morphological data as well as quantitative data on cell metabolism and viability, which were measured over time. Based on this data, a distinction could be made in the ways by which viper and elapid venoms exert their effects on cells. We further made an effort to characterise the bioactive compounds causing these effects, using a combination of liquid chromatography methods followed by bioassaying and protein identification using proteomics. The outcomes of this study might prove valuable for better understanding venom-induced cell- and tissue-damaging pathologies and could be used in the process of developing and improving snakebite treatments.

Medically important snakes and snakebite envenoming in Iran

Snakebite is a common health condition in Iran with a diverse snake fauna, especially in tropical southern and mountainous western areas of the country with plethora of snake species. The list of medically important snakes, circumstances and effects of their bite, and necessary medical care require critical appraisal and should be updated regularly. This study aims to review and map the distributions of medically important snake species of Iran, re-evaluate their taxonomy, review their venomics, describe the clinical effects of envenoming, and discuss medical management and treatment, including the use of antivenom. Nearly 350 published articles and 26 textbooks with information on venomous and mildly venomous snake species and snakebites of Iran, were reviewed, many in Persian (Farsi) language, making them relatively inaccessible to an international readership. This has resulted in a revised updated list of Iran's medically important snake species, with taxonomic revisions of some, compilation of their morphological features, remapping of their geographical distributions, and description of species-specific clinical effects of envenoming. Moreover, the antivenom manufactured in Iran is discussed, together with treatment protocols that have been developed for the hospital management of envenomed patients.

The deep-rooted origin of disulfide-rich spider venom toxins

Spider venoms are a complex concoction of enzymes, polyamines, inorganic salts, and disulfide-rich peptides (DRPs). Although DRPs are widely distributed and abundant, their bevolutionary origin has remained elusive. This knowledge gap stems from the extensive molecular divergence of DRPs and a lack of sequence and structural data from diverse lineages. By evaluating DRPs under a comprehensive phylogenetic, structural and evolutionary framework, we have not only identified 78 novel spider toxin superfamilies but also provided the first evidence for their common origin. We trace the origin of these toxin superfamilies to a primordial knot - which we name 'Adi Shakti', after the creator of the Universe according to Hindu mythology - 375 MYA in the common ancestor of Araneomorphae and Mygalomorphae. As the lineages under evaluation constitute nearly 60% of extant spiders, our findings provide fascinating insights into the early evolution and diversification of the spider venom arsenal. Reliance on a single molecular toxin scaffold by nearly all spiders is in complete contrast to most other venomous animals that have recruited into their venoms diverse toxins with independent origins. By comparatively evaluating the molecular evolutionary histories of araneomorph and mygalomorph spider venom toxins, we highlight their contrasting evolutionary diversification rates. Our results also suggest that venom deployment (e.g. prey capture or self-defense) influences evolutionary diversification of DRP toxin superfamilies.

Killer Knots: Molecular Evolution of Inhibitor Cystine Knot Toxins in Wandering Spiders (Araneae: Ctenidae)

Venom expressed by the nearly 50,000 species of spiders on Earth largely remains an untapped reservoir of a diverse array of biomolecules with potential for pharmacological and agricultural applications. A large fraction of the noxious components of spider venoms are a functionally diverse family of structurally related polypeptides with an inhibitor cystine knot (ICK) motif. The cysteine-rich nature of these toxins makes structural elucidation difficult, and most studies have focused on venom components from the small handful of medically relevant spider species such as the highly aggressive Brazilian wandering spider Phoneutria nigriventer. To alleviate difficulties associated with the study of ICK toxins in spiders, we devised a comprehensive approach to explore the evolutionary patterns that have shaped ICK functional diversification using venom gland transcriptomes and proteomes from phylogenetically distinct lineages of wandering spiders and their close relatives. We identified 626 unique ICK toxins belonging to seven topological elaborations. Phylogenetic tests of episodic diversification revealed distinct regions between cysteine residues that demonstrated differential evidence of positive or negative selection, which may have structural implications towards the specificity and efficacy of these toxins. Increased taxon sampling and whole genome sequencing will provide invaluable insights to further understand the evolutionary processes that have given rise to this diverse class of toxins.

Stings on wings: Proteotranscriptomic and biochemical profiling of the lesser banded hornet (Vespa affinis) venom

Distinct animal lineages have convergently recruited venoms as weaponry for prey capture, anti-predator defence, conspecific competition, or a combination thereof. Most studies, however, have been primarily confined to a narrow taxonomic breadth. The venoms of cone snails, snakes, spiders and scorpions remain particularly well-investigated. Much less explored are the venoms of wasps (Order: Hymenoptera) that are infamous for causing excruciating and throbbing pain, justifying their apex position on Schmidt's pain index, including some that are rated four on four. For example, the lesser banded wasp (V. affinis) is clinically important yet has only been the subject of a few studies, despite being commonly found across tropical and subtropical Asia. Stings from these wasps, especially from multiple individuals of a nest, often lead to clinically severe manifestations, including mastocytosis, myasthenia gravis, optic neuropathy, and life-threatening pathologies such as myocardial infarction and organ failure. However, their venom composition and activity remain unexplored in the Indian subcontinent. Here, we report the proteomic composition, transcriptomic profile, and biochemical and pharmacological activities of V. affinis venom from southern India. Our findings suggest that wasp venoms are rich in diverse toxins that facilitate antipredator defence. Biochemical and pharmacological assessments reveal that these toxins can exhibit significantly higher activities than their homologues in medically important snakes. Their ability to exert potent effects on diverse molecular targets makes them a treasure trove for discovering life-saving therapeutics. Fascinatingly, wasp venoms, being evolutionarily ancient, exhibit a greater degree of compositional and sequence conservation across very distant populations/species, which contrasts with the patterns of venom evolution observed in evolutionarily younger lineages, such as advanced snakes and cone snails.

Membrane-mediated interaction of non-conventional snake three-finger toxins with nicotinic acetylcholine receptors

Nicotinic acetylcholine receptor of α7 type (α7-nAChR) presented in the nervous and immune systems and epithelium is a promising therapeutic target for cognitive disfunctions and cancer treatment. Weak toxin from Naja kaouthia venom (WTX) is a non-conventional three-finger neurotoxin, targeting α7-nAChR with weak affinity. There are no data on interaction mode of non-conventional neurotoxins with nAChRs. Using α-bungarotoxin (classical three-finger neurotoxin with high affinity to α7-nAChR), we showed applicability of cryo-EM to study complexes of α7-nAChR extracellular ligand-binding domain (α7-ECD) with toxins. Using cryo-EM structure of the α7-ECD/WTX complex, together with NMR data on membrane active site in the WTX molecule and mutagenesis data, we reconstruct the structure of α7-nAChR/WTX complex in the membrane environment. WTX interacts at the entrance to the orthosteric site located at the receptor intersubunit interface and simultaneously forms the contacts with the membrane surface. WTX interaction mode with α7-nAChR significantly differs from α-bungarotoxin's one, which does not contact the membrane. Our study reveals the important role of the membrane for interaction of non-conventional neurotoxins with the nicotinic receptors.

The Target Selects the Toxin: Specific Amino Acids in Snake-Prey Nicotinic Acetylcholine Receptors That Are Selectively Bound by King Cobra Venoms

Snake venom is an adaptive ecological trait that has evolved primarily as a form of prey subjugation. Thus, the selection pressure for toxin diversification is exerted by the prey’s physiological targets, with this pressure being particularly acute for specialist feeders, such as the King Cobra species, all of which are snake-prey specialists. However, while extensive research has been undertaken to elucidate key amino acids that guide toxin structure–activity relationships, reciprocal investigations into the specific sites guiding prey-lineage selective effects have been lacking. This has largely been due to the lack of assay systems amenable to systematic amino acid replacements of targeted proteins in the prey’s physiological pathways. To fill this knowledge gap, we used a recently described approach based upon mimotope peptides corresponding to the orthosteric site of nicotinic acetylcholine receptor alpha-1 subunits, a major binding site for snake venom neurotoxins that cause flaccid paralysis. We investigated the venoms of four different types of King Cobra (Cambodian, Javan, Malaysian, and Thai). This approach allowed for the determination of the key amino acid positions in King Cobra snake prey that are selectively bound by the toxins, whereby replacing these amino acids in the snake-prey orthosteric site with those from lizards or rats resulted in a significantly lower level of binding by the venoms, while conversely replacing the lizard or rat amino acids with those from the snake at that position increased the binding. By doing such, we identified three negatively charged amino acids in the snake orthosteric site that are strongly bound by the positively charged neurotoxic three-finger toxins found in King Cobra venom. This study, thus, sheds light on the selection pressures exerted by a specialist prey item for the evolution of lineage-selective toxins.

The structural and functional divergence of a neglected three-finger toxin subfamily in lethal elapids

Bungarus multicinctus is a widely distributed and medically important elapid snake that produces lethal neurotoxic venom. To study and enhance existing antivenom, we explore the complete repertoire of its toxin genes based on de novo chromosome-level assembly and multi-tissue transcriptome data. Comparative genomic analyses suggest that the three-finger toxin family (3FTX) may evolve through the neofunctionalization of flanking LY6E. A long-neglected 3FTX subfamily (i.e., MKA-3FTX) is also investigated. Only one MKA-3FTX gene, which evolves a different protein conformation, is under positive selection and actively transcribed in the venom gland, functioning as a major toxin effector together with MKT-3FTX subfamily homologs. Furthermore, this lethal snake may acquire self-resistance to its β-bungarotoxin via amino acid replacements on fast-evolving KCNA2. This study provides valuable resources for further evolutionary and structure-function studies of snake toxins, which are fundamental for the development of effective antivenoms and drug candidates.

Contextual Constraints: Dynamic Evolution of Snake Venom Phospholipase A2

Venom is a dynamic trait that has contributed to the success of numerous organismal lineages. Predominantly composed of proteins, these complex cocktails are deployed for predation and/or self-defence. Many non-toxic physiological proteins have been convergently and recurrently recruited by venomous animals into their toxin arsenal. Phospholipase A₂ (PLA₂) is one such protein and features in the venoms of many organisms across the animal kingdom, including snakes of the families Elapidae and Viperidae. Understanding the evolutionary history of this superfamily would therefore provide insight into the origin and diversification of venom toxins and the evolution of novelty more broadly. The literature is replete with studies that have identified diversifying selection as the sole influence on PLA₂ evolution. However, these studies have largely neglected the structural/functional constraints on PLA₂s, and the ecology and evolutionary histories of the diverse snake lineages that produce them. By considering these crucial factors and employing evolutionary analyses integrated with a schema for the classification of PLA₂s, we uncovered lineage-specific differences in selection regimes. Thus, our work provides novel insights into the evolution of this major snake venom toxin superfamily and underscores the importance of considering the influence of evolutionary and ecological contexts on molecular evolution.

Therapeutic potential of venom peptides: insights in the nanoparticle-mediated venom formulations

Background

Venoms are the secretions produced by animals, generally for the purpose of self-defense or catching a prey. Biochemically venoms are mainly composed of proteins, lipids, carbohydrates, ions, etc., and classified into three major classes, viz. neurotoxic, hemotoxic and cytotoxic based upon their mode of action. Venoms are composed of different specific peptides/toxins which are responsible for their unique biological actions. Though venoms are generally seen as a source of death, scientifically venom is a complex biochemical substance having a specific pharmacologic action which can be used as agents to diagnose and cure a variety of diseases in humans.

Main body

Many of these venoms have been used since centuries, and their specified therapies can also be found in ancient texts such as Charka Samhita. The modern-day example of such venom therapeutic is captopril, an antihypertensive drug developed from venom of Bothrops jararaca. Nanotechnology is a modern-day science of building materials on a nanoscale with advantages like target specificity, increased therapeutic response and diminished side effects. In the present review we have introduced the venom, sources and related constituents in brief, by highlighting the therapeutic potential of venom peptides and focusing more on the nanoformulations-based approaches. This review is an effort to compile all such report to have an idea about the future direction about the nanoplatforms which should be focused to have more clinically relevant formulations for difficult to treat diseases.

Conclusion

Venom peptides which are fatal in nature if used cautiously and effectively can save life. Several research findings suggested that many of the fatal diseases can be effectively treated with venom peptides. Nanotechnology has emerged as novel strategy in diagnosis, treatment and mitigation of diseases in more effective ways. A variety of nanoformulation approaches have been explored to enhance the therapeutic efficacy and reduce the toxicity and targeted delivery of the venom peptide conjugated with it. We concluded that venom peptides along with nanoparticles can evolve as the new era for potential treatments of ongoing and untreatable diseases.

Graphical Abstract

Taxon-selective venom variation in adult and neonate Daboia russelii (Russell's Viper), and antivenom efficacy

Major variations in venom composition can occur between juvenile and adult venomous snakes. However, due to logistical constraints, antivenoms are produced using adult venoms in immunising mixtures, possibly resulting in limited neutralisation of juvenile snake venoms. Daboia russelii is one of the leading causes of snakebite death across South Asia. Its venom is potently procoagulant, causing stroke in prey animals but causing in humans consumptive coagulopathy-a net anticoagulant state-and sometimes death resulting from hemorrhage. In this in vitro study, we compared the venom activity of-and antivenom efficacy against-six 2-week-old D. russelii relative to that of their parents. Using a coagulation analyser, we quantified the relative coagulotoxicity of these venoms in human, avian, and amphibian plasma. The overall potency on human plasma was similar across all adult and neonate venoms, and SII (Serum Institute of India) antivenom was equipotent in neutralising these coagulotoxic effects. In addition, all venoms were also similar in their action upon avian plasma. In contrast, the neonate venoms were more potent on amphibian plasma, suggesting amphibians make up a larger proportion of neonate diet than adult diet. A similar venom potency in human and avian plasmas but varying selectivity for amphibian plasma suggests ontogenetic differences in toxin isoforms within the factor X or factor V activating classes, thereby providing a testable hypothesis for future transcriptomics work. By providing insights into the functional venom differences between adult and neonate D. russelii venoms, we hope to inform clinical treatment of patients envenomated by this deadly species and to shed new light on the natural history of these extremely medically important snakes.

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.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01208-9.

Antivenom: An immunotherapy for the treatment of snakebite envenoming in sub-Saharan Africa

Snakebite envenoming (SBE) leads to significant morbidity and mortality, resulting in over 90,000 deaths and approximately 400,000 amputations annually. In sub-Saharan Africa (SSA) alone, SBE accounts for over 30,000 deaths per annum. Since 2017, SBE has been classified as a priority Neglected Tropical Disease (NTD) by the World Health Organisation (WHO). The major species responsible for mortality from SBE within SSA are from the Bitis, Dendroaspis, Echis and Naja genera. Pharmacologically active toxins such as metalloproteinases, serine proteinases, 3-finger toxins, kunitz-type toxins, and phospholipase A₂s are the primary snake venom components. These toxins induce cytotoxicity, coagulopathy, hemorrhage, and neurotoxicity in envenomed victims. Antivenom is currently the only available venom-specific treatment for SBE and contains purified equine or ovine polyclonal antibodies, collected from donor animals repeatedly immunized with low doses of adjuvanted venom. The resulting plasma or serum contains a high titre of specific antibodies, which can then be collected and stored until required. The purified antibodies are either whole IgG, monovalent fragment antibody (Fab) or divalent fragment antibody F(ab')₂. Despite pharmacokinetic and pharmacodynamic differences, all three are effective in the treatment of SBE. No antivenom is without adverse reactions but, the level of their impact and severity varies from benign early adverse reactions to the rarely occurring fatal anaphylactic shock. However, the major side effects are largely reversible with immediate administration of adrenaline and corticosteroids. There are 16 different antivenoms marketed within SSA, but the efficacy and safety profiles are only published for less than 50% of these products.

Development of Antibody Detection ELISA Based on Immunoreactive Toxins and Toxin-Derived Peptides to Evaluate the Neutralization Potency of Equine Plasma against Naja atra in Taiwan

Naja atra, also known as Taiwanese cobra, is one of the most prevalent venomous snakes in Taiwan. Clinically, freeze-dried neurotoxic antivenom (FNAV) produced from horses by Taiwan Centers for Disease Control (CDC) has been the only approved treatment for N. atra envenoming for the last few decades. During antivenom production, large numbers of mice are used in the in vivo assay to determine whether the neutralization potency of hyperimmunized equines is satisfactory for large-scale harvesting. However, this in vivo assay is extremely laborious, expensive, and significantly impairs animal welfare. In the present study, we aimed to develop an in vitro ELISA-based system that could serve as an alternative assay to evaluate the neutralization potency of plasma from hyperimmunized equines. We initially obtained 51 plasma samples with known (high or low) neutralization potency assessed in vivo from 9 hyperimmunized equines and subsequently determined their antibody titers against the five major protein components of N. atra venom (neurotoxin (NTX), phospholipase A2 (PLA₂), cytotoxin (CTX), cysteine-rich secretory protein (CRISP), and snake venom metalloproteinase (SVMP)) via ELISA. The antibody titer against NTX was the most effective in discriminating between high and low potency plasma samples. To identify the specific epitope(s) of NTX recognized by neutralization potency-related antibodies, 17 consecutive NTX-derived pentadecapeptides were synthesized and used as antigens to probe the 51 equine plasma samples. Among the 17 peptides, immunoreactive signals for three consecutive peptides (NTX1-8, NTX1-9, and NTX1-10) were significantly higher in the high potency relative to low potency equine plasma groups (p < 0.0001). Our ELISA system based on NTX1-10 peptide (RWRDHRGYRTERGCG) encompassing residues 28–42 of NTX displayed optimal sensitivity (96.88%) and specificity (89.47%) for differentiating between high- and low-potency plasma samples (area under the receiver operating characteristic curve (AUC) = 0.95). The collective data clearly indicate that the antibody titer against NTX protein or derived peptides can be used to efficiently discriminate between high and low neutralization potency of plasma samples from venom-immunized horses. This newly developed antibody detection ELISA based on NTX or its peptide derivatives has good potential to complement or replace the in vivo rodent assay for determining whether the neutralization potency of equine plasma is satisfactory for large-scale harvesting in the antivenom production process against N. atra.

The biology and evolution of spider venoms

Spiders are diverse, predatory arthropods that have inhabited Earth for around 400 million years. They are well known for their complex venom systems that are used to overpower their prey. Spider venoms contain many proteins and peptides with highly specific and potent activities suitable for biomedical or agrochemical applications, but the key role of venoms as an evolutionary innovation is often overlooked, even though this has enabled spiders to emerge as one of the most successful animal lineages. In this review, we discuss these neglected biological aspects of spider venoms. We focus on the morphology of spider venom systems, their major components, biochemical and chemical plasticity, as well as ecological and evolutionary trends. We argue that the effectiveness of spider venoms is due to their unprecedented complexity, with diverse components working synergistically to increase the overall potency. The analysis of spider venoms is difficult to standardize because they are dynamic systems, fine-tuned and modified by factors such as sex, life-history stage and biological role. Finally, we summarize the mechanisms that drive spider venom evolution and highlight the need for genome-based studies to reconstruct the evolutionary history and physiological networks of spider venom compounds with more certainty.

Mutual enlightenment: A toolbox of concepts and methods for integrating evolutionary and clinical toxinology via snake venomics and the contextual stance

Snakebite envenoming is a neglected tropical disease that may claim over 100,000 human lives annually worldwide. Snakebite occurs as the result of an interaction between a human and a snake that elicits either a defensive response from the snake or, more rarely, a feeding response as the result of mistaken identity. Snakebite envenoming is therefore a biological and, more specifically, an ecological problem. Snake venom itself is often described as a "cocktail", as it is a heterogenous mixture of molecules including the toxins (which are typically proteinaceous) responsible for the pathophysiological consequences of envenoming. The primary function of venom in snake ecology is pre-subjugation, with defensive deployment of the secretion typically considered a secondary function. The particular composition of any given venom cocktail is shaped by evolutionary forces that include phylogenetic constraints associated with the snake's lineage and adaptive responses to the snake's ecological context, including the taxa it preys upon and by which it is predated upon. In the present article, we describe how conceptual frameworks from ecology and evolutionary biology can enter into a mutually enlightening relationship with clinical toxinology by enabling the consideration of snakebite envenoming from an "ecological stance". We detail the insights that may emerge from such a perspective and highlight the ways in which the high-fidelity descriptive knowledge emerging from applications of -omics era technologies - "venomics" and "antivenomics" - can combine with evolutionary explanations to deliver a detailed understanding of this multifactorial health crisis.

The evolutionary dynamics of venom toxins made by insects and other animals

Animal venoms are recognised as unique biological systems in which to study molecular evolution. Venom use has evolved numerous times among the insects, and insects today use venom to capture prey, defend themselves from predators, or to subdue and modulate host responses during parasitism. However, little is known about most insect venom toxins or the mode and tempo by which they evolve. Here, I review the evolutionary dynamics of insect venom toxins, and argue that insects offer many opportunities to examine novel aspects of toxin evolution. The key questions addressed are: How do venomous animals evolve from non-venomous animals, and how does this path effect the composition and pharmacology of the venom? What genetic processes (gene duplication, co-option, neofunctionalisation) are most important in toxin evolution? What kinds of selection pressures are acting on toxin-encoding genes and their cognate targets in envenomated animals? The emerging evidence highlights that venom composition and pharmacology adapts quickly in response to changing selection pressures resulting from new ecological interactions, and that such evolution occurs through a stunning variety of genetic mechanisms. Insects offer many opportunities to investigate the evolutionary dynamics of venom toxins due to their evolutionary history rich in venom-related adaptations, and their quick generation time and suitability for culture in the laboratory.

A Meta-Analysis of the Protein Components in Rattlesnake Venom

The specificity and potency of venom components give them a unique advantage in developing various pharmaceutical drugs. Though venom is a cocktail of proteins, rarely are the synergy and association between various venom components studied. Understanding the relationship between various components of venom is critical in medical research. Using meta-analysis, we observed underlying patterns and associations in the appearance of the toxin families. For Crotalus, Dis has the most associations with the following toxins: PDE; BPP; CRL; CRiSP; LAAO; SVMP P-I and LAAO; SVMP P-III and LAAO. In Sistrurus venom, CTL and NGF have the most associations. These associations can predict the presence of proteins in novel venom and understand synergies between venom components for enhanced bioactivity. Using this approach, the need to revisit the classification of proteins as major components or minor components is highlighted. The revised classification of venom components is based on ubiquity, bioactivity, the number of associations, and synergies. The revised classification can be expected to trigger increased research on venom components, such as NGF, which have high biomedical significance. Using hierarchical clustering, we observed that the genera's venom compositions were similar, based on functional characteristics rather than phylogenetic relationships.

Remarkable intrapopulation venom variability in the monocellate cobra (Naja kaouthia) unveils neglected aspects of India's snakebite problem

Interpopulation venom variation has been widely documented in snakes across large geographical distances. This variability is known to markedly influence the effectiveness of snakebite therapy, as antivenoms manufactured against one population may not be effective against others. In contrast, the extent of intrapopulation venom variability, especially at finer geographical scales, remains largely uninvestigated. Moreover, given the historical focus on the 'big four' Indian snakes, our understanding of venom variation in medically important yet neglected snakes, such as the monocellate cobra (Naja kaouthia), remains unclear. To address this shortcoming, we investigated N. kaouthia venoms sampled across a small spatial scale (<50 km) in Eastern India. An interdisciplinary approach employed in this study unveiled considerable intrapopulation differences in the venom proteomic composition, pharmacological and biochemical activities, and toxicity profiles. Documentation of stark differences in venoms at such a finer geographical scale, despite the influence of similar ecological and environmental conditions, is intriguing. Furthermore, evaluation of in vitro and in vivo venom recognition and neutralisation potential of Indian polyvalent 'big four' antivenoms and Thai monovalent N. kaouthia antivenom revealed concerning deficiencies. These results highlight the negative impact of phylogenetic divergence and intrapopulation snake venom variation on the effectiveness of conventional antivenom therapy. SIGNIFICANCE: In contrast to our understanding of snake venom variation across large distances, which is theorised to be shaped by disparities in ecology and environment, intrapopulation variation at finer geographic scales remains scarcely investigated. Assessment of intrapopulation venom variability in Naja kaouthia at a small spatial scale (<50 km) in Eastern India unravelled considerable differences in venom compositions, activities and potencies. While the influence of subtle differences in prey preference and local adaptations cannot be ruled out, these findings, perhaps, also emphasise the role of accelerated molecular evolutionary regimes that rapidly introduce variations in evolutionarily younger lineages, such as advanced snakes. The inability of 'big four' Indian antivenoms and Thai N. kaouthia monovalent antivenom in countering these variations highlights the importance of phylogenetic considerations for the development of efficacious snakebite therapy. Thus, we provide valuable insights into the venoms of one of the most medically important yet neglected Indian snakes.

Production, composition, and mode of action of the painful defensive venom produced by a limacodid caterpillar, Doratifera vulnerans

Venoms have evolved independently several times in Lepidoptera. Limacodidae is a family with worldwide distribution, many of which are venomous in the larval stage, but the composition and mode of action of their venom is unknown. Here, we use imaging technologies, transcriptomics, proteomics, and functional assays to provide a holistic picture of the venom system of a limacodid caterpillar, Doratifera vulnerans Contrary to dogma that defensive venoms are simple in composition, D. vulnerans produces a complex venom containing 151 proteinaceous toxins spanning 59 families, most of which are peptides <10 kDa. Three of the most abundant families of venom peptides (vulnericins) are 1) analogs of the adipokinetic hormone/corazonin-related neuropeptide, some of which are picomolar agonists of the endogenous insect receptor; 2) linear cationic peptides derived from cecropin, an insect innate immune peptide that kills bacteria and parasites by disrupting cell membranes; and 3) disulfide-rich knottins similar to those that dominate spider venoms. Using venom fractionation and a suite of synthetic venom peptides, we demonstrate that the cecropin-like peptides are responsible for the dominant pain effect observed in mammalian in vitro and in vivo nociception assays and therefore are likely to cause pain after natural envenomations by D. vulnerans Our data reveal convergent molecular evolution between limacodids, hymenopterans, and arachnids and demonstrate that lepidopteran venoms are an untapped source of novel bioactive peptides.

Biogeographical venom variation in the Indian spectacled cobra (Naja naja) underscores the pressing need for pan-India efficacious snakebite therapy

Background

Snake venom composition is dictated by various ecological and environmental factors, and can exhibit dramatic variation across geographically disparate populations of the same species. This molecular diversity can undermine the efficacy of snakebite treatments, as antivenoms produced against venom from one population may fail to neutralise others. India is the world’s snakebite hotspot, with 58,000 fatalities and 140,000 morbidities occurring annually. Spectacled cobra (Naja naja) and Russell’s viper (Daboia russelii) are known to cause the majority of these envenomations, in part due to their near country-wide distributions. However, the impact of differing ecologies and environment on their venom compositions has not been comprehensively studied.

Methods

Here, we used a multi-disciplinary approach consisting of venom proteomics, biochemical and pharmacological analyses, and in vivo research to comparatively analyse N. naja venoms across a broad region (>6000 km; seven populations) covering India’s six distinct biogeographical zones.

Findings

By generating the most comprehensive pan-Indian proteomic and toxicity profiles to date, we unveil considerable differences in the composition, pharmacological effects and potencies of geographically-distinct venoms from this species and, through the use of immunological assays and preclinical experiments, demonstrate alarming repercussions on antivenom therapy. We find that commercially-available antivenom fails to effectively neutralise envenomations by the pan-Indian populations of N. naja, including a complete lack of neutralisation against the desert Naja population.

Conclusion

Our findings highlight the significant influence of ecology and environment on snake venom composition and potency, and stress the pressing need to innovate pan-India effective antivenoms to safeguard the lives, limbs and livelihoods of the country’s 200,000 annual snakebite victims.

Annually, India is burdened by the highest number of snake envenomations across the globe, with over 58,000 fatalities and three times the number of morbidities, predominantly affecting the rural agrarian communities. The spectacled cobra (Naja naja) and Russell’s viper (Daboia russelii) are responsible for the vast majority of envenomations in the country, in part, due to their near country-wide distributions. In this study, we unveil the astounding differences in venom composition of N. naja from six different biogeographical zones across the country (>6000 km). We provide a comprehensive account of their disparate venom proteomic profiles, biochemical and pharmacological effects, and the associated potencies. Our study uncovers alarming differences in the efficacy of the marketed polyvalent antivenoms in neutralising these venoms, thereby, emphasising the pressing need to develop dose-efficacious and pan-India effective antivenoms for the treatment of snakebites in the country. This study also highlights the significant influence of ecology and diverse environments on the venom variability, insinuating the necessity for innovating cost-effective and pan-India efficacious solutions to safeguard the lives, limbs and livelihoods of India’s two hundred thousand annual snakebite victims.

Electric Blue: Molecular Evolution of Three-Finger Toxins in the Long-Glanded Coral Snake Species Calliophis bivirgatus

The genus Calliophis is the most basal branch of the family Elapidae and several species in it have developed highly elongated venom glands. Recent research has shown that C. bivirgatus has evolved a seemingly unique toxin (calliotoxin) that produces spastic paralysis in their prey by acting on the voltage-gated sodium (NaV) channels. We assembled a transcriptome from C. bivirgatus to investigate the molecular characteristics of these toxins and the venom as a whole. We find strong confirmation that this genus produces the classic elapid eight-cysteine three-finger toxins, that δδ-elapitoxins (toxins that resemble calliotoxin) are responsible for a substantial portion of the venom composition, and that these toxins form a distinct clade within a larger, more diverse clade of C. bivirgatus three-finger toxins. This broader clade of C. bivirgatus toxins also contains the previously named maticotoxins and is somewhat closely related to cytotoxins from other elapids. However, the toxins from this clade that have been characterized are not themselves cytotoxic. No other toxins show clear relationships to toxins of known function from other species.

A Wolf in Another Wolf's Clothing: Post-Genomic Regulation Dictates Venom Profiles of Medically-Important Cryptic Kraits in India

The Common Krait (Bungarus caeruleus) shares a distribution range with many other ‘phenotypically-similar’ kraits across the Indian subcontinent. Despite several reports of fatal envenomings by other Bungarus species, commercial Indian antivenoms are only manufactured against B. caeruleus. It is, therefore, imperative to understand the distribution of genetically distinct lineages of kraits, the compositional differences in their venoms, and the consequent impact of venom variation on the (pre)clinical effectiveness of antivenom therapy. To address this knowledge gap, we conducted phylogenetic and comparative venomics investigations of kraits in Southern and Western India. Phylogenetic reconstructions using mitochondrial markers revealed a new species of krait, Romulus’ krait (Bungarus romulusi sp. nov.), in Southern India. Additionally, we found that kraits with 17 mid-body dorsal scale rows in Western India do not represent a subspecies of the Sind Krait (B. sindanus walli) as previously believed, but are genetically very similar to B. sindanus in Pakistan. Furthermore, venom proteomics and comparative transcriptomics revealed completely contrasting venom profiles. While the venom gland transcriptomes of all three species were highly similar, venom proteomes and toxicity profiles differed significantly, suggesting the prominent role of post-genomic regulatory mechanisms in shaping the venoms of these cryptic kraits. In vitro venom recognition and in vivo neutralisation experiments revealed a strong negative impact of venom variability on the preclinical performance of commercial antivenoms. While the venom of B. caeruleus was neutralised as per the manufacturer’s claim, performance against the venoms of B. sindanus and B. romulusi was poor, highlighting the need for regionally-effective antivenoms in India.

Novel Bradykinin-Potentiating Peptides and Three-Finger Toxins from Viper Venom: Combined NGS Venom Gland Transcriptomics and Quantitative Venom Proteomics of the Azemiops feae Viper

Feae's viper Azemipos feae belongs to the Azemiopinae subfamily of the Viperidae family. The effects of Viperidae venoms are mostly coagulopathic with limited neurotoxicity manifested by phospholipases A2. From A. feae venom, we have earlier isolated azemiopsin, a novel neurotoxin inhibiting the nicotinic acetylcholine receptor. To characterize other A. feae toxins, we applied label-free quantitative proteomics, which revealed 120 unique proteins, the most abundant being serine proteinases and phospholipases A2. In total, toxins representing 14 families were identified, among which bradykinin-potentiating peptides with unique amino acid sequences possessed biological activity in vivo. The proteomic analysis revealed also basal (commonly known as non-conventional) three-finger toxins belonging to the group of those possessing neurotoxic activity. This is the first indication of the presence of three-finger neurotoxins in viper venom. In parallel, the transcriptomic analysis of venom gland performed by Illumina next-generation sequencing further revealed 206 putative venom transcripts. Together, the study unveiled the venom proteome and venom gland transciptome of A. feae, which in general resemble those of other snakes from the Viperidae family. However, new toxins not found earlier in viper venom and including three-finger toxins and unusual bradykinin-potentiating peptides were discovered.

Snake three-finger α-neurotoxins and nicotinic acetylcholine receptors: molecules, mechanisms and medicine

Snake venom three-finger α-neurotoxins (α-3FNTx) act on postsynaptic nicotinic acetylcholine receptors (nAChRs) at the neuromuscular junction (NMJ) to produce skeletal muscle paralysis. The discovery of the archetypal α-bungarotoxin (α-BgTx), almost six decades ago, exponentially expanded our knowledge of membrane receptors and ion channels. This included the localisation, isolation and characterization of the first receptor (nAChR); and by extension, the pathophysiology and pharmacology of neuromuscular transmission and associated pathologies such as myasthenia gravis, as well as our understanding of the role of α-3FNTxs in snakebite envenomation leading to novel concepts of targeted treatment. Subsequent studies on a variety of animal venoms have yielded a plethora of novel toxins that have revolutionized molecular biomedicine and advanced drug discovery from bench to bedside. This review provides an overview of nAChRs and their subtypes, classification of α-3FNTxs and the challenges of typifying an increasing arsenal of structurally and functionally unique toxins, and the three-finger protein (3FP) fold in the context of the uPAR/Ly6/CD59/snake toxin superfamily. The pharmacology of snake α-3FNTxs including their mechanisms of neuromuscular blockade, variations in reversibility of nAChR interactions, specificity for nAChR subtypes or for distinct ligand-binding interfaces within a subtype and the role of α-3FNTxs in neurotoxic envenomation are also detailed. Lastly, a reconciliation of structure-function relationships between α-3FNTx and nAChRs, derived from historical mutational and biochemical studies and emerging atomic level structures of nAChR models in complex with α-3FNTxs is discussed.

Causes and Consequences of Snake Venom Variation

Snake venoms are mixtures of toxins that vary extensively between and within snake species. This variability has serious consequences for the management of the world’s 1.8 million annual snakebite victims. Advances in ‘omic’ technologies have empowered toxinologists to comprehensively characterize snake venom compositions, unravel the molecular mechanisms that underpin venom variation, and elucidate the ensuing functional consequences. In this review, we describe how such mechanistic processes have resulted in suites of toxin isoforms that cause diverse pathologies in human snakebite victims and we detail how variation in venom composition can result in treatment failure. Finally, we outline current therapeutic approaches designed to circumvent venom variation and deliver next-generation treatments for the world’s most lethal neglected tropical disease.

Toxinology provides multidirectional and multidimensional opportunities: A personal perspective

In nature, toxins have evolved as weapons to capture and subdue the prey or to counter predators or competitors. When they are inadvertently injected into humans, they cause symptoms ranging from mild discomfort to debilitation and death. Toxinology is the science of studying venoms and toxins that are produced by a wide variety of organisms. In the past, the structure, function and mechanisms of most abundant and/or most toxic components were characterized to understand and to develop strategies to neutralize their toxicity. With recent technical advances, we are able to evaluate and determine the toxin profiles using transcriptomes of venom glands and proteomes of tiny amounts of venom. Enormous amounts of data from these studies have opened tremendous opportunities in many directions of basic and applied research. The lower costs for profiling venoms will further fuel the expansion of toxin database, which in turn will provide greater exciting and bright opportunities in toxin research.

Activity of two key toxin groups in Australian elapid venoms show a strong correlation to phylogeny but not to diet

Background

The relative influence of diet and phylogeny on snake venom activity is a poorly understood aspect of snake venom evolution. We measured the activity of two enzyme toxin groups – phospholipase A₂ (PLA₂), and L-amino acid oxidase (LAAO) – in the venom of 39 species of Australian elapids (40% of terrestrial species diversity) and used linear parsimony and BayesTraits to investigate any correlation between enzyme activity and phylogeny or diet.

Results

PLA₂ activity ranged from 0 to 481 nmol/min/mg of venom, and LAAO activity ranged from 0 to 351 nmol/min/mg. Phylogenetic comparative methods, implemented in BayesTraits showed that enzyme activity was strongly correlated with phylogeny, more so for LAAO activity. For example, LAAO activity was absent in both the Vermicella and Pseudonaja/Oxyuranus clade, supporting previously proposed relationships among these disparate taxa. There was no association between broad dietary categories and either enzyme activity. There was strong evidence for faster initial rates of change over evolutionary time for LAAO (delta parameter mean 0.2), but no such pattern in PLA₂ (delta parameter mean 0.64). There were some exceptions to the phylogenetic patterns of enzyme activity: different PLA₂ activity in the ecologically similar sister-species Denisonia devisi and D. maculata; large inter-specific differences in PLA₂ activity in Hoplocephalus and Austrelaps.

Conclusions

We have shown that phylogeny is a stronger influence on venom enzyme activity than diet for two of the four major enzyme families present in snake venoms. PLA₂ and LAAO activities had contrasting evolutionary dynamics with the higher delta value for PLA₂ Some species/individuals lacked activity in one protein family suggesting that the loss of single protein family may not incur a significant fitness cost.

Beyond the 'big four': Venom profiling of the medically important yet neglected Indian snakes reveals disturbing antivenom deficiencies

Background

Snakebite in India causes the highest annual rates of death (46,000) and disability (140,000) than any other country. Antivenom is the mainstay treatment of snakebite, whose manufacturing protocols, in essence, have remained unchanged for over a century. In India, a polyvalent antivenom is produced for the treatment of envenomations from the so called ‘big four’ snakes: the spectacled cobra (Naja naja), common krait (Bungarus caeruleus), Russell’s viper (Daboia russelii), and saw-scaled viper (Echis carinatus). In addition to the ‘big four’, India is abode to many other species of venomous snakes that have the potential to inflict severe clinical or, even, lethal envenomations in their human bite victims. Unfortunately, specific antivenoms are not produced against these species and, instead, the ‘big four’ antivenom is routinely used for the treatment.

Methods

We characterized the venom compositions, biochemical and pharmacological activities and toxicity profiles (mouse model) of the major neglected yet medically important Indian snakes (E. c. sochureki, B. sindanus, B. fasciatus, and two populations of N. kaouthia) and their closest ‘big four’ congeners. By performing WHO recommended in vitro and in vivo preclinical assays, we evaluated the efficiencies of the commercially marketed Indian antivenoms in recognizing venoms and neutralizing envenomations by these neglected species.

Findings

As a consequence of dissimilar ecologies and diet, the medically important snakes investigated exhibited dramatic inter- and intraspecific differences in their venom profiles. Currently marketed antivenoms were found to exhibit poor dose efficacy and venom recognition potential against the ‘neglected many’. Premium Serums antivenom failed to neutralise bites from many of the neglected species and one of the ‘big four’ snakes (North Indian population of B. caeruleus).

Conclusions

This study unravels disturbing deficiencies in dose efficacy and neutralisation capabilities of the currently marketed Indian antivenoms, and emphasises the pressing need to develop region-specific snakebite therapy for the ‘neglected many’.

Snakebite is a ‘neglected tropical disease’ that majorly affects the rural populations in developing countries. India bears the brunt of snakebites with over 46,000 deaths and 140,000 disabilities, annually. A significant number of these bites are attributed to the widely distributed ‘big four’ snakes, namely spectacled cobra (Naja naja), common krait (Bungarus caeruleus), Russell’s viper (Daboia russelii), and saw-scaled viper (Echis carinatus). The commercial antivenoms marketed in India are only manufactured against these four species, while neglecting many other medically relevant snakes with restricted geographic distribution. Snakebite pathology is dependent on the venom composition of the population/species, which can, in turn, vary intra- and inter-specifically. Though this variation severely limits the cross-population/species antivenom efficacy, envenomations by the neglected snakes in India are treated with the ‘big four’ antivenom. Therefore, to unravel the underlying venom variability, we investigated venom proteomic, biochemical/pharmacological and toxicity profiles of the major neglected Indian snakes and their ‘big four’ relatives. To assess the effectiveness of the ‘big four’ antivenom in treating bites from these neglected snakes, we performed preclinical experiments, which revealed alarming inadequacies of the commercial antivenoms. Our findings accentuate the compelling necessity for the innovation of highly efficacious next-generation snakebite therapy in India.

Omics Technologies for Profiling Toxin Diversity and Evolution in Snake Venom: Impacts on the Discovery of Therapeutic and Diagnostic Agents

Snake venoms are primarily composed of proteins and peptides, and these toxins have developed high selectivity to their biological targets. This makes venoms interesting for exploration into protein evolution and structure-function relationships. A single venom protein superfamily can exhibit a variety of pharmacological effects; these variations in activity originate from differences in functional sites, domains, posttranslational modifications, and the formations of toxin complexes. In this review, we discuss examples of how the major venom protein superfamilies have diversified, as well as how newer technologies in the omics fields, such as genomics, transcriptomics, and proteomics, can be used to characterize both known and unknown toxins.Because toxins are bioactive molecules with a rich diversity of activities, they can be useful as therapeutic and diagnostic agents, and successful examples of toxin applications in these areas are also reviewed. With the current rapid pace of technology, snake venom research and its applications will only continue to expand.

Scorpion Toxins: Positive Selection at a Distal Site Modulates Functional Evolution at a Bioactive Site

The bioactive sites of proteins are those that directly interact with their targets. In many immunity- and predation-related proteins, they frequently experience positive selection for dealing with the changes of their targets from competitors. However, some sites that are far away from the interface between proteins and their targets are also identified to evolve under positive selection. Here, we explore the evolutionary implication of such a site in scorpion α-type toxins affecting sodium (Na⁺) channels (abbreviated as α-ScNaTxs) using a combination of experimental and computational approaches. We found that despite no direct involvement in interaction with Na⁺ channels, mutations at this site by different types of amino acids led to toxicity change on both rats and insects in three α-ScNaTxs, accompanying differential effects on their structures. Molecular dynamics simulations indicated that the mutations changed the conformational dynamics of the positively selected bioactive site-containing functional regions by allosteric communication, suggesting a potential evolutionary correlation between these bioactive sites and the distant nonbioactive site. Our results reveal for the first time the cause of fast evolution at nonbioactive sites of scorpion neurotoxins, which is presumably to adapt to the change of their bioactive sites through coevolution to maintain an active conformation for channel binding. This might aid rational design of scorpion Na⁺ channel toxins with improved phyletic selectivity via modification of a distant nonbioactive site.

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.

Sea Anemone Toxins: A Structural Overview

Sea anemones produce venoms of exceptional molecular diversity, with at least 17 different molecular scaffolds reported to date. These venom components have traditionally been classified according to pharmacological activity and amino acid sequence. However, this classification system suffers from vulnerabilities due to functional convergence and functional promiscuity. Furthermore, for most known sea anemone toxins, the exact receptors they target are either unknown, or at best incomplete. In this review, we first provide an overview of the sea anemone venom system and then focus on the venom components. We have organised the venom components by distinguishing firstly between proteins and non-proteinaceous compounds, secondly between enzymes and other proteins without enzymatic activity, then according to the structural scaffold, and finally according to molecular target.

Crystal structure of the unoccupied murine urokinase-type plasminogen activator receptor (uPAR) reveals a tightly packed DII-DIII unit

The urokinase-type plasminogen activator receptor (uPAR) is a cell surface receptor that is capable of binding to a range of extracellular proteins and triggering a series of proteolytic and signaling events. Previous structural studies of uPAR with its ligands uPA and vitronectin revealed that its three domains (DI, DII, and DIII) form a large hydrophobic cavity to accommodate uPA. In the present study, the structure of unoccupied murine uPAR (muPAR) is determined. The structure of DII and DIII of muPAR is well defined and forms a compact globular unit, while DI could not be traced. Molecular dynamic simulations further confirm the rigid binding interface between DII and DIII. This study shows overall structural flexibility of uPAR in the absence of uPA.

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.

Arthropod venoms: Biochemistry, ecology and evolution

Comprising of over a million described species of highly diverse invertebrates, Arthropoda is amongst the most successful animal lineages to have colonized aerial, terrestrial, and aquatic domains. Venom, one of the many fascinating traits to have evolved in various members of this phylum, has underpinned their adaptation to diverse habitats. Over millions of years of evolution, arthropods have evolved ingenious ways of delivering venom in their targets for self-defence and predation. The morphological diversity of venom delivery apparatus in arthropods is astounding, and includes extensively modified pedipalps, tail (telson), mouth parts (hypostome), fangs, appendages (maxillulae), proboscis, ovipositor (stinger), and hair (urticating bristles). Recent investigations have also unravelled an astonishing venom biocomplexity with molecular scaffolds being recruited from a multitude of protein families. Venoms are a remarkable bioresource for discovering lead compounds in targeted therapeutics. Several components with prospective applications in the development of advanced lifesaving drugs and environment friendly bio-insecticides have been discovered from arthropod venoms. Despite these fascinating features, the composition, bioactivity, and molecular evolution of venom in several arthropod lineages remains largely understudied. This review highlights the prevalence of venom, its mode of toxic action, and the evolutionary dynamics of venom in Arthropoda, the most speciose phylum in the animal kingdom.