Proteomic characterization and comparison of venoms from two elapid snakes (Bungarus multicinctus and Naja atra) from China

First Look at the Venoms of Two Sinomicrurus Snakes: Differences in Yield, Proteomic Profiles, and Immunorecognition by Commercial Antivenoms

Chinese coral snakes (Sinomicrurus) are highly neglected regarding their venom profiles and harm to humans, which impedes our ability to deeply understand their biological properties and explore their medicinal potential. In this study, we performed a comparative analysis to reveal the venom profiles of two Chinese coral snakes in terms of their venom yields, proteomic profiles, and immunorecognition by commercial antivenoms. The results showed that Sinomicrurus kelloggi expels more venom (lyophilized venom mass) than Sinomicrurus maccelellandi but possesses a similar solid venom content. These interspecific differences in venom yield were influenced by the snout-vent length. The venoms of these two species varied in their electrophoretic profiles, as well as in the presence or absence and relative abundance of protein families. They exhibited a 3-FTx-predominant phenotype, where the S. maccelellandi venom was dominated by 3-FTx (32.43%), SVMP (23.63%), PLA2 (19.88%), and SVSP (12.61%), while the S. kelloggi venom was dominated by 3-FTx (65.81%), LAAO (11.35%), and AMP (10.09%). While both the commercial Naja atra and Bungarus multicinctus antivenoms could immunorecognize these two Chinese coral snake venoms, the N. atra antivenom possessed a higher neutralization capability than the B. multicinctus antivenom for both species of coral snakes. Our findings show significant interspecific variations in the venom profiles of these Sinomicrurus snakes for the first time. We suggest screening or preparing specific antivenoms with high efficiency for the clinical treatment of envenomation caused by these snakes.

Snake venom weakens neurovascular integrity and promotes vulnerability to neuroinflammation

Snake envenomation poses significant medical challenges, particularly in subtropical and tropical regions, with long-term impacts on neurovascular integrity and neuroinflammation remaining underexplored. This study investigates the effects of venom from four species of venomous snakes in southern China-Zhoushan Cobra (Naja atra, NA), Many-banded Krait (Bungarus multicinctus, BM), Five-paced Pit Viper (Deinagkistrodon acutus, DA), and Chinese Moccasin (Protobothrops mucrosquamatus, PM) - on the blood-brain barrier (BBB) and chronic neuroinflammation. Using mass spectrometry, we analyzed venom protein compositions, while cytotoxic effects on mouse brain endothelial cells (bEND.3) were evaluated to determine IC50 values. In vitro BBB models and in vivo experiments in C57BL/6J mice revealed that NA venom, in particular, significantly compromised BBB integrity without inducing large-scale apoptosis, leading to persistent BBB disruption characterized by increased permeability and selective degradation of extracellular matrix and tight junction proteins. Moreover, to simulate secondary infections that often occur following snakebites, we combined venom exposure with lipopolysaccharide (LPS) treatment, which exacerbated neuroinflammatory responses by intensifying microglial activation and promoting a pro-inflammatory phenotype. These findings highlight the role of snake venom in compromising neurovascular integrity and promoting vulnerability to chronic neuroinflammation, emphasizing the need for further research into venom-induced neuroinflammatory pathways and their potential as therapeutic targets.

Proteomic Profiling of Venoms from Bungarus suzhenae and B. bungaroides: Enzymatic Activities and Toxicity Assessment

Kraits are venomous snakes of the genus Bungarus from the family Elapidae. Their venom typically demonstrates neurotoxicity; however, the toxicity is significantly influenced by the snake's species and geographical origin. Among the Bungarus species, Bungarus suzhenae and B. bungaroides have been poorly studied, with little to no information available regarding their venom composition. In this study, a proteomic approach was employed using LC-MS/MS to identify proteins from trypsin-digested peptides. The analysis revealed 102 venom-related proteins from 18 distinct functional protein families in the venom of B. suzhenae, with the primary components being three-finger toxins (3-FTx, 25.84%), phospholipase A2 (PLA2, 40.29%), L-amino acid oxidase (LAAO, 10.33%), Kunitz-type serine protease inhibitors (KUN, 9.48%), and snake venom metalloproteinases (SVMPs, 6.13%). In the venom of B. bungaroides, 99 proteins from 17 families were identified, with primary components being 3-FTx (33.87%), PLA2 (37.91%), LAAO (4.21%), and KUN (16.60%). Enzymatic activity assays confirmed the presence of key venom enzymes. Additionally, the LD50 values for B. suzhenae and B. bungaroides were 0.0133 μg/g and 0.752 μg/g, respectively, providing a reference for toxicity studies of these two species. This research elucidates the proteomic differences in the venoms of these two species, offering a foundation for developing antivenoms and clinical treatments for envenomation.

Deciphering toxico-proteomics of Asiatic medically significant venomous snake species: A systematic review and interactive data dashboard

Snakebite envenomation (SBE) is a neglected tropical disease (NTD) with an approximate 1.8 million cases annually. The tremendous figure is concerning, and the currently available treatment for snakebite envenomation is antivenom. However, the current antivenom has limited cross-neutralisation activity due to the variations in snake venom composition across species and geographical locations. The proteomics of medically important venomous species is essential as they study the venom compositions within and among different species. The advancement of sophisticated proteomic approaches allows intensive investigation of snake venoms. Nevertheless, there is a need to consolidate the venom proteomics profiles and distribution analysis to examine their variability patterns. This review systematically analysed the proteomics and toxicity profiles of medically important venomous species from Asia across different geographical locations. An interactive dashboard - Asiatic Proteomics Interactive Datasets was curated to consolidate the distribution patterns of the venom compositions, serve as a comprehensive directory for large-scale comparative meta-analyses. The population proteomics demonstrate higher diversities in the predominant venom toxins. Besides, inter-regional differences were also observed in Bungarus sp., Naja sp., Calliophis sp., and Ophiophagus hannah venoms. The elapid venoms are predominated with three-finger toxins (3FTXs) and phospholipase A2 (PLA2). Intra-regional variation is only significantly observed in Naja naja venoms. Proteomics diversity is more prominent in viper venoms, with widespread dominance observed in snake venom metalloproteinase (SVMP) and snake venom serine protease (SVSP). Correlations exist between the proteomics profiles and the toxicity (LD50) of the medically important venomous species. Additionally, the predominant toxins, alongside their pathophysiological effects, were highlighted and discussed as well. The insights of interactive toxico-proteomics datasets provide comprehensive frameworks of venom dynamics and contribute to developing antivenoms for snakebite envenomation. This could reduce misdiagnosis of SBE and accelerate the researchers' data mining process.

Three finger toxins of elapids: structure, function, clinical applications and its inhibitors

The WHO lists snakebite as a "neglected tropical disease". In tropical and subtropical areas, envenoming is an important public health issue. This review article describes the structure, function, chemical composition, natural inhibitors, and clinical applications of Elapids' Three Finger Toxins (3FTX) using scientific research data. The primary venomous substance belonging to Elapidae is 3FTX, that targets nAChR. Three parallel β-sheets combine to create 3FTX, which has four or five disulfide bonds. The three primary types of 3FTX are short-chain, long-chain, and nonconventional 3FTX. The functions of 3FTX depend on the specific toxin subtype and the target receptor or ion channel. The well-known effect of 3FTX is probably neurotoxicity because of the severe consequences of muscular paralysis and respiratory failure in snakebite victims. 3FTX have also been studied for their potential clinical applications. α-bungarotoxin has been used as a molecular probe to study the structure and function of nAChRs (Nicotinic Acetylcholine Receptors). Acid-sensing ion channel (ASIC) isoforms 1a and 1b are inhibited by Mambalgins, derived from Black mamba venom, which hinders their function and provide an analgesic effect. α- Cobra toxin is a neurotoxin purified from Chinese cobra (Naja atra) binds to nAChR at the neuronal junction and causes an analgesic effect for moderate to severe pain. Some of the plants and their compounds have been shown to inhibit the activity of 3FTX, and their mechanisms of action are discussed.

Aligning Post-Column ESI-MS, MALDI-MS, and Coagulation Bioassay Data of Naja spp., Ophiophagus hannah, and Pseudonaja textillis Venoms Chromatographically to Assess MALDI-MS and ESI-MS Complementarity with Correlation of Bioactive Toxins to Mass Spectrometric Data

