Ivan Koludarov

Differential Evolution and Neofunctionalization of Snake Venom Metalloprotease Domains

Snake venom metalloproteases (SVMP) are composed of five domains: signal peptide, propeptide, metalloprotease, disintegrin, and cysteine-rich. Secreted toxins are typically combinatorial variations of the latter three domains. The SVMP-encoding genes of Psammophis mossambicus venom are unique in containing only the signal and propeptide domains. We show that the Psammophis SVMP propeptide evolves rapidly and is subject to a high degree of positive selection. Unlike Psammophis, some species of Echis express both the typical multidomain and the unusual monodomain (propeptide only) SVMP, with the result that a lower level of variation is exerted upon the latter. We showed that most mutations in the multidomain Echis SVMP occurred in the protease domain responsible for proteolytic and hemorrhagic activities. The cysteine-rich and disintegrin-like domains, which are putatively responsible for making the P-III SVMPs more potent than the P-I and P-II forms, accumulate the remaining variation. Thus, the binding sites on the molecules surface are evolving rapidly whereas the core remains relatively conserved. Bioassays conducted on two post-translationally cleaved novel proline-rich peptides from the P. mossambicus propeptide domain showed them to have been neofunctionalized for specific inhibition of mammalian a7 neuronal nicotinic acetylcholine receptors. We show that the proline rich postsynaptic specific neurotoxic peptides from Azemiops feae are the result of convergent evolution within the precursor region of the C-type natriuretic peptide instead of the SVMP. The results of this study reinforce the value of studying obscure venoms for biodiscovery of novel investigational ligands.

Evolution Stings: The Origin and Diversification of Scorpion Toxin Peptide Scaffolds

The episodic nature of natural selection and the accumulation of extreme sequence divergence in venom-encoding genes over long periods of evolutionary time can obscure the signature of positive Darwinian selection. Recognition of the true biocomplexity is further hampered by the limited taxon selection, with easy to obtain or medically important species typically being the subject of intense venom research, relative to the actual taxonomical diversity in nature. This holds true for scorpions, which are one of the most ancient terrestrial venomous animal lineages. The family Buthidae that includes all the medically significant species has been intensely investigated around the globe, while almost completely ignoring the remaining non-buthid families. Australian scorpion lineages, for instance, have been completely neglected, with only a single scorpion species (Urodacus yaschenkoi) having its venom transcriptome sequenced. Hence, the lack of venom composition and toxin sequence information from an entire continent’s worth of scorpions has impeded our understanding of the molecular evolution of scorpion venom. The molecular origin, phylogenetic relationships and evolutionary histories of most scorpion toxin scaffolds remain enigmatic. In this study, we have sequenced venom gland transcriptomes of a wide taxonomical diversity of scorpions from Australia, including buthid and non-buthid representatives. Using state-of-art molecular evolutionary analyses, we show that a majority of CSα/β toxin scaffolds have experienced episodic influence of positive selection, while most non-CSα/β linear toxins evolve under the extreme influence of negative selection. For the first time, we have unraveled the molecular origin of the major scorpion toxin scaffolds, such as scorpion venom single von Willebrand factor C-domain peptides (SV-SVC), inhibitor cystine knot (ICK), disulphide-directed beta-hairpin (DDH), bradykinin potentiating peptides (BPP), linear non-disulphide bridged peptides and antimicrobial peptides (AMP). We have thus demonstrated that even neglected lineages of scorpions are a rich pool of novel biochemical components, which have evolved over millions of years to target specific ion channels in prey animals, and as a result, possess tremendous implications in therapeutics.

Venom proteomic characterization and relative antivenom neutralization of two medically important Pakistani elapid snakes (Bungarus sindanus and Naja naja)

