Next Article in Journal
Degradation of Patulin in Pear Juice and Apple Juice by Ascorbic Acid and the Combination of Ascorbic Acid and Ferrous Iron
Previous Article in Journal
Acute and Subacute Toxicity Studies of Erodium guttatum Extracts by Oral Administration in Rodents
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Toxins
Volume 14
Issue 11
10.3390/toxins14110736
toxins-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

by
Arif Arrahman
1,2,3,
Taline D. Kazandjian
4,
Kristina B. M. Still
1,2,
Julien Slagboom
1,2,
Govert W. Somsen
1,2,
Freek J. Vonk
1,5,
Nicholas R. Casewell
4 and
Jeroen Kool
1,2,*
1
Division of Bioanalytical Chemistry, Department of Chemistry and Pharmaceutical Sciences, Faculty of Sciences, Amsterdam Institute of Molecular and Life Sciences (AIMMS), Vrije Universiteit Amsterdam, De Boelelaan 1083, 1081 HV Amsterdam, The Netherlands
2
Centre for Analytical Sciences Amsterdam (CASA), 1012 WX Amsterdam, The Netherlands
3
Faculty of Pharmacy, Universitas Indonesia, Kampus Baru UI, Depok 16424, Indonesia
4
Centre for Snakebite Research and Interventions. Liverpool School of Tropical Medicine, Pembroke Place, Liverpool L3 5QA, UK
5
Naturalis Biodiversity Centre, Darwinweg 2, 2333 CR Leiden, The Netherlands
*
Author to whom correspondence should be addressed.
Toxins 2022, 14(11), 736; https://doi.org/10.3390/toxins14110736
Submission received: 16 September 2022 / Revised: 17 October 2022 / Accepted: 21 October 2022 / Published: 27 October 2022
(This article belongs to the Special Issue New Trends in Identification and Characterization of Venom Components)
Browse Figures
Versions Notes

Abstract

:
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.
Keywords:
snakebite; coagulopathy; varespladib; marimastat; Dendroaspis; Naja; mass spectrometry
Key Contribution: PLA2s from the cobra venoms inflict anticoagulant effects on bovine plasma, the effects of which are significantly reduced with the addition of varespladib. The key toxin responsible for anticoagulant activity in mamba venoms is still unknown.

1. Introduction

Snakebite is a neglected tropical disease, affecting between 421,000 [1] and 2.7 million individuals [2] each year and resulting in 81,000–138,000 deaths. Estimates of global envenoming rates fluctuate dramatically due to the challenges of collecting consistent and reliable data, which relies on the reporting of snakebite incidences and mortality from both individuals in the community and hospitals [3]. The effects of snakebite are most pronounced in low and middle-income regions [4]. This is due to a combination of factors such as lack of antivenom availability, limitations relating to healthcare infrastructure, accessibly and appropriately trained medical staff, and the ability to afford effective treatment and rehabilitation. Such factors have led to an estimated 93 million people living in locations classed as being vulnerable to snakebite [5], culminating in snakebite being listed by the World Health Organization as a priority neglected tropical disease.
One of the most medically important groups of snakes are the elapids (Serpentes: Elapidae), a snake family that consists of over 300 species, of which the World Health Organization considers around 80 to have medical importance in the context of human snakebite [2,6]. Elapid venoms are complex mixtures of pharmacologically active proteins and polypeptides that play an essential role in the incapacitation, immobilisation, and digestion of prey [7,8,9]. These venoms tend to be dominated by the presence of three-finger toxins (3FTXs) and/or phospholipase A2 (PLA2) toxins [6,10], which often act upon the neuromuscular junction to cause paralysis and respiratory failure [11]. Arguably the most medically important elapids in sub-Saharan Africa are the cobras and near relatives (genus Naja and Hemachatus haemachatus) and the mambas (genus Dendroaspis) [12]. While both groups of these elapid snakes employ a subfamily of venom 3FTXs known as α-neurotoxins to incapacitate prey, cobra venoms also contain a large abundance of cytotoxic 3FTX which, in the context of spitting cobras at least, likely contribute to causing pain and cell damage [13,14], while mamba venoms contain large amounts of dendrotoxins [15], which contribute to neurotoxic symptoms via the blockage of potassium channels [16,17].
Though the predominant symptom of elapid envenoming is neurotoxicity, there is evidence that some of these venoms can additionally display hemotoxic effects. For example, some Australian elapid snakes have venoms that cause potent procoagulation by activating prothrombin [18,19,20,21]. However, certain elapid venoms have anticoagulant effects, which are achieved through several mechanisms: for example, the inhibition of factor X or thrombin, the lysis of fibrin, the activation of plasminogen and the inhibition of platelet aggregation [22,23,24,25,26,27]. Rather than acting directly to incapacitate prey, as observed for the procoagulant toxins of some Australian species, toxins that prevent coagulation could potentially aid in the spread of neurotoxins throughout the body leading to quicker systemic envenomation, though the pathological relevance of such effects in human envenomings remains unclear. Anticoagulant activity has been reported in the venoms of several cobra species [25,28,29]. Cobras contain several venom toxins with the potential for exerting anticoagulant activity, including PLA2s, which can act to inhibit platelet aggregation [30], thrombin [31,32], the tenase and prothrombinase complexes [33,34], and clotting factor Xa (FXa) [25] and 3FTXs, of which cytotoxic 3FTXs (CTXs) have been shown to lyse erythrocytes [35,36,37] and, in the venom of Hemachatus haemachatus, form a complex called hemextin A & B that inhibits the activity of the coagulation factor VIIa. However, recent studies have identified PLA2s as the likely causative agents, both via direct identification following toxin fractionation for Naja nigricollis [38], and via functional inhibition of PLA2 toxins for a much broader range of spitting cobra species [39]. However, there is a lack of such studies focusing on mamba venoms, though prior research suggests they might inhibit platelet aggregation [40].
Current snakebite treatments, known as antivenom and which consist of animal-derived polyclonal antibodies stimulated by snake venoms, are delivered intravenously to snakebite victims. There are several commercially available antivenoms produced for the treatment of cobra and mamba bites in India and Sub-Saharan Africa [41,42,43,44,45], produced using either a single species (monovalent) or covering a range of species (polyvalent). While these are life-saving therapies, they have several limitations associated with their use, including extremely high costs to consumers, limited cross-species efficacy [43,45], cold storage requirements, high incidences of adverse reactions, and a requirement to be delivered intravenously in a healthcare setting [46,47,48,49,50,51]. Additionally, there is a lack of preclinical testing for current antivenoms, whose production methods differ very little today than when the first antivenoms were produced, and the existing preclinical data are not standardised and often not publicly available [45]. Thus, research into refining antivenoms, or finding treatments that can be used synergistically with antivenoms to improve patient outcomes further, are a priority in the field.
One avenue in the pursuit of this goal is research into small molecule inhibitor s [52,53]. Varespladib [54] and marimastat [55] are small molecule drugs that inhibit the activity of PLA2s and matrix metalloproteases, respectively. Originally intended as treatments for inflammatory diseases and cancer, they have since been explored for repurposing as snakebite treatments. Varespladib in particular is a broad-spectrum inhibitor of venom PLA2s that has been shown to delay the lethal effects of venoms from the elapids Bungarus multicinctus, Oxuranus scutellatus, and Micrurus fulvius in mice through the inhibition of presynaptically acting PLA2s [56,57,58], with the injectable form temporarily reversing paralytic effects observed from the prior two species [57]. Furthermore, it confers a significant survival benefit when administered to mice envenomed by lethal doses of venom from certain vipers [59,60,61] and elapids [61], and provides protection against the pathologies of Micrurus dumerilii carinicauda venom in rats when used in combination with antivenom [62]. In terms of its ability to treat haemotoxic pathologies, in vitro studies demonstrate the ability of varespladib to neutralise the anticoagulant activity of PLA2s isolated from the venoms of the vipers Echis carinatus and E. ocellatus, and the elapid Oxuranus scutellatus, and it additionally exhibits some neutralisation of the procoagulant activity of PLA2s from several viper species [63]. Additionally, it has been shown to neutralise PLA2-mediated prothrombinase inhibition in whole venoms from the viper genus Bitis [64] and elapid genus Pseudechis [39], as well as negate [63] the anticoagulant activity of several spitting cobra species [25,38,39], which could not be neutralised by antivenom alone [25,57]. Similarly, marimastat has been found to inhibit the activity of viper SVMP toxins and is capable of conferring preclinical protection against hemotoxic pathology caused by certain envenomings in small animal models, especially when combined with antivenom [53,65,66,67,68,69]. It also prevents haemorrhage in mice from venoms of the vipers Deinagkistrodon acutus and Agkistrodon halys [70], inhibits the fibrinogenolytic activity of CAMP-2 from Crotalus atrox venom [67], and neutralises the SVMP-mediated procoagulant effects of venom from the colubrids Rhabdophis tigrinus [67] and Dispholidus typus [68]. These studies provide enticing evidence to pursue small molecule drugs as (at least adjuvant) treatments for snakebite envenomation.
This study investigates the anticoagulant activity of venoms sourced from medically important cobra and mamba species on bovine plasma and seeks to identify the toxins responsible for such activity. In addition, we evaluated the role of the small molecule inhibitors varespladib and marimastat as potential inhibitors of the anticoagulant toxins found in the elapid venoms under study. The analytical approach involved a workflow consisting of three main steps. The first step was an initial screening of coagulopathic activity of the various crude venoms and then investigating the effect of the small molecules varespladib and marimastat on inhibiting the resulting coagulopathic activity. The second step was the characterisation of nanofractionated toxins in the presence and absence of the same small molecule toxins inhibitors. The third step involved proteomics-based identification of the resulting coagulopathic venom toxins. A schematic overview of the complete analytical and biochemical workflow including workflow explanation is provided in Figure 1. The combination of these experimental approaches demonstrated that small molecule inhibitor varespladib can neutralise the potent anticoagulation effects often observed from cobra venoms.

2. Results

2.1. Initial Screening of Crude Venom Coagulopathic Activity and Inhibition by Varespladib and Marimastat

Crude venom samples from elapids elicited an anticoagulant effect on bovine plasma, with venom from Naja nigricollis and N. pallida eliciting the strongest activity, and Dendroaspis polylepis eliciting the weakest activity (Figure 2 and Supplementary Figure S1). One-way ANOVA revealed significant reductions in anticoagulant activity in cobra venom with all doses of varespladib (p < 0.01), except in the case of N. pallida, which required a minimum dose of 4 µM to show significant reductions (p < 0.0001) (Figure 3, Supplementary Figures S2 and S3). This is in line with recently reported findings [39]. Varespladib caused no discernible effects on the coagulant activity of the mamba venoms. The effects of marimastat on the anticoagulant activity of elapid venoms are more difficult to disentangle. There is no significant effect of marimastat on mamba venoms, however Hemachatus haemachatus, N. haje and N. naja respond significantly to doses of 20 µM, 100 µM, and 4–20 µM, respectively (Supplementary Figures S2 and S4). The significant responses of H. haemachatus and N. naja may be statistical artifacts or errors, as it is unlikely that low doses of marimastat would stimulate a response but not higher doses; however the response of N. haje venom to 100 µM marimastat is comparable to that of varespladib (see Figure 2), suggesting that SVMPs do play a role in the anticoagulant activity of this venom, though seemingly a more minor one than that of PLA2s.

2.2. Nanofractionation, Identification, and Inhibition of Coagulopathic Mamba Venom Toxins

2.2.1. Dendroaspis polylepis

The chromatographic bioassay profile of D. polylepis venom showed a single sharp negative peak at a retention time of 16.5 min (Figure 4D). The lowest venom concentration analysed in which the anticoagulant toxin(s) was/were detected was 0.4 mg/mL (50 μL per injection, 20 μg of venom proteins), with further dilution (4 μg of venom proteins) rendering the anticoagulant activity undetectable. Noticeably, the extent of anticoagulant activity observed following nanofractionation was lower than anticipated based on findings with crude venom (Figure 4B,C). This may result from anticoagulant toxins present in D. polylepis venom (partially) denaturing during the reversed-phase LC separation. One molecular mass, at a retention time of 16.8 min, was assigned with an m/z of 1055.05 (charge +36; mol. mass 37,947.16 Da) for this active anticoagulant peak, falling within the snake venom metalloproteinases mass range (Figure 4E; Table 1). No protein was detected at the retention time of 16.8 min that corresponded to the PLA2 mass range (i.e., 13–15 kDa). So, it can be concluded that anticoagulant activity of D. polylepis was not mediated by PLA2s but by other toxins, perhaps SVMPs. Mascot searches with proteomics data retrieved from tryptic digests of nanofractionated venom toxin fractions collected in wells corresponding to the retention time of the bioactivity peak identified an SVMP toxin (species-specific venom gland transcriptome contig D_polylepis_T0167_SVMP), providing additional evidence that SVMP toxins may be at least partially responsible for the observed anticoagulant bioactivity (Figure 4; Table 1; Supplementary Figure S5).
Further, the resulting database hit shares 83.7% identity with the SVMP toxin “atragin” from Naja atra venom, which is known to exert anticoagulant activity [71]. Perhaps unsurprisingly, given the findings above, the SVMP inhibitor marimastat appeared to neutralise the detected negative anticoagulant peak at higher drug concentrations (Figure 5 and Figure 6; Supplementary Figures S12 and S13), although it is worth pointing out that we observed no evidence of crude venom neutralisation with this drug (Figure 3). However, the PLA2 inhibitor varespladib also abolished the anticoagulant toxin effects in the nanofractionation experiments, though the relatively low anticoagulation signals observed after fractionation in comparison to crude venom analysis make these results difficult to interpret, and thus observed inhibition could be the result of experimental conditions rather than direct toxin neutralisation.

2.2.2. Dendroaspis angusticeps