Snakebite is a serious health issue in tropical and subtropical areas of the world and results in various pathologies, such as hemotoxicity, neurotoxicity, and local swelling, blistering, and tissue necrosis around the bite site. These pathologies may ultimately lead to permanent morbidity and may even be fatal. Understanding the chemical and biological properties of individual snake venom toxins is of great importance when developing a newer generation of safer and more effective snakebite treatments. Two main approaches to ionizing toxins prior to mass spectrometry (MS) analysis are electrospray ionization (ESI) and matrix-assisted laser desorption ionization (MALDI). In the present study, we investigated the use of both ESI-MS and MALDI-MS as complementary techniques for toxin characterization in venom research. We applied nanofractionation analytics to separate crude elapid venoms using reversed-phase liquid chromatography (RPLC) and high-resolution fractionation of the eluting toxins into 384-well plates, followed by online LC-ESI-MS measurements. To acquire clear comparisons between the two ionization approaches, offline MALDI-MS measurements were performed on the nanofractionated toxins. For comparison to the LC-ESI-MS data, we created so-called MALDI-MS chromatograms of each toxin. We also applied plasma coagulation assaying on 384-well plates with nanofractionated toxins to demonstrate parallel biochemical profiling within the workflow. The plotting of post-column acquired MALDI-MS data as so-called plotted MALDI-MS chromatograms to directly align the MALDI-MS data with ESI-MS extracted ion chromatograms allows the efficient correlation of intact mass toxin results from the two MS-based soft ionization approaches with coagulation bioassay chromatograms. This facilitates the efficient correlation of chromatographic bioassay peaks with the MS data. The correlated toxin masses from ESI-MS and/or MALDI-MS were all around 6-8 or 13-14 kDa, with one mass around 20 kDa. Between 24 and 67% of the toxins were observed with good intensity from both ionization methods, depending on the venom analyzed. All Naja venoms analyzed presented anticoagulation activity, whereas pro-coagulation was only observed for the Pseudonaja textillis venom. The data of MALDI-MS can provide complementary identification and characterization power for toxin research on elapid venoms next to ESI-MS.

Proteomic study of localized tissue necrosis by Naja atra venom

Naja atra bites often result in immediate and severe illness. The venom of N. atra contains a complex mixture of toxins that can cause significant damage to the patient's skin tissue. If left untreated, this condition can progress to localized necrosis, potentially resulting in impairment or even amputation in severe cases. Despite the known effects of the venom, the exact mechanisms underlying this tissue necrosis are not fully understood. This study aimed to investigate the protein components responsible for tissue necrosis induced by N. atra venom at both the organism-wide and molecular levels. To achieve this, venom was injected into Bama miniature pigs to cause ulcers, and exudate samples were collected at various time points after injection. Label-free proteomics analysis identified 1119, 1016, 938, 864, and 855 proteins in the exudate at 6, 12, 24, 36, and 48 h post-injection, respectively. Further analysis revealed 431 differentially expressed proteins, with S100A8, MMP-2, MIF, and IDH2 identified as proteins associated with local tissue necrosis. In this study, we established a Bama miniature pig model for N. atra venom injection and performed proteomic analysis of the wound exudate, which provides important insights into the molecular pathology of snakebite-induced tissue necrosis and potential theoretical bases for clinical treatment. Proteomic data from this study can be accessed through ProteomeXchange using the identifier PXD052498.

Specific Amino Acid Residues in the Three Loops of Snake Cytotoxins Determine Their Membrane Activity and Provide a Rationale for a New Classification of These Toxins

Cytotoxins (CTs) are three-finger membrane-active toxins present mainly in cobra venom. Our analysis of the available CT amino acid sequences, literature data on their membrane activity, and conformational equilibria in aqueous solution and detergent micelles allowed us to identify specific amino acid residues which interfere with CT incorporation into membranes. They include Pro9, Ser28, and Asn/Asp45 within the N-terminal, central, and C-terminal loops, respectively. There is a hierarchy in the effect of these residues on membrane activity: Pro9 > Ser28 > Asn/Asp45. Taking into account all the possible combinations of special residues, we propose to divide CTs into eight groups. Group 1 includes toxins containing all of the above residues. Their representatives demonstrated the lowest membrane activity. Group 8 combines CTs that lack these residues. For the toxins from this group, the greatest membrane activity was observed. We predict that when solely membrane activity determines the cytotoxic effects, the activity of CTs from a group with a higher number should exceed that of CTs from a group with a lower number. This classification is supported by the available data on the cytotoxicity and membranotropic properties of CTs. We hypothesize that the special amino acid residues within the loops of the CT molecule may indicate their involvement in the interaction with non-lipid targets.

The effect of Iranian snake, Naja naja oxiana venom on the blood glucose concentration and some biochemical parameters of experimental diabetic rats

Diabetes is a chronic disease resulting from impaired insulin production and function; leading to hyperglycaemia and long-term complications. The treatment for Type I diabetes treatment involves insulin injections while Type II diabetes treatments include drugs such as metformin and sulfonylureas, along with lifestyle changes. These medicines can be expensive and may have adverse effects. Therefore, the search for new therapeutic agents continues. Venoms from various animals yield numerous pharmacologically active compounds. In this study, we investigated the effects of the venom from an Iranian snake, Naja naja oxiana, on blood glucose concentration and certain serum biochemical parameters in male rats with induced diabetes. Diabetes was induced in male rats using either a single injection of streptozotocin (STZ) alone (55 mg/kg i. p.) or STZ (65 mg/kg i. p.) preceded by nicotinamide (230/kg i. p.) administered 15 min earlier. The diabetic rats produced by either method received a single injection of either vehicle or venom (0.2 or 0.4 mg/kg i. p.). In the STZ rats, this was done 13 days after diabetes induction, while in the STZ-nicotinamide rats, venom was injected 3 days after diabetes induction. The venom from Naja naja oxiana significantly reduced blood glucose levels in male rats with diabetes induced by either method. Additionally, the venom decreased serum cholesterol and triglycerides concentrations. However, the venom had no effect on the blood glucose levels of healthy male rats. Pretreatment with the venom did not prevent the induction of diabetes by STZ. These findings suggest that Naja naja oxiana venom exhibits an anti-diabetic effect and could be a potential candidate for effectively controlling diabetes.

Proteomic characteristics of six snake venoms from the Viperidae and Elapidae families in China and their relation to local tissue necrosis

Patients envenomed by snakes from the Viperidae and Elapidae families in China often have varying degrees of local tissue necrosis. Due to the relative clinical characteristics of local tissue necrosis and ulceration following envenoming, this study has analyzed the proteome of six snake venoms from the Viperidae and Elapidae family, and the toxin profiles of each snake were compared and correlated with the clinical manifestations that follow cytotoxic envenoming. Deinagkistrodon acutus and Naja atra envenomation induce severe ulceration, which is absent in Bungarus multicinctus envenomation and mild in the other three vipers. It is interesting to note that the proportion of c-type lectins (CTL) (20.63%) in Deinagkistrodon acutus venom was relatively high, which differs from the venom of other vipers. In addition, three-fingered toxin (3FTx) (2.15%) is present in the venom of Deinagkistrodon acutus, but has not been detected in the remaining three vipers. Snake venom metalloprotease (SVMP) (34.4%-44.7%), phospholipase A2 (PLA2) (9.81%-40.83%), and snake venom serine protease (SVSP) (9.44%-16.2%) represent the most abundant families of toxin in Viperidae venom. The Elapidae venom proteome was mainly composed of neurotoxins and cytotoxins, including 3FTx (39.28%-60.08%) and PLA2 (8.24%-58.95%) toxins, however, the proportion of CRISPS (26.36%) in Naja atra venom was relatively higher compared to Bungarus multicinctus venom. Significant differences in SVMP, SVSP, and 3FTx expression levels exist between the Viperidae and the Elapidae family. The main toxins responsible for the development of tissue necrosis and ulcerations following Viperidae envenoming are hematotoxins (SVSMP, SVSP) and myotoxins (PLA2). Deinagkistrodon acutus venom contains high levels of CTL and traces of 3FTx, leading to more severe local necrosis. However, Naja atra venom can also cause severe local necrosis through the effects of myotoxin (3FTx, CRISP, PLA2). Bungarus multicinctus venom does not contain myotoxins, resulting in pure systemic neurological manifestations no obvious necrosis of local tissue in patients.

Analytical Size Exclusion Chromatography Coupled with Mass Spectrometry in Parallel with High-Throughput Venomics and Bioassaying for Venom Profiling

