Molecular cloning, characterization and evolution of the gene encoding a new group of short-chain alpha-neurotoxins in an Australian elapid, Pseudonaja textilis

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.

Proteomic comparisons of venoms of long-term captive and recently wild-caught Eastern brown snakes (Pseudonaja textilis) indicate venom does not change due to captivity

Snake venom is a highly variable phenotypic character, and its variation and rapid evolution are important because of human health implications. Because much snake antivenom is produced from captive animals, understanding the effects of captivity on venom composition is important. Here, we have evaluated toxin profiles from six long-term (LT) captive and six recently wild-caught (RC) eastern brown snakes, Pseudonaja textilis, utilizing gel electrophoresis, HPLC-MS, and shotgun proteomics. We identified proteins belonging to the three-finger toxins, group C prothrombin activators, Kunitz-type serine protease inhibitors, and phospholipases A2, among others. Although crude venom HPLC analysis showed LT snakes to be higher in some small molecular weight toxins, presence/absence patterns showed no correlation with time in captivity. Shotgun proteomics indicated the presence of similar toxin families among individuals but with variation in protein species. Although no venom sample contained all the phospholipase A2 subunits that form the textilotoxin, all did contain both prothrombin activator subunits. This study indicates that captivity has limited effects on venom composition, that venom variation is high, and that venom composition may be correlated to geographic distribution.

BIOLOGICAL SIGNIFICANCE

Through proteomic comparisons, we show that protein variation within LT and RC groups of snakes (Pseudonaja textilis) is high, thereby resulting in no discernible differences in venom composition between groups. We utilize complementary techniques to characterize the venom proteomes of 12 individual snakes from our study area, and indicate that individuals captured close to one another have more similar venom gel electrophoresis patterns than those captured at more distant locations. These data are important for understanding natural variation in and potential effects of captivity on venom composition.

Venom down under: dynamic evolution of Australian elapid snake toxins

Despite the unparalleled diversity of venomous snakes in Australia, research has concentrated on a handful of medically significant species and even of these very few toxins have been fully sequenced. In this study, venom gland transcriptomes were sequenced from eleven species of small Australian elapid snakes, from eleven genera, spanning a broad phylogenetic range. The particularly large number of sequences obtained for three-finger toxin (3FTx) peptides allowed for robust reconstructions of their dynamic molecular evolutionary histories. We demonstrated that each species preferentially favoured different types of α-neurotoxic 3FTx, probably as a result of differing feeding ecologies. The three forms of α-neurotoxin [Type I (also known as (aka): short-chain), Type II (aka: long-chain) and Type III] not only adopted differential rates of evolution, but have also conserved a diversity of residues, presumably to potentiate prey-specific toxicity. Despite these differences, the different α-neurotoxin types were shown to accumulate mutations in similar regions of the protein, largely in the loops and structurally unimportant regions, highlighting the significant role of focal mutagenesis. We theorize that this phenomenon not only affects toxin potency or specificity, but also generates necessary variation for preventing/delaying prey animals from acquiring venom-resistance. This study also recovered the first full-length sequences for multimeric phospholipase A₂ (PLA₂) ‘taipoxin/paradoxin’ subunits from non-Oxyuranus species, confirming the early recruitment of this extremely potent neurotoxin complex to the venom arsenal of Australian elapid snakes. We also recovered the first natriuretic peptides from an elapid that lack the derived C-terminal tail and resemble the plesiotypic form (ancestral character state) found in viper venoms. This provides supporting evidence for a single early recruitment of natriuretic peptides into snake venoms. Novel forms of kunitz and waprin peptides were recovered, including dual domain kunitz-kunitz precursors and the first kunitz-waprin hybrid precursors from elapid snakes. The novel sequences recovered in this study reveal that the huge diversity of unstudied venomous Australian snakes are of considerable interest not only for the investigation of venom and whole organism evolution but also represent an untapped bioresource in the search for novel compounds for use in drug design and development.

Three-fingered RAVERs: Rapid Accumulation of Variations in Exposed Residues of snake venom toxins

