Snake Venom Gland Organoids.

Wnt dependency and Lgr5 expression define multiple mammalian epithelial stem cell types. Under defined growth factor conditions, such adult stem cells (ASCs) grow as 3D organoids that recapitulate essential features of the pertinent epithelium. Here, we establish long-term expanding venom gland organoids from several snake species. The newly assembled transcriptome of the Cape coral snake reveals that organoids express high levels of toxin transcripts. Single-cell RNA sequencing of both organoids and primary tissue identifies distinct venom-expressing cell types as well as proliferative cells expressing homologs of known mammalian stem cell markers. A hard-wired regional heterogeneity in the expression of individual venom components is maintained in organoid cultures. Harvested venom peptides reflect crude venom composition and display biological activity. This study extends organoid technology to reptilian tissues and describes an experimentally tractable model system representing the snake venom gland.

Recombinant snake antivenoms get closer to the clinic.

Snakebite envenomings kill ~100 000 victims each year and leave many more with permanent sequelae. Antivenoms have been available for more than 125 years but are in need of innovation. A new study by Khalek et al. highlights broadly neutralizing human monoclonal antibodies (mAbs) that might be used to develop recombinant antivenoms with superior therapeutic benefits.

The genetic regulatory architecture and epigenomic basis for age-related changes in rattlesnake venom.

Developmental phenotypic changes can evolve under selection imposed by age- and size-related ecological differences. Many of these changes occur through programmed alterations to gene expression patterns, but the molecular mechanisms and gene-regulatory networks underlying these adaptive changes remain poorly understood. Many venomous snakes, including the eastern diamondback rattlesnake (Crotalus adamanteus), undergo correlated changes in diet and venom expression as snakes grow larger with age, providing models for identifying mechanisms of timed expression changes that underlie adaptive life history traits. By combining a highly contiguous, chromosome-level genome assembly with measures of expression, chromatin accessibility, and histone modifications, we identified cis-regulatory elements and trans-regulatory factors controlling venom ontogeny in the venom glands of C. adamanteus. Ontogenetic expression changes were significantly correlated with epigenomic changes within genes, immediately adjacent to genes (e.g., promoters), and more distant from genes (e.g., enhancers). We identified 37 candidate transcription factors (TFs), with the vast majority being up-regulated in adults. The ontogenetic change is largely driven by an increase in the expression of TFs associated with growth signaling, transcriptional activation, and circadian rhythm/biological timing systems in adults with corresponding epigenomic changes near the differentially expressed venom genes. However, both expression activation and repression contributed to the composition of both adult and juvenile venoms, demonstrating the complexity and potential evolvability of gene regulation for this trait. Overall, given that age-based trait variation is common across the tree of life, we provide a framework for understanding gene-regulatory-network-driven life-history evolution more broadly.

DeTox: a pipeline for the detection of toxins in venomous organisms.

Venomous organisms have independently evolved the ability to produce toxins 101 times during their evolutionary history, resulting in over 200 000 venomous species. Collectively, these species produce millions of toxins, making them a valuable resource for bioprospecting and understanding the evolutionary mechanisms underlying genetic diversification. RNA-seq is the preferred method for characterizing toxin repertoires, but the analysis of the resulting data remains challenging. While early approaches relied on similarity-based mapping to known toxin databases, recent studies have highlighted the importance of structural features for toxin detection. The few existing pipelines lack an integration between these complementary approaches, and tend to be difficult to run for non-experienced users. To address these issues, we developed DeTox, a comprehensive and user-friendly tool for toxin research. It combines fast execution, parallelization and customization of parameters. DeTox was tested on published transcriptomes from gastropod mollusks, cnidarians and snakes, retrieving most putative toxins from the original articles and identifying additional peptides as potential toxins to be confirmed through manual annotation and eventually proteomic analysis. By integrating a structure-based search with similarity-based approaches, DeTox allows the comprehensive characterization of toxin repertoire in poorly-known taxa. The effect of the taxonomic bias in existing databases is minimized in DeTox, as mirrored in the detection of unique and divergent toxins that would have been overlooked by similarity-based methods. DeTox streamlines toxin annotation, providing a valuable tool for efficient identification of venom components that will enhance venom research in neglected taxa.

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.