Modern analytical size exclusion chromatography (SEC) is a suitable technique to separate venom toxin families according to their size characteristics. In this study, a method was developed to separate intact venom toxins from Bungarus multicinctus and Daboia russelii venoms via analytical SEC using volatile, non-salt-containing eluents for post-column mass spectrometry, coagulation bioassaying and high-throughput venomics. Two venoms were used to demonstrate the method developed. While the venom of Bungaurs multicinctus is known to exert anticoagulant effects on plasma, in this study, we showed the existence of both procoagulant toxins and anticoagulant toxins. For Daboia russelii venom, the method revealed characteristic procoagulant effects, with a 90 kDa mass toxin detected and matched with the Factor X-activating procoagulant heterotrimeric glycoprotein named RVV-X. The strong procoagulant effects for this toxin show that it was most likely eluted from size exclusion chromatography non-denatured. In conclusion, the separation of snake venom by size gave the opportunity to separate some specific toxin families from each other non-denatured, test these for functional bioactivities, detect the eluting mass on-line via mass spectrometry and identify the eluted toxins using high-throughput venomics.

Metabolomics analyses of serum metabolites perturbations associated with Naja atra bite

Naja atra bite is one of the most common severe snakebites in emergency departments. Unfortunately, the pathophysiological changes caused by Naja atra bite are unclear due to the lack of good animal models. In this study, an animal model of Naja atra bite in Guangxi Bama miniature pigs was established by intramuscular injection at 2 mg/kg of Naja atra venom, and serum metabolites were systematically analyzed using untargeted metabolomic and targeted metabolomic approaches. Untargeted metabolomic analysis revealed that 5045 chromatographic peaks were obtained in ESI+ and 3871 chromatographic peaks were obtained in ESI-. Screening in ESI+ modes and ESI- modes identified 22 and 36 differential metabolites compared to controls. The presence of 8 core metabolites of glutamine, arginine, proline, leucine, phenylalanine, inosine, thymidine and hippuric acid in the process of Naja atra bite was verified by targeted metabolomics significant difference (P<0.05). At the same time, during the verification process of the serum clinical samples with Naja atra bite, we found that the contents of three metabolites of proline, phenylalanine and inosine in the serum of the patients were significantly different from those of the normal human serum (P<0.05). By conducting functional analysis of core and metabolic pathway analysis, we revealed a potential correlation between changes in key metabolites after the Naja atra bite and the resulting pathophysiological alterations, and our research aims to establish a theoretical foundation for the prompt diagnosis and treatment of Naja atra bite.

Glutamine ameliorates Bungarus multicinctus venom-induced lung and heart injury through HSP70: NF-κB p65 and P53/PUMA signaling pathways involved

Background

Bungarus multicinctus is one of the most dangerous venomous snakes prone to cardiopulmonary damage with extremely high mortality. In our previous work, we found that glutamine (Gln) and glutamine synthetase (GS) in pig serum were significantly reduced after Bungarus multicinctus bite. In the present study, to explore whether there is a link between the pathogenesis of cardiopulmonary injury and Gln metabolic changes induced by Bungarus multicinctus venom. We investigated the effect of Gln supplementation on the lung and heart function after snakebite.

Methods

We supplemented different concentrations of Gln to mice that were envenomated by Bungarus multicinctus to observe the biological behavior, survival rate, hematological and pathological changes. Gln was supplemented immediately or one hour after the venom injection, and then changes in Gln metabolism were analyzed. Subsequently, to further explore the protective mechanism of glutamine on tissue damage, we measured the expression of heat-shock protein70 (HSP70), NF-κB P65, P53/PUMA by western blotting and real-time polymerase in the lung and heart.

Results

Gln supplementation delayed the envenoming symptoms, reduced mortality, and alleviated the histopathological changes in the heart and lung of mice bitten by Bungarus multicinctus. Additionally, Gln increased the activity of glutamine synthetase (GS), glutamate dehydrogenase (GDH) and glutaminase (GLS) in serum. It also balanced the transporter SLC7A11 expression in heart and lung tissues. Bungarus multicinctus venom induced the NF-κB nuclear translocation in the lung, while the HO-1 expression was suppressed. At the same time, venom activated the P53/PUMA signaling pathway and the BAX expression in the heart. Gln treatment reversed the above phenomenon and increased HSP70 expression.

Conclusion

Gln alleviated the glutamine metabolism disorder and cardiopulmonary damage caused by Bungarus multicinctus venom. It may protect lungs and heart against venom by promoting the expression of HSP70, inhibiting the activation of NF-κB and P53/PUMA, thereby delaying the process of snake venom and reducing mortality. The present results indicate that Gln could be a potential treatment for Bungarus multicinctus bite.

Development of sandwich ELISA and lateral flow assay for the detection of Bungarus multicinctus venom

Snakebite envenoming adversely affects human health and life worldwide. Presently, no suitable diagnostic tools for snakebite envenoming are available in China. Therefore, we sought to develop reliable diagnostic tests for snakebite management. We conducted affinity purification experiments to prepare species-specific antivenom antibody (SSAb). In brief, affinity chromatography with an antibody purification column (Protein A) was conducted to purify immunoglobulin G from Bungarus multicinctus (BM) venom hyperimmunized rabbit serum. The cross-reactive antibodies were removed from commercial BM antivenin by immune adsorption on the affinity chromatography columns of the other three venoms, Bungarus Fasciatus (FS), Naja atra (NA), and O. hannah (OH), generating SSAb. The results of western blot analysis and enzyme-linked immunosorbent assay (ELISA) showed the high specificity of the prepared SSAb. The obtained antibodies were then applied to ELISA and lateral flow assay (LFA) to detect BM venom. The resulting ELISA and LFA could specifically and rapidly detect BM venom in various samples with the limits of quantification as 0.1 and 1 ng/ml, respectively. This method could effectively detect snake venom in experimentally envenomed rats (simulating human envenomation), which could distinguish positive and negative samples within 10-15 min. This method also showed promise in serving as a highly useful tool for a rapid clinical distinguishing of BM bites and rational use of antivenom in emergency centers. The study also revealed cross-reactivity between BM and heterogenous venoms, suggesting that they shared common epitopes, which is of great significance for developing detection methods for venoms of the snakes belonging to the same family.

Effect of Several Naja atra Antivenom Injection Methods on the Rabbit Model of Naja naja atra Bite Poisoning

Snakebite is a global public health concern, which often occurs in tropical and subtropical underdeveloped areas, but it is often neglected. In the southern China, Naja naja atra (Chinese cobra) is a common venomous snake that causes swelling and necrosis of local tissues, even amputation and death. Currently, the main therapy is the administration of Naja atra antivenom, which greatly reduces mortality. However, the antivenom is not particularly effective in the improvement of local tissue necrosis. Clinically, antivenom is mainly administered intravenously. We speculated that the method of injection influences the efficacy of antivenom. In this study, the rabbit model was used to explore the effects of different antivenom injection methods on systemic and local poisoning symptoms. If topical injection of antivenom contributes to ameliorate tissue necrosis, then we need to reconsider the use of Naja atra antivenom.

Forensic identification of a fatal snakebite from Bungarus multicinctus (Chinese krait) by pathological and toxicological findings: a case report

Bungarus multicinctus (B. multicinctus) is one of the top ten venomous snakes in China, ranking first in lethality at 26.9-33.3%. However, to our knowledge, no forensic autopsy-related cases of death from B. multicinctus bite poisoning have been reported. There are surprisingly few reported cases of death from poisoning by other species of neurotoxic snakes. Neurotoxic snake venom is often highly toxic, and death can quickly occur when bitten in the wild if victims are not taken to a doctor in time. We presented a case of an adult female in Fujian Province of China who was bitten by a poisonous snake while digging for bamboo shoots in the mountains and died from the bite of B. multicinctus confirmed by enzyme-linked immunosorbent assays (ELISA) results. The autopsy's results, histopathological findings, and ELISA results reported here can be helpful for future forensic practice in B. multicinctus venom poisoning; we also briefly review the pathological changes of neurotoxin poisoning, which may be useful in other types of neurotoxin snake venom poisoning.

Current Insights in the Mechanisms of Cobra Venom Cytotoxins and Their Complexes in Inducing Toxicity: Implications in Antivenom Therapy

Cytotoxins (CTXs), an essential class of the non-enzymatic three-finger toxin family, are ubiquitously present in cobra venoms. These low-molecular-mass toxins, contributing to about 40 to 60% of the cobra venom proteome, play a significant role in cobra venom-induced toxicity, more prominently in dermonecrosis. Structurally, CTXs contain the conserved three-finger hydrophobic loops; however, they also exhibit a certain degree of structural diversity that dictates their biological activities. In their mechanism, CTXs mediate toxicity by affecting cell membrane structures and membrane-bound proteins and activating apoptotic and necrotic cell death pathways. Notably, some CTXs are also responsible for depolarizing neurons and heart muscle membranes, thereby contributing to the cardiac failure frequently observed in cobra-envenomed victims. Consequently, they are also known as cardiotoxins (CdTx). Studies have shown that cobra venom CTXs form cognate complexes with other components that potentiate the toxic effects of the venom's individual component. This review focuses on the pharmacological mechanism of cobra venom CTXs and their complexes, highlighting their significance in cobra venom-induced pathophysiology and toxicity. Furthermore, the potency of commercial antivenoms in reversing the adverse effects of cobra venom CTXs and their complexes in envenomed victims has also been discussed.

