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

Label-free SERS characterization of snake venoms by exploring the cysteine environs with bone-shaped gold nanoparticles

Identification of snake venoms is a vital step in the treatment of fatal snakebites. In this study, we use the gold-thiolate interaction between a cysteine residue and gold nanoparticles to establish a SERS method for the differentiation of the venoms of Trimeresurus stejnegeri and Bungarus multicinctus. We confirm the preference of gold nanoparticles over silver for the SERS study of snake venoms by a binding experiment that also functions to differentiate the two venom samples by colorimetry and UV-vis spectroscopy. We report the SERS spectra of Trimeresurus stejnegeri and Bungarus multicinctus venoms for the first time. The spectra display distinct SERS signatures of the snake venoms on bone-shaped gold nanoparticles made with a house recipe. These signatures correlate to selected segments of the venom proteins due to the anchoring effect of the gold-cysteine bond. The method is quick as it accomplishes in situ isolation of the structure of interest to avoid tedious purification of the samples. The location of the interactive cysteine residue makes a novel characteristic of proteins in general.

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.

Computer modeling of binding of diverse weak toxins to nicotinic acetylcholine receptors

Weak toxins are the "three-fingered" snake venoms toxins grouped together by having an additional disulfide in the N-terminal loop I. In general, weak toxins have low toxicity, and biological targets have been identified for some of them only, recently by detecting the effects on the nicotinic acetylcholine receptors (nAChR). Here the methods of docking and molecular dynamics simulations are used for comparative modeling of the complexes between four weak toxins of known spatial structure (WTX, candoxin, bucandin, gamma-bungarotoxin) and nAChRs. WTX and candoxin are those toxins whose blocking of the neuronal alpha7- and muscle-type nAChR has been earlier shown in binding assays and electrophysiological experiments, while for the other two toxins no such activity has been reported. Only candoxin and WTX are found here to give stable solutions for the toxin-nAChR complexes. These toxins appear to approach the binding site similarly to short alpha-neurotoxins, but their final position resembles that of alpha-cobratoxin, a long alpha-neurotoxin, in the complex with the acetylcholine-binding protein. The final spatial structures of candoxin and WTX complexes with the alpha7 neuronal or muscle-type nAChR are very similar and do not provide immediate answer why candoxin has a much higher affinity than WTX, but both of them share a virtually irreversible mode of binding to one or both these nAChR subtypes. Possible explanation comes from docking and MD simulations which predict fast kinetics of candoxin association with nAChR, no gross changes in the toxin conformation (with smaller toxin flexibility on alpha7 nAChR), while slow WTX binding to nAChR is associated with slow irreversible rearrangement both of the tip of the toxin loop II and of the binding pocket residues locking finally the toxin molecule. Computer modeling showed that the additional disulfide in the loop I is not directly involved in receptor binding of WTX and candoxin, but it stabilizes the structure of loop I which plays an important role in toxin delivery to the binding site. In summary, computer modeling visualized possible modes of binding for those weak toxins which interact with the nAChR, provided no solutions for those weak toxins whose targets are not the nAChRs, and demonstrated that the additional disulfide in loop I cannot be a sound criteria for joining all weak toxins into one group; the conclusion about the diversity of weak toxins made from computer modeling is in accord with the earlier phylogenetic analysis.

Sequences, geographic variations and molecular phylogeny of venom phospholipases and threefinger toxins of eastern India Bungarus fasciatus and kinetic analyses of its Pro31 phospholipases A2

Eight phospholipases A2 (PLAs) and four three-finger-toxins (3FTx) from the pooled venom of Bungarus fasciatus (Bf) were previously studied and sequenced, but their expression pattern in individual Bf venom and possible geographic variations remained to be investigated. We herein analyze the individual venom of two Bf specimens from Kolkata (designated as KBf) to address this question. Seven PLAs and five 3FTx were purified from the KBf venoms, and respective cDNAs were cloned from venom glands of one of the snakes. Comparison of their mass and N-terminal sequence revealed that all the PLAs were conserved in both KBf venoms, but that two of their 3FTx isoforms were variable. When comparing the sequences of these KBf-PLAs with those published, only one was found to be identical to that of Bf Vb-2, and the other five were 94-98% identical to those of Bf II, III, Va, VI and XI-2, respectively. Notably, the most abundant PLA isoforms of Bf and KBf venoms contain Pro31 substitution. They were found to have abnormally low k(cat) values but high affinity for Ca2+. Phylogenetic analysis based on the sequences of venom group IA PLAs showed a close relationship between Bungarus and Australian and marine Elapidae. As the five deduced sequences of KBf-3FTx are only 62-82% identical to the corresponding Bf-3FTx from the pooled venom, the 3FTx apparently have higher degree of individual and geographic variations than the PLAs. None of the KBf-3FTx was found to be neurotoxic or very lethal; phylogenetic analyses of the 3FTx also revealed the unique evolution of Bf as compared with other kraits.