UNLABELLED Intra- and interspecific variation in venom composition has been shown to have a major effect upon the efficacy of antivenoms. Due to the absence of domestically produced antivenoms, Pakistan is wholly reliant upon antivenoms produced in other countries, such as India. However, the efficacy of these antivenoms in neutralising the venoms of Pakistani snakes has not been ascertained. This is symptomatic of the general state of toxicological research in this country, which has a myriad of highly toxic and medically important venomous animals. Thus, there is a dire need for knowledge regarding the fundamental proteomics of these venoms and applied knowledge of the relative efficacy of foreign antivenoms. Here we present the results of our proteomic research on two medically important snakes of Pakistan: Bungarus sindanus and Naja naja. Indian Polyvalent Antivenom (Bharat Serums and Vaccines Ltd), which is currently marketed for use in Pakistan, was completely ineffective against either Pakistani species. In addition to the expected pre- and post-synaptic neurotoxic activity, the venom of the Pakistan population of N. naja was shown to be quite divergent from other populations of this species in being potently myotoxic. These results highlight the importance of studying divergent species and isolated populations, where the same data not only elucidates clinical problems in need of immediate attention, but also uncovers sources for novel toxins with potentially useful activities. BIOLOGICAL SIGNIFICANCE Pakistan Bungarus sindanus and Naja naja venoms are differentially complex. Naja naja is potently myotoxic. Neither venom is neutralized by Indian antivenom. These results have direct implications for the treatment of envenomed patients in Pakistan. The unusually myotoxic effects of Naja naja demonstrates the value of studying remote populations for biodiscovery.

Venom Down Under: Dynamic Evolution of Australian Elapid Snake Toxins

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

Squeezers and leaf-cutters: differential diversification and degeneration of the venom system in toxicoferan reptiles

Although it has been established that all toxicoferan squamates share a common venomous ancestor, it has remained unclear whether the maxillary and mandibular venom glands are evolving on separate gene expression trajectories or if they remain under shared genetic control. We show that identical transcripts are simultaneously expressed not only in the mandibular and maxillary glands, but also in the enigmatic snake rictal gland. Toxin molecular frameworks recovered in this study were three-finger toxin (3FTx), CRiSP, crotamine (beta-defensin), cobra venom factor, cystatin, epididymal secretory protein, kunitz, l-amino acid oxidase, lectin, renin aspartate protease, veficolin, and vespryn. We also discovered a novel low-molecular weight disulfide bridged peptide class in pythonid snake glands. In the iguanian lizards, the most highly expressed are potentially antimicrobial in nature (crotamine (beta-defensin) and cystatin), with crotamine (beta-defensin) also the most diverse. However, a number of proteins characterized from anguimorph lizards and caenophidian snakes with hemotoxic or neurotoxic activities were recruited in the common toxicoferan ancestor and remain expressed, albeit in low levels, even in the iguanian lizards. In contrast, the henophidian snakes express 3FTx and lectin toxins as the dominant transcripts. Even in the constricting pythonid and boid snakes, where the glands are predominantly mucous-secreting, low-levels of toxin transcripts can be detected. Venom thus appears to play little role in feeding behavior of most iguanian lizards or the powerful constricting snakes, and the low levels of expression argue against a defensive role. However, clearly the incipient or secondarily atrophied venom systems of these taxa may be a source of novel compounds useful in drug design and discovery.

Firing the Sting: Chemically Induced Discharge of Cnidae Reveals Novel Proteins and Peptides from Box Jellyfish (Chironex fleckeri) Venom

Cnidarian venom research has lagged behind other toxinological fields due to technical difficulties in recovery of the complex venom from the microscopic nematocysts. Here we report a newly developed rapid, repeatable and cost effective technique of venom preparation, using ethanol to induce nematocyst discharge and to recover venom contents in one step. Our model species was the Australian box jellyfish (Chironex fleckeri), which has a notable impact on public health. By utilizing scanning electron microscopy and light microscopy, we examined nematocyst external morphology before and after ethanol treatment and verified nematocyst discharge. Further, to investigate nematocyst content or “venom” recovery, we utilized both top-down and bottom-up transcriptomics–proteomics approaches and compared the proteome profile of this new ethanol recovery based method to a previously reported high activity and recovery protocol, based upon density purified intact cnidae and pressure induced disruption. In addition to recovering previously characterized box jellyfish toxins, including CfTX-A/B and CfTX-1, we recovered putative metalloproteases and novel expression of a small serine protease inhibitor. This study not only reveals a much more complex toxin profile of Australian box jellyfish venom but also suggests that ethanol extraction method could augment future cnidarian venom proteomics research efforts.

