Molecular evolution and structure-function relationships of crotoxin-like and asparagine-6-containing phospholipases A2 in pit viper venoms

Proteomic Analysis of Heloderma horridum horridum Venom: Assessment to Its Transcriptome and Newfound Proteins

Heloderma horridum horridum, a venomous reptile native to America, has a venom with potential applications in treating type II diabetes. In this work, H. h. horridum venom was extracted, lyophilized, and characterized using enzymatic assays for hyaluronidase, phospholipase, and protease. Proteomic analysis of the venom was conducted employing bottom-up/shotgun approaches, SDS-PAGE, high-pH reversed-phase chromatography, and fractionation of tryptic peptides using nano-LC-MS/MS. The proteins found in H. h. horridum venom were reviewed according to the classification of the transcriptome previously reported. The proteomic approach identified 101 enzymes, 36 other proteins, 15 protein inhibitors, 11 host defense proteins, and 1 toxin, including novel venom components such as calcium-binding proteins, phospholipase A2 inhibitors, serpins, cathepsin, subtilases, carboxypeptidase-like, aminopeptidases, glycoside hydrolases, thioredoxin transferases, acid ceramidase-like, enolase, multicopper oxidases, phosphoglucose isomerase (PGI), fructose-1,6-bisphosphatase class 1, pentraxin-related, peptidylglycine α-hydroxylating monooxygenase/peptidyl-hydroxyglycine α-amidating lyase, carbonic anhydrase, acetylcholinesterase, dipeptidylpeptidase, and lysozymes. These findings contribute to understanding the venomous nature of H. h. horridum and highlight its potential as a source of bioactive compounds. Data are available via PRoteomeXchange with the identifier PXD052417.

Variability in antivenom neutralization of Mexican viperid snake venoms

BACKGROUND

