N-glycome profiling of Bothrops jararaca newborn and adult venoms

Profilings of subproteomes of lectin-binding proteins of nine Bothrops venoms reveal variability driven by different glycan types

Snake venom proteomes have long been investigated to explore a multitude of biologically active components that are used for prey capture and defense, and are involved in the pathological effects observed upon mammalian envenomation. Glycosylation is a major protein post-translational modification in venoms and contributes to the diversification of proteomes. We have shown that Bothrops venoms are markedly defined by their content of glycoproteins, and that most N-glycan structures of eight Bothrops venoms contain sialic acid, while bisected N-acetylglucosamine was identified in Bothrops cotiara venom. To further investigate the mechanisms involved in the generation of different venoms by related snakes, here the glycoproteomes of nine Bothrops venoms (Bothrops atrox, B. cotiara, Bothrops erythromelas, Bothrops fonsecai, B. insularis, Bothrops jararaca, Bothrops jararacussu, Bothrops moojeni and Bothrops neuwiedi) were comparatively analyzed by enrichment with three lectins of different specificities, recognizing bisecting N-acetylglucosamine- and sialic acid-containing glycoproteins, and mass spectrometry. The lectin capture strategy generated venom fractions enriched with several glycoproteins, including metalloprotease, serine protease, and L- amino acid oxidase, in addition to various types of low abundant enzymes. The different contents of lectin-enriched proteins underscore novel aspects of the variability of the glycoprotein subproteomes of Bothrops venoms and point to the role of distinct types of glycan chains in generating different venoms by closely related snake species.

The distinct N-terminomes of Bothrops jararaca newborn and adult venoms

Using approaches of transcriptomics and proteomics we have shown that the phenotype of Bothrops jararaca venom undergoes a significant rearrangement upon neonate to adult transition. Most regulatory processes in biology are intrinsically related to modifications of protein structure, function, and abundance. However, it is unclear to which extent intrinsic proteolysis affects toxins and snake venom phenotypes upon ontogenesis. Here we assessed the natural N-terminome of Bothrops jararaca newborn and adult venoms and explored the degree of N-terminal protein truncation in ontogenetic-based proteome variation. To this end we applied the Terminal Amine Isotopic Labeling of Substrates (TAILS) technology to characterize venom collected in the presence of proteinase inhibitors. We identified natural N-terminal sequences in the newborn (71) and adult (84) venoms, from which only 37 were common to both. However, truncated toxins were found in higher number in the newborn (212) than in the adult (140) venom. Moreover, sequences N-terminally blocked by pyroglutamic acid were identified in the newborn (55) and adult (49) venoms. Most toxin classes identified by their natural N-terminal sequences showed a similar number of unique peptides in the newborn and adult venoms, however, those of serine proteinases and C-type lectins were more abundant in the adult venom. Truncated sequences from at least ten toxin classes were detected, however the catalytic and cysteine-rich domains of metalloproteinases were the most prone to proteolysis, mainly in the newborn venom. Our results underscore the pervasiveness of truncations in most toxin classes and highlight variable post-translational events in newborn and adult venoms.

Venom Profiling of the Insular Species Bothrops alcatraz: Characterization of Proteome, Glycoproteome, and N-Terminome Using Terminal Amine Isotopic Labeling of Substrates

Bothrops alcatraz, a species endemic to Alcatrazes Islands, is regarded as critically endangered due to its small area of occurrence and the declining quality of its habitat. We recently reported the identification of N-glycans attached to toxins of Bothrops species, showing similar compositions in venoms of the B. jararaca complex (B. jararaca, B. insularis, and B. alcatraz). Here, we characterized B. alcatraz venom using electrophoretic, proteomic, and glycoproteomic approaches. Electrophoresis showed that B. alcatraz venom differs from B. jararaca and B. insularis; however, N-glycan removal revealed similarities between them, indicating that the occupation of N-glycosylation sites contributes to interspecies variability in the B. jararaca complex. Metalloproteinase was the major toxin class identified in the B. alcatraz venom proteome followed by serine proteinase and C-type lectin, and overall, the adult B. alcatraz venom resembles that of B. jararaca juvenile specimens. The comparative glycoproteomic analysis of B. alcatraz venom with B. jararaca and B. insularis indicated that there may be differences in the utilization of N-glycosylation motifs among their different toxin classes. Furthermore, we prospected for the first time the N-terminome of a snake venom using the terminal amine isotopic labeling of substrates (TAILS) approach and report the presence of ∼30% of N-termini corresponding to truncated toxin forms and ∼37% N-terminal sequences blocked by pyroglutamic acid in B. alcatraz venom. These findings underscore a low correlation between venom gland transcriptomes and proteomes and support the view that post-translational processes play a major role in shaping venom phenotypes.

