High-Voltage Toxin'Roll: Electrostatic Charge Repulsion as a Dynamic Venom Resistance Trait in Pythonid Snakes

The evolutionary interplay between predator and prey has significantly shaped the development of snake venom, a critical adaptation for subduing prey. This arms race has spurred the diversification of the components of venom and the corresponding emergence of resistance mechanisms in the prey and predators of venomous snakes. Our study investigates the molecular basis of venom resistance in pythons, focusing on electrostatic charge repulsion as a defense against α-neurotoxins binding to the alpha-1 subunit of the postsynaptic nicotinic acetylcholine receptor. Through phylogenetic and bioactivity analyses of orthosteric site sequences from various python species, we explore the prevalence and evolution of amino acid substitutions that confer resistance by electrostatic repulsion, which initially evolved in response to predatory pressure by Naja (cobra) species (which occurs across Africa and Asia). The small African species Python regius retains the two resistance-conferring lysines (positions 189 and 191) of the ancestral Python genus, conferring resistance to sympatric Naja venoms. This differed from the giant African species Python sebae, which has secondarily lost one of these lysines, potentially due to its rapid growth out of the prey size range of sympatric Naja species. In contrast, the two Asian species Python brongersmai (small) and Python bivittatus (giant) share an identical orthosteric site, which exhibits the highest degree of resistance, attributed to three lysine residues in the orthosteric sites. One of these lysines (at orthosteric position 195) evolved in the last common ancestor of these two species, which may reflect an adaptive response to increased predation pressures from the sympatric α-neurotoxic snake-eating genus Ophiophagus (King Cobras) in Asia. All these terrestrial Python species, however, were less neurotoxin-susceptible than pythons in other genera which have evolved under different predatory pressure as: the Asian species Malayopython reticulatus which is arboreal as neonates and juveniles before rapidly reaching sizes as terrestrial adults too large for sympatric Ophiophagus species to consider as prey; and the terrestrial Australian species Aspidites melanocephalus which occupies a niche, devoid of selection pressure from α-neurotoxic predatory snakes. Our findings underline the importance of positive selection in the evolution of venom resistance and suggest a complex evolutionary history involving both conserved traits and secondary evolution. This study enhances our understanding of the molecular adaptations that enable pythons to survive in environments laden with venomous threats and offers insights into the ongoing co-evolution between venomous snakes and their prey.

The Cloning and Characterization of a Three-Finger Toxin Homolog (NXH8) from the Coralsnake Micrurus corallinus That Interacts with Skeletal Muscle Nicotinic Acetylcholine Receptors

Coralsnakes (Micrurus spp.) are the only elapids found throughout the Americas. They are recognized for their highly neurotoxic venom, which is comprised of a wide variety of toxins, including the stable, low-mass toxins known as three-finger toxins (3FTx). Due to difficulties in venom extraction and availability, research on coralsnake venoms is still very limited when compared to that of other Elapidae snakes like cobras, kraits, and mambas. In this study, two previously described 3FTx from the venom of M. corallinus, NXH1 (3SOC1_MICCO), and NXH8 (3NO48_MICCO) were characterized. Using in silico, in vitro, and ex vivo experiments, the biological activities of these toxins were predicted and evaluated. The results showed that only NXH8 was capable of binding to skeletal muscle cells and modulating the activity of nAChRs in nerve-diaphragm preparations. These effects were antagonized by anti-rNXH8 or antielapidic sera. Sequence analysis revealed that the NXH1 toxin possesses eight cysteine residues and four disulfide bonds, while the NXH8 toxin has a primary structure similar to that of non-conventional 3FTx, with an additional disulfide bond on the first loop. These findings add more information related to the structural diversity present within the 3FTx class, while expanding our understanding of the mechanisms of the toxicity of this coralsnake venom and opening new perspectives for developing more effective therapeutic interventions.

Distinct regulatory networks control toxin gene expression in elapid and viperid snakes

BACKGROUND

Venom systems are ideal models to study genetic regulatory mechanisms that underpin evolutionary novelty. Snake venom glands are thought to share a common origin, but there are major distinctions between venom toxins from the medically significant snake families Elapidae and Viperidae, and toxin gene regulatory investigations in elapid snakes have been limited. Here, we used high-throughput RNA-sequencing to profile gene expression and microRNAs between active (milked) and resting (unmilked) venom glands in an elapid (Eastern Brown Snake, Pseudonaja textilis), in addition to comparative genomics, to identify cis- and trans-acting regulation of venom production in an elapid in comparison to viperids (Crotalus viridis and C. tigris).

RESULTS

Although there is conservation in high-level mechanistic pathways regulating venom production (unfolded protein response, Notch signaling and cholesterol homeostasis), there are differences in the regulation of histone methylation enzymes, transcription factors, and microRNAs in venom glands from these two snake families. Histone methyltransferases and transcription factor (TF) specificity protein 1 (Sp1) were highly upregulated in the milked elapid venom gland in comparison to the viperids, whereas nuclear factor I (NFI) TFs were upregulated after viperid venom milking. Sp1 and NFI cis-regulatory elements were common to toxin gene promoter regions, but many unique elements were also present between elapid and viperid toxins. The presence of Sp1 binding sites across multiple elapid toxin gene promoter regions that have been experimentally determined to regulate expression, in addition to upregulation of Sp1 after venom milking, suggests this transcription factor is involved in elapid toxin expression. microRNA profiles were distinctive between milked and unmilked venom glands for both snake families, and microRNAs were predicted to target a diversity of toxin transcripts in the elapid P. textilis venom gland, but only snake venom metalloproteinase transcripts in the viperid C. viridis venom gland. These results suggest differences in toxin gene posttranscriptional regulation between the elapid P. textilis and viperid C. viridis.

CONCLUSIONS

Our comparative transcriptomic and genomic analyses between toxin genes and isoforms in elapid and viperid snakes suggests independent toxin regulation between these two snake families, demonstrating multiple different regulatory mechanisms underpin a venomous phenotype.

Genetic insights: mapping sex-specific loci in Siamese cobra (Naja kaouthia) sheds light on the putative sex determining region

The location of female-specific/linked loci identified in Siamese cobra (Naja kaouthia) previously has been determined through in silico chromosome mapping of the Indian cobra genome (N. naja) as a reference genome. In the present study, we used in silico chromosome mapping to identify sex-specific and linked loci in Siamese cobra. Many sex-specific and sex-linked loci were successfully mapped on the Z sex chromosome, with 227 of the 475 specific loci frequently mapped in a region covering 57 Mb and positioned at 38,992,675-95,561,177 bp of the Indian cobra genome (N. naja). This suggested the existence of a putative sex-determining region (SDR), with one specific locus (PA100000600) homologous to the TOPBP1 gene. The involvement of TOPBP1 gene may lead to abnormal synaptonemal complexes and meiotic chromosomal defects, resulting in male infertility. These findings offer valuable insights into the genetic basis and functional aspects of sex-specific traits in the Siamese cobra, which will contribute to our understanding of snake genetics and evolutionary biology.

miRNAs derived from cobra venom exosomes contribute to the cobra envenomation

Currently, there is an increasing amount of evidence indicating that exosomes and the miRNAs they contain are crucial players in various biological processes. However, the role of exosomes and miRNAs in snake venom during the envenomation process remains largely unknown. In this study, fresh venom from Naja atra of different ages (2-month-old, 1-year-old, and 5-year-old) was collected, and exosomes were isolated through ultracentrifugation. The study found that exosomes with inactivated proteins and enzymes can still cause symptoms similar to cobra envenomation, indicating that substances other than proteins and enzymes in exosomes may also play an essential role in cobra envenomation. Furthermore, the expression profiles of isolated exosome miRNAs were analyzed. The study showed that a large number of miRNAs were co-expressed and abundant in cobra venom exosomes (CV-exosomes) of different ages, including miR-2904, which had high expression abundance and specific sequences. The specific miR-2094 derived from CV-exosomes (CV-exo-miR-2904) was overexpressed both in vitro and in vivo. As a result, CV-exo-miR-2904 induced symptoms similar to cobra envenomation in mice and caused liver damage, demonstrating that it plays a crucial role in cobra envenomation. These results reveal that CV-exosomes and the miRNAs they contain play a significant regulatory role in cobra envenomation. Our findings provide new insights for the treatment of cobra bites and the development of snake venom-based medicines.