Each year, 3,800 cases of snakebite envenomation are reported in Mexico, resulting in 35 fatalities. The only scientifically validated treatment for snakebites in Mexico is the use of antivenoms. Currently, two antivenoms are available in the market, with one in the developmental phase. These antivenoms, produced in horses, consist of F(ab')2 fragments generated using venoms from various species as immunogens. While previous studies primarily focused on neutralizing the venom of the Crotalus species, our study aims to assess the neutralization capacity of different antivenom batches against pit vipers from various genera in Mexico.

METHODOLOGY

We conducted various biological and biochemical tests to characterize the venoms. Additionally, we performed neutralization tests using all three antivenoms to evaluate their effectiveness against lethal activity and their ability to neutralize proteolytic and fibrinogenolytic activities.

RESULTS

Our results reveal significant differences in protein content and neutralizing capacity among different antivenoms and even between different batches of the same product. Notably, the venom of Crotalus atrox is poorly neutralized by all evaluated batches despite being the primary cause of envenomation in the country's northern region. Furthermore, even at the highest tested concentrations, no antivenom could neutralize the lethality of Metlapilcoatlus nummifer and Porthidium yucatanicum venoms. These findings highlight crucial areas for improving existing antivenoms and developing new products.

CONCLUSION

Our research reveals variations in protein content and neutralizing potency among antivenoms, emphasizing the need for consistency in venom characteristics as immunogens. While Birmex neutralizes more LD50 per vial, Antivipmyn excels in specific neutralization. The inability of antivenoms to neutralize certain venoms, especially M. nummifer and P. yucatanicum, highlights crucial improvement opportunities, given the medical significance of these species.

A Review of Rattlesnake Venoms

Venom components are invaluable in biomedical research owing to their specificity and potency. Many of these components exist in two genera of rattlesnakes, Crotalus and Sistrurus, with high toxicity and proteolytic activity variation. This review focuses on venom components within rattlesnakes, and offers a comparison and itemized list of factors dictating venom composition, as well as presenting their known characteristics, activities, and significant applications in biosciences. There are 64 families and subfamilies of proteins present in Crotalus and Sistrurus venom. Snake venom serine proteases (SVSP), snake venom metalloproteases (SVMP), and phospholipases A2 (PLA2) are the standard components in Crotalus and Sistrurus venom. Through this review, we highlight gaps in the knowledge of rattlesnake venom; there needs to be more information on the venom composition of three Crotalus species and one Sistrurus subspecies. We discuss the activity and importance of both major and minor components in biomedical research and drug development.

Comparative Analysis of Alpha-1 Orthosteric-Site Binding by a Clade of Central American Pit Vipers (Genera Atropoides, Cerrophidion, Metlapilcoatlus, and Porthidium)

The distribution and relative potency of post-synaptic neurotoxic activity within Crotalinae venoms has been the subject of less investigation in comparison with Elapidae snake venoms. No previous studies have investigated post-synaptic neurotoxic activity within the Atropoides, Metlapilcoatlus, Cerrophidion, and Porthidium clade. Given the specificity of neurotoxins to relevant prey types, we aimed to uncover any activity present within this clade of snakes that may have been overlooked due to lower potency upon humans and thus not appearing as a clinical feature. Using biolayer interferometry, we assessed the relative binding of crude venoms to amphibian, lizard, bird, rodent and human α-1 nAChR orthosteric sites. We report potent alpha-1 orthosteric site binding in venoms from Atropoides picadoi, Metlapilcoatlus occiduus, M. olmec, M. mexicanus, M. nummifer. Lower levels of binding, but still notable, were evident for Cerrophidion godmani, C. tzotzilorum and C. wilsoni venoms. No activity was observed for Porthidium venoms, which is consistent with significant alpha-1 orthosteric site neurotoxicity being a trait that was amplified in the last common ancestor of Atropoides/Cerrophidion/Metlapilcoatlus subsequent to the split by Porthidium. We also observed potent taxon-selective activity, with strong selection for non-mammalian targets (amphibian, lizard, and bird). As these are poorly studied snakes, much of what is known about them is from clinical reports. The lack of affinity towards mammalian targets may explain the knowledge gap in neurotoxic activity within these species, since symptoms would not appear in bite reports. This study reports novel venom activity, which was previously unreported, indicating toxins that bind to post-synaptic receptors may be more widespread in pit vipers than previously considered. While these effects appear to not be clinically significant due to lineage-specific effects, they are of significant evolutionary novelty and of biodiscovery interest. This work sets the stage for future research directions, such as the use of in vitro and in vivo models to determine whether the alpha-1 orthosteric site binding observed within this study confers neurotoxic venom activity.

Sequence Divergence in Venom Genes Within and Between Montane Pitviper (Viperidae: Crotalinae: Cerrophidion) Species is Driven by Mutation-Drift Equilibrium

Snake venom can vary both among and within species. While some groups of New World pitvipers-such as rattlesnakes-have been well studied, very little is known about the venom of montane pitvipers (Cerrophidion) found across the Mesoamerican highlands. Compared to most well-studied rattlesnakes, which are widely distributed, the isolated montane populations of Cerrophidion may facilitate unique evolutionary trajectories and venom differentiation. Here, we describe the venom gland transcriptomes for populations of C. petlalcalensis, C. tzotzilorum, and C. godmani from Mexico, and a single individual of C. sasai from Costa Rica. We explore gene expression variation in Cerrophidion and sequence evolution of toxins within C. godmani specifically. Cerrophidion venom gland transcriptomes are composed primarily of snake venom metalloproteinases, phospholipase A[Formula: see text]s (PLA[Formula: see text]s), and snake venom serine proteases. Cerrophidion petlalcalensis shows little intraspecific variation; however, C. godmani and C. tzotzilorum differ significantly between geographically isolated populations. Interestingly, intraspecific variation was mostly attributed to expression variation as we did not detect signals of selection within C. godmani toxins. Additionally, we found PLA[Formula: see text]-like myotoxins in all species except C. petlalcalensis, and crotoxin-like PLA[Formula: see text]s in the southern population of C. godmani. Our results demonstrate significant intraspecific venom variation within C. godmani and C. tzotzilorum. The toxins of C. godmani show little evidence of directional selection where variation in toxin sequence is consistent with evolution under a model of mutation-drift equilibrium. Cerrophidion godmani individuals from the southern population may exhibit neurotoxic venom activity given the presence of crotoxin-like PLA[Formula: see text]s; however, further research is required to confirm this hypothesis.

Individual variations in the protein profiles and functional activities of the eyelash palm pit-viper (Bothriechis schlegelii) venom from the Colombian southwest region

Bothriechis schlegelii is a venomous snake found in Central and South America, mainly sighted in regions devoted to agriculture. However, in Colombia, little is known about its contribution to the total envenoming cases. Furthermore, there are no reports of the biochemical and functional activities of venoms from the southwest populations, and the differences respecting other populations are unknown. This study analyzed the protein profiles of venom samples obtained from three specimens originating from this region of Colombia using electrophoresis and chromatography. The lethality, edema-induction, hemorrhagic, defibrinating, coagulant, and indirect hemolytic activities were also evaluated. As a result, venoms were composed of proteins with a wide range of molecular weights, most of them below <37 kDa, with differences between male and female electrophoretic and chromatographic profiles. These variations were also observed in the evaluation of venom functional activities such as pro-coagulant, indirect hemolytic, and edema-inducing activities, whereas neither hemorrhagic nor defibrinating activities were detected. These results are also different considering reports with venom samples from other geographical locations, restating the existence of high intraspecific variability in B. schlegelii venoms, which could have relevant pathophysiological and therapeutic implications.

Dissection of phospholipases A2 reveals multifaceted peptides targeting cancer cells, Leishmania and bacteria

Cationic peptides bio-inspired by natural toxins have been recognized as an efficient strategy for the treatment of different health problems. Due to the specific interaction with substrates from biological membranes, snake venom phospholipases (PLA2s) represent valuable scaffolds for the research and development of short peptides targeting parasites, bacteria, and cancer cells. Considering this, we evaluated the in vitro therapeutic potential of three biomimetic peptides (pCergo, pBmTxJ and pBmje) based on three different amino acid sequences from Asp49 PLA2s. First, short amino acid sequences (12-17 in length) derived from these membranolytic toxins were selected using a combination of bioinformatics tools, including AntiCP, AMPA, PepDraw, ToxinPred, and HemoPI. The peptide, from each polypeptide sequence, with the greatest average antimicrobial index, no toxicity, and no hemolysis predicted was synthesized, purified, and characterized. According to in vitro assays performed, pBmje showed moderate cytotoxicity specifically against MCF-7 (breast cancer cells) with an EC50 of 464.85 µM, whereas pBmTxJ showed an antimicrobial effect against Staphylococcus aureus (ATCC 25923) with an MIC of 37.5 µM, and pCergo against E. coli (ATCC 25922) with an MIC of 75 µM. In addition, pCergo showed antileishmanial activity with an EC50 of 93.69 µM and 110.40 µM against promastigotes of Leishmania braziliensis and L. amazonensis, respectively. Altogether, these results confirmed the versatility of PLA2-derived synthetic peptides, highlighting the relevance of the use of these membrane-interacting toxins as specific archetypes for drug design focused on public health problems.

Omics Technologies for Profiling Toxin Diversity and Evolution in Snake Venom: Impacts on the Discovery of Therapeutic and Diagnostic Agents

Snake venoms are primarily composed of proteins and peptides, and these toxins have developed high selectivity to their biological targets. This makes venoms interesting for exploration into protein evolution and structure-function relationships. A single venom protein superfamily can exhibit a variety of pharmacological effects; these variations in activity originate from differences in functional sites, domains, posttranslational modifications, and the formations of toxin complexes. In this review, we discuss examples of how the major venom protein superfamilies have diversified, as well as how newer technologies in the omics fields, such as genomics, transcriptomics, and proteomics, can be used to characterize both known and unknown toxins.Because toxins are bioactive molecules with a rich diversity of activities, they can be useful as therapeutic and diagnostic agents, and successful examples of toxin applications in these areas are also reviewed. With the current rapid pace of technology, snake venom research and its applications will only continue to expand.

Biochemical and immunochemical characterization of venoms from snakes of the genus Agkistrodon

In the present work, venoms from five species of the genus Agkistrodon were evaluated in terms of their enzymatic (Phospholipase A2 and caseinolytic) and biological (edema forming, hemorrhagic, procoagulant and lethal) effects. Horses were used to produce monovalent hyperimmune sera against each of three venoms (A. bilineatus, A. contortrix and A. piscivorus) and their neutralizing potency, expressed as Median Effective Dose (ED50), was determined against the venoms of all five species. In terms of PLA2 and caseinolytic activities, all venoms are extremely homogeneous. PLA2 activity is high, while caseinolytic activity is low when in contrast with that of the rattlesnake Crotalus simus. On the other hand, biological activities showed marked interspecific differences, particularly between the species from Mexico and those from the United States. Mexican species displayed higher edema-forming, hemorrhagic and lethal effects than US species, while none of the species studied presented procoagulant activity. All three monovalent hyperimmune sera showed good neutralizing potency against the analyzed venoms. Nonetheless, we observed relevant immunochemical differences among the venoms using ELISA and Western Blot assays. We conclude that the venoms of A. piscivorus (USA) and A. bilineatus would be ideal to use as immunogens for the production of a polyvalent antivenom with good neutralizing potency against the venoms of all the species of the genus.

Detection and quantification of a β-neurotoxin (crotoxin homologs) in the venom of the rattlesnakes Crotalus simus, C. culminatus and C. tzabcan from Mexico

Snake venom may vary in composition and toxicity across the geographic distribution of a species. In the case of the three species of the Neotropical rattlesnakes Crotalus simus, C. culminatus and C. tzabcan recent research has revealed that their venoms can contain a neurotoxic component (crotoxin homologs), but is not always the case. In the present work, we detected and quantified crotoxin homologs in venom samples from three species distributed across Mexico, to describe variation at the individual and subspecific levels, using slot blot and ELISA immunoassays. We found that all C. simus individuals analyzed had substantial percentages of crotoxin homologs in their venoms (7.6–44.3%). In contrast, C. culminatus lacked them completely and six of ten individuals of the species C. tzabcan had low percentages (3.0–7.7%). We also found a direct relationship between the lethality of a venom and the percentage of crotoxin homologs it contained, indicating that the quantity of this component influences venom lethality in the rattlesnake C. simus.

•

Monoclonal antibodies were produced that specifically recognized crotoxin homologs in venoms of Crotalus species.

•

Crotoxin homologs were quantified in three species of Crotalus: C. simus, C. culminatus and C. tzabcan.

•

All specimens of C. simus contained crotoxin homologs at different levels, while C. culminatus venoms lacked them completely.

•

In C. tzabcan, some venoms possess and other lack crotoxin homologs.

Intravascular hemolysis induced by phospholipases A2 from the venom of the Eastern coral snake, Micrurus fulvius: Functional profiles of hemolytic and non-hemolytic isoforms

A unique feature of the venom of Micrurus fulvius (Eastern coral snake) is its ability to induce severe intravascular hemolysis in particular species, such as dogs or mice. This effect was previously shown to be induced by distinct phospholipase A₂ (PLA₂) isoforms which cause direct hemolysis in vitro, an uncommon finding for such enzymes. The functional profiles of PLA₂-17, a direct hemolytic enzyme, and PLA₂-12, a co-existing venom isoform lacking such effect, were compared. The enzymes differed not only in their ability to cause intravascular hemolysis: PLA₂-17 additionally displayed lethal, myotoxic, and anticoagulant actions, whereas PLA₂-12 lacked these effects. PLA₂-12 was much more active in hydrolyzing a monodisperse synthetic substrate than PLA₂-17, but the catalytic activity of latter was notably higher on a micellar substrate, or towards pure phospholipid artificial monolayers under controlled lateral pressures. Interestingly, PLA₂-17 could hydrolyze substrate at a pressure of 20 mN m^-1, in contrast to PLA₂-12 or the non-toxic pancreatic PLA₂. This suggests important differences in the monolayer penetrating power, which could be related to differences in toxicity. Comparative examination of primary structures and predicted three-dimensional folding of PLA₂-12 and PLA₂-17, revealed that differences concentrate in their N-terminal and central regions, leading to variations of the surface properties at the membrane interacting interface. PLA₂-17 presents a less basic interfacial surface than PLA₂-12, but more bulky aromatic residues, which could be associated to its higher membrane-penetrating strength. Altogether, these structural and functional comparative observations suggest that the ability of PLA₂s to penetrate substrate interfaces could be a major determinant of toxicity, perhaps more important than protein surface charge.

Tracing the evolution of venom phospholipases A2 in Gloydius strauchii and related pitvipers: A tale of two acidic isozymes

Two acidic Asp⁴⁹-PLA₂s with Glu⁶ substitution and a neutral Lys⁴⁹-PLA (designated Gst-K49) were cloned from G. strauchii venom glands, their full amino acid sequences were deduced. The predominant acidic PLA₂ (designated Gst-E6a) contains 124 residues and the M¹⁸W³⁰ substitutions, while the minor acidic PLA₂ (designated Gst-E6b) contains 122 residues and the V¹⁸A³⁰ substitutions. Their sequences are most similar to those of the respective orthologous PLA₂s of G. intermedius venom. Gst-E6a and Gst-E6b appear to be paralogs and possibly have different predatory targets or functions. The LC-MS/MS results indicate the presence of only three PLA₂ gene products in the crude venom, the relative expression levels were in the order of Gst-E6a ≫ Gst-E6b > Gst-K49, as confirmed by qPCR results. In contrast to other Gloydius, G. strauchii venom does not contain neurotoxic or basic anticoagulant Asp⁴⁹-PLA₂s, but Gst-K49 is the first Lys⁴⁹-PLA₂ identified in Gloydius venoms. However, its venom content is relatively low and its pI value 7.3 is much lower than those of other Lys⁴⁹-PLA₂s and. The Lys⁴⁹-PLA₂ genes appear to regress in the venom of most of Gloydius and related rattlesnake, and this evolutionary regression occurred before the dispersal of Asian pitvipers to the New World.

Coralsnake Venomics: Analyses of Venom Gland Transcriptomes and Proteomes of Six Brazilian Taxa

Venom gland transcriptomes and proteomes of six Micrurus taxa (M. corallinus, M. lemniscatus carvalhoi, M. lemniscatus lemniscatus, M. paraensis, M. spixii spixii, and M. surinamensis) were investigated, providing the most comprehensive, quantitative data on Micrurus venom composition to date, and more than tripling the number of Micrurus venom protein sequences previously available. The six venomes differ dramatically. All are dominated by 2-6 toxin classes that account for 91-99% of the toxin transcripts. The M. s. spixii venome is compositionally the simplest. In it, three-finger toxins (3FTxs) and phospholipases A₂ (PLA₂s) comprise >99% of the toxin transcripts, which include only four additional toxin families at levels ≥0.1%. Micrurus l. lemniscatus venom is the most complex, with at least 17 toxin families. However, in each venome, multiple structural subclasses of 3FTXs and PLA₂s are present. These almost certainly differ in pharmacology as well. All venoms also contain phospholipase B and vascular endothelial growth factors. Minor components (0.1-2.0%) are found in all venoms except that of M. s. spixii. Other toxin families are present in all six venoms at trace levels (<0.005%). Minor and trace venom components differ in each venom. Numerous novel toxin chemistries include 3FTxs with previously unknown 8- and 10-cysteine arrangements, resulting in new 3D structures and target specificities. 9-cysteine toxins raise the possibility of covalent, homodimeric 3FTxs or heterodimeric toxins with unknown pharmacologies. Probable muscarinic sequences may be reptile-specific homologs that promote hypotension via vascular mAChRs. The first complete sequences are presented for 3FTxs putatively responsible for liberating glutamate from rat brain synaptosomes. Micrurus C-type lectin-like proteins may have 6-9 cysteine residues and may be monomers, or homo- or heterodimers of unknown pharmacology. Novel KSPIs, 3× longer than any seen previously, appear to have arisen in three species by gene duplication and fusion. Four species have transcripts homologous to the nociceptive toxin, (MitTx) α-subunit, but all six species had homologs to the β-subunit. The first non-neurotoxic, non-catalytic elapid phospholipase A₂s are reported. All are probably myonecrotic. Phylogenetic analysis indicates that the six taxa diverged 15-35 million years ago and that they split from their last common ancestor with Old World elapines nearly 55 million years ago. Given their early diversification, many cryptic micrurine taxa are anticipated.

Full-Length Venom Protein cDNA Sequences from Venom-Derived mRNA: Exploring Compositional Variation and Adaptive Multigene Evolution

Envenomation of humans by snakes is a complex and continuously evolving medical emergency, and treatment is made that much more difficult by the diverse biochemical composition of many venoms. Venomous snakes and their venoms also provide models for the study of molecular evolutionary processes leading to adaptation and genotype-phenotype relationships. To compare venom complexity and protein sequences, venom gland transcriptomes are assembled, which usually requires the sacrifice of snakes for tissue. However, toxin transcripts are also present in venoms, offering the possibility of obtaining cDNA sequences directly from venom. This study provides evidence that unknown full-length venom protein transcripts can be obtained from the venoms of multiple species from all major venomous snake families. These unknown venom protein cDNAs are obtained by the use of primers designed from conserved signal peptide sequences within each venom protein superfamily. This technique was used to assemble a partial venom gland transcriptome for the Middle American Rattlesnake (Crotalus simus tzabcan) by amplifying sequences for phospholipases A₂, serine proteases, C-lectins, and metalloproteinases from within venom. Phospholipase A₂ sequences were also recovered from the venoms of several rattlesnakes and an elapid snake (Pseudechis porphyriacus), and three-finger toxin sequences were recovered from multiple rear-fanged snake species, demonstrating that the three major clades of advanced snakes (Elapidae, Viperidae, Colubridae) have stable mRNA present in their venoms. These cDNA sequences from venom were then used to explore potential activities derived from protein sequence similarities and evolutionary histories within these large multigene superfamilies. Venom-derived sequences can also be used to aid in characterizing venoms that lack proteomic profiles and identify sequence characteristics indicating specific envenomation profiles. This approach, requiring only venom, provides access to cDNA sequences in the absence of living specimens, even from commercial venom sources, to evaluate important regional differences in venom composition and to study snake venom protein evolution.

This work demonstrates that full-length venom protein messenger RNAs are present in secreted venoms and can be used to acquire full-length protein sequences of toxins from both front-fanged (Elapidae, Viperidae) and rear-fanged (Colubridae) snake venoms, eliminating the need to use venom glands. Full-length transcripts were obtained from venom samples that were fresh, newly lyophilized, old, field desiccated or commercially prepared, representing a significant advance over previous attempts which produced only partial sequence transcripts. Transcripts for all major venom protein families (metalloproteinases, serine proteases, C-type lectins, phospholipases A₂ and three-finger toxins) responsible for clinically significant snakebite symptoms were obtained from venoms. These sequences aid in the identification and characterization of venom proteome profiles, allowing for the identification of peptide sequences, specific isoforms, and novel venom proteins. The application of this technique will help to provide venom protein sequences for many snake species, including understudied rear-fanged snakes. Venom protein transcripts offer important insights into potential snakebite envenomation profiles and the molecular evolution of venom protein multigene families. By requiring only venom to obtain venom protein cDNAs, the approach detailed here will provide access to cDNA-based protein sequences from commercial and other venom sources, facilitating study of snake venom protein composition and evolution.

Complete clinical course of envenomation by Protobothrops mangshanensis: delayed coagulopathy and response to Trimeresurus albolabris antivenom

INTRODUCTION

Protobothrops mangshanensis, the Mangshan pit viper, is a rare pit viper native to the area surrounding Mount Mang in China's Hunan province. Toxicity from envenomation is not well characterized.

CASE DETAILS

A 33-year-old male presented to an emergency department (ED) after being bitten on the forearm by his P. mangshanensis. He complained of mild swelling and pain at the bite site. He was admitted for observation and toxicology consultation. Following initially normal coagulation studies including platelets, prothrombin time (PT), activated partial thromboplastin time (aPTT), fibrinogen and D-dimer, fibrinogen decreased to 121 mg/dL and D-dimer concurrently rose to 377 ng/mL over 24 h. On hospital day 2 fibrinogen stabilized at 109 mg/dL and he was discharged with outpatient laboratory monitoring. Three days later, he returned with bruising to the contralateral arm. Fibrinogen was undetectable (<40 mg/dL) and PT was 14.6 s. He declined admission but returned 2 d later with bruising to the nose. Bloodwork revealed immeasurably prolonged PT, aPTT, and thrombin time, but he eloped. Late that evening he returned and was treated with three vials of Green pit viper (Trimeresurus albolabris) antivenom. Within 24 h coagulopathy improved markedly; at five days, coagulation abnormalities resolved.

DISCUSSION

Mangshan pit viper envenomations may cause isolated hemotoxicity, despite molecular studies suggesting additional neurotoxicity and myotoxicity. T. albolabris antivenom appears effective in treating the resultant coagulopathy.

CONCLUSION

We report the natural history of envenomation by the Mangshan pit viper. A delayed coagulopathy, apparently fibrinolytic in nature, is unaccompanied by local tissue destruction and responsive to Green pit viper antivenom.

Molecular cloning and pharmacological properties of an acidic PLA2 from Bothrops pauloensis snake venom

In this work, we describe the molecular cloning and pharmacological properties of an acidic phospholipase A(2) (PLA(2)) isolated from Bothrops pauloensis snake venom. This enzyme, denominated BpPLA(2)-TXI, was purified by four chromatographic steps and represents 2.4% of the total snake venom protein content. BpPLA(2)-TXI is a monomeric protein with a molecular mass of 13.6 kDa, as demonstrated by Matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF) analysis and its theoretical isoelectric point was 4.98. BpPLA(2)-TXI was catalytically active and showed some pharmacological effects such as inhibition of platelet aggregation induced by collagen or ADP and also induced edema and myotoxicity. BpPLA(2)-TXI displayed low cytotoxicity on TG-180 (CCRF S 180 II) and Ovarian Carcinoma (OVCAR-3), whereas no cytotoxicity was found in regard to MEF (Mouse Embryonic Fibroblast) and Sarcoma 180 (TIB-66). The N-terminal sequence of forty-eight amino acid residues was determined by Edman degradation. In addition, the complete primary structure of 122 amino acids was deduced by cDNA from the total RNA of the venom gland using specific primers, and it was significantly similar to other acidic D49 PLA(2)s. The phylogenetic analyses showed that BpPLA(2)-TXI forms a group with other acidic D49 PLA(2)s from the gender Bothrops, which are characterized by a catalytic activity associated with anti-platelet effects.

