Yamileth Angulo

Venoms, venomics, antivenomics

Venoms comprise mixtures of peptides and proteins tailored by Natural Selection to act on vital systems of the prey or victim. Here we review our proteomic protocols for uncoiling the composition, immunological profile, and evolution of snake venoms. Our long‐term goal is to gain a deep insight of all viperid venom proteomes. Knowledge of the inter‐ and intraspecies ontogenetic, individual, and geographic venom variability has applied importance for the design of immunization protocols aimed at producing more effective polyspecific antivenoms. A practical consequence of assessing the cross‐reactivity of heterologous antivenoms is the possibility of circumventing the restricted availability of species‐specific antivenoms in some regions. Further, the high degree of target specificity makes toxins valuable scaffolds for drug development.

Snake venomics of the Central American rattlesnake Crotalus simus and the South American Crotalus durissus complex points to neurotoxicity as an adaptive paedomorphic trend along Crotalus dispersal in South America.

We report a comparative venomic and antivenomic characterization of the venoms of newborn and adult specimens of the Central American rattlesnake, Crotalus simus, and of the subspecies cumanensis, durissus, ruruima, and terrificus of South American Crotalus durissus. Neonate and adult C. simus share about 50% of their venom proteome. The venom proteome of 6-week-old C. simus is predominantly made of the neurotoxic heterodimeric phospholipase A(2) (PLA(2) crotoxin) (55.9%) and serine proteinases (36%), whereas snake venom Zn(2+)-metalloproteinases (SVMPs), exclusively of class PIII, represent only 2% of the total venom proteins. In marked contrast, venom from adult C. simus comprises toxins from 7 protein families. A large proportion (71.7%) of these toxins are SVMPs, two-thirds of which belong to the PIII class. These toxin profiles correlate well with the overall biochemical and pharmacological features of venoms from adult (hemorrhagic) and newborn (neurotoxic) C. simus specimens. The venoms of the South American Crotalus subspecies belong to one of two distinct phenotypes. C. d. cumanensis exhibits high levels of SVMPs and low lethal potency (LD(50)), whereas C. d. subspecies terrificus, ruruima, and durissus have low SVMP activity and high neurotoxicity to mice. Their overall toxin compositions explain the outcome of envenomation by these species. Further, in all C. simus and C. durissus venoms, the concentration of neurotoxins (crotoxin and crotamine) is directly related with lethal activity, whereas lethality and metalloproteinase activity show an inverse relationship. The similar venom toxin profiles of newborn C. simus and adult C. durissus terrificus, ruruima, and durissus subspecies strongly suggests that the South American taxa have retained juvenile venom characteristics in the adult form (paedomorphism) along their North-South stepping-stone dispersal. The driving force behind paedomorphism is often competition or predation pressure. The increased concentration of the neurotoxins crotoxin and crotamine in South American rattlesnake venoms strongly argues that the gain of neurotoxicity and lethal venom activities to mammals may have represented the key axis along which overall venom toxicity has evolved during Crotalus durissus invasion of South America. The paedomorphic trend is supported by a decreasing LNC (lethal neurotoxicity coefficient, defined as the ratio between the average LD(50) of the venom and the crotoxin + crotamine concentration) along the North-South axis, coincident with the evolutionary dispersal pattern of the Neotropical rattlesnakes. The indistinguisable immunoreactivity patterns of Costa Rican and Venezuelan polyvalent antivenoms toward C. simus and C. durissus venoms strongly suggest the possibility of using these antivenoms indistinctly for the management of snakebites by adult C. simus and by certain C. d. cumanensis populations exhibiting a hemorrhagic venom phenotype. The antivenomic results also explain why the antivenoms effectively neutralize the hemorrhagic activity of adult C. simus venoms but does not protect against adult C. durissus sp. and newborn C. simus envenomations. The identification of evolutionary trends among tropical Crotalus, as reported here, may have an impact in defining the mixture of venoms for immunization to produce an effective pan-American anti-Crotalus antivenom.