A Combined Bioassay and Nanofractionation Approach to Investigate the Anticoagulant Toxins of Mamba and Cobra Venoms and Their Inhibition by Varespladib

Envenomation by elapid snakes primarily results in neurotoxic symptoms and, consequently, are the primary focus of therapeutic research concerning such venoms. However, mounting evidence suggests these venoms can additionally cause coagulopathic symptoms, as demonstrated by some Asian elapids and African spitting cobras. This study sought to investigate the coagulopathic potential of venoms from medically important elapids of the genera Naja (true cobras), Hemachatus (rinkhals), and Dendroaspis (mambas). Crude venoms were bioassayed for coagulant effects using a plasma coagulation assay before RPLC/MS was used to separate and identify venom toxins in parallel with a nanofractionation module. Subsequently, coagulation bioassays were performed on the nanofractionated toxins, along with in-solution tryptic digestion and proteomics analysis. These experiments were then repeated on both crude venoms and on the nanofractionated venom toxins with the addition of either the phospholipase A2 (PLA2) inhibitor varespladib or the snake venom metalloproteinase (SVMP) inhibitor marimastat. Our results demonstrate that various African elapid venoms have an anticoagulant effect, and that this activity is significantly reduced for cobra venoms by the addition of varespladib, though this inhibitor had no effect against anticoagulation caused by mamba venoms. Marimastat showed limited capacity to reduce anticoagulation in elapids, affecting only N. haje and H. haemachatus venom at higher doses. Proteomic analysis of nanofractionated toxins revealed that the anticoagulant toxins in cobra venoms were both acidic and basic PLA2s, while the causative toxins in mamba venoms remain uncertain. This implies that while PLA2 inhibitors such as varespladib and metalloproteinase inhibitors such as marimastat are viable candidates for novel snakebite treatments, they are not likely to be effective against mamba envenomings.

Venom Variation of Neonate and Adult Chinese Cobras in Captivity Concerning Their Foraging Strategies

The venom and transcriptome profile of the captive Chinese cobra (Naja atra) is not characterized until now. Here, LC-MS/MS and illumine technology were used to unveil the venom and trascriptome of neonates and adults N. atra specimens. In captive Chinese cobra, 98 co-existing transcripts for venom-related proteins was contained. A total of 127 proteins belong to 21 protein families were found in the profile of venom. The main components of snake venom were three finger toxins (3-FTx), snake venom metalloproteinase (SVMP), cysteine-rich secretory protein (CRISP), cobra venom factor (CVF), and phosphodiesterase (PDE). During the ontogenesis of captive Chinese cobra, the rearrangement of snake venom composition occurred and with obscure gender difference. CVF, 3-FTx, PDE, phospholipase A2 (PLA2) in adults were more abundant than neonates, while SVMP and CRISP in the neonates was richer than the adults. Ontogenetic changes in the proteome of Chinese cobra venom reveals different strategies for handling prey. The levels of different types of toxin families were dramatically altered in the wild and captive specimens. Therefore, we speculate that the captive process could reshape the snake venom composition vigorously. The clear comprehension of the composition of Chinese cobra venom facilitates the understanding of the mechanism of snakebite intoxication and guides the preparation and administration of traditional antivenom and next-generation drugs for snakebite.

Differences between Two Groups of Burmese Vipers (Viperidae: Azemiops) in the Proteomic Profiles, Immunoreactivity and Biochemical Functions of Their Venoms

Two recently revised Azemiops snakes with apparent differences in their external appearances and skeletal morphologies but unclear genetic boundaries have been proposed. Some researchers have refrained from using the newly proposed taxonomy because these two “species” might be two clades corresponding to different geographical populations of Azemiops feae. To improve the understanding of the kinship of these two Burmese viper groups, more of their characteristics should be explored in depth. We performed a comparative analysis of the proteomic profiles and biochemical activities of snake venoms from these two groups (Sichuan A. feae and Zhejiang A. feae) and evaluated the immunorecognition capacity of commercial antivenoms toward them. Eight protein families were identified in venoms from these two groups, while phospholipase B was only detected in venom from Sichuan A. feae. These protein families displayed varying degrees of differences in relative abundance between venoms, and phospholipase A₂ (Sichuan A. feae: 57.15%; Zhejiang A. feae: 65.94%) was the predominated component. Gloydius brevicaudus antivenom exhibited the strongest capacity to immunologically recognize these two venoms, but this was mainly limited to components with high molecular masses, some of which differed between venoms. Additionally, Zhejiang A. feae venom was more toxic than Sichuan A. feae venom, and the venoms expressed remarkable differences in enzymatic activities, probably resulting from the variation in the relative abundance of specific protein families. Our findings unveil differences between the two Burmese viper groups in terms of proteomic profiles, immunoreactivity, and the biochemical functions of their venoms. This information will facilitate the management of snakebites caused by these snakes.

Bibliometric Analysis of Literature in Snake Venom-Related Research Worldwide (1933-2022)

Snake envenomation is a severe economic and health concern affecting countries worldwide. Snake venom carries a wide variety of small peptides and proteins with various immunological and pharmacological properties. A few key research areas related to snake venom, including its applications in treating cancer and eradicating antibiotic-resistant bacteria, have been gaining significant attention in recent years. The goal of the current study was to analyze the global profile of literature in snake venom research. This study presents a bibliometric review of snake venom-related research documents indexed in the Scopus database between 1933 and 2022. The overall number of documents published on a global scale was 2999, with an average annual production of 34 documents. Brazil produced the highest number of documents (n = 729), followed by the United States (n = 548), Australia (n = 240), and Costa Rica (n = 235). Since 1963, the number of publications has been steadily increasing globally. At a worldwide level, antivenom, proteomics, and transcriptomics are growing hot issues for research in this field. The current research provides a unique overview of snake venom research at global level from 1933 through 2022, and it may be beneficial in guiding future research.

Equatorial Spitting Cobra (Naja sumatrana) from Malaysia (Negeri Sembilan and Penang), Southern Thailand, and Sumatra: Comparative Venom Proteomics, Immunoreactivity and Cross-Neutralization by Antivenom

The Equatorial Spitting Cobra (Naja sumatrana) is a medically important venomous snake species in Southeast Asia. Its wide geographical distribution implies potential intra-specific venom variation, while there is no species-specific antivenom available to treat its envenoming. Applying a protein-decomplexing proteomic approach, the study showed that three-finger toxins (3FTX), followed by phospholipases A₂ (PLA₂), were the major proteins well-conserved across N. sumatrana venoms of different locales. Variations were noted in the subtypes and relative abundances of venom proteins. Of note, alpha-neurotoxins (belonging to 3FTX) are the least in the Penang specimen (Ns-PG, 5.41% of total venom proteins), compared with geographical specimens from Negeri Sembilan (Ns-NS, 14.84%), southern Thailand (Ns-TH, 16.05%) and Sumatra (Ns-SU, 10.81%). The alpha-neurotoxin abundance, in general, correlates with the venom’s lethal potency. The Thai Naja kaouthia Monovalent Antivenom (NkMAV) was found to be immunoreactive toward the N. sumatrana venoms and is capable of cross-neutralizing N. sumatrana venom lethality to varying degrees (potency = 0.49–0.92 mg/mL, interpreted as the amount of venom completely neutralized per milliliter of antivenom). The potency was lowest against NS-SU venom, implying variable antigenicity of its lethal alpha-neurotoxins. Together, the findings suggest the para-specific and geographical utility of NkMAV as treatment for N. sumatrana envenoming in Southeast Asia.

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.

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.

On characterizing the Red-headed Krait (Bungarus flaviceps) venom: Decomplexation proteomics, immunoreactivity and toxicity cross-neutralization by hetero-specific antivenoms

The Red-headed Krait (Bungarus flaviceps) is a medically important venomous snake species in Southeast Asia, while there is no specific antivenom available for its envenoming. This study investigated the venom composition through a decomplexation proteomic approach, and examined the immunoreactivity as well as cross-neutralization efficacy of two hetero-specific krait antivenoms, Bungarus candidus Monovalent Antivenom (BcMAV) and Bungarus fasciatus Monovalent Antivenom (BfMAV), against the venom of B. flaviceps from Peninsular Malaysia. A total of 43 non-redundant proteoforms belonging to 10 toxin families were identified in the venom proteome, which is dominated by phospholipases A₂ including beta-bungarotoxin lethal subunit (56.20 % of total venom proteins), Kunitz-type serine protease inhibitors (19.40 %), metalloproteinases (12.85 %) and three-finger toxins (7.73 %). The proteome varied in quantitative aspect from the earlier reported Indonesian (Sumatran) sample, suggesting geographical venom variation. BcMAV and BfMAV were immunoreactive toward the B. flaviceps venom, with BcMAV being more efficacious in immunological binding. Both antivenoms cross-neutralized the venom lethality with varying efficacy, where BcMAV was more potent than BfMAV by ~13 times (normalized potency: 38.04 mg/g vs. 2.73 mg/g, defined as the venom amount completely neutralized by one-gram antivenom protein), supporting the potential utility of BcMAV for para-specific neutralization against B. flaviceps venom.