Domain loss enabled evolution of novel functions in the snake three-finger toxin gene superfamily

Three-finger toxins (3FTXs) are a functionally diverse family of toxins, apparently unique to venoms of caenophidian snakes. Although the ancestral function of 3FTXs is antagonism of nicotinic acetylcholine receptors, redundancy conferred by the accumulation of duplicate genes has facilitated extensive neofunctionalization, such that derived members of the family interact with a range of targets. 3FTXs are members of the LY6/UPAR family, but their non-toxin ancestor remains unknown. Combining traditional phylogenetic approaches, manual synteny analysis, and machine learning techniques (including AlphaFold2 and ProtT5), we have reconstructed a detailed evolutionary history of 3FTXs. We identify their immediate ancestor as a non-secretory LY6, unique to squamate reptiles, and propose that changes in molecular ecology resulting from loss of a membrane-anchoring domain and changes in gene expression, paved the way for the evolution of one of the most important families of snake toxins.

Highly Evolvable: Investigating Interspecific and Intraspecific Venom Variation in Taipans (Oxyuranus spp.) and Brown Snakes (Pseudonaja spp.)

Snake venoms are complex mixtures of toxins that differ on interspecific (between species) and intraspecific (within species) levels. Whether venom variation within a group of closely related species is explained by the presence, absence and/or relative abundances of venom toxins remains largely unknown. Taipans (Oxyuranus spp.) and brown snakes (Pseudonaja spp.) represent medically relevant species of snakes across the Australasian region and provide an excellent model clade for studying interspecific and intraspecific venom variation. Using liquid chromatography with ultraviolet and mass spectrometry detection, we analyzed a total of 31 venoms covering all species of this monophyletic clade, including widespread localities. Our results reveal major interspecific and intraspecific venom variation in Oxyuranus and Pseudonaja species, partially corresponding with their geographical regions and phylogenetic relationships. This extensive venom variability is generated by a combination of the absence/presence and differential abundance of venom toxins. Our study highlights that venom systems can be highly dynamical on the interspecific and intraspecific levels and underscores that the rapid toxin evolvability potentially causes major impacts on neglected tropical snakebites.

Elucidating the Venom Diversity in Sri Lankan Spectacled Cobra (Naja naja) through De Novo Venom Gland Transcriptomics, Venom Proteomics and Toxicity Neutralization

Inadequate effectiveness of Indian antivenoms in treating envenomation caused by the Spectacled Cobra/Indian Cobra (Naja naja) in Sri Lanka has been attributed to geographical variations in the venom composition. This study investigated the de novo venom-gland transcriptomics and venom proteomics of the Sri Lankan N. naja (NN-SL) to elucidate its toxin gene diversity and venom variability. The neutralization efficacy of a commonly used Indian antivenom product in Sri Lanka was examined against the lethality induced by NN-SL venom in mice. The transcriptomic study revealed high expression of 22 toxin genes families in NN-SL, constituting 46.55% of total transcript abundance. Three-finger toxins (3FTX) were the most diversely and abundantly expressed (87.54% of toxin gene expression), consistent with the dominance of 3FTX in the venom proteome (72.19% of total venom proteins). The 3FTX were predominantly S-type cytotoxins/cardiotoxins (CTX) and α-neurotoxins of long-chain or short-chain subtypes (α-NTX). CTX and α-NTX are implicated in local tissue necrosis and fatal neuromuscular paralysis, respectively, in envenomation caused by NN-SL. Intra-species variations in the toxin gene sequences and expression levels were apparent between NN-SL and other geographical specimens of N. naja, suggesting potential antigenic diversity that impacts antivenom effectiveness. This was demonstrated by limited potency (0.74 mg venom/ml antivenom) of the Indian polyvalent antivenom (VPAV) in neutralizing the NN-SL venom. A pan-regional antivenom with improved efficacy to treat N. naja envenomation is needed.

Venom of the Annulated Sea Snake Hydrophis cyanocinctus: A Biochemically Simple but Genetically Complex Weapon

Given that the venom system in sea snakes has a role in enhancing their secondary adaption to the marine environment, it follows that elucidating the diversity and function of venom toxins will help to understand the adaptive radiation of sea snakes. We performed proteomic and de novo NGS analyses to explore the diversity of venom toxins in the annulated sea snake (Hydrophis cyanocinctus) and estimated the adaptive molecular evolution of the toxin-coding unigenes and the toxicity of the major components. We found three-finger toxins (3-FTxs), phospholipase A2 (PLA2) and cysteine-rich secretory protein (CRISP) in the venom proteome and 59 toxin-coding unigenes belonging to 24 protein families in the venom-gland transcriptome; 3-FTx and PLA2 were the most abundant families. Nearly half of the toxin-coding unigenes had undergone positive selection. The short- (i.p. 0.09 μg/g) and long-chain neurotoxin (i.p. 0.14 μg/g) presented fairly high toxicity, whereas both basic and acidic PLA2s expressed low toxicity. The toxicity of H. cyanocinctus venom was largely determined by the 3-FTxs. Our data show the venom is used by H. cyanocinctus as a biochemically simple but genetically complex weapon and venom evolution in H. cyanocinctus is presumably driven by natural selection to deal with fast-moving prey and enemies in the marine environment.

De Novo Venom-Gland Transcriptomics of Spine-Bellied Sea Snake (Hydrophis curtus) from Penang, Malaysia-Next-Generation Sequencing, Functional Annotation and Toxinological Correlation

Envenomation resulted from sea snake bite is a highly lethal health hazard in Southeast Asia. Although commonly caused by sea snakes of Hydrophiinae, each species is evolutionarily distinct and thus, unveiling the toxin gene diversity within individual species is important. Applying next-generation sequencing, this study investigated the venom-gland transcriptome of Hydrophis curtus (spine-bellied sea snake) from Penang, West Malaysia. The transcriptome was de novo assembled, followed by gene annotation and sequence analyses. Transcripts with toxin annotation were only 96 in number but highly expressed, constituting 48.18% of total FPKM in the overall transcriptome. Of the 21 toxin families, three-finger toxins (3FTX) were the most abundantly expressed and functionally diverse, followed by phospholipases A₂. Lh_FTX001 (short neurotoxin) and Lh_FTX013 (long neurotoxin) were the most dominant 3FTXs expressed, consistent with the pathophysiology of envenomation. Lh_FTX001 and Lh_FTX013 were variable in amino acid compositions and predicted epitopes, while Lh_FTX001 showed high sequence similarity with the short neurotoxin from Hydrophis schistosus, supporting cross-neutralization effect of Sea Snake Antivenom. Other toxins of low gene expression, for example, snake venom metalloproteinases and L-amino acid oxidases not commonly studied in sea snake venom were also identified, enriching the knowledgebase of sea snake toxins for future study.

Electric Blue: Molecular Evolution of Three-Finger Toxins in the Long-Glanded Coral Snake Species Calliophis bivirgatus

The genus Calliophis is the most basal branch of the family Elapidae and several species in it have developed highly elongated venom glands. Recent research has shown that C. bivirgatus has evolved a seemingly unique toxin (calliotoxin) that produces spastic paralysis in their prey by acting on the voltage-gated sodium (NaV) channels. We assembled a transcriptome from C. bivirgatus to investigate the molecular characteristics of these toxins and the venom as a whole. We find strong confirmation that this genus produces the classic elapid eight-cysteine three-finger toxins, that δδ-elapitoxins (toxins that resemble calliotoxin) are responsible for a substantial portion of the venom composition, and that these toxins form a distinct clade within a larger, more diverse clade of C. bivirgatus three-finger toxins. This broader clade of C. bivirgatus toxins also contains the previously named maticotoxins and is somewhat closely related to cytotoxins from other elapids. However, the toxins from this clade that have been characterized are not themselves cytotoxic. No other toxins show clear relationships to toxins of known function from other species.