The sweet side of venom: Glycosylated prothrombin activating metalloproteases from Dispholidus typus (boomslang) and Thelotornis mossambicanus (twig snake)

Dispholidus typus and Thelotornis mossambicanus are closely related rear-fanged colubrid snakes that both possess strongly procoagulant venoms. However, despite similarities in overall venom biochemistry and resulting clinical manifestations, the underlying venom composition differs significantly between the two species. As a result, the only available antivenom-which is a monovalent antivenom for D. typus-has minimal cross reactivity with T. mossambicanus and is not a clinically viable option. It was hypothesised that this lack of cross reactivity is due to the additional large metalloprotease protein within T. mossambicanus venom, which may also be responsible for faster coagulation times. In this study, we found that T. mossambicanus venom is a more powerful activator of prothrombin than that of D. typus and that the SVMP transcripts from T. mossambicanus form a clade with those from D. typus. The sequences from D. typus and T. mossambicanus were highly similar in length, with the calculated molecular weights of the T. mossambicanus transcripts being significantly less than the molecular weights of some isoforms on the 1D SDS-PAGE gels. Analyses utilising degylcosylating enzymes revealed that T. mossambicanus SVMPs are glycosylated during post-translational modification, but that this does not lead to the different molecular weight bands observed in 1D SDS-PAGE gels. However, differences in glycosylation patterns may still explain some of the difference between the enzymatic activities and neutralization by antivenom that have been observed in these venoms. The results of this study provide new information regarding the treatment options for patients envenomated by T. mossambicanus as well as the evolution of these dangerous snakes.

Naja annulifera Snake: New insights into the venom components and pathogenesis of envenomation

Background

Naja annulifera is a medically important venomous snake occurring in some of the countries in Sub-Saharan Africa. Accidental bites result in severe coagulation disturbances, systemic inflammation and heart damage, as reported in dogs, and death, by respiratory arrest, in humans. Despite the medical importance of N. annulifera, little is known about its venom composition and the pathogenesis of envenomation. In this paper, the toxic, inflammatory and immunogenic properties of N. annulifera venom were analyzed.

Methodology/Principal findings

Venom proteomic analysis identified 79 different proteins, including Three Finger Toxins, Cysteine Rich Secretory Proteins, Metalloproteinases, Phospholipases A₂ (PLA₂), Hyaluronidase, L-amino-acid oxidase, Cobra Venom Factor and Serine Proteinase. The presence of PLA₂, hyaluronidase, fibrinogenolytic and anticoagulant activities was detected using functional assays. The venom was cytotoxic to human keratinocytes. In an experimental murine model of envenomation, it was found that the venom induced local changes, such as swelling, which was controlled by anti-inflammatory drugs. Moreover, the venom caused death, which was preceded by systemic inflammation and pulmonary hemorrhage. The venom was shown to be immunogenic, inducing a strong humoral immune response, with the production of antibodies able to recognize venom components with high molecular weight and to neutralize its lethal activity.

Conclusions/Significance

The results obtained in this study demonstrate that N. annulifera venom contains toxins able to induce local and systemic inflammation, which can contribute to lung damage and death. Moreover, the venom is immunogenic, an important feature that must be considered during the production of a therapeutic anti-N. annulifera antivenom.

