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

In vitro immunoreactivity and in vivo neutralization of Trimeresurus gracilis venom with antivenoms targeting four pit viper species

Snakebite envenomation is a significant global health issue that requires specific antivenom treatments. In Taiwan, available antivenoms target a variety of snakes, but none specifically target Trimeresurus gracilis, an endemic and protected species found in the high mountain areas of Taiwan. This study evaluated the effectiveness of existing antivenoms against T. gracilis venom, focusing on a bivalent antivenom developed for Trimeresurus stejnegeri and Protobothrops mucrosquamatus (TsPmAV), as well as monovalent antivenoms for Deinagkistrodon acutus (DaAV) and Gloydius brevicaudus (GbAV). Our research involved in vivo toxicity testing in mice and in vitro immunobinding experiments using (chaotropic) enzyme-linked immunosorbent assays, comparing venoms from four pit viper species (T. gracilis, T. stejnegeri, P. mucrosquamatus, and D. acutus) with three types of antivenoms. These findings indicate that TsPmAV partially neutralized T. gracilis venom, marginally surpassing the efficacy of DaAV. In vitro tests revealed that GbAV displayed higher binding capacities toward T. gracilis venom than TsPmAV or DaAV. Comparisons of electrophoretic profiles also reveal that T. gracilis venom has fewer snake venom C-type lectin like proteins than D. acutus, and has more P-I snake venom metalloproteases or fewer phospholipase A2 than G. brevicaudus, T. stejnegeri, or P. mucrosquamatus. This study highlights the need for antivenoms that specifically target T. gracilis, as current treatments using TsPmAV show limited effectiveness in neutralizing local effects in patients. These findings provide crucial insights into clinical treatment protocols and contribute to the understanding of the evolutionary adaptation of snake venom, aiding in the development of more effective antivenoms for human health.

Malaysian and Chinese King Cobra Venom Cytotoxicity in Melanoma and Neonatal Foreskin Fibroblasts Is Mediated by Age and Geography

Snake venoms constitute a complex, rapidly evolving trait, whose composition varies between and within populations depending on geographical location, age and preys (diets). These factors have determined the adaptive evolution for predatory success and link venom heterogeneity with prey specificity. Moreover, understanding the evolutionary drivers of animal venoms has streamlined the biodiscovery of venom-derived compounds as drug candidates in biomedicine and biotechnology. The king cobra (Ophiophagus hannah; Cantor, 1836) is distributed in diverse habitats, forming independent populations, which confer differing scale markings, including between hatchlings and adults. Furthermore, king cobra venoms possess unique cytotoxic properties that are used as a defensive trait, but their toxins may also have utility as promising anticancer-agent candidates. However, the impact of geographical distribution and age on these potential venom applications has been typically neglected. In this study, we hypothesised that ontogenetic venom variation accompanies the morphological distinction between hatchlings and adults. We used non-transformed neonatal foreskin (NFF) fibroblasts to examine and compare the variability of venom cytotoxicity between adult captive breeding pairs from Malaysian and Chinese lineages, along with that of their progeny upon hatching. In parallel, we assessed the anticancer potential of these venoms in human-melanoma-patient-derived cells (MM96L). We found that in a geographical distribution and gender-independent manner, venoms from hatchlings were significantly less cytotoxic than those from adults (NFF; ~Log EC50: 0.5-0.6 vs. 0.2-0.35 mg/mL). This is consistent with neonates occupying a semifossorial habitat, while adults inhabit more above-ground habitats and are therefore more conspicuous to potential predators. We also observed that Malaysian venoms exhibited a slightly higher cytotoxicity than those from the Chinese cobra cohorts (NFF; Log EC50: 0.1-0.3 vs. 0.3-0.4 mg/mL), which is consistent with Malaysian king cobras being more strongly aposematically marked. These variations are therefore suggestive of differential anti-predator strategies associated with the occupation of distinct niches. However, all cobra venoms were similarly cytotoxic in both melanoma cells and fibroblasts, limiting their potential medical applications in their native forms.

Hydrogen peroxide from L-amino acid oxidase of king cobra (Ophiophagus hannah) venom attenuates Pseudomonas biofilms

Because of the high incidence of Pseudomonas aeruginosa biofilms-related nosocomial infections, venoms from common Thai snakes were tested. Although venoms from king cobra (Ophiophagus hannah; OH) and green pit viper (Trimeresurus albolabris) showed the broadest antibacterial spectrum, OH venom demonstrated more profound anti-biofilm activities against P. aeruginosa. Additionally, purified L-amino acid oxidase from OH venom (OH-LAAO), using a three-step chromatography and protein identification, reduced biofilm mass as indicated by the downregulation of several genes, including the genes for biofilm synthesis (algD and pslB) and biofilm regulators (algU, gacA, and siaD). Moreover, OH-LAAO disrupted Pseudomonas-preformed biofilms via upregulation of several genes for biofilm dispersion (nbdA, bdlA, and dipA) and biofilm degradation (endA and pslG), resulting in a reduction of the biofilm biomass. Due to the antimicrobial effects and anti-biofilm activities (reduced production plus increased dispersion) neutralized by catalase, a hydrogen peroxide (H₂O₂)-degrading enzyme, the enhanced H₂O₂ by OH venom might be one of the anti-biofilm mechanisms. Hence, OH-LAAO was proposed as a novel agent against Pseudomonas biofilms for either treatment or prevention. More studies are interesting.

An immunodetection assay developed using cobra cytotoxin-specific antibodies: Potential diagnostics for cobra envenoming