Three-finger toxins (3FTx) represent one of the most abundantly secreted and potently toxic components of colubrid (Colubridae), elapid (Elapidae) and psammophid (Psammophiinae subfamily of the Lamprophidae) snake venom arsenal. Despite their conserved structural similarity, they perform a diversity of biological functions. Although they are theorised to undergo adaptive evolution, the underlying diversification mechanisms remain elusive. Here, we report the molecular evolution of different 3FTx functional forms and show that positively selected point mutations have driven the rapid evolution and diversification of 3FTx. These diversification events not only correlate with the evolution of advanced venom delivery systems (VDS) in Caenophidia, but in particular the explosive diversification of the clade subsequent to the evolution of a high pressure, hollow-fanged VDS in elapids, highlighting the significant role of these toxins in the evolution of advanced snakes. We show that Type I, II and III α-neurotoxins have evolved with extreme rapidity under the influence of positive selection. We also show that novel Oxyuranus/Pseudonaja Type II forms lacking the apotypic loop-2 stabilising cysteine doublet characteristic of Type II forms are not phylogenetically basal in relation to other Type IIs as previously thought, but are the result of secondary loss of these apotypic cysteines on at least three separate occasions. Not all 3FTxs have evolved rapidly: κ-neurotoxins, which form non-covalently associated heterodimers, have experienced a relatively weaker influence of diversifying selection; while cytotoxic 3FTx, with their functional sites, dispersed over 40% of the molecular surface, have been extremely constrained by negative selection. We show that the a previous theory of 3FTx molecular evolution (termed ASSET) is evolutionarily implausible and cannot account for the considerable variation observed in very short segments of 3FTx. Instead, we propose a theory of Rapid Accumulation of Variations in Exposed Residues (RAVER) to illustrate the significance of point mutations, guided by focal mutagenesis and positive selection in the evolution and diversification of 3FTx.

Transcriptomic analysis of the venom gland of the red-headed krait (Bungarus flaviceps) using expressed sequence tags

Background

The Red-headed krait (Bungarus flaviceps, Squamata: Serpentes: Elapidae) is a medically important venomous snake that inhabits South-East Asia. Although the venoms of most species of the snake genus Bungarus have been well characterized, a detailed compositional analysis of B. flaviceps is currently lacking.

Results

Here, we have sequenced 845 expressed sequence tags (ESTs) from the venom gland of a B. flaviceps. Of the transcripts, 74.8% were putative toxins; 20.6% were cellular; and 4.6% were unknown. The main venom protein families identified were three-finger toxins (3FTxs), Kunitz-type serine protease inhibitors (including chain B of β-bungarotoxin), phospholipase A₂ (including chain A of β-bungarotoxin), natriuretic peptide (NP), CRISPs, and C-type lectin.

Conclusion

The 3FTxs were found to be the major component of the venom (39%). We found eight groups of unique 3FTxs and most of them were different from the well-characterized 3FTxs. We found three groups of Kunitz-type serine protease inhibitors (SPIs); one group was comparable to the classical SPIs and the other two groups to chain B of β-bungarotoxins (with or without the extra cysteine) based on sequence identity. The latter group may be functional equivalents of dendrotoxins in Bungarus venoms. The natriuretic peptide (NP) found is the first NP for any Asian elapid, and distantly related to Australian elapid NPs. Our study identifies several unique toxins in B. flaviceps venom, which may help in understanding the evolution of venom toxins and the pathophysiological symptoms induced after envenomation.

Molecular diversity in venom proteins of the Russell's viper (Daboia russellii russellii) and the Indian cobra (Naja naja) in Sri Lanka

To examine the molecular diversity of the venom proteins of the Russell's viper (Daboia russellii russellii) and the Indian cobra (Naja naja) in Sri Lanka, we isolated 38 venom proteins through a combination of anion exchange chromatography followed by reversed-phase high performance liquid chromatography. From the venom of D. r. russellii we isolated 15 proteins: 5 isozymes of phospholipase A(2) (PLA(2)), 4 serine proteases, 2 C-type lectin-like proteins, 2 L-amino acid oxidases, 1 cysteine-rich secretory protein (CRISP), and 1 metalloproteinase. From the venom of N. naja we isolated 23 proteins: 10 isoforms of cytotoxins (CTX), 7 PLA(2) isozymes, 2 muscarinic toxinlike proteins, 2 CRISPs, 1 nerve growth factor, and 1 new thrombin-like serine protease. Most of these proteins contained new amino acid sequences for each species, indicating molecular diversity in venom proteins. The entire amino acid sequences of PLA(2)3 from D. r. russellii and CTX7 from N. naja were determined. Additionally, the polymorphic amino acid residues of PLA(2)3 were preferentially localized on the potential antigenic sites. While 2 types of PLA(2) (N and S types) were found in D. r. russellii (India) and D. r. siamensis (Java), all the PLA(2)s from D. r. siamensis (Burma) were N type, and those from D. r. russellii (Sri Lanka) were primarily S type.

