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

Acetylcholine-Binding Protein Affinity Profiling of Neurotoxins in Snake Venoms with Parallel Toxin Identification

Snakebite is considered a concerning issue and a neglected tropical disease. Three-finger toxins (3FTxs) in snake venoms primarily cause neurotoxic effects since they have high affinity for nicotinic acetylcholine receptors (nAChRs). Their small molecular size makes 3FTxs weakly immunogenic and therefore not appropriately targeted by current antivenoms. This study aims at presenting and applying an analytical method for investigating the therapeutic potential of the acetylcholine-binding protein (AChBP), an efficient nAChR mimic that can capture 3FTxs, for alternative treatment of elapid snakebites. In this analytical methodology, snake venom toxins were separated and characterised using high-performance liquid chromatography coupled with mass spectrometry (HPLC-MS) and high-throughput venomics. By subsequent nanofractionation analytics, binding profiling of toxins to the AChBP was achieved with a post-column plate reader-based fluorescence-enhancement ligand displacement bioassay. The integrated method was established and applied to profiling venoms of six elapid snakes (Naja mossambica, Ophiophagus hannah, Dendroaspis polylepis, Naja kaouthia, Naja haje and Bungarus multicinctus). The methodology demonstrated that the AChBP is able to effectively bind long-chain 3FTxs with relatively high affinity, but has low or no binding affinity towards short-chain 3FTxs, and as such provides an efficient analytical platform to investigate binding affinity of 3FTxs to the AChBP and mutants thereof and to rapidly identify bound toxins.

Unveiling the Venom Composition of the Colombian Coral Snakes Micrurus helleri, M. medemi, and M. sangilensis

Little is known of the biochemical composition and functional features of the venoms of poorly known Colombian coral snakes. Here, we provide a preliminary characterization of the venom of two Colombian endemic coral snake species, Micrurus medemi and M. sangilensis, as well as Colombian populations of M. helleri. Electrophoresis and RP-HPLC techniques were used to identify venom components, and assays were conducted to detect enzyme activities, including phospholipase A2, hyaluronidase, and protease activities. The median lethal dose was determined using murine models. Cytotoxic activities in primary cultures from hippocampal neurons and cancer cell lines were evaluated. The venom profiles revealed similarities in electrophoretic separation among proteins under 20 kDa. The differences in chromatographic profiles were significant, mainly between the fractions containing medium-/large-sized and hydrophobic proteins; this was corroborated by a proteomic analysis which showed the expected composition of neurotoxins from the PLA2 (~38%) and 3FTx (~17%) families; however, a considerable quantity of metalloproteinases (~12%) was detected. PLA2 activity and protease activity were higher in M. helleri venom according to qualitative and quantitative assays. M. medemi venom had the highest lethality. All venoms decreased cell viability when tested on tumoral cell cultures, and M. helleri venom had the highest activity in neuronal primary culture. These preliminary studies shed light on the venoms of understudied coral snakes and broaden the range of sources that could be used for subsequent investigations of components with applications to specific diseases. Our findings also have implications for the clinical manifestations of snake envenoming and improvements in its medical management.

Developing a computational pharmacokinetic model of systemic snakebite envenomation and antivenom treatment