Alignment-free distance measure based on return time distribution for sequence analysis: applications to clustering, molecular phylogeny and subtyping

The data deluge in post-genomic era demands development of novel data mining tools. Existing molecular phylogeny analyses (MPAs) developed for individual gene/protein sequences are alignment-based. However, the size of genomic data and uncertainties associated with alignments, necessitate development of alignment-free methods for MPA. Derivation of distances between sequences is an important step in both, alignment-dependant and alignment-free methods. Various alignment-free distance measures based on oligo-nucleotide frequencies, information content, compression techniques, etc. have been proposed. However, these distance measures do not account for relative order of components viz. nucleotides or amino acids. A new distance measure, based on the concept of 'return time distribution' (RTD) of k-mers is proposed, which accounts for the sequence composition and their relative orders. Statistical parameters of RTDs are used to derive a distance function. The resultant distance matrix is used for clustering and phylogeny using Neighbor-joining. Its performance for MPA and subtyping was evaluated using simulated data generated by block-bootstrap, receiver operating characteristics and leave-one-out cross validation methods. The proposed method was successfully applied for MPA of family Flaviviridae and subtyping of Dengue viruses. It is observed that method retains resolution for classification and subtyping of viruses at varying levels of sequence similarity and taxonomic hierarchy.

Venom variability and envenoming severity outcomes of the Crotalus scutulatus scutulatus (Mojave rattlesnake) from Southern Arizona