Role of accelerated segment switch in exons to alter targeting (ASSET) in the molecular evolution of snake venom proteins

Background

Snake venom toxins evolve more rapidly than other proteins through accelerated changes in the protein coding regions. Previously we have shown that accelerated segment switch in exons to alter targeting (ASSET) might play an important role in its functional evolution of viperid three-finger toxins. In this phenomenon, short sequences in exons are radically changed to unrelated sequences and hence affect the folding and functional properties of the toxins.

Results

Here we analyzed other snake venom protein families to elucidate the role of ASSET in their functional evolution. ASSET appears to be involved in the functional evolution of three-finger toxins to a greater extent than in several other venom protein families. ASSET leads to replacement of some of the critical amino acid residues that affect the biological function in three-finger toxins as well as change the conformation of the loop that is involved in binding to specific target sites.

Conclusion

ASSET could lead to novel functions in snake venom proteins. Among snake venom serine proteases, ASSET contributes to changes in three surface segments. One of these segments near the substrate binding region is known to affect substrate specificity, and its exchange may have significant implications for differences in isoform catalytic activity on specific target protein substrates. ASSET therefore plays an important role in functional diversification of snake venom proteins, in addition to accelerated point mutations in the protein coding regions. Accelerated point mutations lead to fine-tuning of target specificity, whereas ASSET leads to large-scale replacement of multiple functionally important residues, resulting in change or gain of functions.

Evolution of soldier-specific venomous protease in social aphids

In social aphids of the genus Tuberaphis a cysteine protease gene of the family cathepsin B exhibits soldier-specific expression and intestinal protease production. The product is orally excreted and injected by soldier nymphs into natural enemies, thereby exerting an insecticidal activity. In an attempt to gain insights into when and how the novel venomous protease for the altruistic caste has evolved, we investigated the soldier-specific type (S-type) and nonspecific type (N-type) cathepsin B genes from social and nonsocial aphids. All the social aphids examined, representing the genera Tuberaphis, Astegopteryx, and Cerataphis, possessed both the S-type and N-type genes. Phylogenetically distant nonsocial aphids also possessed cathepsin B genes allied to the S-type and the N-type, indicating the evolutionary origin of these genes in the common ancestor of extant aphids. In Tuberaphis species the S-type genes exhibited significant soldier-specific expression and accelerated molecular evolution whereas the N-type genes did not. In Astegopteryx and Cerataphis species, meanwhile, both the S-type and N-type genes exhibited neither remarkable soldier-specific expression nor accelerated molecular evolution. These results suggest that the S-type gene acquired the soldier-specific expression and the venom function after divergence of the genus Tuberaphis. On the structural model of the S-type protease of Tuberaphis styraci the accelerated molecular evolution was associated with the molecular surface rather than the catalytic cleft, suggesting that the venom activity was probably acquired by relatively minor modifications on the molecular surface rather than by generation of a novel active site. In Cerataphis jamuritsu the S-type gene was, although containing a stop codon, structurally almost intact and still transcribed, suggesting recent pseudogenization of the gene copy and possible relevance to relaxed functional constraint in the highly multiplied protease gene family. On the basis of these results we suggest that the massive amplification in aphid cathepsin B genes might have predisposed the evolution of venomous protease in the social aphid lineage and argue that gene duplication, accelerated molecular evolution, and acquisition of novel gene function must have played considerable roles in the evolution of complex biological systems including insect sociality.

Accelerated exchange of exon segments in Viperid three-finger toxin genes (Sistrurus catenatus edwardsii; Desert Massasauga)

Background

Snake venoms consist primarily of proteins and peptides showing a myriad of potent biological activities which have been shaped by both adaptive and neutral selective forces. Venom proteins are encoded by multigene families that have evolved through a process of gene duplication followed by accelerated evolution in the protein coding region.

Results