Horizontal gene transfer underlies the painful stings of asp caterpillars (Lepidoptera: Megalopygidae).

Larvae of the genus Megalopyge (Lepidoptera: Zygaenoidea: Megalopygidae), known as asp or puss caterpillars, produce defensive venoms that cause severe pain. Here, we present the anatomy, chemistry, and mode of action of the venom systems of caterpillars of two megalopygid species, the Southern flannel moth Megalopyge opercularis and the black-waved flannel moth Megalopyge crispata. We show that megalopygid venom is produced in secretory cells that lie beneath the cuticle and are connected to the venom spines by canals. Megalopygid venoms consist of large aerolysin-like pore-forming toxins, which we have named megalysins, and a small number of peptides. The venom system differs markedly from those of previously studied venomous zygaenoids of the family Limacodidae, suggestive of an independent origin. Megalopygid venom potently activates mammalian sensory neurons via membrane permeabilization and induces sustained spontaneous pain behavior and paw swelling in mice. These bioactivities are ablated by treatment with heat, organic solvents, or proteases, indicating that they are mediated by larger proteins such as the megalysins. We show that the megalysins were recruited as venom toxins in the Megalopygidae following horizontal transfer of genes from bacteria to the ancestors of ditrysian Lepidoptera. Megalopygids have recruited aerolysin-like proteins as venom toxins convergently with centipedes, cnidarians, and fish. This study highlights the role of horizontal gene transfer in venom evolution.

Mast cell-derived prostaglandin D2 limits the subcutaneous absorption of honey bee venom in mice.

Mast cells play pivotal roles in innate host defenses against venom. Activated mast cells release large amounts of prostaglandin D2 (PGD2). However, the role of PGD2 in such host defense remains unclear. We found that c-kit-dependent and c-kit-independent mast cell-specific hematopoietic prostaglandin D synthase (H-pgds) deficiency significantly exacerbated honey bee venom (BV)-induced hypothermia and increased mortality rates in mice. BV absorption via postcapillary venules in the skin was accelerated upon endothelial barrier disruption resulting in increased plasma venom concentrations. These results suggest that mast cell-derived PGD2 may enhance host defense against BV and save lives by inhibiting BV absorption into circulation.

Peptide toxins that target vertebrate voltage-gated sodium channels underly the painful stings of harvester ants.

Harvester ants (genus Pogonomyrmex) are renowned for their stings which cause intense, long-lasting pain, and other neurotoxic symptoms in vertebrates. Here, we show that harvester ant venoms are relatively simple and composed largely of peptide toxins. One class of peptides is primarily responsible for the long-lasting local pain of envenomation via activation of peripheral sensory neurons. These hydrophobic, cysteine-free peptides potently modulate mammalian voltage-gated sodium (NaV) channels, reducing the voltage threshold for activation and inhibiting channel inactivation. These toxins appear to have evolved specifically to deter vertebrates.

Alpha-1-B glycoprotein (A1BG) inhibits sterol-binding and export by CRISP2.

Proteins belonging to the CAP superfamily are present in all kingdoms of life and have been implicated in various processes, including sperm maturation and cancer progression. They are mostly secreted glycoproteins and share a unique conserved CAP domain. The precise mode of action of these proteins, however, has remained elusive. Saccharomyces cerevisiae expresses three members of this protein family, which bind sterols in vitro and promote sterol secretion from cells. This sterol-binding and export function of yeast Pry proteins is conserved in the mammalian CRISP proteins and other CAP superfamily members. CRISP3 is an abundant protein of the human seminal plasma and interacts with alpha-1-B glycoprotein (A1BG), a human plasma glycoprotein that is upregulated in different types of cancers. Here we examined whether the interaction between CRISP proteins and A1BG affects the sterol-binding function of CAP family members. Co-expression of A1BG with CAP proteins abolished their sterol export function in yeast and their interaction inhibits sterol-binding in vitro. We map the interaction between A1BG and CRISP2 to the third of five repeated immunoglobulin-like (Ig) domains within A1BG. Interestingly, the interaction between A1BG and CRISP2 requires magnesium, suggesting that coordination of Mg2+ by the highly conserved tetrad residues within the CAP domain is essential for a stable interaction between the two proteins. The observation that A1BG modulates the sterol-binding function of CRISP2, has potential implications for the role of A1BG and related Ig domain containing proteins in cancer progression and the toxicity of reptile venoms containing CRISP proteins.