Twenty-one Mojave rattlesnakes, Crotalus scutulatus scutulatus (C. s. scutulatus), were collected from Arizona and New Mexico U.S.A. Venom proteome of each specimen was analyzed using reverse-phase HPLC and SDS-PAGE. The toxicity of venoms was analyzed using lethal dose 50 (LD(50)). Health severity outcomes between two Arizona counties U.S.A., Pima and Cochise, were determined by retrospective chart review of the Arizona Poison and Drug Information Center (APDIC) database between the years of 2002 and 2009. Six phenotypes (A-F) were identified based on three venom protein families; Mojave toxin, snake venom metalloproteinases PI and PIII (SVMP), and myotoxin-A. Venom changed geographically from SVMP-rich to Mojave toxin-rich phenotypes as you move from south central to southeastern Arizona. Phenotypes containing myotoxin-A were only found in the transitional zone between the SVMP and Mojave toxin phenotypes. Venom samples containing the largest amounts of SVMP or Mojave toxin had the highest and lowest LD(50s), respectively. There was a significant difference when comparing the presence of neurotoxic effects between Pima and Cochise counties (p=0.001). No significant difference was found when comparing severity (p=0.32), number of antivenom vials administered (p=0.17), days spent in a health care facility (p=0.23) or envenomation per 100,000 population (p=0.06). Although not part of the original data to be collected, death and intubations, were also noted. There is a 10× increased risk of death and a 50× increased risk of intubations if envenomated in Cochise County.