Rapid Radiations and the Race to Redundancy: An Investigation of the Evolution of Australian Elapid Snake Venoms.

Australia is the stronghold of the front-fanged venomous snake family Elapidae. The Australasian elapid snake radiation, which includes approximately 100 terrestrial species in Australia, as well as Melanesian species and all the world’s true sea snakes, may be less than 12 million years old. The incredible phenotypic and ecological diversity of the clade is matched by considerable diversity in venom composition. The clade’s evolutionary youth and dynamic evolution should make it of particular interest to toxinologists, however, the majority of species, which are small, typically inoffensive, and seldom encountered by non-herpetologists, have been almost completely neglected by researchers. The present study investigates the venom composition of 28 species proteomically, revealing several interesting trends in venom composition, and reports, for the first time in elapid snakes, the existence of an ontogenetic shift in the venom composition and activity of brown snakes (Pseudonaja sp.). Trends in venom composition are compared to the snakes’ feeding ecology and the paper concludes with an extended discussion of the selection pressures shaping the evolution of snake venom.

Proteomic comparison of Hypnale hypnale (Hump-Nosed Pit-Viper) and Calloselasma rhodostoma (Malayan Pit-Viper) venoms

UNLABELLED Treatment of Hypnale hypnale bites with commercial antivenoms, even those raised against its sister taxon Calloselasma rhodostoma, has never been clinically successful. As these two genera have been separated for 20million years, we tested to see whether significant variations in venom had accumulated during this long period of evolutionary divergence, and thus could be responsible for the failure of antivenom. Proteomic analyses of C. rhodostoma and H. hypnale venom were performed using 1D and 2D PAGE as well as 2D-DIGE. C. rhodostoma venom was diverse containing large amounts of Disintegrin, Kallikrein, l-amino acid oxidase, Lectin, phospholipase A2 (acidic, basic and neutral) and Snake Venom Metalloprotease. In contrast, while H. hypnale also contained a wide range of toxin types, the venom was overwhelmingly dominated by two molecular weight forms of basic PLA2. 2D-DIGE (2-D Fluorescence Difference Gel Electrophoresis analysis) showed that even when a particular toxin class was shared between the two venoms, there were significant molecular weights or isoelectric point differences. This proteomic difference explains the past treatment failures with C. rhodostoma antivenom and highlights the need for a H. hypnale specific antivenom. BIOLOGICAL SIGNIFICANCE These results have direct implications for the treatment of envenomed patients in Sri Lanka. The unusual venom profile of Hypnale hypnale underscores the biodiscovery potential of novel snake venoms.

Correlation between ontogenetic dietary shifts and venom variation in Australian brown snakes ( Pseudonaja )

Venom is a key evolutionary trait, as evidenced by its widespread convergent evolution across the animal kingdom. In an escalating prey-predator arms race, venoms evolve rapidly to guarantee predatory or defensive success. Variation in venom composition is ubiquitous among snakes. Here, we tested variation in venom activity on substrates relevant to blood coagulation among Pseudonaja (brown snake) species, Australian elapids responsible for the majority of medically important human envenomations in Australia. A functional approach was employed to elucidate interspecific variation in venom activity in all nine currently recognised species of Pseudonaja. Fluorometric enzymatic activity assays were performed to test variation in whole venom procoagulant activity among species. Analyses confirmed the previously documented ontogenetic shift from non-coagulopathic venom in juveniles to coagulopathic venom as adults, except for the case of P. modesta, which retains non-coagulopathic venom as an adult. These shifts in venom activity correlate with documented ontogenetic shifts in diet among brown snakes from specialisation on reptilian prey as juveniles (and throughout the life cycle of P. modesta), to a more generalised diet in adults that includes mammals. The results of this study bring to light findings relevant to both clinical and evolutionary toxinology.

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.