Snakebite envenomation is responsible for over 100,000 deaths and 400,000 cases of disability annually, most of which are preventable through access to safe and effective antivenoms. Snake venom toxins span a wide molecular weight range, influencing their absorption, distribution, and elimination within the body. In recent years, a range of scaffolds have been applied to antivenom development. These scaffolds similarly span a wide molecular weight range and subsequently display diverse pharmacokinetic behaviours. Computational simulations represent a powerful tool to explore the interplay between these varied antivenom scaffolds and venoms, to assess whether a pharmacokinetically optimal antivenom exists. The purpose of this study was to establish a computational model of systemic snakebite envenomation and treatment, for the quantitative assessment and comparison of conventional and next-generation antivenoms. A two-compartment mathematical model of envenomation and treatment was defined and the system was parameterised using existing data from rabbits. Elimination and biodistribution parameters were regressed against molecular weight to predict the dynamics of IgG, F(ab')₂, Fab, scFv, and nanobody antivenoms, spanning a size range of 15-150 kDa. As a case study, intramuscular envenomation by Naja sumatrana (equatorial spitting cobra) and its treatment using Fab, F(ab')₂, and IgG antivenoms was simulated. Variable venom dose tests were applied to visualise effective antivenom dose levels. Comparisons to existing antivenoms and experimental rescue studies highlight the large dose reductions that could result from recombinant antivenom use. This study represents the first comparative in silico model of snakebite envenomation and treatment.

Strategies for Heterologous Expression, Synthesis, and Purification of Animal Venom Toxins

Animal venoms are complex mixtures containing peptides and proteins known as toxins, which are responsible for the deleterious effect of envenomations. Across the animal Kingdom, toxin diversity is enormous, and the ability to understand the biochemical mechanisms governing toxicity is not only relevant for the development of better envenomation therapies, but also for exploiting toxin bioactivities for therapeutic or biotechnological purposes. Most of toxinology research has relied on obtaining the toxins from crude venoms; however, some toxins are difficult to obtain because the venomous animal is endangered, does not thrive in captivity, produces only a small amount of venom, is difficult to milk, or only produces low amounts of the toxin of interest. Heterologous expression of toxins enables the production of sufficient amounts to unlock the biotechnological potential of these bioactive proteins. Moreover, heterologous expression ensures homogeneity, avoids cross-contamination with other venom components, and circumvents the use of crude venom. Heterologous expression is also not only restricted to natural toxins, but allows for the design of toxins with special properties or can take advantage of the increasing amount of transcriptomics and genomics data, enabling the expression of dormant toxin genes. The main challenge when producing toxins is obtaining properly folded proteins with a correct disulfide pattern that ensures the activity of the toxin of interest. This review presents the strategies that can be used to express toxins in bacteria, yeast, insect cells, or mammalian cells, as well as synthetic approaches that do not involve cells, such as cell-free biosynthesis and peptide synthesis. This is accompanied by an overview of the main advantages and drawbacks of these different systems for producing toxins, as well as a discussion of the biosafety considerations that need to be made when working with highly bioactive proteins.

Development of a Broad-Spectrum Antiserum against Cobra Venoms Using Recombinant Three-Finger Toxins

Three-finger toxins (3FTXs) are the most clinically relevant components in cobra (genus Naja) venoms. Administration of the antivenom is the recommended treatment for the snakebite envenomings, while the efficacy to cross-neutralize the different cobra species is typically limited, which is presumably due to intra-specific variation of the 3FTXs composition in cobra venoms. Targeting the clinically relevant venom components has been considered as an important factor for novel antivenom design. Here, we used the recombinant type of long-chain α-neurotoxins (P01391), short-chain α-neurotoxins (P60770), and cardiotoxin A3 (P60301) to generate a new immunogen formulation and investigated the potency of the resulting antiserum against the venom lethality of three medially important cobras in Asia, including the Thai monocled cobra (Naja kaouthia), the Taiwan cobra (Naja atra), and the Thai spitting cobra (Naja Siamensis) snake species. With the fusion of protein disulfide isomerase and the low-temperature settings, the correct disulfide bonds were built on these recombinant 3FTXs (r3FTXs), which were confirmed by the circular dichroism spectra and tandem mass spectrometry. Immunization with r3FTX was able to induce the specific antibody response to the native 3FTXs in cobra venoms. Furthermore, the horse and rabbit antiserum raised by the r3FTX mixture is able to neutralize the venom lethality of the selected three medically important cobras. Thus, the study demonstrated that the r3FTXs are potential immunogens in the development of novel antivenom with broad neutralization activity for the therapeutic treatment of victims involving cobra snakes in countries.

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.