Snake venomics of Crotalus tigris: the minimalist toxin arsenal of the deadliest Nearctic rattlesnake venom. Evolutionary Clues for generating a pan-specific antivenom against crotalid type II venoms [corrected]

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.

Enzymatic toxins from snake venom: structural characterization and mechanism of catalysis

Snake venoms are cocktails of enzymes and non-enzymatic proteins used for both the immobilization and digestion of prey. The most common snake venom enzymes include acetylcholinesterases, l-amino acid oxidases, serine proteinases, metalloproteinases and phospholipases A(2) . Higher catalytic efficiency, thermal stability and resistance to proteolysis make these enzymes attractive models for biochemists, enzymologists and structural biologists. Here, we review the structures of these enzymes and describe their structure-based mechanisms of catalysis and inhibition. Some of the enzymes exist as protein complexes in the venom. Thus we also discuss the functional role of non-enzymatic subunits and the pharmacological effects of such protein complexes. The structures of inhibitor-enzyme complexes provide ideal platforms for the design of potent inhibitors which are useful in the development of prototypes and lead compounds with potential therapeutic applications.

Unique structural characteristics and evolution of a cluster of venom phospholipase A2 isozyme genes of Protobothrops flavoviridis snake

Protobothrops flavoviridis (Crotalinae) venom gland phospholipase A(2) (PLA(2)) isozyme genes have evolved in an accelerated manner to acquire diverse physiological activities in their products. For elucidation of the multiplication mechanism of PLA(2) genes, a 25,026 bp genome segment harboring five PLA(2) isozyme genes was obtained from Amami-Oshima P. flavoviridis liver and sequenced. The gene PfPLA 2 encoded [Lys(49)]PLA(2) called BPII, the gene PfPLA 4 neurotoxic [Asp(49)]PLA(2) called PLA-N, the gene PfPLA 5 basic [Asp(49)]PLA(2) called PLA-B, and PfPLA 1(psi) and PfPLA 3(psi) were the inactivated genes. The 5' truncated reverse transcriptase (RT) elements, whose intact forms constitute long interspersed nuclear elements (LINEs), were found in close proximity to the 3' end of PLA(2) genes and named PLA(2) gene-coupled RT fragments (PcRTFs). The facts that PcRTFs have the stem-loop and repetitive sequence in the 3' untranslated region (UTR) which is characteristic of CR1 LINEs suggest that PcRTFs are the debris of P. flavoviridis ancestral CR1 LINEs, denoted as PfCR1s. Since the associated pairs of PLA(2) genes and PcRTFs are arranged in tandem in the 25,026 bp segment, it is thought that an ancestral PLA(2) gene-PfCR1 unit (PfPLA-PfCR1) which was produced by retrotransposition of PfCR1 by itself to the 3' end of PLA(2) gene duplicated several times to form a multimer of PfPLA-PfCR1, a cluster of PLA(2) genes, in the period after Crotalinae and Viperinae snakes branched off. Recombinational hot spot of a 37bp segment, named Scomb, was found in the region 548 bp upstream from the TATA box of PLA(2) genes. Thus, it could be assumed that multiplication of PfPLA-PfCR1 occurred by unequal crossing over of the segment, -Scomb-PfPLA-PfCR1-Scomb-. The PfCR1 moieties were afterward disrupted in the 5' portion to PcRTFs. The detection of two types of PcRTFs different in length which were produced by elimination of two definitive sequences in PfCR1 moiety possibly by gene conversion clearly supports such process but not multiplication of the PLA(2) gene-PcRTF unit.