Convergent evolution of pain-inducing defensive venom components in spitting cobras

Convergent evolution provides insights into the selective drivers underlying evolutionary change. Snake venoms, with a direct genetic basis and clearly defined functional phenotype, provide a model system for exploring the repeated evolution of adaptations. While snakes use venom primarily for predation, and venom composition often reflects diet specificity, three lineages of cobras have independently evolved the ability to spit venom at adversaries. Using gene, protein, and functional analyses, we show that the three spitting lineages possess venoms characterized by an up-regulation of phospholipase A2 (PLA2) toxins, which potentiate the action of preexisting venom cytotoxins to activate mammalian sensory neurons and cause enhanced pain. These repeated independent changes provide a fascinating example of convergent evolution across multiple phenotypic levels driven by selection for defense.

Protein composition and biochemical characterization of venom from Sonoran Coral Snakes (Micruroides euryxanthus)

The elapid genus, Micruroides, is considered the sister clade of all New World coral snakes (Genus Micrurus), is monotypic, and is represented by Sonoran Coral Snakes, Micruroides euryxanthus. Coral snakes of the genus Micrurus have been reported to have venoms that are predominantly composed of phospholipases A₂ (PLA₂) or three finger toxins (3FTx), but the venoms of the genus Micruroides are almost completely unstudied. Here, we present the first description of the venom of M. euryxanthus including identification of some proteins as well as transcriptomic, and biological activity assays. The most abundant components within M. euryxanthus venom are 3FTxs (62.3%) and there was relatively low proportion of PLA₂s (14.2%). The venom phenotype supports the hypothesis that the common ancestor of Micrurus and Micruroides had a 3FTx-dominated venom. Within the venom, there were two nearly identical α-neurotoxins (α-Ntx), one of which was designated Eurytoxin, that account for approximately 60% of the venom's lethality to mice. Eurytoxin was cloned, expressed in a soluble and active form, and used to produce rabbit hyperimmune serum. This allowed the analysis of its immunochemical properties, showing them to be different from the recombinant αNTx D.H., present in the venoms of some species of Micrurus. Finally, we observed that the commercial antivenom produced in Mexico for coral snake envenomation is unable to neutralize the lethality from M. euryxanthus venom. This work allowed the classification of Micruroides venom into the 3FTx-predominant group and identified the main components responsible for toxicity to mice.

Horse immunization with short-chain consensus α-neurotoxin generates antibodies against broad spectrum of elapid venomous species

Antivenoms are fundamental in the therapy for snakebites. In elapid venoms, there are toxins, e.g. short-chain α-neurotoxins, which are quite abundant, highly toxic, and consequently play a major role in envenomation processes. The core problem is that such α-neurotoxins are weakly immunogenic, and many current elapid antivenoms show low reactivity towards them. We have previously developed a recombinant consensus short-chain α-neurotoxin (ScNtx) based on sequences from the most lethal elapid venoms from America, Africa, Asia, and Oceania. Here we report that an antivenom generated by immunizing horses with ScNtx can successfully neutralize the lethality of pure recombinant and native short-chain α-neurotoxins, as well as whole neurotoxic elapid venoms from diverse genera such as Micrurus, Dendroaspis, Naja, Walterinnesia, Ophiophagus and Hydrophis. These results provide a proof-of-principle for using recombinant proteins with rationally designed consensus sequences as universal immunogens for developing next-generation antivenoms with higher effectiveness and broader neutralizing capacity.

Proteomic and functional variation within black snake venoms (Elapidae: Pseudechis)

Pseudechis (black snakes) is an Australasian elapid snake genus that inhabits much of mainland Australia, with two representatives confined to Papua New Guinea. The present study is the first to analyse the venom of all 9 described Pseudechis species (plus one undescribed species) to investigate the evolution of venom composition and functional activity. Proteomic results demonstrated that the typical Pseudechis venom profile is dominated by phospholipase A₂ toxins. Strong cytotoxicity was the dominant function for most species. P. porphyriacus, the most basal member of the genus, also exhibited the most divergent venom composition, being the only species with appreciable amounts of procoagulant toxins. The relatively high presence of factor Xa recovered in P. porphyriacus venom may be related to a predominantly amphibian diet. Results of this study provide important insights to guide future ecological and toxinological investigations.

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.

Synthetic peptide antigens derived from long-chain alpha-neurotoxins: Immunogenicity effect against elapid venoms

Three-finger toxins (3FTXs), especially α-neurotoxins, are the most poorly neutralized elapid snake toxins by current antivenoms. In this work, the conserved structural similarity and motif arrangements of long-chain α-neurotoxins led us to design peptides with consensus sequences. Eight long-chain α-neurotoxins (also known as Type II) were used to generate a consensus sequence from which two peptides were chemically synthesized, LCP1 and LCP2. Rabbit sera raised against them were able to generate partially-neutralizing antibodies, which delayed mice mortality in neutralization assays against Naja haje, Dendrospis polylepis and Ophiophagus hannah venoms.

Amplification of lipotoxic cardiomyopathy in the VDR gene knockout mouse

Previous studies demonstrated that the liganded vitamin D receptor (VDR) plays an important role in controlling cardiovascular homeostasis. Both the whole animal VDR gene knockout (VDR^-/-) and the myocyte-specific VDR gene deletion result in changes in cardiac structure and function. Clinical states associated with cardiac steatosis (obesity and diabetes mellitus) are also associated with low circulating 25 OH vitamin D levels. We, therefore, examined the effects of VDR deficiency (VDR^-/- mouse) in a murine model of cardiac steatosis that expresses the terminal enzyme involved in triglyceride synthesis, diacylglycerol acyltransferase 1 (DGAT1), selectively in the cardiac myocyte. These mice display early cardiac dysfunction and late cardiomyopathy and heart failure. In the present study, we demonstrate that mice harboring both genetic modifications (i.e., MHC-DGAT1 Tg and VDR^-/-) exhibit an increase in myocyte size, heart weight/body weight ratio and natriuretic peptide gene expression, all markers of cardiac hypertrophy, that exceed that seen in either VDR^-/- or the MHC-DGAT1 Tg mice alone. This was accompanied by a dramatic increase in interstitial fibrosis and increased expression of collagen 1a1 and collagen 3a1, as well as the osteopontin and matrix metalloproteinase 2, genes. At a functional level, this resulted in a 37% reduction in ejection fraction and 55% reduction in fractional shortening in the DGAT1; VDR^-/- mice relative to the controls. Collectively, these data demonstrate that deficiency in the vitamin D signaling system enhances the pathological phenotype in this experimental cardiomyopathy and suggest an important role for vitamin D in modulating disease severity in common cardiovascular disorders.

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.

Comparison of the primary structures, cytotoxicities, and affinities to phospholipids of five kinds of cytotoxins from the venom of Indian cobra, Naja naja

The molecular mechanism underlying the hemolytic and cytolytic processes of cobra cytotoxins (CTXs) is not yet fully elucidated. To examine this, we analyzed the amino acid sequences, hemolytic and cytotoxic activities, and affinities to phospholipids of the five major CTXs purified from the venom of Indian cobra, Naja naja. CTX2, CTX7, and CTX8 belonged to S-type, and CTX9 and CTX10 to P-type. Comparisons of CTX7 with CTX8 and CTX9 with CTX10 revealed similar primary structures and hemolytic and cytolytic activities. CTX2, whose primary structure was rather different from the others, showed several times weaker hemolytic and cytolytic biological activities than the others. The comparison of CTX2 with CTX7 suggested the importance of Lys30 in loop II for the strong hemolytic and cytolytic activities of S-type CTXs. Cloning of 12 CTX cDNAs from the Naja naja venom cDNA library revealed that 18 out of 23 substitutions found in CTX cDNAs were nonsynonymous. This clearly indicated the accelerated evolution of CTX genes. Multiple sequence alignment of 51 kinds of CTX cDNAs and calculations of nonsynonymous and synonymous substitutions indicated that the codons coding the three loops' regions, which may interact with the hydrophobic tails of phospholipids, have undergone an accelerated evolution. In contrast, the codons coding for amino acid residues considered to participate in the recognition and binding of the hydrophilic head groups of phospholipids, eight Cys residues, and those likely stabilizing β core structure, were all conserved.