King Cobra and snakebite envenomation: on the natural history, human-snake relationship and medical importance of Ophiophagus hannah

King Cobra (Ophiophagus hannah) has a significant place in many cultures, and is a medically important venomous snake in the world. Envenomation by this snake is highly lethal, manifested mainly by neurotoxicity and local tissue damage. King Cobra may be part of a larger species complex, and is widely distributed across Southeast Asia, southern China, northern and eastern regions as well as the Western Ghats of India, indicating potential geographical variation in venom composition. There is, however, only one species-specific King Cobra antivenom available worldwide that is produced in Thailand, using venom from the snake of Thai origin. Issues relating to the management of King Cobra envenomation (e.g., variation in the composition and toxicity of the venom, limited availability and efficacy of antivenom), and challenges faced in the research of venom (in particular proteomics), are rarely addressed. This article reviews the natural history and sociocultural importance of King Cobra, cases of snakebite envenomation caused by this species, current practice of management (preclinical and clinical), and major toxinological studies of the venom with a focus on venom proteomics, toxicity and neutralization. Unfortunately, epidemiological data of King Cobra bite is scarce, and venom proteomes reported in various studies revealed marked discrepancies in details. Challenges, such as inconsistency in snake venom sampling, varying methodology of proteomic analysis, lack of mechanistic and antivenomic studies, and controversy surrounding antivenom use in treating King Cobra envenomation are herein discussed. Future directions are proposed, including the effort to establish a standard, comprehensive Pan-Asian proteomic database of King Cobra venom, from which the venom variation can be determined. Research should be undertaken to characterize the toxin antigenicity, and to develop an antivenom with improved efficacy and wider geographical utility. The endeavors are aligned with the WHO´s roadmap that aims to reduce the disease burden of snakebite by 50% before 2030.

Peptide Fraction pOh2 Exerts Antiadipogenic Activity through Inhibition of C/EBP-α and PPAR-γ Expression in 3T3-L1 Adipocytes

Many studies have comprehensively examined the venom of Ophiophagus hannah snake. Its venom comprises different compounds exhibiting a wide range of pharmacological activities. In this investigation, four peptide fractions (PFs), ranging from 3 kDa to 10 kDa, isolated from the Vietnamese snake venom of O. hannah were separated by HPLC and investigated for their inhibitory activity on adipogenesis in 3T3-L1 adipocytes. The most effective PF was then further purified, generating two peptides, pOh1 and pOh2. Upon investigation of these two peptides on 3T3-L1 adipocytes, it was revealed that, at 10 μg/mL, pOh2 was able to inhibit the lipid accumulation in 3T3-L1 adipocytes by up to 56%, without affecting cell viability. Furthermore, the pOh2 downregulated the gene expression of important transcription factors C/EBP-α and PPAR-γ. In addition, aP2 and GPDH adipocyte-specific markers were also significantly reduced compared to untreated differentiated cells. Taken together, pOh2 inhibited the expression of key transcription factors C/EBP-α and PPAR-γ and their target genes, aP2 and GPDH, thereby blocking the adipocyte differentiation. In conclusion, this novel class of peptide might have potential for in vivo antiobesity effects.

The Snake with the Scorpion's Sting: Novel Three-Finger Toxin Sodium Channel Activators from the Venom of the Long-Glanded Blue Coral Snake (Calliophis bivirgatus)

Millions of years of evolution have fine-tuned the ability of venom peptides to rapidly incapacitate both prey and potential predators. Toxicofera reptiles are characterized by serous-secreting mandibular or maxillary glands with heightened levels of protein expression. These glands are the core anatomical components of the toxicoferan venom system, which exists in myriad points along an evolutionary continuum. Neofunctionalisation of toxins is facilitated by positive selection at functional hotspots on the ancestral protein and venom proteins have undergone dynamic diversification in helodermatid and varanid lizards as well as advanced snakes. A spectacular point on the venom system continuum is the long-glanded blue coral snake (Calliophis bivirgatus), a specialist feeder that preys on fast moving, venomous snakes which have both a high likelihood of prey escape but also represent significant danger to the predator itself. The maxillary venom glands of C. bivirgatus extend one quarter of the snake's body length and nestle within the rib cavity. Despite the snake's notoriety its venom has remained largely unstudied. Here we show that the venom uniquely produces spastic paralysis, in contrast to the flaccid paralysis typically produced by neurotoxic snake venoms. The toxin responsible, which we have called calliotoxin (δ-elapitoxin-Cb1a), is a three-finger toxin (3FTx). Calliotoxin shifts the voltage-dependence of NaV1.4 activation to more hyperpolarised potentials, inhibits inactivation, and produces large ramp currents, consistent with its profound effects on contractile force in an isolated skeletal muscle preparation. Voltage-gated sodium channels (NaV) are a particularly attractive pharmacological target as they are involved in almost all physiological processes including action potential generation and conduction. Accordingly, venom peptides that interfere with NaV function provide a key defensive and predatory advantage to a range of invertebrate venomous species including cone snails, scorpions, spiders, and anemones. Enhanced activation or delayed inactivation of sodium channels by toxins is associated with the extremely rapid onset of tetanic/excitatory paralysis in envenomed prey animals. A strong selection pressure exists for the evolution of such toxins where there is a high chance of prey escape. However, despite their prevalence in other venomous species, toxins causing delay of sodium channel inhibition have never previously been described in vertebrate venoms. Here we show that NaV modulators, convergent with those of invertebrates, have evolved in the venom of the long-glanded coral snake. Calliotoxin represents a functionally novel class of 3FTx and a structurally novel class of NaV toxins that will provide significant insights into the pharmacology and physiology of NaV. The toxin represents a remarkable case of functional convergence between invertebrate and vertebrate venom systems in response to similar selection pressures. These results underscore the dynamic evolution of the Toxicofera reptile system and reinforces the value of using evolution as a roadmap for biodiscovery.

A Distinct Functional Site in Ω-Neurotoxins: Novel Antagonists of Nicotinic Acetylcholine Receptors from Snake Venom

Snake venom α-neurotoxins from the three-finger toxin (3FTx) family are competitive antagonists with nanomolar affinity and high selectivity for nicotinic acetylcholine receptors (nAChR). Here, we report the characterization of a new group of competitive nAChR antagonists: Ω-neurotoxins. Although they belong to the 3FTx family, the characteristic functional residues of α-neurotoxins are not conserved. We evaluated the subtype specificity and structure-function relationships of Oh9-1, an Ω-neurotoxin from Ophiophagus hannah venom. Recombinant Oh9-1 showed reversible postsynaptic neurotoxicity in the micromolar range. Experiments with different nAChR subtypes expressed in Xenopus oocytes indicated Oh9-1 is selective for rat muscle type α1β1εδ (adult) and α1β1γδ (fetal) and rat neuronal α3β2 subtypes. However, Oh9-1 showed low or no affinity for other human and rat neuronal subtypes. Twelve individual alanine-scan mutants encompassing all three loops of Oh9-1 were evaluated for binding to α1β1εδ and α3β2 subtypes. Oh9-1's loop-II residues (M25, F27) were the most critical for interactions and formed the common binding core. Mutations at T23 and F26 caused a significant loss in activity at α1β1εδ receptors but had no effect on the interaction with the α3β2 subtype. Similarly, mutations at loop-II (H7, K22, H30) and -III (K45) of Oh9-1 had a distinctly different impact on its activity with these subtypes. Thus, Oh9-1 interacts with these nAChRs via distinct residues. Unlike α-neurotoxins, the tip of loop-II is not involved. We reveal a novel mode of interaction, where both sides of the β-strand of Oh9-1's loop-II interact with α1β1εδ, but only one side interacts with α3β2. Phylogenetic analysis revealed functional organization of the Ω-neurotoxins independent of α-neurotoxins. Thus, Ω-neurotoxin: Oh9-1 may be a new, structurally distinct class of 3FTxs that, like α-neurotoxins, antagonize nAChRs. However, Oh9-1 binds to the ACh binding pocket via a different set of functional residues.