Prey Shifts Drive Venom Evolution in Cone Snails.

Venom systems are complex traits that have independently emerged multiple times in diverse plant and animal phyla. Within each venomous lineage there typically exists interspecific variation in venom composition where several factors have been proposed as drivers of variation, including phylogeny and diet. Understanding these factors is of broad biological interest and has implications for the development of antivenom therapies and venom-based drug discovery. Because of their high species richness and the presence of several major evolutionary prey shifts, venomous marine cone snails (genus Conus) provide an ideal system to investigate drivers of interspecific venom variation. Here, by analyzing the venom gland expression profiles of ∼3,000 toxin genes from 42 species of cone snail, we elucidate the role of prey-specific selection pressures in shaping venom variation. By analyzing overall venom composition and individual toxin structures, we demonstrate that the shifts from vermivory to piscivory in Conus are complemented by distinct changes in venom composition independent of phylogeny. In vivo injections of venom from piscivorous cone snails in fish further showed a higher potency compared with venom of nonpiscivores demonstrating a selective advantage. Together, our findings provide compelling evidence for the role of prey shifts in directing the venom composition of cone snails and expand our understanding of the mechanisms of venom variation and diversification.

Sea Anemone Membrane Attack Complex/Perforin Superfamily Demonstrates an Evolutionary Transitional State between Venomous and Developmental Functions.

Gene duplication is a major force driving evolutionary innovation. A classic example is generating new animal toxins via duplication of physiological protein-encoding genes and recruitment into venom. While this process drives the innovation of many animal venoms, reverse recruitment of toxins into nonvenomous cells remains unresolved. Using comparative genomics, we find members of the Membrane Attack Complex and Perforin Family (MAC) have been recruited into venom-injecting cells (cnidocytes), in soft and stony corals and sea anemones, suggesting that the ancestral MAC was a cnidocyte expressed toxin. Further investigation into the model sea anemone Nematostella vectensis reveals that three members have undergone Nematostella-specific duplications leading to their reverse recruitment into endomesodermal cells. Furthermore, simultaneous knockdown of all three endomesodermally expressed MACs leads to mis-development, supporting that these paralogs have nonvenomous function. By resolving the evolutionary history and function of MACs in Nematostella, we provide the first proof for reverse recruitment from venom to organismal development.

Whole transcriptome sequencing reveals the activity of the PLA2 family members in Androctonus crassicauda (Scorpionida: Buthidae) venom gland.

Phospholipase A2 is the most abundant venom gland enzyme, whose activity leads to the activation of the inflammatory response by accumulating lipid mediators. This study aimed to identify, classify, and investigate the properties of venom PLA2 isoforms. Then, the present findings were confirmed by chemically measuring the activity of PLA2. The sequences representing PLA2 annotation were extracted from the Androctonus crassicauda transcriptome dataset using BLAS searches against the local PLA2 database. We found several cDNA sequences of PLA2 classified and named by conducting multiple searches as platelet-activating factor acetylhydrolases, calcium-dependent PLA2s, calcium-independent PLA2s, and secreted PLA2s. The largest and smallest isoforms of these proteins range between approximately 70.34 kDa (iPLA2) and 17.75 kDa (cPLA2). Among sPLA2 isoforms, sPLA2GXIIA and sPLA2G3 with ORF encoding 169 and 299 amino acids are the smallest and largest secreted PLA2, respectively. These results collectively suggested that A. crassicauda venom has PLA2 activity, and the members of this protein family may have important biological roles in lipid metabolism. This study also revealed the interaction between members of PLA2s in the PPI network. The results of this study would greatly help with the classification, evolutionary relationships, and interactions between PLA2 family proteins in the gene network.