The venom-gland transcriptome of the eastern coral snake (Micrurus fulvius) reveals high venom complexity in the intragenomic evolution of venoms

BACKGROUND

Snake venom is shaped by the ecology and evolution of venomous species, and signals of positive selection in toxins have been consistently documented, reflecting the role of venoms as an ecologically critical phenotype. New World coral snakes (Elapidae) are represented by three genera and over 120 species and subspecies that are capable of causing significant human morbidity and mortality, yet coral-snake venom composition is poorly understood in comparison to that of Old World elapids. High-throughput sequencing is capable of identifying thousands of loci, while providing characterizations of expression patterns and the molecular evolutionary forces acting within the venom gland.

RESULTS

We describe the de novo assembly and analysis of the venom-gland transcriptome of the eastern coral snake (Micrurus fulvius). We identified 1,950 nontoxin transcripts and 116 toxin transcripts. These transcripts accounted for 57.1% of the total reads, with toxins accounting for 45.8% of the total reads. Phospholipases A(2) and three-finger toxins dominated expression, accounting for 86.0% of the toxin reads. A total of 15 toxin families were identified, revealing venom complexity previously unknown from New World coral snakes. Toxins exhibited high levels of heterozygosity relative to nontoxins, and overdominance may favor gene duplication leading to the fixation of advantageous alleles. Phospholipase A(2) expression was uniformly distributed throughout the class while three-finger toxin expression was dominated by a handful of transcripts, and phylogenetic analyses indicate that toxin divergence may have occurred following speciation. Positive selection was detected in three of the four most diverse toxin classes, suggesting that venom diversification is driven by recurrent directional selection.

CONCLUSIONS

We describe the most complete characterization of an elapid venom gland to date. Toxin gene duplication may be driven by heterozygote advantage, as the frequency of polymorphic toxin loci was significantly higher than that of nontoxins. Diversification among toxins appeared to follow speciation reflecting species-specific adaptation, and this divergence may be directly related to dietary shifts and is suggestive of a coevolutionary arms race.

Venom proteomic characterization and relative antivenom neutralization of two medically important Pakistani elapid snakes (Bungarus sindanus and Naja naja)

Intra- and interspecific variation in venom composition has been shown to have a major effect upon the efficacy of antivenoms. Due to the absence of domestically produced antivenoms, Pakistan is wholly reliant upon antivenoms produced in other countries, such as India. However, the efficacy of these antivenoms in neutralising the venoms of Pakistani snakes has not been ascertained. This is symptomatic of the general state of toxicological research in this country, which has a myriad of highly toxic and medically important venomous animals. Thus, there is a dire need for knowledge regarding the fundamental proteomics of these venoms and applied knowledge of the relative efficacy of foreign antivenoms. Here we present the results of our proteomic research on two medically important snakes of Pakistan: Bungarus sindanus and Naja naja. Indian Polyvalent Antivenom (Bharat Serums and Vaccines Ltd), which is currently marketed for use in Pakistan, was completely ineffective against either Pakistani species. In addition to the expected pre- and post-synaptic neurotoxic activity, the venom of the Pakistan population of N. naja was shown to be quite divergent from other populations of this species in being potently myotoxic. These results highlight the importance of studying divergent species and isolated populations, where the same data not only elucidates clinical problems in need of immediate attention, but also uncovers sources for novel toxins with potentially useful activities.

BIOLOGICAL SIGNIFICANCE

Pakistan Bungarus sindanus and Naja naja venoms are differentially complex. Naja naja is potently myotoxic. Neither venom is neutralized by Indian antivenom. These results have direct implications for the treatment of envenomed patients in Pakistan. The unusually myotoxic effects of Naja naja demonstrates the value of studying remote populations for biodiscovery.

Trypsin and chymotrypsin inhibitor peptides from the venom of Chinese Daboia russellii siamensis

Two trypsin inhibitors and one chymotrypsin inhibitor from Chinese Daboia russellii siamensis venom, denoted as CBPTI-1, CBPTI-2 and CBPTI-3 were purified, characterized and cloned from lyophilized venom-derived cDNA libraries. The N-terminus of CBPTI-1 was modified and not amenable to Edman degradation sequencing, however an internal partial sequence was found to be SGRCRGHLRRIYYNPDSNKCE. The N-termini of CBPTI-2 and CBPTI-3 were unmodified and their partial sequences were established as HDRPTFCNLAPESGRCRAH and HDRPKFCYLPADPGECMAYIRSFYYDS respectively. From cloning studies CBPTI-1 was found to consist of 66 amino acid residues, while CBPTI-2 and CBPTI-3 precursors consist of 60 amino acid residues, including 6 cysteine residues. Another cDNA sequence (CBPTI-4) was also obtained. Alignment of cDNA sequences showed that CBPTI-3 exhibited similar sequence homology to CBPTI-4 cDNA except for an 8 nucleotide deletion in the open-reading frame. CBPTI-1 and CBPTI-2 were demonstrated to be potent trypsin inhibitors, but were also shown to be effectively potent in chymotrypsin inhibition. The K(i) values of CBPTI-1 and CBPTI-2 for trypsin inhibition were 4.07 × 10(-7) M and 6.66 × 10(-7) M, respectively, and they were non-competitive in their activity. CBPTI-3 showed chymotrypsin inhibition activity with a K(i) value of 2.55 × 10(-9) M, but did not show trypsin inhibitor activity.

Isolation, characterization, cloning and expression of an alpha-neurotoxin from the venom of the Mexican coral snake Micrurus laticollaris (Squamata: Elapidae)

A new member of short chain α-neurotoxic protein family from venom of the Mexican coral snake, Micrurus laticollaris, was characterized. This protein, named MlatA1, possesses 61 amino acids with 8 conserved cysteine residues, sharing 30-91% sequence identity with other fully sequenced Micrurus toxins. MlatA1 (LD50i.v. = 0.064 mg/kg) antagonizes with both fetal and adult nicotinic acetylcholine receptor (nAChR) as well as α-7 neuronal nAChR in a dose-dependent way. Specific rabbit anti-Mlat serum (titer higher than 18,000) does not show any protective ability against this toxin, nevertheless it was able to recognize protein bands in six out of twelve Micrurus venoms showing the existence of two distinct antigenic groups for α-neurotoxins in North American coral snakes species. The MlatA1 gene was cloned and used to produce recombinant toxin (rMlatA1) that was recognized by rabbit anti-native toxin but was depleted of toxic activity.

Conservative mutations in the C2 domains of factor VIII and factor V alter phospholipid binding and cofactor activity

Factor VIII and factor V share structural homology and bind to phospholipid membranes via tandem, lectin-like C domains. Their respective C2 domains bind via 2 pairs of hydrophobic amino acids and an amphipathic cluster. In contrast, the factor V-like, homologous subunit (Pt-FV) of a prothrombin activator from Pseudonaja textilis venom is reported to function without membrane binding. We hypothesized that the distinct membrane-interactive amino acids of these proteins contribute to the differing membrane-dependent properties. We prepared mutants in which the C2 domain hydrophobic amino acid pairs were changed to the homologous residues of the other protein and a factor V mutant with 5 amino acids changed to those from Pt-FV (FV(MTTS/Y)). Factor VIII mutants were active on additional membrane sites and had altered apparent affinities for factor X. Some factor V mutants, including FV(MTTS/Y), had increased membrane interaction and apparent membrane-independent activity that was the result of phospholipid retained during purification. Phospholipid-free FV(MTTS/Y) showed increased activity, particularly a 10-fold increase in activity on membranes lacking phosphatidylserine. The reduced phosphatidylserine requirement correlated to increased activity on resting and stimulated platelets. We hypothesize that altered membrane binding contributes to toxicity of Pt-FV.