Absence of phospholipase A(2) in most Crotalus horridus venom due to translation blockage: comparison with Crotalus horridus atricaudatus venom

To investigate the peculiar absence of phospholipases A(2) (PLA(2)s) in most Crotalus horridus (CH) venom, we cloned and sequenced the venom PLA(2)s of three CH specimens from different regions. The results revealed that all the venom glands contained mRNAs that encoded an acidic PLA(2) (designated as either CH-E6 or CH-E6'). The predicted CH-E6 from the Iowan CH and CH-E6' from the South Carolinian CH differed by only one amino acid residue, while the PLA(2) cDNA cloned from the Kentuckian CH contained an early stop codon instead of a Tyr(22) codon. Only the individual South Carolinian CH venom was found to contain the CH-E6' protein whose mass was confirmed by MALDI-TOF analysis. Our results suggest that low PLA(2) expression levels in most CH venom can be attributed to translation blockage. We also purified two acidic PLA(2)s and canebrake toxin from the pooled venom of Crotalus horridus atricaudatus (neurotoxic CH subspecies). One of the acidic PLA(2)s was identical to CH-E6 and showed high lipolytic activity and weak anti-platelet activities. The possibility that C. h. atricaudatus could be a hybrid between CH and Crotalus scutulatus is discussed.

Protein complexes in snake venom

Snake venom contains mixture of bioactive proteins and polypeptides. Most of these proteins and polypeptides exist as monomers, but some of them form complexes in the venom. These complexes exhibit much higher levels of pharmacological activity compared to individual components and play an important role in pathophysiological effects during envenomation. They are formed through covalent and/or non-covalent interactions. The subunits of the complexes are either identical (homodimers) or dissimilar (heterodimers; in some cases subunits belong to different families of proteins). The formation of complexes, at times, eliminates the non-specific binding and enhances the binding to the target molecule. On several occasions, it also leads to recognition of new targets as protein-protein interaction in complexes exposes the critical amino acid residues buried in the monomers. Here, we describe the structure and function of various protein complexes of snake venoms and their role in snake venom toxicity.

Crotalus durissus collilineatus venom gland transcriptome: analysis of gene expression profile

Crotalus durissus rattlesnakes are responsible for the most lethal cases of snakebites in Brazil. Crotalus durissus collilineatus subspecies is related to a great number of accidents in Southeast and Central West regions, but few studies on its venom composition have been carried out to date. In an attempt to describe the transcriptional profile of the C. durissus collilineatus venom gland, we generated a cDNA library and the sequences obtained could be identified by similarity searches on existing databases. Out of 673 expressed sequence tags (ESTs) 489 produced readable sequences comprising 201 singletons and 47 clusters of two or more ESTs. One hundred and fifty reads (60.5%) produced significant hits to known sequences. The results showed a predominance of toxin-coding ESTs instead of transcripts coding for proteins involved in all cellular functions. The most frequent toxin was crotoxin, comprising 88% of toxin-coding sequences. Crotoxin B, a basic phospholipase A(2) (PLA(2)) subunit of crotoxin, was represented in more variable forms comparing to the non-enzymatic subunit (crotoxin A), and most sequences coding this molecule were identified as CB1 isoform from Crotalus durissus terrificus venom. Four percent of toxin-related sequences in this study were identified as growth factors, comprising five sequences for vascular endothelial growth factor (VEGF) and one for nerve growth factor (NGF) that showed 100% of identity with C. durissus terrificus NGF. We also identified two clusters for metalloprotease from PII class comprising 3% of the toxins, and two for serine proteases, including gyroxin (2.5%). The remaining 2.5% of toxin-coding ESTs represent singletons identified as homologue sequences to cardiotoxin, convulxin, angiotensin-converting enzyme inhibitor and C-type natriuretic peptide, Ohanin, crotamin and PLA(2) inhibitor. These results allowed the identification of the most common classes of toxins in C. durissus collilineatus snake venom, also showing some unknown classes for this subspecies and even for C. durissus species, such as cardiotoxins and VEGF.

Cloning, characterization and phylogenetic analyses of members of three major venom families from a single specimen of Walterinnesia aegyptia

Walterinnesia aegyptia is a monotypic elapid snake inhabiting in Africa and Mideast. Although its envenoming is known to cause rapid deaths and paralysis, structural data of its venom proteins are rather limited. Using gel filtration and reverse-phase HPLC, phospholipases A(2) (PLAs), three-fingered toxins (3FTxs), and Kunitz-type protease inhibitors (KIns) were purified from the venom of a single specimen of this species caught in northern Egypt. In addition, specific primers were designed and PCR was carried out to amplify the cDNAs encoding members of the three venom families, respectively, using total cDNA prepared from its venom glands. Complete amino acid sequences of two acidic PLAs, three short chain 3FTxs, and four KIns of this venom species were thus deduced after their cDNAs were cloned and sequenced. They are all novel sequences and match the mass data of purified proteins. For members of each toxin family, protein sequences were aligned and subjected to molecular phylogenetic analyses. The results indicated that the PLAs and a Kunitz inhibitor of W. aegyptia are most similar to those of king cobra venom, and its 3FTxs belongs to either Type I alpha-neurotoxins or weak toxins of orphan-II subtype. It is remarkable that both king cobra and W. aegyptia cause rapid deaths of the victims, and a close evolutionary relationship between them is speculated.

Venom phospholipases of Russell's vipers from Myanmar and eastern India--cloning, characterization and phylogeographic analysis

Venoms of Russell's vipers (genus Daboia) are known for their deadly coagulopathic and other effects. We herein studied various isoforms of venom phospholipases A(2) (PLAs) from two Daboia species at their geographic boundary. From Myanmar Daboia siamensis venom (designated as DsM), four PLAs (designated DsM-aI, aI', aII' and bI') were purified, and the cDNAs encoding two acidic (DsM-aI and aII) and two basic PLAs (DsM-bI and S1) were also cloned from its venom-glands. DsM-S1 is identical to the major venom PLA of southern India Daboia russelii, but the protein is absent from the venom. Additionally, four PLAs (designated DrK-aI, aII, bI and bII) were cloned from cDNA obtained from venom glands of a Kolkata D. russelii, and the PLAs were purified from the pooled venom (designated as DrK). The acidic DrK-aI is the most neurotoxic and lethal among these PLAs; DsM-aI which differs from DrK-aI by only the Phe2 substitution shows greatly reduced enzymatic activity and lethality. Both acidic PLAs do not form dimeric complex with basic PLAs in the same venoms. DsM-bI' is neurotoxic and lethal but its orthologous DrK-bI (97% identical to DsM-bI') is a much weaker toxin. Given the fact that most of the orthologous PLAs of DrK and DsM share 97-100% sequence identity, Daboia vipers of Myanmar and Kolkata must be closely related. Molecular phylogenetic analyses on 30 venom PLAs of Eurasian vipers' revealed co-evolution of five subtypes of venom PLAs in both Daboia and Vipera genera. Our results shed light on the intra- and inter-species variations and structure-function relationships of viperid venom PLAs.