Here we report five gene structures of three-finger toxins from a viperid snake, Sistrurus catenatus edwardsii. These toxin genes are structured similarly to elapid and hydrophiid three-finger toxin genes, with two introns and three exons. Both introns and exons show distinct patterns of segmentation, and the insertion/deletion of segments may define their evolutionary history. The segments in introns, when present, are highly similar to their corresponding segments in other members of the gene family. In contrast, some segments in the exons show high similarity, while others are often distinctly different among corresponding regions of the isoforms.

Conclusion

Ordered, conserved exon structure strongly suggests that segments in corresponding regions in exons have been exchanged with distinctly different ones during the evolution of these genes. Such a "switching" of segments in exons may result in drastically altering the molecular surface topology and charge, and hence the molecular targets of these three-finger toxins. Thus the phenomenon of accelerated segment switch in exons to alter targeting (ASSET) may play an important role in the evolution of three-finger toxins, resulting in a family of toxins with a highly conserved structural fold but widely varying biological activities.

Molecular evolution of toxin genes in Elapidae snakes

The venom of the sea krait, Laticauda semifasciata, consists primarily of two toxic proteins, phospholipase A(2) (PLA(2)) and a three-finger-structure toxin. We have cloned both toxic protein genes, including the upstream region. PLA(2) genes contain three types of inserted sequences: an AG-rich region, a chicken repeat 1-like long interspersed nucleotide element sequence and an intron II 3' side repeat sequence. The molecular divergence of L. semifasciata PLA(2) genes was defined on the basis of the inserted sequences and their sequence homology. The length of intron I in the three-finger-structure toxin genes differs from species to species. The alignment analysis of intron I of the three-finger-structure toxin genes revealed that the intron I sequence of the ancestral gene comprised ten genetic regions. A hypothetical evolutionary process for the three-finger-structure toxin genes has also been developed.

Venom Proteomes of Closely RelatedSistrurusRattlesnakes with Divergent Diets

The protein composition of the venoms of the three subspecies of Sistrurus catenatus (S. c. catenatus, tergeminus, and edwardsii) and a basal species, Sistrurus miliarius barbouri, were analyzed by RP-HPLC, N-terminal sequencing, MALDI-TOF peptide mass fingerprinting, and CID-MS/MS. The venoms of the four Sistrurus taxa contain proteins from 11 families. The protein family profile and the relative abundance of each protein group in the different venoms are not conserved. Myotoxins and 2-chain PLA2s were detected only in S.c. catenatus and S.c. tergeminus, whereas C-type BPP and Kunitz-type inhibitors were exclusively found in S.c. edwardsii and Sistrurus miliarius barbouri. Among major protein families, taxa were most similar in their metalloproteases (protein similarity coefficient value: 34%) and most divergent in PLA2s (12%), with values for disintegrins and serine proteases lying between these extremes (25 and 20%, respectively). The patterns of venom diversity points to either a gain in complexity in S. catenatus taxa or a loss of venom diversity occurring early on in the evolution of the group involving the lineage connecting S. milarius to the other taxa. The high degree of differentiation in the venom proteome among recently evolved congeneric taxa emphasizes the uniqueness of the venom composition of even closely related species that have different diets. Comparative proteomic analysis of Sistrurus venoms provides a comprehensive catalog of secreted proteins, which may contribute to a deeper understanding of the biology and ecology of these North American snakes and may also serve as a starting point for studying structure-function correlations of individual toxins.

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 purification of α-neurotoxins from king cobra (Ophiophagus hannah)

Thirteen complete and three partial cDNA sequences were cloned from the constructed king cobra (Ophiophagus hannah) venom gland cDNA library. Phylogenetic analysis of nucleotide sequences of king cobra with those from other snake venoms revealed that obtained cDNAs are highly homologous to snake venom alpha-neurotoxins. Alignment of deduced mature peptide sequences of the obtained clones with those of other reported alpha-neurotoxins from the king cobra venom indicates that our obtained 16 clones belong to long-chain neurotoxins (seven), short-chain neurotoxins (seven), weak toxin (one) and variant (one), respectively. Up to now, two out of 16 newly cloned king cobra alpha-neurotoxins have identical amino acid sequences with CM-11 and Oh-6A/6B, which have been characterized from the same venom. Furthermore, five long-chain alpha-neurotoxins and two short-chain alpha-neurotoxins were purified from crude venom and their N-terminal amino acid sequences were determined. The cDNAs encoding the putative precursors of the purified native peptide were also determined based on the N-terminal amino acid sequencing. The purified alpha-neurotoxins showed different lethal activities on mice.

