Metabolome-Based Classification of Snake Venoms by Bioinformatic Tools

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

High-Throughput Venomics

In this study, we present high-throughput (HT) venomics, a novel analytical strategy capable of performing a full proteomic analysis of a snake venom within 3 days. This methodology comprises a combination of RP-HPLC-nanofractionation analytics, mass spectrometry analysis, automated in-solution tryptic digestion, and high-throughput proteomics. In-house written scripts were developed to process all the obtained proteomics data by first compiling all Mascot search results for a single venom into a single Excel sheet. Then, a second script plots each of the identified toxins in so-called Protein Score Chromatograms (PSCs). For this, for each toxin, identified protein scores are plotted on the y-axis versus retention times of adjacent series of wells in which a toxin was fractionated on the x-axis. These PSCs allow correlation with parallel acquired intact toxin MS data. This same script integrates the PSC peaks from these chromatograms for semiquantitation purposes. This new HT venomics strategy was performed on venoms from diverse medically important biting species; Calloselasma rhodostoma, Echis ocellatus, Naja pallida, Bothrops asper, Bungarus multicinctus, Crotalus atrox, Daboia russelii, Naja naja, Naja nigricollis, Naja mossambica, and Ophiophagus hannah. Our data suggest that high-throughput venomics represents a valuable new analytical tool for increasing the throughput by which we can define venom variation and should greatly aid in the future development of new snakebite treatments by defining toxin composition

Varespladib Inhibits the Phospholipase A 2 and Coagulopathic Activities of Venom Components from Hemotoxic Snakes

Phospholipase A2 (PLA2) enzymes are important toxins found in many snake venoms, and they can exhibit a variety of toxic activities including causing hemolysis and/or anticoagulation. In this study, the inhibiting effects of the small molecule PLA2 inhibitor varespladib on snake venom PLA2s was investigated by nanofractionation analytics, which combined chromatography, mass spectrometry (MS), and bioassays. The venoms of the medically important snake species Bothrops asper, Calloselasma rhodostoma, Deinagkistrodon acutus, Daboia russelii, Echis carinatus, Echis ocellatus, and Oxyuranus scutellatus were separated by liquid chromatography (LC) followed by nanofractionation and interrogation of the fractions by a coagulation assay and a PLA2 assay. Next, we assessed the ability of varespladib to inhibit the activity of enzymatic PLA2s and the coagulopathic toxicities induced by fractionated snake venom toxins, and identified these bioactive venom toxins and those inhibited by varespladib by using parallel recorded LC-MS data and proteomics analysis. We demonstrated here that varespladib was not only capable of inhibiting the PLA2 activities of hemotoxic snake venoms, but can also effectively neutralize the coagulopathic toxicities (most profoundly anticoagulation) induced by venom toxins. While varespladib effectively inhibited PLA2 toxins responsible for anticoagulant effects, we also found some evidence that this inhibitory molecule can partially abrogate procoagulant venom effects caused by different toxin families. These findings further emphasize the potential clinical utility of varespladib in mitigating the toxic effects of certain snakebites

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

Studying Venom Toxin Variation Using Accurate Masses from Liquid Chromatography–Mass Spectrometry Coupled with Bioinformatic Tools

This study provides a new methodology for the rapid analysis of numerous venom samples in an automated fashion. Here, we use LC-MS (Liquid Chromatography–Mass Spectrometry) for venom separation and toxin analysis at the accurate mass level combined with new in-house written bioinformatic scripts to obtain high-throughput results. This analytical methodology was validated using 31 venoms from all members of a monophyletic clade of Australian elapids: brown snakes (Pseudonaja spp.) and taipans (Oxyuranus spp.). In a previous study, we revealed extensive venom variation within this clade, but the data was manually processed and MS peaks were integrated into a time-consuming and labour-intensive approach. By comparing the manual approach to our new automated approach, we now present a faster and more efficient pipeline for analysing venom variation. Pooled venom separations with post-column toxin fractionations were performed for subsequent high-throughput venomics to obtain toxin IDs correlating to accurate masses for all fractionated toxins. This workflow adds another dimension to the field of venom analysis by providing opportunities to rapidly perform in-depth studies on venom variation. Our pipeline opens new possibilities for studying animal venoms as evolutionary model systems and investigating venom variation to aid in the development of better antivenoms

A Combined Bioassay and Nanofractionation Approach to Investigate the Anticoagulant Toxins of Mamba and Cobra Venoms and Their Inhibition by Varespladib.