N. annulifera is a dangerous snake that belongs to the Elapidae family. It is found in some of the countries in Sub-Saharan Africa and has caused accidents in humans and dogs. In this study, we characterized some of the biochemical, toxic and immunogenic properties of N. annulifera venom. We showed that the venom is composed of several proteins, some of which display enzymatic activities, such as phospholipase A2, hyaluronidase, metalloproteinases and serine proteinases. The venom promoted disturbances in the human coagulation system and was cytotoxic to human epidermal cells. Using a mouse model, we showed that the venom promotes local reactions that were reduced with anti-inflammatory drugs. The venom caused systemic inflammation, lung hemorrhage and death. Further, the venom stimulated production of high antibody titers when injected into mice and the antiserum produced was able to inhibit venom-induced death. This study demonstrated that N. annulifera venom contains toxins that trigger inflammatory process, which may contribute to the envenomation pathology. Moreover, the venom is immunogenic, an important aspect for the production of an efficient N. annulifera antivenom.

Organic and Peptidyl Constituents of Snake Venoms: The Picture Is Vastly More Complex Than We Imagined

Small metabolites and peptides in 17 snake venoms (Elapidae, Viperinae, and Crotalinae), were quantified using liquid chromatography-mass spectrometry. Each venom contains >900 metabolites and peptides. Many small organic compounds are present at levels that are probably significant in prey envenomation, given that their known pharmacologies are consistent with snake envenomation strategies. Metabolites included purine nucleosides and their bases, neurotransmitters, neuromodulators, guanidino compounds, carboxylic acids, amines, mono- and disaccharides, and amino acids. Peptides of 2⁻15 amino acids are also present in significant quantities, particularly in crotaline and viperine venoms. Some constituents are specific to individual taxa, while others are broadly distributed. Some of the latter appear to support high anabolic activity in the gland, rather than having toxic functions. Overall, the most abundant organic metabolite was citric acid, owing to its predominance in viperine and crotaline venoms, where it chelates divalent cations to prevent venom degradation by venom metalloproteases and damage to glandular tissue by phospholipases. However, in terms of their concentrations in individual venoms, adenosine, adenine, were most abundant, owing to their high titers in Dendroaspis polylepis venom, although hypoxanthine, guanosine, inosine, and guanine all numbered among the 50 most abundant organic constituents. A purine not previously reported in venoms, ethyl adenosine carboxylate, was discovered in D. polylepis venom, where it probably contributes to the profound hypotension caused by this venom. Acetylcholine was present in significant quantities only in this highly excitotoxic venom, while 4-guanidinobutyric acid and 5-guanidino-2-oxopentanoic acid were present in all venoms.

Structures of N-Glycans of Bothrops Venoms Revealed as Molecular Signatures that Contribute to Venom Phenotype in Viperid Snakes

The complexity of snake venoms has long been investigated to explore a myriad of biologically active proteins and peptides that are used for immobilizing or killing prey, and are responsible for the pathological effects observed on envenomation. Glycosylation is the main post-translational modification (PTM) of viperid venoms but currently there is little understanding of how protein glycosylation impacts the variation of venom proteomes. We have previously reported that Bothrops venom glycoproteomes contain a core of components that markedly define their composition and parallel their phylogenetic classification. Here we extend those observations to eight Bothrops species evaluating the N-glycomes by LC-MS as assigned cartoon structures and detailing those structures separately as methylated analogs using ion-trap mass spectrometry (MSⁿ). Following ion disassembly through multiple steps provided sequence and linkage isomeric details that characterized 52 unique compositions in Bothrops venoms. These occurred as 60 structures, of which 26 were identified in the venoms of the Jararaca Complex (B. alcatraz, B. insularis, and B. jararaca), 20 in B. erythromelas, B. jararacussu, B. moojeni and B. neuwiedi venoms, and 22 in B. cotiara venom. Further, quantitative analysis of these N-glycans showed variable relative abundances in the venoms. For the first time a comprehensive set of N-glycan structures present in snake venoms are defined. Despite the fact that glycosylation is not template-defined, the N-glycomes of these venoms mirror the phylogeny cladograms of South American bothropoid snakes reported in studies on morphological, molecular data and feeding habits, exhibiting distinct molecular signatures for each venom. Considering the complexity of N-glycan moieties generally found in glycoproteins, characterized by different degrees of branching, isomer structures, and variable abundances, our findings point to these factors as another level of complexity in Bothrops venoms, features that could dramatically contribute to their distinct biological activities.