Sequence analysis and phylogenetic relationship of genes encoding heterodimeric phospholipases A2 from the venom of the scorpion Anuroctonus phaiodactylus

Some scorpion venom contain heterodimeric phospholipases A2. They were shown to be toxic to insects and to cause edema and/or hemolysis of mammalian erythrocytes. This manuscript describes the results of cDNA cloning of five different heterodimeric phospholipases from the venomous glands of the Mexican scorpion Anuroctonus phaiodactylus. The amino acid sequence deduced from the heterodimeric phospholipases open reading frames corresponds in each case to a different isoform. The nucleotide sequences corresponding to two of these genes were also obtained by directly sequencing genomic DNA. The cDNA isoforms show high similarity with the heterodimeric phospholipase Phaiodactylipin purified from the same scorpion. However, similar phospholipases were also found in scorpions from other species and the sequences available were used to construct a phylogenetic tree. In order to understand better the gene structure and phylogeny of these enzymes we analyzed their sequences and compared them with secretory phospholipases of other sources from groups I, II and III. The genomic DNA sequence of a similar phospholipase from bee venomous glands was also cloned. The information available on a Drosophila phospholipase was included in this analysis. The phospholipases of groups I and II contain a conserved exon-intron structure (four or five exons of the mature segment of the enzyme are separated by three or four introns). Also, the gene structure of the phospholipases from A. phaiodactylus and that of the bee venom, belonging to group III phospholipases, are interrupted by three introns. The mature peptide of the bee enzyme is a single polypeptide chain, coded by four exons, whereas those from the scorpion studied here although having four exons, showed the presence of two different polypeptides in its native state. The mature protein is processed after synthesis, producing the heterodimeric structure: a long and a short-peptide chain, linked by a disulfide bridge. The small subunit is the one coded by the fourth exon. The human phospholipase A2 and that of Drosophila, also classified into the group III phospholipases, have a considerably different exon-intron organization.

