Proteomic analysis and antivenomics study of Western India Naja naja venom: correlation between venom composition and clinical manifestations of cobra bite in this region

Assessing Target Specificity of the Small Molecule Inhibitor MARIMASTAT to Snake Venom Toxins: A Novel Application of Thermal Proteome Profiling

New treatments that circumvent the pitfalls of traditional antivenom therapies are critical to address the problem of snakebite globally. Numerous snake venom toxin inhibitors have shown promising cross-species neutralization of medically significant venom toxins in vivo and in vitro. The development of high-throughput approaches for the screening of such inhibitors could accelerate their identification, testing, and implementation and thus holds exciting potential for improving the treatments and outcomes of snakebite envenomation worldwide. Energetics-based proteomic approaches, including thermal proteome profiling and proteome integral solubility alteration (PISA) assays, represent "deep proteomics" methods for high throughput, proteome-wide identification of drug targets and ligands. In the following study, we apply thermal proteome profiling and PISA methods to characterize the interactions between venom toxin proteoforms in Crotalus atrox (Western Diamondback Rattlesnake) and the snake venom metalloprotease (SVMP) inhibitor marimastat. We investigate its venom proteome-wide effects and characterize its interactions with specific SVMP proteoforms, as well as its potential targeting of non-SVMP venom toxin families. We also compare the performance of PISA thermal window and soluble supernatant with insoluble precipitate using two inhibitor concentrations, providing the first demonstration of the utility of a sensitive high-throughput PISA-based approach to assess the direct targets of small molecule inhibitors for snake venom.

The Need for Next-Generation Antivenom for Snakebite Envenomation in India

The limitations posed by currently available antivenoms have emphasized the need for alternative treatments to counteract snakebite envenomation. Even though exact epidemiological data are lacking, reports have indicated that most global snakebite deaths are reported in India. Among the many problems associated with snakebite envenomation, issues related to the availability of safer and more efficient antivenoms are of primary concern. Since India has the highest number of global snakebite deaths, efforts should be made to reduce the burden associated with snakebite envenoming. Alternative methods, including aptamers, camel antivenoms, phage display techniques for generating high-affinity antibodies and antibody fragments, small-molecule inhibitors, and natural products, are currently being investigated for their effectiveness. These alternative methods have shown promise in vitro, but their in vivo effectiveness should also be evaluated. In this review, the issues associated with Indian polyvalent antivenoms in neutralizing venom components from geographically distant species are discussed in detail. In a nutshell, this review gives an overview of the current drawbacks of using animal-derived antivenoms and several alternative strategies that are currently being widely explored.

The concept of Big Four: Road map from snakebite epidemiology to antivenom efficacy

Snake envenomation is a life-threatening disease caused by the injection of venom toxins from the venomous snake bite. Snakebite is often defined as the occupational or domestic hazard mostly affecting the rural population. India experiences a high number of envenoming cases and fatality due to the nation's diversity in inhabiting venomous snakes. The Indian Big Four snakes namely Russell's viper (Daboia russelii), spectacled cobra (Naja naja), common krait (Bungarus caeruleus), and saw-scaled viper (Echis carinatus) are responsible for majority of the snake envenoming cases and death. The demographic characteristics including occupation, stringent snake habitat management, poor healthcare facilities and ignorance of the rural victims are the primary influencers of high mortality. Biogeographic venom variation greatly influences the clinical pathologies of snake envenomation. The current antivenoms against the Big Four snakes are found to be less immunogenic against the venom toxins emphasizing the necessity of alternative approaches for antivenom generation. This review summarizes the burden of snake envenomation in India by the Big Four snakes including the geographic distribution of snake species and biogeographic venom variation. We have provided comprehensive information on snake venom proteomics that has aided the better understanding of venom induced pathological features, summarized the impact of current polyvalent antivenom therapy highlighting the need for potential antivenom treatment for the effective management of snakebites.

Application of an Extracellular Matrix-Mimicking Fluorescent Polymer for the Detection of Proteolytic Venom Toxins