Molecular cloning, expression and characterization of three short chain alpha-neurotoxins from the venom of sea snake--Hydrophiinae Hydrophis cyanocinctus Daudin

Three different genes named sn311, sn316 and sn285 were discovered by large-scale randomly sequencing the high quality cDNA library of the venom glands from Hydrophiinae Hydrophis cyanocinctus Daudin. Sequence analysis showed that these three genes encoded three different short chain alpha-neurotoxins of 81 amino acids, which contained a signal peptide of 21 amino acids and followed by a mature peptide of 60 amino acids. Amino acid comparison reveals that mature peptides of sn311 and sn316 are highly homologous, with the only variance at position 46, which is Lys46 and Ser46, respectively. Whereas the mature peptide of sn285 lacks the most conserved amino acids in short chain alpha-neurotoxins, Asp31 and Arg33. The coding sequences of three neurotoxins were cloned into a thioredoxin (TRX) fusion expression vector (pTRX) and expressed as soluble recombinant fusion proteins in E. coli. After purification, approximately 10 mg/l recombinant proteins with the purity up to 95% were obtained. These three recombinant proteins are designated as rSN311, rSN316 and rSN285, they have a molecular weight of 6.963, 6.920 and 6.756 kDa, respectively, which are similar to those predicted from amino acid sequences. LD50 values of rSN311, rSN316 and rSN285 are 0.0827, 0.095, and 0.0647 mg/kg to mice, respectively. Studies on effects of these recombinant proteins on neuromuscular transmission were carried out, and results indicate that they all can produce prompt blockade of neuromuscular transmission, but display distinct biological activity characteristic individually. The results from UV-circular dichroism (CD) spectra indicate that they share similar secondary structure compared to other identified alpha-neurotoxins, and no significant structural differences in these recombinant proteins are observed.

Structure and chromosomal localization of the gene for crotamine, a toxin from the South American rattlesnake, Crotalus durissus terrificus

Crotamine is a 42 amino acid-long basic polypeptide, one of the major components of the South American rattlesnake, Crotalus durissus terrificus, venom. The mRNA has about 340 nucleotides and codifies a pre-crotamine, including the signal peptide, the mature crotamine, and a final lysine. In this report, we describe the crotamine gene with 1.8 kb organized into three exons separated by a long phase-1 (900 bp) and a short phase-2 (140 bp) introns. Exon 1 includes the 5'-untranslated region and codifies the first 19 amino acids of the signal peptide. Exon 2 codifies 42 amino acids, three belonging to the signal peptide and 39 to the mature crotamine. Exon 3 codifies the last three amino acids of the mature toxin and the terminal lysine. The crotamine gene was mapped by in situ hybridization to the end of the long arm of chromosome 2, the intensity of signals differing between the two homologues. This may reflect a difference in gene copy numbers between chromosomes, a possible explanation for the variable amounts of crotamine found in the venom.

Snake postsynaptic neurotoxins: gene structure, phylogeny and applications in research and therapy

Snake venoms are complex mixtures of biologically active polypeptides that target a variety of vital physiological functions in mammals. alpha-Neurotoxins, toxins that cause paralysis by binding to the nicotinic receptors at the postsynaptic region of the neuromuscular junction have been widely studied in terms of their structure-function relationships as well as gene structure, organization and expression. In this review, we describe the structure of alpha-neurotoxin genes and discuss their evolutionary relationships. Almost all members of neurotoxins have been found to exhibit a common evolutionary origin. The importance of alpha-neurotoxins in therapy and research has also been discussed to highlight their potential applications especially in the area of drug discovery.

Group IB phospholipase A2 from Pseudonaja textilis

Pseudonaja textilis, an Australian Elapid, is known to produce a highly toxic venom. Both protein profiling and N-terminal sequence analysis showed the presence of four new phospholipases A(2) in this venom. Besides being non-lethal, the phospholipase A(2) proteins were found to be moderately active enzymes and they showed procoagulant property. cDNA cloning and characterization indicated the presence of two isoforms of PLA(2) proteins in a single snake, each containing the "pancreatic loop," characteristic of group IB phospholipase A(2). The genomic cloning also confirmed the presence of two genes each containing four exons that are interrupted by three introns. Phylogenetic analysis showed that the venom group IB PLA(2) gene is primitive and could have evolved from the same ancestor as the mammalian and venom group IA PLA(2) genes. In the present study, we report that the Pt-PLA2 gene could be responsible for the production of PL1, 2, and 3 possibly via RNA editing process.