Venom gland transcriptomes of two elapid snakes (Bungarus multicinctus and Naja atra) and evolution of toxin genes

BACKGROUND

Kraits (genus Bungarus) and cobras (genus Naja) are two representative toxic genera of elapids in the old world. Although they are closely related genera and both of their venoms are very toxic, the compositions of their venoms are very different. To unveil their detailed venoms and their evolutionary patterns, we constructed venom gland cDNA libraries and genomic bacterial artificial chromosome (BAC) libraries for Bungarus multicinctus and Naja atra, respectively. We sequenced about 1500 cDNA clones for each of the venom cDNA libraries and screened BAC libraries of the two snakes by blot analysis using four kinds of toxin probes; i.e., three-finger toxin (3FTx), phospholipase A2 (PLA2), kunitz-type protease inhibitor (Kunitz), and natriuretic peptide (NP).

RESULTS

In total, 1092 valid expressed sequences tags (ESTs) for B. multicinctus and 1166 ESTs for N. atra were generated. About 70% of these ESTs can be annotated as snake toxin transcripts. 3FTx (64.5%) and β bungarotoxin (25.1%) comprise the main toxin classes in B. multicinctus, while 3FTx (95.8%) is the dominant toxin in N. atra. We also observed several less abundant venom families in B. multicinctus and N. atra, such as PLA2, C-type lectins, and Kunitz. Peculiarly a cluster of NP precursors with tandem NPs was detected in B. multicinctus. A total of 71 positive toxin BAC clones in B. multicinctus and N. atra were identified using four kinds of toxin probes (3FTx, PLA2, Kunitz, and NP), among which 39 3FTx-positive BACs were sequenced to reveal gene structures of 3FTx toxin genes.

CONCLUSIONS

Based on the toxin ESTs and 3FTx gene sequences, the major components of B. multicinctus venom transcriptome are neurotoxins, including long chain alpha neurotoxins (α-ntx) and the recently originated β bungarotoxin, whereas the N. atra venom transcriptome mainly contains 3FTxs with cytotoxicity and neurotoxicity (short chain α-ntx). The data also revealed that tandem duplications contributed the most to the expansion of toxin multigene families. Analysis of nonsynonymous to synonymous nucleotide substitution rate ratios (dN/dS) indicates that not only multigene toxin families but also other less abundant toxins might have been under rapid diversifying evolution.

Characterisation of a mannose-binding C-type lectin from Oxyuranus scutellatus snake venom

C-type lectins are calcium-dependent sugar binding proteins and are distributed ubiquitously amongst vertebrate organisms. As part of a wider study on Australian snake venom components, we have identified and characterised a C-type lectin from the venom of Oxyuranus scutellatus (Australian coastal taipan) with mannose-binding activity. This protein exhibited a subunit molecular mass of 15 kDa and was found to bind mannose and also bind to and agglutinate erythrocytes in a Ca(2+)-dependent manner. cDNA transcripts coding for C-lectin proteins were cloned and sequenced from six Australian elapid snake species and an antibody generated against the O. scutellatus mannose-binding C-lectin identified C-lectin proteins in the venom of 13 Australian elapid snakes by immunoblotting. Experimental evidence and molecular modelling also suggest that this protein exhibits a unique dimeric structure. This is the first confirmed example of a snake venom C-lectin with mannose-binding activity.

Osmoregulation of natriuretic peptide receptors in bromoethylamine-treated rat kidney

Extracellular osmolarity is known as an important factor for the regulation of natriuretic peptide receptors (NPRs). We investigated the intra-renal osmoregulation of NPRs using renal medullectomized rats with bromoethylamine hydrobromide (BEA, 200mg/kg). The administration of BEA caused the decreased food intake and body weight. Water intake was decreased on the first day and then increased from the second day. Urine volume was persistently increased from the first day and free water clearance was also increased from the second day. Urinary excretions of sodium and potassium were decreased on the second day and then recovered to control level. Plasma levels of atrial natriuretic peptide (ANP) and Dendroaspis natriuretic peptide (DNP) in BEA-treated rats were not different from control rats. The inactive renin was increased. The maximum binding capacities of (125)I-ANP as well as (125)I-DNP decreased in glomeruli and medulla of BEA-treated rat kidneys but the binding affinity was not changed. In renal cortex, the gene expressions of ANP, NPR-A, and NPR-B were not changed but that of NPR-C decreased. In renal medulla, the gene expressions of NPR-A, -B, and -C decreased without change in ANP mRNA. Both renal medullary osmolarity and sodium concentration by BEA treatment were lower than those in control kidney. The cGMP concentrations in renal medulla and urine in BEA-treated rats were higher than those in control rats. These results suggest that the increased cGMP production may be partly involved in the decrease in NPRs mRNA expression and their binding capacities by BEA-induced medullectomy.

Transcriptomic basis for an antiserum against Micrurus corallinus (coral snake) venom

BACKGROUND

Micrurus corallinus (coral snake) is a tropical forest snake belonging to the family Elapidae. Its venom shows a high neurotoxicity associated with pre- and post-synaptic toxins, causing diaphragm paralysis, which may result in death. In spite of a relatively small incidence of accidents, serum therapy is crucial for those bitten. However, the adequate production of antiserum is hampered by the difficulty in obtaining sufficient amounts of venom from a small snake with demanding breeding conditions. In order to elucidate the molecular basis of this venom and to uncover possible immunogens for an antiserum, we generated expressed sequences tags (ESTs) from its venom glands and analyzed the transcriptomic profile. In addition, their immunogenicity was tested using DNA immunization.

RESULTS

A total of 1438 ESTs were generated and grouped into 611 clusters. Toxin transcripts represented 46% of the total ESTs. The two main toxin classes consisted of three-finger toxins (3FTx) (24%) and phospholipases A(2) (PLA(2)s) (15%). However, 8 other classes of toxins were present, including C-type lectins, natriuretic peptide precursors and even high-molecular mass components such as metalloproteases and L-amino acid oxidases. Each class included an assortment of isoforms, some showing evidence of alternative splicing and domain deletions. Five antigenic candidates were selected (four 3FTx and one PLA(2)) and used for a preliminary study of DNA immunization. The immunological response showed that the sera from the immunized animals were able to recognize the recombinant antigens.

CONCLUSION

Besides an improvement in our knowledge of the composition of coral snake venoms, which are very poorly known when compared to Old World elapids, the expression profile suggests abundant and diversified components that may be used in future antiserum formulation. As recombinant production of venom antigens frequently fails due to complex disulfide arrangements, DNA immunization may be a viable alternative. In fact, the selected candidates provided an initial evidence of the feasibility of this approach, which is less costly and not dependent on the availability of the venom.

Functional characterization of human C3/cobra venom factor hybrid proteins for therapeutic complement depletion

Cobra venom factor (CVF) is a structural and functional analog of complement C3 isolated from cobra venom. Both CVF and C3b can bind factor B and subsequently form the bimolecular C3/C5 convertases CVF,Bb or C3b,Bb, respectively. The two homologous enzymes exhibit several differences of which the difference in physico-chemical stability is most important, allowing continuous activation of C3 and C5 by CVF,Bb, leading to serum complement depletion. Here we describe the detailed functional properties of two hybrid proteins in which the 113 or 315 C-terminal residues of C3 were replaced with corresponding CVF sequences. Both hybrid proteins formed stable convertases that exhibited C3-cleaving activity, although at different rates. Neither convertase cleaved C5. Both convertases showed partial resistance to inactivation by factors H and I, allowing them to deplete complement in human serum. These data demonstrate that functionally important structural differences between CVF and C3 are located in the very C-terminal region of both homologous proteins, and that small substitutions in human C3 with homologous CVF sequence result in C3 derivatives with CVF-like functions. Such hybrid proteins are important tools to study the structure/function relationships in both C3 and CVF, and these "humanized CVF" proteins may become reagents for therapeutic complement depletion.