The cytotoxicity caused by snake venoms is a serious medical problem that greatly contributes to the morbidity observed in snakebite patients. The cytotoxic components found in snake venoms belong to a variety of toxin classes and may cause cytotoxic effects by targeting a range of molecular structures, including cellular membranes, the extracellular matrix (ECM) and the cytoskeleton. Here, we present a high-throughput assay (384-well plate) that monitors ECM degradation by snake venom toxins via the application of fluorescent versions of model ECM substrates, specifically gelatin and collagen type I. Both crude venoms and fractionated toxins of a selection of medically relevant viperid and elapid species, separated via size-exclusion chromatography, were studied using the self-quenching, fluorescently labelled ECM-polymer substrates. The viperid venoms showed significantly higher proteolytic degradation when compared to elapid venoms, although the venoms with higher snake venom metalloproteinase content did not necessarily exhibit stronger substrate degradation than those with a lower one. Gelatin was generally more readily cleaved than collagen type I. In the viperid venoms, which were subjected to fractionation by SEC, two (B. jararaca and C. rhodostoma, respectively) or three (E. ocellatus) active proteases were identified. Therefore, the assay allows the study of proteolytic activity towards the ECM in vitro for crude and fractionated venoms.

Multilevel Comparison of Indian Naja Venoms and Their Cross-Reactivity with Indian Polyvalent Antivenoms

Snake envenoming is caused by many biological species, rather than a single infectious agent, each with a multiplicity of toxins in their venom. Hence, developing effective treatments is challenging, especially in biodiverse and biogeographically complex countries such as India. The present study represents the first genus-wide proteomics analysis of venom composition across Naja species (N. naja, N. oxiana, and N. kaouthia) found in mainland India. Venom proteomes were consistent between individuals from the same localities in terms of the toxin families present, but not in the relative abundance of those in the venom. There appears to be more compositional variation among N. naja from different locations than among N. kaouthia. Immunoblotting and in vitro neutralization assays indicated cross-reactivity with Indian polyvalent antivenom, in which antibodies raised against N. naja are present. However, we observed ineffective neutralization of PLA2 activities of N. naja venoms from locations distant from the source of immunizing venoms. Antivenom immunoprofiling by antivenomics revealed differential antigenicity of venoms from N. kaouthia and N. oxiana, and poor reactivity towards 3FTxs and PLA2s. Moreover, there was considerable variation between antivenoms from different manufacturers. These data indicate that improvements to antivenom manufacturing in India are highly desirable.

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.

Pharmacological re-assessment of traditional medicinal plants-derived inhibitors as antidotes against snakebite envenoming: A critical review

ETHNOPHARMACOLOGICAL RELEVANCE

Traditional healers have used medicinal plants to treat snakebite envenomation worldwide; however, mostly without scientific validation. There have been many studies on the therapeutic potential of the natural products against snake envenomation.

AIM OF THE STUDY

This review has highlighted snake venom inhibitory activity of bioactive compounds and peptides from plants that have found a traditional use in treating snakebite envenomation. We have systematically reviewed the scenario of different phases of natural snake venom inhibitors characterization covering a period from 1994 until the present and critically analysed the lacuna of the studies if any, and further scope for their translation from bench to bedside.

MATERIALS AND METHODS

The medicinal plant-derived compounds used against snakebite therapy were reviewed from the available literature in public databases (Scopus, MEDLINE) from 1994 till 2020. The search words used were 'natural inhibitors against snakebite,' 'natural products as therapeutics against snakebite,' 'natural products as antidote against snake envenomation,' ' snake venom toxin natural inhibitors,' 'snake venom herbal inhibitors'. However, the scope of this review does not include computational (in silico) predictions without any wet laboratory validation and snake venom inhibitory activity of the crude plant extracts. In addition, we have also predicted the ADMET properties of the identified snake venom inhibitors to highlight their valuable pharmacokinetics for future clinical studies.

RESULTS