Selective block of human Kv1.1 channels and an epilepsy-associated gain-of-function mutation by AETX-K peptide.

Dysfunction of the human voltage-gated K+ channel Kv1.1 has been associated with epilepsy, multiple sclerosis, episodic ataxia, myokymia, and cardiorespiratory dysregulation. We report here that AETX-K, a sea anemone type I (SAK1) peptide toxin we isolated from a phage display library, blocks Kv1.1 with high affinity (Ki ~ 1.6 pM) and notable specificity, inhibiting other Kv channels we tested a million-fold less well. Nuclear magnetic resonance (NMR) was employed both to determine the three-dimensional structure of AETX-K, showing it to employ a classic SAK1 scaffold while exhibiting a unique electrostatic potential surface, and to visualize AETX-K bound to the Kv1.1 pore domain embedded in lipoprotein nanodiscs. Study of Kv1.1 in Xenopus oocytes with AETX-K and point variants using electrophysiology demonstrated the blocking mechanism to employ a toxin-channel configuration we have described before whereby AETX-K Lys23 , two positions away on the toxin interaction surface from the classical blocking residue, enters the pore deeply enough to interact with K+ ions traversing the pathway from the opposite side of the membrane. The mutant channel Kv1.1-L296 F is associated with pharmaco-resistant multifocal epilepsy in infants because it significantly increases K+ currents by facilitating opening and slowing closure of the channels. Consistent with the therapeutic potential of AETX-K for Kv1.1 gain-of-function-associated diseases, AETX-K at 4 pM decreased Kv1.1-L296 F currents to wild-type levels; further, populations of heteromeric channels formed by co-expression Kv1.1 and Kv1.2, as found in many neurons, showed a Ki of ~10 nM even though homomeric Kv1.2 channels were insensitive to the toxin (Ki > 2000 nM).

Revealing molecular determinants governing mambalgin-3 pharmacology at acid-sensing ion channel 1 variants.

Acid-sensing ion channels (ASICs) are trimeric proton-gated cation channels that play a role in neurotransmission and pain sensation. The snake venom-derived peptides, mambalgins, exhibit potent analgesic effects in rodents by inhibiting central ASIC1a and peripheral ASIC1b. Despite their distinct species- and subtype-dependent pharmacology, previous structure-function studies have focussed on the mambalgin interaction with ASIC1a. Currently, the specific channel residues responsible for this pharmacological profile, and the mambalgin pharmacophore at ASIC1b remain unknown. Here we identify non-conserved residues at the ASIC1 subunit interface that drive differences in the mambalgin pharmacology from rat ASIC1a to ASIC1b, some of which likely do not make peptide binding interactions. Additionally, an amino acid variation below the core binding site explains potency differences between rat and human ASIC1. Two regions within the palm domain, which contribute to subtype-dependent effects for mambalgins, play key roles in ASIC gating, consistent with subtype-specific differences in the peptides mechanism. Lastly, there is a shared primary mambalgin pharmacophore for ASIC1a and ASIC1b activity, with certain peripheral peptide residues showing variant-specific significance for potency. Through our broad mutagenesis studies across various species and subtype variants, we gain a more comprehensive understanding of the pharmacophore and the intricate molecular interactions that underlie ligand specificity. These insights pave the way for the development of more potent and targeted peptide analogues required to advance our understating of human ASIC1 function and its role in disease.

Exploring oak processionary caterpillar induced lepidopterism (Part 1): unveiling molecular insights through transcriptomics and proteomics.

Lepidopterism, a skin inflammation condition caused by direct or airborne exposure to irritating hairs (setae) from processionary caterpillars, is becoming a significant public health concern. Recent outbreaks of the oak processionary caterpillar (Thaumetopoea processionea) have caused noteworthy health and economic consequences, with a rising frequency expected in the future, exacerbated by global warming promoting the survival of the caterpillar. Current medical treatments focus on symptom relief due to the lack of an effective therapy. While the source is known, understanding the precise causes of symptoms remain incomplete understood. In this study, we employed an advanced method to extract venom from the setae and identify the venom components through high-quality de novo transcriptomics, venom proteomics, and bioinformatic analysis. A total of 171 venom components were identified, including allergens, odorant binding proteins, small peptides, enzymes, enzyme inhibitors, and chitin biosynthesis products, potentially responsible for inflammatory and allergic reactions. This work presents the first comprehensive proteotranscriptomic database of T. processionea, contributing to understanding the complexity of lepidopterism. Furthermore, these findings hold promise for advancing therapeutic approaches to mitigate the global health impact of T. processionea and related caterpillars.