Isolation, expression and characterization of a novel dual serine protease inhibitor, OH-TCI, from king cobra venom

Snake venom Kunitz/BPTI members are good tools for understanding of structure-functional relationship between serine proteases and their inhibitors. A novel dual Kunitz/BPTI serine proteinase inhibitor named OH-TCI (trypsin- and chymotrypsin-dual inhibitor from Ophiophagus hannah) was isolated from king cobra venom by three chromatographic steps of gel filtration, trypsin affinity and reverse phase HPLC. OH-TCI is composed of 58 amino acid residues with a molecular mass of 6339Da. Successful expression of OH-TCI was performed as the maltose-binding fusion protein in E. coli DH5alpha. Much different from Oh11-1, the purified native and recombinant OH-TCI both had strong inhibitory activities against trypsin and chymotrypsin although the sequence identity (74.1%) between them is very high. The inhibitor constants (K(i)) of recombinant OH-TCI were 3.91 x 10(-7) and 8.46 x10(-8)M for trypsin and chymotrypsin, respectively. To our knowledge, it was the first report of Kunitz/BPTI serine proteinase inhibitor from snake venom that had equivalent trypsin and chymotrypsin inhibitory activities.

A novel serine protease inhibitor from Bungarus fasciatus venom

By Sephadex G-50 gel filtration, cation-exchange CM-Sephadex C-25 chromatography and reversed phase high-performance liquid chromatography (HPLC), a novel serine protease inhibitor named bungaruskunin was purified and characterized from venom of Bungarus fasciatus. Its cDNA was also cloned from the cDNA library of B. fasciatus venomous glands. The predicted precursor is composed of 83 amino acid (aa) residues including a 24-aa signal peptide and a 59-aa mature bungaruskunin. Bungaruskunin showed maximal similarity (64%) with the predicted serine protease inhibitor blackelin deduced from the cDNA sequence of the red-bellied black snake Pseudechis porphyriacus. Bungaruskunin is a Kunitz protease inhibitor with a conserved Kunitz domain and could exert inhibitory activity against trypsin, chymotrypsin, and elastase. By screening the cDNA library, two new B chains of beta-bungarotoxin are also identified. The overall structures of bungaruskunin and beta-bungarotoxin B chains are similar; especially they have highly conserved signal peptide sequences. These findings strongly suggest that snake Kunitz/BPTI protease inhibitors and neurotoxic homologs may have originated from a common ancestor.

Distinct activities of novel neurotoxins from Australian venomous snakes for nicotinic acetylcholine receptors

Envenomation from Australian elapid snakes results in a myriad of neurological effects due to post-synaptic neurotoxins that bind and inhibit nicotinic acetylcholine receptors (nAChRs) of neurons and muscle fibres. However, despite the significant physiological effects of these toxins, they have remained largely undercharacterised at the molecular level. This study describes the identification and comparative analysis of multiple neurotoxin isoforms from ten Australian snakes, including functional characterisation of two of these isoforms, Os SNTX-1 from Oxyuranus scutellatus and the more potent Pt LNTX-1 from Pseudonaja textilis. Electrophysiological recordings from adrenal chromaffin cells demonstrate that both neurotoxins act as competitive antagonists of nAChRs in a concentration-dependent manner. Their effects upon spontaneous and nerve-evoked membrane responses at the amphibian neuromuscular junction provide further evidence that both toxins bind muscle nAChRs in an irreversible manner. This study represents one of the most comprehensive descriptions to date of the sequences and activity of individual Australian elapid neurotoxins.

Expression pattern of three-finger toxin and phospholipase A2 genes in the venom glands of two sea snakes, Lapemis curtus and Acalyptophis peronii: comparison of evolution of these toxins in land snakes, sea kraits and sea snakes

Background

Snake venom composition varies widely both among closely related species and within the same species, based on ecological variables. In terrestrial snakes, such variation has been proposed to be due to snakes' diet. Land snakes target various prey species including insects (arthropods), lizards (reptiles), frogs and toads (amphibians), birds (aves), and rodents (mammals), whereas sea snakes target a single vertebrate class (fishes) and often specialize on specific types of fish. It is therefore interesting to examine the evolution of toxins in sea snake venoms compared to that of land snakes.

Results

Here we describe the expression of toxin genes in the venom glands of two sea snakes, Lapemis curtus (Spine-bellied Sea Snake) and Acalyptophis peronii (Horned Sea Snake), two members of a large adaptive radiation which occupy very different ecological niches. We constructed cDNA libraries from their venom glands and sequenced 214 and 192 clones, respectively. Our data show that despite their explosive evolutionary radiation, there is very little variability in the three-finger toxin (3FTx) as well as the phospholipase A₂ (PLA₂) enzymes, the two main constituents of Lapemis curtus and Acalyptophis peronii venom. To understand the evolutionary trends among land snakes, sea snakes and sea kraits, pairwise genetic distances (intraspecific and interspecific) of 3FTx and PLA₂ sequences were calculated. Results show that these proteins appear to be highly conserved in sea snakes in contrast to land snakes or sea kraits, despite their extremely divergent and adaptive ecological radiation.

Conclusion

Based on these results, we suggest that streamlining in habitat and diet in sea snakes has possibly kept their toxin genes conserved, suggesting the idea that prey composition and diet breadth may contribute to the diversity and evolution of venom components.

Identification and characterization of Harobin, a novel fibrino(geno)lytic serine protease from a sea snake (Lapemis hardwickii)

A gene encoding a novel serine protease designated as Harobin is cloned and identified from a sea snake venom gland bacteriophage T7 library. It has 265 amino acids and shares 50-70% similarity to terrestrial snake serine proteases. In addition to the 12 conservative Cys, it has three more Cys residues that may contribute to its higher enzymatic stability. Harobin is expressed in Pichia pastoris and purified. Recombinant Harobin exhibits an amidolytic activity, and specifically degrades Aalpha, Bbeta-chain of fibrinogen. It functions as a defibrase both in vitro and in vivo, and reduces thrombosis. Harobin prolongs the coagulation time and the bleeding time of mice and reduces the fibrinogen levels of rats as well. Meanwhile, intravenous injection of Harobin leads to the reduction of blood pressure in SHR rats. It results from the ability of Harobin that cleaves angiotensin I and release bradykinin from plasma kininogen in vitro and in vivo. These data suggest that Harobin is a novel defibrase and has a potential to be an agent for the therapy of thrombosis and hypertension.

Antimicrobial activity of omwaprin, a new member of the waprin family of snake venom proteins

We have isolated and characterized omwaprin, a 50-amino-acid cationic protein from the venom of inland taipan (Oxyuranus microlepidotus). It is a new member of the waprin family of snake venom proteins. A synthetic gene was designed and constructed for expressing the recombinant protein in Escherichia coli. Recombinant omwaprin was used for carrying out functional analyses. The protein is non-toxic to Swiss albino mice at doses of up to 10 mg/kg when administered intraperitoneally. However, it shows selective and dose-dependant antibacterial activity against Gram-positive bacteria. The minimum inhibitory doses were in the range 2-10 microg for selected species of bacteria in radial diffusion assays. The antibacterial activity is salt-tolerant up to 350 mM NaCl. However, omwaprin lost its antibacterial activity upon reduction and alkylation of its cysteine residues, or upon deletion of six N-terminal amino acid residues, four of which are positively charged. These observations indicate that the three-dimensional structure constrained by four disulfide bonds and the N-terminal residues are essential for its activity. The mechanism of action is via membrane disruption, as shown by scanning electron microscopy. Importantly, omwaprin lacks haemolytic activity on human erythrocytes. This demonstrates the specificity of omwaprin for bacterial membranes. Unlike other reported WAP (whey acidic protein) domain-containing antibacterial proteins, including elafin, EPPIN (epididymal proteinase inhibitor), SWAM1 and SWAM2 [single WAP (whey acidic protein) motif proteins 1 and 2] and SLPI (secretory leucocyte proteinase inhibitor), omwaprin shows species-specific activity on the Gram-positive bacteria tested.