The therapeutic application of plant-derived natural inhibitors to treat snakebite envenomation as an auxiliary to antivenom therapy has been gaining significant momentum. Pharmacological reassessment of the natural compounds derived from traditional medicinal plants has demonstrated inhibition of the principal toxic enzymes of snake venoms at various extents to curb the lethal and/or deleterious effects of venomous snakebite. Nevertheless, such molecules are yet to be commercialized for clinical application in the treatment of snakebite. There are many obstacles in the marketability of the plant-derived natural products as snake envenomation antidote and strategies must be explored for the translation of these compounds from drug candidates to their clinical application.

CONCLUSION

In order to minimize the adverse implications of snake envenomation, strategies must be developed for the smooth transition of these plant-derived small molecule inhibitors from bench to bedside. In this article we have presented an inclusive review and have critically analysed natural products for their therapeutic potential against snake envenomation, and have proposed a road map for use of natural products as antidote against snakebite.

The chemistry of snake venom and its medicinal potential

The fascination and fear of snakes dates back to time immemorial, with the first scientific treatise on snakebite envenoming, the Brooklyn Medical Papyrus, dating from ancient Egypt. Owing to their lethality, snakes have often been associated with images of perfidy, treachery and death. However, snakes did not always have such negative connotations. The curative capacity of venom has been known since antiquity, also making the snake a symbol of pharmacy and medicine. Today, there is renewed interest in pursuing snake-venom-based therapies. This Review focuses on the chemistry of snake venom and the potential for venom to be exploited for medicinal purposes in the development of drugs. The mixture of toxins that constitute snake venom is examined, focusing on the molecular structure, chemical reactivity and target recognition of the most bioactive toxins, from which bioactive drugs might be developed. The design and working mechanisms of snake-venom-derived drugs are illustrated, and the strategies by which toxins are transformed into therapeutics are analysed. Finally, the challenges in realizing the immense curative potential of snake venom are discussed, and chemical strategies by which a plethora of new drugs could be derived from snake venom are proposed.

A comparison of two different analytical workflows to determine the venom proteome composition of Naja kaouthia from North-East India and immunological profiling of venom against commercial antivenoms

The Indian monocled cobra (Naja kaouthia) is one of the most prevalent venomous snakes in northeast India (NEI) and is the cause of many fatalities. The composition of NEI N. kaouthia venom (NkV) was deciphered using two different proteomic approaches: (i) 1D SDS-PAGE coupled to label-free quantification of protein bands using stringent identification criteria and (ii) reversed-phase high-performance liquid chromatography (RP-HPLC) followed by quantification based on area under the RP-HPLC peaks. The proteomic data from both strategies were compared. Proteomic analyses from both workflows identified 32 proteins (toxins) distributed over 10-14 snake venom protein families in NEI NkV. The relative abundances of the venom proteins determined from the analytical workflows coincided with the densitometry band intensities of the NEI NkV. Phospholipase A₂ (13.1-16.0%) and three-finger toxins (58.5-64.2%) represented the most abundant enzymatic and non-enzymatic proteins in NEI NkV, respectively. Immuno-cross-reactivity studies by enzyme-linked immunoassay and immunoblot analyses pointed to the poor efficacy of commercial PAVs in recognizing the low molecular mass (<15 kDa) toxins of NEI NkV. Spectrofluorometric titration determined the presence of NEI NkV-specific antibodies in commercial PAV, at a level that was higher than that previously reported for eastern India NkV-specific antibodies in commercial antivenom.

Interpopulational variation and ontogenetic shift in the venom composition of Lataste's viper (Vipera latastei, Boscá 1878) from northern Portugal