Cobra (Naja spp.) envenoming is a life-threatening medical emergency, and a correct diagnosis is crucial to initiating timely and appropriate antivenom treatment. However, snakebite diagnostics remain unavailable in Southeast Asia. This study, therefore, developed an immunodetection assay with a potential diagnostic application for cobra envenoming. The cytotoxin of Naja kaouthia (Thai Monocled Cobra) (Nk-CTX) was purified from its venom to produce CTX-specific antibodies in rabbits and chickens. A double-antibody sandwich enzyme-linked immunosorbent assay was developed using the purified anti-Nk-CTX antibodies (immunoglobulin G and immunoglobulin Y), and its selectivity, specificity, and sensitivity for the venoms of five major cobra species in Southeast Asia (N. kaouthia, Naja sumatrana, Naja sputatrix, and Naja siamensis, Naja philippinensis) were studied. The results showed the immunoassay discriminates cobra venoms from other species commonly implicated in snakebites in Southeast Asia, i.e., the Malayan Krait, Many-banded Krait, King Cobra, Eastern Russell's Viper, Malayan Pit Viper and White-lipped Pit Viper. The immunoassay has a high sensitivity for the five cobra venoms, with detection limits (LoD) ranging from 0.6 to 2.6 ng/ml. Together, the findings suggest the potential diagnostic application of the cytotoxin immunoassay for cobra envenoming. The immunoassay was found to exhibit high immunoreactivity toward ten Asiatic cobra venoms (absorbance>1.5), in contrast to African cobra venoms with low immunoreactivity (absorbance<0.9). Considering the varying CTX antigenicity between Asiatic and African cobras, the immunoassay for African cobras should utilize antibodies produced specifically from the cytotoxins of African cobra venoms.

Snake Venomics: Fundamentals, Recent Updates, and a Look to the Next Decade

Venomic research, powered by techniques adapted from proteomics, transcriptomics, and genomics, seeks to unravel the diversity and complexity of venom through which knowledge can be applied in the treatment of envenoming, biodiscovery, and conservation. Snake venom proteomics is most extensively studied, but the methods varied widely, creating a massive amount of information which complicates data comparison and interpretation. Advancement in mass spectrometry technology, accompanied by growing databases and sophisticated bioinformatic tools, has overcome earlier limitations of protein identification. The progress, however, remains challenged by limited accessibility to samples, non-standardized quantitative methods, and biased interpretation of -omic data. Next-generation sequencing (NGS) technologies enable high-throughput venom-gland transcriptomics and genomics, complementing venom proteomics by providing deeper insights into the structural diversity, differential expression, regulation and functional interaction of the toxin genes. Venomic tissue sampling is, however, difficult due to strict regulations on wildlife use and transfer of biological materials in some countries. Limited resources for techniques and funding are among other pertinent issues that impede the progress of venomics, particularly in less developed regions and for neglected species. Genuine collaboration between international researchers, due recognition of regional experts by global organizations (e.g., WHO), and improved distribution of research support, should be embraced.

Addressing the global snakebite crisis with geo-spatial analyses - Recent advances and future direction

Venomous snakebite is a neglected tropical disease that annually leads to hundreds of thousands of deaths or long-term physical and mental ailments across the developing world. Insufficient data on spatial variation in snakebite risk, incidence, human vulnerability, and accessibility of medical treatment contribute substantially to ineffective on-ground management. There is an urgent need to collect data, fill knowledge gaps and address on-ground management problems. The use of novel, and transdisciplinary approaches that take advantage of recent advances in spatio-temporal models, 'big data', high performance computing, and fine-scale spatial information can add value to snakebite management by strategically improving our understanding and mitigation capacity of snakebite. We review the background and recent advances on the topic of snakebite related geospatial analyses and suggest avenues for priority research that will have practical on-ground applications for snakebite management and mitigation. These include streamlined, targeted data collection on snake distributions, snakebites, envenomings, venom composition, health infrastructure, and antivenom accessibility along with fine-scale models of spatio-temporal variation in snakebite risk and incidence, intraspecific venom variation, and environmental change modifying human exposure. These measures could improve and 'future-proof' antivenom production methods, antivenom distribution and stockpiling systems, and human-wildlife conflict management practices, while simultaneously feeding into research on venom evolution, snake taxonomy, ecology, biogeography, and conservation.

Assessment of snakebite burdens, clinical features of envenomation, and strategies to improve snakebite management in Vietnam

The sheer paucity of scientific documentation of herpetofauna in Vietnam and the rudimentary healthcare response to snakebite have stimulated this review. Over six decades of data culled from public data bases and search engines, have been used to assess snakebite burdens, clinical features of envenomation, and strategies for snakebite management in Vietnam. In addition, biochemical and proteomic analyses to decipher venom composition, rapid analytical techniques to be used for clinical diagnosis of snakebite in Vietnam have been discussed in detail. The assessment of efficacy, safety, and quality of commercial antivenom produced in Vietnam and improvement of antivenom production to meet the national requirement has been critically examined. It is apparent that snake bite incidence in Vietnam is exacerbated by mismatch in demand and supply of antivenom therapy, insufficient medical facilities, preference for traditional healers and poor management of clinical records. The impediments arising from geographical and species-specific variation in venom composition can be overcome by the 'Omics approach', and scientific documentation of pathophysiological manifestations post envenomation. The development of next generation of therapeutics, encouraging clinical research, novel approaches and social awareness against snakebite and its treatments have been suggested to significantly reduce the snakebite mortality and morbidity in this region.

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.