Convergent evolution of venom gland transcriptomes across Metazoa.

Animals have repeatedly evolved specialized organs and anatomical structures to produce and deliver a mixture of potent bioactive molecules to subdue prey or predators-venom. This makes it one of the most widespread, convergent functions in the animal kingdom. Whether animals have adopted the same genetic toolkit to evolved venom systems is a fascinating question that still eludes us. Here, we performed a comparative analysis of venom gland transcriptomes from 20 venomous species spanning the main Metazoan lineages to test whether different animals have independently adopted similar molecular mechanisms to perform the same function. We found a strong convergence in gene expression profiles, with venom glands being more similar to each other than to any other tissue from the same species, and their differences closely mirroring the species phylogeny. Although venom glands secrete some of the fastest evolving molecules (toxins), their gene expression does not evolve faster than evolutionarily older tissues. We found 15 venom gland-specific gene modules enriched in endoplasmic reticulum stress and unfolded protein response pathways, indicating that animals have independently adopted stress response mechanisms to cope with mass production of toxins. This, in turn, activates regulatory networks for epithelial development, cell turnover, and maintenance, which seem composed of both convergent and lineage-specific factors, possibly reflecting the different developmental origins of venom glands. This study represents a first step toward an understanding of the molecular mechanisms underlying the repeated evolution of one of the most successful adaptive traits in the animal kingdom.

A novel broad spectrum venom metalloproteinase autoinhibitor in the rattlesnake Crotalus atrox evolved via a shift in paralog function.

The complexity of snake venom composition reflects adaptation to the diversity of prey and may be driven at times by a coevolutionary arms race between snakes and venom-resistant prey. However, many snakes are also resistant to their own venom due to serum-borne inhibitors of venom toxins, which raises the question of how snake autoinhibitors maintain their efficacy as venom proteins evolve. To investigate this potential three-way arms race among venom, prey, and autoinhibitors, we have identified and traced the evolutionary origin of serum inhibitors of snake venom metalloproteinases (SVMPs) in the Western Diamondback rattlesnake Crotalus atrox which possesses the largest known battery of SVMP genes among crotalids examined. We found that C. atrox expresses five members of a Fetuin A-related metalloproteinase inhibitor family but that one family member, FETUA-3, is the major SVMP inhibitor that binds to approximately 20 different C. atrox SVMPs and inhibits activities of all three SVMP classes. We show that the fetua-3 gene arose deep within crotalid evolution before the origin of New World species but, surprisingly, fetua-3 belongs to a different paralog group than previously identified SVMP inhibitors in Asian and South American crotalids. Conversely, the C. atrox FETUA-2 ortholog of previously characterized crotalid SVMP inhibitors shows limited activity against C. atrox SVMPs. These results reveal that there has been a functional evolutionary shift in the major SVMP inhibitor in the C. atrox lineage as the SVMP family expanded and diversified in the Crotalus lineage. This broad-spectrum inhibitor may be of potential therapeutic interest.

Breaking muscle: neurotoxic and myotoxic effects of Central American snake venoms and the relative efficacies of antivenom and varespladib.