Lataste's viper (Vipera latastei) is a venomous European viper endemic to the Iberian Peninsula, recognised as medically important by the World Health Organization. To date, no comprehensive characterisation of this species' venom has been reported. Here, we analysed the venoms of juvenile and adult specimens of V. latastei from two environmentally different populations from northern Portugal. Using bottom-up venomics, we produced six venom proteomes (three per population) from vipers belonging to both age classes (i.e., two juveniles and four adults), and RP-HPLC profiles of 54 venoms collected from wild specimens. Venoms from juveniles and adults differed in their chromatographic profiles and relative abundances of their toxins, suggesting the occurrence of ontogenetic changes in venom composition. Specifically, snake venom metalloproteinase (SVMP) was the most abundant toxin family in juvenile venoms, while snake venom serine proteinases (SVSPs), phospholipases A₂ (PLA₂s), and C-type lectin-like (CTLs) proteins were the main toxins comprising adult venoms. The RP-HPLC venom profiles were found to vary significantly between the two sampled localities, indicating geographic variability. Furthermore, the presence/absence of certain peaks in the venom chromatographic profiles appeared to be significantly correlated also to factors like body size and sex of the vipers. Our findings show that V. latastei venom is a variable phenotype. The intraspecific differences we detected in its composition likely mirror changes in the feeding ecology of this species, taking place during different life stages and under different environmental pressures. SIGNIFICANCE: Lataste's viper (Vipera latastei) is a medically important viper endemic to the Iberian Peninsula, inhabiting different habitats and undergoing a marked ontogenetic dietary shift. In the current study, we report the first proteomic analysis of V. latastei venom from two environmentally different localities in northern Portugal. Our bottom-up venomic analyses show that snake venom serine proteinases (SVSPs), phospholipases A2 (PLA2s), and C-type lectin-like (CTLs) proteins are the major components of adult V. latastei venom. The comparative analysis of young and adult venoms suggests the occurrence of ontogenetic shift in toxin abundances, with snake venom metalloproteinases (SVMPs) being the predominant toxins in juvenile venoms. Moreover, geographic venom variation between the two studied populations is also detected, with our statistical analyses suggesting that factors like body size and sex of the vipers are possibly at play in its determination. Our work represents the first assessment of the composition of V. latastei venom, and the first step towards a better understanding of the drivers behind its variability.

Lipid-induced monokine cyclophilin-A promotes adipose tissue dysfunction implementing insulin resistance and type 2 diabetes in zebrafish and mice models of obesity

Several studies have implicated obesity-induced macrophage-adipocyte cross-talk in adipose tissue dysfunction and insulin resistance. However, the molecular cues involved in the cross-talk of macrophage and adipocyte causing insulin resistance are currently unknown. Here, we found that a lipid-induced monokine cyclophilin-A (CyPA) significantly attenuates adipocyte functions and insulin sensitivity. Targeted inhibition of CyPA in diet-induced obese zebrafish notably reduced adipose tissue inflammation and restored adipocyte function resulting in improvement of insulin sensitivity. Silencing of macrophage CyPA or pharmacological inhibition of CyPA by TMN355 effectively restored adipocytes' functions and insulin sensitivity. Interestingly, CyPA incubation markedly increased adipocyte inflammation along with an impairment of adipogenesis, however, mutation of its cognate receptor CD147 at P309A and G310A significantly waived CyPA's effect on adipocyte inflammation and its differentiation. Mechanistically, CyPA-CD147 interaction activates NF-κB signaling which promotes adipocyte inflammation by upregulating various pro-inflammatory cytokines gene expression and attenuates adipocyte differentiation by inhibiting PPARγ and C/EBPβ expression via LZTS2-mediated downregulation of β-catenin. Moreover, inhibition of CyPA or its receptor CD147 notably restored palmitate or CyPA-induced adipose tissue dysfunctions and insulin sensitivity. All these results indicate that obesity-induced macrophage-adipocyte cross-talk involving CyPA-CD147 could be a novel target for the management of insulin resistance and type 2 diabetes.

Venomics and antivenomics of Indian spectacled cobra (Naja naja) from the Western Ghats