The chromatographic bioassay profile of D. angusticeps showed a broad negative peak at the retention time region of 16.5–18.3 min at the highest (1 mg/mL; 50 μL per injection, 50 μg of venom proteins) concentration injected (Figure 5D; Supplementary Figure S6), correlating with our crude venom experiments that revealed more substantial anticoagulant venom activity in comparison with the congener D. polylepis (Figure 2). At reduced venom concentrations, only a sharper peak at around 16.5–17.3 min remained, with further dilution resulting in a significant loss of anticoagulant activity. At the retention time of 16.5–17.3 min, there were many coeluting venom toxins detected. The most noticeable traces observed by MS were from toxins with masses in the range of 6–7 kDa (which are probably 3FTxs and/or Kunitz-type serine protease inhibitors) and from a toxin with a mass around 60 kDa (possibly an SVMP) measured with relatively low intensity (Figure 5E; Table 1). A 4 μM concentration of varespladib seemingly inhibited the anticoagulant activity of the nanofractionated venom (at 1 mg/mL concentration), while as low as a 0.16 μM concentration of marimastat exerted a similar effect (Figure 5B,C). Marimastat thus inhibited the anticoagulant effect of D. angusticeps venom at lower doses, whereas the higher doses required for varespladib inhibition hint at non-specific effects. However, as with D. polylepis, neither drug substantially inhibited crude venom anticoagulant activity, even at high doses (Figure 3). From the resulting proteomics data, D. angusticeps venom toxins detected in the bioactivity region included two 3FTxs Toxin Tx7335 (3NOJ_DENAN; rt 17.1 min), and Rho-elapitoxin-Da1b (3SI1B_DENAN; rt 17.1 min) and the Kunitz-type serine protease inhibitor long epsilon-dendrotoxin His55 (VKTHE_DENAN; rt 16.5 min) (Figure 5F,G; Table 1; Supplementary Figure S13). Other minor components detected included a venom nerve growth factor (vNGF; rt 16.7 min) and SVMPs (rt 17.1 min, 61,599.77 Da, D_angusticeps_T0082), which could at least partially explain the inhibitory profile of marimastat in the nanofractionation experiments.

2.3. Nanofractionation, Identification, and Inhibition of Coagulopathic Cobra and Rinkhals Venom Toxins

2.3.1. Naja naja

A broad negative peak defines the chromatographic bioassay profile of N. naja at a retention time region spanning 16.3–19.1 min (Figure 6C; Supplementary Figure S7). Venom dilution experiments showed a decrease of this broad negative peak into a single sharp peak with a retention time of 16.5 min, with anticoagulant activity still observed at a venom concentration of 0.04 mg/mL (50 μL per injection; 2 μg of venom proteins). The 0.2 mg/mL (50 μL per injection; 10 μg of venom proteins) venom concentration was selected for further use in inhibition experiments with varespladib. In contrast to the nanofractionation findings observed with the Dendroaspis spp., for all Naja species studied herein, marimastat failed to exert any inhibitory effect on anticoagulant venom activities at any of the tested concentrations (as was also the case for the crude venom analyses; see Figure 3), and thus will not be further discussed. Contrastingly, the anticoagulant activity of N. naja venom decreased markedly with increasing concentrations of varespladib, with complete neutralisation observed at a concentration of 0.8 µM (Figure 6B; Supplementary Figure S7). From the LC-MS data for N. naja (Figure 6D; Table 1), in the retention time region of 16.3–19.1 min, several proteins within the mass range of 6–7 kDa (i.e., 7803.47 Da; 7819.46 Da; 6750.33 Da; and 6733.32 Da) were detected, indicating the likely presence of 3FTxs. In addition, a protein with a mass of 13,321.59 Da, likely a PLA2, was also detected and eluted at the same retention time as the sharp bioactivity peak (~16.5 min). Proteomics data (Figure 6E,F; Table 1; Supplementary Figure S14) confirmed the presence of an acidic PLA2 (PA2A2_NAJNA, 13,321.59 Da) at the retention time of 16.5 min, which in combination with the varespladib inhibition data, strongly suggests that this PLA2 toxin is at least partially responsible for the observed anticoagulant venom effects.

2.3.2. Naja pallida

Similar to that of the congener N. naja, the chromatographic bioassay profile of N. pallida venom was defined by a broad negative peak at a retention time of 16.0–19.0 min (Figure 7C; Table 1; Supplementary Figure S8). Upon diluted venom analyses, the venom still retained anticoagulant properties down to a concentration of 0.04 mg/mL, though only a large peak at around 16.0 to 16.7 min remained. A 0.2 mg/mL concentration was chosen as the optimal concentration for subsequent assessment of anticoagulation inhibition by varespladib (Figure 7B; Supplementary Figure S8). From the LC-MS data (Figure 7D; Table 1), several toxins were found within the retention time region of 16.0–16.7 min and exhibited mass ranges between 6–7 kDa (6782.34 Da; 6815.32 Da), which indicated likely 3FTxs, and 13–14 kDa (i.e., 13,290.79 Da at rt 16.2 min and 13,205.74 Da and 13,320.73 Da at rt 16.4 min), indicating likely PLA2s. Although varespladib reduced the anticoagulant venom toxin activities upon increasing concentrations, with the broad peak observed in the venom only control significantly reduced by even the lowest concentration (0.032 μM) of varespladib tested, these experiments revealed a smaller anticoagulant peak that was not neutralised by varespladib even at the highest concentration tested (20 μM) (Figure 7B). From the proteomics data generated from the anticoagulant bioactivity region (Figure 7; Table 1; Supplementary Figure S15), both the basic phospholipase A2 (13290.79 Da, N_pallida_T0443) and basic phospholipase A2 nigexine (13,320.73 Da, PA2B_NAJPA) were detected.

2.3.3. Naja nigricollis

The chromatographic bioassay profile of N. nigricollis (Tanzania) can be described as a broad negative peak at a retention time region of 16.2–18.1 min (Figure 8C). The venom still exhibited anticoagulant properties down to a venom concentration of 0.08 mg/mL. Nevertheless, at lowering venom concentrations, only two large sharp negative peaks at around 16.1 and 16.8 min are left in Figure 8C. These correspond clearly to the first two sharp peaks in the LC-UV chromatogram (Figure 8A). The LC-MS data (Figure 8D; Table 1) revealed a toxin found with a mass of 13,245.77 Da at a retention time of 16.1 min, while at 16.8 min, an exact mass of 13,279.66 Da was found. Varespladib was found to effectively neutralise the anticoagulant activity of N. nigricollis venom (Figure 8B; Table 1; Supplementary Figure S9). The broad anticoagulation peak changed into two individual anticoagulation bioactivity peaks which became smaller upon increasing varespladib concentrations. A 20 μM dose of varespladib effectively inhibited all the anticoagulation activities of the venom and the proteomics data (Figure 8E,F; Table 1; Supplementary Figure S16) detected the presence of a basic PLA2 (PA2B4_NAJNG) and an acidic PLA2 (Acidic Phospholipase A2, N_nigricollis_T3199_T0053) at rt. 16.1 min and 16.8 min, respectively. These findings are identical to those we recently reported from the venom of Nigerian Naja nigricollis using the same methodology [38], where again a basic and acidic PLA2 with rts. of 16.1 and 16.8 min were responsible for the observed anticoagulant venom activity.

2.3.4. Naja haje

The chromatographic bioassay profile of N. haje shows a broad negative peak at the retention time region of 15.5–18.5 min (Figure 9C). The very broad anticoagulant peak observed indicates the potential presence of closely eluting anticoagulant toxins. However, after performing experiments with diluted venom concentrations, a sharp bioactive peak was observed at retention time 16.3–16.7 min, and anticoagulant activity was detected down to venom concentrations of 0.08 mg/mL, thus indicating lower activity compared with congeners N. pallida and N. naja. Based on the LC-MS data, a protein detected with a mass of 13,963.38 Da at a retention time of 16.5 min, suggests a PLA2 toxin is likely associated with anticoagulant activity, though some 3FTxs were also detected around the same retention time (Figure 9D; Table 1). The superimposed bioassay chromatograms and LC-UV chromatograms from the varespladib inhibition experiments are shown in Figure 9B (Supplementary Figure S10). Figure 9B shows that the observed anticoagulant activity of N. haje venom decreased with an increasing concentration of varespladib. However, the anticoagulation activity of the venom was only fully neutralised by the highest concentration (20 μM) of varespladib, despite reduced venom activity in comparison with some other cobra venoms tested. While this data might hint at the presence of other toxins contributing to the anticoagulant effects observed, the processed proteomics data revealed a toxin eluting at a retention time of 16.5 min, which was annotated as an acidic PLA2 (PA2A3_NAJSG), suggesting once again that PLA2 toxins likely underpin anticoagulant venom activity in this genus (Figure 9E,F; Table 1; Supplementary Figure S17).

2.3.5. Hemachatus haemachatus

The chromatographic bioassay profile of H. haemachatus was defined by a broad negative peak at a retention time region of 15.8–17.9 min, similar to that of related cobra species (Naja spp.) (Figure 10C). Venom dilution analyses showed anticoagulant activity was retained down to concentrations of 0.08 mg/mL, where a sharp, rather than broad, negative peak with a retention time of 16.0 min was visible. In line with data from Naja spp., the resulting mass data (Figure 10D; Table 1) detected several proteins within the mass range of 6–7 kDa (6788.43 Da, rt. 16.4 min; 6832.45 Da, rt. 16.7 min; and 6779.36 Da, rt. 17.8) which are most probably 3FTxs. A protein with a mass of 13,498.76 Da was also detected, and likely to be a PLA2. This protein was eluted with essentially the same retention time as the bioactivity peak at 15.8 min. Moreover, varespladib dose-dependently inhibited the anticoagulation activity observed in Hemachatus haemachatus venom, though only the high 20 μM concentration was capable of completely inhibiting the anticoagulation activity (Figure 10B; Supplementary Figure S11). From the proteomics data (Figure 10E,F; Table 1; Supplementary Figure S18), Hemachatus haemachatus venom mainly consisted of 3FTxs which eluted within the retention time profile of 16.4–18.2 min, though a PLA2 toxin was identified (15.8 min, PA2B1_HEMHA, basic phospholipase A2) at the exact retention time of the anticoagulation bioactivity peak.

3. Discussion