Dynamic changes of serum protein in rats with acute intoxication of Chinese cobra snake venom by proteomic analysis

To elucidate the toxic mechanism of snake venom at the protein level, proteomics technology was applied to investigate the effect of venom on circulation in the mammalian body. Temporal proteomic analysis was performed to profile the dynamic changes in the sera of Sprague–Dawley rats administered with Chinese cobra venom or saline. Using 8-plex iTRAQ analysis, 392 and 636 serum proteins were identified to be linearly upregulated or downregulated over time in the low-dose group and high-dose group, respectively. These proteins were mainly associated with the acute phase response pathway, complement system, and liver X receptor (LXR)/retinoid X receptor (RXR) and farnesoid X receptor (FXR)/RXR activation pathways. Compared with the low-dose group, the immune response and integrin pathways were inhibited in the high-dose group, although no obvious effect was observed. With consistently higher or lower expression in the high-dose group compared to the low-dose group throughout the whole process of venom poisoning, two proteins, Kininogen-1 (KNG1) and orosomucoid 1 (ORM1), which are involved in metabolism and immune response, occupied a core position in the pathway network and are considered venom dose-dependent biomarker candidates.

Suppression of cardiomyocyte functions by β-CTX isolated from the Thai king cobra (Ophiophagus hannah) venom via an alternative method

Background

Beta-cardiotoxin (β-CTX), the three-finger toxin isolated from king cobra (Ophiophagus hannah) venom, possesses β-blocker activity as indicated by its negative chronotropy and its binding property to both β-1 and β-2 adrenergic receptors and has been proposed as a novel β-blocker candidate. Previously, β-CTX was isolated and purified by FPLC. Here, we present an alternative method to purify this toxin. In addition, we tested its cytotoxicity against different mammalian muscle cell types and determined the impact on cardiac function in isolated cardiac myocyte so as to provide insights into the pharmacological action of this protein.

Methods

β-CTX was isolated from the crude venom of the Thai king cobra using reverse-phased and cation exchange HPLC. In vitro cellular viability MTT assays were performed on mouse myoblast (C2C12), rat smooth muscle (A7r5), and rat cardiac myoblast (H9c2) cells. Cell shortening and calcium transient dynamics were recorded on isolated rat cardiac myocytes over a range of β-CTX concentration.

Results

Purified β-CTX was recovered from crude venom (0.53% w/w). MTT assays revealed 50% cytotoxicity on A7r5 cells at 9.41 ± 1.14 µM (n = 3), but no cytotoxicity on C2C12 and H9c2 cells up to 114.09 µM. β-CTX suppressed the extend of rat cardiac cell shortening in a dose-dependent manner; the half-maximal inhibition concentration was 95.97 ± 50.10 nM (n = 3). In addition, the rates of cell shortening and re-lengthening were decreased in β-CTX treated myocytes concomitant with a prolongation of the intracellular calcium transient decay, indicating depression of cardiac contractility secondary to altered cardiac calcium homeostasis.

Conclusion

We present an alternative purification method for β-CTX from king cobra venom. We reveal cytotoxicity towards smooth muscle and depression of cardiac contractility by this protein. These data are useful to aid future development of pharmacological agents derived from β-CTX.

Interspecific and intraspecific venom enzymatic variation among cobras (Naja sp. and Ophiophagus hannah)

The genera Ophiophagus and Naja comprise part of a clade of snakes referred to as cobras, dangerously venomous front-fanged snakes in the family Elapidae responsible for significant human mortality and morbidity throughout Asia and Africa. We evaluated venom enzyme variation for eleven cobra species and three N. kaouthia populations using SDS-PAGE venom fingerprinting and numerous enzyme assays. Acetylcholinesterase and PLA₂ activities were the most variable between species, and PLA₂ activity was significantly different between Malaysian and Thailand N. kaouthia populations. Venom metalloproteinase activity was low and significantly different among most species, but levels were identical for N. kaouthia populations; minor variation in venom L-amino acid oxidase and phosphodiesterase activities were seen between cobra species. Naja siamensis venom lacked the α-fibrinogenolytic activity common to other cobra venoms. In addition, venom from N. siamensis had no detectable metalloproteinase activity and exhibited an SDS-PAGE profile with reduced abundance of higher mass proteins. Venom profiles from spitting cobras (N. siamensis, N. pallida, and N. mossambica) exhibited similar reductions in higher mass proteins, suggesting the evolution of venoms of reduced complexity and decreased enzymatic activity among spitting cobras. Generally, the venom proteomes of cobras show highly abundant three-finger toxin diversity, followed by large quantities of PLA₂s. However, PLA₂ bands and activity were very reduced for N. haje, N. annulifera and N. nivea. Venom compositionalenzy analysis provides insight into the evolution, diversification and distribution of different venom phenotypes that complements venomic data, and this information is critical for the development of effective antivenoms and snakebite treatment.

Geographical variations in king cobra (Ophiophagus hannah) venom from Thailand, Malaysia, Indonesia and China: On venom lethality, antivenom immunoreactivity and in vivo neutralization