Myotoxic phospholipases A2 in Bothrops snake venoms: Effect of chemical modifications on the enzymatic and pharmacological properties of bothropstoxins from Bothrops jararacussu

Venoms from eight Bothrops spp. were fractionated by ion-exchange chromatography on CM-Sepharose at pH 8.0 for the purification of myotoxins. Chromatographic profiles showed differences regarding myotoxic components among these venoms. B. alternatus, B. atrox and B. jararaca venoms did not show the major basic myotoxic fractions identified in the other venoms. Polyacrylamide gel electrophoresis for basic proteins also showed distinct patterns for these toxins. In vivo, all the isolated myotoxins induced release of creatine kinase due to necrosis of muscle fibers, accompanied by polymorphonuclear cell infiltration, and edema in the mouse paw. In addition, the toxins showed cytotoxic and liposome-disrupting activities in vitro. B. jararacussu bothropstoxins-I (BthTX-I) and II (BthTX-II) were submitted to chemical modifications of: His, by 4-bromophenacyl bromide (BPB) or photooxidation by Rose Bengal (RB); Tyr, by 2-nitrobenzenesulphonyl fluoride (NBSF); and Trp, by o-nitrophenylsulphenyl chloride (NPSC). The myotoxic and cytotoxic activities of BthTX-I, a Lys49 PLA(2) homologue, after modification by BPB, RB, NBSF and NPSC, were reduced to 50%, 20%, 75%, 65% and 13%, 0.5%, 76%, 58%, respectively. However, the edema-inducing and liposome-disrupting activities were not significantly reduced by the above modifications. BPB-treated BthTX-II, an Asp49 PLA(2) homologue, lost most of its catalytic, indirect hemolytic, anticoagulant, myotoxic and cytotoxic activities. The edema-inducing and liposome-disrupting activities were reduced to 50% and 80%, respectively. Lethality caused by BthTX-I and -II was strongly reduced after treatment with BPB or RB, but only partially with NBSF or NPSC. BthTX-I and -II, both native or modified, migrated similarly in a charge-shift electrophoresis. Antibodies raised against BthTX-I or -II, B. asper Basp-II and the C-terminal 115-129 peptide from Basp-II did not show significant differences in their cross-reactivity with the modified toxins, except with RB photooxidized toxins.

Comparative study of the cytolytic activity of myotoxic phospholipases A2 on mouse endothelial (tEnd) and skeletal muscle (C2C12) cells in vitro

A rapid in vitro cytolytic effect of some myotoxic phospholipases A2 (PLA2s) isolated from the venoms of Viperidae snakes has been previously described. This study was undertaken to investigate if cytolytic activity is a common property of the myotoxic proteins from this group. Murine endothelial cells (tEnd) and skeletal muscle myotubes (C2C12) were utilized as targets. The release of lactic dehydrogenase was quantified as a measure of cell damage, 3 h after exposure of cells to the different PLA2s, including representatives from the genera Bothrops, Agkistrodon, Trimeresurus, Crotalus (family Viperidae), and Notechis (family Elapidae). All of the group II myotoxic PLA2s tested displayed rapid cytolytic activity when tested in the micromolar range of concentrations (8-32 microM). In contrast, the group I myotoxic PLA2 notexin was devoid of this activity. Aspartate-49 and lysine-49 PLA2 group II variants showed a comparable cytolytic effect. Skeletal muscle myotubes, obtained after fusion and differentiation of C2C12 myoblasts, were significantly more susceptible to the cytolytic action of myotoxins than endothelial cells, previously reported to be more susceptible than undifferentiated myoblasts under the same assay conditions. Cytolytic activity appears to be a common characteristic of group II myotoxic PLA2s of the Viperidae. Bee venom PLA2, a group III enzyme of known myotoxicity, also displayed cytotoxic activity on C2C12 myotubes, being devoid of activity on endothelial cells. These results suggest that in vitro differentiated skeletal muscle myotubes may represent a suitable model target for the study of myotoxic PLA2s of the structural group II found in snake venoms.

Snake Venomics and Antivenomics of the Arboreal Neotropical Pitvipers Bothriechis lateralis and Bothriechis schlegelii

We report the comparative proteomic characterization of the venoms of two related neotropical arboreal pitvipers from Costa Rica of the genus Bothriechis, B. lateralis (side-striped palm pit viper) and B. schlegelii (eyelash pit viper). The crude venoms were fractionated by reverse-phase HPLC, followed by analysis of each chromatographic fraction by SDS-PAGE, N-terminal sequencing, MALDI-TOF mass fingerprinting, and collision-induced dissociation tandem mass spectrometry of tryptic peptides. The venom proteomes of B. lateralis and B. schlegelii comprise similar number of distinct proteins belonging, respectively, to 8 and 7 protein families. The two Bothriechis venoms contain bradykinin-potentiating peptides (BPPs), and proteins from the phospholipase A 2 (PLA 2), serine proteinase, l-amino acid oxidase (LAO), cysteine-rich secretory protein (CRISP), and Zn (2+)-dependent metalloproteinase (SVMP) families, albeit each species exhibit different relative abundances. Each venom also contains unique components, for example, snake venom vascular endothelial growth factor (svVEGF) and C-type lectin-like molecules in B. lateralis, and Kazal-type serine proteinase inhibitor-like proteins in B. schlegelii. Using a similarity coefficient, we estimate that the similarity of the venom proteins between the two Bothriechis taxa may be <10%, indicating a high divergence in their venom compositions, in spite of the fact that both species have evolved to adapt to arboreal habits. The major toxin families of B. lateralis and B. schlegelii are SVMP (55% of the total venom proteins) and PLA 2 (44%), respectively. Their different venom toxin compositions provide clues for rationalizing the distinct signs of envenomation caused by B. schlegelii and B. lateralis. An antivenomic study of the immunoreactivity of the Instituto Clodomiro Picado (ICP) polyvalent antivenom toward Bothriechis venoms revealed that l-amino acid oxidase and SVMPs represent the major antigenic protein species in both venoms. Our results provide a ground for rationalizing the reported protection of the ICP polyvalent antivenom against the hemorrhagic, coagulant, defibrinating, caseinolytic and fibrin(ogen)olytic activities of Bothriechis ( schlegelii, lateralis) venoms. However, these analyses also evidenced the limited recognition capability of the polyvalent antivenom toward a number of Bothriechis venom components, predominantly BPPs, svVEGF, Kazal-type inhibitors, some PLA 2 proteins, some serine proteinases, and CRISP molecules.

Snake population venomics and antivenomics of Bothrops atrox: Paedomorphism along its transamazonian dispersal and implications of geographic venom variability on snakebite management

We describe two geographically differentiated venom phenotypes across the wide distribution range of Bothrops atrox, from the Colombian Magdalena Medio Valley through Puerto Ayacucho and El Paují, in the Venezuelan States of Amazonas and Orinoquia, respectively, and São Bento in the Brazilian State of Maranhão. Colombian and Venezuelan venoms show an ontogenetic toxin profile phenotype whereas Brazilian venoms exhibit paedomorphic phenotypes. Venoms from each of the 16 localities sampled contain both population-specific toxins and proteins shared by neighboring B. atrox populations. Mapping the molecular similarity between conspecific populations onto a physical map of B. atrox range provides clues for tracing dispersal routes that account for the current biogeographic distribution of the species. The proteomic pattern is consistent with a model of southeast and southwest dispersal and allopatric fragmentation northern of the Amazon Basin, and trans-Amazonian expansion through the Andean Corridor and across the Amazon river between Monte Alegre and Santarém. An antivenomic approach applied to assess the efficacy towards B. atrox venoms of two antivenoms raised in Costa Rica and Brazil using Bothrops venoms different than B. atrox in the immunization mixtures showed that both antivenoms immunodepleted very efficiently the major toxins (PIII-SVMPs, serine proteinases, CRISP, LAO) of paedomorphic venoms from Puerto Ayacucho (Venezuelan Amazonia) through São Bento, but had impaired reactivity towards PLA(2) and P-I SVMP molecules abundantly present in ontogenetic venoms. The degree of immunodepletion achieved suggests that each of these antivenoms may be effective against envenomations by paedomorphic, and some ontogenetic, B. atrox venoms.