The myth of cobra venom cytotoxin: More than just direct cytolytic actions

Cobra venom cytotoxin (CTX) is a non-enzymatic three-finger toxin that constitutes 40-60% of cobra venom. Thus, it plays an important role in the pathophysiology of cobra envenomation, especially in local dermonecrosis. The three-finger hydrophobic loops of CTX determine the cytotoxicity. Nevertheless, the actual mechanisms of cytotoxicity are not fully elucidated as they involve not only cytolytic actions but also intracellular signalling-mediated cell death pathways. Furthermore, the possible transition cell death pattern remains to be explored. The actual molecular mechanisms require further studies to unveil the relationship between different CTXs from different cobra species and cell types which may result in differential cell death patterns. Here, we discuss the biophysical interaction of CTX with the cell membrane involving four binding modes: electrostatic interaction, hydrophobic partitioning, isotropic phase, and oligomerisation. Oligomerisation of CTX causes pore formation in the membrane lipid bilayer. Additionally, the CTX-induced apoptotic pathway can be executed via death receptor-mediated extrinsic pathways and mitochondrial-mediated intrinsic pathways. We also discuss lysosomal-mediated necrosis and the occurrence of necroptosis following CTX action. Collectively, we provided an insight into concentration-dependent transition of cell death pattern which involves different mechanistic actions. This contributes a new direction for further investigation of cytotoxic pathways activated by the CTXs for future development of biotherapeutics targeting pathological effects caused by CTX.

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.

Analysis of the Composition of Deinagkistrodon acutus Snake Venom Based on Proteomics, and Its Antithrombotic Activity and Toxicity Studies

There is a strong correlation between the composition of Deinagkistrodon acutus venom proteins and their potential pharmacological effects. The proteomic analysis revealed 103 proteins identified through label-free proteomics from 30 different snake venom families. Phospholipase A2 (30.0%), snaclec (21.0%), antithrombin (17.8%), thrombin (8.1%) and metalloproteinases (4.2%) were the most abundant proteins. The main toxicity of Deinagkistrodon acutus venom is hematotoxicity and neurotoxicity, and it acts on the lung. Deinagkistrodon acutus venom may have anticoagulant and antithrombotic effects. In summary, the protein profile and related toxicity and pharmacological activity of Deinagkistrodon acutus venom from southwest China were put forward for the first time. In addition, we revealed the relationship between the main toxicity, pharmacological effects, and the protein components of snake venom.

The development of surgical risk score and evaluation of necrotizing soft tissue infection in 161 Naja atra envenomed patients

BACKGROUND

Naja atra bites cause wound necrosis, secondary infection, and necrotizing soft tissue infection (NSTI) requiring repetitive surgeries. Little information is known about the predictors for surgery after these bites.

MATERIALS AND METHODS

We retrospectively evaluated 161 patients envenomed by N. atra, 80 of whom underwent surgery because of wound necrosis and infection. We compared the patients' variables between surgical and non-surgical groups. To construct a surgical risk score, we converted the regression coefficients of the significant factors in the multivariate logistic regression into integers. We also examined the deep tissue cultures and pathological findings of the debrided tissue.

RESULTS

A lower limb as the bite site, a ≥3 swelling grade, bullae or blister formation, gastrointestinal (GI) effects, and fever were significantly associated with surgery in the multivariate logistic regression analysis. The surgical risk scores for these variables were 1, 1, 2, 1, and 2, respectively. At a ≥3-point cutoff value, the model has 71.8% sensitivity and 88.5% specificity for predicting surgery, with an area under the receiver operating characteristic curve of 0.88. The histopathological examinations of the debrided tissues supported the diagnosis of snakebite-induced NSTI. Twelve bacterial species were isolated during the initial surgery and eleven during subsequent surgeries.

DISCUSSION AND CONCLUSIONS

From the clinical perspective, swelling, bullae or blister formation, GI effects, and fever appeared quickly after the bite and before surgery. The predictive value of these factors for surgery was acceptable, with a ≥3-point risk score. The common laboratory parameters did not always predict the outcomes of N. atra bites without proper wound examination. Our study supported the diagnosis of NSTI and demonstrated the changes in bacteriology during the surgeries, which can have therapeutic implications for N. atra bites.

Snake Venom Proteomics of Samar Cobra (Naja samarensis) from the Southern Philippines: Short Alpha-Neurotoxins as the Dominant Lethal Component Weakly Cross-Neutralized by the Philippine Cobra Antivenom

The Samar Cobra, Naja samarensis, is endemic to the southern Philippines and is a WHO-listed Category 1 venomous snake species of medical importance. Envenomation caused by N. samarensis results in neurotoxicity, while there is no species-specific antivenom available for its treatment. The composition and neutralization of N. samarensis venom remain largely unknown to date. This study thus aimed to investigate the venom proteome of N. samarensis for a comprehensive profiling of the venom composition, and to examine the immunorecognition as well as neutralization of its toxins by a hetero-specific antivenom. Applying C₁₈ reverse-phase high-performance liquid chromatography (RP-HPLC) and tandem mass spectrometry (LC-MS/MS), three-finger toxins (3FTx) were shown to dominate the venom proteome by 90.48% of total venom proteins. Other proteins in the venom comprised snake venom metalloproteinases, phospholipases A_2, cysteine-rich secretory proteins, venom nerve growth factors, L-amino acid oxidases and vespryn, which were present at much lower abundances. Among all, short-chain alpha-neurotoxins (SαNTX) were the most highly expressed toxin within 3FTx family, constituting 65.87% of the total venom proteins. The SαNTX is the sole neurotoxic component of the venom and has an intravenous median lethal dose (LD₅₀) of 0.18 μg/g in mice. The high abundance and low LD₅₀ support the potent lethal activity of N. samarensis venom. The hetero-specific antivenom, Philippine Cobra Antivenom (PCAV, raised against Naja philippinensis) were immunoreactive toward the venom and its protein fractions, including the principal SαNTX. In efficacy study, PCAV was able to cross-neutralize the lethality of SαNTX albeit the effect was weak with a low potency of 0.20 mg/ml (defined as the amount of toxin completely neutralized per milliliter of the antivenom). With a volume of 5 ml, each vial of PCAV may cross-neutralize approximately 1 mg of the toxin in vivo. The findings support the potential para-specific use of PCAV in treating envenomation caused by N. samarensis while underscoring the need to improve the potency of its neutralization activity, especially against the highly lethal alpha-neurotoxins.

Compositional and toxicological investigation of pooled venom from farm-raised Naja atra

Background:

Naja atra is a venomous snake species medically relevant in China. In the current study, we evaluated the composition and toxicological profile of venom collected from farm-raised N. atra.

Methods:

Venom was collected from third-generation captive bred N. atra on a snake farm in Hunan Province, China. The venom was analyzed using sodium dodecyl sulfate polyacrylamide gel electrophoresis and nano-liquid chromatography with electrospray ionization tandem mass spectrometry. In addition, hemolytic activity, median lethal dose, serum biochemical and histopathological parameters were accessed.

Results:

N. atra venom proteome was dominated by phospholipase A₂ (46.5%) and three-finger toxins (41.4 %), and a set of common low relative abundance proteins, including cysteine-rich secretory proteins (4.7%), NGF-beta (2.4%), snake venom metalloproteinase (1.5%), glutathione peroxidase (0.6%), vespryn (0.3%), and 5ʹ-nucleotidases (0.2%) were also found. Furthermore, the venom exhibited direct hemolytic activity, neurotoxicity, myotoxicity, and high lethal potency in mice, with a subcutaneous median lethal dose of 1.02 mg/kg. Histopathological analysis and serum biochemical tests revealed that venom caused acute hepatic, pulmonary and renal injury in mice.

Conclusion:

This study revealed the composition and toxicity of venom collected from farm-raised N. atra, thereby providing a reference for the analysis of venom samples collected from captive-born venomous snakes in the future.

A strategy for the enrichment and characterization of disulfide bond-contained proteins from Chinese cobra (Naja atra) venom