Venom proteome profiling of Naja naja from the Western Ghats region in Kerala was achieved through SDS-PAGE and RP-HPLC followed by Q-TOF LC-MS/MS analysis, incorporating PEAKS and Novor assisted de novo sequencing methodologies. A total of 115 proteins distributed across 17 different enzymatic and non-enzymatic venom protein families were identified through conventional and 39 peptides through homology-driven proteomics approaches. Fourteen peptides derived through de novo complements the Mascot data indicating the importance of homology-driven approaches in improving protein sequence information. Among the protein families identified, glutathione peroxidase and endonuclease were reported for the first time in the Indian cobra venom. Immunological cross-reactivity assessed using Indian polyvalent antivenoms suggested that VINS showed better EC₅₀ (2.48 µg/mL) value than that of PSAV (6.04 µg/mL) and Virchow (6.03 µg/mL) antivenoms. Western blotting experiments indicated that all the antivenoms elicited poor binding specificities, especially towards low molecular mass proteins. Second-generation antivenomics studies revealed that VINS antivenom was less efficient to detect many low molecular mass proteins such as three-finger toxins and Kunitz-type serine protease Inhibitors. Taken together, the present study enabled a large-scale characterization of the venom proteome of Naja naja from the Western Ghats and emphasized the need for developing more efficient antivenoms.

Immunorecognition capacity of Indian polyvalent antivenom against venom toxins from two populations of Echis carinatus

Clinicians report low efficacy of Indian polyvalent antivenom (PAV), with >20 vials required for treatment of a snakebite envenoming. We hypothesize that the antivenom efficacy could be reduced due to insufficient antibodies against some venom toxins. To test this, we used third-generation antivenomics to reveal bound and unbound venom toxins of Echis carinatus venom from Goa (ECVGO) and Tamil Nadu (ECVTN). We used 60, 120, 180, 240, 300, and 360 μg of venom and passed through mini-columns containing ~5 mg Antivenom bound to CNBr beads. The non-retained (unbound) and retained (bound) toxins were identified using reverse-phase HPLC and tandem mass spectrometry. Low molecular weight toxins - Short disintegrins (5.3 kDa) and DIS domain of P-II SVMP from ECVGO and ECVTN showed poor binding with antivenom. The immunorecognition sites of antivenom saturated at the lower antivenom-venom ratio for ECVGO than for ECVTN. The immunoretained capacity of antivenom against ECVTN was 140.6 μg and ECVGO was 125.1 μg. The amount of immunoretained toxins quantified can further be used to estimate the efficacy of antivenom by correlating it with in-vivo studies. The unbound toxins identified from this study could be targeted to improve the effectiveness of antivenom.

Advances in the Therapeutic Application of Small-Molecule Inhibitors and Repurposed Drugs against Snakebite

The World Health Organization has declared snakebite as a neglected tropical disease. Antivenom administration is the sole therapy against venomous snakebite; however, several limitations of this therapy reinforce the dire need for an alternative and/or additional treatment against envenomation. Inhibitors against snake venoms have been explored from natural resources and are synthesized in the laboratory; however, repurposing of small-molecule therapeutics (SMTs) against the principal toxins of snake venoms to inhibit their lethality and/or obnoxious effect of envenomation has been garnering greater attention owing to their established pharmacokinetic properties, low-risk attributes, cost-effectiveness, ease of administration, and storage stability. Nevertheless, SMTs are yet to be approved and commercialized for snakebite treatment. Therefore, we have systematically reviewed and critically analyzed the scenario of small synthetic inhibitors and repurposed drugs against snake envenomation from 2005 to date and proposed novel approaches and commercialization strategies for the development of efficacious therapies against snake envenomation.

Assessment of quality and pre-clinical efficacy of a newly developed polyvalent antivenom against the medically important snakes of Sri Lanka