Enzymatic and Pro-Inflammatory Activities of Bothrops lanceolatus Venom: Relevance for Envenomation

Bothrops lanceolatus, commonly named ‘Fer-de-Lance’, is an endemic snake of the French Caribbean Island of Martinique. Envenomations by B. lanceolatus present clinical aspects characterized by systemic thrombotic syndrome and important local inflammation, involving edema and pain but limited hemorrhage. To investigate mechanisms of venom-induced inflammation, B. lanceolatus venom was characterized, its cross-reactivity with bothropic antivenom explored, its cytotoxicity on human keratinocytes and vascular cells, and the production of cytokines and chemokines were analyzed. We used electrophoretic separation, zymography, colorimetric or fluorimetric enzymatic assays, and immunochemical assays. Therapeutic South American bothropic antivenom cross-reacted with B. lanceolatus venom and completely or partially abolished its PLA2, hyaluronidase, and proteolytic activities, as well as its cytotoxicity for keratinocytes. The substrate specificity of B. lanceolatus venom proteases was emphasized. B. lanceolatus venom cytotoxicity was compared to the B. jararaca venom. Both venoms were highly cytotoxic for keratinocytes (HaCaT), whereas B. lanceolatus venom showed particularly low toxicity for endothelial cells (EAhy926). Patterns of cytokine and chemokine production by cells exposed to the venoms were highly pro-inflammatory. Thus, the results presented here show that B. lanceolatus venom toxins share important antigenic similarities with South American Bothrops species toxins, although their proteases have acquired particular substrate specificity. Moreover, the venom displays important cytotoxic and pro-inflammatory action on human cell types such as keratinocytes and endothelial cells, which are important players in the local and systemic compartments affected by the envenomation.

Processing of Snake Venom Metalloproteinases: Generation of Toxin Diversity and Enzyme Inactivation

Snake venom metalloproteinases (SVMPs) are abundant in the venoms of vipers and rattlesnakes, playing important roles for the snake adaptation to different environments, and are related to most of the pathological effects of these venoms in human victims. The effectiveness of SVMPs is greatly due to their functional diversity, targeting important physiological proteins or receptors in different tissues and in the coagulation system. Functional diversity is often related to the genetic diversification of the snake venom. In this review, we discuss some published evidence that posit that processing and post-translational modifications are great contributors for the generation of functional diversity and for maintaining latency or inactivation of enzymes belonging to this relevant family of venom toxins.

Proteomic characterization and comparison of venoms from two elapid snakes (Bungarus multicinctus and Naja atra) from China

Bungarus multicinctus (many-banded krait) and Naja atra (Chinese cobra) are widely distributed and medically important venomous snakes in China; however, their venom proteomic profiles have not been fully compared. Here, we fractionated crude venoms and analyzed them using a combination of proteomic techniques. Three-finger toxins (3-FTx) and phospholipase A2 (PLA2) were most abundant in both species, respectively accounting for 32.6% and 66.4% of total B. multicinctus venom, and 84.3% and 12.2% of total N. atra venom. Venoms from these two species contained one common protein family and six less abundant species-specific protein families. The proteomic profiles of B. multicinctus and N. atra venoms and analysis of toxicological activity in mice suggested that 3-FTx and PLA2 are the major contributors to clinical symptoms caused by envenomation. The venoms differed in enzymatic activity, likely the result of inter-specific variation in the amount of related venom components. Antivenomics assessment revealed that a small number of venom components (3-FTxs and PLA2s in B. multicinctus, and 3-FTxs in N. atra) could not be immunocaptured completely, suggesting that we should pay attention to enhancing the immune response of these components in designing commercial antivenoms for B. multicinctus and N. atra.

BIOLOGICAL SIGNIFICANCE

The proteomic profiles of venoms from two medically important snake species - B. multicinctus and N. atra - have been explored. Quantitative and qualitative differences are evident in both venoms when proteomic profiles and transcriptomic results are compared; this is a reminder that combined approaches are needed to explore the precise composition of snake venom. Two protein families (3-FTx and PLA2) of high abundance in these snake venoms are major players in the biochemical and pharmacological effects of envenomation. Elucidation of the proteomic profiles of these snake venoms is helpful in understanding composition-function relationships and will facilitate the clinical application of antivenoms.

