Snake Venomics and Antivenomics of Bothrops diporus, a Medically Important Pitviper in Northeastern Argentina

Comparative in vitro and in silico analysis of the ability of basic Asp49 phospholipase A2 and Lys49-phospholipase A2-like myotoxins from Bothrops diporus venom to inhibit the metastatic potential of murine mammary tumor cells and endothelial cell tubulogenesis: Asp49 vs Lys49 phospholipases A2: Inhibition of metastasis and angiogenesis

Snake venoms are a complex mixture of proteins and polypeptides that represent a valuable source of potential molecular tools for understanding physiological processes for the development of new drugs. In this study two major PLA2s, named PLA2-I (Asp49) and PLA2-II (Lys49), isolated from the venom of Bothrops diporus from Northeastern Argentina, have shown cytotoxic effects on LM3 murine mammary tumor cells, with PLA2-II-like exhibiting a stronger effect compared to PLA2-I. At sub-cytotoxic levels, both PLA2s inhibited adhesion, migration, and invasion of these adenocarcinoma cells. Moreover, these toxins hindered tubulogenesis in endothelial cells, implicating a potential role in inhibiting tumor angiogenesis. All these inhibitory effects were more pronounced for the catalytically-inactive toxin. Additionally, in silico studies strongly suggest that this PLA2-II-like myotoxin could effectively block fibronectin binding to the integrin receptor, offering a dual advantage over PLA2-I in interacting with the αVβ3 integrin. In conclusion, this study reports for the first time, integrating both in vitro and in silico approaches, a comparative analysis of the antimetastatic and antiangiogenic potential effects of two isoforms, an Asp49 PLA2-I and a Lys49 PLA2-II-like, both isolated from Bothrops diporus venom.

Analysis of the genus B othrops snake venom: An inter and intraspecific comparative study

The genus Bothrops are considered Category 1 of medical importance by the World Health Organization, responsible for approximately 85 % of snakebites occurring throughout Brazil. Main factors determining snake venom variations can be genetics, diet, gender, geographic distribution, age, or even seasonality. In this study, we compared the composition of protein profile, biochemical activities, and immunorecognition of toxins present in the venom of eight adults of Bothrops species (B. alternatus, B. atrox, B. jararaca, B. jararacussu, B. leucurus, B. moojeni, B. neuwiedi and B. pauloensis). The following methods were used to analyze the venoms: protein dosage; electrophoresis in polyacrylamide gel containing SDS; High Performance Liquid Chromatography - Reverse Phase; enzymatic activities, western blotting and Enzyme Linked Immuno Sorbent Assay. The results show inter and intraspecific differences in the electrophoretic profile. LAAO and PLA2 activities, in general, were higher in males than females and proteolytic activity was higher in females than males. The bothropic antivenom produced by Instituto Butantan recognized most of the protein bands in all Bothrops species analyzed, with only the regions between 37 and 25 kDa presenting lower intensity. A notable variability in the chromatograms was observed. Bothrops venom demonstrated inter-intraspecific disparities in protein composition and biochemical activity.

Exploring the antibacterial potential of venoms from Argentinian animals

The resistance to antimicrobials developed by several bacterial species has become one of the main health problems in recent decades. It has been widely reported that natural products are important sources of antimicrobial compounds. Considering that animal venoms are under-explored in this line of research, in this study, we screened the antibacterial activity of venoms of eight snake and five lepidopteran species from northeastern Argentina. Twofold serial dilutions of venoms were tested by the agar well-diffusion method and the minimum inhibitory concentration (MIC) determination against seven bacterial strains. We studied the comparative protein profile of the venoms showing antibacterial activity. Only the viperid and elapid venoms showed remarkable dose-dependent antibacterial activity towards most of the strains tested. Bothrops diporus venom showed the lowest MIC values against all the strains, and S. aureus ATCC 25923 was the most sensitive strain for all the active venoms. Micrurus baliocoryphus venom was unable to inhibit the growth of Enterococcus faecalis. Neither colubrid snake nor lepidopteran venoms exhibited activity on any bacterial strain tested. The snake venoms exhibiting antibacterial activity showed distinctive protein profiles by SDS-PAGE, highlighting that we could reveal for the first time the main protein families which may be thought to contribute to the antibacterial activity of M. baliocoryphus venom. This study paves the way to search for new antibacterial agents from Argentinian snake venoms, which may be a further opportunity to give an added value to the local biodiversity.

Venom characterization of the Brazilian Pampa snake Bothrops pubescens by top-down and bottom-up proteomics