Molecular isoforms of cobra venom factor-like proteins in the venom of Austrelaps superbus

Cobra venom factor (CVF) is characteristic of the elapid cobras and has not been reported from venoms of any other families of snakes. During our search for novel proteins, we isolated a polypeptide from the venom of the snake Austrelaps superbus (Lowland Copperhead) that showed structural similarity to C-terminal segment of the alpha-chain of CVF and hence named as AVFalphac (AVF-A. superbus venom factor). cDNA sequence of AVFalphac and its precursor indicated the presence of two isoforms of CVF-like proteins in A. superbus venom gland. This is the first report of molecular isoforms of CVF-like proteins in the venom of an Australian elapid snake. We have determined the complete cDNA sequence of both the isoforms (AVF-1 and AVF-2). They differ in their potential glycosylation sites and the characteristic thioester bond sequence. They display the overall domain structure of CVF and complement C3 proteins. By real-time quantitative analysis, we show that there is a 140-fold difference in the mRNA expression levels of the two isoforms in the venom gland of A. superbus. We also show the presence of AVF-1 and its variant (not AVF-2) in A. superbus venom by partial purification, dot blots, Western blots and peptide mapping using mass spectrometry. Partially purified proteins activate human Factor B in the presence of Factor D and Mg(2+), and deplete the complement activity in human and guinea pig serum. The bimolecular complex (AVFBb) formed activates complement C3 but not complement C5. Thus, AVF proteins may serve as potential candidates for therapeutic complement depletion without side effects. Thus, the discovery of CVF-like proteins in the venom of this Australian elapid snake provides an alternative source of research tools, and contributes to our understanding of the structure-function relationships and evolution of new members of CVF-like proteins.

A model study of interactions between TcAChE peripheral site segment Tyr70Val71 and loop 1 of Fasciculin 2

The continuous chain of residues (Thr7 to Ala12) of Loop1 of Fas2 (F1) and its interaction with the peripheral binding sites (Tyr70-Val71) of AChE (P1) has been studied. Our results suggest that the flexibility of Loop1 might be caused by either the partially protonated guanidine group of Arg11 under experimental conditions or by the interaction with the negatively charged center of substrates. The binding energy of F1-P1 is predicted to be -16.6 kcal/mol at the B3LYP/6-311G(d,p) level, which is assumed to originate from one isolated O7...HN10 H-bond, one possible O10...HC71 unconventional O...HC type H-bonding, and the improved pi-bonding cooperativity around the peptide group of the AChE segment Tyr70-Val71. The classical Kitaura-Morokuma energy decomposition analysis, the NPA charge analysis, and the AIM analysis consistently reveal that the peptide group in segment P1 is more polarizable, which might play the key role in the interactions between F1 and P1. The PCM solvent effect corrected results reveal decrease of the interaction energy of the considered model. The importance of Thr8 of Fas2 in the P-site binding of AChE is also concluded. Site-directed mutations on either the Fas2 residue of Thr8 or the AChE residue of Tyr70 are expected to alter the binding behavior of the Loop1 of Fas2 with AChE.

Post-translational modification accounts for the presence of varied forms of nerve growth factor in Australian elapid snake venoms

The Australian elapid snakes are amongst the most venomous snakes in the world, but much less is known about the overall venom composition in comparison to Asian and American snakes. We have used a combined approach of cDNA cloning and 2-DE with MS to identify nerve growth factor (NGF) in venoms of the Australian elapid snakes and demonstrate its neurite outgrowth activity. While a single 730 nucleotide ORF, coding for a 243 amino acid precursor protein was detected in all snakes, use of 2-DE identified NGF proteins with considerable variation in molecular size within and between the different snakes. The variation in size can be explained at least in part by N-linked glycosylation. It is possible that these modifications alter the stability, activity and other characteristics of the snake NGFs. Further characterisation is necessary to delineate the function of the individual NGF isoforms.

Electrostatic contribution to the binding stability of protein-protein complexes

To investigate roles of electrostatic interactions in protein binding stability, electrostatic calculations were carried out on a set of 64 mutations over six protein-protein complexes. These mutations alter polar interactions across the interface and were selected for putative dominance of electrostatic contributions to the binding stability. Three protocols of implementing the Poisson-Boltzmann model were tested. In vdW4 the dielectric boundary between the protein low dielectric and the solvent high dielectric is defined as the protein van der Waals surface and the protein dielectric constant is set to 4. In SE4 and SE20, the dielectric boundary is defined as the surface of the protein interior inaccessible to a 1.4-A solvent probe, and the protein dielectric constant is set to 4 and 20, respectively. In line with earlier studies on the barnase-barstar complex, the vdW4 results on the large set of mutations showed the closest agreement with experimental data. The agreement between vdW4 and experiment supports the contention of dominant electrostatic contributions for the mutations, but their differences also suggest van der Waals and hydrophobic contributions. The results presented here will serve as a guide for future refinement in electrostatic calculation and inclusion of nonelectrostatic effects.

Novel genes encoding six kinds of three-finger toxins in Ophiophagus hannah (king cobra) and function characterization of two recombinant long-chain neurotoxins

Three-finger toxins are a family of low-molecular-mass toxins (<10 kDa) having very similar three-dimensional structures. In the present study, 19 novel cDNAs coding three-finger toxins were cloned from the venom gland of Ophiophagus hannah (king cobra). Alignment analysis showed that the putative peptides could be divided into six kinds of three-finger toxins: LNTXs (long-chain neurotoxins), short-chain neurotoxins, cardiotoxins (CTXs), weak neurotoxins, muscarinic toxins and a toxin with a free SH group. Furthermore, a phylogenetic tree was established on the basis of the toxin cDNAs and the previously reported similar nucleotide sequences from the same source venom. It indicated that three-finger-toxin genes in O. hannah diverged early in the course of evolution by long- and short-type pathways. Two LNTXs, namely rLNTX1 (recombinant LNTX1) and rLNTX3, were expressed and showed cytolytic activity in addition to their neurotoxic function. By comparing the functional residues, we offer some possible explanations for the differences in their neurotoxic function. Moreover, a plausible elucidation of the additonal cytolytic activity was achieved by hydropathy-profile analysis. This, to our knowledge, is the first observation that recombinant long chain alpha-neurotoxins have a CTX-like cytolytic activity.

Cloning and characterisation of natriuretic peptides from the venom glands of Australian elapids

The venom from Australian elapid snakes contains a complex mixture of polypeptide toxins that adversely affect multiple homeostatic systems within their prey in a highly specific and targeted manner. Included in these toxin families are the recently described venom natriuretic peptides, which display similar structure and vasoactive functions to mammalian natriuretic peptides. This paper describes the identification and detailed comparative analysis of the cDNA transcripts coding for the mature natriuretic peptide from a total of nine Australian elapid snake species. Multiple isoforms were identified in a number of species and represent the first description of a natriuretic peptide from the venom gland for most of these snakes. Two distinct natriuretic peptide isoforms were selected from the common brown snake (Pseudonaja textilis), PtNP-a, and the mulga (Pseudechis australis), PaNP-c, for recombinant protein expression and functional analysis. Only one of these peptides, PtNP-a, displayed cGMP stimulation indicative of normal natriuretic peptide activity. Interestingly, both recombinant peptides demonstrated a dose-dependent inhibition of angiotensin converting enzyme (ACE) activity, which is predictive of the vasoactive effects of the toxin. The natriuretic peptides, however, did not possess any coagulopathic activity, nor did they inhibit or potentiate thrombin, adenosine diphosphate or arachidonic acid induced platelet aggregation. The data presented in this study represent a significant resource for understanding the role of various natriuretic peptides isoforms during the envenomation process by Australian elapid snakes.

Molecular evolution caught in action: gene duplication and evolution of molecular isoforms of prothrombin activators in Pseudonaja textilis (brown snake)

BACKGROUND