With the exception of the two mamba species studied, the elapid venoms in this study showed strong anticoagulant activity [20], particularly in the venoms of spitting cobras [25,72], but has, until now, been unverified in mambas. Both the crude coagulation assays with inhibitors and the nanofractionation results described here implicate PLA2s as the main anticoagulant agent in cobra venoms. This is unsurprising; cobra venom PLA2s have been shown to impede the process of coagulation through the inhibition of platelet aggregation [30], by non-enzymatic binding to the clotting factor Xa [25], the disruption of the tenase and prothrombinase complexes [33,34] and the lysis of fibrin.
Contrasting with the cobras, the identity of anticoagulant venom toxins in mamba venoms remains uncertain. Although inhibition of anticoagulation by varespladib was shown for nanofractionated mamba venoms at higher drug concentrations, no inhibitory effects were observed on crude mamba venoms in the plasma bioassay. The nanofractionation results are challenging to interpret due to the reduced anticoagulation levels observed compared to crude venom analyses. It may be that anticoagulant venom toxins in mamba venoms are prone to denaturation during the reversed-phase separation process. The inhibition of the remaining lower activity observed at higher drug concentrations may be a non-specific effect, thereby rendering these results inconclusive.
Further, proteomic data has previously demonstrated that mamba venoms lack any abundant PLA2s [15], which are known to be stable during reversed-phase LC separations in the proteomes of these species. Thus, other toxin families are highly likely to be responsible for the anticoagulation effects observed in mamba venoms. Our attempts to identify the toxin(s) causing the small anticoagulant peak observed in the D. polylepis bioassay chromatograms (at a retention time of 16.5 min) resulted in the detection of a bioactive venom metalloproteinase (SVMP) but the effect can also be caused by another protein that coeluted at the same retention time. Although strong matches to any elapid toxin present in the UniProt database were lacking, when searched against a D. polylepis venom gland transcriptome database [15], besides several 3FTXs eluting within the bioactivity range, one SVMP was identified from proteomic analysis, and this had a retention time consistent with the mass range of SVMPs. Despite this, the lack of inhibition of venom activity from all doses of marimastat suggests that SVMPs do not play a significant role in the anticoagulant activity in mamba venoms unless these SVMPs are unaffected by this broad-spectrum metalloproteinase-inhibiting drug. The potential roles of 3FTxs and Kunitz-type serine protease inhibitors in the observed anticoagulant activity of mamba venoms remain to be seen, and it may be worth isolating these toxins to purity followed by performing coagulation and other bioassays in the future. However, it is also possible that these toxins are coincidentally co-eluting with the actual bioactive toxin(s) in the experiments described herein.
Our findings show that the small molecule inhibitor varespladib can neutralise the potent anticoagulation effects often observed from cobra venoms. However, for mamba venoms, neither varespladib nor marimastat provided any inhibitory effects against venom-induced anticoagulation. These findings suggest that distinct toxin types are likely responsible for the effects observed in the mambas, and that anticoagulant venom effects likely evolved independently on multiple occasions in African elapids. Previous work investigated the effectiveness of antivenom on neutralising modulation of coagulation effects by snake venom using a similar analytical approach [53,63,73]. This study investigates how the small molecule inhibitors varespladib and marimastat can neutralise the same coagulant modulating effects caused by toxins in cobra and mamba venoms. Although not elaborately discussed, this research includes for all venom under study a comprehensive set of venom proteomics data. The summarised venomics data are given in this paper in the figures, while the complete detailed proteomics datasets are presented in the supporting information.
Marimastat failed to reduce the anticoagulant activity of the cobra venoms used in this study, with the exception of Naja haje, in which the highest dose (100 µM) elicited a reduction in anticoagulant activity comparable with that of varespladib at doses of as little as 0.8 µM (see Supplementary Figure S2). This conflicts with previous work performed on spitting cobras which showed a significant reduction in anticoagulant activity upon the application of marimastat and the related matrix metalloproteinase inhibitor prinomastat in most species tested [39]. The lack of inhibitory effect of marimastat observed in most species tested here could reflect a denaturation process during nanofractionation, which SVMP toxins would be subjected to, although this does not explain the lack of inhibition observed in our crude venom experiments. Thus, we conclude that, although marimastat exhibits considerable promise as a lead repurposed drug for snakebite envenoming caused by SVMP-rich viper venoms [37,74], its potential utility against elapid-induced coagulopathy appears limited, though it may be effective in cases of snakebite by N. haje in higher doses. The proportion of N. haje venom that is made up of SVMPs is debatable, ranging from 0 to 10% of the toxic component [14,75]; however, PLA2s also make up a considerably small amount of the venom; 0.1–4% [14,75], which begs the question of why varespladib is still effective at neutralising this venom at all and more so than marimastat. Marimastat has been shown to work at neutralising SVMP-mediated toxicity and lethality, especially when used in conjunction with antivenom [53,65,66,67,68,69], however this study finds no clear justification for its use as a treatment for incidents of elapid envenomation, likely due to lack of medically-important SVMPs in elapid venoms compared to those of vipers. Thus, marimastat was not included in the part of the study dealing with separated and fractionated toxins as it had not much effects on the anticoagulation evoked by the elapid venoms under study when measured as crude venoms.
Varespladib has shown considerable potential as a treatment for snakebite envenoming in previous laboratory studies, including abrogating the enzymatic effects of PLA2s from several vipers and the elapid Oxyuranus scutellatus, and Asian cobras [25,37,58,63], and can neutralise the haemotoxic effects from many elapid species [25,38,39,63]. However, the drug concentrations required to completely inhibit the anticoagulant activity of PLA2 toxins in vitro were found to vary extensively; between 0.8 μM to 20 μM in some studies [63,64]. In this study, varespladib demonstrated potent inhibition of in vitro anticoagulation activities for crude cobra venoms at doses of 4 μM and above, and for the nanofractionated toxins from these venoms. For all cobras under investigation, one to a few PLA2s could be identified as the likely bioactive toxins causing the anticoagulation activities observed, and which could all be effectively inhibited by varespladib. Varespladib is a drug candidate that has undergone phase I and II trials for the treatment of cardiovascular diseases [76,77,78] and, more recently, has acquired significant interest as a candidate drug for the treatment of snakebite, together with its orally available derivative, methyl-varespladib. Following extensive in vitro and preclinical in vivo successes, methyl-varespladib has recently entered into phase II clinical trials for snakebite indication in the USA (https://clinicaltrials.gov/ct2/show/NCT04996264, accessed on 8 October 2021). Thus, these PLA2 inhibitors are currently lead candidates for the treatment of PLA2-mediated toxicities (Clare et al., 2021), and here we provide additional evidence of their potential utility against the anticoagulant activities of the cobra, though not mamba, venoms. This supports a growing amount of evidence that varespladib is a good candidate as a snakebite therapeutic drug, alongside conventional antivenom treatment, that may aid in improving the outcomes for victims of elapid envenomation.
The question arises of why typically neurotoxic elapids have anticoagulant components in their venoms. Perhaps the anticoagulant effects of said toxins aid in spreading other neurotoxic components, such as 3FTXs, to encourage systemic envenomation. This may be evidenced by studies that show that lethal doses of venoms from the elapids Micrurus fulvius and Oxyuranus scutellatus, and neurotoxically-acting viper Vipera berus nikolskii can be neutralised by the addition of varespladib to inhibit the action of PLA2s [57,58,59,79,80]. Further, varespladib can delay the onset of lethality caused by metalloproteinase-rich venoms in vivo, again suggesting a role for PLA2 inhibition in mitigating toxin spread or dampening the severity of pathology [53]. The difficulty with this hypothesis is that PLA2s are highly flexible in function, being able to act directly to cause presynaptic neurotoxicity or inhibit coagulation or contribute to cytotoxic effects [81]. Indeed, it may be that this multifunctionality of PLA2s, which can act directly and indirectly, could be a contributing factor to the potent neurotoxicity of certain elapid snakes. However, what is clear is that African elapids have independently evolved distinct anticoagulant toxins. While these are unquestionably PLA2s in cobras (and which are particularly abundant and bioactive in spitting cobras) [20,38,72], the identity of those in mambas is much less clear, though based on findings herein, they seem likely to be toxins distinct from PLA2s [15]. The multiple independent origins of anticoagulants in these elapid snakes hint at a biologically relevant role in this activity for prey capture. Future clinical research will hopefully also contextualise the relevance of such findings in human snakebite envenomings.

4. Conclusions

This study investigated the anticoagulant activities of mamba and cobra venoms. These effects could be reduced dose-dependently with the addition of varespladib in the case of cobras and rinkhals. For all cobras, strong anticoagulation effects were observed. In case of the mambas, of the two mambas under study, Dendroaspis polylepis and Dendroaspis angusticeps, weaker anticoagulation effects were observed, of which the Dendroaspis angusticeps venom was the more potent. Neither varespladib nor marimastat neutralised the coagulation modulation evoked by mamba venom fractions tested in this study. These observations, combined with nanofractionation data identifying fractions likely to be involved in such activity, implicates the role of PLA2s from cobra, but not mamba, venoms in inhibiting blood coagulation pathways. The toxins responsible for anticoagulant effects caused by mamba venoms remain unclear. However, collectively these findings suggest that these two groups of medically important elapid snakes have independently evolved anticoagulant venom activities. Although not elaborately discussed, this research includes for all venoms under study a comprehensive set of venom proteomics data. The summarised venomics data are given in this paper in the figures, while the selected retention time detailing proteomics datasets are presented in the supporting information.

5. Materials and Methods

5.1. Chemical and Biological Reagents

Water was purified with a Milli-Q Plus system (Millipore, Amsterdam, The Netherlands). Dimethyl sulfoxide (DMSO) was supplied by Riedel-de-Haen (Zwijndrecht, The Netherlands). Acetonitrile (ACN) LC/MS grade was obtained from Concord, NC, USA. Formic acid (FA, MS grade) and acetic acid (AA) were purchased from Biosolve (Valkenswaard, The Netherlands), while CaCl2, K2HPO4, KH2PO4, NH4HCO3, phosphate buffer saline (PBS) tablets, and other salts used for buffer preparation were analytical grade and purchased from standard suppliers (Merck, Kenilworth, UK; Fluka, Bucharest, Romania; or Sigma-Aldrich, Darmstadt, Germany). The PLA2 inhibitor varespladib, the SVMP inhibitor marimastat, alkylating agent iodoacetamide, and argatroban control were purchased from Sigma-Aldrich, Darmstadt, Germany. A 0.1 mg/mL stock solution of argatroban was used for aligning the bioassay chromatogram with UV traces according to previously-defined methods [38]. The 1000 µM stock solution of varespladib and marimastat were made in DMSO and kept at −20 °C for short term storage and at −80 °C for longer-term storage. Bovine plasma was sourced from Biowest (Nuaille, France). The plasma was defrosted in a water bath at 37 °C, centrifuged, and aliquoted (approximately 10 mL aliquots in 15 mL falcon tubes) and kept at −20 °C until use. Lyophilised trypsin Gold (mass spectrometry grade) enzyme was purchased from Promega Corporation, Madison, USA. The trypsin was reconstituted in 50 mM acetic acid to obtain 1 µg/µL concentration and was aliquoted (10 µL in 1 mL Eppendorf tubes). The aliquoted solution was kept at −80 °C.
Lyophilised venom from black mamba (Dendroaspis polylepis, Tanzania), Eastern green mamba (Dendroaspis angusticeps, Tanzania), Indian cobra (Naja naja, captive-bred), red spitting cobra (Naja pallida, Tanzania), black-necked spitting cobra (Naja nigricollis, Tanzania), Egyptian cobra (Naja haje, Uganda) and rinkhals (Haemachatus haemachatus, South Africa) were sourced from snakes held in the Herpetarium of the Centre for Snakebite Research and Interventions at the Liverpool School of Tropical Medicine, UK. Venoms were stored at 4 °C until reconstituted, whereafter storage was at −80 °C.

5.2. Crude venom Coagulation Bioassaying and Inhibition with Varespladib and Marimastat

Plasma clotting assays were performed on citrated bovine plasma, following protocols previously outlined [82]. Plasma was stored at −80 °C and defrosted in a water bath at 37 °C, before being centrifuged at 2000 rpm for 4 min to remove the precipitate. The supernatant was then pipetted into a 50 mL falcon tube and kept on ice until use. For crude venom assays, venom samples were diluted with PBS, resulting in the following final assay doses: 1, 10, 50, and 100 ng. Samples were pipetted in 10 µL volumes, in replicates of 7, onto a Greiner Bio-One transparent 384-well microplate. Negative controls contained 10 µL PBS (i.e., no venom). To all wells, 20 µL of 20 mM calcium chloride was added using a Thermo Scientific Multidrop™ 384 Labsystems multidrop pipette. The plate was spun down for 5–10 s on a Benchmark Scientific PlateFuge™ MicroPlate MicroCentrifuge to remove bubbles. The multidrop was then flushed with deionised water before using it to add 20 µL of plasma into each well. Plates were read kinetically at a wavelength of 595 nm and a temperature of 25 °C using a FLUOstar Omega plate reader, with a run time set to 40 cycles at the minimum possible cycle time for the plate occupancy (141 s). Clotting activity was measured as the mean absorbance curve (area under the curve (AUC)) of each venom concentration in the relevant clotting period (determined by the clotting curve of the control to be between 24 and 71 min), standardised as a percentage of the mean AUC of the negative control in the same period. The crude clotting curves determined that 100 ng of venom in each plate well should be for the small molecule inhibition assays to follow.
To test the efficacy of small molecule inhibitors at preventing the anticoagulant activity of elapid venoms, dilutions of the SVMP inhibitor marimastat or the PLA2 inhibitor varespladib were made up using DMSO, to make final assay reaction concentrations of 0.8 µM, 4 µM, 20 µM, and 100 µM. Inhibitor-containing samples consisted of 6.25 µL venom (concentration 10 µg/mL) and 3.75 µL inhibitor. Positive controls consisted of 100 ng of crude venom and negative controls of 10 µL PBS only. All samples were pipetted into a 384-well plate in quadruplicate and the plate was incubated at 37 °C for 30 min to allow the inhibitors to take effect. As previously described, 20 µL of 20 mM calcium chloride was then added to the plate via the multidrop, centrifuged in a microcentrifuge, and then 20 µL of thawed bovine plasma was added. The plate was read in the FLUOstar at 595 nm for 45 cycles at the minimum possible cycle time for the plate occupancy (147 s). Clotting activity was measured as the mean area under the curve (AUC) of each sample type in the relevant clotting period (determined by the clotting curve of the control to be between 22 and 71 min), standardised by the mean AUC of the negative control in the same period. Readings were captured in at least triplicate. One-way analysis of variance (ANOVA) analysis was performed on clotting activity in GraphPad Prism 5.0, with Dunnett’s Multiple Comparison post-hoc tests used to determine the significance of the reduction of anticoagulant activity between the control well and the inhibitor for each species. These data were used to generate dose–response data for each drug and venom.

5.3. Coagulation of Nanofractionated Elapid Venom

5.3.1. An Overview of the Workflow

The workflow of the analysis of nanofractionated venom toxins on coagulation activity, the parallel acquired LC-MS-based accurate mass determination of the separated toxins, and the inhibition of their anticoagulation activity by varespladib and/or marimastat is represented in the lower part (ii) of Figure 1. Nanofractionation of venom (50 µL injections) was performed into 384-well plates for subsequent plasma coagulation assaying of all high-resolution collected fractions. This bioassay measures differences in coagulation velocity between control wells, which produce a typical coagulation profile, and wells containing bioactive venom toxins [38,82]. Plotting end-point measurement absorbances measured in each well against the time of fractionation results in negative peaks for anticoagulant toxins in these wells, where clotting is thus reduced or even fully inhibited. By plotting the slopes of the coagulation curves kinetically measured for each well in time against fractionation time, for procoagulant venom toxins, positive peaks are observed in these so-called bioactivity chromatograms. The parallel obtained LC-UV, and LC-MS data are plotted as superimposed chromatograms for matching UV peaks and exact masses (using the Total Ion Current, TIC, and Extracted Ion Currents, XICs) with negative bioactivity peaks observed. For analysis of the small molecule inhibitors varespladib and marimastat on their capability of neutralising elapid venom toxin-induced coagulation activities, as observed in the chromatographic coagulation bioassay profiles, each venom was nanofractionated at its optimal concentration (i.e., that concentration resulting in good bioactivity signals, but not resulting in peaks that are too broad). Nanofractionated venom toxins in the wells of the (vacuum centrifuge freeze-dried) well plates were then pre-incubated using different concentrations of varespladib (for cobras and mambas) or marimastat (for mambas), before the coagulation activity assay was performed. The resulting bioassay chromatograms of the serially diluted varespladib and marimastat experiments per venom were analysed and presented as superimposed chromatograms.

5.3.2. Instrumental Setup for HPLC/MS-Nanofractionation

LC-MS carried out the separation of venom toxins for the post-column assay with high-resolution nanofractionation. The samples (50 µL) were injected into a Shimadzu SIL-20A autosampler, and LC separation was performed with an LC system controlled via Shimadzu Lab Solution software. the gradient was set using a binary Shimadzu LC-20AB pump (pump A and B) at a total flow rate of 0.5 mL/min. Mobile phase A was water-ACN-FA (98:2:0.1, v/v/v) and mobile phase B was water-ACN-FA (2:98:0.1, v/v/v). The following gradient was used: 0% to 10% B (10 min), 10% to 95% B (20 min), 95% B (2 min), 90% to 0% B (7 min), 0% B (2 min). A 100 × 4.6 mm ID analytical column with XbridgeTM BEH300 reversed-phase C18 material (particle diameter 3.5 µm and pore size 300 Å) was used for separation. The column eluate was split in a 1:9 ratio using a low-dead-volume flow splitter. The smaller flow (0.05 mL/min) was led to a Shimadzu SPD 20A UV-Vis detector and Bruker Maxis mass spectrometer (MS). For mass analysis, the Maxis HD mass spectrometer was equipped with an electrospray ionisation source (ESI) and operated in positive ion mode. The parameters of the ESI source were source temperature 200 °C; capillary voltage 4500 V; dry gas flow 4.0 L/min; mass range 500–3000 m/z with a data-sampling time of 1 s. Protein mass (Da) were calculated using Data Analysis 5.0 (Bruker, Darmstadt, Germany). The larger eluate flow was nanofractionated (6 s/well) onto transparent 384-well plates (Greiner Bio-One, Alpen aan den Rijn, The Netherlands) controlled by Ariadne software. After nanofractionation, the plates were overnight vacuum centrifuged to dryness at room temperature using a Christ Rotational Vacuum Concentrator RVC 2–23 CD Plus (Salm en Kipp, Breukelen, The Netherlands), with a cooling trap at −80 °C. The plates were then stored at −80 °C until use.

5.3.3. Coagulation Activities of Nanofractionated Venom an Inhibition by Varespladib and Marimastat

All measurements of coagulation activity were carried out on vacuum centrifuged freeze-dried 384-well plates and performed in duplicate. The plates were first placed at room temperature for 30 min toe defrost, and assay preparation was performed at room temperature. For measurements, a working solution of 20 mM CaCl2 was prepared (2 weeks expiration time when kept at 4 °C) and frozen aliquots of citrated bovine plasma were rapidly defrosted in a water bath at 37 °C until thoroughly defrosted. After defrosting, the plasma was centrifuged (2000 RPM for 4 min) briefly to remove any particulate matter. The mixing of plasma with calcium chloride solution was carried out by first transferring 20 μL of 20 mM calcium chloride solution to the wells of the well plate followed by adding 20 µL citrated bovine plasma, using a Thermo Scientific Multidrop™ pipetting robot. The preparation and mixing process was carried out within 5 min.
To neutralise potential PLA2s and/or SVMP activities, the respective inhibitors varespladib or marimastat were added to the well plate wells prior to adding calcium chloride and plasma. For this, serial dilutions (0.16 µM; 0.8 µM; 4 µM; 20 µM; and 100 µM) of varespladib and marimastat prepared in PBS were used. With these solutions, 10 µL of inhibitor solution was added to rows 8–18 of well plates with nanofractionated venom. Additionally, PBS was used as a negative control in rows 4–7 and in rows 19–22. Plates were then spun down by centrifuging at 2000 RPM for 2 min and then incubating at room temperature for 30 min with gentle shaking (using a shaker at 60 RPM). The resulting final assay concentrations of varespladib and marimastat were 0.032 µM; 0.16 µM; 0.8 µM; 4 µM; and 20 µM. Absorbance was measured at a wavelength of 595 nm by a Thermo Scientific Varioskan Lux™ Platereader using SkanIt 4.1 software. Measurement was performed at room temperature in one kinetic loop consisting of 80 readings within 1.5 h. Two data-processing methods were carried out to produce: (1) single reading at 80th reading for anticoagulation (i.e., endpoint read), (2) slope of a reading range for procoagulation, in case any procoagulation was observed. No significant procoagulation was observed in any of the experiments conducted in this study.

5.4. Proteomics Analysis

5.4.1. In-Solution Tryptic Digestions from 384-Well Plates Containing Nanofractionated Venom Toxins

Vacuum centrifuged well-plates containing nanofractionated venom toxins were removed from the freezer and defrosted at room temperature in a closed flow cabinet. Wells containing bioactive of interest were reconstituted in 40 μL milli-Q water. Plates were then centrifuged gently at 1000 RPM for 30 sec using a well plate centrifuge (Eppendorf® Centrifuge 5810R, Nijmegen, The Netherlands), followed by careful shaking at 60 RPM for 30 min on a plate shaker (IKA® KS 4000 IC Control, Staufen, Germany). Digestion buffer was prepared by mixing a reducing agent (0.5% β-mercaptoethanol in milli-Q water) and freshly prepared AmmBi buffer pH 8.2 (25 mM ammonium bicarbonate) at a 1:9 ratio. Thereafter, 50 μL of digestion buffer and 30 μL of reconstituted venom proteins were transferred to an Eppendorf tube, followed by vortexing at 1000 RPM for 10 sec using a tube centrifuge (Prism® mini labnet). The samples were then reduced at 95 °C for 10 min using a dry block heating thermostat (Biosan® Bio TDB-100) for incubation. Next, 9 μL of alkylation agent (100 mM Iodoacetamide) was added followed by incubation for 30 min at room temperature in a closed cabinet to prevent light degradation of the alkylating agent. Following incubation, trypsin stock solution (1 μg/μL in 50 mM acetic acid solution) was diluted by AmmBi buffer pH 8.2 to obtain a 0.1 μg/μL concentration. Then, 5 μL of trypsin solution was added to each Eppendorf tube, vortexed, and subsequently incubated at 37 °C for 1.5–3.0 h using an incubator (Memmert, Schwabach, Germany). After that, an additional 5 μL of trypsin solution was added to the Eppendorf tubes and the samples were incubated at 37 °C overnight inside the incubator. After this incubation step, 10 μL of 5% formic acid solution was added to each tube to quench the digestion process. The resulting samples were finally transferred to autosampler vials and analysed by nanoLC–ESI-MS/MS.

5.4.2. Instrument Setup for Proteomics

An Ultimate 3000 nano HPLC module (Thermo Scientific, Waltham, MA, USA) coupled to a Bruker TIMS-TOF Mass Spectrometer (Bruker Daltonics, Bremen, Germany) was used. Sample (1 µL) was injected with a WPS-3000(RS) autosampler, and nanoLC separation was performed with a nanoLC system controlled via Chromeleon 7.2 SR4 MUb software. The gradient was set using a nanoLC binary pump (A and B) at a total flow rate of 0.5 µL/min. Mobile phase A was water-FA (100:0.1, v/v), and mobile phase B was water-ACN-FA (20:80:0.05, v/v/v). The system was also equipped with a loading pump, for which solvent water-ACN-FA (99:1:0.05, v/v/v) was used. The following gradient was used: 1% B (10 min), 1–20% B (5 min), 20–50% B (30 min), 50–85% B (1 min), 85% B (5 min), 85–1% B (0.5 min), and 1% B (9.5 min). For sample trapping, an Acclaim PepMap 100 reversed-phase C18 trapping column (particle diameter 5 µm; pore size 100 Å; and column dimensions of 5 × 0.3 mm) was used. An Acclaim PepMap 100 reversed-phase C18 analytical column (particle diameter 2 µm; pore size 100 Å; and column dimensions of 150 × 0.75 mm), was used to separate the peptides in the samples subsequently. Both the analytical column and the trapping column were placed in a column oven, of which the temperature was set at 45 °C. The column eluate was transferred to a Bruker TIMS-TOF Mass Spectrometer (Bruker Daltonics, Bremen, Germany) equipped with a captive spray ionisation source in positive ion mode. The parameters of the CSI source were source temperature, 150 °C; desolvation temperature, 180 °C; capillary voltage, 1300 V; gas flow, 3 L/min. The monitored mass range was m/z 300–3000 with a data-sampling time of 0.5 s. The collision energy was 10 eV with prepulse storage of 10 µs. For MS/MS, the following parameters were used: precursor ion list 300–3000 m/z; number of precursor ions 3; threshold 250; preferred charge states +1 up to +3; fixed MS/MS acquisition 2 Hz (0.5 s).

5.4.3. Proteomics Data Interpretation and Analysis

The tryptic digests were analysed with nanoLC–MS/MS, and the obtained MS/MS data for each of the samples were converted into mascot generic format files (MGF) by using Bruker DataAnalysis software. The obtained MGF files were subjected to Mascot database searching by using two different databases: the Swiss-Prot database and a species-specific database (MassIVE accession numbers MSV000081885 and MSV000080491, NCBI sequence read archive numbers SRX5026268, SRX7700833, SRX7700829, SRX5026269, SRX2768426, and SRX2768427). A detailed method for performing the in-solution tryptic digests is described in Section 5.4.1. The following parameters were used for the Mascot searches: (1) since iodoacetamide was used as an alkylating agent, fixed modification: carbamidomethyl (C) was chosen (adding 34 Da to methionine residues); (2) variable modifications: amidation (C-term) and methionine oxidation (M), (3) peptide tolerance of ±0.1% and MS/MS tolerance of ±0.05 Da, (4) peptide charges of +1, +2, and +3 allowed.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/toxins14110736/s1, Supplementary Data contains: (a) raw mass spectrometry data of seven elapid venom species. The raw data were generated from Bruker mass spectrometer and can be accessed using Compass Data Analysis 5.0 software. (b) The Mascot Generic Files (MGF) were created from deconvoluting the raw mass spectrometry data using nanoLC–MS/MS. (c) an LC-UV raw chromatogram data that was generated using Shimadzu Lab Solution software. Supplementary Figures consist of figures that are not shown in the main text: Supplementary Figure S1. Elapid venoms demonstrate varying degrees of anticoagulant activity. Supplementary Figure S2. The PLA2 inhibitor varespladib, but not the SVMP inhibitor marimastat, cause significant reductions in anticoagulant activity induced by cobra venom. Supplementary Figure S3. The anticoagulant curves of venom with the addition of the PLA2 inhibitor varespladib show a reduction in the anticoagulant activity of cobra venoms with increasing doses of varespladib. Supplementary Figure S4. The coagulation curves of venom with the addition of the SVMP inhibitor marimastat show no reduction in the anticoagulant activity of cobra venoms but do inhibit the procoagulant toxins of Hemachatus haemachatus, causing the venom to become anticoagulant. Supplementary Figures S5–S11 Duplication results of superimposed pro- and anticoagulation bioassay chromatogram of Dendroaspis angusticeps, D. angusticeps, Naja naja, Naja pallida, Naja nigricollis, Naja haje and, Hemachatus haemachatus venom respectively. Figures S12–S18 Detailed results of proteomics searches using Mascot software against the Species-specific database and Swiss-Prot database to identify anticoagulant venom proteins of Dendroaspis angusticeps, D. angusticeps, Naja naja, Naja pallida, Naja nigricollis, Naja haje and, Hemachatus haemachatus venom respectively.

Author Contributions

Conceptualisation, J.K., G.W.S., N.R.C., and A.A.; methodology, A.A., T.D.K., K.B.M.S., and J.S.; software, A.A. and K.B.M.S.; validation, J.K., G.W.S., and N.R.C.; writing—original draft preparation, A.A., T.D.K., and writing—review and editing, J.K., F.J.V., K.B.M.S., N.R.C., and G.W.S.; visualisation, A.A. and T.D.K.; supervision, J.K., N.R.C., and funding acquisition, A.A. and N.R.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by an Indonesian Endowment Fund for Education (LPDP) Ph.D. research studentship [S-5622/LPDP.4/2020] to A.A. and a UK Medical Research Council grant (MR/S00016X/1) and a Wellcome Trust and Royal Society Sir Henry Dale Fellowship (200517/Z/16/Z) to N.R.C.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The underpinning raw MS data of seven elapid venom species, MGF files, and LC-UV raw chromatogram can be found in FigShare using the following links. MGF Files of cobra venom DOI: doi.org/10.6084/m9.figshare.19397903. MGF Files of mamba venom DOI: doi.org/10.6084/m9.figshare.19397852. LC-UV raw chromatogram DOI: doi.org/10.6084/m9.figshare.19447583.

Acknowledgments

The authors thank Paul Rowley and Edouard Crittenden for maintenance of snakes and provision of venom.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kasturiratne, A.; Wickremasinghe, A.R.; De Silva, N.; Gunawardena, N.K.; Pathmeswaran, A.; Premaratna, R.; Savioli, L.; Lalloo, D.G.; De Silva, H.J. The Global Burden of Snakebite: A Literature Analysis and Modelling Based on Regional Estimates of Envenoming and Deaths. PLoS Med. 2008, 5, 1591–1604. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Gutiérrez, J.M.; Calvete, J.J.; Habib, A.G.; Harrison, R.A.; Williams, D.J.; Warrell, D.A. Snakebite Envenoming. Nat. Rev. Dis. Primers 2017, 3, 17063. [Google Scholar] [CrossRef] [PubMed]
  3. Fox, S.; Rathuwithana, A.C.; Kasturiratne, A.; Lalloo, D.G.; de Silva, H.J. Underestimation of Snakebite Mortality by Hospital Statistics in the Monaragala District of Sri Lanka. Trans. R. Soc. Trop. Med. Hyg. 2006, 100, 693–695. [Google Scholar] [CrossRef]
  4. Harrison, R.A.; Hargreaves, A.; Wagstaff, S.C.; Faragher, B.; Lalloo, D.G. Snake Envenoming: A Disease of Poverty. PLoS Negl. Trop. Dis. 2009, 3, e569. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Longbottom, J.; Shearer, F.M.; Devine, M.; Alcoba, G.; Chappuis, F.; Weiss, D.J.; Ray, S.E.; Ray, N.; Warrell, D.A.; Bill, F.; et al. Vulnerability to Snakebite Envenoming: A Global Mapping of Hotspots. Lancet 2018, 6736, 673–684. [Google Scholar] [CrossRef] [Green Version]
  6. Tasoulis, T.; Isbister, G.K. A Review and Database of Snake Venom Proteomes. Toxins 2017, 9, 290. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Aird, S.D. Ophidian Envenomation Strategies and the Role of Purines. Toxicon 2002, 40, 335–393. [Google Scholar] [CrossRef]
  8. Kumar, T.K.S.; Pandian, S.K.; Srisailam, S.; Yu, C. Structure and Function of Snake Venom Cardiotoxins. J. Toxicol. Toxin Rev. 1998, 17, 183–211. [Google Scholar] [CrossRef]
  9. Thwin, M.M.; Gopalakrishnakone, P. Snake Envenomation and Protective Natural Endogenous Proteins: A Mini Review of the Recent Developments (1991–1997). Toxicon 1998, 36, 1471–1482. [Google Scholar] [CrossRef]
  10. Barua, A.; Mikheyev, A.S.; Russo, C. Many Options, Few Solutions: Over 60 My Snakes Converged on a Few Optimal Venom Formulations. Mol. Biol. Evol. 2019, 36, 1964–1974. [Google Scholar] [CrossRef]
  11. Ponnappa, K.C.; Saviour, P.; Ramachandra, N.B.; Kini, R.M.; Gowda, T.V. INN-Toxin, a Highly Lethal Peptide from the Venom of Indian Cobra (Naja Naja) Venom-Isolation, Characterization and Pharmacological Actions. Peptides 2008, 29, 1893–1900. [Google Scholar] [CrossRef] [PubMed]
  12. Warrell, D.A. Clinical Toxicology of Snakebite in Asia. In Handbook of Clinical Toxicology of Animal Venoms and Poisons; Meier, J., White, J., Eds.; Informa Healthcare USA Inc: New York, NY, USA, 2008; pp. 496–588. ISBN 978-0-8493-4489-3. [Google Scholar]
  13. Panagides, N.; Jackson, T.; Ikonomopoulou, M.; Arbuckle, K.; Pretzler, R.; Yang, D.; Ali, S.; Koludarov, I.; Dobson, J.; Sanker, B.; et al. How the Cobra Got Its Flesh-Eating Venom: Cytotoxicity as a Defensive Innovation and Its Co-Evolution with Hooding, Aposematic Marking, and Spitting. Toxins 2017, 9, 103. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Kazandjian, T.D.; Robinson, S.D.; Greene, H.W.; Carter, D.A.; Wouters, R.M.; Wagstaff, S.C.; Arias, A.S.; Albulescu, L.-O.; McCabe, C.V.; da Silva, R.R.; et al. Convergent Evolution of Pain-Inducing Defensive Venom Components in Spitting Cobras. bioRxiv 2020. [Google Scholar] [CrossRef]
  15. Ainsworth, S.; Petras, D.; Engmark, M.; Süssmuth, R.D.; Whiteley, G.; Albulescu, L.O.; Kazandjian, T.D.; Wagstaff, S.C.; Rowley, P.; Wüster, W.; et al. The Medical Threat of Mamba Envenoming in Sub-Saharan Africa Revealed by Genus-Wide Analysis of Venom Composition, Toxicity and Antivenomics Profiling of Available Antivenoms. J. Proteom. 2018, 172, 173–189. [Google Scholar] [CrossRef] [PubMed]
  16. Harvey, A.L. Twenty Years of Dendrotoxins. Toxicon 2001, 39, 15–26. [Google Scholar] [CrossRef]
  17. Harvey, A.L.; Robertson, B. Dendrotoxins: Structure-Activity Relationships and Effects on Potassium Ion Channels. Curr. Med. Chem. 2004, 11, 3065–3072. [Google Scholar] [CrossRef]
  18. Marshall, L.R.; Herrmann, R.P. Coagulant and Anticoagulant Actions of Australian Snake Venoms. Thromb. Haemost. 1983, 50, 707–711. [Google Scholar] [CrossRef]
  19. Lister, C.; Arbuckle, K.; Jackson, T.N.W.; Debono, J.; Zdenek, C.N.; Dashevsky, D.; Dunstan, N.; Allen, L.; Hay, C.; Bush, B.; et al. Catch a Tiger Snake by Its Tail: Differential Toxicity, Co-Factor Dependence and Antivenom Efficacy in a Procoagulant Clade of Australian Venomous Snakes. Comp. Biochem. Physiol. Part C. Toxicol. Pharmacol. 2017, 202, 39–54. [Google Scholar] [CrossRef] [Green Version]
  20. Bittenbinder, M.A.; Dobson, J.S.; Zdenek, C.N.; op den Brouw, B.; Naude, A.; Vonk, F.J.; Fry, B.G. Differential Destructive (Non-Clotting) Fibrinogenolytic Activity in Afro-Asian Elapid Snake Venoms and the Links to Defensive Hooding Behavior. Toxicol. In Vitro 2019, 60, 330–335. [Google Scholar] [CrossRef]
  21. Youngman, N.J.; Zdenek, C.N.; Dobson, J.S.; Bittenbinder, M.A.; Gillett, A.; Hamilton, B.; Dunstan, N.; Allen, L.; Veary, A.; Veary, E.; et al. Mud in the Blood: Novel Potent Anticoagulant Coagulotoxicity in the Venoms of the Australian Elapid Snake Genus Denisonia (Mud Adders) and Relative Antivenom Efficacy. Toxicol. Lett. 2019, 302, 1–6. [Google Scholar] [CrossRef]
  22. Slagboom, J.; Kool, J.; Harrison, R.A.; Casewell, N.R. Haemotoxic Snake Venoms: Their Functional Activity, Impact on Snakebite Victims and Pharmaceutical Promise. Br. J. Haematol. 2017, 177, 947–959. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. McCleary, R.J.R.; Kini, R.M. Snake Bites and Hemostasis/Thrombosis. Thromb. Res. 2013, 132, 642–646. [Google Scholar] [CrossRef] [PubMed]
  24. Ali, S.A.; Jackson, T.N.W.; Casewell, N.R.; Low, D.H.W.; Rossi, S.; Baumann, K.; Fathinia, B.; Visser, J.; Nouwens, A.; Hendrikx, I.; et al. Extreme Venom Variation in Middle Eastern Vipers: A Proteomics Comparison of Eristicophis Macmahonii, Pseudocerastes Fieldi and Pseudocerastes Persicus. J. Proteom. 2015, 116, 106–113. [Google Scholar] [CrossRef]
  25. Bittenbinder, M.A.; Zdenek, C.N.; Op Den Brouw, B.; Youngman, N.J.; Dobson, J.S.; Naude, A.; Vonk, F.J.; Fry, B.G. Coagulotoxic Cobras: Clinical Implications of Strong Anticoagulant Actions of African Spitting Naja Venoms That Are Not Neutralised by Antivenom but Are by LY315920 (Varespladib). Toxins 2018, 10, 516. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Mantok, S.; Sekhon, H.; Sahi, G.; Jones, P. Cronfa-Swansea University Open Access Repository. J. Small Bus. Enterp. Dev. 2019, 92–105. [Google Scholar]
  27. Seneci, L.; Zdenek, C.N.; Bourke, L.A.; Cochran, C.; Elda, E.; Neri-castro, E.; Frank, N.; Fry, B.G.; Lucia, S.; Sciences, B.; et al. A symphony of destruction: Dynamic differential fibrinogenolytic toxicity by rattlesnake (Crotalus and Sistrurus) venoms. Comp. Biochem. Physiol. Part C Toxicol. Pharmacol. 2021, 245, 109034. [Google Scholar] [CrossRef] [PubMed]
  28. Asad, M.H.H.B.; Razi, M.T.; Khan, T.; Najam-Us-Saqib, Q.; Murtaza, G.; Hussain, M.S.; Hussain, M.S.; Karim, S.; Hussain, I. Coagulopathies in Naja Naja Karachiensis (Black Pakistan Cobra) Bites and Its Effect on Coagulation Tests upon Storage of Platelet-Poor Plasma. Acta Pol. Pharm.-Drug Res. 2012, 69, 1031–1034. [Google Scholar]
  29. MacKay, N.; Ferguson, J.C.; McNicol, G.P. Effects of Three Cobra Venoms on Blood Coagulation, Platelet Aggregation, and Fibrinolysis. J. Clin. Pathol. 1969, 22, 304–311. [Google Scholar] [CrossRef] [Green Version]
  30. Sundell, I.B.; Rånby, M.; Zuzel, M.; Robinson, K.A.; Theakston, R.D.G. In Vitro Procoagulant and Anticoagulant Properties of Naja Naja Naja Venom. Toxicon 2003, 42, 239–247. [Google Scholar] [CrossRef]
  31. Osipov, A.V.; Filkin, S.Y.; Makarova, Y.V.; Tsetlin, V.I.; Utkin, Y.N. A New Type of Thrombin Inhibitor, Noncytotoxic Phospholipase A2, from the Naja Haje Cobra Venom. Toxicon 2010, 55, 186–194. [Google Scholar] [CrossRef]
  32. Dutta, S.; Gogoi, D.; Mukherjee, A.K. Anticoagulant Mechanism and Platelet Deaggregation Property of a Non-Cytotoxic, Acidic Phospholipase A2 Purified from Indian Cobra (Naja Naja) Venom: Inhibition of Anticoagulant Activity by Low Molecular Weight Heparin. Biochimie 2015, 110, 93–106. [Google Scholar] [CrossRef]
  33. Kini, R.M. Structure-Function Relationships and Mechanism of Anticoagulant Phospholipase A2 Enzymes from Snake Venoms. Toxicon 2005, 45, 1147–1161. [Google Scholar] [CrossRef] [PubMed]
  34. Stefansson, S.; Kini, R.M.; Evans, H.J. The Inhibition of Clotting Complexes from the Extrinsic Coagulation Cascade by the Phospholipase A2 Isoenzymes from Naja Nigricollis Venom. Thromb. Res. 1989, 55, 481–491. [Google Scholar] [CrossRef]
  35. Boffa, M.C.; Barbier, D.; de Angulo, M. Anticoagulant Effect of Cardiotoxins. Thromb. Res. 1983, 32, 635–640. [Google Scholar] [CrossRef]
  36. Jiang, M.S.; Fletcher, J.E.; Smith, L.A. Factors Influencing the Hemolysis of Human Erythrocytes by Cardiotoxins from Naja Naja Kaouthia and Naja Naja Atra Venoms and a Phospholipase A2with Cardiotoxin-like Activities from Bungarus Fasciatus Venom. Toxicon 1989, 27, 247–257. [Google Scholar] [CrossRef]
  37. Xie, C.; Bittenbinder, M.A.; Slagboom, J.; Arrahman, A.; Kool, J. Erythrocyte Haemotoxicity Profiling of Snake Venom Toxins after Nanofractionation. J. Chromatogr. B 2021, 1176, 122586. [Google Scholar] [CrossRef]
  38. Kazandjian, T.; Arrahman, A.; Still, K.B.M.; Somsen, G.W.; Vonk, F.J.; Nicholas, R.; Wilkinson, M.C.; Kool, J. Anticoagulant Activity of Naja Nigricollis Venom Is Mediated by Phospholipase A2 Toxins and Inhibited by Varespladib. Toxins 2021, 13, 302. [Google Scholar] [CrossRef]
  39. Chowdhury, A.; Lewin, M.R.; Zdenek, C.N.; Carter, R.; Fry, B.G. The Relative Efficacy of Chemically Diverse Small-Molecule Enzyme-Inhibitors Against Anticoagulant Activities of African Spitting Cobra (Naja Species) Venoms. Front. Immunol. 2021, 12, 4215. [Google Scholar] [CrossRef]
  40. MacKay, N.; Ferguson, J.C.; McNicol, G.P. Effects of Three Mamba Venoms on the Haemostatic Mechanism. Br. J. Haematol. 1968, 15, 549–560. [Google Scholar] [CrossRef]
  41. Kaur, N.; Iyer, A.; Sunagar, K. Evolution Bites—Timeworn Inefficacious Snakebite Therapy in the Era of Recombinant Vaccines. Indian. Pediatr. 2021, 58, 219–223. [Google Scholar] [CrossRef]
  42. Attarde, S.; Iyer, A.; Khochare, S.; Shaligram, U.; Vikharankar, M.; Sunagar, K. The Preclinical Evaluation of a Second-Generation Antivenom for Treating Snake Envenoming in India. Toxins 2022, 14, 168. [Google Scholar] [CrossRef] [PubMed]
  43. Harrison, R.A.; Oluoch, G.O.; Ainsworth, S.; Alsolaiss, J.; Bolton, F.; Arias, A.S.; Gutiérrez, J.M.; Rowley, P.; Kalya, S.; Ozwara, H.; et al. Preclinical Antivenom-Efficacy Testing Reveals Potentially Disturbing Deficiencies of Snakebite Treatment Capability in East Africa. PLoS Negl. Trop. Dis. 2017, 11, e0005969. [Google Scholar] [CrossRef] [PubMed]
  44. Menzies, S.K.; Litschka-Koen, T.; Edge, R.J.; Alsolaiss, J.; Crittenden, E.; Hall, S.R.; Westhorpe, A.; Thomas, B.; Murray, J.; Shongwe, N.; et al. Two Snakebite Antivenoms Have Potential to Reduce Eswatini’s Dependency upon a Single, Increasingly Unavailable Product: Results of Preclinical Efficacy Testing. PLoS Negl. Trop. Dis. 2022, 16, e0010496. [Google Scholar] [CrossRef] [PubMed]
  45. Ainsworth, S.; Menzies, S.K.; Casewell, N.R.; Harrison, R.A. An Analysis of Preclinical Efficacy Testing of Antivenoms for Sub-Saharan Africa: Inadequate Independent Scrutiny and Poor-Quality Reporting Are Barriers to Improving Snakebite Treatment and Management. PLoS Negl. Trop. Dis. 2020, 14, e0008579. [Google Scholar] [CrossRef]
  46. Kalogeropoulos, K.; Treschow, A.F.; Auf Dem Keller, U.; Escalante, T.; Rucavado, A.; Gutiérrez, J.M.; Laustsen, A.H.; Workman, C.T. Protease Activity Profiling of Snake Venoms Using High-Throughput Peptide Screening. Toxins 2019, 11, 170. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Harrison, R.A.; Cook, D.A.; Renjifo, C.; Casewell, N.R.; Currier, R.B.; Wagstaff, S.C. Research Strategies to Improve Snakebite Treatment: Challenges and Progress. J. Proteom. 2011, 74, 1768–1780. [Google Scholar] [CrossRef]
  48. Kini, R.M.; Sidhu, S.S.; Laustsen, A.H. Biosynthetic Oligoclonal Antivenom (BOA) for Snakebite and next-Generation Treatments for Snakebite Victims. Toxins 2018, 10, 534. [Google Scholar] [CrossRef] [Green Version]
  49. Habib, A.G.; Brown, N.I. The Snakebite Problem and Antivenom Crisis from a Health-Economic Perspective. Toxicon 2018, 150, 115–123. [Google Scholar] [CrossRef]
  50. Gutiérrez, J.M. Global Availability of Antivenoms: The Relevance of Public Manufacturing Laboratories. Toxins 2019, 11, 5. [Google Scholar] [CrossRef] [Green Version]
  51. Habibid, A.G.; Musa, B.M.; Iliyasuid, G.; Hamza, M.; Kuznik, A.; Chippauxid, J.P. Challenges and Prospects of Snake Antivenom Supply in Sub-Saharan Africa. PLoS Negl. Trop. Dis. 2020, 14, e0008374. [Google Scholar] [CrossRef]
  52. Bulfone, T.C.; Samuel, S.P.; Bickler, P.E.; Lewin, M.R. Developing Small Molecule Therapeutics for the Initial and Adjunctive Treatment of Snakebite. J. Trop. Med. 2018, 2018, 4320175. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Albulescu, L.O.; Xie, C.; Ainsworth, S.; Alsolaiss, J.; Crittenden, E.; Dawson, C.A.; Softley, R.; Bartlett, K.E.; Harrison, R.A.; Kool, J.; et al. A Therapeutic Combination of Two Small Molecule Toxin Inhibitors Provides Broad Preclinical Efficacy against Viper Snakebite. Nat. Commun. 2020, 11, 6094. [Google Scholar] [CrossRef] [PubMed]
  54. Varespladib. Am. J. Cardiovasc. Drugs 2011, 11, 137–143. [CrossRef]
  55. Thomas, A.L.; Steward, W.P. Marimastat: The Clinical Development of a Matrix Metalloproteinase Inhibitor. Expert Opin. Investig. Drugs 2000, 9, 2913–2922. [Google Scholar] [CrossRef]
  56. Fontana Oliveira, I.C.; Gutiérrez, J.M.; Lewin, M.R.; Oshima-Franco, Y. Varespladib (LY315920) Inhibits Neuromuscular Blockade Induced by Oxyuranus Scutellatus Venom in a Nerve-Muscle Preparation. Toxicon 2020, 187, 101–104. [Google Scholar] [CrossRef] [PubMed]
  57. Gutiérrez, J.M.; Lewin, M.R.; Williams, D.J.; Lomonte, B. Varespladib (LY315920) and Methyl Varespladib (LY333013) Abrogate or Delay Lethality Induced by Presynaptically Acting Neurotoxic Snake Venoms. Toxins 2020, 12, 131. [Google Scholar] [CrossRef] [Green Version]
  58. Lewin, M.R.; Samue, S.P.; Merkel, J.; Bickler, P. Varespladib (LY315920) Appears to Be a Very Potent, Broad-Spectrum, Inhibitor of Snake Venom PLA2s from Six Continents. Toxicon 2016, 117, 103–104. [Google Scholar] [CrossRef]
  59. Zinenko, O.; Tovstukha, I.; Korniyenko, Y. PLA2 Inhibitor Varespladib as an Alternative to the Antivenom Treatment for Bites from Nikolsky’s Viper Vipera Berus Nikolskii. Toxins 2020, 12, 356. [Google Scholar] [CrossRef]
  60. Wang, Y.; Zhang, J.; Zhang, D.; Xiao, H.; Xiong, S.; Huang, C. Exploration of the Inhibitory Potential of Varespladib for Snakebite Envenomation. Molecules 2018, 23, 391. [Google Scholar] [CrossRef] [Green Version]
  61. Lewin, M.; Samuel, S.; Merkel, J.; Bickler, P. Varespladib (LY315920) Appears to Be a Potent, Broad-Spectrum, Inhibitor of Snake Venom Phospholipase A2 and a Possible Pre-Referral Treatment for Envenomation. Toxins 2016, 8, 248. [Google Scholar] [CrossRef] [Green Version]
  62. Silva-Carvalho, R.; Gaspar, M.Z.; Quadros, L.H.B.; Lobo, L.G.G.; Giuffrida, R.; Santarém, C.L.; Silva, E.O.; Gerez, J.R.; Silva, N.J.; Hyslop, S.; et al. Partial Efficacy of a Brazilian Coralsnake Antivenom and Varespladib in Neutralizing Distinct Toxic Effects Induced by Sublethal Micrurus Dumerilii Carinicauda Envenoming in Rats. Toxicon 2022, 213, 99–104. [Google Scholar] [CrossRef] [PubMed]
  63. Xie, C.; Albulescu, L.-O.; Still, K.; Slagboom, J.; Zhao, Y.; Jiang, Z.; Somsen, G.; Vonk, F.; Casewell, N.; Kool, J. Varespladib Inhibits the Phospholipase A2 and Coagulopathic Activities of Venom Components from Haemotoxic Snakes. Biomedicines 2020, 8, 165. [Google Scholar] [CrossRef] [PubMed]
  64. Youngman, N.J.; Walker, A.; Naude, A.; Coster, K.; Sundman, E.; Fry, B.G. Varespladib (LY315920) Neutralises Phospholipase A2 Mediated Prothrombinase-Inhibition Induced by Bitis Snake Venoms. Comp. Biochem. Physiol. Part-C. Toxicol. Pharmacol. 2020, 236, 108818. [Google Scholar] [CrossRef] [PubMed]
  65. Howes, J.M.; Theakston, R.D.G.; Laing, G.D. Neutralization of the Haemorrhagic Activities of Viperine Snake Venoms and Venom Metalloproteinases Using Synthetic Peptide Inhibitors and Chelators. Toxicon 2007, 49, 734–739. [Google Scholar] [CrossRef] [PubMed]
  66. Slagboom, J.; Mladić, M.; Xie, C.; Kazandjian, T.D.; Vonk, F.; Somsen, G.W.; Casewell, N.R.; Kool, J. High Throughput Screening and Identification of Coagulopathic Snake Venom Proteins and Peptides Using Nanofractionation and Proteomics Approaches. PLoS Negl. Trop. Dis. 2020, 14, e0007802. [Google Scholar] [CrossRef] [Green Version]
  67. Layfield, H.J.; Williams, H.F.; Ravishankar, D.; Mehmi, A.; Sonavane, M.; Salim, A.; Vaiyapuri, R.; Lakshminarayanan, K.; Vallance, T.M.; Bicknell, A.B.; et al. Repurposing Cancer Drugs Batimastat and Marimastat to Inhibit the Activity of a Group i Metalloprotease from the Venom of the Western Diamondback Rattlesnake, Crotalus Atrox. Toxins 2020, 12, 309. [Google Scholar] [CrossRef]
  68. Menzies, S.K.; Clare, R.H.; Xie, C.; Westhorpe, A.; Hall, S.R.; Edge, R.J.; Alsolaiss, J.; Crittenden, E.; Marriott, A.E.; Harrison, R.A.; et al. In Vitro and in Vivo Preclinical Venom Inhibition Assays Identify Metalloproteinase Inhibiting Drugs as Potential Future Treatments for Snakebite Envenoming by Dispholidus Typus. Toxicon X 2022, 14, 100118. [Google Scholar] [CrossRef]
  69. Arias, A.S.; Rucavado, A.; Gutiérrez, J.M. Peptidomimetic Hydroxamate Metalloproteinase Inhibitors Abrogate Local and Systemic Toxicity Induced by Echis Ocellatus (Saw-Scaled) Snake Venom. Toxicon 2017, 132, 40–49. [Google Scholar] [CrossRef]
  70. Shan, L.L.; Gao, J.F.; Zhang, Y.X.; Shen, S.S.; He, Y.; Wang, J.; Ma, X.M.; Ji, X. Proteomic Characterization and Comparison of Venoms from Two Elapid Snakes (Bungarus Multicinctus and Naja Atra) from China. J. Proteom. 2016, 138, 83–94. [Google Scholar] [CrossRef]
  71. Kini, R.M. Anticoagulant Proteins from Snake Venoms: Structure, Function and Mechanism. Biochem. J. 2006, 397, 377–387. [Google Scholar] [CrossRef] [Green Version]
  72. Choudhury, M.; McCleary, R.J.R.; Kesherwani, M.; Kini, R.M.; Velmurugan, D. Comparison of Proteomic Profiles of the Venoms of Two of the ‘Big Four’ Snakes of India, the Indian Cobra (Naja Naja) and the Common Krait (Bungarus Caeruleus), and Analyses of Their Toxins. Toxicon 2017, 135, 33–42. [Google Scholar] [CrossRef] [PubMed]
  73. Xie, C.; Slagboom, J.; Albulescu, L.O.; Bruyneel, B.; Still, K.B.M.; Vonk, F.J.; Somsen, G.W.; Casewell, N.R.; Kool, J. Antivenom Neutralization of Coagulopathic Snake Venom Toxins Assessed by Bioactivity Profiling Using Nanofractionation Analytics. Toxins 2020, 12, 53. [Google Scholar] [CrossRef] [PubMed]
  74. Albulescu, L.O.; Hale, M.S.; Ainsworth, S.; Alsolaiss, J.; Crittenden, E.; Calvete, J.J.; Evans, C.; Wilkinson, M.C.; Harrison, R.A.; Kool, J.; et al. Preclinical Validation of a Repurposed Metal Chelator as an Early-Intervention Therapeutic for Hemotoxic Snakebite. Sci. Transl. Med. 2020, 12, eaay8314. [Google Scholar] [CrossRef] [PubMed]
  75. Malih, I.; Ahmad rusmili, M.R.; Tee, T.Y.; Saile, R.; Ghalim, N.; Othman, I. Proteomic Analysis of Moroccan Cobra Naja Haje Legionis Venom Using Tandem Mass Spectrometry. J. Proteom. 2014, 96, 240–252. [Google Scholar] [CrossRef] [PubMed]
  76. Rosenson, R.S.; Elliott, M.; Stasiv, Y.; Hislop, C. Randomized Trial of an Inhibitor of Secretory Phospholipase A2 on Atherogenic Lipoprotein Subclasses in Statin-Treated Patients with Coronary Heart Disease. Eur. Heart J. 2011, 32, 999–1005. [Google Scholar] [CrossRef] [Green Version]
  77. Nicholls, S.J.; Cavender, M.A.; Kastelein, J.J.P.; Schwartz, G.; Waters, D.D.; Rosenson, R.S.; Bash, D.; Hislop, C. Inhibition of Secretory Phospholipase A 2 in Patients with Acute Coronary Syndromes: Rationale and Design of the Vascular Inflammation Suppression to Treat Acute Coronary Syndrome for 16 Weeks (VISTA-16) Trial. Cardiovasc. Drugs Ther. 2012, 26, 71–75. [Google Scholar] [CrossRef]
  78. Nicholls, S.J.; Kastelein, J.J.P.; Schwartz, G.G.; Bash, D.; Rosenson, R.S.; Cavender, M.A.; Brennan, D.M.; Koenig, W.; Jukema, J.W.; Nambi, V.; et al. Varespladib and Cardiovascular Events in Patients with an Acute Coronary Syndrome: The VISTA-16 Randomized Clinical Trial. JAMA 2014, 311, 252–262. [Google Scholar] [CrossRef] [Green Version]
  79. Lewin, M.R.; Gilliam, L.L.; Gilliam, J.; Samuel, S.P.; Bulfone, T.C.; Bickler, P.E.; Gutiérrez, J.M. Delayed LY333013 (Oral) and LY315920 (Intravenous) Reverse Severe Neurotoxicity and Rescue Juvenile Pigs from Lethal Doses of Micrurus Fulvius (Eastern Coral Snake) Venom. Toxins 2018, 10, 479. [Google Scholar] [CrossRef] [Green Version]
  80. Lewin, M.R.; María Gutiérrez, J.; Samuel, S.P.; Herrera, M.; Bryan-Quirós, W.; Lomonte, B.; Bickler, P.E.; Bulfone, T.C.; Williams, D.J. Delayed Oral LY333013 Rescues Mice from Highly Neurotoxic, Lethal Doses of Papuan Taipan (Oxyuranus Scutellatus) Venom. Toxins 2018, 10, 380. [Google Scholar] [CrossRef] [Green Version]
  81. Ferraz, C.R.; Arrahman, A.; Xie, C.; Casewell, N.R.; Lewis, R.J.; Kool, J.; Cardoso, F.C. Multifunctional Toxins in Snake Venoms and Therapeutic Implications: From Pain to Hemorrhage and Necrosis. Front. Ecol. Evol. 2019, 7, 218. [Google Scholar] [CrossRef] [Green Version]
  82. Still, K.B.M.; Nandlal, R.S.S.; Slagboom, J.; Somsen, G.W.; Casewell, N.R.; Kool, J. Multipurpose HTS Coagulation Analysis: Assay Development and Assessment of Coagulopathic Snake Venoms. Toxins 2017, 9, 382. [Google Scholar] [CrossRef] [PubMed]
Toxins 14 00736 g001 550
Figure 1. A schematic overview of the complete analytical and biochemical workflow. There are two main experiments shown in Figure 1 that run simultaneously in this study: (i) crude snake venom plate reader-based bioassaying (for assessing inhibition potential of anticoagulant activity by the small molecule inhibitors varespladib and marimastat), (ii) venom separation (using reversed-phase high-performance liquid chromatography; RP-HPLC or LC in short) coupled to UV detection followed by mass spectrometry (MS) with parallel nanofractionation for high-resolution fraction collection of separated toxins onto 384-well plates. Well plates with nanofractionated venom toxins are then vacuum centrifuged to dryness overnight, followed by bioassaying to assess anticoagulant activity. This data obtained is then processed to deliver bioassay chromatograms. Dried well plates can also be used for subsequent toxin identification by proteomics. Key: TIMS-TOF-MS, trapped ion mobility spectrometry time of flight mass spectrometry; HRMS, high-resolution mass spectrometry (HRMS); RP-nanoLC, reversed-phase nanoflow liquid chromatography.
Figure 1. A schematic overview of the complete analytical and biochemical workflow. There are two main experiments shown in Figure 1 that run simultaneously in this study: (i) crude snake venom plate reader-based bioassaying (for assessing inhibition potential of anticoagulant activity by the small molecule inhibitors varespladib and marimastat), (ii) venom separation (using reversed-phase high-performance liquid chromatography; RP-HPLC or LC in short) coupled to UV detection followed by mass spectrometry (MS) with parallel nanofractionation for high-resolution fraction collection of separated toxins onto 384-well plates. Well plates with nanofractionated venom toxins are then vacuum centrifuged to dryness overnight, followed by bioassaying to assess anticoagulant activity. This data obtained is then processed to deliver bioassay chromatograms. Dried well plates can also be used for subsequent toxin identification by proteomics. Key: TIMS-TOF-MS, trapped ion mobility spectrometry time of flight mass spectrometry; HRMS, high-resolution mass spectrometry (HRMS); RP-nanoLC, reversed-phase nanoflow liquid chromatography.
Toxins 14 00736 g001
Toxins 14 00736 g002 550
Figure 2. Venoms from elapids demonstrate anticoagulant activity on bovine plasma. Crude venom activity was measured as the mean area under the curve (AUC) of each venom concentration in the relevant clotting period (24–71 min), standardised by the mean AUC of the negative control in the same period. Lines were generated using a linear model and error bars represent the standard error of the mean (SEM) of 3–7 replicates (wells containing excess bubbles caused reading errors and had to be removed from analysis).
Figure 2. Venoms from elapids demonstrate anticoagulant activity on bovine plasma. Crude venom activity was measured as the mean area under the curve (AUC) of each venom concentration in the relevant clotting period (24–71 min), standardised by the mean AUC of the negative control in the same period. Lines were generated using a linear model and error bars represent the standard error of the mean (SEM) of 3–7 replicates (wells containing excess bubbles caused reading errors and had to be removed from analysis).
Toxins 14 00736 g002
Toxins 14 00736 g003 550
Figure 3. Anticoagulant venom activity of cobra venoms are inhibited by the PLA2 inhibitor varespladib. H. haemachatus, N. haje, and N. naja additionally show significant inhibition by some marimastat doses. Clotting activity was measured as the mean area under the curve (AUC) of each sample type in the relevant clotting period (22–71 min), standardised by the mean AUC of the negative control in the same period. Error bars represent the standard error of the mean (SEM) of 3–4 measurements (some sample wells needed to be removed due to disruption in readings caused by bubbles). Asterisks represent significant increases in coagulant activity as measured by one-way ANOVA and Dunnett’s Multiple Comparison post-hoc testing; *—p < 0.05, **—p < 0.01, ***—p < 0.0001.
Figure 3. Anticoagulant venom activity of cobra venoms are inhibited by the PLA2 inhibitor varespladib. H. haemachatus, N. haje, and N. naja additionally show significant inhibition by some marimastat doses. Clotting activity was measured as the mean area under the curve (AUC) of each sample type in the relevant clotting period (22–71 min), standardised by the mean AUC of the negative control in the same period. Error bars represent the standard error of the mean (SEM) of 3–4 measurements (some sample wells needed to be removed due to disruption in readings caused by bubbles). Asterisks represent significant increases in coagulant activity as measured by one-way ANOVA and Dunnett’s Multiple Comparison post-hoc testing; *—p < 0.05, **—p < 0.01, ***—p < 0.0001.
Toxins 14 00736 g003
Toxins 14 00736 g004 550
Figure 4. Correlation of LC-UV and MS chromatogram and bioassay chromatogram of D. polylepis venom. (A) UV chromatogram of D. polylepis venom (1 mg/mL, 50 μL injection volume, post-column split into 1:9 ratio. The smaller portion went to UV (220 and 254 nm recorded), and then to MS detection. Superimposed bioassay chromatograms resulting from analyses of D. polylepis venom (1 mg/mL, 50 µL injection volume) in the presence of different concentrations of (B) varespladib and (C) marimastat. (D) Superimposed anticoagulation bioassay chromatograms resulting from analyses of several diluted D. polylepis venom samples ranging from 1 mg/mL to 0.2 mg/mL (50 µL per injection). (E) Base peak chromatogram (BPC) and extracted ion currents (XIC) mass spectrometry. Bioassay chromatograms are correlated to UV and MS. Proteomic annotation of D. polylepis venom using Mascot software against (F) the Swiss-Prot database and (G) species-specific venom gland transcriptomic-derived database. The protein score represents the probability that designated proteins, defined by elution time, are present in the sample.
Figure 4. Correlation of LC-UV and MS chromatogram and bioassay chromatogram of D. polylepis venom. (A) UV chromatogram of D. polylepis venom (1 mg/mL, 50 μL injection volume, post-column split into 1:9 ratio. The smaller portion went to UV (220 and 254 nm recorded), and then to MS detection. Superimposed bioassay chromatograms resulting from analyses of D. polylepis venom (1 mg/mL, 50 µL injection volume) in the presence of different concentrations of (B) varespladib and (C) marimastat. (D) Superimposed anticoagulation bioassay chromatograms resulting from analyses of several diluted D. polylepis venom samples ranging from 1 mg/mL to 0.2 mg/mL (50 µL per injection). (E) Base peak chromatogram (BPC) and extracted ion currents (XIC) mass spectrometry. Bioassay chromatograms are correlated to UV and MS. Proteomic annotation of D. polylepis venom using Mascot software against (F) the Swiss-Prot database and (G) species-specific venom gland transcriptomic-derived database. The protein score represents the probability that designated proteins, defined by elution time, are present in the sample.
Toxins 14 00736 g004
Toxins 14 00736 g005 550
Figure 5. Correlation of LC-UV and MS chromatogram and bioassay chromatogram of D. angusticeps venom. (A) UV chromatogram of D. angusticeps venom (1 mg/mL, 50 μL injection volume, post-column split into 1:9 ratio. The smaller portion went to UV (220 and 254 nm recorded); and then to MS detection. Superimposed bioassay chromatograms resulting from analyses of D. angusticeps venom (1 mg/mL, 50 µL injection volume) in the presence of different concentrations of (B) varespladib and (C) marimastat. (D) Superimposed anticoagulation bioassay chromatograms resulting from analyses of several diluted D. angusticeps venom samples ranging from 1 mg/mL to 0.2 mg/mL (50 µL per injection). (E) Base peak chromatogram (BPC) and extracted ion currents (XIC) mass spectrometry. Bioassay chromatograms are correlated to UV and MS. Proteomic annotation of D. angusticeps venom using Mascot software against (F) the Swiss-Prot database and (G) species-specific venom gland transcriptomic-derived database. The protein score represents the probability that designated proteins, defined by elution time, are present in the sample.
Figure 5. Correlation of LC-UV and MS chromatogram and bioassay chromatogram of D. angusticeps venom. (A) UV chromatogram of D. angusticeps venom (1 mg/mL, 50 μL injection volume, post-column split into 1:9 ratio. The smaller portion went to UV (220 and 254 nm recorded); and then to MS detection. Superimposed bioassay chromatograms resulting from analyses of D. angusticeps venom (1 mg/mL, 50 µL injection volume) in the presence of different concentrations of (B) varespladib and (C) marimastat. (D) Superimposed anticoagulation bioassay chromatograms resulting from analyses of several diluted D. angusticeps venom samples ranging from 1 mg/mL to 0.2 mg/mL (50 µL per injection). (E) Base peak chromatogram (BPC) and extracted ion currents (XIC) mass spectrometry. Bioassay chromatograms are correlated to UV and MS. Proteomic annotation of D. angusticeps venom using Mascot software against (F) the Swiss-Prot database and (G) species-specific venom gland transcriptomic-derived database. The protein score represents the probability that designated proteins, defined by elution time, are present in the sample.
Toxins 14 00736 g005
Toxins 14 00736 g006 550
Figure 6. Correlation of LC-UV and MS chromatogram and bioassay chromatogram of N. naja venom. (A) UV chromatogram of N. naja venom (1 mg/mL, 50 μL injection volume, post-column split into 1:9 ratio of which the smaller portion went to UV (220 and 254 nm recorded); and then to MS detection. Superimposed bioassay chromatograms resulting from analyses of N. naja venom (0.2 mg/mL, 50 µL injection volume) in the presence of different concentrations of (B) varespladib. (C) Superimposed anticoagulation bioassay chromatograms resulting from analyses of several diluted N. naja venom samples ranging from 1 mg/mL to 0.008 mg/mL (50 µL per injection). (D) Base peak chromatogram (BPC) and extracted ion currents (XIC) mass spectrometry. Bioassay chromatograms are correlated to UV and MS. Proteomic annotation of N. naja venom using Mascot software against (E) the Swiss-Prot database and (F) species-specific venom gland transcriptomic-derived database.
Figure 6. Correlation of LC-UV and MS chromatogram and bioassay chromatogram of N. naja venom. (A) UV chromatogram of N. naja venom (1 mg/mL, 50 μL injection volume, post-column split into 1:9 ratio of which the smaller portion went to UV (220 and 254 nm recorded); and then to MS detection. Superimposed bioassay chromatograms resulting from analyses of N. naja venom (0.2 mg/mL, 50 µL injection volume) in the presence of different concentrations of (B) varespladib. (C) Superimposed anticoagulation bioassay chromatograms resulting from analyses of several diluted N. naja venom samples ranging from 1 mg/mL to 0.008 mg/mL (50 µL per injection). (D) Base peak chromatogram (BPC) and extracted ion currents (XIC) mass spectrometry. Bioassay chromatograms are correlated to UV and MS. Proteomic annotation of N. naja venom using Mascot software against (E) the Swiss-Prot database and (F) species-specific venom gland transcriptomic-derived database.
Toxins 14 00736 g006
Toxins 14 00736 g007 550
Figure 7. Correlation of LC-UV and MS chromatogram and bioassay chromatogram of N. pallida venom. (A) UV chromatogram of N. pallida venom (1 mg/mL, 50 μL injection volume, post-column split into 1:9 ratio of which the smaller portion went to UV (220 and 254 nm recorded); and then to MS detection. Superimposed bioassay chromatograms resulting from analyses of N. pallida venom (0.2 mg/mL, 50 µL injection volume) in the presence of different concentrations of (B) varespladib. (C) Superimposed anticoagulation bioassay chromatograms resulting from analyses of several diluted N. pallida venom samples ranging from 1 mg/mL to 0.008 mg/mL (50 µL per injection). (D) Base peak chromatogram (BPC) and extracted ion currents (XIC) mass spectrometry. Bioassay chromatograms are correlated to UV and MS. Proteomic annotation of N. pallida venom using Mascot software against (E) the Swiss-Prot database and (F) species-specific venom gland transcriptomic-derived database. The protein score represents the probability that designated proteins, defined by elution time, are present in the sample.
Figure 7. Correlation of LC-UV and MS chromatogram and bioassay chromatogram of N. pallida venom. (A) UV chromatogram of N. pallida venom (1 mg/mL, 50 μL injection volume, post-column split into 1:9 ratio of which the smaller portion went to UV (220 and 254 nm recorded); and then to MS detection. Superimposed bioassay chromatograms resulting from analyses of N. pallida venom (0.2 mg/mL, 50 µL injection volume) in the presence of different concentrations of (B) varespladib. (C) Superimposed anticoagulation bioassay chromatograms resulting from analyses of several diluted N. pallida venom samples ranging from 1 mg/mL to 0.008 mg/mL (50 µL per injection). (D) Base peak chromatogram (BPC) and extracted ion currents (XIC) mass spectrometry. Bioassay chromatograms are correlated to UV and MS. Proteomic annotation of N. pallida venom using Mascot software against (E) the Swiss-Prot database and (F) species-specific venom gland transcriptomic-derived database. The protein score represents the probability that designated proteins, defined by elution time, are present in the sample.
Toxins 14 00736 g007
Toxins 14 00736 g008 550
Figure 8. Correlation of LC-UV and MS chromatogram and bioassay chromatogram of N. nigricollis venom. (A) UV chromatogram of N. nigricollis venom (1 mg/mL, 50 μL injection volume, post-column split into 1:9 ratio and a smaller portion went to UV (220 and 254 nm recorded); and then to MS detection. Superimposed bioassay chromatograms resulting from analyses of N. nigricollis venom (0.2 mg/mL, 50 µL injection volume) in the presence of different concentrations of (B) varespladib. (C) Superimposed anticoagulation bioassay chromatograms resulting from analyses of several diluted N. nigricollis venom samples ranging from 1 mg/mL to 0.008 mg/mL (50 µL per injection). (D) Base peak chromatogram (BPC) and extracted ion currents (XIC) mass spectrometry. Bioassay chromatograms are correlated to UV and MS. Proteomic annotation of N. nigricollis venom using Mascot software against (E) the Swiss-Prot database and (F) species-specific venom gland transcriptomic-derived database. The protein score represents the probability that designated proteins, defined by elution time, are present in the sample.
Figure 8. Correlation of LC-UV and MS chromatogram and bioassay chromatogram of N. nigricollis venom. (A) UV chromatogram of N. nigricollis venom (1 mg/mL, 50 μL injection volume, post-column split into 1:9 ratio and a smaller portion went to UV (220 and 254 nm recorded); and then to MS detection. Superimposed bioassay chromatograms resulting from analyses of N. nigricollis venom (0.2 mg/mL, 50 µL injection volume) in the presence of different concentrations of (B) varespladib. (C) Superimposed anticoagulation bioassay chromatograms resulting from analyses of several diluted N. nigricollis venom samples ranging from 1 mg/mL to 0.008 mg/mL (50 µL per injection). (D) Base peak chromatogram (BPC) and extracted ion currents (XIC) mass spectrometry. Bioassay chromatograms are correlated to UV and MS. Proteomic annotation of N. nigricollis venom using Mascot software against (E) the Swiss-Prot database and (F) species-specific venom gland transcriptomic-derived database. The protein score represents the probability that designated proteins, defined by elution time, are present in the sample.
Toxins 14 00736 g008
Toxins 14 00736 g009 550
Figure 9. Correlation of LC-UV and MS chromatogram and bioassay chromatogram of N. haje venom. (A) UV chromatogram of N. haje venom (1 mg/mL, 50 μL injection volume, post-column split into 1:9 ratio of which the smaller portion went to UV (220 and 254 nm recorded); and then to MS detection. Superimposed bioassay chromatograms resulting from analyses of N. haje venom (0.2 mg/mL, 50 µL injection volume) in the presence of different concentrations of (B) varespladib. (C) Superimposed anticoagulation bioassay chromatograms resulting from analyses of several diluted N. haje venom samples ranging from 1 mg/mL to 0.008 mg/mL (50 µL per injection). (D) Base peak chromatogram (BPC) and extracted ion currents (XIC) mass spectrometry. Bioassay chromatograms are correlated to UV and MS. Proteomic annotation of N. haje venom using Mascot software against (E) the Swiss-Prot database and (F) species-specific venom gland transcriptomic-derived database. The protein score represents the probability that designated proteins, defined by elution time, are present in the sample.
Figure 9. Correlation of LC-UV and MS chromatogram and bioassay chromatogram of N. haje venom. (A) UV chromatogram of N. haje venom (1 mg/mL, 50 μL injection volume, post-column split into 1:9 ratio of which the smaller portion went to UV (220 and 254 nm recorded); and then to MS detection. Superimposed bioassay chromatograms resulting from analyses of N. haje venom (0.2 mg/mL, 50 µL injection volume) in the presence of different concentrations of (B) varespladib. (C) Superimposed anticoagulation bioassay chromatograms resulting from analyses of several diluted N. haje venom samples ranging from 1 mg/mL to 0.008 mg/mL (50 µL per injection). (D) Base peak chromatogram (BPC) and extracted ion currents (XIC) mass spectrometry. Bioassay chromatograms are correlated to UV and MS. Proteomic annotation of N. haje venom using Mascot software against (E) the Swiss-Prot database and (F) species-specific venom gland transcriptomic-derived database. The protein score represents the probability that designated proteins, defined by elution time, are present in the sample.
Toxins 14 00736 g009
Toxins 14 00736 g010 550
Figure 10. Correlation of LC-UV and MS chromatogram and bioassay chromatogram of H. Haemachatus venom. (A) UV chromatogram of H. haemachatus venom (1 mg/mL, 50 μL injection volume, post-column split into 1:9 ratio. The smaller portion went to UV (220 and 254 nm recorded), and then to MS detection. Superimposed bioassay chromatograms resulting from analyses of H. haemachatus venom (0.2 mg/mL, 50 µL injection volume) in the presence of different concentrations of (B) varespladib. (C) Superimposed anticoagulation bioassay chromatograms resulting from analyses of several diluted H. haemachatus venom samples ranging from 1 mg/mL to 0.008 mg/mL (50 µL per injection). (D) Base peak chromatogram (BPC) and extracted ion currents (XIC) mass spectrometry. Bioassay chromatograms are correlated to UV and MS. Proteomic annotation of H. haemachatus venom using Mascot software against (E) the Swiss-Prot database and (F) species-specific venom gland transcriptomic-derived database. The protein score represents the probability that designated proteins, defined by elution time, are present in the sample. Note: after nanofractionation analytics, no procoagulant toxin peaks were observed (not shown in the figure). SVSPs and other mainly larger toxins can denature during the chromatographic separation thereby losing their biological activity (which we have observed in other ongoing research projects).
Figure 10. Correlation of LC-UV and MS chromatogram and bioassay chromatogram of H. Haemachatus venom. (A) UV chromatogram of H. haemachatus venom (1 mg/mL, 50 μL injection volume, post-column split into 1:9 ratio. The smaller portion went to UV (220 and 254 nm recorded), and then to MS detection. Superimposed bioassay chromatograms resulting from analyses of H. haemachatus venom (0.2 mg/mL, 50 µL injection volume) in the presence of different concentrations of (B) varespladib. (C) Superimposed anticoagulation bioassay chromatograms resulting from analyses of several diluted H. haemachatus venom samples ranging from 1 mg/mL to 0.008 mg/mL (50 µL per injection). (D) Base peak chromatogram (BPC) and extracted ion currents (XIC) mass spectrometry. Bioassay chromatograms are correlated to UV and MS. Proteomic annotation of H. haemachatus venom using Mascot software against (E) the Swiss-Prot database and (F) species-specific venom gland transcriptomic-derived database. The protein score represents the probability that designated proteins, defined by elution time, are present in the sample. Note: after nanofractionation analytics, no procoagulant toxin peaks were observed (not shown in the figure). SVSPs and other mainly larger toxins can denature during the chromatographic separation thereby losing their biological activity (which we have observed in other ongoing research projects).
Toxins 14 00736 g010
Table
Table 1. Overview of tentatively assigned venom toxins associated with the observed anticoagulation peaks. The table shows m/z values, including the charge state measured for each m/z value of intact toxins observed in MS, their retention times, and the calculated accurate mass for each toxin. Moreover, proteomics-derived Mascot (Swiss-Prot) protein identification and Mascot-derived corresponding masses are shown, as a match to the species-specific venom gland transcriptome databases. Hits were obtained by database searches of nanoLC–MS/MS data measured from tryptic digests of the respective venom toxins after nanofractionation and collecting these venom toxins from their respective wells in which they were fractionated. As many toxins, especially the larger toxins such as SVMPs, can have post-translational modifications (PTMs) such as glycosylation, the molecular mass (mol. mass (in Da)) found by high resolution mass spectrometry does not always match the Mascot exact mass (in Da).
Table 1. Overview of tentatively assigned venom toxins associated with the observed anticoagulation peaks. The table shows m/z values, including the charge state measured for each m/z value of intact toxins observed in MS, their retention times, and the calculated accurate mass for each toxin. Moreover, proteomics-derived Mascot (Swiss-Prot) protein identification and Mascot-derived corresponding masses are shown, as a match to the species-specific venom gland transcriptome databases. Hits were obtained by database searches of nanoLC–MS/MS data measured from tryptic digests of the respective venom toxins after nanofractionation and collecting these venom toxins from their respective wells in which they were fractionated. As many toxins, especially the larger toxins such as SVMPs, can have post-translational modifications (PTMs) such as glycosylation, the molecular mass (mol. mass (in Da)) found by high resolution mass spectrometry does not always match the Mascot exact mass (in Da).
SpeciesR. Time (min)m/zChargeMol. Mass (Da)Protein Identification/Mascot HitsMascot Exact Mass (Da)Species-specific Database HitsToxin Class
Dendroaspis polylepis14.41019.937+7129.47Three-finger toxinN.A.N.A.3FTx
15.31007.317+7040.16Three-finger toxinN.A.N.A3FTx
15.81006.037+7031.173SLS_DENPO (Calciseptin) Dendroaspis polylepis7031.21D_polylepis_T00103FTx
16.4999.348+7982.683L24_DENPO (Alpha-elapitoxin-Dpp2d) Dendroaspis polylepis7986.68N.A.3FTx
16.4999.34+2625,957.30Snake venom metalloprotease23433.41 (non-glycosylated)D_polylepis_T3990_T0439SVMP
16.81055.347+7377.363SIM3_DENAN (Muscarinic toxin 3) Dendroaspis angusticeps7392.54D_polylepis_T23203FTx
16.81055.0536+37,947.16Snake venom metalloproteinase27,106.10 (non-glycosylated)D_polylepis_T0167SVMP
17.71081.836+6481.97Three finger toxinsN.A.D_polylepis_T23313FTx
Dendroaspis angusticeps15.1965.147+6745.913SE2_DENAN (Fasciculin-2) Dendroaspis angusticeps6744.96D_angusticeps_T1380_T26423FTx
15.5971.997+6793.933SE1_DENAN (Fasciculin-1) Dendroaspis angusticeps6793.98D_angusticeps_T27373FTx
16.1943.447+6593.05Three-finger toxinN.A.D_angusticeps_T10113FTx
16.4945.447+6607.07VKTHE_DENAN (Kunitz-type serine protease inhibitor long epsilon-dendrotoxin His55) Dendroaspis angusticeps6609.17D_angusticeps_T3547Kunitz-type serine protease
17.11073.637+7504.033SI1B_DENAN (Rho-elapitoxin-Da1b) Dendroaspis angusticeps7507.48D_angusticeps_T44053FTx
1048.767+7330.293NOJ_DENAN (Toxin Tx7335) Dendroaspis angusticeps7330.33D_angusticeps_T37773FTx
17.11045.0559+61,599.77Snake venom metalloproteinase62,965.82D_angusticeps_T0082SVMP
18.2951.357+6648.353SA5_NAJKA (Cytotoxin 5) Naja kaouthia6846.95N.A.3FTx
Naja naja15.31116.367+7803.473L22_NAJNA (Long neurotoxin 2) Naja naja7805.48N.A.3FTx
15.71118.797+7819.463L23_NAJNA (Long neurotoxin 3) Naja naja7817.54N.A.3FTx
16.31663.708+13,293.52PA2B3_NAJMO (Basic phospholipase A2 CM-III) Naja mossambica13,288.82N.A.PLA2
1678.078+13,408.51PA2A2_NAJME (Acidic phospholipase A2 DE-II) Naja melanoleuca13,413.07N.A.PLA2
16.4972.197+6795.273SA4_NAJHA (Cytotoxin 4) Naja haje6794.58N.A.3FTx
16.51482.079+13,321.59PA2A2_NAJNA (Acidic phospholipase A2 2) Naja naja13,322.64N.A.PLA2
1474.189+13,249.27PA2B4_NAJNG (Basic phospholipase A2) Naja nigricollis13,245.22N_naja_T0142_T1042_T2728_T2640_
T2414_T1496_T0668_T2421		PLA2
1447.129+13,007.02NGFV_NAJNA (Venom nerve growth factor) Naja naja13,008.10N.A.vNGF
1681.588+13,435.58Phospholipase A2N.A.N.A.PLA2
16.61649.458+13,179.51PA2A_NAJAT (Acidic phospholipase A2 natratoxin) Naja atra13,174.25N.A.PLA2
1672.588+13,364.56Phospholipase A2 N.A.N.A.PLA2
16.62217.7612+26,598.04Snake venom metalloproteinaseN.A.N.A.SVMP
17.21001.787+7002.443SOFL_NAJNA (Cytotoxin like basic protein) Naja naja7006.56N_naja_T09173FTx
17.4965.917+6750.333SAA_NAJNA (Cytotoxin 10) Naja naja
3SA2_NAJNA (Cytotoxin 2) Naja naja		6751.34
6750.36		N_naja_T2400_T1263_T0382_T23243FTx
3FTx
963.337+6733.323SA3_NAJNA (Cytotoxin 3) Naja naja6732.34N.A.3FTx
17.6969.7+7+6778.333SA1_NAJNA (Cytotoxin 1) Naja naja
3SA7_NAJNA (Cytotoxin 7) Naja naja
3SA8_NAJNA (Cytotoxin 8) Naja naja		6778.32
6779.30
6780.28		N_naja_T04813FTx
3FTx
3FTx
18.31113.056+6669.273SA9_NAJNA (Cytotoxin 9) Naja naja6661.77N_naja_T2420_T2418_T2687_T16723FTx
Naja pallida12.5970.747+6785.083SA4_NAJHA (Cytotoxin 4) Naja haje6789.37N_pallida_T1076
N_pallida_T0738		3FTx
14.9987.357+6901.42Three finger toxinsN.A.N.A.3FTx
16.21478.769+13,290.79(Basic phospholipase A2) Naja pallida13,266.83N_pallida_T0443PLA2
16.41496.209+13,205.74Phospholipase A2N.A.N.A.PLA2
1666.728+13,320.73PA2B_NAJPA (Basic phospholipase A2 nigexine) Naja pallida13,316.93N.A.PLA2
963.197+6732.29Three finger toxins6731.85N_pallida_T2193_T16903FTx
17.3994.067+6948.363S1CB_NAJNA (Cobrotoxin homolog) Naja naja6943.98N.A.3FTx
17.5970.347+6782.343S11_NAJPA (Short neurotoxin 1) Naja pallida6782.10N_pallida_T09543FTx
17.8984.777+6883.393SA3_NAJMO (Cytotoxin 3) Naja mossambica 6881.42N_pallida_T1038_T13233FTx
18.2975.057+6815.313SA1_NAJPA (Cytotoxin 1) Naja pallida6814.31N_pallida_T1011_T2149_T1010_T08893FTx
Naja nigricollis15.61036.757+7247.19Three-finger toxinN.A.N.A.3FTx
16.11473.659+13,245.77PA2B4_NAJNG (Phospholipase A2 Basic) Naja nigricollis13,244.90N_nigricollis_T2086_T1848_T3198PLA2
16.81477.429+13,279.66Acidic phospholipase A213,279.08N_nigricollis_T3199_T0053PLA2
17.0959.337+6705.273SA4_NAJMO (Cytotoxin 4) Naja mossambica6702.34N.A.3FTx
975.337+6816.293SA1_NAJPA (Cytotoxin 1) Naja pallida6814.31N_nigricollis_VG_T3171_T16483FTx
17.2965.907+6751.263SA9_NAJNA (Cytotoxin 9) Naja naja6570.36N.A.3FTx
1656.979+13,167.58PA2A1_NAJMO (Acidic phospholipase A2 CM-I13,195.75N.A.PLA2
17.4956.177+6683.203SAN_NAJNG (Naniproin) Naja nigricollis6682.41N_nigricollis_VG_T27713FTx
17.7947.487+6811.323SA5_NAJAT (Cytotoxin 5) Naja atra6810.35N.A.3FTx
18.0984.777+6883.383SA3_NAJMO (Cytotoxin 3) Naja mossambica6881.42N.A.3FTx
18.4975.057+6815.313SA1_NAJMO (Cytotoxin 1) Naja mossambica6813.33N_nigricollis_VG_T1821_T11783FTx
Naja haje14.4963.218+7692.63Three finger toxinsN.A.N.A.3FTx
16.11031.167+7208.093SI3_NAJMO (Short neurotoxin 3) Naja mossambica7210.58N.A.3FTx
16.51552.389+13,963.36PA2A3_NAJSG (Acidic phospholipase A2 3) Naja sagitifera 13,955.23N.A.PLA2
16.9981.897+6863.213SA2_NAJME (Cytotoxin 2) Naja melanoleuca6863.54N_haje_T2831_T1704_T3906_T3905_
T0971		3FTx
17.2967.317+6761.183SA5_NAJHH (Cytotoxin 5) Naja haje6765.18N.A.3FTx
17.6974.737+6813.293SA1_NAJMO (Cytotoxin 1) Naja mossambica6813.33N.A.3FTx
17.9979.327+6845.243SA2_NAJHA (Cytotoxin 2) Naja haje
3SA5_NAJHA (Cytotoxin 5) Naja haje
3SA6_NAJHA (Cytotoxin 6) Naja haje
3SA7_NAJHA (Cytotoxin 7) Naja haje		6845.24
6843.27
6844.25
6843.27		N.A.3FTx
3FTx
3FTx
3FTx
Hemachatus haemachatus15.31014.457+7091.103SUB_DENAN (Muscarinic toxin) Dendroaspis angusticeps7092.35N.A.3FTx
15.81501.539+13,498.76PA2B1_HEMHA (Basic phospholipase A2 DE-1) Hemachatus haemachatus13,495.94H_haemachatus_T1146PLA2
15.81031.227+7208.50Three finger toxinsN.A.N.A.3FTx
16.4971.217+6788.433S11_NAJPA (Short neurotoxin 1) Naja pallida
3SA4_NAJHA (Cytotoxin 4) Naja haje
3SB2_HEMHA (Cytotoxin 2) Hemachatus haemachatus		6782.10
6789.37
6787.44		H_haemachatus_T12353FTx
3FTx
3FTx
16.7977.507+6832.453S11_NAJHA (Short neurotoxin 1) Naja haje
3SBH_HEMHA (Three-finger hemachatoxin) Hemachatus haemachatus
3SB1_HEMHA (Cytotoxin 1) Hemachatus haemachatus		6831.05
6381.45
6831.45		H_haemachatus_T1274_T1588_T1866_
T1092_T1440_T1175		3FTx
3FTx
3FTx
18.0972.347+6795.353SA8_NAJHA (Cytotoxin 8) Naja haje6799.32N.A.3FTx
18.2969.917+6779.363SB3_HEMHA (Cytotoxin 3) Hemachatus haemachatus6780.32N.A.3FTx
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Arrahman, A.; Kazandjian, T.D.; Still, K.B.M.; Slagboom, J.; Somsen, G.W.; Vonk, F.J.; Casewell, N.R.; Kool, J. A Combined Bioassay and Nanofractionation Approach to Investigate the Anticoagulant Toxins of Mamba and Cobra Venoms and Their Inhibition by Varespladib. Toxins 2022, 14, 736. https://doi.org/10.3390/toxins14110736

AMA Style

Arrahman A, Kazandjian TD, Still KBM, Slagboom J, Somsen GW, Vonk FJ, Casewell NR, Kool J. A Combined Bioassay and Nanofractionation Approach to Investigate the Anticoagulant Toxins of Mamba and Cobra Venoms and Their Inhibition by Varespladib. Toxins. 2022; 14(11):736. https://doi.org/10.3390/toxins14110736

Chicago/Turabian Style

Arrahman, Arif, Taline D. Kazandjian, Kristina B. M. Still, Julien Slagboom, Govert W. Somsen, Freek J. Vonk, Nicholas R. Casewell, and Jeroen Kool. 2022. "A Combined Bioassay and Nanofractionation Approach to Investigate the Anticoagulant Toxins of Mamba and Cobra Venoms and Their Inhibition by Varespladib" Toxins 14, no. 11: 736. https://doi.org/10.3390/toxins14110736

APA Style

Arrahman, A., Kazandjian, T. D., Still, K. B. M., Slagboom, J., Somsen, G. W., Vonk, F. J., Casewell, N. R., & Kool, J. (2022). A Combined Bioassay and Nanofractionation Approach to Investigate the Anticoagulant Toxins of Mamba and Cobra Venoms and Their Inhibition by Varespladib. Toxins, 14(11), 736. https://doi.org/10.3390/toxins14110736

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop