Alicia Pérez

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.

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.

Snake venomics and antivenomics of Crotalus durissus subspecies from Brazil: Assessment of geographic variation and its implication on snakebite management

We report the comparative proteomic and antivenomic characterization of the venoms of subspecies cascavella and collilineatus of the Brazilian tropical rattlesnake Crotalus durissus. The venom proteomes of C. d. collilineatus and C. d. cascavella comprise proteins in the range of 4-115 kDa belonging to 9 and 8 toxin families, respectively. Collilineatus and cascavella venoms contain 20-25 main toxins belonging to the following protein families: disintegrin, PLA(2), serine proteinase, cysteine-rich secretory protein (CRISP), vascular endothelial growth factor-like (VEGF), L-amino acid oxidase, C-type lectin-like, and snake venom metalloproteinase (SVMP). As judged by reverse-phase HPLC and mass spectrometry, cascavella and collilineatus share about 90% of their venom proteome. However, the relative occurrence of the toxin families departs among the two C. durissus subspecies venoms. The most notable difference is the presence of the myotoxin crotamine in some C. d. collilineatus specimens (averaging 20.8% of the total proteins of pooled venom), which is absent in the venom of C. d. cascavella. On the other hand, the neurotoxic PLA(2) crotoxin represents the most abundant protein in both C. durissus venoms, comprising 67.4% of the toxin proteome in C. d. collilineatus and 72.5% in C. d. cascavella. Myotoxic PLA(2)s are also present in the two venoms albeit in different relative concentrations (18.1% in C. d. cascavella vs. 4.6% in C. d. collilineatus). The venom composition accounts for the clinical manifestations caused by C. durissus envenomations: systemic neurotoxicity and myalgic symptoms and coagulation disturbances, frequently accompanied by myoglobinuria and acute renal failure. The overall compositions of C. d. subspecies cascavella and collilineatus venoms closely resemble that of C. d. terrificus, supporting the view that these taxa can be considered geographical variations of the same species. Pooled venom from adult C.d. cascavella and neonate C.d. terrificus lack crotamine, whereas this skeletal muscle cell membrane depolarizing inducing myotoxin accounts for approximately 20% of the total toxins of venom pooled from C.d. collilineatus and C.d. terrificus from Southern Brazil. The possible relevance of the observed venom variability among the tropical rattlesnake subspecies was assessed by antivenomics using anti-crotalic antivenoms produced at Instituto Butantan and Instituto Vital Brazil. The results revealed that both antivenoms exhibit impaired immunoreactivity towards crotamine and display restricted ( approximately 60%) recognition of PLA(2) molecules (crotoxin and D49-myotoxins) from C. d. cascavella and C. d. terrificus venoms. This poor reactivity of the antivenoms may be due to a combination of factors: on the one hand, an inappropriate choice of the mixture of venoms for immunization and, on the other hand, the documented low immunogenicity of PLA(2) molecules. C. durissus causes most of the lethal snakebite accidents in Brazil. The implication of the geographic variation of venom composition for the treatment of bites by different C. durissus subspecies populations is discussed.

Snake Venomics of Crotalus tigris: The Minimalist Toxin Arsenal of the Deadliest Neartic Rattlesnake Venom. Evolutionary Clues for Generating a Pan-Specific Antivenom against Crotalid Type II Venoms

We report the proteomic and antivenomic characterization of Crotalus tigris venom. This venom exhibits the highest lethality for mice among rattlesnakes and the simplest toxin proteome reported to date. The venom proteome of C. tigris comprises 7-8 gene products from 6 toxin families; the presynaptic β-neurotoxic heterodimeric PLA(2), Mojave toxin, and two serine proteinases comprise, respectively, 66 and 27% of the C. tigris toxin arsenal, whereas a VEGF-like protein, a CRISP molecule, a medium-sized disintegrin, and 1-2 PIII-SVMPs each represent 0.1-5% of the total venom proteome. This toxin profile really explains the systemic neuro- and myotoxic effects observed in envenomated animals. In addition, we found that venom lethality of C. tigris and other North American rattlesnake type II venoms correlates with the concentration of Mojave toxin A-subunit, supporting the view that the neurotoxic venom phenotype of crotalid type II venoms may be described as a single-allele adaptation. Our data suggest that the evolutionary trend toward neurotoxicity, which has been also reported for the South American rattlesnakes, may have resulted by pedomorphism. The ability of an experimental antivenom to effectively immunodeplete proteins from the type II venoms of C. tigris, Crotalus horridus , Crotalus oreganus helleri, Crotalus scutulatus scutulatus, and Sistrurus catenatus catenatus indicated the feasibility of generating a pan-American anti-Crotalus type II antivenom, suggested by the identification of shared evolutionary trends among South and North American Crotalus species.

Molecular cloning of disintegrins from Cerastes vipera and Macrovipera lebetina transmediterranea venom gland cDNA libraries: insight into the evolution of the snake venom integrin-inhibition system

We report the cloning and sequence analysis of Cerastes vipera and Macrovipera lebetina transmediterranea cDNAs coding for short non-RGD (Arg-Gly-Asp) disintegrins and for dimeric disintegrin subunits. The mRNAs belong to the short-coding class, suggesting that these disintegrin mRNAs may be more widely distributed than previously thought. Our data also argue for a common ancestry of the mRNAs of short disintegrins and those coding for precursors of dimeric disintegrin chains. The Macrovipera lebetina transmediterranea dimeric disintegrin reported to inhibit the laminin-binding integrins alpha3beta1, alpha6beta1 and alpha7beta1 was analysed using a proteomic approach and was shown to bear MLD (Met-Leu-Asp) and VGD (Val-Gly-Asp) motifs. The results highlight the fact that disintegrins have evolved a restricted panel of integrin-blocking sequences that segregate with defined branches of the phylogenetic tree of the integrin alpha-chains, providing novel insights into the evolutionary adaptation of the snake venom antagonists to the ligand-binding sites of their target integrin receptors.

cDNA Cloning and Functional Expression of Jerdostatin, a Novel RTS-disintegrin from Trimeresurus jerdonii and a Specific Antagonist of the α1β1 Integrin

Jerdostatin represents a novel RTS-containing short disintegrin cloned by reverse transcriptase-PCR from the venom gland mRNA of the Chinese Jerdons pit viper Trimeresurus jerdonii. The jerdostatins precursor cDNA contained a 333-bp open reading frame encoding a signal peptide, a pre-peptide, and a 43-amino acid disintegrin domain, whose amino acid sequence displayed 80% identity with that of the KTS-disintegrins obtustatin and viperistatin. The jerdostatin cDNA structure represents the first complete open reading frame of a short disintegrin and points to the emergence of jerdostatin from a short-coding gene. The different residues between jerdostatin and obtustatin/viperistatin are segregated within the integrin-recognition loop and the C-terminal tail. Native jerdostatin (r-jerdostatin-R21) and a R21K mutant (r-jerdostatin-K21) were produced in Escherichia coli. In each case, two conformers were isolated. One-dimensional 1H NMR showed that conformers 1 and 2 of r-jerdostatin-R21 represent, respectively, well folded and unfolded proteins. The two conformers of the wild-type and the R21K mutant inhibited the adhesion of α1-K562 cells to collagen IV with IC50 values of 180 and 703 nm, respectively. The IC50 values of conformers 2 of r-jerdostatin-R21 and r-jerdostatin-K21 were, respectively, 5.95 and 12.5 μm. Neither r-jerdostatin-R21 nor r-jerdostatin-K21 showed inhibitory activity toward other integrins, including αIIbβ3, αvβ3, α2β1, α5β1, α4β1, α6β1, and α9β1 up to a concentration of 24 μm. Although the RTS motif appears to be more potent than KTS inhibiting the α1β1 integrin, r-jerdostatin-R21 is less active than the KTS-disintegrins, strongly suggesting that substitutions outside the integrin-binding motif and/or C-terminal proteolytic processing are responsible for the decreased inhibitory activity.

Crystal structure and statistical coupling analysis of highly glycosylated peroxidase from royal palm tree (Roystonea regia)

Royal palm tree peroxidase (RPTP) is a very stable enzyme in regards to acidity, temperature, H(2)O(2), and organic solvents. Thus, RPTP is a promising candidate for developing H(2)O(2)-sensitive biosensors for diverse applications in industry and analytical chemistry. RPTP belongs to the family of class III secretory plant peroxidases, which include horseradish peroxidase isozyme C, soybean and peanut peroxidases. Here we report the X-ray structure of native RPTP isolated from royal palm tree (Roystonea regia) refined to a resolution of 1.85A. RPTP has the same overall folding pattern of the plant peroxidase superfamily, and it contains one heme group and two calcium-binding sites in similar locations. The three-dimensional structure of RPTP was solved for a hydroperoxide complex state, and it revealed a bound 2-(N-morpholino) ethanesulfonic acid molecule (MES) positioned at a putative substrate-binding secondary site. Nine N-glycosylation sites are clearly defined in the RPTP electron-density maps, revealing for the first time conformations of the glycan chains of this highly glycosylated enzyme. Furthermore, statistical coupling analysis (SCA) of the plant peroxidase superfamily was performed. This sequence-based method identified a set of evolutionarily conserved sites that mapped to regions surrounding the heme prosthetic group. The SCA matrix also predicted a set of energetically coupled residues that are involved in the maintenance of the structural folding of plant peroxidases. The combination of crystallographic data and SCA analysis provides information about the key structural elements that could contribute to explaining the unique stability of RPTP.

Insights Into the Structural Basis of the Ph- Dependent Dimer-Tetramer Equilibrium Through Crystallographic Analysis of Recombinant Diocleinae Lectins.

The structural ground underlying the pH-dependency of the dimer-tetramer transition of Diocleinae lectins was investigated by equilibrium sedimentation and X-ray crystal structure determination of wild-type and site-directed mutants of recombinant lectins. Synthetic genes coding for the full-length alpha-chains of the seed lectins of Dioclea guianensis (termed r-alphaDguia) and Dioclea grandiflora (termed r-alphaDGL) were designed and expressed in Escherichia coli. This pioneering approach, which will be described in detail in the present paper, yielded recombinant lectins displaying carbohydrate-binding activity, dimer-tetramer equilibria and crystal structures indistinguishable from their natural homologues. Conversion of the pH-stable tetrameric r-alphaDGL into a structure exhibiting pH-dependent dimer-tetramer transition was accomplished through mutations that abolished the interdimeric interactions at the central cavity of the tetrameric lectins. Both the central and the peripheral interacting regions bear structural information for formation of the canonical legume lectin tetramer. We hypothesize that the strength of the ionic contacts at these sites may be modulated by the pH, leading to dissociation of those lectin structures that are not locked into a pH-stable tetramer through interdimeric contacts networking the central cavity loops.

Characterization of pea (Pisum sativum) seed protein fractions

BACKGROUND Legume seed proteins have to be chemically characterized in order to properly link their nutritional effects with their chemical structure. RESULTS Vicilin and albumin fractions devoid of cross-contamination, as assessed by mass peptide fingerprinting analysis, were obtained from defatted pea (Pisum sativum cv. Bilbo) meal. The extracted protein fractions contained 56.7-67.7 g non-starch polysaccharides kg⁻¹. The vicilin fraction was higher than legumins in arginine, isoleucine, leucine, phenylalanine and lysine. The most abundant amino acids in the albumin fraction were aspartic acid, glutamic acid, lysine and arginine, and the amounts of methionine were more than double than those in legumins and vicilins. The pea albumin fraction showed a clear enrichment of protease inhibitory activity when compared with the seed meal. In vitro digestibility values for pea proteins were 0.63 ±  0.04, 0.88 ±  0.04 and 0.41 ±  0.23 for legumins, vicilins and albumins respectively. CONCLUSION Vicilin and albumin fractions devoid of cross-contamination with other proteins were obtained from pea seed meal. The vicilin fraction also contained low amounts of soluble non-starch polysaccharides and was enriched in isoleucine, leucine, phenylalanine and lysine. In vitro digestibility values for pea proteins were similar or even numerically higher than those for control proteins.

Venomic Analysis of the Poorly Studied Desert Coral Snake, Micrurus tschudii tschudii, Supports the 3FTx/PLA2 Dichotomy across Micrurus Venoms

The venom proteome of the poorly studied desert coral snake Micrurus tschudii tschudii was unveiled using a venomic approach, which identified ≥38 proteins belonging to only four snake venom protein families. The three-finger toxins (3FTxs) constitute, both in number of isoforms (~30) and total abundance (93.6% of the venom proteome), the major protein family of the desert coral snake venom. Phospholipases A2 (PLA2s; seven isoforms, 4.1% of the venom proteome), 1–3 Kunitz-type proteins (1.6%), and 1–2 l-amino acid oxidases (LAO, 0.7%) complete the toxin arsenal of M. t. tschudii. Our results add to the growing evidence that the occurrence of two divergent venom phenotypes, i.e., 3FTx- and PLA2-predominant venom proteomes, may constitute a general trend across the cladogenesis of Micrurus. The occurrence of a similar pattern of venom phenotypic variability among true sea snake (Hydrophiinae) venoms suggests that the 3FTx/PLA2 dichotomy may be widely distributed among Elapidae venoms.