The evolution of structurally and functionally similar proteins with highly diverse physiological roles within a single organism is of great interest. Australian elapid snakes offer an excellent opportunity to study the molecular evolution of prothrombin activators. Venom from Pseudonaja textilis contains pseutarin C, a group C prothrombin activator. Its enzymatic subunit is structurally and functionally similar to mammalian factor (F) Xa, whereas its non-enzymatic subunit is similar to FVa. As vertebrates, the snakes also contain a system to activate prothrombin in their own blood during injury. These hemostatic factors are produced in the liver.

RESULTS

Here we describe the presence of two molecular forms of FX expressed in the liver of P. textilis. Both isoforms have molecular signatures and domain architecture of FX. However, one isoform shows approximately 94% sequence identity with the snake FX from Tropidechis carinatus, whereas the other is much closer (90% identity) to the catalytic subunit of pseutarin C (PCCS). Real-time polymerase chain reaction reveals that the latter isoform is expressed approximately 56 000 times lower in the liver of P. textilis. However, the isoforms are not expressed in the venom gland.

CONCLUSION

A detailed analysis of deletions and insertions along with the sequence indicates that the second isoform is an intermediate caught in the evolution of venom prothrombin activator from the blood coagulation FX. Thus, this isoform represents a 'molecular fossil' and reveals the likely evolutionary path of recruitment of FX in the venom gland.

Mutagenesis Studies on the N-terminus and Thr54 of Naja naja atra (Taiwan cobra) Chymotrypsin Inhibitor

Ala-screening mutagenesis studies on Arg1, Pro2, Arg3, Phe4 and Thr54 of Naja naja atra (Taiwan cobra) chymotrypsin inhibitor showed that inhibitory potency and gross conformation of the mutants were not significantly different from those of wild-type inhibitor. Nevertheless, the R1A mutant had an appreciable decrease in the structural stability underlying thermal unfolding and urea-induced denaturation. Alternatively, deleting the first three residues at the N-terminus caused a reduction in structural stability as well as inhibitory potency. In sharp contrast to wild-type and other mutated inhibitors, R1A mutant and truncated mutant completely lost their inhibitory activity when the inhibitors were incubated with chymotrypsin for periods of up to 3 h. The loss of activity correlated with chymotryptic cleavage of inhibitors as evidenced by SDA-PAGE. Taken together, these results reflect that the globally structural rigidity of N. naja atra chymotrypsin inhibitor functionally affects the sustainable period in inhibiting chymotrypsin activity, and that the intact N-terminus might contribute to this event.

Lachesis muta (Viperidae) cDNAs reveal diverging pit viper molecules and scaffolds typical of cobra (Elapidae) venoms: implications for snake toxin repertoire evolution

Efforts to describe toxins from the two major families of venomous snakes (Viperidae and Elapidae) usually reveal proteins belonging to few structural types, particular of each family. Here we carried on an effort to determine uncommon cDNAs that represent possible new toxins from Lachesis muta (Viperidae). In addition to nine classes of typical toxins, atypical molecules never observed in the hundreds of Viperidae snakes studied so far are highly expressed: a diverging C-type lectin that is related to Viperidae toxins but appears to be independently originated; an ohanin-like toxin, which would be the third member of the most recently described class of Elapidae toxins, related to human butyrophilin and B30.2 proteins; and a 3FTx-like toxin, a new member of the widely studied three-finger family of proteins, which includes major Elapidae neurotoxins and CD59 antigen. The presence of these common and uncommon molecules suggests that the repertoire of toxins could be more conserved between families than has been considered, and their features indicate a dynamic process of venom evolution through molecular mechanisms, such as multiple recruitments of important scaffolds and domain exchange between paralogs, always keeping a minimalist nature in most toxin structures in opposition to their nontoxin counterparts.

Bayesian segmental models with multiple sequence alignment profiles for protein secondary structure and contact map prediction

In this paper, we develop a segmental semi-Markov model (SSMM) for protein secondary structure prediction which incorporates multiple sequence alignment profiles with the purpose of improving the predictive performance. The segmental model is a generalization of the hidden Markov model where a hidden state generates segments of various length and secondary structure type. A novel parameterized model is proposed for the likelihood function that explicitly represents multiple sequence alignment profiles to capture the segmental conformation. Numerical results on benchmark data sets show that incorporating the profiles results in substantial improvements and the generalization performance is promising. By incorporating the information from long range interactions in beta-sheets, this model is also capable of carrying out inference on contact maps. This is an important advantage of probabilistic generative models over the traditional discriminative approach to protein secondary structure prediction. The Web server of our algorithm and supplementary materials are available at http://public.kgi.edu/-wild/bsm.html.

Ohanin, a novel protein from king cobra venom: its cDNA and genomic organization

Ohanin, from king cobra venom, is a novel protein which induces hypolocomotion and hyperalgesia in mice [Pung, Y.F., Wong, P.T.H., Kumar, P.P., Hodgson W.C., Kini, R.M., 2005. Ohanin, a novel protein from king cobra venom induces hypolocomotion and hyperalgesia in mice. J. Biol. Chem. 280, 13137-13147.]. It is weakly similar to PRY-SPRY domains (B30.2-like domain). Here we report the complete cDNA and genomic organization of ohanin. Interestingly, cDNA sequence does not show significant sequence similarity to any known sequences, including those of B30.2-like domain-containing proteins. Its full-length cDNA sequence of 1558 bp encodes for prepro-ohanin with a propeptide segment at the C-terminal. Ohanin is the first member of a new subfamily of proteins containing B30.2-like domain with short N-terminal segment. We named this subfamily as vespryns. There are two mRNA subtypes differing in their 5'-untranslated regions. Southern hybridization study shows that ohanin is encoded by a single gene. Its genomic sequence is 7086 bp with five exons and four introns, and the two types of mRNAs are generated by alternative splicing of exon 2. Our results indicate that ohanin and vespryns may have evolved from the same ancestral gene as B30.2 domain.

Identification and analysis of venom gland-specific genes from the coastal taipan (Oxyuranus scutellatus) and related species

Australian terrestrial elapid snakes contain amongst the most potently toxic venoms known. However, despite the well-documented clinical effects of snake bite, little research has focussed on individual venom components at the molecular level. To further characterise the components of Australian elapid venoms, a complementary (cDNA) microarray was produced from the venom gland of the coastal taipan (Oxyuranus scutellatus) and subsequently screened for venom gland-specific transcripts. A number of putative toxin genes were identified, including neurotoxins, phospholipases, a pseudechetoxin-like gene, a venom natriuretic peptide and a nerve growth factor together with other genes involved in cellular maintenance. Venom gland-specific components also included a calglandulin-like protein implicated in the secretion of toxins from the gland into the venom. These toxin transcripts were subsequently identified in seven other related snake species, producing a detailed comparative analysis at the cDNA and protein levels. This study represents the most detailed description to date of the cloning and characterisation of different genes associated with envenomation from Australian snakes.

Cloning and functional expression of venom prothrombin activator protease from Pseudonaja textilis with whole blood procoagulant activity

The snake venom group C prothrombin activators contain a number of components that enhance the rate of prothrombin activation. The cloning and expression of full-length cDNA for one of these components, an activated factor X (factor Xa)-like protease from Pseudonaja textilis as well as the generation of functional chimeric constructs with procoagulant activity were described. The complete cDNA codes for a propeptide, light chain, activation peptide (AP) and heavy chain related in sequence to mammalian factor X. Efficient expression of the protease was achieved with constructs where the AP was deleted and the cleavage sites between the heavy and light chains modified, or where the AP was replaced with a peptide involved in insulin receptor processing. In human kidney cells (H293F) transfected with these constructs, up to 80% of the pro-form was processed to heavy and light chains. Binding of the protease to barium citrate and use of specific antibodies demonstrated that gamma-carboxylation of glutamic acid residues had occurred on the light chain in both cases, as observed in human factor Xa and the native P. textilis protease. The recombinant protease caused efficient coagulation of whole citrated blood and citrated plasma that was enhanced by the presence of Ca2+. This study identified the complete cDNA sequence of a factor Xa-like protease from P. textilis and demonstrated for the first time the expression of a recombinant form of P. textilis protease capable of blood coagulation.