The wide distribution of king cobra (Ophiophagus hannah), a medically important venomous snake in Asia could be associated with geographical variation in the toxicity and antigenicity of the venom. This study investigated the lethality of king cobra venoms (KCV) from four geographical locales (Malaysia, Thailand, Indonesia, China), and the immunological binding as well as in vivo neutralization activities of three antivenom products (Thai Ophiophagus hannah monovalent antivenom, OHMAV; Indonesian Serum Anti Bisa Ular, SABU; Chinese Naja atra monovalent antivenom, NAMAV) toward the venoms. The Indonesian and Chinese KCV were more lethal (median lethal dose, LD₅₀ ~0.5 μg/g) than those from Malaysia and Thailand (LD₅₀ ~1.0 μg/g). The antivenoms, composed of F(ab)'₂, were variably immunoreactive toward the KCV from all locales, with OHMAV exhibited the highest immunological binding activity. In mice, OHMAV neutralized the neurotoxic lethality of Thai KCV most effectively (normalized potency = 118 mg venom neutralized per g antivenom) followed by Malaysian, Indonesian and Chinese KCV. In comparison, the hetero-specific SABU was remarkably less potent by at least 6 to10 folds, whereas NAMAV appeared to be non-effective. The finding supports that a specific king cobra antivenom is needed for the effective treatment of king cobra envenomation in each region.

Naja atra venom-spit ophthalmia in Taiwan: An epidemiological survey from 1990 to 2016

BACKGROUND

Venomous snakebites are common during hot seasons in Taiwan. However, rarely is venom spat directly into the subject's eyes, causing eye injury. Despite being uncommon, analytical data regarding venom-spit ophthalmia in Taiwan have been lacking. This study thus aimed to conduct an epidemiological survey on Naja atra venom-spit ophthalmia in Taiwan during the past decades to improve future care of such patients.

METHODS

Registered records of cases with snake venom injuries at the Taiwan National Poison Control Center from 1990 to 2016 were retrospectively reviewed, enrolling those with records of cobra venom-spit eye injuries and excluding exotic species. Demographic data, clinical symptoms, snake species, ocular conditions, management, and outcomes were recorded and analyzed.

RESULTS

A total of 39 cases suffering from Naja atra venom-spit ophthalmia were enrolled. The overall incidence rate was 1.6 cases per million people from 1990 to 2016. Among the included cases, most were unilaterally injured (82%), male (95%), aged 18 to 59 years (90%), injured during catching (51%), and injured in spring and summer (92%). Ocular symptoms occurred in 90% of the cases, majority of which included ocular pain (90%) and redness (85%). Conjunctivitis (67%) and corneal injury (59%) of involved eyes were common. Immediate water irrigations were done in all cases, most of whom received further topical eye drops, including antibiotics, corticosteroids, and vasoconstrictors. Although topical or intravenous antivenoms were administered in 11 cases, no obviously superior outcome was observed. Most cases (77%) were symptom free after the acute stage.

CONCLUSION

Although Naja atra venom-spit ophthalmia in Taiwan is uncommon, the risk for transient ocular symptoms and corneal/conjunctival injury does exists. Prompt ocular irrigation and professional ophthalmic care after envenomation help prevent serious ocular sequelae. Moreover, superior outcomes were not achieved with the use of antivenom. Nonetheless, further studies are required to clarify the role antivenoms play in venom-spit ophthalmia.

A first step towards building spectral libraries as complementary tools for snake venom proteome/peptidome studies

Snake venoms are complex mixtures of a large number of distinct proteins and peptides with biological activity. Peptide spectral libraries are compilations of previously identified MS/MS spectra obtained from proteomics experiments. Here we present the generation and use of a Venom Peptidome and a Venom Proteome spectral library for the analysis of venom proteomes and peptidomes from distinct snake species.

Adaptive evolution of distinct prey-specific toxin genes in rear-fanged snake venom

Venom proteins evolve rapidly, and as a trophic adaptation are excellent models for predator-prey evolutionary studies. The key to a deeper understanding of venom evolution is an integrated approach, combining prey assays with analysis of venom gene expression and venom phenotype. Here, we use such an approach to study venom evolution in the Amazon puffing snake, Spilotes sulphureus, a generalist feeder. We identify two novel three-finger toxins: sulditoxin and sulmotoxin 1. These new toxins are not only two of the most abundant venom proteins, but are also functionally intriguing, displaying distinct prey-specific toxicities. Sulditoxin is highly toxic towards lizard prey, but is non-toxic towards mammalian prey, even at greater than 22-fold higher dosage. By contrast, sulmotoxin 1 exhibits the reverse trend. Furthermore, evolutionary analysis and structural modelling show highest sequence variability in the central loop of these proteins, probably driving taxon-specific toxicity. This is, to our knowledge, the first case in which a bimodal and contrasting pattern of toxicity has been shown for proteins in the venom of a single snake in relation to diet. Our study is an example of how toxin gene neofunctionalization can result in a venom system dominated by one protein superfamily and still exhibit flexibility in prey capture efficacy.

Analysis of the efficacy of Taiwanese freeze-dried neurotoxic antivenom against Naja kaouthia, Naja siamensis and Ophiophagus hannah through proteomics and animal model approaches