Venom-gland transcriptome and venom proteome of the Malaysian king cobra (Ophiophagus hannah)

Background

The king cobra (Ophiophagus hannah) is widely distributed throughout many parts of Asia. This study aims to investigate the complexity of Malaysian Ophiophagus hannah (MOh) venom for a better understanding of king cobra venom variation and its envenoming pathophysiology. The venom gland transcriptome was investigated using the Illumina HiSeq™ platform, while the venom proteome was profiled by 1D-SDS-PAGE-nano-ESI-LCMS/MS.

Results

Transcriptomic results reveal high redundancy of toxin transcripts (3357.36 FPKM/transcript) despite small cluster numbers, implying gene duplication and diversification within restricted protein families. Among the 23 toxin families identified, three-finger toxins (3FTxs) and snake-venom metalloproteases (SVMPs) have the most diverse isoforms. These 2 toxin families are also the most abundantly transcribed, followed in descending order by phospholipases A₂ (PLA₂s), cysteine-rich secretory proteins (CRISPs), Kunitz-type inhibitors (KUNs), and L-amino acid oxidases (LAAOs). Seventeen toxin families exhibited low mRNA expression, including hyaluronidase, DPP-IV and 5’-nucleotidase that were not previously reported in the venom-gland transcriptome of a Balinese O. hannah. On the other hand, the MOh proteome includes 3FTxs, the most abundantly expressed proteins in the venom (43 % toxin sbundance). Within this toxin family, there are 6 long-chain, 5 short-chain and 2 non-conventional 3FTx. Neurotoxins comprise the major 3FTxs in the MOh venom, consistent with rapid neuromuscular paralysis reported in systemic envenoming. The presence of toxic enzymes such as LAAOs, SVMPs and PLA₂ would explain tissue inflammation and necrotising destruction in local envenoming. Dissimilarities in the subtypes and sequences between the neurotoxins of MOh and Naja kaouthia (monocled cobra) are in agreement with the poor cross-neutralization activity of N. kaouthia antivenom used against MOh venom. Besides, the presence of cobra venom factor, nerve growth factors, phosphodiesterase, 5’-nucleotidase, and DPP-IV in the venom proteome suggests its probable hypotensive action in subduing prey.

Conclusion

This study reports the diversity and abundance of toxins in the venom of the Malaysian king cobra (MOh). The results correlate with the pathophysiological actions of MOh venom, and dispute the use of Naja cobra antivenoms to treat MOh envenomation. The findings also provide a deeper insight into venom variations due to geography, which is crucial for the development of a useful pan-regional antivenom.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1828-2) contains supplementary material, which is available to authorized users.

Ophiophagus hannah venom: proteome, components bound by Naja kaouthia antivenin and neutralization by N. kaouthia neurotoxin-specific human ScFv