BACKGROUND: The snake genera Atropoides, Cerrophidion, and Metlapilcoatlus form a clade of neotropical pit vipers distributed across Mexico and Central America. This study evaluated the myotoxic and neurotoxic effects of nine species of Atropoides, Cerrophidion, and Metlapilcoatlus, and the neutralising efficacy of the ICP antivenom from Costa Rica against these effects, in the chick biventer cervicis nerve-muscle preparation. Given the prominence of PLA2s within the venom proteomes of these species, we also aimed to determine the neutralising potency of the PLA2 inhibitor, varespladib. RESULTS: All venoms showed myotoxic and potential neurotoxic effects, with differential intra-genera and inter-genera potency. This variation was also seen in the antivenom ability to neutralise the muscle damaging pathophysiological effects observed. Variation was also seen in the relative response to the PLA2 inhibitor varespladib. While the myotoxic effects of M. mexicanus and M. nummifer venoms were effectively neutralised by varespladib, indicating myotoxicity is PLA2 mediated, those of C. godmani and M. olmec venoms were not, revealing that the myotoxicity is driven by non-PLA2 toxin types. CONCLUSIONS: This study characterises the myotoxic and neurotoxic venom activity, as well as neutralisation of venom effects from the Atropoides, Cerrophidion, and Metlapilcoatlus clade of American crotalids. Our findings contribute significant clinical and evolutionary knowledge to a clade of poorly researched snakes. In addition, these results provide a platform for future research into the reciprocal interaction between ecological niche specialisation and venom evolution, as well as highlighting the need to test purified toxins to accurately evaluate the potential effects observed in these venoms.

Host and venom evolution in parasitoid wasps: does independently adapting to the same host shape the evolution of the venom gland transcriptome?

BACKGROUND: Venoms have repeatedly evolved over 100 occasions throughout the animal tree of life, making them excellent systems for exploring convergent evolutionary novelty. Growing evidence supports that venom evolution is predominantly driven by prey or host-related selection pressures, and the expression patterns of venom glands reflect adaptive evolution. However, it remains elusive whether the evolution of expression patterns in venom glands is likewise a convergent evolution driven by their prey/host species. RESULTS: We utilized parasitoid wasps that had independently adapted to Drosophila hosts as models to investigate the convergent evolution of venom gland transcriptomes in 19 hymenopteran species spanning ~ 200 million years of evolution. Comparative transcriptome analysis reveals that the global expression patterns among the venom glands of Drosophila parasitoid wasps do not achieve higher similarity compared to non-Drosophila parasitoid wasps. Further evolutionary analyses of expression patterns at the single gene, orthogroup, and Gene Ontology (GO) term levels indicate that some orthogroups/GO terms show correlation with the Drosophila parasitoid wasps. However, these groups rarely include genes highly expressed in venom glands or putative venom genes in the Drosophila parasitoid wasps. CONCLUSIONS: Our study suggests that convergent evolution may not play a predominant force shaping gene expression levels in the venom gland of the Drosophila parasitoid wasps, offering novel insights into the co-evolution between venom and prey/host.

From birth to bite: the evolutionary ecology of India's medically most important snake venoms.

BACKGROUND: Snake venoms can exhibit remarkable inter- and intraspecific variation. While diverse ecological and environmental factors are theorised to explain this variation, only a handful of studies have attempted to unravel their precise roles. This knowledge gap not only impedes our understanding of venom evolution but may also have dire consequences on snakebite treatment. To address this shortcoming, we investigated the evolutionary ecology of venoms of Russell's viper (Daboia russelii) and spectacled cobra (Naja naja), India's two clinically most important snakes responsible for an alarming number of human deaths and disabilities. METHODOLOGY: Several individuals (n = 226) of D. russelii and N. naja belonging to multiple clutches (n = 9) and their mothers were maintained in captivity to source ontogenetic stage-specific venoms. Using various in vitro and in vivo assays, we assessed the significance of prey, ontogeny and sex in driving venom composition, function, and potency. RESULTS: Considerable ontogenetic shifts in venom profiles were observed in D. russelii, with the venoms of newborns being many times as potent as juveniles and adults against mammalian (2.3-2.5 ×) and reptilian (2-10 ×) prey. This is the first documentation of the ontogenetic shift in viperine snakes. In stark contrast, N. naja, which shares a biogeographic distribution similar to D. russelii, deployed identical biochemical cocktails across development. Furthermore, the binding kinetics of cobra venom toxins against synthetic target receptors from various prey and predators shed light on the evolutionary arms race. CONCLUSIONS: Our findings, therefore, provide fascinating insights into the roles of ecology and life history traits in shaping snake venoms.