The envenomation from the Bothrops genus is characterized by systemic and local effects caused by the main toxin families in the venom. In Bothrops pubescens venom we were able to identify 89 protein groups belonging to 13 toxin families with the bottom-up proteomics approach and 40 unique proteoforms belonging to 6 toxin families with the top-down proteomics approach. We also identified multi-proteoform complexes of dimeric L-amino acid oxidase using native top-down mass spectrometry.

Bibliometric Analysis of Literature in Snake Venom-Related Research Worldwide (1933-2022)

Snake envenomation is a severe economic and health concern affecting countries worldwide. Snake venom carries a wide variety of small peptides and proteins with various immunological and pharmacological properties. A few key research areas related to snake venom, including its applications in treating cancer and eradicating antibiotic-resistant bacteria, have been gaining significant attention in recent years. The goal of the current study was to analyze the global profile of literature in snake venom research. This study presents a bibliometric review of snake venom-related research documents indexed in the Scopus database between 1933 and 2022. The overall number of documents published on a global scale was 2999, with an average annual production of 34 documents. Brazil produced the highest number of documents (n = 729), followed by the United States (n = 548), Australia (n = 240), and Costa Rica (n = 235). Since 1963, the number of publications has been steadily increasing globally. At a worldwide level, antivenom, proteomics, and transcriptomics are growing hot issues for research in this field. The current research provides a unique overview of snake venom research at global level from 1933 through 2022, and it may be beneficial in guiding future research.

Immunorecognition and Neutralization of Crotalus durissus cumanensis Venom by a Commercial Antivenom Produced in Colombia

In Colombia, on average 2.9% of the nearly 5600 snakebite events that occur annually involve the rattlesnake Crotalus durissus cumanensis. The envenomation by this snake is mainly characterized by neurotoxicity and the main toxin is crotoxin (~64.7% of the total venom). The Instituto Nacional de Salud (INS) produces a polyvalent antivenom aimed at the treatment of bothropic, crotalid, and lachesic envenomations; nonetheless, its immune reactivity profile and neutralizing capacity over biological activities of the C. d. cumanensis venom has been poorly evaluated. In this sense, the study aims: (1) to describe an in-depth exploration of its immunoreactivity through second-generation antivenomics and HPLC fraction-specific ELISA immunoprofiles; and (2) to evaluate the neutralization pattern of the rattlesnake venom in vitro and in vivo biological activities. The results obtained showed a variable recognition of crotoxin subunits, in addition to a molecular mass-dependent immunoreactivity pattern in which the disintegrins were not recognized, and snake venom metalloproteinases and L-amino acid oxidases were the most recognized. Additionally, a high neutralization of proteolytic and coagulant activities was observed, but not over the PLA₂ activity. Further, the median effective dose against C. d. cumanensis venom lethality was 962 μL of antivenom per mg of venom. In conclusion, (1) the antivenom recognition over the crotoxin and the disintegrins of the C. d. cumanensis should be improved, thus aiming upcoming efforts for the exploration of new techniques and approaches in antivenom production in Colombia, and (2) the neutralization activity of the antivenom seems to follow the molecular mass-dependent recognition pattern, although other explanations should be explored.

Bothrops Jararaca Snake Venom Modulates Key Cancer-Related Proteins in Breast Tumor Cell Lines

Cancer is characterized by the development of abnormal cells that divide in an uncontrolled way and may spread into other tissues where they may infiltrate and destroy normal body tissue. Several previous reports have described biochemical anti-tumorigenic properties of crude snake venom or its components, including their capability of inhibiting cell proliferation and promoting cell death. However, to the best of our knowledge, there is no work describing cancer cell proteomic changes following treatment with snake venoms. In this work we describe the quantitative changes in proteomics of MCF7 and MDA-MB-231 breast tumor cell lines following treatment with Bothrops jararaca snake venom, as well as the functional implications of the proteomic changes. Cell lines were treated with sub-toxic doses at either 0.63 μg/mL (low) or 2.5 μg/mL (high) of B. jararaca venom for 24 h, conditions that cause no cell death per se. Proteomics analysis was conducted on a nano-scale liquid chromatography coupled on-line with mass spectrometry (nLC-MS/MS). More than 1000 proteins were identified and evaluated from each cell line treated with either the low or high dose of the snake venom. Protein profiling upon venom treatment showed differential expression of several proteins related to cancer cell metabolism, immune response, and inflammation. Among the identified proteins we highlight histone H3, SNX3, HEL-S-156an, MTCH2, RPS, MCC2, IGF2BP1, and GSTM3. These data suggest that sub-toxic doses of B. jararaca venom have potential to modulate cancer-development related protein targets in cancer cells. This work illustrates a novel biochemical strategy to identify therapeutic targets against cancer cell growth and survival.

From birth to adulthood: An analysis of the Brazilian lancehead (Bothrops moojeni) venom at different life stages

The Brazilian lancehead (Bothrops moojeni) has a wide distribution in Brazil and represents a serious public health hazard. Previous works reported that the symptoms of snakebites caused by B. moojeni juveniles’ bites were mainly related to coagulation, while those caused by adults’ bites had a more prominent local damage. In this work, we analyzed the venoms of B. moojeni at different life stages to better understand the ontogeny shift in this species. Snakes were grouped by age and sex, and venom pools were formed accordingly. Compositional analyses by one-dimensional electrophoresis (1-DE), chromatography, and mass spectrometry revealed that ontogenetic changes might be mostly related to phospholipase A₂ (PLA₂) and metalloproteases. Regarding the venoms functional aspect, proteolytic, L-amino acid oxidase, PLA₂, and coagulant in vitro activities were assayed, but only the first and the last ones showed age-related changes, with the venom of snakes up to 1 year-old displaying lower proteolytic and higher coagulant activities, while those from 2 years-old onward presented the opposite relation. The venoms of 3 years-old snakes were exceptions to the compositional and functional pattern of adults as both venoms presented profiles similar to neonates. Sex-related differences were observed in specific groups and were not age-related. In vivo experiments (median lethal dose and hemorrhagic activity) were statistically similar between neonates and adults, however we verified that the adult venom killed mice faster comparing to the neonates. All venoms were mostly recognized by the antibothropic serum and displayed similar profiles to 1-DE in western blotting. In conclusion, the Brazilian lancehead venom showed ontogenetic shift in its composition and activities. Furthermore, this change occurred in snakes from 1 to 2 years-old, and interestingly the venom pools from 3 years-old snakes had particular characteristics, which highlights the importance of comprehensive studies to better understand venom variability.

Antimalarial activity of basic phospholipases A2 isolated from Paraguayan Bothrops diporus venom against Plasmodium falciparum

Malaria is a parasitic infectious disease and was responsible for 400.000 deaths in 2018. Plasmodium falciparum represents the species that causes most human deaths due to severe malaria. In addition, studies prove the resistance of P. falciparum to drugs used to treat malaria, making the search for new drugs with antiplasmodial potential necessary. In this context, the literature describes snake venoms as a rich source of molecules with microbicidal potential, including phospholipases A2 (PLA2s). In this sense, the present study aimed to isolate basic PLA2s from Paraguayan Bothrops diporus venom and evaluate their antiplasmodial potential. Basic PLA2s were obtained using two chromatographic steps. Initially, B. diporus venom was subjected to ion exchange chromatography (IEC). The electrophoretic profile of the fractions from the IEC permitted the selection of 3 basic fractions, which were subjected to reverse phase chromatography, resulting in the isolation of the PLA2s. The toxins were tested for enzymatic activity using a chromogenic substrate and finally, the antiplasmodial, cytotoxic potential and hemolytic activity of the isolated toxins were evaluated. The electrophoretic profile of the fractions from the IEC permitted the selection of 3 basic fractions, which were subjected to reverse phase chromatography, resulting in the isolation of the two enzymatically active PLA2s, BdTX-I and BdTX-II and the PLA2 homologue BdTX-III. The antiplasmodial potential was evaluated and the toxins showed IC50 values of: 2.44, 0.0153 and 0.59 μg/mL respectively, presenting PLA2 selectivity according to the selectivity index results (SI) calculated against HepG2 cells. The results show that the 3 basic phospholipases isolated in this study have a potent antiparasitic effect against the W2 strain of P. falciparum. In view of the results obtained in this work, further research are necessary to determine the mechanism of action by which these toxins cause cell death in parasites.

Geographic variation of individual venom profile of Crotalus durissus snakes

Background:

South American rattlesnakes are represented in Brazil by a single species, Crotalus durissus, which has public health importance due to the severity of its envenomation and to its wide geographical distribution. The species is subdivided into several subspecies, but the current classification is controversial. In Brazil, the venoms of C. d. terrificus and C. d. collilineatus are used for hyperimmunization of horses for antivenom production, even though the distinction of these two subspecies are mostly by their geographical distribution. In this context, we described a comparative compositional and functional characterization of individual C. d. collilineatus and C. d. terrificus venoms from three Brazilian states.

Methods:

We compared the compositional patterns of C. d. terrificus and C. d. collilineatus individual venoms by 1-DE and RP-HPLC. For functional analyzes, the enzymatic activities of PLA₂, LAAO, and coagulant activity were evaluated. Finally, the immunorecognition of venom toxins by the crotalic antivenom produced at Butantan Institute was evaluated using Western blotting.