Snake venomics and antivenomics of Bothrops atrox venoms from Colombia and the Amazon regions of Brazil, Perú and Ecuador suggest the occurrence of geographic variation of venom phenotype by a trend towards paedomorphism

The venom proteomes of Bothrops atrox from Colombia, Brazil, Ecuador, and Perú were characterized using venomic and antivenomic strategies. Our results evidence the existence of two geographically differentiated venom phenotypes. The venom from Colombia comprises at least 26 different proteins belonging to 9 different groups of toxins. PI-metalloproteinases and K49-PLA(2) molecules represent the most abundant toxins. On the other hand, the venoms from Brazilian, Ecuadorian, and Peruvian B. atrox contain predominantly PIII-metalloproteinases. These toxin profiles correlate with the venom phenotypes of adult and juvenile B. asper from Costa Rica, respectively, suggesting that paedomorphism represented a selective trend during the trans-Amazonian southward expansion of B. atrox through the Andean Corridor. The high degree of crossreactivity of a Costa Rican polyvalent (Bothrops asper, Lachesis stenophrys, Crotalus simus) antivenom against B. atrox venoms further evidenced the close evolutionary kinship between B. asper and B. atrox. This antivenom was more efficient immunodepleting proteins from the venoms of B. atrox from Brazil, Ecuador, and Perú than from Colombia. Such behaviour may be rationalized taking into account the lower content of poorly immunogenic toxins, such as PLA(2) molecules and PI-SVMPs in the paedomorphic venoms. The immunological profile of the Costa Rican antivenom strongly suggests the possibility of using this antivenom for the management of snakebites by B. atrox in Colombia and the Amazon regions of Ecuador, Perú and Brazil.

Snake venomics of the South and Central American Bushmasters. Comparison of the toxin composition of Lachesis muta gathered from proteomic versus transcriptomic analysis

We report the proteomic characterization of the venoms of two closely related pit vipers of the genus Lachesis, L. muta (South American Bushmaster) and L. stenophrys (Central American Bushmaster), and compare the toxin repertoire of the former revealed through a proteomic versus a transcriptomic approach. The protein composition of the venoms of Lachesis muta and L. stenophrys were analyzed by RP-HPLC, N-terminal sequencing, MALDI-TOF peptide mass fingerprinting and CID-MS/MS. Around 30-40 proteins of molecular masses in the range of 13-110 kDa and belonging, respectively, to only 8 and 7 toxin families were identified in L. muta and L. stenophrys venoms. In addition, both venoms contained a large number of bradykinin-potentiating peptides (BPP) and a C-type natriuretic peptide (C-NP). BPPs and C-NP comprised around 15% of the total venom proteins. In both species, the most abundant proteins were Zn(2+)-metalloproteinases (32-38%) and serine proteinases (25-31%), followed by PLA(2)s (9-12%), galactose-specific C-type lectin (4-8%), l-amino acid oxidase (LAO, 3-5%), CRISP (1.8%; found in L. muta but not in L. stenophrys), and NGF (0.6%). On the other hand, only six L. muta venom-secreted proteins matched any of the previously reported 11 partial or full-length venom gland transcripts, and venom proteome and transcriptome depart in their relative abundances of different toxin families. As expected from their close phylogenetic relationship, the venoms of L. muta and L. stenophrys share (or contain highly similar) proteins, in particular BPPs, serine proteinases, a galactose-specific C-type lectin, and LAO. However, they dramatically depart in their respective PLA(2) complement. Intraspecific quantitative and qualitative differences in the expression of PLA(2) molecules were found when the venoms of five L. muta specimens (3 from Bolivia and 2 from Peru) and the venom of the same species purchased from Sigma were compared. These observations indicate that these class of toxins represents a rapidly-evolving gene family, and suggests that functional differences due to structural changes in PLA(2)s molecules among these snakes may have been a hallmark during speciation and adaptation of diverging snake populations to new ecological niches, or competition for resources in existing ones. Our data may contribute to a deeper understanding of the biology and ecology of these snakes, and may also serve as a starting point for studying structure-function correlations of individual toxins.