Snake envenomation is a severe problem in Sri Lanka (SL) and Indian polyvalent antivenom (PAV) is mostly used for treating snakebite albeit due to geographical variation in venom composition, Indian PAV shows poor efficacy in neutralizing the lethality and toxicity of venom from the same species of snakes in SL. Therefore, the quality and in vivo venom neutralization potency of a country-specific PAV produced against the venom of the five most medically important snakes of SL (Daboia russelii, Echis carinatus, Hypnale hypnale, Naja naja, Bungarus caeruleus) was assessed. LC-MS/MS analysis of two batches of PAV showed the presence of 88.7-97.2% IgG and traces of other plasma proteins. The tested PAVs contained minor amounts of undigested IgG and F(ab')2 aggregates, showed complement activation, were devoid of IgE, endotoxin, and content of preservative was below the threshold level. Immunological cross-reactivity and in vitro neutralization of enzymatic activities, pharmacological properties demonstrated superior efficacy of SL PAV compared to Indian PAV against SL snake venoms. The in vivo neutralization study showed that the tested PAVs are potent to neutralize the lethality and venom-induced toxicity of SL snake venoms. Therefore, our study suggests that introduction of SL-specific PAV will improve snakebite management in SL.

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.

The application of laboratory-based analytical tools and techniques for the quality assessment and improvement of commercial antivenoms used in the treatment of snakebite envenomation

Snakebite envenomation is a public health problem of high impact, particularly for the developing world. Antivenom, which contains whole or protease-digested immunoglobulin G, purified from the plasma of hyper-immunized animals (mainly horses), is the mainstay for the treatment of snakebite envenomation. The success of antivenom therapy depends upon its ability to abrogate or reduce the local and systemic toxicity of envenomation. In addition, antivenom administration must be safe for the patients. Therefore, antivenom manufacturers must ensure that these products are effective and safe in the treatment of envenomations. Antivenom efficacy and safety are determined by the physicochemical characteristics of formulations, purity of the immunoglobulin fragments and antibodies, presence of protein aggregates, endotoxin burden, preservative load, and batch to batch variation, as well as on the ability to neutralize the most important toxins of the venoms against which the antivenom is designed. In this context, recent studies have shown that laboratory-based simple analytical techniques, for example, size exclusion chromatography, sodium dodecyl sulphate polyacrylamide gel electrophoresis, mass spectrometry, immunological profiling including immuno-turbidimetry and enzyme-linked immunosorbent assays, Western blotting, immune-chromatographic technique coupled to mass spectrometry analysis, reverse-phase high performance liquid chromatography, spectrofluorometric analysis, in vitro neutralization of venom enzymatic activities, and other methodologies, can be applied for the assessment of antivenom quality, safety, stability, and efficacy. This article reviews the usefulness of different analytical techniques for the quality assessment of commercial antivenoms. It is suggested that these tests should be applied for screening the quality of commercial antivenoms before their preclinical and clinical assessment.

Biogeographic venom variation in Russell's viper (Daboia russelii) and the preclinical inefficacy of antivenom therapy in snakebite hotspots

BACKGROUND

Snakebite in India results in over 58,000 fatalities and a vast number of morbidities annually. The majority of these clinically severe envenomings are attributed to Russell's viper (Daboia russelii), which has a near pan-India distribution. Unfortunately, despite its medical significance, the influence of biogeography on the composition and potency of venom from disparate D. russelii populations, and the repercussions of venom variation on the neutralisation efficacy of marketed Indian antivenoms, remain elusive.

METHODS

Here, we employ an integrative approach comprising proteomic characterisation, biochemical analyses, pharmacological assessment, and venom toxicity profiling to elucidate the influence of varying ecology and environment on the pan-Indian populations of D. russelii. We then conducted in vitro venom recognition experiments and in vivo neutralisation assays to evaluate the efficacy of the commercial Indian antivenoms against the geographically disparate D. russelii populations.

FINDINGS

We reveal significant intraspecific variation in the composition, biochemical and pharmacological activities and potencies of D. russelii venoms sourced from five distinct biogeographic zones across India. Contrary to our understanding of the consequences of venom variation on the effectiveness of snakebite therapy, commercial antivenom exhibited surprisingly similar neutralisation potencies against the majority of the investigated populations, with the exception of low preclinical efficacy against the semi-arid population from northern India. However, the ability of Indian antivenoms to counter the severe morbid effects of Daboia envenoming remains to be evaluated.

CONCLUSION

The concerning lack of antivenom efficacy against the north Indian population of D. russelii, as well as against two other 'big four' snake species in nearby locations, underscores the pressing need to develop pan-India effective antivenoms with improved efficacy in high snakebite burden locales.

Biogeographical venom variation in the Indian spectacled cobra (Naja naja) underscores the pressing need for pan-India efficacious snakebite therapy

Background

Snake venom composition is dictated by various ecological and environmental factors, and can exhibit dramatic variation across geographically disparate populations of the same species. This molecular diversity can undermine the efficacy of snakebite treatments, as antivenoms produced against venom from one population may fail to neutralise others. India is the world’s snakebite hotspot, with 58,000 fatalities and 140,000 morbidities occurring annually. Spectacled cobra (Naja naja) and Russell’s viper (Daboia russelii) are known to cause the majority of these envenomations, in part due to their near country-wide distributions. However, the impact of differing ecologies and environment on their venom compositions has not been comprehensively studied.

Methods

Here, we used a multi-disciplinary approach consisting of venom proteomics, biochemical and pharmacological analyses, and in vivo research to comparatively analyse N. naja venoms across a broad region (>6000 km; seven populations) covering India’s six distinct biogeographical zones.

Findings

By generating the most comprehensive pan-Indian proteomic and toxicity profiles to date, we unveil considerable differences in the composition, pharmacological effects and potencies of geographically-distinct venoms from this species and, through the use of immunological assays and preclinical experiments, demonstrate alarming repercussions on antivenom therapy. We find that commercially-available antivenom fails to effectively neutralise envenomations by the pan-Indian populations of N. naja, including a complete lack of neutralisation against the desert Naja population.

Conclusion

Our findings highlight the significant influence of ecology and environment on snake venom composition and potency, and stress the pressing need to innovate pan-India effective antivenoms to safeguard the lives, limbs and livelihoods of the country’s 200,000 annual snakebite victims.

Annually, India is burdened by the highest number of snake envenomations across the globe, with over 58,000 fatalities and three times the number of morbidities, predominantly affecting the rural agrarian communities. The spectacled cobra (Naja naja) and Russell’s viper (Daboia russelii) are responsible for the vast majority of envenomations in the country, in part, due to their near country-wide distributions. In this study, we unveil the astounding differences in venom composition of N. naja from six different biogeographical zones across the country (>6000 km). We provide a comprehensive account of their disparate venom proteomic profiles, biochemical and pharmacological effects, and the associated potencies. Our study uncovers alarming differences in the efficacy of the marketed polyvalent antivenoms in neutralising these venoms, thereby, emphasising the pressing need to develop dose-efficacious and pan-India effective antivenoms for the treatment of snakebites in the country. This study also highlights the significant influence of ecology and diverse environments on the venom variability, insinuating the necessity for innovating cost-effective and pan-India efficacious solutions to safeguard the lives, limbs and livelihoods of India’s two hundred thousand annual snakebite victims.

A Wolf in Another Wolf's Clothing: Post-Genomic Regulation Dictates Venom Profiles of Medically-Important Cryptic Kraits in India

The Common Krait (Bungarus caeruleus) shares a distribution range with many other ‘phenotypically-similar’ kraits across the Indian subcontinent. Despite several reports of fatal envenomings by other Bungarus species, commercial Indian antivenoms are only manufactured against B. caeruleus. It is, therefore, imperative to understand the distribution of genetically distinct lineages of kraits, the compositional differences in their venoms, and the consequent impact of venom variation on the (pre)clinical effectiveness of antivenom therapy. To address this knowledge gap, we conducted phylogenetic and comparative venomics investigations of kraits in Southern and Western India. Phylogenetic reconstructions using mitochondrial markers revealed a new species of krait, Romulus’ krait (Bungarus romulusi sp. nov.), in Southern India. Additionally, we found that kraits with 17 mid-body dorsal scale rows in Western India do not represent a subspecies of the Sind Krait (B. sindanus walli) as previously believed, but are genetically very similar to B. sindanus in Pakistan. Furthermore, venom proteomics and comparative transcriptomics revealed completely contrasting venom profiles. While the venom gland transcriptomes of all three species were highly similar, venom proteomes and toxicity profiles differed significantly, suggesting the prominent role of post-genomic regulatory mechanisms in shaping the venoms of these cryptic kraits. In vitro venom recognition and in vivo neutralisation experiments revealed a strong negative impact of venom variability on the preclinical performance of commercial antivenoms. While the venom of B. caeruleus was neutralised as per the manufacturer’s claim, performance against the venoms of B. sindanus and B. romulusi was poor, highlighting the need for regionally-effective antivenoms in India.

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.

Functional venomics of the Big-4 snakes of Pakistan

In South Asia, the "Big-4" venomous snakes Naja naja, Bungarus caeruleus, Daboia russelii, and Echis carinatus are so-called because they are the most medically important snakes in the region. Antivenom is the only effective treatment option for snakebite envenoming but antivenom is not produced domestically in Pakistan making the country reliant on polyvalent products imported from India and Saudi Arabia. The present study investigated the toxin composition and activity of the venoms of Pakistani specimens by means of proteomic and physio/pharmacological experiments. To evaluate the composition of venoms, 1D/2D-PAGE of crude venoms and RP-HPLC followed by SDS-PAGE were performed. Enzymatic, hemolytic, coagulant and platelet aggregating activities of crude venoms were assayed and were concordant with expectations based on the abundance of protein species in each. Neutralization assays were performed using Bharat polyvalent antivenom (BPAV), a product raised against venoms from Big-4 specimens from southern India. BPAV exhibited cross-reactivity against the Pakistani venoms, however, neutralization of clinically relevant activities was variable and rarely complete. Cumulatively, the presented data not only highlight geographical variations present in the venoms of the Big-4 snakes of South Asia, but also demonstrate the neutralization potential of Indian polyvalent against the venom of Pakistani specimens. Given the partial neutralization observed, it is clear that whilst BPAV is a life-saving product in Pakistan, in future it is hoped that a region-specific product might be manufactured domestically, using venoms of local snakes in the immunising mixture.

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.

The action of Echis carinatus and Naja naja venoms on human neutrophils; an emphasis on NETosis

BACKGROUND

Neutrophils are the first line defense cells of the innate immunity. As a final defense, they discharge their de-condensed chromatin/DNA fibers, the NETs (Neutrophil Extracellular Traps), by a process called NETosis. Two types of NETosis have been currently described: the suicidal/delayed/classical-type, which is ROS dependent that results in the ejection of nuclear DNA, and the vital/rapid/early-type, which may or may not require ROS but, eject nuclear/mitochondrial DNA or both. Thus, Echis carinatus and Naja naja venoms are comparatively studied for their NET inducing property.

METHODS

Formation of NETs, cell viability, ROS, and Ca²⁺ levels are estimated. An in vivo toxicity study and possible cellular signaling have been addressed using immunoblots and pharmacological inhibitors.

RESULTS

E. carinatus and N. naja venoms respectively induce suicidal and vital NETosis. E. carinatus venom induces NETosis by activating NOX and PAD-4 enzymes in a ROS dependent manner via PKC/ERK/JNK signaling axis, while N. naja venom does it by activating PAD-4 enzyme, but independent of ROS requirement and as well as PKC/ERK/JNK activation.

CONCLUSION

For the first time our study demonstrates the distinct action of E. carinatus and N. naja venoms on the process of NETosis. NETosis being a newly explored area in snake venom pharmacodynamics, it is important to study its impact on the various pathophysiological properties induced by snake venoms.

SIGNIFICANCE

Understanding the varied actions of snake venoms on neutrophils/blood cells and the role of DNase are likely to provide insights for better management of snakebite pathophysiology.