In Southeast Asia, envenoming resulting from cobra snakebites is an important public health issue in many regions, and antivenom therapy is the standard treatment for the snakebite. Because these cobras share a close evolutionary history, the amino acid sequences of major venom components in different snakes are very similar. Therefore, either monovalent or polyvalent antivenoms may offer paraspecific protection against envenomation of humans by several different snakes. In Taiwan, a bivalent antivenom—freeze-dried neurotoxic antivenom (FNAV)—against Bungarus multicinctus and Naja atra is available. However, whether this antivenom is also capable of neutralizing the venom of other species of snakes is not known. Here, to expand the clinical application of Taiwanese FNAV, we used an animal model to evaluate the neutralizing ability of FNAV against the venoms of three common snakes in Southeast Asia, including two ‘true’ cobras Naja kaouthia (Thailand) and Naja siamensis (Thailand), and the king cobra Ophiophagus hannah (Indonesia). We further applied mass spectrometry (MS)-based proteomic techniques to characterize venom proteomes and identify FNAV-recognizable antigens in the venoms of these Asian snakes. Neutralization assays in a mouse model showed that FNAV effectively neutralized the lethality of N. kaouthia and N. siamensis venoms, but not O. hannah venom. MS-based venom protein identification results further revealed that FNAV strongly recognized three-finger toxin and phospholipase A2, the major protein components of N. kaouthia and N. siamensis venoms. The characterization of venom proteomes and identification of FNAV-recognizable venom antigens may help researchers to further develop more effective antivenom designed to block the toxicity of dominant toxic proteins, with the ultimate goal of achieving broadly therapeutic effects against these cobra snakebites.

Cobra envenomation is a public health issue in Southeast Asia. Currently, antivenom therapy is the standard treatment for snakebite. However, antivenoms are not available in many rural countries and communities or have only limited effectiveness. Taiwan has wealth of experience in producing antivenoms, including the bivalent freeze-dried neurotoxic antivenom (FNAV), which is raised against venom proteins from Bungarus multicinctus and Naja atra. Our results showed that FNAV effectively neutralized the lethality of Naja kaouthia(Thailand) and Naja siamensis (Thailand) venoms, but not Ophiophagus hannah (Indonesia) venom, in an animal model. We further characterized the venom proteome profiles of the four cobras and identified three abundant proteins—neurotoxin, cytotoxin and phospholipase A2—in the venom of N. atra, N. kaouthia and N. siamensisas the major antigens recognized by FNAV. In contrast, we found that β-cardiotoxin and phospholipase A2, common toxin proteins in all king cobra venom samples, are weakly or not recognized by FNAV. Our data provide evidence suggesting the potential use of Taiwan’s FNAV to treat envenomation by other cobra species (N. kaouthia and N. siamensis) in Southeast Asia. Moreover, our findings support the previous recommendation and current experimental approach that major cobra toxins are used as antigens to generate more efficient antivenoms than those currently available.

How the Cobra Got Its Flesh-Eating Venom: Cytotoxicity as a Defensive Innovation and Its Co-Evolution with Hooding, Aposematic Marking, and Spitting

The cytotoxicity of the venom of 25 species of Old World elapid snake was tested and compared with the morphological and behavioural adaptations of hooding and spitting. We determined that, contrary to previous assumptions, the venoms of spitting species are not consistently more cytotoxic than those of closely related non-spitting species. While this correlation between spitting and non-spitting was found among African cobras, it was not present among Asian cobras. On the other hand, a consistent positive correlation was observed between cytotoxicity and utilisation of the defensive hooding display that cobras are famous for. Hooding and spitting are widely regarded as defensive adaptations, but it has hitherto been uncertain whether cytotoxicity serves a defensive purpose or is somehow useful in prey subjugation. The results of this study suggest that cytotoxicity evolved primarily as a defensive innovation and that it has co-evolved twice alongside hooding behavior: once in the Hemachatus + Naja and again independently in the king cobras (Ophiophagus). There was a significant increase of cytotoxicity in the Asian Naja linked to the evolution of bold aposematic hood markings, reinforcing the link between hooding and the evolution of defensive cytotoxic venoms. In parallel, lineages with increased cytotoxicity but lacking bold hood patterns evolved aposematic markers in the form of high contrast body banding. The results also indicate that, secondary to the evolution of venom rich in cytotoxins, spitting has evolved three times independently: once within the African Naja, once within the Asian Naja, and once in the Hemachatus genus. The evolution of cytotoxic venom thus appears to facilitate the evolution of defensive spitting behaviour. In contrast, a secondary loss of cytotoxicity and reduction of the hood occurred in the water cobra Naja annulata, which possesses streamlined neurotoxic venom similar to that of other aquatic elapid snakes (e.g., hydrophiine sea snakes). The results of this study make an important contribution to our growing understanding of the selection pressures shaping the evolution of snake venom and its constituent toxins. The data also aid in elucidating the relationship between these selection pressures and the medical impact of human snakebite in the developing world, as cytotoxic cobras cause considerable morbidity including loss-of-function injuries that result in economic and social burdens in the tropics of Asia and sub-Saharan Africa.

Mapping Proteoforms and Protein Complexes From King Cobra Venom Using Both Denaturing and Native Top-down Proteomics

Characterizing whole proteins by top-down proteomics avoids a step of inference encountered in the dominant bottom-up methodology when peptides are assembled computationally into proteins for identification. The direct interrogation of whole proteins and protein complexes from the venom of Ophiophagus hannah (king cobra) provides a sharply clarified view of toxin sequence variation, transit peptide cleavage sites and post-translational modifications (PTMs) likely critical for venom lethality. A tube-gel format for electrophoresis (called GELFrEE) and solution isoelectric focusing were used for protein fractionation prior to LC-MS/MS analysis resulting in 131 protein identifications (18 more than bottom-up) and a total of 184 proteoforms characterized from 14 protein toxin families. Operating both GELFrEE and mass spectrometry to preserve non-covalent interactions generated detailed information about two of the largest venom glycoprotein complexes: the homodimeric l-amino acid oxidase (∼130 kDa) and the multichain toxin cobra venom factor (∼147 kDa). The l-amino acid oxidase complex exhibited two clusters of multiproteoform complexes corresponding to the presence of 5 or 6 N-glycans moieties, each consistent with a distribution of N-acetyl hexosamines. Employing top-down proteomics in both native and denaturing modes provides unprecedented characterization of venom proteoforms and their complexes. A precise molecular inventory of venom proteins will propel the study of snake toxin variation and the targeted development of new antivenoms or other biotherapeutics.