Venomous snakebites are an important health problem in tropical and subtropical countries. King cobra (Ophiophagushannah) is the largest venomous snake found in South and Southeast Asia. In this study, the O. hannah venom proteome and the venom components cross-reactive to N. kaouthia monospecific antivenin were studied. O. hannah venom consisted of 14 different protein families, including three finger toxins, phospholipases, cysteine-rich secretory proteins, cobra venom factor, muscarinic toxin, L-amino acid oxidase, hypothetical proteins, low cysteine protein, phosphodiesterase, proteases, vespryn toxin, Kunitz, growth factor activators and others (coagulation factor, endonuclease, 5’-nucleotidase). N. kaouthia antivenin recognized several functionally different O. hannah venom proteins and mediated paratherapeutic efficacy by rescuing the O. hannah envenomed mice from lethality. An engineered human ScFv specific to N. kaouthia long neurotoxin (NkLN-HuScFv) cross-neutralized the O. hannah venom and extricated the O. hannah envenomed mice from death in a dose escalation manner. Homology modeling and molecular docking revealed that NkLN-HuScFv interacted with residues in loops 2 and 3 of the neurotoxins of both snake species, which are important for neuronal acetylcholine receptor binding. The data of this study are useful for snakebite treatment when and where the polyspecific antivenin is not available. Because the supply of horse-derived antivenin is limited and the preparation may cause some adverse effects in recipients, a cocktail of recombinant human ScFvs for various toxic venom components shared by different venomous snakes, exemplified by the in vitro produced NkLN-HuScFv in this study, should contribute to a possible future route for an improved alternative to the antivenins.

Functional proteomic approach to discover geographic variations of king cobra venoms from Southeast Asia and China

This study deciphers the geographic variations of king cobra (Ophiophagus hannah) venom using functional proteomics. Pooled samples of king cobra venom (abbreviated as Ohv) were obtained from Indonesia, Malaysia, Thailand, and two provinces of China, namely Guangxi and Hainan. Using two animal models to test and compare the lethal effects, we found that the Chinese Ohvs were more fatal to mice, while the Southeast Asian Ohvs were more fatal to lizards (Eutropis multifasciata). Various phospholipases A2 (PLA2s), three-finger toxins (3FTxs) and Kunitz-type inhibitors were purified from these Ohvs and compared. Besides the two Chinese Ohv PLA2s with known sequences, eight novel PLA2s were identified from the five Ohv samples and their antiplatelet activities were compared. While two 3FTxs (namely oh-55 and oh-27) were common in all the Ohvs, different sets of 3FTx markers were present in the Chinese and Southeast Asian Ohvs. All the Ohvs contain the Kunitz inhibitor, OH-TCI, while only the Chinese Ohvs contain the inhibitor variant, Oh11-1. Relative to the Chinese Ohvs which contained more phospholipases, the Southeast Asian Ohvs had higher metalloproteinase, acetylcholine esterase, and alkaline phosphatase activities.

BIOLOGICAL SIGNIFICANCE

Remarkable variations in five king cobra geographic samples reveal fast evolution and dynamic translational regulation of the venom which probably adapted to different prey ecology as testified by the lethal tests on mice and lizards. Our results predict possible variations of the king cobra envenoming to human and the importance of using local antivenin for snakebite treatment.

Selective targeting of G-protein-coupled receptor subtypes with venom peptides

The G-protein-coupled receptor (GPCR) family is one of the largest gene superfamilies with approx. 370 members responding to endogenous ligands in humans and a roughly equal amount of receptors sensitive to external stimuli from the surrounding. A number of receptors from this superfamily are well recognized targets for medical treatment of various disease conditions, whereas for many others the potential medical benefit of interference is still obscure. A general problem associated with GPCR research and therapeutics is the insufficient specificity of available ligands to differentiate between closely homologous receptor subtypes. In this context, venom peptides could make a significant contribution to the development of more specific drugs. Venoms from certain animals specialized in biochemical hunting contain a mixture of molecules that are directed towards a variety of membrane proteins. Peptide toxins isolated from these mixtures usually exhibit high specificity for their targets. Muscarinic toxins found from mamba snakes attracted much attention during the 1990s. These are 65-66 amino acid long peptides with a structural three-finger folding similar to the α-neurotoxins and they target the muscarinic acetylcholine receptors in a subtype-selective manner. Recently, several members of the three-finger toxins from mamba snakes as well as conotoxins from marine cone snails have been shown to selectively interact with subtypes of adrenergic receptors. In this review, we will discuss the GPCR-directed peptide toxins found from different venoms and how some of these can be useful in exploring specific roles of receptor subtypes.