Novel neurotoxins from Taiwan banded krait (Bungarus multicinctus) venom: purification, characterization and gene organization

Two novel neurotoxins BM10-1 and BM10-2 were isolated from Bungarus multicinctus (Taiwan banded krait) venom using the combinations of chromatography on a SP-Sephadex C-25 column and a reverse phase HPLC column. BM10-1 contained 66 amino acid residues including 10 Cys residues, while BM10-2 consisted of 65 amino acid residues with 8 Cys residues. The secondary structure of both BM10-1 and BM10-2 was dominated with beta-sheet, but their gross conformation differed as evidenced by CD spectra and acrylamide quenching studies. BM10-1 inhibited carbachol-induced muscle contraction in a reversible manner and the dose for achieving 50% inhibition was approximately fourfold that of alpha-bungarotoxin. BM10-2 exhibited an irreversible but weak inhibition on carbachol-induced muscle contraction. Sequence alignment of neurotoxins with BM10-1 and BM10-2 suggested that the manner in the manifestation of their activity could be partly elucidated by the residues at toxin second loop. The genomic DNAs encoding BM10-1 and BM10-1-like protein (BM10-1L) were amplified by PCR. The two genes shared virtually identical structural organization and high degree of sequence identity with B. multicinctus neurotoxin genes. Compared to intron sequences of these genes, the protein-coding regions were highly variable. The difference between BM10-1 gene and BM10-1L gene notably arose from the third exon. These results suggest the evolution of B. multicinctus neurotoxins via the path of gene duplication.

Structurally conserved α-neurotoxin genes encode functionally diverse proteins in the venom of Naja sputatrix 1

The structure and organization of the genes encoding the long-chain neurotoxins and four other isoforms of weak neurotoxins in the venom of Naja sputatrix are reported. The genes contained three exons interrupted by two introns, a structure similar to other members of the three-finger toxin family. The proteins encoded by these genes, however, show varied affinity towards nicotinic acetylcholine receptors. Phylogenetic analysis of these genes showed that the weak neurotoxin gene is confined to a distinct group. We also observe that specific mutations of the gene provide the diversity in function in these toxins while maintaining a common structural scaffold. This forms the first report where the molecular basis of evolution of postsynaptic neurotoxins from an ancestral gene can be demonstrated using the same species of snake.

Adaptive evolution in the snake venom Kunitz/BPTI protein family

Snake venoms are rich sources of serine proteinase inhibitors that are members of the Kunitz/BPTI (bovine pancreatic trypsin inhibitor) family. However, only a few of their gene sequences have been determined from snakes. We therefore cloned the cDNAs for the trypsin and chymotrypsin inhibitors from a Vipera ammodytes venom gland cDNA library. Phylogenetic analysis of these and other snake Kunitz/BPTI homologs shows the presence of three clusters, where sequences cluster by functional role. Analysis of the nucleotide sequences from the snake Kunitz/BPTI family shows that positive Darwinian selection was operating on the highly conserved BPTI fold, indicating that this family evolved by gene duplication and rapid diversification.

Muscarinic toxin-like proteins from Taiwan banded krait (Bungarus multicinctus) venom: purification, characterization and gene organization

Two novel proteins, BM8 and BM14, were isolated from Bungarus multicinctus (Taiwan banded krait) venom using the combination of chromatography on a SP-Sephadex C-25 column and a reverse-phase HPLC column. Both proteins contained 82 amino acid residues including 10 cysteine residues, but there were two amino acid substitutions at positions 37 and 38 (Glu37-Ala38 in BM8; Lys37-Lys38 in BM14). CD spectra and acrylamide quenching studies revealed that the gross conformation of BM8 and BM14 differed. In contrast to BM8, BM14 inhibited the binding of [3H]quinuclidinyl benzilate to the M2 muscarinic acetylcholine (mAchR) receptor subtype. Trinitrophenylation of Lys residues abolished the mAchR-binding activity of BM14, indicating that the Lys substitutions at positions 37 and 38 played a crucial role in the activity of BM14. The genomic DNA encoding the precursor of BM14 was amplified by PCR. The gene shared virtually identical structural organization with alpha-neurotoxin and cardiotoxin genes. The intron sequences of these genes shared a sequence identity up to 84%, but the protein-coding regions were highly variable. These results suggest that BM8, BM14, neurotoxins and cardiotoxins may have originated from a common ancestor, and the evolution of snake venom proteins shows a tendency to diversify their functions.