Purification, characterization and cytokine release function of a novel Arg-49 phospholipase A(2) from the venom of Protobothrops mucrosquamatus

Group IIA phospholipase A(2) (PLA(2)) are major components in Viperidae/Crotalidae venom. In the present study, a novel PLA(2) named promutoxin with Arg at the site 49 has been purified from the venom of Protobothrops mucrosquamatus by chromatography. It consists of 122 amino acid residues with a molecular mass of 13,656 Da assessed by MALDI-TOF. It has the structural features of snake venom group IIA PLA(2)s, but has no PLA(2) enzymatic activity. Promutoxin shows higher amino acid sequence identity to the K49 PLA(2)s (72-95%) than to D49 PLA(2)s (52-58%). Promutoxin exhibits potent myotoxicity in the animal model with as little as 1 microg of promutoxin causing myonecrosis and myoedema in the gastrocnemius muscle of mice. Promutoxin is also able to stimulate the release of IL-12, TNFalpha, IL-6 and IL-1beta from human monocytes, and induce IL-2, TNFalpha and IL-6 release from T cells, indicating that this snake venom group IIA PLA(2) is actively involved in the inflammatory process in man caused by snake venom poisoning.

Comparative proteomics and subtyping of venom phospholipases A2 and disintegrins of Protobothrops pit vipers

To explore the venom diversity and systematics of pit vipers under the genus Protobothrops, the venom phospholipases A2 (PLA2s) of P. mangshanensis, P. elegans and P. tokarensis were purified and characterized for the first time. The results were compared with the corresponding venom data of other co-generic species including P. mucrosquamatus, P. flavoviridis and P. jerdonii. Based on sequence features at the N-terminal regions, we identified five PLA2 subtypes, i.e., the Asp49-PLA2s with N6, E6 or R6 substitution and the Lys49-PLA2. However, not all subtypes were expressed in each of the species. Venom N6-PLA2s from P. mangshanensis and P. tokarensis venom were weakly neurotoxic toward chick biventer cervicis tissue preparations. The venoms of P. tokarensis and P. flavoviridis contained identical PLA2 isoforms. In most Protobothrop disintegrins, sequences flanking the RGD-motif are conserved. Phylogenetic analyses based on amino acid sequences of both families of the acidic PLA2s and the disintegrins clarify that these species could belong to a monophyletic group.