A new strategy combined gold-coated magnetic nanocomposites assisted enrichment with mass spectrometry was developed for the characterization of disulfide bond-contained proteins from Chinese cobra (Naja atra) venom. In this work, core-shell nanocomposites were synthesized by the seed-mediated growth method and used for the enrichment of snake venom proteins containing disulfide bonds. A total of 3545 tryptic digested peptides derived from 96 venom proteins in Naja atra venom were identified. The venom proteins comprised 14 toxin families including three-finger toxins, phospholipase A₂ , snake venom metalloproteinase, cobra venom factor, and so forth. Extra 16 venom proteins were detected exclusively in the nanocomposites set, among which 11 venom proteins were from the three-finger toxins family. In the present study, the proposed simple and efficient protocol replaced the tedious and laborious technologies commonly used for pre-separating crude snake venom, suggesting widely implementation in low-abundance or trace disulfide bond-contained proteins or peptides characterization.

Proteomics and neutralization of Bungarus multicinctus (Many-banded Krait) venom: Intra-specific comparisons between specimens from China and Taiwan

The Many-banded Krait (Bungarus multicinctus) is a medically important venomous snake in East Asia. This study investigated the venom proteomes of B. multicinctus from Guangdong, southern China (BM-China) and insular Taiwan (BM-Taiwan), and the neutralization activities of two antivenom products (produced separately in China and Taiwan) against the lethal effect of the venoms. The venom proteomes of both specimens contained similar toxin families, notwithstanding small variations in the subtypes and abundances of minor components. More than 90% of the total venom proteins belong to three-finger toxins (3FTx, including alpha-neurotoxins) and phospholipases A₂ (PLA₂, including beta-bungarotoxins), supporting their key involvement in the pathophysiology of krait envenomation which manifests as pre- and post-synaptic neurotoxicity. The venoms exhibited potent neurotoxic and lethal effects with extremely low i.v. LD₅₀ of 0.027 μg/g (Bm-China) and 0.087 μg/g (Bm-Taiwan), respectively, in mice. Bungarus multicinctus monovalent antivenom (BMMAV) produced in China and Neuro bivalent antivenom (NBAV) produced in Taiwan were immunoreactive toward both venoms and their toxin fractions. The antivenoms neutralized the venom lethality variably, with BMMAV being more efficacious than NBAV by approximately two-fold. Findings suggest that the monovalent antivenom has a higher potency presumably due to its species-specificity toward the krait venom.

Elucidating the Venom Diversity in Sri Lankan Spectacled Cobra (Naja naja) through De Novo Venom Gland Transcriptomics, Venom Proteomics and Toxicity Neutralization

Inadequate effectiveness of Indian antivenoms in treating envenomation caused by the Spectacled Cobra/Indian Cobra (Naja naja) in Sri Lanka has been attributed to geographical variations in the venom composition. This study investigated the de novo venom-gland transcriptomics and venom proteomics of the Sri Lankan N. naja (NN-SL) to elucidate its toxin gene diversity and venom variability. The neutralization efficacy of a commonly used Indian antivenom product in Sri Lanka was examined against the lethality induced by NN-SL venom in mice. The transcriptomic study revealed high expression of 22 toxin genes families in NN-SL, constituting 46.55% of total transcript abundance. Three-finger toxins (3FTX) were the most diversely and abundantly expressed (87.54% of toxin gene expression), consistent with the dominance of 3FTX in the venom proteome (72.19% of total venom proteins). The 3FTX were predominantly S-type cytotoxins/cardiotoxins (CTX) and α-neurotoxins of long-chain or short-chain subtypes (α-NTX). CTX and α-NTX are implicated in local tissue necrosis and fatal neuromuscular paralysis, respectively, in envenomation caused by NN-SL. Intra-species variations in the toxin gene sequences and expression levels were apparent between NN-SL and other geographical specimens of N. naja, suggesting potential antigenic diversity that impacts antivenom effectiveness. This was demonstrated by limited potency (0.74 mg venom/ml antivenom) of the Indian polyvalent antivenom (VPAV) in neutralizing the NN-SL venom. A pan-regional antivenom with improved efficacy to treat N. naja envenomation is needed.

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

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

Venom of the Annulated Sea Snake Hydrophis cyanocinctus: A Biochemically Simple but Genetically Complex Weapon

Given that the venom system in sea snakes has a role in enhancing their secondary adaption to the marine environment, it follows that elucidating the diversity and function of venom toxins will help to understand the adaptive radiation of sea snakes. We performed proteomic and de novo NGS analyses to explore the diversity of venom toxins in the annulated sea snake (Hydrophis cyanocinctus) and estimated the adaptive molecular evolution of the toxin-coding unigenes and the toxicity of the major components. We found three-finger toxins (3-FTxs), phospholipase A2 (PLA2) and cysteine-rich secretory protein (CRISP) in the venom proteome and 59 toxin-coding unigenes belonging to 24 protein families in the venom-gland transcriptome; 3-FTx and PLA2 were the most abundant families. Nearly half of the toxin-coding unigenes had undergone positive selection. The short- (i.p. 0.09 μg/g) and long-chain neurotoxin (i.p. 0.14 μg/g) presented fairly high toxicity, whereas both basic and acidic PLA2s expressed low toxicity. The toxicity of H. cyanocinctus venom was largely determined by the 3-FTxs. Our data show the venom is used by H. cyanocinctus as a biochemically simple but genetically complex weapon and venom evolution in H. cyanocinctus is presumably driven by natural selection to deal with fast-moving prey and enemies in the marine environment.

Serum metabolomics of Bama miniature pigs bitten by Bungarus multicinctus

Bungarus multicinctus is one of the top ten venomous snakes in China, and its bite causes acute and severe diseases, but its pathophysiology remains poorly elucidated. Thus, an animal model of Bungarus multicinctus bite was established by intramuscular injection of 30μg/kg of Bungarus multicinctus venom, and then the serum metabolites were subsequently screened, identified and validated by ultra-performance liquid chromatography-quadrupole-time of flight-mass spectrometry (UPLC-Q-TOF-MS) methods to explore the potential biomakers and possible metabolic pathways. Untargeted metabolomics analysis showed that 36 and 38 endogenous metabolites levels changed in ESI+ and ESI-, respectively, KEGG pathway analysis showed that 5 metabolic pathways, including mineral absorption, central carbon metabolism in cancer, protein digestion and absorption, aminoacyl-tRNA biosynthesis and ABC transporters might be closely related to Bungarus multicinctus bite. Targeted metabolomics analysis showed that there were significant differences in serum D-proline, L-leucine and L-glutamine after Bungarus multicinctus bite (P < 0.05). In addition, receiver operating characteristic (ROC) analysis showed that the diagnostic efficiency of L-Glutamine was superior to other potential biomarkers and the AUC value was 0.944. Moreover, we found evidence for differences in the pathophysiology of glutamine between Bungarus multicinctus bite group and normal group, specifically with the content of glutamine synthetase (GS) and glutaminase (GLS). Taken together, the current study has successfully established an animal model of Bungarus multicinctus bite, and further identified the links between the metabolic perturbations and the pathophysiology and the potential diagnostic biomakers of Bungarus multicinctus bite, which provided valuable insights for studying the mechanism of Bungarus multicinctus bite.

Venom-gland transcriptomic, venomic, and antivenomic profiles of the spine-bellied sea snake (Hydrophis curtus) from the South China Sea

Background

A comprehensive evaluation of the -omic profiles of venom is important for understanding the potential function and evolution of snake venom. Here, we conducted an integrated multi-omics-analysis to unveil the venom-transcriptomic and venomic profiles in a same group of spine-bellied sea snakes (Hydrophis curtus) from the South China Sea, where the snake is a widespread species and might generate regionally-specific venom potentially harmful to human activities. The capacity of two heterologous antivenoms to immunocapture the H. curtus venom was determined for an in-depth evaluation of their rationality in treatment of H. curtus envenomation. In addition, a phylogenetic analysis by maximum likelihood was used to detect the adaptive molecular evolution of full-length toxin-coding unigenes.

Results

A total of 90,909,384 pairs of clean reads were generated via Illumina sequencing from a pooled cDNA library of six specimens, and yielding 148,121 unigenes through de novo assembly. Sequence similarity searching harvested 63,845 valid annotations, including 63,789 non-toxin-coding and 56 toxin-coding unigenes belonging to 22 protein families. Three protein families, three-finger toxins (3-FTx), phospholipase A₂ (PLA₂), and cysteine-rich secretory protein, were detected in the venom proteome. 3-FTx (27.15% in the transcriptome/41.94% in the proteome) and PLA₂ (59.71%/49.36%) were identified as the most abundant families in the venom-gland transcriptome and venom proteome. In addition, 24 unigenes from 11 protein families were shown to have experienced positive selection in their evolutionary history, whereas four were relatively conserved throughout evolution. Commercial Naja atra antivenom exhibited a stronger capacity than Bungarus multicinctus antivenom to immunocapture H. curtus venom components, especially short neurotoxins, with the capacity of both antivenoms to immunocapture short neurotoxins being weaker than that for PLA₂s.