An analysis of venom ontogeny and prey-specific toxicity in the Monocled Cobra (Naja kaouthia)

Venoms of snakes of the family Elapidae (cobras, kraits, mambas, and relatives) are predominantly composed of numerous phospholipases A2 (PLA2s) and three-finger toxins (3FTxs), some of which are lethal while others are not significantly toxic. Currently, the only identified prey-specific toxins are several nonconventional 3FTxs, and given the large diversity of 3FTxs within Monocled Cobra (Naja kaouthia) venom, it was hypothesized that several 3FTxs, previously found to be non-toxic or weakly toxic 3FTxs in murine models, could potentially be toxic towards non-murine prey. Additionally, it was hypothesized that ontogenetic dietary shifts will be correlated with observable changes in specific 3FTx isoform abundance. Adult and juvenile N. kaouthia venom composition was investigated using ion-exchange FPLC, 1D and 2D SDS-PAGE, mass spectrometry, and various enzymatic and LD50 assays. Alpha-cobratoxin (α-elapitoxin) was the only significantly toxic (LD50 < 1 μg/g) 3FTx found in N. kaouthia venom and was equally toxic toward both lizard and mouse models. The abundance and diversity of 3FTxs and most enzyme activities did not vary between adult and juvenile cobra venoms; however, total venom PLA2 activity and specific PLA2 isoforms did vary, with juveniles lacking several of the least acidic PLA2s, and these differences could have both biological (related to predation) and clinical (antivenom efficacy) implications. Nevertheless, the ubiquitous presence of α-cobratoxin in both adult and juvenile cobra venoms, with high toxicity toward both reptiles and mammals, represents a venom compositional strategy wherein a single potent toxin effectively immobilizes a variety of prey types encountered across life history stages.

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.

Immunological cross-reactivity and neutralization of the principal toxins of Naja sumatrana and related cobra venoms by a Thai polyvalent antivenom (Neuro Polyvalent Snake Antivenom)

The low potency of cobra antivenom has been an area of concern in immunotherapy for cobra envenomation. This study sought to investigate factors limiting the neutralizing potency of cobra antivenom, using a murine model. We examined the immunological reactivity and neutralizing potency of a Thai polyvalent antivenom against the principal toxins of Naja sumatrana (Equatorial spitting cobra) venom and two related Asiatic cobra venom α-neurotoxins. The antivenom possesses moderate neutralizing potency against phospholipases A2 (P, potency of 0.98mg/mL) and moderately weak neutralizing potency against long-chain α-neurotoxins (0.26-0.42mg/mL) but was only weakly effective in neutralizing the short-chain α-neurotoxins and cardiotoxins (0.05-0.08mg/mL). The poor neutralizing potency of the antivenom on the low molecular mass short-chain neurotoxins and cardiotoxins is presumably the main limiting factor of the efficacy of the cobra antivenom. Our results also showed that phospholipase A2, which exhibited the highest ELISA reactivity and avidity, was most effectively neutralized, whereas N. sumatrana short-chain neurotoxin, which exhibited the lowest ELISA reactivity and avidity, was least effectively neutralized by the antivenom. These observations suggest that low immunoreactivity (low ELISA reactivity and avidity) is one of the reasons for poor neutralization of the cobra venom low molecular mass toxins. Nevertheless, the overall results show that there is a lack of congruence between the immunological reactivity of the toxins toward antivenom and the effectiveness of toxin neutralization by the antivenom, indicating that there are other factors that also contribute to the weak neutralization capacity of the antivenom. Several suggestions have been put forward to overcome the low efficacy of the cobra antivenom. The use of a 'proper-mix' formulation of cobra venoms as immunogen, whereby the immunogen mixture used for hyperimmunization contains a mix of various types of α-neurotoxins and cardiotoxins in sufficient amount, may also help to improve the efficacy and broaden the neutralization spectrum of the antivenom.

Are ticks venomous animals?

INTRODUCTION

As an ecological adaptation venoms have evolved independently in several species of Metazoa. As haematophagous arthropods ticks are mainly considered as ectoparasites due to directly feeding on the skin of animal hosts. Ticks are of major importance since they serve as vectors for several diseases affecting humans and livestock animals. Ticks are rarely considered as venomous animals despite that tick saliva contains several protein families present in venomous taxa and that many Ixodida genera can induce paralysis and other types of toxicoses. Tick saliva was previously proposed as a special kind of venom since tick venom is used for blood feeding that counteracts host defense mechanisms. As a result, the present study provides evidence to reconsider the venomous properties of tick saliva.

RESULTS

Based on our extensive literature mining and in silico research, we demonstrate that ticks share several similarities with other venomous taxa. Many tick salivary protein families and their previously described functions are homologous to proteins found in scorpion, spider, snake, platypus and bee venoms. This infers that there is a structural and functional convergence between several molecular components in tick saliva and the venoms from other recognized venomous taxa. We also highlight the fact that the immune response against tick saliva and venoms (from recognized venomous taxa) are both dominated by an allergic immunity background. Furthermore, by comparing the major molecular components of human saliva, as an example of a non-venomous animal, with that of ticks we find evidence that ticks resemble more venomous than non-venomous animals. Finally, we introduce our considerations regarding the evolution of venoms in Arachnida.

