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

The structural and functional divergence of a neglected three-finger toxin subfamily in lethal elapids

Bungarus multicinctus is a widely distributed and medically important elapid snake that produces lethal neurotoxic venom. To study and enhance existing antivenom, we explore the complete repertoire of its toxin genes based on de novo chromosome-level assembly and multi-tissue transcriptome data. Comparative genomic analyses suggest that the three-finger toxin family (3FTX) may evolve through the neofunctionalization of flanking LY6E. A long-neglected 3FTX subfamily (i.e., MKA-3FTX) is also investigated. Only one MKA-3FTX gene, which evolves a different protein conformation, is under positive selection and actively transcribed in the venom gland, functioning as a major toxin effector together with MKT-3FTX subfamily homologs. Furthermore, this lethal snake may acquire self-resistance to its β-bungarotoxin via amino acid replacements on fast-evolving KCNA2. This study provides valuable resources for further evolutionary and structure-function studies of snake toxins, which are fundamental for the development of effective antivenoms and drug candidates.

Three-finger proteins from snakes and humans acting on nicotinic receptors: Old and new

The first toxin to give rise to the three-finger protein (TFP) family was α-bungarotoxin (α-Bgt) from Bungarus multicinctus krait venom. α-Bgt was crucial for research on nicotinic acetylcholine receptors (nAChRs), and in this Review article we focus on present data for snake venom TFPs and those of the Ly6/uPAR family from mammalians (membrane-bound Lynx1 and secreted SLURP-1) interacting with nAChRs. Recently isolated from Bungarus candidus venom, αδ-bungarotoxins differ from α-Bgt: they bind more reversibly and distinguish two binding sites in Torpedo californica nAChR. Naja kaouthia α-cobratoxin, classical blocker of nAChRs, was shown to inhibit certain GABA-A receptor subtypes, whereas α-cobratoxin dimer with 2 intermolecular disulfides has a novel type of 3D structure. Non-conventional toxin WTX has additional 5th disulfide not in the central loop, as α-Bgt, but in the N-terminal loop, like all Ly6/uPAR proteins, and inhibits α7 and Torpedo nAChRs. A water-soluble form of Lynx1, ws-Lynx1, was expressed in E. coli, its ¹ H-NMR structure and binding to several nAChRs determined. For SLURP-1, similar information was obtained with its recombinant analogue rSLURP-1. A common feature of ws-Lynx1, rSLURP-1, and WTX is their activity against nAChRs and muscarinic acetylcholine receptors. Synthetic SLURP-1, identical to the natural protein, demonstrated some differences from rSLURP-1 in distinguishing nAChR subtypes. The loop II fragment of the Lynx1 was synthesized having the same µM affinity for the Torpedo nAChR as ws-Lynx1. This review illustrates the productivity of parallel research of nAChR interactions with the two TFP groups.

Snake venom three-finger toxins and their potential in drug development targeting cardiovascular diseases

Cardiovascular diseases such as coronary and peripheral artery diseases, venous thrombosis, stroke, hypertension, and heart failure are enormous burden to health and economy globally. Snake venoms have been the sources of discovery of successful therapeutics targeting cardiovascular diseases. For example, the first-in-class angiotensin-converting enzyme inhibitor captopril was designed largely based on bradykinin-potentiating peptides from Bothrops jararaca venom. In the recent years, sensitive and high throughput approaches drive discovery and cataloging of new snake venom toxins. As one of the largest class of snake venom toxin, there are now>700 sequences of three-finger toxins (3FTxs) available, many of which are yet to be studied. While the function of 3FTxs are normally associated with neurotoxicity, increasingly more 3FTxs have been characterized to have pharmacological effects on cardiovascular systems. Here we focus on this family of snake venom toxins and their potential in developing therapeutics against cardiovascular diseases.

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

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

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.

Nature-Derived Peptides: A Growing Niche for GPCR Ligand Discovery

G protein-coupled receptors (GPCRs) represent important drug targets, as they regulate pivotal physiological processes and they have proved to be readily druggable. Natural products have been and continue to be amongst the most valuable sources for drug discovery and development. Here, we surveyed small molecules and (poly-)peptides derived from plants, animals, fungi, and bacteria, which modulate GPCR signaling. Among naturally occurring compounds, peptides from plants, cone-snails, snakes, spiders, scorpions, fungi, and bacteria are of particular interest as lead compounds for the development of GPCR ligands, since they cover a chemical space, which differs from that of synthetic small molecules. Peptides, however, face challenges, some of which can be overcome by studying plant-derived compounds. We argue here that the opportunities outweigh the challenges.

Proteomic analysis to unravel the complex venom proteome of eastern India Naja naja: Correlation of venom composition with its biochemical and pharmacological properties