Conclusions

Our study clarified the venom-gland transcriptomic and venomic profiles along with the within-group divergence of a H. curtus population from the South China Sea. Adaptive evolution of most venom components driven by natural selection appeared to occur rapidly during evolutionary history. Notably, the utility of commercial N. atra and B. multicinctus antivenoms against H. curtus toxins was not comprehensive; thus, the development of species-specific antivenom is urgently needed.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07824-7.

Convergent evolution of pain-inducing defensive venom components in spitting cobras

Convergent evolution provides insights into the selective drivers underlying evolutionary change. Snake venoms, with a direct genetic basis and clearly defined functional phenotype, provide a model system for exploring the repeated evolution of adaptations. While snakes use venom primarily for predation, and venom composition often reflects diet specificity, three lineages of cobras have independently evolved the ability to spit venom at adversaries. Using gene, protein, and functional analyses, we show that the three spitting lineages possess venoms characterized by an up-regulation of phospholipase A2 (PLA2) toxins, which potentiate the action of preexisting venom cytotoxins to activate mammalian sensory neurons and cause enhanced pain. These repeated independent changes provide a fascinating example of convergent evolution across multiple phenotypic levels driven by selection for defense.

In Vitro Neurotoxicity of Chinese Krait (Bungarus multicinctus) Venom and Neutralization by Antivenoms

Bungarus multicinctus, the Chinese krait, is a highly venomous elapid snake which causes considerable morbidity and mortality in southern China. B. multicinctus venom contains pre-synaptic PLA₂ neurotoxins (i.e., β-bungarotoxins) and post-synaptic neurotoxins (i.e., α-bungarotoxins). We examined the in vitro neurotoxicity of B. multicinctus venom, and the efficacy of specific monovalent Chinese B. multicinctus antivenom, and Australian polyvalent elapid snake antivenom, against venom-induced neurotoxicity. B. multicinctus venom (1-10 μg/mL) abolished indirect twitches in the chick biventer cervicis nerve-muscle preparation as well as attenuating contractile responses to exogenous ACh and CCh, but not KCl. This indicates a post-synaptic neurotoxic action but myotoxicity was not evident. Given that post-synaptic α-neurotoxins have a more rapid onset than pre-synaptic neurotoxins, the activity of the latter in the whole venom will be masked. The prior addition of Chinese B. multicinctus antivenom (12 U/mL) or Australian polyvalent snake antivenom (15 U/mL), markedly attenuated the neurotoxic actions of B. multicinctus venom (3 μg/mL) and prevented the inhibition of contractile responses to ACh and CCh. The addition of B. multicinctus antivenom (60 U/mL), or Australian polyvalent snake antivenom (50 U/mL), at the t₉₀ time point after the addition of B. multicinctus venom (3 μg/mL), did not restore the twitch height over 180 min. The earlier addition of B. multicinctus antivenom (60 U/mL), at the t₂₀ or t₅₀ time points, also failed to prevent the neurotoxic effects of the venom but did delay the time to abolish twitches based on a comparison of t₉₀ values. Repeated washing of the preparation with physiological salt solution, commencing at the t₂₀ time point, failed to reverse the neurotoxic effects of venom or delay the time to abolish twitches. This study showed that B. multicinctus venom displays marked in vitro neurotoxicity in a skeletal muscle preparation which is not reversed by antivenom. This does not appear to be related to antivenom efficacy, but due to the irreversible/pseudo-irreversible nature of the neurotoxins.

Immunoreactivity and neutralization study of Chinese Bungarus multicinctus antivenin and lab-prepared anti-bungarotoxin antisera towards purified bungarotoxins and snake venoms

Bungarus multicinctus is the most venomous snake distributed in China and neighboring countries of Myanmar, Laos, north Vietnam and Thailand. The high mortality rate of B. multicinctus envenomation is attributed to the lethal components of α-, β-, γ- and κ- bungarotoxins contained in the venom. Although anti-B. multicinctus sera were produced in Shanghai, Taiwan and Vietnam, the most widely clinic used product was term as B. multicinctus antivenin and manufactured by Shanghai Serum Bio-technology Co. Ltd. In the present investigation, high purity α-, β- and γ-bungarotoxins were separately isolated from B. multicinctus crude venom. Rabbit anti- α-, β- and γ-bungarotoxin antisera were prepared by common methods, respectively. LD50 values of α-, β- and γ-bungarotoxins were systematically determined via three administration pathways (intraperitoneal, intramuscular and intravenous injections) in Kunming mice. LD50 values of β-bungarotoxin were closely related with injection routines but those of both α- and γ-bungarotoxins were not dependent on the injection routines. Commercial B. multicinctus antivenin showed strong immunoreaction with high molecular weight fractions of the B. multicinctus but weakly recognized low molecular weight fractions like α- and γ-bungarotoxins. Although B. multicinctus antivenin showed immunoreaction with high molecular weight fractions of Bungarus fasciatus, Naja atra, Ophiophagus hannah venoms but the antivenin only demonstrated animal protection efficacy against O. hannah venom. These results indicated that the high molecular weight fractions of the O. hannah played an important role in venom lethality but those of B. fasciatus and N. atra did not have such a role.

Proteomic Investigations of Two Pakistani Naja Snake Venoms Species Unravel the Venom Complexity, Posttranslational Modifications, and Presence of Extracellular Vesicles

Latest advancement of omics technologies allows in-depth characterization of venom compositions. In the present work we present a proteomic study of two snake venoms of the genus Naja i.e., Naja naja (black cobra) and Naja oxiana (brown cobra) of Pakistani origin. The present study has shown that these snake venoms consist of a highly diversified proteome. Furthermore, the data also revealed variation among closely related species. High throughput mass spectrometric analysis of the venom proteome allowed to identify for the N. naja venom 34 protein families and for the N. oxiana 24 protein families. The comparative evaluation of the two venoms showed that N. naja consists of a more complex venom proteome than N. oxiana venom. Analysis also showed N-terminal acetylation (N-ace) of a few proteins in both venoms. To the best of our knowledge, this is the first study revealing this posttranslational modification in snake venom. N-ace can shed light on the mechanism of regulation of venom proteins inside the venom gland. Furthermore, our data showed the presence of other body proteins, e.g., ankyrin repeats, leucine repeats, zinc finger, cobra serum albumin, transferrin, insulin, deoxyribonuclease-2-alpha, and other regulatory proteins in these venoms. Interestingly, our data identified Ras-GTpase type of proteins, which indicate the presence of extracellular vesicles in the venom. The data can support the production of distinct and specific anti-venoms and also allow a better understanding of the envenomation and mechanism of distribution of toxins. Data are available via ProteomeXchange with identifier PXD018726.

An Examination of the Neutralization of In Vitro Toxicity of Chinese Cobra (Naja atra) Venom by Different Antivenoms

The Chinese Cobra (Naja atra) is an elapid snake of major medical importance in southern China. We describe the in vitro neurotoxic, myotoxic, and cytotoxic effects of N. atra venom, as well as examining the efficacy of three Chinese monovalent antivenoms (N. atra antivenom, Gloydius brevicaudus antivenom and Deinagkistrodon acutus antivenom) and an Australian polyvalent snake antivenom. In the chick biventer cervicis nerve-muscle preparation, N. atra venom (1–10 µg/mL) abolished indirect twitches in a concentration-dependent manner, as well as abolishing contractile responses to exogenous acetylcholine chloride (ACh) and carbamylcholine chloride (CCh), indicative of post-synaptic neurotoxicity. Contractile responses to potassium chloride (KCl) were also significantly inhibited by venom indicating myotoxicity. The prior addition of Chinese N. atra antivenom (0.75 U/mL) or Australian polyvalent snake antivenom (3 U/mL), markedly attenuated the neurotoxic actions of venom (3 µg/mL) and prevented the inhibition of contractile responses to ACh, CCh, and KCl. The addition of Chinese antivenom (0.75 U/mL) or Australian polyvalent antivenom (3 U/mL) at the t₉₀ time point after the addition of venom (3 µg/mL), partially reversed the inhibition of twitches and significantly reversed the venom-induced inhibition of responses to ACh and CCh, but had no significant effect on the response to KCl. Venom (30 µg/mL) also abolished direct twitches in the chick biventer cervicis nerve-muscle preparation and caused a significant increase in baseline tension, further indicative of myotoxicity. N. atra antivenom (4 U/mL) prevented the myotoxic effects of venom (30 µg/mL). However, G. brevicaudus antivenom (24 U/mL), D. acutus antivenom (8 U/mL) and Australian polyvalent snake antivenom (33 U/mL) were unable to prevent venom (30 µg/mL) induced myotoxicity. In the L6 rat skeletal muscle myoblast cell line, N. atra venom caused concentration-dependent inhibition of cell viability, with a half maximal inhibitory concentration (IC₅₀) of 2.8 ± 0.48 μg/mL. N. atra antivenom significantly attenuated the cytotoxic effect of the venom, whereas Australian polyvalent snake antivenom was less effective but still attenuated the cytotoxic effects at lower venom concentrations. Neither G. brevicaudus antivenom or D. acutus antivenom were able to prevent the cytotoxicity. This study indicates that Chinese N. atra monovalent antivenom is efficacious against the neurotoxic, myotoxic and cytotoxic effects of N. atra venom but the clinical effectiveness of the antivenom is likely to be diminished, even if given early after envenoming. The use of Chinese viper antivenoms (i.e., G. brevicaudus and D. acutus antivenoms) in cases of envenoming by the Chinese cobra is not supported by the results of the current study.