Trends in Snakebite Envenomation Therapy: Scientific, Technological and Public Health Considerations

The therapy of snakebite envenomation has been based on the parenteral administration of animal-derived antivenoms. Despite the success of this treatment at reducing the impact of snakebite mortality and morbidity, mostly due to their capacity to neutralize systemically-acting toxins, antivenoms are of relatively low efficacy in the prevention of venom-induced local tissue damage, which often leads to permanent disability. The issue of safety also remains a concern, particularly for some antivenoms which induce a relatively high incidence of adverse reactions. Consequently, there is a need to improve the therapy of snakebite envenomations on the following lines: (a) the technologies to produce antivenoms require improvements aimed at obtaining more refined preparations of higher efficacy and safety, while being affordable for the public health systems of developing countries. (b) The growing knowledge on the biochemistry and toxicology of snake venoms should pave the way for the identification of natural and synthetic inhibitors of venom toxins, particularly of those involved in local tissue pathology. Such inhibitors might become a highly effective therapeutic tool for the abrogation of venom-induced local tissue damage. (c) A better knowledge of the inflammatory events secondary to venom actions may open the possibility of modulating such response, in order to prevent further tissue damage and to promote successful tissue repair and regeneration. A global partnership, involving many participants and combining scientific, technological and public health actions, is required to achieve a leap forward in the treatment of snakebite envenomations world-wide.

Integrated “omics” profiling indicates that miRNAs are modulators of the ontogenetic venom composition shift in the Central American rattlesnake, Crotalus simus simus

BackgroundUnderstanding the processes that drive the evolution of snake venom is a topic of great research interest in molecular and evolutionary toxinology. Recent studies suggest that ontogenetic changes in venom composition are genetically controlled rather than environmentally induced. However, the molecular mechanisms underlying these changes remain elusive. Here we have explored the basis and level of regulation of the ontogenetic shift in the venom composition of the Central American rattlesnake, Crotalus s. simus using a combined proteomics and transcriptomics approach.ResultsProteomic analysis showed that the ontogenetic shift in the venom composition of C. s. simus is essentially characterized by a gradual reduction in the expression of serine proteinases and PLA2 molecules, particularly crotoxin, a β-neurotoxic heterodimeric PLA2, concominantly with an increment of PI and PIII metalloproteinases at age 9–18 months. Comparison of the transcriptional activity of the venom glands of neonate and adult C. s. simus specimens indicated that their transcriptomes exhibit indistinguisable toxin family profiles, suggesting that the elusive mechanism by which shared transcriptomes generate divergent venom phenotypes may operate post-transcriptionally. Specifically, miRNAs with frequency count of 1000 or greater exhibited an uneven distribution between the newborn and adult datasets. Of note, 590 copies of a miRNA targeting crotoxin B-subunit was exclusively found in the transcriptome of the adult snake, whereas 1185 copies of a miRNA complementary to a PIII-SVMP mRNA was uniquely present in the newborn dataset. These results support the view that age-dependent changes in the concentration of miRNA modulating the transition from a crotoxin-rich to a SVMP-rich venom from birth through adulhood can potentially explain what is observed in the proteomic analysis of the ontogenetic changes in the venom composition of C. s. simus.ConclusionsExisting snake venom toxins are the result of early recruitment events in the Toxicofera clade of reptiles by which ordinary genes were duplicated, and the new genes selectively expressed in the venom gland and amplified to multigene families with extensive neofunctionalization throughout the approximately 112–125 million years of ophidian evolution. Our findings support the view that understanding the phenotypic diversity of snake venoms requires a deep knowledge of the mechanisms regulating the transcriptional and translational activity of the venom gland. Our results suggest a functional role for miRNAs. The impact of specific miRNAs in the modulation of venom composition, and the integration of the mechanisms responsible for the generation of these miRNAs in the evolutionary landscape of the snakes venom gland, are further challenges for future research.