The complex venom proteome of the eastern India (EI) spectacled cobra (Naja naja) was analyzed using tandem mass spectrometry of cation-exchange venom fractions. About 75% of EI N. naja venom proteins were <18kDa and cationic at physiological pH of blood. SDS-PAGE (non-reduced) analysis indicated that in the native state venom proteins either interacted with each-other or self-aggregated resulting in the formation of higher molecular mass complexes. Proteomic analysis revealed that 43 enzymatic and non-enzymatic proteins in EI N. naja venom with a percent composition of about 28.4% and 71.6% respectively were distributed over 15 venom protein families. The three finger toxins (63.8%) and phospholipase A₂s (11.4%) were the most abundant families of non-enzymatic and enzymatic proteins, respectively. nanoLC-ESI-MS/MS analysis demonstrated the occurrence of acetylcholinesterase, phosphodiesterase, cholinesterase and snake venom serine proteases in N. naja venom previously not detected by proteomic analysis. ATPase, ADPase, hyaluronidase, TAME, and BAEE-esterase activities were detected by biochemical analysis; however, due to a limitation in the protein database depository they were not identified in EI N. naja venom by proteomic analysis. The proteome composition of EI N. naja venom was well correlated with its in vitro and in vivo pharmacological properties in experimental animals and envenomed human.

BIOLOGICAL SIGNIFICANCE

Proteomic analysis reveals the complex and diverse protein profile of EI N. naja venom which collectively contributes to the severe pathophysiological manifestation upon cobra envenomation. The study has also aided in comprehending the compositional variation in venom proteins of N. naja within the Indian sub-continent. In addition, this study has also identified several enzymes in EI N. naja venom which were previously uncharacterized by proteomic analysis of Naja venom.

Nonconventional three-finger toxin BMLCL from krait Bungarus multicinctus venom with high affinity interacts with nicotinic acetylcholine receptors

Nonconventional three-finger toxin BMLCL was isolated from B. multicinctus venom, and its interaction with different subtypes of nicotinic acetylcholine receptor (nAChR) was studied. It was found that BMLCL is able to interact with high efficiency with both α7 and muscle type nAChRs.

Muscarinic toxins

Muscarinic toxins isolated from the venom of Dendroaspis snakes may interact with a high affinity, large selectivity and various functional properties with muscarinic receptors. Therefore, these toxins are invaluable tools for studying the physiological role, molecular functioning and structural organization of the five subtypes of these G-Protein Coupled Receptors. We review the data on the most relevant results dealing with the isolation/identification, mode of action, structure/function and exploitation of these toxins and finally highlight the unresolved issues related to their pharmacological studies.

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.

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.

Muscarinic toxins: tools for the study of the pharmacological and functional properties of muscarinic receptors

Muscarinic receptors mediate metabotropic actions of acetylcholine in the CNS and PNS and autocrine functions of acetylcholine in non-neuronal systems. Because of the lack of highly selective muscarinic ligands, the precise location, functional role, and roles in various diseases of the five muscarinic receptor subtypes remain unclear. Muscarinic toxins isolated from the venom of Dendroaspis snakes have a natural high affinity and selectivity, associated with roles as competitive antagonists, allosteric modulators, and potential agonists. These toxins may therefore be invaluable tools for studying muscarinic receptors. We review data on the structural and pharmacological characterization of the muscarinic toxins, focusing on recent structure-function studies on toxin-receptor interactions. We discuss the potential benefits of using these toxins for investigating muscarinic function in vivo.

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.

A catalog for transcripts in the venom gland of the Agkistrodon acutus: identification of the toxins potentially involved in coagulopathy

Agkistrodon acutus is a special agkistrodon halys, only distributed in Southern China, with a few exceptions in Vietnam. It is a cherished element used in traditional Chinese medicine. In order to produce a global panorama of gene expression in the Agkistrodon acutus venom gland, a non-normalized cDNA library was constructed, and 8696 high quality 5' end expressed sequenced tags (ESTs) were sequenced and analyzed. The initial sequences were assembled into 2855 clusters. Of these clusters, only 45.60% clusters matched known sequence and 54.40% had no match to any known sequence in GenBank. Except for putative cellular proteins (1184 clusters), the remaining 118 clusters (40.16% of all ESTs) corresponded to sequences associated with diverse toxin function. According to expression abundance, the major toxin components were metalloproteinases (32.08%) and C-type lectin (5.22%), and other components including bradykinin-potentiating peptide (0.90%), serine proteases (0.51%), nucleotidase and nuclease (0.41%), phospholipase A2 (0.30%), disintegrin (0.05%), cytokine-like molecules (0.06%), and other proteins (0.63%). The majority of these components are thought to be responsible for coagulopathy after A. acutus bites. We have therefore generated a comprehensive catalog of the A. acutus venom gland described so far. Gene expression from the very specialized secretory tissue, especially for those involved in coagulopathy, can be surveyed and provide important information in finding novel toxins.

Purification and characterization of Taiwan cobra venom proteins with weak toxicity

Two proteins G2a and G2b with molecular masses of approximately 24 kDa were isolated from Naja naja atra (Taiwan cobra) venom using sequential chromatography on gel filtration, ion-exchanger and reverse phase columns. The results of Edman degradation and mass analysis revealed that G2a is a cysteine-rich protein reported previously, and G2b is a novel polypeptide. CD spectra showed that the gross conformation of G2a and G2b notably differed. G2a exhibited an activity higher than that noted with G2b on inhibiting carbachol-induced muscle contraction. However, the two proteins weakly blocked muscle contraction evoked by K+. The observations that the two proteins exhibit the toxic activity in the concentration of micromolar range suggest that they are inherently weak toxins as other snake venom cysteine-rich proteins.

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.