CONCLUSIONS

Taking into account the composition of tick saliva, the venomous functions that ticks have while interacting with their hosts, and the distinguishable differences between human (non-venomous) and tick salivary proteins, we consider that ticks should be referred to as venomous ectoparasites.

Correlation between the glycan variations and defibrinogenating activities of acutobin and its recombinant glycoforms

Acutobin isolated from Deinagkistrodon acutus venom has been used to prevent or treat stroke in patients. This defibrinogenating serine protease is a 39 kDa glycoprotein containing terminal disialyl-capped N-glycans. After sialidase treatment, the enzyme showed similar catalytic activities toward chromogenic substrate, and cleaved the Aα chain of fibrinogen as efficiently as the native acutobin did. However, the level of fibrinogen degradation products in mice after i.p.-injection of desialylated-acutobin was significantly lower than the level after acutobin injection, suggesting that the disialyl moieties may improve or prolong the half-life of acutobin. Two recombinant enzymes with identical protein structures and similar amidolytic activities to those of native acutobin were expressed from HEK293T and SW1353 cells and designated as HKATB and SWATB, respectively. Mass spectrometric profiling showed that their glycans differed from those of acutobin. In contrast to acutobin, HKATB cleaved not only the Aα chain but also the Bβ and γ chains of human fibrinogens, while SWATB showed a reduced α-fibrinogenase activity. Non-denaturing deglycosylation of these proteases by peptide N-glycosidase F significantly reduced their fibrinogenolytic activities and thermal stabilities. The in vivo defibrinogenating effect of HKATB was inferior to that of acutobin in mice. Taken together, our results suggest that the conjugated glycans of acutobin are involved in its interaction with fibrinogen, and that the selection of cells optimally expressing efficient glycoforms and further glycosylation engineering are desirable before a recombinant product can replace the native enzyme for clinical use.

Pharmacokinetics of Naja sumatrana (equatorial spitting cobra) venom and its major toxins in experimentally envenomed rabbits

BACKGROUND

The optimization of snakebite management and the use of antivenom depend greatly on the knowledge of the venom's composition as well as its pharmacokinetics. To date, however, pharmacokinetic reports on cobra venoms and their toxins are still relatively limited. In the present study, we investigated the pharmacokinetics of Naja sumatrana (Equatorial spitting cobra) venom and its major toxins (phospholipase A2, neurotoxin and cardiotoxin), following intravenous and intramuscular administration into rabbits.

PRINCIPAL FINDINGS

The serum antigen concentration-time profile of the N. sumatrana venom and its major toxins injected intravenously fitted a two-compartment model of pharmacokinetics. The systemic clearance (91.3 ml/h), terminal phase half-life (13.6 h) and systemic bioavailability (41.9%) of N. sumatrana venom injected intramuscularly were similar to those of N. sputatrix venom determined in an earlier study. The venom neurotoxin and cardiotoxin reached their peak concentrations within 30 min following intramuscular injection, relatively faster than the phospholipase A2 and whole venom (Tmax=2 h and 1 h, respectively). Rapid absorption of the neurotoxin and cardiotoxin from the injection site into systemic circulation indicates fast onsets of action of these principal toxins that are responsible for the early systemic manifestation of envenoming. The more prominent role of the neurotoxin in N. sumatrana systemic envenoming is further supported by its significantly higher intramuscular bioavailability (Fi.m.=81.5%) compared to that of the phospholipase A2 (Fi.m.=68.6%) or cardiotoxin (Fi.m.=45.6%). The incomplete absorption of the phospholipase A2 and cardiotoxin may infer the toxins' affinities for tissues at the injection site and their pathological roles in local tissue damages through synergistic interactions.

CONCLUSION/SIGNIFICANCE

Our results suggest that the venom neurotoxin is absorbed very rapidly and has the highest bioavailability following intramuscular injection, supporting its role as the principal toxin in systemic envenoming.

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.

Two color morphs of the pelagic yellow-bellied sea snake, Pelamis platura, from different locations of Costa Rica: snake venomics, toxicity, and neutralization by antivenom

The yellow-bellied sea snake, Pelamis platura, is the most broadly distributed snake species. Despite being endowed with a highly lethal venom, a proteomic analysis of its toxin composition was unavailable. The venoms of specimens collected in Golfo de Papagayo and Golfo Dulce (Costa Rica), where two distinctive color morphs occur, were chromatographically compared. The latter inhabits a fjord-like gulf where the transit of oceanic sea snakes into and from the basin is restricted, thus possibly affecting gene flow. RP-HPLC evidenced a conserved venom protein profile in both populations, despite their divergent color phenotypes. Following a trend observed in other sea snakes, P. platura venom is relatively simple, being composed of proteins of the three-finger toxin (3FTx), phospholipase A2 (PLA2), cysteine-rich secretory protein (CRISP), 5'-nucleotidase, and metalloproteinase families. The first three groups represent 49.9%, 32.9%, and 9.1% of total venom protein, respectively. The most abundant component (~26%) is pelamitoxin (P62388), a short-chain 3FTx, followed by a major basic PLA2 (~20%) and a group of three isoforms of CRISPs (~9%). Whereas isolated pelamitoxin was highly lethal to mice, neither the PLA2 nor the CRISP fraction caused death. However, the PLA2 rapidly increased plasma creatine kinase activity after intramuscular injection, indicating its myotoxic action. Differing from myotoxic PLA2s of viperids, this PLA2 was not cytolytic to murine myogenic cells in vitro, suggesting possible differences in its mechanism of action. The median lethal dose (LD50) estimates for P. platura crude venom in mice and in three species of fishes did not differ significantly. The sea snake antivenom manufactured by CSL Ltd. (Australia), which uses Enhydrina schistosa as immunogen, cross-recognized the three major components of P. platura venom and, accordingly, neutralized the lethal activity of crude venom and pelamitoxin, therefore being of potential usefulness in the treatment of envenomations by this species.