Envenomation by elapid snakes primarily results in neurotoxic symptoms and, consequently, are the primary focus of therapeutic research concerning such venoms. However, mounting evidence suggests these venoms can additionally cause coagulopathic symptoms, as demonstrated by some Asian elapids and African spitting cobras. This study sought to investigate the coagulopathic potential of venoms from medically important elapids of the genera Naja (true cobras), Hemachatus (rinkhals), and Dendroaspis (mambas). Crude venoms were bioassayed for coagulant effects using a plasma coagulation assay before RPLC/MS was used to separate and identify venom toxins in parallel with a nanofractionation module. Subsequently, coagulation bioassays were performed on the nanofractionated toxins, along with in-solution tryptic digestion and proteomics analysis. These experiments were then repeated on both crude venoms and on the nanofractionated venom toxins with the addition of either the phospholipase A2 (PLA2) inhibitor varespladib or the snake venom metalloproteinase (SVMP) inhibitor marimastat. Our results demonstrate that various African elapid venoms have an anticoagulant effect, and that this activity is significantly reduced for cobra venoms by the addition of varespladib, though this inhibitor had no effect against anticoagulation caused by mamba venoms. Marimastat showed limited capacity to reduce anticoagulation in elapids, affecting only N. haje and H. haemachatus venom at higher doses. Proteomic analysis of nanofractionated toxins revealed that the anticoagulant toxins in cobra venoms were both acidic and basic PLA2s, while the causative toxins in mamba venoms remain uncertain. This implies that while PLA2 inhibitors such as varespladib and metalloproteinase inhibitors such as marimastat are viable candidates for novel snakebite treatments, they are not likely to be effective against mamba envenomings

Solenodon genome reveals convergent evolution of venom in eulipotyphlan mammals

Venom systems are key adaptations that have evolved throughout the tree of life and typically facilitate predation or defense. Despite venoms being model systems for studying a variety of evolutionary and physiological processes, many taxonomic groups remain understudied, including venomous mammals. Within the order Eulipotyphla, multiple shrew species and solenodons have oral venom systems. Despite morphological variation of their delivery systems, it remains unclear whether venom represents the ancestral state in this group or is the result of multiple independent origins. We investigated the origin and evolution of venom in eulipotyphlans by characterizing the venom system of the endangered Hispaniolan solenodon (Solenodon paradoxus). We constructed a genome to underpin proteomic identifications of solenodon venom toxins, before undertaking evolutionary analyses of those constituents, and functional assessments of the secreted venom. Our findings show that solenodon venom consists of multiple paralogous kallikrein 1 (KLK1) serine proteases, which cause hypotensive effects in vivo, and seem likely to have evolved to facilitate vertebrate prey capture. Comparative analyses provide convincing evidence that the oral venom systems of solenodons and shrews have evolved convergently, with the 4 independent origins of venom in eulipotyphlans outnumbering all other venom origins in mammals. We find that KLK1s have been independently coopted into the venom of shrews and solenodons following their divergence during the late Cretaceous, suggesting that evolutionary constraints may be acting on these genes. Consequently, our findings represent a striking example of convergent molecular evolution and demonstrate that distinct structural backgrounds can yield equivalent functions

Highly Evolvable: Investigating Interspecific and Intraspecific Venom Variation in Taipans (Oxyuranus spp.) and Brown Snakes (Pseudonaja spp.)

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

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

Monitoring Snake Venom-Induced Extracellular Matrix Degradation and Identifying Proteolytically Active Venom Toxins Using Fluorescently Labeled Substrates

Snakebite envenoming is an important public health issue with devastating consequences and annual mortality rates that range between 81,000 and 138,000. Snake venoms may cause a range of pathophysiological effects affecting the nervous system and the cardiovascular system. Moreover, snake venom may have tissue-damaging activities that result in lifelong morbidities such as amputations, muscle degeneration, and organ malfunctioning. The tissue-damaging components in snake venoms comprise multiple toxin classes with various molecular targets including cellular membranes and the extracellular matrix (ECM). In this study, we present multiple assay formats that enable investigation of snake venom-induced ECM degradation using a variety of (dye-quenched) fluorescently labeled ECM components. Using a combinatorial approach, we were able to characterise different proteolytic profiles for different medically relevant snake venoms, followed by identification of the responsible components within the snake venoms. This workflow could provide valuable insights into the key mechanisms by which proteolytic venom components exert their effects and could therefore prove useful for the development of effective snakebite treatments against this severe pathology