Structure and organization of the cardiotoxin genes in Naja naja sputatrix

Origins, genomic structure and copy number variation of snake venom myotoxins

Crotamine, myotoxin a and homologs are short peptides that often comprise major fractions of rattlesnake venoms and have been extensively studied for their bioactive properties. These toxins are thought to be important for rapidly immobilizing mammalian prey and are implicated in serious, and sometimes fatal, responses to envenomation in humans. While high quality reference genomes for multiple venomous snakes are available, the loci that encode myotoxins have not been successfully assembled in any existing genome assembly. Here, we integrate new and existing genomic and transcriptomic data from the Prairie Rattlesnake (Crotalus viridis viridis) to reconstruct, characterize, and infer the chromosomal locations of myotoxin-encoding loci. We integrate long-read transcriptomics (Pacific Bioscience's Iso-Seq) and short-read RNA-seq to infer gene sequence diversity and characterize patterns of myotoxin and paralogous β-defensin expression across multiple tissues. We also identify two long non-coding RNA sequences which both encode functional myotoxins, demonstrating a newly discovered source of venom coding sequence diversity. We also integrate long-range mate-pair chromatin contact data and linked-read sequencing to infer the structure and chromosomal locations of the three myotoxin-like loci. Further, we conclude that the venom-associated myotoxin is located on chromosome 1 and is adjacent to non-venom paralogs. Consistent with this locus contributing to venom composition, we find evidence that the promoter of this gene is selectively open in venom gland tissue and contains transcription factor binding sites implicated in broad trans-regulatory pathways that regulate snake venoms. This study provides the best genomic reconstruction of myotoxin loci to date and raises questions about the physiological roles and interplay between myotoxin and related genes, as well as the genomic origins of snake venom variation.

Toxinology provides multidirectional and multidimensional opportunities: A personal perspective

In nature, toxins have evolved as weapons to capture and subdue the prey or to counter predators or competitors. When they are inadvertently injected into humans, they cause symptoms ranging from mild discomfort to debilitation and death. Toxinology is the science of studying venoms and toxins that are produced by a wide variety of organisms. In the past, the structure, function and mechanisms of most abundant and/or most toxic components were characterized to understand and to develop strategies to neutralize their toxicity. With recent technical advances, we are able to evaluate and determine the toxin profiles using transcriptomes of venom glands and proteomes of tiny amounts of venom. Enormous amounts of data from these studies have opened tremendous opportunities in many directions of basic and applied research. The lower costs for profiling venoms will further fuel the expansion of toxin database, which in turn will provide greater exciting and bright opportunities in toxin research.

Novel long-chain neurotoxins from Bungarus candidus distinguish the two binding sites in muscle-type nicotinic acetylcholine receptors

αδ-Bungarotoxins, a novel group of long-chain α-neurotoxins, manifest different affinity to two agonist/competitive antagonist binding sites of muscle-type nicotinic acetylcholine receptors (nAChRs), being more active at the interface of α-δ subunits. Three isoforms (αδ-BgTx-1-3) were identified in Malayan Krait (Bungarus candidus) from Thailand by genomic DNA analysis; two of them (αδ-BgTx-1 and 2) were isolated from its venom. The toxins comprise 73 amino acid residues and 5 disulfide bridges, being homologous to α-bungarotoxin (α-BgTx), a classical blocker of muscle-type and neuronal α7, α8, and α9α10 nAChRs. The toxicity of αδ-BgTx-1 (LD₅₀ = 0.17-0.28 µg/g mouse, i.p. injection) is essentially as high as that of α-BgTx. In the chick biventer cervicis nerve-muscle preparation, αδ-BgTx-1 completely abolished acetylcholine response, but in contrast with the block by α-BgTx, acetylcholine response was fully reversible by washing. αδ-BgTxs, similar to α-BgTx, bind with high affinity to α7 and muscle-type nAChRs. However, the major difference of αδ-BgTxs from α-BgTx and other naturally occurring α-neurotoxins is that αδ-BgTxs discriminate the two binding sites in the Torpedo californica and mouse muscle nAChRs showing up to two orders of magnitude higher affinity for the α-δ site as compared with α-ε or α-γ binding site interfaces. Molecular modeling and analysis of the literature provided possible explanations for these differences in binding mode; one of the probable reasons being the lower content of positively charged residues in αδ-BgTxs. Thus, αδ-BgTxs are new tools for studies on nAChRs.

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.

Purification and characterization of a novel antinociceptive toxin from Cobra venom (Naja naja atra)

Snake venoms have demonstrated antinociceptive activity, and certain isolated neurotoxins have demonstrated significant analgesia in animal models. Here we report a novel analgesic toxin which was isolated from Naja naja atra and was given the name 'najanalgesin'. The LD(50) of the crude venom and najanalgesin were 0.89mg/kg and 2.69mg/kg, respectively. We used the writhing test and hot plate test to evaluate the antinociceptive properties of the crude venom and najanalgesin after intraperitoneal (ip) administration. The analgesic mechanism of najanalgesin was also studied. The response latency time was significantly prolonged in the hot plate test after ip administration of the crude venom of Naja naja atra (0.111-0.445mg/kg) in a dose-dependent manner. Najanalgesin (1mg/kg) elicited almost the same antinociceptive effect as that of the crude venom of Naja naja atra at the dose of 0.445mg/kg and remained for 6h after intraperitoneal injection, shown by hot plate test. The percentage of increase in the latency time for the venom and the najanalgesin 3h after drug administration was 96.2% and 112%, respectively. The number of writhes decreased to almost 1/3, 1/6, and 1/12 of the NS (physiological saline) group after intraperitoneal administration of najanalgesin at 0.25, 0.5, and 1.0mg/kg, respectively. Pretreatment with atropine (1mg/kg) or naloxone (3mg/kg) blocked the antinociception of najanalgesin in the hot plate test. Based on the sequence information, najanalgesin is found to be highly homologous with the conventional CTXs (cardiotoxins). To our knowledge, no study had previously reported that a toxin which was homologous with CTXs possessed the antinociceptive activity. Thus, this is the first report that the antinociceptive effect induced by najanalgesin is mediated by cholinergic and opioidergic mechanisms.

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.

Genetic organization of Bungarus multicinctus protease inhibitor-like proteins

The structural organization of the genes encoding Bungarus multicinctus protease inhibitor-like proteins (PILPs), PILP-1, PILP-2 and PILP-3, are reported in this study. Unlike PILP-2 and PILP-3, recombinant PILP-1 exhibited inhibitory activity on trypsin. PILP genes and B chain genes shared identical organization with three exons interrupted by two introns in similar positions. On the contrary, intron 1 of these genes had a similar size, a notable variation with the size of intron 2 was observed. It was found that two regions at the second intron of B1 chain and B2 chain genes were absent in that of PILP genes. Noticeably, intronic insertion in the second intron of B chain genes appeared in the promoter region of PILP-1 gene, but not in that of PILP-2 and PILP-3 genes. Comparative analyses of PILP genes and B chain genes showed that the protein-coding regions of the exons are more diverse than introns, except for in the signal peptide domain. These results suggest that PILP genes and B chain genes originate from a common ancestor, and that accelerated evolution may diversify PILP and B chain genes as that proposed for snake venom phospholipase A(2), neurotoxin and cardiotoxin genes.

Divergence of genes encoding B chains of beta-bungarotoxins

The structural organization of the genes encoding B2, B4, B5 and B6 chains of beta-bungarotoxins are reported in this study. These genes shared virtually identical overall organization with three exons interrupted by two introns in similar positions. On the contrary, intron 1 of these genes had a similar size, a notable variation with the size of intron 2 was observed. It was found that two regions at the second intron of B1 and B2 chains were absent in that of B4, B5 and B6 chains. RT-PCR analyses indicated that Bungarus multicinctus venom gland, heart, liver and muscle expressed the RNA transcripts showing sequence similarity with the intronic segment being deleted in B4, B5 and B6 chain genes. This reflects that the ancestral gene of the intronic segment might insert in multiple loci of B. multicinctus genome. Comparative analyses of B chain genes showed that the protein-coding regions of the exons are more diverse than introns, except for in the signal peptide domain. These results suggest that intron insertions or deletions occur with the evolution of B chains, and that accelerated evolution may diversify the protein-coding sequence of B chain genes same as snake phospholipase A2, neurotoxin and cardiotoxin genes.

Taiwan cobra chymotrypsin inhibitor: cloning, functional expression and gene organization

A cDNA encoding chymotrypsin inhibitor was constructed from the cellular RNA isolated from the venom glands of Naja atra (Taiwan cobra). The resultant amino acid sequence showed that the mature protein is comprised of 57 amino acid residues with six cysteine residues. Cloned protein was expressed and isolated from the inclusion bodies of E. coli and refolded into a functional protein in vitro. Deleting the first three residues at its N-terminus caused a moderate increase in the inhibitory constant (K(i)) against chymotrypsin. The genomic DNA encoding the chymotrypsin inhibitor was amplified by PCR. The gene shares virtually an identical structural organization with the beta-bungarotoxin B1 chain (a snake Kunitz/BPTI neurotoxic homolog) gene. Moreover, the overall sequence identity of the N. atra chymotrypsin inhibitor and beta-bungarotoxin B1 chain genes was up to 83%. These findings strongly suggest that snake Kunitz/BPTI protease inhibitors and neurotoxic homologs may have originated from a common ancestor.

The speciation of conger eel galectins by rapid adaptive evolution

Many cases of accelerated evolution driven by positive Darwinian selection are identified in the genes of venomous and reproductive proteins. This evolutional phenomenon might have important consequences in their gene-products' functions, such as multiple specific toxins for quick immobilization of the prey and the establishment of barriers to fertilization that might lead to speciation, and in the molecular evolution of novel genes. Recently, we analyzed the molecular evolution of two galectins isolated from the skin mucus of conger eel (Conger myriaster), named congerins I and II, by cDNA cloning and X-ray structural analysis, and we found that they have evolved in the rapid adaptive manner to emergence of a new structure including strand-swapping and a unique new ligand-binding site. In this review article we summarize and discuss the molecular evolution, especially the rapid adaptive evolution, and the structure-function relationships of conger eel galectins.

Identification of crotasin, a crotamine-related gene of Crotalus durissus terrificus

Crotamine is a cationic peptide (4.9 kDa, pI 9.5) of South American rattlesnake, Crotalus durissus terrificus' venom. Its presence varies according to the subspecies or the geographical locality of a given species. At the genomic level, we observed the presence of 1.8 kb gene, Crt-p1, in crotamine-positive specimens and its absence in crotamine-negative ones. In this work, we described a crotamine-related 2.5 kb gene, crotasin (Cts-p2), isolated from crotamine-negative specimens. Reverse transcription coupled to polymerase chain reaction indicates that Cts-p2 is abundantly expressed in several snake tissues, but scarcely expressed in the venom gland. The genome of crotamine-positive specimen contains both Crt-p1 and Cts-p2 genes. The present data suggest that both crotamine and crotasin have evolved by duplication of a common ancestor gene, and the conservation of their three disulfide bonds indicates that they might adopt the same fold as beta-defensin. The physiological function of the crotasin is not yet known.

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.

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.

Identification of alpha-bungarotoxin (A31) as the major postsynaptic neurotoxin, and complete nucleotide identity of a genomic DNA of Bungarus candidus from Java with exons of the Bungarus multicinctus alpha-bungarotoxin (A31) gene

The Malayan krait (Bungarus candidus) is one of the most medically significant snake species in Southeast Asia. No specific antivenom exists to treat envenoming by this species. Death within 30 min after its bite has been reported from Java, suggesting the presence of highly lethal postsynaptic neurotoxins in the venom of these snakes. We purified and identified the major postsynaptic toxin in the venom of B. candidus from Java. The toxin was indistinguishable from alpha-bungarotoxin (A31), a toxin originally isolated from Bungarus multicinctus, in its mass (7983.75 Da), LD50 (0.23 microg/g in mice i.p.), affinity to nicotinic acetylcholine receptors, and by its 40 N-terminal amino acid residues as determined by Edman degradation. Identity with alpha-bungarotoxin was confirmed by cloning and sequencing a genomic DNA from B. candidus which encodes the 74 amino acid sequence of alpha-bungarotoxin (A31) and part of its signal peptide, revealing complete identity to the alpha-bungarotoxin (A31) gene in exon and 98.9% identity in intron sequences. The entire mitochondrial cytochrome b gene of the krait species B. candidus from Java and B. multicinctus from Taiwan was sequenced for comparison, suggesting that these snakes are phylogenetically closely related. alpha-Bungarotoxin appears to be widely present and conserved in Southeast and East Asian black-and-white kraits across populations and taxa.

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.

Purification and characterization of a neurotoxin from the venom of Ophiophagus hannah (king cobra)

A neurotoxin, Oh9-1, from the venom of Ophiophagus hannah was isolated by a combination of ion-exchange chromatography and reverse phase HPLC. Amino acid sequence analysis revealed that Oh9-1 consists of 57 amino acids and eight cysteine residues. This protein was mainly constituted with beta-sheet as evidenced by CD spectrum. Oh9-1 inhibited carbachol-induced muscle contraction in an irreversible manner and the dose for achieving 50% inhibition was approximately fourfold that of alpha-bungarotoxin. Since the residues in alpha-neurotoxins closely involve in the binding with acetylcholine receptors are not highly conserved in this toxin molecule, Oh9-1 represents a novel type of neurotoxin structurally distinct from alpha-neurotoxins.

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.

Purification and Characterization of a Chymotrypsin Inhibitor from the Venom of Ophiophagus hannah (King Cobra)

A chymotrypsin inhibitor from the venom of Ophiophagus hannah was isolated by a combination of ion-exchange chromatography and reverse phase HPLC. Amino acid sequence analysis revealed that this protein consists of 58 amino acids, six of these being cysteine residues and is highly homologous to Kunitz-type protease inhibitors. ESI-mass spectrum showed that the protein had a mass of 6493, which is in agreement with that predicted from its primary structure. In contrast to P1 Leu, Met, Phe, Trp, and Tyr appearing in other chymotrypsin inhibitors, a P1 Asn in the novel inhibitor may cause a weak binding (Ki = 3.52 microM) with chymotrypsin. Phylogenetic analysis suggests that the functional variations of the chymotrypsin inhibitor and other Kunitz-type inhibitors probably distinguish from dendrotoxins by accelerated evolution.

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.

The multiplicity of cardiotoxins from Naja naja atra (Taiwan cobra) venom

Four novel cardiotoxins were isolated from Naja naja atra (Taiwan cobra) venom by successive separation on a SP-Sephadex C-25 column and a reverse phase column. Amino acid sequences of the cardiotoxins were determined by Edman degradation and carboxypeptidase digestion. It shows that these cardiotoxins comprise 60 amino acid residues. Comparative analyses on the amino acid sequences of cardiotoxins from the venoms of N. naja atra and other Naja species indicated that amino acid substitutions of cardiotoxin isoforms frequently occurred at positions 7-11, 27-32 and 45-47. The hypervariable segments encoded by the second and third exon of cardiotoxin genes are located at or near the tips of loop structure of cardiotoxin molecules. These results, together with the suggestions that the residues at the tips of cardiotoxins' loop structure were involved in the manifestation of the biological activities of cardiotoxins, reflect that the preferential mutations may contribute to alterations in the function of cardiotoxin molecules. Analysis on the secondary structure of pre-mRNAs of N. naja atra cardiotoxin 4 gene and N. naja sputatrix cardiotoxin 3 gene has shown that the hypervariable regions of the exon 2 pertain to form intra-exon pairings and are not involved in the formation of intron-exon pairings. Since the pairings of splice sites and gene architecture were supposed to be associated with intron-exon recognition, it is likely that the preferred loci of mutations occurring with the evolution of cardiotoxin genes would not affect the processing of cardiotoxin precursors.

Structure-function studies on Taiwan cobra long neurotoxin homolog

A novel long neurotoxin homolog was purified from Naja naja atra (Taiwan cobra) venom using the combination of ion exchange chromatography and reverse phase high performance liquid chromatography. The determined protein sequence was essentially the same as that deduced from the cDNA amplified by reverse transcriptase-polymerase chain reaction. The long neurotoxin homolog exhibited an activity that inhibited acetylcholine-induced muscle contractions, as with N. naja atra cobrotoxin. The degree of inhibition caused by the addition of long neurotoxin homolog was approximately 70% of that observed with the addition of cobrotoxin. Unlike the well-known short and long neurotoxins, this neurotoxin homolog contained two additional cysteine residues forming a disulfide linkage in the N-terminal region. Circular dichroism measurement and computer models of the neurotoxin reveal that its secondary structure was not abundant in beta-sheet as noted with short and long neurotoxins. This less ordered structure may be associated with the lower activity noted with the long neurotoxin homolog. Together with the finding that the known long neurotoxin homologs exclusively appear in the venoms of the Naja and Bungarus genera, the long neurotoxin homologs should represent an evolutionary branch from the long and short neurotoxins in the Elapidae family.

Structure and phylogeny of the venom group I phospholipase A(2) gene

Phospholipases A(2) (PLA(2)s) catalyzing the hydrolysis of phospholipids form a family of proteins with diverse physiological and pharmacological properties. While there have been several reports on the cloning of PLA(2) cDNAs, very few studies have been carried out on the PLA(2) genes and, most importantly, no information has been available on the gene structure and function of group I venom PLA(2). This study, on the PLA(2) gene from a spitting cobra, besides being the very first report on any venom group I PLA(2) gene, constitutes the missing link in the biology and evolution of phospholipases. The 4-kb gene consists of four exons and three introns and resembles the human pancreatic PLA(2) gene. However, the size of intron 3 in particular is much smaller than that in the pancreatic gene. Interestingly, the information for the toxic and most of the pharmacological properties of the venom PLA(2) can be attributed to the end of exon 3 and the whole of exon 4 of the gene. This functional delineation fits in well with the theory of adaptive evolution exhibited by the venom PLA(2)s. We also show that the mammalian pancreatic and elapid PLA(2)s have similar paths of evolution (probably following gene duplication) from a common ancestral gene. Venom group II phospholipases, although evolved from the same ancestor, diverged early in evolution from the group I PLA(2) genes. Intriguingly, CAT reporter gene assays and DNase 1 footprinting studies on the promoter and its deletion constructs using CHO and HepG2 cell lines identified the possible involvement of cis elements such as Sp1, AP2, gamma-IRE, and (TG)(12) repeats in the expression of the gene in a tissue-specific manner.

Postsynaptic short-chain neurotoxins from Pseudonaja textilis. cDNA cloning, expression and protein characterization

Two lethal proteins, which specifically bind to the nAChR from Torpedo californica, were isolated from the venom of Pseudonaja textilis, the common brown snake from Australia. The isolated proteins have masses of 6236 and 6345 Da and are structurally related to short-chain neurotoxins from other elapids. Six cDNAs encoding isoforms of related neurotoxins were cloned using the RT-PCR of the venom gland mRNAs. The sequences of the corresponding proteins consist of 57-58 amino acid residues and display several unique features when compared with all known short-chain neurotoxins. Accordingly, they grouped separately in phylogenetic analysis. The six cDNAs were expressed in Escherichia coli and the recombinant proteins were characterized. They have similar masses and display similar toxicities and binding constants to the nAChR as the native toxins isolated from the venom. Thus, a new group of short-chain postsynaptic neurotoxins from the venom of an Australian elapid has been characterized.

Postsynaptic α-Neurotoxin Gene of the Spitting Cobra, Naja naja sputatrix: Structure, Organization, and Phylogenetic Analysis

The venom of the spitting cobra, Naja naja sputatrixcontains highly potent α-neurotoxins (NTXs) in addition to phospholipase A2 (PLA2) and cardiotoxin (CTX). In this study, we report the complete characterization of three genes that are responsible for the synthesis of three isoforms of α-NTX in the venom of a single spitting cobra. DNA amplification by long-distance polymerase chain reaction (LD-PCR) and genome walking have provided information on the gene structure including their promoter and 5′ and 3′ UTRs. Each NTX isoform is ∼4 kb in size and contains three exons and two introns. The sequence homology among these isoforms was found to be 99%. Two possible transcription sites were identified by primer extension analysis and they corresponded to the adenine (A) nucleotide at positions +1 and −45. The promoter also contains two TATA boxes and a CCAAT box. Putative binding sites for transcriptional factors AP-2 and GATA are also present. The high percentage of similarity observed among the NTX gene isoforms of N. n. sputatrix as well as with the α-NTX and κ-NTX genes from other land snakes suggests that the NTX gene has probably evolved from a common ancestral gene.

[The genomic DNA sequences reported in this paper have been submitted to GenBank databases under accession nos. AF096999 to AF097001.]