Results:

The protein profile of individual venoms from C. d. collilineatus and C. d. terrificus showed a comparable overall composition, despite some intraspecific variation, especially regarding crotamine and LAAO. Interestingly, HPLC analysis showed a geographic pattern concerning PLA₂. In addition, a remarkable intraspecific variation was also observed in PLA₂, LAAO and coagulant activities. The immunorecognition pattern of individual venoms from C. d. collilineatus and C. d. terrificus by crotalic antivenom produced at Butantan Institute was similar.

Conclusions:

The results highlighted the individual variability among the venoms of C. durissus ssp. specimens. Importantly, our data point to a geographical variation of C. durissus ssp. venom profile, regardless of the subspecies, as evidenced by PLA₂ isoforms complexity, which may explain the increase in venom neurotoxicity from Northeastern through Southern Brazil reported for the species.

Sexual and ontogenetic variation of Bothrops leucurus venom

Various factors, such as geographical origin, climate, sex, age and diet can influence the composition and pathophysiological activities of snake venoms. In this study, we examined the sexual and ontogenetic variations in the venom of Bothrops leucurus, a pitviper responsible for more than 80% of the snakebites in the state of Bahia, northeastern Brazilian. The venoms of 31 snakes were pooled according to sex and age (young, adult and old) and screened by SDS-PAGE (in reducing and non-reducing conditions), reverse-phase high performance liquid chromatography (RP-HPLC), gelatin zymography, and immunoblotting with therapeutic bothropic antivenom (BAV) from the Instituto Butantan. The electrophoretic and chromatographic profiles showed intraspecific ontogenetic variation, whereas sexual variations were less evident. All venoms showed gelatinolytic activity associated with 50-75 kDa protein bands. In addition, all venoms, regardless of the snakes' sex and age, cross-reacted to similar extents with BAV. Our findings show that B. leucurus venom changes during ontogenetic development and demonstrate sexual differences in its composition, indicating differences in biological activity.

Venoms and Isolated Toxins from Snakes of Medical Impact in the Northeast Argentina: State of the Art. Potential Pharmacological Applications

Among the ophidians that inhabit the Northeast of Argentina, the genus Bothrops such as B. alternatus and B. diporus species (also known as yararás) and Crotalus durisus terrificus (named cascabel), represent the most studied snake venom for more than thirty years. These two genera of venomous snakes account for the majority of poisonous snake envenomations and therefore, constitute a medical emergency in this region. This review presents a broad description of the compiled knowledge about venomous snakebite: its pathophysiological action, protein composition, isolated toxins, toxin synergism, toxin-antitoxin cross-reaction assays. Properties of some isolated toxins support a potential pharmacological application.

Preliminary Biochemical and Venomic Characterization of the Venom of Phalotris lemniscatus (Serpentes, Colubridae)

BACKGROUND

For many decades, research on snake venom toxinology focused mainly on the venoms of Viperidae and Elapidae species, which were traditionally the only ones considered as venomous. However, much less interest has been given to the venom produced by opisthoglyphous colubrid snakes, since they were typically considered of no clinical relevance.

OBJECTIVE

The aim of this work is to perform a preliminary biochemical and venomic characterization of the venom of the colubrid snake Phalotris lemniscatus, a species that has been responsible for two relevant cases of envenomation in Uruguay.

METHODS

We extracted venom from collected specimens and performed different biochemical and proteomic assays to understand its toxin composition.

RESULTS

We found that the venom of P. lemniscatus is composed of protein families typically present in snake venoms, such as metallo and serine preoteases, L-amino acid oxidases, phospholipases A2s, Ctype lectines-like, Kunitz-type proteins and three-finger toxins. Activity assays demonstrated a highly active gelatinolytic component as well as a potent capability to induce blood coagulation.

CONCLUSION

The results indicate that the venom of P. lemniscatus contains hemotoxic activities and components that resemble those found in Viperidae (Bothrops) snakes and that can induce a clinically relevant accident. Further studies are needed to better understand the venom composition of this colubrid snake and its most active compounds.

Antibody Cross-Reactivity in Antivenom Research