Combined venomics, venom gland transcriptomics, bioactivities, and antivenomics of two Bothrops jararaca populations from geographic isolated regions within the Brazilian Atlantic rainforest

Bothrops jararaca is a slender and semi-arboreal medically relevant pit viper species endemic to tropical and subtropical forests in southern Brazil, Paraguay, and northern Argentina (Misiones). Within its geographic range, it is often abundant and is an important cause of snakebite. Although no subspecies are currently recognized, geographic analyses have revealed the existence of two well-supported B. jararaca clades that diverged during the Pliocene ~3.8Mya and currently display a southeastern (SE) and a southern (S) Atlantic rainforest (Mata Atlântica) distribution. The spectrum, geographic variability, and ontogenetic changes of the venom proteomes of snakes from these two B. jararaca phylogroups were investigated applying a combined venom gland transcriptomic and venomic analysis. Comparisons of the venom proteomes and transcriptomes of B. jararaca from the SE and S geographic regions revealed notable interpopulational variability that may be due to the different levels of population-specific transcriptional regulation, including, in the case of the southern population, a marked ontogenetic venom compositional change involving the upregulation of the myotoxic PLA2 homolog, bothropstoxin-I. This population-specific marker can be used to estimate the proportion of venom from the southern population present in the B. jararaca venom pool used for the Brazilian soro antibotrópico (SAB) antivenom production. On the other hand, the southeastern population-specific D49-PLA2 molecules, BinTX-I and BinTX-II, lend support to the notion that the mainland ancestor of Bothrops insularis was originated within the same population that gave rise to the current SE B. jararaca phylogroup, and that this insular species endemic to Queimada Grande Island (Brazil) expresses a pedomorphic venom phenotype. Mirroring their compositional divergence, the two geographic B. jararaca venom pools showed distinct bioactivity profiles. However, the SAB antivenom manufactured in Vital Brazil Institute neutralized the lethal effect of both venoms to a similar extent. In addition, immobilized SAB antivenom immunocaptured most of the venom components of the venoms of both B. jararaca populations, but did not show immunoreactivity against vasoactive peptides. The Costa Rican bothropic-crotalic-lachesic (BCL) antivenom showed the same lack of reactivity against vasoactive peptides but, in addition, was less efficient immunocapturing PI- and PIII-SVMPs from the SE venom, and bothropstoxin-I, a CRISP molecule, and a D49-PLA2 from the venom of the southern B. jararaca phylogroup. The remarkable paraspecificity exhibited by the Brazilian and the Costa Rican antivenoms indicates large immunoreactive epitope conservation across the natural history of Bothrops, a genus that has its roots in the middle Miocene. This article is part of a Special Issue entitled: Omics Evolutionary Ecolog.

Biochemical and functional characterization of Parawixia bistriata spider venom with potential proteolytic and larvicidal activities

Toxins purified from the venom of spiders have high potential to be studied pharmacologically and biochemically. These biomolecules may have biotechnological and therapeutic applications. This study aimed to evaluate the protein content of Parawixia bistriata venom and functionally characterize its proteins that have potential for biotechnological applications. The crude venom showed no phospholipase, hemorrhagic, or anti-Leishmania activities attesting to low genotoxicity and discrete antifungal activity for C. albicans. However the following activities were observed: anticoagulation, edema, myotoxicity and proteolysis on casein, azo-collagen, and fibrinogen. The chromatographic and electrophoretic profiles of the proteins revealed a predominance of acidic, neutral, and polar proteins, highlighting the presence of proteins with high molecular masses. Five fractions were collected using cation exchange chromatography, with the P4 fraction standing out as that of the highest purity. All fractions showed proteolytic activity. The crude venom and fractions P1, P2, and P3 showed larvicidal effects on A. aegypti. Fraction P4 showed the presence of a possible metalloprotease (60 kDa) that has high proteolytic activity on azo-collagen and was inhibited by EDTA. The results presented in this study demonstrate the presence of proteins in the venom of P. bistriata with potential for biotechnological applications.