Alpha-neurotoxin gene expression in Naja sputatrix: identification of a silencer element in the promoter region

Alpha-neurotoxin (alpha-NTX) from the venom of cobra, Naja sputatrix, is a highly lethal post-synaptic toxin that is responsible for the lethality caused by the venom. However, this toxin is found at low levels (3%) in the crude venom. The expression of its gene is determined by a promoter which is 90% similar to the promoter of another three-fingered toxin, cardiotoxin (CTX), which is produced in large amounts (60%) in the same venom. Functional analysis of the NTX-2 gene promoter demonstrated the presence of a silencer element of 24 nucleotides (nt -678 to -655) at its 5(') flanking region. This element has been found to play a major role in the down-regulation of NTX-2 gene expression. A point mutation on this silencer appears to attenuate its repressive property in CTX-2 gene.

Characterization and gene organization of Taiwan banded krait (Bungarus multicinctus) gamma-bungarotoxin

gamma-Bungarotoxin was isolated from Bungarus multicinctus (Taiwan banded krait) venom using a combination of chromatography on a SP-Sephadex C-25 column and a reverse-phase high-performance liquid chromatography column. Circular dichroism (CD) measurement revealed that its secondary structure was dominant with beta-sheet structure as is that of snake venom alpha-neurotoxins and cardiotoxins. gamma-Bungarotoxin exhibits activity on inhibiting the binding of [3H]quinuclidinyl benzilate to the M2 muscarinic acetylcholine receptor subtype, and competes weakly with radioiodinated alpha-bungarotoxin for binding to the Torpedo nicotinic acetylcholine receptor. Moreover, the toxin inhibits collagen-induced platelet aggregation, with an IC50 of approximately 200 nM. The genomic DNA encoding the gamma-bungarotoxin precursor is amplified by polymerase chain reaction (PCR). The gene is organized with three exons separated by two introns, and shares virtually identical overall organization with those reported for alpha-neurotoxin and cardiotoxin genes, including similar intron insertions. The intron sequences of these genes share sequence identity up to 85%, but the exon sequences are highly variable. These observations suggest that gamma-bungarotoxin, alpha-neurotoxins, and cardiotoxins originate from a common ancestor, and the evolution of these genes shows a tendency to diversify the functions of snake venom proteins.

Cloning and characterization of the pseudonajatoxin b precursor

An Australian common brown snake, Pseudonaja textilis, is known to contain highly lethal neurotoxins. Among them, a long-chain alpha-neurotoxin, pseudonajatoxin b, has been identified. In this report, while presenting evidence for the presence of at least four such long-chain alpha-neurotoxins in the venom of P. textilis, we describe the characteristics of both the mRNA and the gene responsible for the synthesis of these neurotoxins. A precursor toxin synthesized from the gene has been identified as being capable of producing the isoforms possibly by post-translational modifications at its C-terminal end. Recombinant toxins corresponding to the precursor and its product have been found to possess similar binding affinities for muscular acetylcholine receptors (IC(50)=3x10(-8) M) and a lethality, LD(50), of 0.15 microg/g in mice.

Adaptive evolution of animal toxin multigene families

Animal toxins comprise a diverse array of proteins that have a variety of biochemical and pharmacological functions. A large number of animal toxins are encoded by multigene families. From studies of several toxin multigene families at the gene level the picture is emerging that most have been functionally diversified by gene duplication and adaptive evolution. The number of pharmacological activities in most toxin multigene families results from their adaptive evolution. The molecular evolution of animal toxins has been analysed in some multigene families, at both the intraspecies and interspecies levels. In most toxin multigene families, the rate of non-synonymous to synonymous substitutions (dN/dS) is higher than one. Thus natural selection has acted to diversify coding sequences and consequently the toxin functions. The selection pressure for the rapid adaptive evolution of animal toxins is the need for quick immobilization of the prey in classical predator and prey interactions. Currently available evidence for adaptive evolution in animal toxin multigene families will be considered in this review.