Antivenom cross-reactivity has been investigated for decades to determine which antivenoms can be used to treat snakebite envenomings from different snake species. Traditionally, the methods used for analyzing cross-reactivity have been immunodiffusion, immunoblotting, enzyme-linked immunosorbent assay (ELISA), enzymatic assays, and in vivo neutralization studies. In recent years, new methods for determination of cross-reactivity have emerged, including surface plasmon resonance, antivenomics, and high-density peptide microarray technology. Antivenomics involves a top-down assessment of the toxin-binding capacities of antivenoms, whereas high-density peptide microarray technology may be harnessed to provide in-depth knowledge on which toxin epitopes are recognized by antivenoms. This review provides an overview of both the classical and new methods used to investigate antivenom cross-reactivity, the advantages and disadvantages of each method, and examples of studies using the methods. A special focus is given to antivenomics and high-density peptide microarray technology as these high-throughput methods have recently been introduced in this field and may enable more detailed assessments of antivenom cross-reactivity.

Toxin-resolved antivenomics-guided assessment of the immunorecognition landscape of antivenoms

Snakebite envenoming represents a major issue in rural areas of tropical and subtropical regions across sub-Saharan Africa, South to Southeast Asia, Latin America and Oceania. Antivenoms constitute the only scientifically validated therapy for snakebite envenomings, provided they are safe, effective, affordable, accessible and administered appropriately. However, the lack of financial incentives in a technology that has remained relatively unchanged for more than a century, has contributed to some manufacturers leaving the market and others downscaling production or increasing the prices, leading to a decline in the availability and accessibility for these life-saving antidotes to millions of rural poor most at risk from snakebites in low income countries. The shortage of antivenoms can be significantly alleviated by optimizing the use of current antivenoms (through the assessment of their specific and paraspecific efficacy against the different medically relevant homologous and heterologous snake venoms) and by generating novel polyspecific antivenoms exhibiting broad clinical spectrum and wide geographic distribution range. Research on venoms has been continuously enhanced by advances in technology. Particularly, the last decade has witnessed the development of omics strategies for unravelling the toxin composition of venoms ("venomics") and to assess the immunorecognition profile of antivenoms ("antivenomics"). Here, we review recent developments and reflect on near future innovations that promise to revolutionize the mutually enlightening relationship between evolutionary and translational venomics.

A Review and Database of Snake Venom Proteomes

Advances in the last decade combining transcriptomics with established proteomics methods have made possible rapid identification and quantification of protein families in snake venoms. Although over 100 studies have been published, the value of this information is increased when it is collated, allowing rapid assimilation and evaluation of evolutionary trends, geographical variation, and possible medical implications. This review brings together all compositional studies of snake venom proteomes published in the last decade. Compositional studies were identified for 132 snake species: 42 from 360 (12%) Elapidae (elapids), 20 from 101 (20%) Viperinae (true vipers), 65 from 239 (27%) Crotalinae (pit vipers), and five species of non-front-fanged snakes. Approximately 90% of their total venom composition consisted of eight protein families for elapids, 11 protein families for viperines and ten protein families for crotalines. There were four dominant protein families: phospholipase A₂s (the most common across all front-fanged snakes), metalloproteases, serine proteases and three-finger toxins. There were six secondary protein families: cysteine-rich secretory proteins, l-amino acid oxidases, kunitz peptides, C-type lectins/snaclecs, disintegrins and natriuretic peptides. Elapid venoms contained mostly three-finger toxins and phospholipase A₂s and viper venoms metalloproteases, phospholipase A₂s and serine proteases. Although 63 protein families were identified, more than half were present in <5% of snake species studied and always in low abundance. The importance of these minor component proteins remains unknown.

Preclinical Evaluation of the Efficacy of Antivenoms for Snakebite Envenoming: State-of-the-Art and Challenges Ahead

Animal-derived antivenoms constitute the mainstay in the therapy of snakebite envenoming. The efficacy of antivenoms to neutralize toxicity of medically-relevant snake venoms has to be demonstrated through meticulous preclinical testing before their introduction into the clinical setting. The gold standard in the preclinical assessment and quality control of antivenoms is the neutralization of venom-induced lethality. In addition, depending on the pathophysiological profile of snake venoms, the neutralization of other toxic activities has to be evaluated, such as hemorrhagic, myotoxic, edema-forming, dermonecrotic, in vitro coagulant, and defibrinogenating effects. There is a need to develop laboratory assays to evaluate neutralization of other relevant venom activities. The concept of the 3Rs (Replacement, Reduction, and Refinement) in Toxinology is of utmost importance, and some advances have been performed in their implementation. A significant leap forward in the study of the immunological reactivity of antivenoms against venoms has been the development of “antivenomics”, which brings the analytical power of mass spectrometry to the evaluation of antivenoms. International partnerships are required to assess the preclinical efficacy of antivenoms against snake venoms in different regions of the world in order to have a detailed knowledge on the neutralizing profile of these immunotherapeutics.