Directed evolution of a three-finger neurotoxin by using cDNA display yields antagonists as well as agonists of interleukin-6 receptor signaling

Background

Directed evolution of biomolecules such as DNA, RNA and proteins containing high diversity has emerged as an effective method to obtain molecules for various purposes. In the recent past, proteins from non-immunoglobulins have attracted attention as they mimic antibodies with respect to binding potential and provide further potential advantages. In this regard, we have attempted to explore a three-finger neurotoxin protein (3F). 3F proteins are small (~7 kDa), structurally well defined, thermally stable and resistant to proteolysis that presents them as promising candidates for directed evolution.

Results

We have engineered a snake α-neurotoxin that belongs to the 3F family by randomizing the residues in the loops involved in binding with acetylcholine receptors and employing cDNA display to obtain modulators of interleukin-6 receptor (IL-6R). Selected candidates were highly specific for IL-6R with dissociation constants and IC50s in the nanomolar range. Antagonists as well as agonists were identified in an IL-6 dependent cell proliferation assay. Size minimization yielded peptides of about one-third the molecular mass of the original proteins, without significant loss of activities and, additionally, lead to the identification of the loops responsible for function.

Conclusions

This study shows 3F protein is amenable to introduce amino acid changes in the loops that enable preparation of a high diversity library that can be utilized to obtain ligands against macromolecules. We believe this is the first report of protein engineering to convert a neurotoxin to receptor ligands other than the parent receptor, the identification of an agonist from non-immunoglobulin proteins, the construction of peptide mimic of IL-6, and the successful size reduction of a single-chain protein.

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.

Structure, function and evolution of three-finger toxins: mini proteins with multiple targets

Snake venoms are complex mixtures of pharmacologically active peptides and proteins. These protein toxins belong to a small number of superfamilies of proteins. Three-finger toxins belong to a superfamily of non-enzymatic proteins found in all families of snakes. They have a common structure of three beta-stranded loops extending from a central core containing all four conserved disulphide bonds. Despite the common scaffold, they bind to different receptors/acceptors and exhibit a wide variety of biological effects. Thus, the structure-function relationships of this group of toxins are complicated and challenging. Studies have shown that the functional sites in these 'sibling' toxins are located on various segments of the molecular surface. Targeting to a wide variety of receptors and ion channels and hence distinct functions in this group of mini proteins is achieved through a combination of accelerated rate of exchange of segments as well as point mutations in exons. In this review, we describe the structural and functional diversity, structure-function relationships and evolution of this group of snake venom toxins.

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

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

Beta-cardiotoxin: a new three-finger toxin from Ophiophagus hannah (king cobra) venom with beta-blocker activity

Snake venoms have provided a number of novel ligands with therapeutic potential. We have constructed a partial cDNA library from the mRNA of Ophiophagus hannah (king cobra) venom gland tissue and identified five new genes encoding proteins belonging to the three-finger toxin family of snake venom proteins. We have isolated and characterized one of these beta-sheet containing proteins with a mass of 7012.43 +/- 0.91 Da from the venom. The protein was nonlethal up to a dose of 10 mg/kg when injected intraperitoneally into Swiss albino mice. However, it induces labored breathing and death at a dose of 100 mg/kg. It does not show any hemolytic or anticoagulant activity. It caused a dose-dependent decrease of heart rate in vivo (anesthetized Sprague-Dawley rats) and also ex vivo (Langendorff isolated rat heart). This is in contrast to classical cardiotoxins from snake venom that increase the heart rate in animals. Radioligand displacement studies showed that this protein targets beta-adrenergic receptors with a binding affinity (Ki) of 5.3 and 2.3 microM toward beta1 and beta2 subtypes, respectively, to bring about its effect, and hence, it was named as beta-cardiotoxin. This is the first report of a natural exogenous beta-blocker.