Kinetic analysis of effects of temperature and time on the regulation of venom expression in Bungarus multicinctus

Venom gland is a highly efficient venom production system to maintain their predatory arsenal. Venom toxins mRNA has been shown to increase abruptly in snake after venom expenditure, while the dynamics of venom accumulation during synthesis are poorly understood. Here, PacBio long-read sequencing, Illumina RNA sequencing (RNA-seq), and label-free proteome quantification were used to investigate the composition landscape and time- and temperature-dependent dynamics changes of the Bungarus multicinctus venom gland system. Transcriptome data (19.5223 Gb) from six adult B. multicinctus tissues were sequenced using PacBio RS II to generate a reference assembly, and average 7.28 Gb of clean RNA-seq data was obtained from venom glands by Illumina sequencing. Differentially expressed genes (DEGs) mainly were protein processing rather than venom toxins. Post-translational modifications provided the evidence of the significantly different proportions of toxins in the venom proteome with the changing of replenishment time and temperature, but constant of venom toxins mRNA in the venom gland transcriptome. Dynamic of toxins and genes involved in venom gland contraction suggesting the formation of the mature venom gland system would take at least 9 days. In addition, 59 toxin processing genes were identified, peptidylprolyl isomerase B of which underwent positive selection in Toxicofera. These results provide a reference for determining the extraction time of venom, production of polyclonal and monoclonal antibody for precise treatment plans of venomous snakebites, and construction of an in vitro synthesis system for snake venom protein.

Preparation and detection of sea snake antisera raised in rabbits

Antivenoms are currently the most effective medication used in the treatment of snakebites. However, there were relatively few studies on preparation of antivenoms targeting sea snakes, especially common sea snakes in China. In this study, we sought to prepare and detect mono- and bispecific antisera raised in rabbits against venoms of two sea snakes, Hydrophis cyanocinctus and H. curtus. The results of enzyme-linked immunosorbent assay showed that the rabbit antisera generally showed clearly detectable immunological cross-reactions after the third immunization and indicated that the strength of cross-reactions increased with an increase in the immunizing dose. Proteins within the H. cyanocinctus and H. curtus venoms showed similar profiles and were mainly concentrated in the low-molecular-weight region (8-25 kDa). Western blotting results revealed that the bands of these low-molecular weight proteins were dense and showed strong immunogenicity. Although we detected comparatively few bands of the high-molecular-weight proteins, these also showed strong immunogenicity. Our results indicate that both mono- and bispecific antisera both can neutralize H. cyanocinctus and H. curtus venoms, and in this regard, the monospecific H. curtus and bispecific antiserum were found to be superior to the H. cyanocinctus antiserum. Given the increasing frequency of snakebites worldwide, we believe that the findings of this study will have high practical applicability.

An in vitro α-neurotoxin-nAChR binding assay correlates with lethality and in vivo neutralization of a large number of elapid neurotoxic snake venoms from four continents

The aim of this study was to develop an in vitro assay for use in place of in vivo assays of snake venom lethality and antivenom neutralizing potency. A novel in vitro assay has been developed based on the binding of post-synaptically acting α-neurotoxins to nicotinic acetylcholine receptor (nAChR), and the ability of antivenoms to prevent this binding. The assay gave high correlation in previous studies with the in vivo murine lethality tests (Median Lethal Dose, LD50), and the neutralization of lethality assays (Median Effective Dose, ED50) by antisera against Naja kaouthia, Naja naja and Bungarus candidus venoms. Here we show that, for the neurotoxic venoms of 20 elapid snake species from eight genera and four continents, the in vitro median inhibitory concentrations (IC50s) for α-neurotoxin binding to purified nAChR correlated well with the in vivo LD50s of the venoms (R2 = 0.8526, p < 0.001). Furthermore, using this assay, the in vitro ED50s of a horse pan-specific antiserum against these venoms correlated significantly with the corresponding in vivo murine ED50s, with R2 = 0.6896 (p < 0.01). In the case of four elapid venoms devoid or having a very low concentration of α-neurotoxins, no inhibition of nAChR binding was observed. Within the philosophy of 3Rs (Replacement, Reduction and Refinement) in animal testing, the in vitro α-neurotoxin-nAChR binding assay can effectively substitute the mouse lethality test for toxicity and antivenom potency evaluation for neurotoxic venoms in which α-neurotoxins predominate. This will greatly reduce the number of mice used in toxicological research and antivenom production laboratories. The simpler, faster, cheaper and less variable in vitro assay should also expedite the development of pan-specific antivenoms against various medically important snakes in many parts of the world.

A pan-specific antiserum produced by a novel immunization strategy shows a high spectrum of neutralization against neurotoxic snake venoms

Snakebite envenomation is a neglected tropical disease of high mortality and morbidity largely due to insufficient supply of effective and affordable antivenoms. Snake antivenoms are mostly effective against the venoms used in their production. It is thus crucial that effective and affordable antivenom(s) with wide para-specificity, capable of neutralizing the venoms of a large number of snakes, be produced. Here we studied the pan-specific antiserum prepared previously by a novel immunization strategy involving the exposure of horses to a ‘diverse toxin repertoire’ consisting of 12 neurotoxic Asian snake toxin fractions/ venoms from six species. This antiserum was previously shown to exhibit wide para-specificity by neutralizing 11 homologous and 16 heterologous venoms from Asia and Africa. We now show that the antiserum can neutralize 9 out of 10 additional neurotoxic venoms. Altogether, 36 snake venoms belonging to 10 genera from 4 continents were neutralized by the antiserum. Toxin profiles previously generated using proteomic techniques of these 36 venoms identified α-neurotoxins, β-neurotoxins, and cytotoxins as predominant toxins presumably neutralized by the antiserum. The bases for the wide para-specificity of the antiserum are discussed. These findings indicate that it is feasible to generate antivenoms of wide para-specificity against elapid neurotoxic venoms from different regions in the world and raises the possibility of a universal neurotoxic antivenom. This should reduce the mortality resulting from neurotoxic snakebite envenomation.

Comparative venom proteomics of banded krait (Bungarus fasciatus) from five geographical locales: Correlation of venom lethality, immunoreactivity and antivenom neutralization

The banded krait, Bungarus fasciatus is a medically important venomous snake in Asia. The wide distribution of this species in Southeast Asia and southern China indicates potential geographical variation of the venom which may impact the clinical management of snakebite envenomation. This study investigated the intraspecific venom variation of B. fasciatus from five geographical locales through a venom decomplexing proteomic approach, followed by toxinological and immunological studies. The venom proteomes composed of a total of 9 toxin families, comprising 22 to 31 proteoforms at varying abundances. The predominant proteins were phospholipase A₂ (including beta-bungarotoxin), Kunitz-type serine protease inhibitor (KSPI) and three-finger toxins (3FTx), which are toxins that cause neurotoxicity and lethality. The venom lethality varied with geographical origins of the snake, with intravenous median lethal doses (LD₅₀) ranging from 0.45-2.55 µg/g in mice. The Thai Bungarus fasciatus monovalent antivenom (BFMAV) demonstrated a dose-dependent increasing immunological binding activity toward all venoms; however, its in vivo neutralization efficacy varied vastly with normalized potency values ranging from 3 to 28 mg/g, presumably due to the compositional differences of dominant proteins in the different venoms. The findings support that antivenom use should be optimized in different geographical areas. The development of a pan-regional antivenom may be a more sustainable solution for the treatment of snakebite envenomation.

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.