Snake venoms are complex mixtures of toxins that differ on interspecific (between species) and intraspecific (within species) levels. Whether venom variation within a group of closely related species is explained by the presence, absence and/or relative abundances of venom toxins remains largely unknown. Taipans (Oxyuranus spp.) and brown snakes (Pseudonaja spp.) represent medically relevant species of snakes across the Australasian region and provide an excellent model clade for studying interspecific and intraspecific venom variation. Using liquid chromatography with ultraviolet and mass spectrometry detection, we analyzed a total of 31 venoms covering all species of this monophyletic clade, including widespread localities. Our results reveal major interspecific and intraspecific venom variation in Oxyuranus and Pseudonaja species, partially corresponding with their geographical regions and phylogenetic relationships. This extensive venom variability is generated by a combination of the absence/presence and differential abundance of venom toxins. Our study highlights that venom systems can be highly dynamical on the interspecific and intraspecific levels and underscores that the rapid toxin evolvability potentially causes major impacts on neglected tropical snakebites.

Snakebite is considered a neglected tropical disease, and it is one of the most intricate ones. The variability found in snake venom is what makes it immensely complex to study. These variations are present both in the big and the small molecules found in snake venom. This study focused on examining the variability found in the venom’s small molecules (i.e., mass range of 100–1000 Da) between two main families of venomous snakes—Elapidae and Viperidae—managing to create a model able to classify unknown samples by means of specific features, which can be extracted from their LC–MS data and output in a comprehensive list. The developed model also allowed further insight into the composition of snake venom by highlighting the most relevant metabolites of each group by clustering similarly composed venoms. The model was created by means of support vector machines and used 20 features, which were merged into 10 principal components. All samples from the first and second validation data subsets were correctly classified. Biological hypotheses relevant to the variation regarding the metabolites that were identified are also given.

Envenoming resulting from snakebites is recognized as a priority neglected tropical disease by The World Health Organization. The Bothrops genus, consisting of different pitviper species, is considered the most medically significant taxa in Central and South America. Further research into Bothrops venom composition is important to aid in the development of safer and more effective snakebite treatments. In addition, the discovery of Bothrops toxins that could potentially be used for medical or diagnostic purposes is of interest to the pharmaceutical industry. This study aimed to employ high-throughput (HT) venomics to qualitatively analyze venom composition while utilizing coagulation bioassays for identifying coagulopathic toxins and characterizing coagulopathic activity in various Bothrops venoms. Using the recently demonstrated HT venomics workflow in combination with post-column coagulopathic bioassaying, focus was placed at anticoagulant toxins. Well-known procoagulant toxins were also investigated, taking into account that using the HT venomics workflow, procoagulant toxins are especially prone to denaturation during the reversed-phase chromatographic separations performed in the workflow. The findings revealed that the venoms of B. atrox and B. jararaca harbored procoagulant toxins, whereas those of B. alternatus and B. neuwiedi contained both procoagulant and anticoagulant toxins. In general, anticoagulation was associated with phospholipases A₂s, while procoagulation was associated with snake venom metalloproteinases and snake venom serine proteases. These results showed the identification of coagulopathic venom toxins in the Bothrops venoms analyzed using multiple analytical methods that complement each other. Additionally, each venom underwent qualitative characterization of its composition.

The variation in venom between and within snake species has significant implications for snakebite treatment. This highlights the critical importance of studying venom composition and its variations, not only for medical purposes but also from an evolutionary perspective. This study explores analytics for characterizing venom variability, focusing on venom toxin accurate masses, and emphasizes how the complexity of studying snake venom variability can be addressed by using liquid chromatography mass spectrometry (LC-MS) analysis with bioinformatics tools. This was demonstrated by investigating LC-MS data obtained from the venoms of 15 true cobras (Naja spp.), 5 mambas (Dendroaspis spp.) and 28 vipers (Crotalus and Bothrops spp.; total of 20 Elapidae and 28 Viperidae venoms), with newly developed bioinformatics tools. The measured LC-MS data was processed in an automated fashion and sorted based on the monoisotopic accurate masses of all toxins found, their peak intensities, and their retention times in LC. The data was then investigated using bioinformatic tools, before the toxin data available in open-source databases was used to predict the class of a toxin by means of its mass. This study highlights the importance of studying venom variability, which is performed by our combinatorial approach of intact-toxin analysis and toxin grouping by accurate mass.

The cytotoxicity caused by snake venoms is a serious medical problem that greatly contributes to the morbidity observed in snakebite patients. The cytotoxic components found in snake venoms belong to a variety of toxin classes and may cause cytotoxic effects by targeting a range of molecular structures, including cellular membranes, the extracellular matrix (ECM) and the cytoskeleton. Here, we present a high-throughput assay (384-well plate) that monitors ECM degradation by snake venom toxins via the application of fluorescent versions of model ECM substrates, specifically gelatin and collagen type I. Both crude venoms and fractionated toxins of a selection of medically relevant viperid and elapid species, separated via size-exclusion chromatography, were studied using the self-quenching, fluorescently labelled ECM–polymer substrates. The viperid venoms showed significantly higher proteolytic degradation when compared to elapid venoms, although the venoms with higher snake venom metalloproteinase content did not necessarily exhibit stronger substrate degradation than those with a lower one. Gelatin was generally more readily cleaved than collagen type I. In the viperid venoms, which were subjected to fractionation by SEC, two (B. jararaca and C. rhodostoma, respectively) or three (E. ocellatus) active proteases were identified. Therefore, the assay allows the study of proteolytic activity towards the ECM in vitro for crude and fractionated venoms.

Snakebite is considered a concerning issue and a neglected tropical disease. Three-finger toxins (3FTxs) in snake venoms primarily cause neurotoxic effects since they have high affinity for nicotinic acetylcholine receptors (nAChRs). Their small molecular size makes 3FTxs weakly immunogenic and therefore not appropriately targeted by current antivenoms. This study aims at presenting and applying an analytical method for investigating the therapeutic potential of the acetylcholine-binding protein (AChBP), an efficient nAChR mimic that can capture 3FTxs, for alternative treatment of elapid snakebites. In this analytical methodology, snake venom toxins were separated and characterised using high-performance liquid chromatography coupled with mass spectrometry (HPLC-MS) and high-throughput venomics. By subsequent nanofractionation analytics, binding profiling of toxins to the AChBP was achieved with a post-column plate reader-based fluorescence-enhancement ligand displacement bioassay. The integrated method was established and applied to profiling venoms of six elapid snakes (Naja mossambica, Ophiophagus hannah, Dendroaspis polylepis, Naja kaouthia, Naja haje and Bungarus multicinctus). The methodology demonstrated that the AChBP is able to effectively bind long-chain 3FTxs with relatively high affinity, but has low or no binding affinity towards short-chain 3FTxs, and as such provides an efficient analytical platform to investigate binding affinity of 3FTxs to the AChBP and mutants thereof and to rapidly identify bound toxins.