BIOLOGICAL SIGNIFICANCE

Integrative analyses of animal venoms that combine the power of proteomics (venomics) with the characterization of their functional and immunological properties are significantly expanding knowledge on these remarkable bioweapons, both from a basic and a medical perspective. Costa Rica harbors a unique population of the yellow-bellied sea snake, Pelamis platura, that is restricted to a fjord-like gulf (Golfo Dulce). This population differs markedly from oceanic populations found elsewhere along the Pacific coast of this country, by presenting a patternless bright yellow coloration, instead of the typical bicolored or tricolored pattern of this species. It has been suggested that the dominance of this yellow-morph in Golfo Dulce might reflect gene flow restrictions, caused by the oceanographic conditions at this location. The present study demonstrates that the remarkable phenotypic variation between the two color morphs inhabiting Golfo Dulce and Golfo de Papagayo, respectively, is not associated with differences in the expression of venom components, as shown by their conserved RP-HPLC profiles. Proteomic analysis revealed the relatively simple toxin composition of P. platura venom, which contains three predominant types of proteins: three-finger toxins (protein abundance: 49.9%), phospholipases A2 (32.9%), and cysteine-rich secretory proteins (9.1%), together with few minor components. Further, the involvement of these most abundant proteins in the toxic effects of the venom, and their cross-recognition and neutralization by a sea snake antivenom produced against the venom of Enhydrina schistosa, were analyzed.

In-vitro neurotoxicity of two Malaysian krait species (Bungarus candidus and Bungarus fasciatus) venoms: neutralization by monovalent and polyvalent antivenoms from Thailand

Bungarus candidus and Bungarus fasciatus are two species of krait found in Southeast Asia. Envenoming by these snakes is often characterized by neurotoxicity and, without treatment, causes considerable morbidity and mortality. In this study, the in vitro neurotoxicity of each species, and the effectiveness of two monovalent antivenoms and a polyvalent antivenom, against the neurotoxic effects of the venoms, were examined in a skeletal muscle preparation. Both venoms caused concentration-dependent inhibition of indirect twitches, and attenuated responses to exogenous nicotinic receptor agonists, in the chick biventer preparation, with B. candidus venom being more potent than B. fasciatus venom. SDS-PAGE and western blot analysis indicated different profiles between the venoms. Despite these differences, most proteins bands were recognized by all three antivenoms. Antivenom, added prior to the venoms, attenuated the neurotoxic effect of the venoms. Interestingly, the respective monovalent antivenoms did not neutralize the effects of the venom from the other Bungarus species indicating a relative absence of cross-neutralization. Addition of a high concentration of polyvalent antivenom, at the t₉₀ time point after addition of venom, partially reversed the neurotoxicity of B. fasciatus venom but not B. candidus venom. The monovalent antivenoms had no significant effect when added at the t₉₀ time point. This study showed that B. candidus and B. fasciatus venoms display marked in vitro neurotoxicity in the chick biventer preparation and administration of antivenoms at high dose is necessary to prevent or reverse neurotoxicity.

Proteomic and biochemical analyses of short-tailed pit viper (Gloydius brevicaudus) venom: age-related variation and composition-activity correlation

We conducted an in-depth analysis of the proteomic and biochemical profiles of the venom of neonate and adult short-tailed pit vipers (Gloydius brevicaudus). Identified proteins were assigned to a few main toxin families. Disintegrin, phospholipase A2 (PLA2), serine proteinase, cysteine-rich secretory protein, C-type lectin-like protein, l-amino acid oxidase and snake venom metalloproteinase (SVMP) were detected in both venoms, while 5'-nucleotidase was detected only in the adult venom. SVMP was the predominant protein family in both venoms (neonate: 65.7%; adult: 64.4%), followed by PLA2 (neonate: 13.4%; adult: 25.0%). Antivenomic analysis revealed that commercial G. brevicaudus antivenom almost neutralized the chromatographic peaks with medium and high molecular masses in both venoms, but did not completely recognize peaks with low molecular mass. Toxicological and enzymatic activities show remarkable age-related variation in G. brevicaudus venom, probably resulting from variation in venom composition. Our data demonstrate age-related variation across venomics, antivenomics and biochemical profiles of G. brevicaudus venom, and have implications for the management of G. brevicaudus bites, including improving antivenom preparation by combining both venoms.

BIOLOGICAL SIGNIFICANCE

This study investigates the composition and biochemical activity of neonate and adult Gloydius brevicaudus venoms. We found remarkable age-related variation in venom biological activity, likely the result of variation in venom composition. Antivenomics analysis was used to explore difference in neonate and adult G. brevicaudus venoms. Our findings have implications for the diagnosis and clinical management of G. brevicaudus bites, and the design of venom mixtures that will increase the efficacy of commercial antivenom. This article is part of a Special Issue entitled: Proteomics of non-model organisms.