Next Article in Journal
Impact of Microcystin-LR on Liver Function Varies by Dose and Sex in Mice
Previous Article in Journal
Current Status and Future Opportunities of Omics Tools in Mycotoxin Research
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Toxins
Volume 10
Issue 11
10.3390/toxins10110434
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

Proteomic Characterization of Two Medically Important Malaysian Snake Venoms, Calloselasma rhodostoma (Malayan Pit Viper) and Ophiophagus hannah (King Cobra)

by
Sugita Kunalan
1,
Iekhsan Othman
1,
Sharifah Syed Hassan
1 and
Wayne C. Hodgson
2,*
1
Jeffrey Cheah School of Medicine and Health Sciences, Monash University Malaysia, Jalan Lagoon Selatan, Bandar Sunway, Selangor Darul Ehsan 47500, Malaysia
2
Monash Venom Group, Department of Pharmacology, Faculty of Medicine, Nursing and Health Sciences, Monash University, Victoria 3800, Australia
*
Author to whom correspondence should be addressed.
Toxins 2018, 10(11), 434; https://doi.org/10.3390/toxins10110434
Submission received: 18 October 2018 / Revised: 23 October 2018 / Accepted: 23 October 2018 / Published: 26 October 2018
(This article belongs to the Section Animal Venoms)
Browse Figures
Versions Notes

Abstract

:
Calloselasma rhodostoma (CR) and Ophiophagus hannah (OH) are two medically important snakes found in Malaysia. While some studies have described the biological properties of these venoms, feeding and environmental conditions also influence the concentration and distribution of snake venom toxins, resulting in variations in venom composition. Therefore, a combined proteomic approach using shotgun and gel filtration chromatography, analyzed by tandem mass spectrometry, was used to examine the composition of venoms from these Malaysian snakes. The analysis revealed 114 proteins (15 toxin families) and 176 proteins (20 toxin families) in Malaysian Calloselasma rhodostoma and Ophiophagus hannah species, respectively. Flavin monoamine oxidase, phospholipase A2, phosphodiesterase, snake venom metalloproteinase, and serine protease toxin families were identified in both venoms. Aminopeptidase, glutaminyl-peptide cyclotransferase along with ankyrin repeats were identified for the first time in CR venom, and insulin, c-type lectins/snaclecs, hepatocyte growth factor, and macrophage colony-stimulating factor together with tumor necrosis factor were identified in OH venom for the first time. Our combined proteomic approach has identified a comprehensive arsenal of toxins in CR and OH venoms. These data may be utilized for improved antivenom production, understanding pathological effects of envenoming, and the discovery of biologically active peptides with medical and/or biotechnological value.
Keywords:
Calloselasma rhodostoma; Ophiophagus hannah; toxin; proteomics; venomics; venom proteome
Key Contribution: Though earlier investigations have provided insight into the major families of toxins in these Malaysian venoms, our work has identified several yet undiscovered proteins: Aminopeptidase; glutaminyl-peptide cyclotransferase and ankyrin repeats in CR venom as well as insulin; and c-type lectins along with several notable growth factors in OH venom.

1. Introduction

Approximately 5.5 million people are victims of snake bites every year, causing about 400,000 amputations, and 20,000 to 125,000 fatalities [1], representing a global public health burden. Despite its destructive propensity, snake venom has been studied to identify lead compounds for use as therapeutic agents. Promising reports of antimicrobial (bacteria, virus, fungi, and protozoa), anti-inflammatory, antinociceptive, anti-cancer, anti-hypertensive, anti-coagulation, and anti-fibrinolytic activities of snake venoms [2,3,4,5,6,7] have initiated more intensive investigations of venom and venom glands using venomics [8,9]. However, it is widely hypothesized that feeding and environmental conditions influence the concentration and distribution of snake venom components, causing variation between individual, species, genus, and families [10]. Such variation can have major medical implications for the treatment of human snakebite victims. Thus, the reason for the conception of venomics and optimism is that characterization of the large protein variability within snake venoms may contribute to a deeper understanding of their biological as well as pathological effects.
In this study, we investigated the proteomic composition of Calloselasma rhodostoma and Ophiohagus hannah venoms, two medically important snakes found in Malaysia [11]. Both venoms are a rich source of proteins that have been reported to contain biologically active peptides with medical and biotechnological value. One such example is Ancrod (or Viprinex®), a snake venom serine protease purified from C. rhodostoma venom. Ancrod is a thrombin-like enzyme that breaks down fibrinogen and was developed for the treatment of ischaemic stroke and myocardial infarction as well as deep-vein thrombosis [12,13]. However, Ancrod has been suspended from use due to inconclusive clinical trial results [14,15]. Another example is the α-neurotoxin, hannalgesin, that was derived from a small peptide called prohanin, originating from the venom of O. hannah [16]. Despite being an α-neurotoxin, hannalgesin does not paralyze skeletal muscles, but provides analgesic effects. However, the development of the drug was terminated in the preclinical phase [17,18].
Until recently, there have been relatively few reports on the proteomic profile of these two medically important Malaysian snakes. In a comprehensive study by Tang et al., a total of 96 proteins were identified for C. rhodostoma venom [19]. The study employed a series of chromatographic separations, followed by in-gel tryptic digestion of one-dimension (1D) gel electrophoresis proteins and, lastly, identification based on mass spectrometry. In earlier investigations, C. rhodostoma venom proteins were studied using two-dimension (2D) gel electrophoresis together with mass spectrometry. Ali et al. (2013) reported the presence of disintegrin, kallikrein (thrombin-like protein), L-amino acid oxidase (LAAO), lectin, and phospholipase A2 (acidic, basic, and neutral) as well as snake venom metalloproteinases (SVMPs) belonging to six protein families, whereas Vejayan et al., (2014) presented evidence for 26 proteins that could be grouped into six protein families, as demonstrated by Ali et al. [20,21]. Despite the positive identification of toxins, there were 89 unidentified spots that indicates a substantial gap in the complete venom profiling of the C. rhodostoma venom. Apart from the relatively small number of identified proteins, both studies did not provide information on the relative abundance of the identified proteins.
An elaborate venom-gland transcriptomic and venom proteomic study of O. hannah venom revealed 128 transcripts and 45 proteins, which were classified into 16 individual protein families [22]. The major protein families reported were three-finger toxin (3FTx), SVMP, cysteine-rich secretory protein (CRiSP), LAAO, Vespryn, phospholipase A2 (PLA2), and cobra venom factor (CVF). In addition, nine low abundance (≤2%) protein families were detected. The results of the study were compared with the Balinese (Indonesian) O. hannah venom investigation by Vonk et al. (2013) [23]. In contrast, they reported 73 proteins using a combination of shotgun and 2D gel electrophoresis. The major protein families described for Balinese O. hannah venom were highly similar to the Malaysian O. hannah venom, in which the most abundant toxins were 3FTxs and SVMPs. Although similarities were detected, we know that intraspecies toxin composition variations exist between venoms. As mentioned above, Vejayan et al. (2014) also investigated the proteomic profile of Malaysian O. hannah venom together with C. rhodostoma venom. Only 15 proteins belonging to five toxin families could be deduced from 50 spots that were isolated from the 2D gel. The large number of uncharacterized proteins, coupled with the discrepancy of the number of identified proteins between the Malaysian and Indonesian species, indicates the need for further proteomic exploration of Malaysian O. hannah venom.
In this work, we employed a combination of two proteomic approaches to provide a more comprehensive proteomic profile of two medically significant venomous snakes from Malaysia. The first approach involves a direct tryptic digest of freeze-dried crude venom proteins followed by ESI-LC-MS/MS, whereas the second approach subjects the crude venom proteins to a molecular-weight based separation prior to tryptic digestion and subsequent ESI-LC-MS/MS. Using this strategy, we report a total of 114 proteins for C. rhodostoma and 176 proteins for O. hannah venoms, successfully identifying more proteins and peptides than previously reported.

2. Results

We characterized the venom composition of two Malaysian snakes, i.e., Calloselasma rhodostoma and Ophiophagus hannah, using two proteomic approaches: Direct tryptic digest of freeze-dried crude venom proteins followed by liquid chromatography coupled with electrospray tandem mass spectrometry (shotgun-ESI-LC-MS/MS) and tryptic digest of gel-filtration chromatography venom protein fractions and subsequent liquid chromatography coupled with electrospray tandem mass spectrometry (GF-ESI-LC-MS/MS). The combined approaches revealed a diverse set of data that has been summarized in Table 1, Table 2, Table 3 and Table 4 and is discussed in detail below. A complementary visual representation of the crude venom composition is also presented.

2.1. SDS-PAGE of Crude Venom Proteins

The overall protein composition of C. rhodostoma and O. hannah venoms was assessed by 1D SDS-PAGE (Figure 1). The Coomassie Brilliant blue staining showed varying electrophoretic (molecular weight-distribution) protein profiles. This is not surprising as the venom proteins belong to different snake families, one a Viperidae (C. rhodostoma) and the other an Elapidae (O. hannah). The separations revealed that both venoms were composed of highly heterogeneous proteins varying in terms of band intensity as well as migration. Although staining was observed in the entirety of the gel background for both venoms, O. hannah venom appears to have more large molecular weight proteins compared to C. rhodostoma venom. C. rhodostoma venom proteins were generally <60 kDa whereas O. hannah venom proteins appeared to be concentrated in several regions, namely 70 kDa, 25 kDa, and 10 kDa.
The crude venom protein patterns observed in both gels show strong similarity with electrophoretic profiles from previous studies using Malaysian species [22,24]. A visual representation of the crude venom protein pattern provides an indication of possible geographic and species variation as well as illustrates the complexity of the sample being investigated.

2.2. Gel Filtration Separation of Crude Snake Venom Using Superdex G-75

For the GF-ESI-LC-MS/MS based approach, C. rhodostoma and O. hannah venoms were subjected to a molecular weight separation according to Section 4.2.4. The proteins were eluted based on their molecular sizes, beginning with large molecular weight proteins. C. rhodostoma separation yielded six distinctive peaks and two minor fractions, i.e., a total of eight fractions (Figure 2). Similarly, O. hannah venom separation yielded six peaks and one minor fraction, with a total of seven fractions (Figure 3). All resulting fractions were pooled according to their individual peaks and subjected to in-solution tryptic digestion prior to characterization.

2.3. Snake Venom Proteomic Characterization Using Shotgun Approach (Shotgun-ESI-LC-MS/MS)

Using the shotgun-ESI-LC-MS/MS approach, a total of 47 proteins belonging to 12 snake venom protein families were identified in C. rhodostoma venom (Table 1). The major protein families categorized were c-type lectin/snaclec (CTL/snaclec), phosphodiesterase (PDE), PLA2, SVMP, 5′-nucleotidase (5′-NTD), and snake venom serine protease (SVSP), followed by minor flavin monoamine oxidase (FMO), phospholipase B (PLB), CRiSP, glutaminyl-peptide cyclotransferase (QPCT), Trysinogen, and Ankyrin protein families. The latter three proteins were solely identified using this approach. For O. hannah venom, 76 proteins were identified from 19 toxin families. Quantitatively, there were two prominent families that made up the majority of the total venom, 3FTx with 23 peptides, followed by SVMP with 21 peptides (Table 2). The minority of venom was made up of Kunitz-type, CRiSP, FMO, PDE, 5′-NTD, PLA2, nerve-growth factor (NGF), acetylcholinesterase (AChE), complement C3 homolog (CC3H), endonuclease, neprilysin, ohanin/vespryn, platelet-derived growth factor/vascular endothelial growth factor (PDGF/VEGF), and SVSP protein families. Some notable proteins identified included the presence of insulin growth factors, hepatocyte growth factor activator, and CTL/snaclec peptides in O. hannah venom, which will be discussed in detail below. 5′-NTD and neprilysin proteins were exclusively detected using this approach.

2.4. Snake Venom Proteomic Characterization Using Gel Filtration Protein Fractions (GF-ESI-LC-MS/MS)

All resulting peptides from digested individual chromatographic fractions were classified according to protein families as per the previous section. From eight chromatographic fractions obtained for C. rhodostoma venom, a total of 99 proteins belonging to 12 protein families were detected (Figure 2, Table 3). The C. rhodostoma protein families were predominantly comprised of SVMPs as well as SVSPs and PDEs, while others families such as FMO, CTL/snaclec, PLA2, PLB, NGF, 5′-NTD, endonuclease, and CRiSP, were represented by lesser numbers of peptides. Of interest, the aminopeptidase protein family was newly discovered in C. rhodostoma venom using this proteomic approach. Other proteins not identified using the shotgun technique included endonucleases and NGFs (Table 3). For O. hannah venom, 150 proteins from 18 major families were observed for the seven chromatographic fractions (Figure 3, Table 4). The OH venom families consisted of 3FTx, SVMP, PLA2, CRiSP, PDE, Complement C3, Kunitz-type, Insulin, SVGF, SVSP, ohanin/vespryn, FMO, endonuclease, CTL/snaclec, AChE, NGF, PDGF/VEGF, and PLB. Interestingly, the presence of the single PLB peptide in the venom was detected only using this proteomic approach.

2.5. Summary of Identified C. rhodostoma and O. hannah Snake Venom Proteins

To obtain an overview of the major protein classes in Malaysian C. rhodostoma and O. hannah snake venoms, we combined two proteomic approaches that have made it possible to detect a higher number of proteins than previously reported. Our proteomic analysis revealed 114 proteins, belonging to 15 major protein families, for C. rhodostoma venom, and 176 proteins, belonging to 20 protein families, for O. hannah venom (Figure 4, Figure 5 and Figure 6, Table 5). Amongst these identified protein families are several new, previously undiscovered proteins for these venoms. These include aminopeptidase, glutaminyl-peptide cyclotransferase, and Ankyrin repeat for C. rhodostoma venom and insulin, CTL/snaclec, and SVGF (hepatocyte growth factor activator, macrophage colony stimulating factor, tumor necrosis factor receptor) for O. hannah venom (Table 5).
Using the shotgun approach, a total of 47 and 76 proteins were detected in CR and OH venoms, whereby 15 and 26 proteins were exclusively identified using the shotgun approach. In comparison, 99 and 150 proteins were identified using the gel-filtration approach, with 67 and 100 proteins being detected exclusively with this approach (Figure 4). Although the most common and widely used liquid chromatography separation for venom proteome profiling is reverse-phase C18 HPLC, we have utilized a molecular weight based separation to try to overcome some potential difficulties in the elution of low abundant large molecular weight proteins [25,26]. It appears that prior protein separation, together with the use of the unconventional molecular weight based separation, does indeed allow for decomplexation of crude venom proteins that has enabled a higher number of proteins to be detected.
We were also interested in the identification of low abundant proteins that could be of biological interest as a potential source of active compounds or play an important role in snake envenoming. Thus, the shotgun approach was employed to target potentially indiscernible small molecular weight proteins, even by SDS-PAGE [27]. Looking at the migration and strong intensity of the protein bands in the lower molecular weight range (Figure 1), we believe that the shotgun approach provided the necessary complementary analysis. Thus, combined, the results of both techniques have facilitated a detailed proteomic view of the entire venom components (Table 5, Figure 5 and Figure 6).

2.5.1. Snake Venom Protein Classes in Both C. rhodostoma and O. hannah Venoms

Generally, there are some protein families that are commonly present in the venoms of all snake families, i.e., Viperidae, Atractaspididae; Elapidae, and Colubridae, while others are restricted to certain families [28]. For instance, viper venoms usually induce hemorrhagic, hypotensive, and inflammatory effects due to high concentrations of metalloproteases, serine proteases, and C-type lectins, whereas elapid venoms usually induce neurotoxic symptoms owing to large amounts of pre- or post-synaptic neurotoxins. While these venoms may display different venom presentations, we identified a total of 11 families that were shared between these two snake families (Table 5, Figure 5 and Figure 6). The toxin families between the two snakes are compared and discussed below.

SVMP

This class of proteins represents a major toxin family in viper venoms, including Vipernaes and Crotalinaes. In Elapidaes, the proteins make up smaller proportions of the venom. The snake venom metalloproteinases (SVMPs) are known to display a wide range of physiological activates, such as hemorrhagic, fibrinolytic, prothrombin activating, blood coagulation factor X activating, apoptotic, platelet aggregation inhibition, proinflammatory, and inactivation of blood serine proteinase inhibitors [29]. Based on their molecular weight differences and domain structure, the SVMPs are divided into three main classes: P-I, P-II, and P-III. Comparatively, P-I SVMPs are relatively small at 20–30 kDa and contain a single catalytic metalloproteinase (MP) domain. P-II SVMPs, with molecular weight between 30–60 kDa, consist of an MP domain and an additional disintegrin domain. P-III proteins are the largest at 60–100 kDa, composed of both the MP and disintegrin domains as well as a cysteine-rich domain. In general, the P-IIIs are considered the most potent of SVMPs, followed by the P-IIs and the least hemorrhagic P-Is. Interestingly, all classes of SVMPs are responsible for the classic Viperid envenoming effect, which is the ability to cause hemorrhage at the bite site [30,31]. Paradoxically, novel human therapeutics may be achieved from the bioactive peptides of SVMPs because of their haemorrhage, coagulopathy, and inflammatory responses. Several examples of SVMPs with clinical utility include Alfimeprase (A. contortrix contortrix), Integrelin (S. miliarius barbouri), Aggrastat (E. carinatus), and NN-PF3 (N. naja) [6]. In this study, we report 20% SVMPs for the venom of the viperid C. rhodostoma and 25% for the venom of the elapid O. hannah (Table 5, Figure 5 and Figure 6). Although the results appear contrasting, previous studies have shown significant amount of SVMPS in O. hannah with 12% (Balinese specimen) and 24% (Malaysian specimen) [22,23]. However, it is important to note that despite the high amounts of SVMP detected, it is not the primary toxin family identified in the Elapidae venom. In contrast, even though only 20% were detected in C. rhodostoma venom, this was the largest toxin family represented in this Viperidae venom.

SVSP

Our study revealed a distinct difference in the abundance of snake venom serine proteinases (SVSPs) identified in C. rhodostoma venom (12%) compared to O. hannah venom (1%) (Table 5, Figure 5 and Figure 6). This is not unexpected, as SVSPs are commonly found in viperid venoms, but only in relatively low amounts in some Elapid venoms [32]. Together with SVMPs, the SVSPs are proteolytic enzymes that play important roles during snake envenoming. This family of proteins is known to exert hemostatic effects on prey by acting on various components of the coagulation cascade, on the fibrinolytic and kallikrein–kinin systems and cells to produce an imbalance of their systems [33,34]. The SVSPs that have been biologically characterized have both, or either, fibrinogenolytic and fibrinolytic activities and based on that, they either promote or inhibit coagulation. Some SVSPs are also known as ‘thrombin-like’ enzymes or abbreviated to SVTLE (snake venom thrombin like enzymes) due to their ability to mimic the fibrinogenolytic activities of thrombin [35]. Thrombin-like enzymes, ancrod and ancrod-2, were discovered in C. rhodostoma venom in high abundance using both shotgun and gel-filtration methods. First purified in 1967, ancrod was described as a potential anticoagulant and thoroughly researched for many years due to its potential clinical utility [36]. In O. hannah venom, only two SVSPs were identified, Alpha and beta fibrinogenase OhS1 and Neuroendocrine convertase 1. The Alpha and beta fibrinogenase OhS1 was also identified by Tan (2015) in the Malaysian king cobra species and is purported to possess potent fibrinogenolytic and amidolytic activities without any haemorrhagic consequences. In venomous animals, the neuroendocrine convertase 1 protein has only been identified in spider and gila monster venoms and the transcript of cone snail venom-glands [37,38,39,40]. The function of this protein is to convert inactive proteins into bioactive components by cleaving after Lys-Arg or after a tetrabasic sequence (Arg-X-Arg/Lys-Arg) [41,42]. Findings postulate that ‘venom proprotein convertase’ is most likely involved in the activation of exendins in gila monster, protoxins in spider venom, and conotoxins in cone snails (cleavage after dibasic residues) [38,43]. It is suggested that proforms of peptide toxins are activated by proteolytic cleavage by venom proprotein convertase, such as neuroendocrine 1 [37].

CRiSP

Cysteine-rich secretory proteins (CRiSPs) are relatively small glycoproteins (molecular weight 20 to 30 kDa) with a high degree of amino acid sequence similarity and highly conserved specific pattern of 16 cysteine residues [44,45,46]. This family of proteins are widely distributed and have been established in mammals, Drosophila, Xenopus, Heloderma horridum horridum, Caenorhabditis elegans, and Conus textile. Some of their proposed biological functions are gamete fusion, sperm maturation, anti-microbial and matrix-degradation, sperm chemoattractant, Ca2+ and K+ channel blockers, longevity and stress resistance regulators, and protease-like activity [44]. In snakes, they are found in most families, including Viperidae, Elipidae, and Colubridaes, demonstrating Ca2+ channel blocker-like properties and cyclic nucleotide-gated channel blocker activities [47,48]. This study revealed 15 CRiSP proteins in O. hannah venom (9%), but only one CRiSP protein in C. rhodostoma venom (1%) (Table 5, Figure 5 and Figure 6). This finding coincides with previous studies, whereby individual acidic and basic fragments of the same protein were identified. The 15 CRiSPs in O. hannah venom showed sequence similarities to 12 different species (Table 2 and Table 4). The ophanin from O. hannah venom was first purified and characterized in 2003 and compared to nine other snake venom CRiSP proteins. Based on the phylogenetic tree constructed in that study, it was suggested that O. hannah is more closely related to the Viperidae snakes rather than Elapidae snakes owing to the similarity of the nucleotide sequence of ophanin with the Viperidae snakes [49].

Phospholipase A2

Phospholipase A2’s activates the hydrolysis of glycerophospholipids at the sn-2 position of the glycerol backbone to dissociate fatty acids and corresponding 1-acyl lysophospholipids [2,50]. Over the years, the toxic and pharmacological effects presented by these enzymes have been well-documented. Effects include neurotoxic, myotoxic, cardiotoxic, cytotoxic, hemolytic, hypotensive, platelet aggregation, anticoagulant, pro-inflammatory, edematogenic, and bactericidal activities [50,51]. Widespread in nature, PLA2 has been identified in venoms from all snake families; i.e., Colubridae, Elapidae, Viperidae, and Hydrophilae [32]. Here, we report relatively high amounts of PLA2 enzymes in both C. rhodostoma and O. hannah venoms with 10% and 5%, respectively (Table 5, Figure 5 and Figure 6). The number of identified PLA2 peptides in this study is noticeably higher than reported in previous studies [19,20,21,22,52,53]. Ali et al. (2013), Vejayan et al. (2014), and Tsai et al. (2001) reported the presence of both acidic and basic isoforms of PLA2 in C. rhodotoma venom, whereas Tang et al. (2016) primarily identified the basic form of the enzyme. Additionally, Tang et al. identified two PLA2 peptides in minute amounts with sequence similarity to N. kaouthia venom, one of which was an acidic fraction found in less than 0.01% abundance. In our study, we identified a total of 11 PLA2 peptides (both acidic and basic isoforms), discovering more than two times the number of PLA2 compared to the latest proteomic analysis [19]. PLA2 peptides from C. rhodostoma venom had similar sequences to PLA2 peptides from O. hannah, T. stejnegeri, C. atrox, T. sabahi, and C. horridus venoms. Tsai et al. and Tan et al. both reported that the major PLA2 isolated from the Malaysian king cobra is the acidic isoform of the enzyme [22,52]. Based on the protein score obtained in our study, we are in agreement with their findings. The additional PLA2 peptides from O. hannah venom showed similarity to sequences from C. rhodostoma, B. fasciatus, M. fulvius, and T. biscutatus venoms. However, the higher number of PLA2s identified for both venoms here provides evidence that intragenic variations may indeed be influenced by ontogenic factors as suggested by Tan et al. (2015), which in turn will have significant impacts on the outcome of snake envenoming in Malaysia.

PLB

Phospholipase type B (PLB) is another member of the phospholipase superfamily. However, compared to PLA2, information on this enzyme is relatively scarce. The lack of data could be due to the limited identification of this enzyme in snake venomics analyses and/or the minute amounts of the enzyme that have been detected in snake venoms [54]. PLB hydrolyzes monoacyl phosphatides by liberating free fatty acid and forming glycerophosphoryl derivatives. To date, phospholipase B activity has been reported in snakes by Doery and Pearson (1964), Van Deenen and Haas (1963), Mohamed (1969), Shiloah (1973), and Bernheimer (1987) [55,56,57,58,59]. One study presented hemolytic and cytotoxic activities in Pseudechis colletti venom, indicating possible biological effects of this enzyme. In mice, PLB causes myoglobinuria, a muscle destruction occurrence also exerted by PLA2. Interestingly, Shiloah (1973) also observed that purified N. naja PLA2 isoforms presented PLB activity at specific experimental conditions. This finding led to the suggestion that PLB activity found in snake venoms may be due to a dual action by PLA2 [5]. In this study, we discovered 6% of PLB peptides in C. rhodostoma venom and 1% in O. hannah venom (Table 5, Figure 5 and Figure 6). A previous proteome study for the Malaysian C. rhodostoma venom reported less than 1% of this enzyme, whereas the Malaysian O. hannah venom presented only one PLB peptide, which is in accord with our study. Although, PLB identification and quantification has been limited thus far, Tan (2015) reported that the enzyme is present in some vipers and pit vipers [54,60,61,62]. In a recent proteomic study by Zainal Abidin et al. (2017) this finding was confirmed with a total of 8% and 2% of PLB enzymes reported for Malaysian T. wagleri and C. purpureomaculatus vipers, respectively. Although the exact pathophysiological role of PLB in venom has not been elucidated, the recent increasing abundance discovered in viper venoms may suggest a more significant role in snake envenoming than first predicted.

CTL/snaclec

Classically, CTLs bind to calcium and sugar residues (hence, “C” calcium and lectin). However, the C-type lectin-like proteins of snake venom do not contain the calcium and sugar binding loop and, thus, do not have lectin activity [63,64]. To avoid confusion between classic C-type lectins and C-type lectin-like proteins, the latter was termed snaclec (snake venom c-type lectins). Snaclecs exhibit numerous biological activities, including anticoagulation, pro-coagulation, and platelet modulating (agonist and antagonist) activities [63]. The phenomenon of contrasting platelet activity in snake venom has been reported in several crotalids, such as H. hypnale, D. acutus, and C. rhodostoma species [19,65,66]. Nevertheless, CTLs and snaclecs are considered to be abundant components of snake venom, particularly in vipers. In elapids, only minor amounts (<2%) have been reported in about one third of the family. Here, we report the abundance of CTL/snaclecs in C. rhodostoma with 12% compared to 1% in O. hannah venom (Table 5, Figure 5 and Figure 6). In C. rhodostoma venom, we identified all four subunits of rhodocetin, the platelet aggregation inhibitor snaclec. We also identified both alpha and beta subunits of rhodocytin, the platelet aggregation inducer snaclec that was previously known as aggretin and rhodoaggretin. Previous studies have shown evidence of these proteins grouped as one of the major components in C. rhodostoma venom [19,21]. In contrast, Tan et al. (2015) reported that although CTL venom-gland transcripts from the Malaysian and Balinese O. hannah venom were discovered, the transcribed genes were not translated into proteins. In this study, we present the first report of snaclec proteins found in O. hannah venom. The two CTL/snaclec peptides discovered in O. hannah venom showed sequence similarity to the alpha and beta subunits of the platelet aggregation inhibitor, rhodocetin, from C. rhodostoma venom.

FMO (l-amino Acid Oxidase)

L-amino acid oxidases (LAAOs) belong to the flavin monoamine oxidase (FMO) family and are considered one of the best studied protein families from snake venom. Since the discovery of snake venom LAAOs, many reports have been published on various species of the enzyme. Generally known as flavoenzymes, LAAO catalyzes the stereospecific oxidative deamination of an l-amino acid substrate to produce an α-ketoacid equivalent as well as ammonia and hydrogen peroxide [67,68,69]. Although the pathophysiological role in snake venom has not been defined, researchers have taken advantage of this well-characterized enzyme to study the biological and pharmacological activities of LAAOs. Some of the reported biological effects include hemorrhage, edema, coagulation, and platelet aggregation, whereas pharmacological effects include antimicrobial (virus, bacteria, and parasite) and anti-cancer activities [70]. In a recent snake venom database review, it was reported that 91% of all crotalinae species contained LAAO enzymes and the maximum amount identified thus far is from the crotalinae R. cotiara, constituting 20% of the venom. Although important in elapidaes and viperinaes, the enzymes only make up approximately 6% of their venoms [32]. In the current study, we report 7% LAAOs in C. rhodostoma venom and 2% in O. hannah venom (Table 5, Figure 5 and Figure 6). Although, the relative abundance reported here seems to be much lesser than reported by Ponnudurai et al. (1994), this value is in agreement with Tang et al. (2016), who reported a little less than 7% of LAAO in C. rhodostoma venom [71]. The study suggested that the discrepancy between the reported quantities is probably due to the lack of sensitivity in detecting other proteins in the whole venom. As for O. hannah venom, Tan et al. (2015) reported approximately 6% of LAAO from only two peptides [22], which is the same number of peptides detected in the current study. The difference in abundance is probably due to different quantitative analyses employed to determine the relative abundance. The sequence similarity of O. hannah LAAO peptides detected was observed in C. rhodostoma and O. hannah venoms, and the similarity with C. rhodostoma LAAO peptides was seen in C. rhodostoma, B. schelegelii, T. stejnegeri, L. muta, and B. fasciatus venoms. Several amine oxidase peptides were also observed in both venoms, with three peptides found in C. rhodostoma venom and one peptide in O. hannah venom.

NGF

This proteome analysis revealed 4% and 1% of nerve growth factors (NGFs) in C. rhodostoma and O. hannah venoms, respectively (Table 5, Figure 5 and Figure 6). In C. rhodostoma venom, a total of four peptides were discovered showing sequence similarity to venoms from B. jararacussu, P. flavoviridis, and B. atrox, whereas two peptides were found in O. hannah venom, with similarity to venoms from W. aegyptia and O. hannah. Snake venom nerve growth factors were first discovered serendipitously by Cohen and Levi-Montalcini [72]. In an effort to break down the nucleic acid portion of a purified component from a growth promoting secretion, Agkistrodon piscivorus venom was shown to be a significant inducer of nerve growth [72,73]. Since its discovery, numerous snake venom nerve growth factors have been purified and characterized from various snake venoms [74]. Reports have revealed that sv-NGFs are approximately 13 to 35 kDa, consisting of 10 to 20% of carbohydrates and are of great therapeutic interest for the treatment of neurodegenerative disease [75]. Apart from its growth promoting activity, sv-NGF has been proposed to possess vasculotoxin-like behavior, potentially to render the venom injection site more susceptible to facilitate infiltration of venom components into target tissues [76]. While sv-NGF has been established in Viperidae, Crotalinae, and Elapidae families, it is shown to occur in minute amounts (i.e., 0.1–0.5%) in several snake venoms [76,77]. Our findings here show a relatively higher amount than previously reported for both venoms [19,22].

PDE

Although PDEs are considered ubiquitously present in snake venoms, reports have described lower levels of PDE activity in elapids compared to crotalids and viperids [78]. In accordance, our analysis revealed 16% and 6% of phosphodiesterase (PDE) proteins in C. rhodostoma and O. hannah venoms, respectively (Table 5, Figure 5 and Figure 6). This class of toxins, also termed exonucleases, catalyzes the hydrolysis of phosphodiester bonds to release 5′-mononucleotides from the 3′ end of polynucleotides to provide nucleotide substrates to facilitate the activity of other toxins, such as 5′-nucleotidases [5,78]. In 1963, Russell et al. observed mean arterial pressure reduction and locomotor depression using PDE fractions from several snake venoms. Apart from this evidence of their potential function and as well as their hydrolytic properties, there is little information on the other pharmacological activities of these toxins [79]. The PDE peptides identified from C. rhodostoma venom showed similarity to peptides from C. adamateus, C. horridus, and S. miliarius barbouri venoms (Table 1 and Table 3). Likewise, O. hannah PDE peptides displayed similarity to C. adamanteus and C. horridus, as well as O. okinavensis, M. fulvius, and M. tener PDE peptides (Table 2 and Table 4). Interestingly, only minor amounts of PDEs were observed in the other proteomic analysis of King cobra and Malayan pit viper venoms [19,22,23,80].

Endonuclease

It is important to clarify that denominating endonucleases from nucleases has been a somewhat ambiguous endeavor. This is because endonucleolytic activity is an inherent property of venom PDEs (exonucleases). Endonucleases are made up of DNases that hydrolyze DNA, as well as RNases that hydrolyze RNA. Whereas, exonucleases are made up of PDEs that have the ability to hydrolyze both DNA and RNA [81]. To facilitate differentiation, parameters, such as the pH optima and metal ion requirement, were considered. To date, although this approach has provided some clarity, describing exclusive nuclease (endo and exo) activity is still considered a challenge. In our study, endonucleases were categorized from both venoms. Albeit in fairly low amounts, there was 3% in C. rhodostoma venom and 1% in O. hannah venom (Table 5, Figure 5 and Figure 6). Two peptides showing sequence similarity with O. hannah venom were discovered for both venoms, with an additional DNase-like protein showing sequence similarity with C. adamanteus found for C. rhodostoma venom (Table 2, Table 3 and Table 4). Thus far, no endonuclease proteins have been identified for these Malaysian specimens.

Nucleotidase

5′-Nucleotidase (5′-NTD) are high molecular weight glycoproteins and metalloenzymes ranging from 70 kDa to 100 kDa [77,79]. Primarily, 5′-NTDs endogenously liberates purines by selectively hydrolyzing 5′-nucleotidase to nucleosides [78,82]. Recent studies have indicated that 5′-NTDs play a crucial role in snake envenoming strategies of prey immobilization (via hypotension and paralysis) and prey digestion. This is because purine is the central component in all these strategies, whereby injection of exogenous nucleosides via venom and release of endogenous nucleosides via prey tissue, exert synchronous effects to bring about death. Thus, the ability of 5′-NTDs to provide a steady supply of purines has labelled them a multitoxin [79]. Though their biological functions have not been exhaustively studied, these enzymes are responsible for anticoagulant effects and inhibition of platelet aggregation in several snake venoms [83,84,85,86]. While ubiquitously distributed in snake venoms, 5′-NTD activity is more prominent in viperid venoms than elapid venoms. Our study reports 5% of 5′-NTDs in viperid C. rhodostoma venom and 2% in elapid O. hannah venom (Table 5, Figure 5 and Figure 6). Interestingly, the 5′-NTDs found in O. hannah venom were mainly ecto-nucleotidases showing sequence similarity to M. tener and M. fulvius (Table 2). Ecto-enzymes, such as ecto-5′-nucleotidases, are essentially membrane-embedded enzymes, with their active sites positioned on the exterior surface of the cell [78].

2.5.2. Snake Venom Protein Classes in C. rhodostoma Venom

In our study, we have shown that O. hannah venom contains a relatively high abundance of disintegrins (SVMP) and minor amounts of CTLs, protein components that occur in almost all viper venoms. However, there were several protein families that were detected exclusively in the viper, C. rhodostoma venom, albeit in low abundance. These toxin families are described in the following section.

Low Abundance Proteins Only (<2%)

A total of three protein families with relative abundance of less than 2% each were discovered in C. rhodostoma venom. These proteins were from the aminopeptidase, glutaminyl-peptide cyclotransferase (QPCT), and Ankyrin (ANK; designated as ‘Others’) families (Table 5, Figure 5). The aminopeptidase protein family is considered one of several less well understood enzyme constituents of snake venom. Reports of aminopeptidase activity have been presented in elapidae and viperidae crude venoms as well as a crotalinae venom fraction [28]. However, although there has been evidence of aminopeptidase transcripts in several species of snake venom glands, aminopeptidase proteins have only been identified in three snake venom proteome studies [25,87,88]. Here, we report the proteomic identification of aminopeptidase peptides in the venom of crotalinae, C. rhodostoma for the first time (Table 3). The two identified peptides sequences showed similarity to the aminopeptidase sequence from C. horridus venom. Thus far, the role of aminopeptidases in snake venom has not been fully elucidated. However, analogous comparisons have been made with mammalian aminopeptidases, proposing their importance in regulating brain function and blood pressure [28,89]. Further suggestions include disrupting physiological processes though cleavage of oligopeptide N-terminals, degradation of host tissues, increasing permeability to venom components, or processing toxins components within the venom [28].
One glutaminyl-peptide cyclotransferase (QPCT) was detected in the C. rhodostoma proteome, constituting only 1% of the entire crude venom (Table 5, Figure 5). The peptide identified for this protein showed sequence similarity to QPCTs from the venom of another crotalinae, A. contortrix contortrix (Table 3). Remarkably, glutaminyl cyclases were found to be notably present in crotalines, but not viperine venoms. Studies that have identified this protein in snake venom proteomes include other crotalinaes, C. durissus terrifcus, C. atrox, and C. purpureomaculatus [32]. Interestingly, Zainal Abidin et al. (2016) reported the absence of this protein in another crotalinae, T. wagleri. However, the proteome of T. wagleri was reported to be somewhat atypical compared to other viperid species [32,62]. Additionally, transcripts have been reported in the venom glands of more crotalinaes, S. catenus edwardsii and B. jararaca [90,91,92,93]. QPCTs are considered to be another poorly understood protein component in snake venoms. Very few QPCTs have been purified or isolated and therefore not thoroughly evaluated [94]. Nevertheless, it has been suggested that QPCTs play a role in post-translational modification (PTM), whereby they catalyze the N-terminal pyroglutamate (pGlu) portion of toxins to release ammonia or water molecules. Atypical to toxins, QPCTs may be involved in toxin maturation, protection from exopeptidase degradation, and aiding accurate protein conformation [95].
One of the low abundance peptides identified in C. rhodostoma venom identified an unusual protein domain designated as Ankyrin repeat-containing protein showing sequence similarity to A. contortrix contortrix venom (Table 1, Figure 5). Interestingly, the Ankyrin protein was initially described as a 33 amino acid repeat found amongst a few proteins in Drosophila melanogaster, a common fruit fly. Later, it was termed a cytoskeletal protein that was made up almost entirely of these short repeats [96]. Since then, the Ankyrin (ANK) repeats have been found to occur in thousands of proteins in viruses, prokaryotes, and eukaryotes. No particular function has been attributed to the ANK repeats, however, there are suggestions of their involvement in mediating protein-protein interactions [97]. In venomous animals, the most prominent work on Ankyrin repeats has been reported in a black widow spider, Latrodectus tredecimguttatus, transcriptomics study. The Ankyrin repeats were grouped as a superfamily consisting of five subfamilies (mainly α-latrotoxins), which differed in the number and distribution of ANK domains present in individual toxin families. As the ANK repeats were predominantly abundant in the main neurotoxic proteins identified in the venom, ANKs were included as part of the neurotoxin family [98]. However, looking at the low protein score (Table S1) coupled with the identification of only a single ANK peptide, we believe that the protein does not play a significant role in the venom and may serve in other general non-toxic functions.

2.5.3. Snake Venom Protein Classes in O. hannah Venom

One of the key differences between C. rhodostoma venom and O. hannah venom was the absence of 3FTx in the viper venom. This toxin family was identified as the most abundant protein family in the elapid venom. Other major protein families that were found exclusively in O. hannah venom were complement C3 homolog, Kunitz-type, and insulin, while several other low abundance families were also discovered in the venom, notably SVGF, hepatocyte growth factor (Table 5). The nine protein families only present in O. hannah venom are discussed below.

3FTx

Three-finger toxins (3FTx) are a class of proteins that belong to the snake venom superfamily, with a conserved three finger like appearance in three dimensional structures. 3FTx are commonly found in venoms of elapid snakes and are known to be a major contributing lethal factor to Elapidae and also Hydrophidae venoms [99,100,101]. As expected, we found a total of 55 proteins (31%) belonging to this family from the Malaysian O. hannah venom, but not C. rhodostoma venom (Table 5, Figure 6). The majority of the 3FTxs detected belonged to the long-chain subfamily (65%), whereas the rest were appropriated to the short-chain and non-conventional subfamilies. All identified peptides showed similar sequences to Elapidae venoms; N. kaouthia, N. oxiana, N. annulata annulata, N. haje haje, N. nivea, N. melanoleuca, N. naja, N. sputatrix, D. polylepis polylepis, D. viridis, B. multicinctus, B. flaviceps, A. superbus, and A. labialis (Table 2 and Table 4). Interestingly, there was a sequence similarity detected with a Viperidae venom, C. rhombeatus, of a short chain 3FTx, neurotoxin-like protein 1. Evolution has allowed snakes to adapt to a variety of prey and although previously thought to be unique to elapidaes, studies have shown the presence of elapidae components across snake families, including a non-venomous C. radiatus gland secretion [102].

Complement C3 Homolog

The complement C3 homolog family consists of proteins interacting with components of the serum complement system that modulate the immune system. In venom, this group is a structural and functional analog of complement component C3b, the activated form of C3. The cobra venom factor (CVF) is one of the best studied proteins in this family [103]. It is a glycoprotein, ~150 kDa and is frequently composed of alpha, beta, and gamma chains. When introduced into the bloodstream, it stimulates complement activation and results in the consumption of the complement activity. The CVF protein is known to be specific to the cobra species of the Elapidae family. Although the function of the purified CVF has been demonstrated in vivo, the exact biological role of CVF in cobra venom is still not fully understood [103,104]. In our findings, we have categorized seven proteins (4%) in this family consisting of complement-depleting factor, Ophiophagus venom factor, and venom factors (Table 2 and Table 4, Figure 6). Despite being annotated as a venom factor of a cobra, the purified protein is known to be non-toxic [103].

Kunitz-Type

Kunitz-type is a family of toxins consisting of a Kunitz domain with a serine protease inhibiting function. The characteristic domains are made up of a conserved fold of about 60 amino acids and three disulfide bridges. The family is also known as the bovine pancreatic trypsin inhibitor (BPTI) and is naturally present across the animal taxa [105,106]. Similar to 3FTx, this toxin family has evolved due to gene duplication and positive selection. Originally, the role of Kunitz type toxins was to inhibit the function of serine protease, but, over the years, some have evolved to block ion channels, such as K+ and Ca2+. However, a few of these toxins have the novel ability of reserving their protease inhibitory function as well as ion channel modulatory function, and therefore acquiring a bifunctional status [107]. Here, we have identified six ‘original’ Kunitz-type proteins appearing to have retained their protease inhibitory function. The Kunitz-type protein made up 3% of the total venom (Table 2, Table 4 and Table 5, Figure 6). Although previous studies have isolated these serine protease inhibitors from Elapidae and Viperidae venoms [108], there were no peptides that matched the identity of a Kunitz-type protein for Malaysian C. rhodostoma venom.

Insulin

Over the years, the advancement of investigational techniques in protein chemistry have propelled the recognition of insulin as being a part of a large protein superfamily that is present in various types of organisms. This superfamily consists of insulin and related peptides, such as insulin-like peptides and insulin-like growth factors as well as relaxins [109]. One very exciting discovery is the first example of insulin in a venom and furthermore an invertebrate of a cone snail, Conus geographus [110]. In snake venom, particularly O. hannah, there has been evidence at the transcriptional level of venom glands [22,23] and proteomic level [80]. However, due to the lack of evidence in other snake venoms, the transcriptomic findings were classified as rare and thus, in need of more thorough or detailed investigation. Here, we would like to document the presence of five insulin-like growth factor 1 proteins identified for the first time at the proteomic level of Malaysian O. hannah venom, making up 3% of the total venom (Table 2 and Table 4, Figure 6). Previously, the lack of identification at the proteomic level was attributed to extremely minute amounts of proteins in the venom. Therefore, we believe that our discovery is noteworthy as almost all except one protein was identified through the concentrated molecular weight separation fractions compared to the crude shotgun method.

Low Abundant Proteins in O. hannah Venom (<2%)

Snake venom acetylcholinesterase (AChE) was first reported in 1938 when strong choline esterase activity was observed in three cobra species. Together with butyrylcholinesterase (BuChE), AChE makes up the snake venom cholinesterase family and functions to hydrolyze choline esters [111]. Studies have shown that Bungarus fasciatus venom contains high concentrations of AChE at approximately 8 mg/g of dry venom with strong activity reported at around 750,000 Ellman units/mg, designating this venom as the richest source of this enzyme. Interestingly, AChE activity has only been demonstrated in Elapidae venoms, but not in Viperidae or Crotalidae venoms [112]. Not surprisingly, we detected only two peptides belonging to the AChE family, one showing sequence similarity with B. fasciatus and the other a carboxylic ester hydrolase from S. fasciata (Table 2 and Table 4). No AChE peptides were recovered from C. rhodostoma venom. Neprilysin is a zinc-metalloendopeptidase with specificity to the bacterial zinc peptidase, thermolysin. Originally located in the rabbit kidney, the first observation of the protein in snake was in the form of a transcript in the saw scaled viper venom of E. pyramidum leakeyi [113,114]. Other neprilysin transcripts were reported in the Balinese and Malaysian king cobra. The latter study also reported the first proteomic evidence of neprilysin in Malaysian O. hannah venom [22]. In this study, we report the presence of one neprilysin (Table 2) compared to two identified by the other Malaysian sample. In mammals, neprilysin has been implicated in neuropeptide regulation, however, its role in venom has not been elucidated [115].
Another low abundant, but noteworthy, toxin family in O. hannah venom is the Ohanin/Vespryn family. This family is a relatively new family of snake venom proteins established in 2005. Ohanin, a novel protein with unusual molecular mass, was isolated from King cobra venom that displayed an uncharacteristically low content of cysteine residues (1%). Comparison of amino acid sequence revealed a 93% sequence similarity with Thai cobrin, an isoform from the monocled cobra, Naja kaouthia. However, to date, there is no literature available on Thai cobrin [116]. The proteins possess the PRY and SPRY domain (B30.2-like domain) sequences that are the trademark for vespryn proteins, but do not share sequence similarity to other proteins in the snake venom database. This feature, coupled with the low cysteine residue content, makes both proteins exclusive members of this toxin family. Ohanin has been synthesized as a precursor named prohanin, with a propeptide segment at the C-terminal segment that completes the SPRY domain [117]. Prohanin (THA 903) has shown strong pharmacological potential and is currently undergoing clinical trials for the treatment of chronic pain [118]. In this study, we identified ohanin from O. hannah venom and its isoform, Thaicobrin, from the N. kaouthia venom (Table 2 and Table 4), revealing a total of 1% of this protein family in the overall venom composition (Table 4, Figure 6).
The proteomic analysis of Malaysian king cobra in this study also discovered single peptides of snake venom growth factors (SVGFs): Hepatocyte-growth factor activator, macrophages colony-stimulating factor (MCSF), tumor necrosis factor receptor (grouped as SVGF), and one other peptide, PDGF/VEGF growth factor (Table 2 and Table 4, Figure 6). Historically, growth factors represent growth and cell proliferation. Since then, the attributes have evolved to include neutrophins, such as cell differentiation, transformation, synthesis secretions, motility, and death [119]. Over the years, hundreds of snake venom growth factors (SVGFs) have been isolated, studied, and deposited in the UniProt databank. This includes the abovementioned platelet-derived growth factor (PDGF), vascular endothelial growth factor (VEGF), tumor necrosis factor (TNF), hepatocyte growth factor-tyrosine kinase (HGF), insulin-like growth factor-binding protein (IGF), and nerve growth factor (NGF) that has been discussed throughout this proteomic analysis [77]. Remarkably, prior to this study, the hepatocyte growth factor activator protein was only identified at the proteomic level in Thai O. hannah [80]. The investigations of O. hannah species from Malaysia and Bali did not report the presence of this peptide in their proteomic analyses [21,22,23]. In this study, the hepatocyte growth factor protein was discovered using both proteomic approaches and constitutes approximately 1% of the entire venom (Table 2 and Table 4, Figure 6). The hepatocyte growth factor is known to stimulate development and migration of vascular endothelial cells. This particular growth factor has demonstrated accelerated glomerular restoration via capillary endothelial cell growth and capillary regeneration in experimental progressive glomerulonephritis [120]. The other two growth factors (TNF and MCSF) were also not identified in the Malaysian, Thai, or Balinese specimens [21,22,23,80]. However, the presence of VEGF proteins was reported in the Balinese and Thai proteomes [23,80].
The focus of this study was to provide comprehensive proteomic profiles of Malaysian CR and OH snakes by maximally identifying the proteins present within their venom composition. We utilized venom samples from three separate pools (species) from the northern region of Malaysia, which provides a robust statistical framework for general venom composition identification. Nonetheless, further work, including a broader area of sampling around Malaysia, is necessary for a more detailed appreciation of individual differences between species.
Finally, this venomic investigation has provided evidence of major venom proteins, including SVMP, PDE, SVSP, CTL/snaclec, PLA2, FMO, PLB, 5′-NTD, NGF, Endonuclease, and CRiSP, as well as minor venom proteins consisting of Aminopeptidase, QPCT, Trypsinogen, and Ankyrin repeats in C. rhodostoma venom. Likewise, major protein components from 3FTx, SVMP, CRiSP, PDE, PLA2, Complement C3, Kunitz-type, and Insulin, along with minor components from SVGF, FMO, 5′-NTD, CTL/snaclec, SVSP, AChE, Ohanin/Vespryn, NGF, Endonuclease, Neprilysin, PDGF/VEGF, and PLB were determined in O. hannah venom. At this point of the study, it is difficult to distinguish if the small numbers of low abundance proteins found in the venom are indeed novel species-specific proteins or rare proteins that are yet to be discovered in other snake species. Nevertheless, the data obtained from our global proteomic analysis has revealed the presence of inherent medically and pharmacologically important components within Malaysian C. rhodostoma and O. hannah crude venoms, possessing potential clinical utility as tools or as prototypes for drug design studies.

3. Conclusions

Using a combination of two proteomic approaches, shotgun and gel-filtration fractionation techniques, we have provided a holistic view of the complex Malaysian C. rhodostoma and O. hannah crude venom protein compositions. As a forethought, we would like to highlight that proteins from the insulin, CTL/snaclec, and SVGFs families were identified for the first time in O. hannah venom in a proteomics study. Proteins from the Aminopeptidase, QPCT, and Ankyrin families in C. rhodostoma venom were also reported for the first time. In addition, the total number of proteins identified in both venoms is higher than reported in previous proteomic studies. Looking at the discrepancies between previous studies and our study, this finding is perhaps not exhaustive, but definitely supplementary, to previous findings and data collection. Nonetheless, we reiterate that the discovery of these new proteins will significantly impact the knowledge of the venom components in both these species, which may be used as a reference for improved antivenom production, understanding pathological effects of envenoming, and the discovery of biologically active peptides with medical and biotechnological values.

4. Materials and Methods

4.1. Materials

4.1.1. Snake Venom

Calloselasma rhodostoma and Ophiophagus hannah snake venoms were milked from adult snakes originating from the northern part of the west coast of Peninsular Malaysia (Perlis). Crude venoms were obtained from Mr. Zainuddin Ismail, a private snake enthusiast from Perlis. Manual milking was performed by placing the snake’s fangs on a parafilm-covered sterile container. Extracted venom samples were transported on ice to Monash University Malaysia (Sunway Campus). The venom samples were immediately frozen at −80 °C, lyophilized via a freeze-dryer, and stored at −20 °C until required. Quantification of the crude protein was carried out by performing a Bicinchoninic acid (BCA) assay as outlined by the manufacturer (Pierce Biotechnology, Rockford, IL, USA), using bovine serum albumin (BSA) as the standard.

4.1.2. Columns, Drugs and Chemicals

The Superdex G-75 (10 mm × 300 mm) and Hi-trap G-25 (16 mm × 25 mm) were purchased from GE-Healthcare (Uppsala, Sweden), whereas the Zorbax 300SB-C18 (75 μm × 150 mm, 5 μm) was purchased from Agilent (Santa Clara, CA, USA). The chemicals and reagents used were of analytical grade, except LCMS grade materials from Fisher Scientific (Leicestershire, UK) for all LCMS work. Water used in this study was ultrapure.

4.2. Methods

4.2.1. Sodium Dodedyl Sulfate-Polyacrylamide Gel Electrophoresis (SDS-PAGE)

Protein profiles of C. rhodostoma and O. hannah crude venom samples were assessed by SDS-PAGE according to Laemmli, 1970 [121]. Electrophoresis runs were carried out under reducing conditions in 12% resolving gels and 4% stacking gels in room temperature. Fractions were boiled at 95 °C for 5 min and allowed to cool to room temperature before being loaded onto the gel and run for approximately 60 min at 200 V in 1× SDS running buffer. The gels were fixed for 10 min in fixing solution and visualized through staining with Coomassie Blue R-350 (1 h or overnight) and de-staining with Coomassie de-stain (overnight). The molecular weights of protein bands were estimated using broad range SDS-PAGE molecular weight standards from Thermo Scientific (Rockford, AZ, USA).

4.2.2. In-solution Protein Digestion

Crude venom solutions (Section 4.2.3) and gel-filtration protein fractions (Section 4.2.4) were prepared for identification using an in-solution digestion procedure according to the manufacturer’s instructions (Agilent Technologies, Santa Clara, CA, USA). Briefly, 25 μL of 100 mM ammonium bicarbonate (ABC), 25 μL of trifluoroethanol (TFE), and 1 μL of dithiothreitol (DTT) were added to concentrated protein crude venom samples and fractions. The mixture was briefly vortexed and incubated at 60 °C for 60 min, followed by incubation with 4 μL of iodoacetamide at 37 °C for another 60 min. 1 μL of DTT was added to the solution and incubated again for another 60 min. 100 mM ABC and trypsin [1:50 (w/w) enzyme-to-substrate ratio; approximately 500 ng] were then added to the solution and incubated at 37 °C for 18 h. At the end of the incubation, trypsin activity was terminated by adding formic acid and finally subjected to centrifugal vacuum evaporation to remove solvents and concentrate the sample simultaneously. Samples were monitored for 1 to 2 days until completely dry and subsequently stored at −20 °C until LC-MS/MS analysis.

4.2.3. Identification of Crude Venom Proteins Using Shotgun LCMS/MS

Previously frozen lyophilized C. rhodostoma and O. hannah venoms were adjusted to 200 μg protein samples and subjected to in-solution tryptic digest protocol as described in Section 4.2.2. The crude venom samples from three separately pooled stocks were utilized for this study.

4.2.4. Fractionation of C. rhodostoma and O. hannah Crude Venoms by Superdex G-75 Gel Filtration Chromatography

Freeze-dried crude venoms (10 mg) were dissolved in 0.01 M ammonium acetate (250 μL; pH 7.0) and centrifuged at 5000× g for 1 min to remove insoluble material. The clarified supernatant was subjected to gel filtration chromatography via Superdex G-75 pre-equilibrated and eluted with the same buffer at a flow rate of 0.5 mL/min. The fractionation was monitored using a fast protein liquid chromatography system (AKTA Prime Plus, GE Healthcare, Uppsala, Sweden) and carried out at room temperature. The chromatographic data were acquired and processed using Unicorn 5.20 Workstation software (GE Healthcare, Uppsala, Sweden). All fractions were pooled and subjected to buffer exchange using a Hi-Trap G-25 desalting column. Desalted fractions were concentrated and stored at −20 °C until required for mass spectrometry identification analysis. Samples from three separately pooled stocks were adjusted to 200 μg before being subjected to in-solution tryptic digest.

4.2.5. LC-MS/MS Analysis of Crude Venom and Gel-Filtration Protein Fractions

Tryptic-digested peptides stored in −20 °C were reconstituted in 20 μL of 0.1% formic acid in water. 1 μL of peptide sample was injected on a nano C-18 column (see below) using the Agilent 1200 capillary nanoflow liquid chromatography system coupled to a 6550 iFunnel Q-TOF LC/MS mass spectrometer (Agilent) via a nano-electrospray (Nano-ESI) source. Digested peptide products were separated using the C18 300 Å Large Capacity Chip column (Zorbax 300SB-C18) that was pre-equilibrated with 0.1% formic acid in water (A) and eluted in a multi-gradient approach using 90% acetonitrile in 0.1% formic acid in water (B) at a flow rate of 4 μL/min. The method utilized a 47 min gradient starting with a 30 min ramping from 3% to 50% solution B, followed by a further increase of 50% to 95% solution B for 2 min, 95% solution B for 7 min, and finally a decrease from 95% to 3% solution B for another 8 min. BSA was used as a standard for reference. Mass spectrometry (MS/MS) was performed using these following settings: Gas temperature 300 °C; drying gas 5 L/min; capillary voltage of 2050 V; and fragmenter voltage 300 V. The auto MS/MS mode was used, with positive ion polarity setup. Acquisition ranges were 110–3000 m/z for MS scan and 50–3000 m/z for MS/MS scan.

4.2.6. Data Analysis

Data analysis and protein identification by automated de novo sequencing was performed using PEAKS 8.0 software (Bioinformatics Solutions Inc., Waterloo, ON, Canada). MS/MS spectra were searched against the Swiss-Prot database (July 2017) by selecting all proteins from the Serpentes taxa subjected to the following parameters: Digestion Trypsin; maximum number of missed cleavages 3; fixed modification: carbamidomethylation; mass error tolerance (parent and fragment) 0.1 Da; false discovery rate (FDR) 1%; unique peptides 2; and peptide hit threshold (−10lgP) 20. Venom components were classified according to protein families and their relative abundances were calculated according to Zainal Abidin et al. (2016) [62] as follows:
number   of   proteins   ( protein   family ) total   proteins   detected   using   LC − MS / MS × 100 %
Relative abundance was expressed as percentages and represented in pie charts. A complete list of identified peptides and proteins is included as Supplementary Material (Tables S1–S4).

Supplementary Materials

The following are available online at https://www.mdpi.com/2072-6651/10/11/434/s1, Table S1: Complete mass spectrophotometry data of Malaysian Calloselasma rhodostoma (Malayan pit viper) crude venom via shotgun-ESI-LC-MS/MS analysis, Table S2: Complete mass spectrophotometry data of Malaysian Ophiophagus hannah (King cobra) crude venom via shotgun-ESI-LC-MS/MS analysis, Table S3: Complete mass spectrophotometry data of Malaysian Calloselasma rhodostoma (Malayan pit viper) crude venom via GF-ESI-LC-MS/MS analysis, Table S4: Complete mass spectrophotometry data of Malaysian Ophiophagus hannah (King cobra) crude venom via GF-ESI-LC-MS/MS analysis.

Author Contributions

S.K. was responsible for the design and execution of the experiments, data analysis and writing of the manuscript. I.O. conceived the experiments and provided critical feedback. S.S.H. provided critical feedback. W.C.H. analyzed the data, wrote and provided critical feedback for the manuscript.

Funding

W.C.H. is funded by a National Health and Medical Research Council (NHMRC)—Australia, Centres for Research Excellence Grant ID: 1110343.

Acknowledgments

The authors would like to acknowledge Jeffrey Cheah School of Medicine and Health Sciences, Monash University Malaysia for research facilities and support. We would also like to thank Nurziana Sharmilla, Tee Ting Yee and Adillah Akhasan from Jeffrey Cheah School of Medicine and Health Sciences Proteomics laboratory for assistance with operation of the LC-MS systems.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Williams, D.; Gutiérrez, J.M.; Harrison, R.; Warrell, D.A.; White, J.; Winkel, K.D.; Gopalakrishnakone, P. The global snake bite initiative: An antidote for snake bite. Lancet 2010, 375, 89–91. [Google Scholar] [CrossRef]
  2. Koh, D.; Armugam, A.; Jeyaseelan, K. Snake venom components and their applications in biomedicine. Cell. Mol. Life Sci. CMLS 2006, 63, 3030–3041. [Google Scholar] [CrossRef] [PubMed]
  3. Kasai, K.; Ishikawa, T.; Nakamura, T.; Miura, T. Antibacterial properties of l-amino acid oxidase: Mechanisms of action and perspectives for therapeutic applications. Appl. Microbiol. Biotechnol. 2015, 99, 7847–7857. [Google Scholar] [CrossRef] [PubMed]
  4. Gomes, A.; Bhattacharjee, P.; Mishra, R.; Biswas, A.K.; Dasgupta, S.C.; Giri, B.; Debnath, A.; Gupta, S.D.; Das, T.; Gomes, A. Anticancer potential of animal venoms and toxins. Indian J. Exp. Boil. 2010, 48, 93–103. [Google Scholar]
  5. Iwanaga, S.; Suzuki, T. Enzymes in snake venom. In Snake Venoms; Springer: Berlin/Heidelberg, Germany, 1979; pp. 61–158. [Google Scholar]
  6. Chellapandi, P. Structural, functional and therapeutic aspects of snake venom metalloproteinases. Mini-Rev. Org. Chem. 2014, 11, 28–44. [Google Scholar] [CrossRef]
  7. King, G.F. Venoms as a platform for human drugs: Translating toxins into therapeutics. Expert Opin. Boil. Ther. 2011, 11, 1469–1484. [Google Scholar] [CrossRef] [PubMed]
  8. Oldrati, V.; Arrell, M.; Violette, A.; Perret, F.; Sprüngli, X.; Wolfender, J.-L.; Stöcklin, R. Advances in venomics. Mol. Biosyst. 2016, 12, 3530–3543. [Google Scholar] [CrossRef] [PubMed]
  9. Calvete, J.J. Snake venomics: From the inventory of toxins to biology. Toxicon 2013, 75, 44–62. [Google Scholar] [CrossRef] [PubMed]
  10. De Lima, D.C.; Alvarez Abreu, P.; de Freitas, C.C.; Santos, D.O.; Borges, R.O.; Dos Santos, T.C.; Mendes Cabral, L.; Rodrigues, C.R.; Castro, H.C. Snake venom: Any clue for antibiotics and CAM? Evid.-Based Complement. Altern. Med. 2005, 2, 39–47. [Google Scholar] [CrossRef] [PubMed]
  11. Das, I.; Ahmed, N.; Liat, L.B. Venomous terrestrial snakes of Malaysia: Their identity and biology. In Clinical Toxinology in Asia Pacific and Africa; Gopalakrishnakone, P., Faiz, M.A., Fernando, R., Gnanathasan, C.A., Habib, A.G., Yang, C.C., Eds.; Springer: Dordrecht, The Netherlands, 2015. [Google Scholar]
  12. Nolan, C.; Hall, L.; Barlow, G. [17] Ancrod, the coagulating enzyme from Malayan pit viper (Agkistrodon rhodostoma) venom. Methods Enzymol. 1976, 45, 205–213. [Google Scholar] [PubMed]
  13. Marsh, N.; Williams, V. Practical applications of snake venom toxins in haemostasis. Toxicon 2005, 45, 1171–1181. [Google Scholar] [CrossRef] [PubMed]
  14. McCleary, R.J.; Kang, T.S.; Kini, R.M. Reptile venoms as a platform for drug development. In Venoms to Drugs: Venom as a Source for the Development of Human Therapeutics, Chapter: Reptile Venoms as a Platform for Drug Development; Royal Society of Chemistry: Cambridge, UK, 2015; pp. 129–162. [Google Scholar]
  15. 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]
  16. Pu, X.; Wong, P.; Gopalakrishnakone, P. A novel analgesic toxin (hannalgesin) from the venom of king cobra (Ophiophagus hannah). Toxicon 1995, 33, 1425–1431. [Google Scholar] [CrossRef]
  17. Koh, C.Y.; Kini, R.M. From snake venom toxins to therapeutics–cardiovascular examples. Toxicon 2012, 59, 497–506. [Google Scholar] [CrossRef] [PubMed]
  18. Harvey, A.L. Toxins and drug discovery. Toxicon 2014, 92, 193–200. [Google Scholar] [CrossRef] [PubMed]
  19. Tang, E.L.; Tan, C.H.; Fung, S.Y.; Tan, N.H. Venomics of Calloselasma rhodostoma, the Malayan pit viper: A complex toxin arsenal unraveled. J. Proteom. 2016, 148, 44–56. [Google Scholar] [CrossRef] [PubMed]
  20. Ali, S.A.; Baumann, K.; Jackson, T.N.; Wood, K.; Mason, S.; Undheim, E.A.; Nouwens, A.; Koludarov, I.; Hendrikx, I.; Jones, A.; et al. Proteomic comparison of Hypnale hypnale (Hump-Nosed Pit-Viper) and Calloselasma rhodostoma (Malayan Pit-Viper) venoms. J. Proteom. 2013, 91, 338–343. [Google Scholar] [CrossRef] [PubMed]
  21. Vejayan, J.; Khoon, T.L.; Ibrahim, H. Comparative analysis of the venom proteome of four important Malaysian snake species. J. Venom. Anim. Toxins Incl. Trop. Dis. 2014, 20, 6. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Tan, C.H.; Tan, K.Y.; Fung, S.Y.; Tan, N.H. Venom-gland transcriptome and venom proteome of the Malaysian king cobra (Ophiophagus hannah). BMC Genom. 2015, 16, 687. [Google Scholar] [CrossRef] [PubMed]
  23. Vonk, F.J.; Casewell, N.R.; Henkel, C.V.; Heimberg, A.M.; Jansen, H.J.; McCleary, R.J.; Kerkkamp, H.M.; Vos, R.A.; Guerreiro, I.; Calvete, J.J.; et al. The king cobra genome reveals dynamic gene evolution and adaptation in the snake venom system. Proc. Natl. Acad. Sci. USA 2013, 110, 20651–20656. [Google Scholar] [CrossRef] [PubMed]
  24. Khow, O.; Chulasugandha, P.; Pakmanee, N. Venom protein of the haematotoxic snakes Cryptelytrops albolabris, Calloselasma rhodostoma, and Daboia russelii siamensis. Sci. Asia 2011, 37, 377–381. [Google Scholar] [CrossRef]
  25. Choudhury, M.; McCleary, R.J.; 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]
  26. Shuting, L.; Jingqiang, W.; Zhang, X.; Yan, R.; Ning, W.; Kang, Z.; Xishu, C.; Caifeng, Z.; Xiaolei, L.; Jianmin, S.; et al. Proteomic characterization of two snake venoms: Naja naja atra and Agkistrodon halys. Biochem. J. 2004, 384, 119–127. [Google Scholar]
  27. De Oliveira, I.S.; Cardoso, I.A.; Bordon, K.D.C.F.; Carone, S.E.I.; Boldrini-França, J.; Pucca, M.B.; Zoccal, K.F.; Faccioli, L.H.; Sampaio, S.V.; Rosa, J.C.; et al. Global proteomic and functional analysis of Crotalus durissus collilineatus individual venom variation and its impact on envenoming. J. Proteom. 2018, 2018. [Google Scholar] [CrossRef] [PubMed]
  28. Vaiyapuri, S.; Wagstaff, S.C.; Watson, K.A.; Harrison, R.A.; Gibbins, J.M.; Hutchinson, E.G. Purification and functional characterisation of rhiminopeptidase A, a novel aminopeptidase from the venom of Bitis gabonica rhinoceros. PLoS Negl. Trop. Dis. 2010, 4, e796. [Google Scholar] [CrossRef] [PubMed]
  29. Markland, F.S., Jr.; Swenson, S. Snake venom metalloproteinases. Toxicon 2013, 62, 3–18. [Google Scholar] [CrossRef] [PubMed]
  30. Fox, J.W.; Serrano, S.M. Timeline of key events in snake venom metalloproteinase research. J. Proteom. 2009, 72, 200–209. [Google Scholar] [CrossRef] [PubMed]
  31. Fox, J.W.; Serrano, S.M. Structural considerations of the snake venom metalloproteinases, key members of the M12 reprolysin family of metalloproteinases. Toxicon 2005, 45, 969–985. [Google Scholar] [CrossRef] [PubMed]
  32. Tasoulis, T.; Isbister, G.K. A review and database of snake venom proteomes. Toxins 2017, 9, 290. [Google Scholar] [CrossRef] [PubMed]
  33. Serrano, S.M. The long road of research on snake venom serine proteinases. Toxicon 2013, 62, 19–26. [Google Scholar] [CrossRef] [PubMed]
  34. Serrano, S.M.; Maroun, R.C. Snake venom serine proteinases: Sequence homology vs. Substrate specificity, a paradox to be solved. Toxicon 2005, 45, 1115–1132. [Google Scholar] [CrossRef] [PubMed]
  35. Stubbs, M.T.; Bode, W. A player of many parts: The spotlight falls on thrombin’s structure. Thromb. Res. 1993, 69, 1–58. [Google Scholar] [CrossRef]
  36. Esnouf, M.; Tunnah, G. The isolation and properties of the thrombin-like activity from Ancistrodon rhodostoma venom. Br. J. Haematol. 1967, 13, 581–590. [Google Scholar] [CrossRef] [PubMed]
  37. Sanggaard, K.W.; Dyrlund, T.F.; Thomsen, L.R.; Nielsen, T.A.; Brøndum, L.; Wang, T.; Thøgersen, I.B.; Enghild, J.J. Characterization of the gila monster (Heloderma suspectum suspectum) venom proteome. J. Proteom. 2015, 117, 1–11. [Google Scholar] [CrossRef] [PubMed]
  38. Sanggaard, K.W.; Bechsgaard, J.S.; Fang, X.; Duan, J.; Dyrlund, T.F.; Gupta, V.; Jiang, X.; Cheng, L.; Fan, D.; Feng, Y.; et al. Spider genomes provide insight into composition and evolution of venom and silk. Nat. Commun. 2014, 5, 3765. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Safavi-Hemami, H.; Hu, H.; Gorasia, D.G.; Bandyopadhyay, P.K.; Veith, P.D.; Young, N.D.; Reynolds, E.C.; Yandell, M.; Olivera, B.M.; Purcell, A.W. Combined proteomic and transcriptomic interrogation of the venom gland of Conus geographus uncovers novel components and functional compartmentalization. Mol. Cell. Proteom. 2014, 13, 938–953. [Google Scholar] [CrossRef] [PubMed]
  40. Robinson, S.D.; Safavi-Hemami, H.; McIntosh, L.D.; Purcell, A.W.; Norton, R.S.; Papenfuss, A.T. Diversity of conotoxin gene superfamilies in the venomous snail, Conus victoriae. PLoS ONE 2014, 9, e87648. [Google Scholar] [CrossRef] [PubMed]
  41. Hook, V.; Yasothornsrikul, S.; Greenbaum, D.; Medzihradszky, K.F.; Troutner, K.; Toneff, T.; Bundey, R.; Logrinova, A.; Reinheckel, T.; Peters, C.; et al. Cathepsin L and Arg/Lys aminopeptidase: A distinct prohormone processing pathway for the biosynthesis of peptide neurotransmitters and hormones. Boil. Chem. 2004, 385, 473–480. [Google Scholar] [CrossRef] [PubMed]
  42. Jackson, R.S.; Creemers, J.W.; Ohagi, S.; Raffin-Sanson, M.-L.; Sanders, L.; Montague, C.T.; Hutton, J.C.; O’Rahilly, S. Obesity and impaired prohormone processing associated with mutations in the human prohormone convertase 1 gene. Nat. Genet. 1997, 16, 303. [Google Scholar] [CrossRef] [PubMed]
  43. Milne, T.J.; Abbenante, G.; Tyndall, J.D.; Halliday, J.; Lewis, R.J. Isolation and characterization of a cone snail protease with homology to crisp proteins of the pathogenesis-related protein superfamily. J. Boil. Chem. 2003, 278, 31105–31110. [Google Scholar] [CrossRef] [PubMed]
  44. Yamazaki, Y.; Morita, T. Structure and function of snake venom cysteine-rich secretory proteins. Toxicon 2004, 44, 227–231. [Google Scholar] [CrossRef] [PubMed]
  45. Yamazaki, Y.; Morita, T. Snake venom components affecting blood coagulation and the vascular system: Structural similarities and marked diversity. Curr. Pharm. Des. 2007, 13, 2872–2886. [Google Scholar] [CrossRef] [PubMed]
  46. Sunagar, K.; Johnson, W.E.; O’Brien, S.J.; Vasconcelos, V.; Antunes, A. Evolution of crisps associated with toxicoferan-reptilian venom and mammalian reproduction. Mol. Boil. Evol. 2012, 29, 1807–1822. [Google Scholar] [CrossRef] [PubMed]
  47. Matsunaga, Y.; Yamazaki, Y.; Hyodo, F.; Sugiyama, Y.; Nozaki, M.; Morita, T. Structural divergence of cysteine-rich secretory proteins in snake venoms. J. Biochem. 2008, 145, 365–375. [Google Scholar] [CrossRef] [PubMed]
  48. Brown, R.L.; Haley, T.L.; West, K.A.; Crabb, J.W. Pseudechetoxin: A peptide blocker of cyclic nucleotide-gated ion channels. Proc. Natl. Acad. Sci. USA 1999, 96, 754–759. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Yamazaki, Y.; Hyodo, F.; Morita, T. Wide distribution of cysteine-rich secretory proteins in snake venoms: Isolation and cloning of novel snake venom cysteine-rich secretory proteins. Arch. Biochem. Biophys. 2003, 412, 133–141. [Google Scholar] [CrossRef]
  50. Kini, R.M. Excitement ahead: Structure, function and mechanism of snake venom phospholipase A2 enzymes. Toxicon 2003, 42, 827–840. [Google Scholar] [CrossRef] [PubMed]
  51. Gutiérrez, J.M.; Lomonte, B. Phospholipases A2: Unveiling the secrets of a functionally versatile group of snake venom toxins. Toxicon 2013, 62, 27–39. [Google Scholar] [CrossRef] [PubMed]
  52. Tsai, I.H.; Wang, Y.M.; Au, L.C.; Ko, T.P.; Chen, Y.H.; Chu, Y.F. Phospholipases A2 from Callosellasma rhodostoma venom gland. FEBS J. 2000, 267, 6684–6691. [Google Scholar] [CrossRef]
  53. Tsai, S.-H.; Chen, Y.-C.; Chen, L.; Wang, Y.-M.; Tsai, I.-H. Binding of a venom Lys-49 phospholipase A2 to LPS and suppression of its effects on mouse macrophages. Toxicon 2007, 50, 914–922. [Google Scholar] [CrossRef] [PubMed]
  54. Tan, N.H.; Fung, S.Y.; Tan, K.Y.; Yap, M.K.K.; Gnanathasan, C.A.; Tan, C.H. Functional venomics of the Sri Lankan Russell’s viper (Daboia russelii) and its toxinological correlations. J. Proteom. 2015, 128, 403–423. [Google Scholar] [CrossRef] [PubMed]
  55. Doery, H.M.; Pearson, J.E. Phospholipase B in snake venoms and bee venom. Biochem. J. 1964, 92, 599. [Google Scholar] [CrossRef] [PubMed]
  56. Van Deenen, L.L.M.; de Haas, G.H. The substrate specificity of phospholipase A. Biochim. Biophys. Acta (BBA) Spéc. Sect. Lipids Relat. Subj. 1963, 70, 538–553. [Google Scholar]
  57. Mohamed, A.; Kamel, A.; Ayobe, M. Studies of phospholipase A and B activities of Egyptian snake venoms and a scorpion toxin. Toxicon 1969, 6, 293–298. [Google Scholar] [CrossRef]
  58. Shiloah, J.; Klibansky, C.; De Vries, A.; Berger, A. Phospholipase B activity of a purified phospholipase A from Vipera palestinae venom. J. Lipid Res. 1973, 14, 267–278. [Google Scholar] [PubMed]
  59. Bernheimer, A.; Linder, R.; Weinstein, S.; Kim, K.-S. Isolation and characterization of a phospholipase B from venom of collett’s snake, Pseudechis colletti. Toxicon 1987, 25, 547–554. [Google Scholar] [CrossRef]
  60. Kovalchuk, S.I.; Ziganshin, R.H.; Starkov, V.G.; Tsetlin, V.I.; Utkin, Y.N. Quantitative proteomic analysis of venoms from Russian vipers of Pelias group: Phospholipases A2 are the main venom components. Toxins 2016, 8, 105. [Google Scholar] [CrossRef] [PubMed]
  61. Ziganshin, R.H.; Kovalchuk, S.I.; Arapidi, G.P.; Starkov, V.G.; Hoang, A.N.; Nguyen, T.T.T.; Nguyen, K.C.; Shoibonov, B.B.; Tsetlin, V.I.; Utkin, Y.N. Quantitative proteomic analysis of Vietnamese krait venoms: Neurotoxins are the major components in Bungarus multicinctus and phospholipases A2 in Bungarus fasciatus. Toxicon 2015, 107, 197–209. [Google Scholar] [CrossRef] [PubMed]
  62. Zainal Abidin, S.A.; Rajadurai, P.; Chowdhury, M.E.H.; Ahmad Rusmili, M.R.; Othman, I.; Naidu, R. Proteomic characterization and comparison of malaysian Tropidolaemus wagleri and Cryptelytrops purpureomaculatus venom using shotgun-proteomics. Toxins 2016, 8, 299. [Google Scholar] [CrossRef] [PubMed]
  63. Morita, T. Structures and functions of snake venom CLPs (C-type lectin-like proteins) with anticoagulant-, procoagulant-, and platelet-modulating activities. Toxicon 2005, 45, 1099–1114. [Google Scholar] [CrossRef] [PubMed]
  64. Clemetson, K.J. Snaclecs (snake C-type lectins) that inhibit or activate platelets by binding to receptors. Toxicon 2010, 56, 1236–1246. [Google Scholar] [CrossRef] [PubMed]
  65. Qinghua, L.; Xiaowei, Z.; Wei, Y.; Chenji, L.; Yijun, H.; Pengxin, Q.; Xingwen, S.; Songnian, H.; Guangmei, Y. A catalog for transcripts in the venom gland of the Agkistrodon acutus: Identification of the toxins potentially involved in coagulopathy. Biochem. Biophys. Res. Commun. 2006, 341, 522–531. [Google Scholar] [CrossRef] [PubMed]
  66. Tan, C.H.; Tan, N.H.; Sim, S.M.; Fung, S.Y.; Gnanathasan, C.A. Proteomic investigation of Sri Lankan hump-nosed pit viper (Hypnale hypnale) venom. Toxicon 2015, 93, 164–170. [Google Scholar] [CrossRef] [PubMed]
  67. Du, X.-Y.; Clemetson, K.J. Snake venom l-amino acid oxidases. Toxicon 2002, 40, 659–665. [Google Scholar] [CrossRef]
  68. Fung, S.; Tan, N. Snake venom l-amino acid oxidases. In Handbook of Venoms and Toxins of Reptiles; CRC Press Inc.: London, UK, 2009; pp. 221–235. [Google Scholar]
  69. Baby, J.; Rajan Sheeja, S.; Jeevitha, M.; Ajisha, S. Pharmacological effects of snake venom l-amino acid oxidases. Int. J. Res. Ayurveda Pharm. 2011, 2, 114–120. [Google Scholar]
  70. Fry, B. Venomous Reptiles and Their Toxins: Evolution, Pathophysiology and Biodiscovery; Oxford University Press: Oxford, UK, 2015. [Google Scholar]
  71. Ponnudurai, G.; Chung, M.C.; Tan, N.-H. Purification and properties of the l-amino acid oxidase from Malayan pit viper (Calloselasma rhodostoma) venom. Arch. Biochem. Biophys. 1994, 313, 373–378. [Google Scholar] [CrossRef] [PubMed]
  72. Cohen, S.; Levi-Montalcini, R. A nerve growth-stimulating factor isolated from snake venom. Proc. Natl. Acad. Sci. USA 1956, 42, 571–574. [Google Scholar] [CrossRef] [PubMed]
  73. McCleary, R.J.; Kini, R.M. Non-enzymatic proteins from snake venoms: A gold mine of pharmacological tools and drug leads. Toxicon 2013, 62, 56–74. [Google Scholar] [CrossRef] [PubMed]
  74. Trummal, K.; Tõnismägi, K.; Paalme, V.; Järvekülg, L.; Siigur, J.; Siigur, E. Molecular diversity of snake venom nerve growth factors. Toxicon 2011, 58, 363–368. [Google Scholar] [CrossRef] [PubMed]
  75. Earl, S.T.; Birrell, G.W.; Wallis, T.P.; St Pierre, L.D.; Masci, P.P.; de Jersey, J.; Gorman, J.J.; Lavin, M.F. Post-translational modification accounts for the presence of varied forms of nerve growth factor in Australian elapid snake venoms. Proteomics 2006, 6, 6554–6565. [Google Scholar] [CrossRef] [PubMed]
  76. Kostiza, T.; Meier, J. Nerve growth factors from snake venoms: Chemical properties, mode of action and biological significance. Toxicon 1996, 34, 787–806. [Google Scholar] [CrossRef]
  77. Boldrini-França, J.; Cologna, C.T.; Pucca, M.B.; Bordon, K.D.C.F.; Amorim, F.G.; Anjolette, F.A.P.; Cordeiro, F.A.; Wiezel, G.A.; Cerni, F.A.; Pinheiro-Junior, E.L.; et al. Minor snake venom proteins: Structure, function and potential applications. Biochim. Biophys. Acta (BBA) Gen. Subj. 2017, 1861, 824–838. [Google Scholar] [CrossRef] [PubMed]
  78. Aird, S.D. Ophidian envenomation strategies and the role of purines. Toxicon 2002, 40, 335–393. [Google Scholar] [CrossRef]
  79. Dhananjaya, B.L.; D’souza, C.J. The pharmacological role of nucleotidases in snake venoms. Cell Biochem. Funct. 2010, 28, 171–177. [Google Scholar] [CrossRef] [PubMed]
  80. Danpaiboon, W.; Reamtong, O.; Sookrung, N.; Seesuay, W.; Sakolvaree, Y.; Thanongsaksrikul, J.; Dong-din-on, F.; Srimanote, P.; Thueng-in, K.; Chaicumpa, W. Ophiophagus hannah venom: Proteome, components bound by Naja kaouthia antivenin and neutralization by N. kaouthia neurotoxin-specific human ScFv. Toxins 2014, 6, 1526–1558. [Google Scholar] [CrossRef] [PubMed]
  81. Dhananjaya, B.; D’souza, C.J. An overview on nucleases (DNase, RNase, and phosphodiesterase) in snake venoms. Biochemistry. 2010, 75, 1–6. [Google Scholar] [CrossRef] [PubMed]
  82. Aird, S.D. Taxonomic distribution and quantitative analysis of free purine and pyrimidine nucleosides in snake venoms. Comp. Biochem. Physiol. Part B Biochem. Mol. Boil. 2005, 140, 109–126. [Google Scholar] [CrossRef] [PubMed]
  83. Ouyang, C.; Huang, T.-F. Inhibition of platelet aggregation by 5′-nucleotidase purified from Trimeresurus gramineus snake venom. Toxicon 1983, 21, 491–501. [Google Scholar] [CrossRef]
  84. Ouyang, C.; Huang, T.-F. Platelet aggregation inhibitors from Agkistrodon acutus snake venom. Toxicon 1986, 24, 1099–1106. [Google Scholar] [CrossRef]
  85. Boffa, M.-C.; Boffa, G.A. Correlations between the enzymatic activities and the factors active on blood coagulation and platelet aggregation from the venom of Vipera aspis. Biochim. Biophys. Acta (BBA) Gen. Subj. 1974, 354, 275–290. [Google Scholar] [CrossRef]
  86. Dhananjaya, B.; Nataraju, A.; Rajesh, R.; Gowda, C.R.; Sharath, B.; Vishwanath, B.; D’Souza, C.J. Anticoagulant effect of Naja naja venom 5′ nucleotidase: Demonstration through the use of novel specific inhibitor, vanillic acid. Toxicon 2006, 48, 411–421. [Google Scholar] [CrossRef] [PubMed]
  87. Yap, M.K.; Fung, S.Y.; Tan, K.Y.; Tan, N.H. Proteomic characterization of venom of the medically important Southeast Asian Naja sumatrana (Equatorial spitting cobra). Acta Trop. 2014, 133, 15–25. [Google Scholar] [CrossRef] [PubMed]
  88. Melani, R.D.; Araujo, G.D.; Carvalho, P.C.; Goto, L.; Nogueira, F.C.; Junqueira, M.; Domont, G.B. Seeing beyond the tip of the iceberg: A deep analysis of the venome of the Brazilian Rattlesnake, Crotalus durissus terrificus. EuPA Open Proteom. 2015, 8, 144–156. [Google Scholar] [CrossRef]
  89. Aird, S.D.; Watanabe, Y.; Villar-Briones, A.; Roy, M.C.; Terada, K.; Mikheyev, A.S. Quantitative high-throughput profiling of snake venom gland transcriptomes and proteomes (Ovophis okinavensis and Protobothrops flavoviridis). BMC Genom. 2013, 14, 790. [Google Scholar] [CrossRef] [PubMed]
  90. Georgieva, D.; Öhler, M.; Seifert, J.; Von Bergen, M.; Arni, R.K.; Genov, N.; Betzel, C. Snake venomic of Crotalus durissus terrificus—Correlation with pharmacological activities. J. Proteome Res. 2010, 9, 2302–2316. [Google Scholar] [CrossRef] [PubMed]
  91. Fox, J.W.; Ma, L.; Nelson, K.; Sherman, N.E.; Serrano, S.M. Comparison of indirect and direct approaches using ion-trap and Fourier transform ion cyclotron resonance mass spectrometry for exploring viperid venom proteomes. Toxicon 2006, 47, 700–714. [Google Scholar] [CrossRef] [PubMed]
  92. Chapeaurouge, A.; Reza, M.A.; Mackessy, S.P.; Carvalho, P.C.; Valente, R.H.; Teixeira-Ferreira, A.; Perales, J.; Lin, Q.; Kini, R.M. Interrogating the venom of the viperid snake Sistrurus catenatus edwardsii by a combined approach of electrospray and MALDI mass spectrometry. PLoS ONE 2015, 10, e0092091. [Google Scholar] [CrossRef] [PubMed]
  93. Gonçalves-Machado, L.; Pla, D.; Sanz, L.; Jorge, R.J.B.; Leitão-De-Araújo, M.; Alves, M.L.M.; Alvares, D.J.; De Miranda, J.; Nowatzki, J.; de Morais-Zani, K.; et al. Combined venomics, venom gland transcriptomics, bioactivities, and antivenomics of two Bothrops jararaca populations from geographic isolated regions within the Brazilian Atlantic rainforest. J. Proteom. 2016, 135, 73–89. [Google Scholar] [CrossRef] [PubMed]
  94. Wang, Y.-M.; Huang, K.-F.; Tsai, I.-H. Snake venom glutaminyl cyclases: Purification, cloning, kinetic study, recombinant expression, and comparison with the human enzyme. Toxicon 2014, 86, 40–50. [Google Scholar] [CrossRef] [PubMed]
  95. Pawlak, J.; Kini, R.M. Snake venom glutaminyl cyclase. Toxicon 2006, 48, 278–286. [Google Scholar] [CrossRef] [PubMed]
  96. Breeden, L.; Nasmyth, K. Similarity between cell-cycle genes of budding yeast and fission yeast and the Notch gene of Drosophila. Nature 1987, 329, 651–654. [Google Scholar] [CrossRef] [PubMed]
  97. Foord, R.; Taylor, I.A.; Sedgwick, S.G.; Smerdon, S.J. X-ray structural analysis of the yeast cell cycle regulator Swi6 reveals variations of the ankyrin fold and has implications for Swi6 function. Nat. Struct. Mol. Boil. 1999, 6, 157–165. [Google Scholar]
  98. He, Q.; Duan, Z.; Yu, Y.; Liu, Z.; Liu, Z.; Liang, S. The venom gland transcriptome of Latrodectus tredecimguttatus revealed by deep sequencing and cDNA library analysis. PLoS ONE 2013, 8, e81357. [Google Scholar] [CrossRef] [PubMed]
  99. Kini, R.M. Molecular moulds with multiple missions: Functional sites in three-finger toxins. Clin. Exp. Pharmacol. Physiol. 2002, 29, 815–822. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  100. Kini, R.M.; Doley, R. Structure, function and evolution of three-finger toxins: Mini proteins with multiple targets. Toxicon 2010, 56, 855–867. [Google Scholar] [CrossRef] [PubMed]
  101. Dufton, M.; Hider, R. Structure and pharmacology of elapid cytotoxins. Pharmacol. Ther. 1988, 36, 1–40. [Google Scholar] [CrossRef]
  102. Fry, B.G.; Lumsden, N.G.; Wüster, W.; Wickramaratna, J.C.; Hodgson, W.C.; Kini, R.M. Isolation of a neurotoxin (α-colubritoxin) from a nonvenomous colubrid: Evidence for early origin of venom in snakes. J. Mol. Evol. 2003, 57, 446–452. [Google Scholar] [CrossRef] [PubMed]
  103. Vogel, C.-W.; Bredehorst, R.; Fritzinger, D.C.; Grunwald, T.; Ziegelmüller, P.; Kock, M.A. Structure and function of cobra venom factor, the complement-activating protein in cobra venom. In Natural Toxins 2; Springer: Boston, MA, USA, 1996; pp. 97–114. [Google Scholar]
  104. Vogel, C.-W.; Fritzinger, D.C. Cobra venom factor: Structure, function, and humanization for therapeutic complement depletion. Toxicon 2010, 56, 1198–1222. [Google Scholar] [CrossRef] [PubMed]
  105. Ranasinghe, S.; McManus, D.P. Structure and function of invertebrate Kunitz serine protease inhibitors. Dev. Comp. Immunol. 2013, 39, 219–227. [Google Scholar] [CrossRef] [PubMed]
  106. Župunski, V.; Kordiš, D.; Gubenšek, F. Adaptive evolution in the snake venom Kunitz/BPTI protein family. FEBS Lett. 2003, 547, 131–136. [Google Scholar] [CrossRef] [Green Version]
  107. Yang, W.; Feng, J.; Wang, B.; Cao, Z.; Li, W.; Wu, Y.; Chen, Z. BF9, the first functionally characterized snake toxin peptide with Kunitz-type protease and potassium channel inhibiting properties. J. Biochem. Mol. Toxicol. 2014, 28, 76–83. [Google Scholar] [CrossRef] [PubMed]
  108. Shafqat, J.; Beg, O.U.; Yin, S.J.; Zaidi, Z.H.; Jörnvall, H. Primary structure and functional properties of cobra (Naja naja naja) venom Kunitz-type trypsin inhibitor. FEBS J. 1990, 194, 337–341. [Google Scholar] [CrossRef]
  109. Shabanpoor, F.; Separovic, F.; Wade, J.D. The human insulin superfamily of polypeptide hormones. Vitam. Horm. 2009, 80, 1–31. [Google Scholar] [PubMed]
  110. Safavi-Hemami, H.; Gajewiak, J.; Karanth, S.; Robinson, S.D.; Ueberheide, B.; Douglass, A.D.; Schlegel, A.; Imperial, J.S.; Watkins, M.; Bandyopadhyay, P.K.; et al. Specialized insulin is used for chemical warfare by fish-hunting cone snails. Proc. Natl. Acad. Sci. USA 2015, 112, 1743–1748. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  111. Iyengar, N.; Sehra, K.; Mukerji, B.; Chopra, R. Choline esterase in cobra venom. Curr. Sci. 1938, 7, 51–53. [Google Scholar]
  112. Frobert, Y.; Créminon, C.; Cousin, X.; Rémy, M.-H.; Chatel, J.-M.; Bon, S.; Bon, C.; Grassi, J. Acetylcholinesterases from Elapidae snake venoms: Biochemical, immunological and enzymatic characterization. Biochim. Biophys. Acta (BBA) Protein Struct. Mol. Enzym. 1997, 1339, 253–267. [Google Scholar] [CrossRef]
  113. Kerr, M.; Kenny, A. The purification and specificity of a neutral endopeptidase from rabbit kidney brush border. Biochem. J. 1974, 137, 477–488. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  114. Casewell, N.R.; Harrison, R.A.; Wüster, W.; Wagstaff, S.C. Comparative venom gland transcriptome surveys of the saw-scaled vipers (Viperidae: Echis) reveal substantial intra-family gene diversity and novel venom transcripts. BMC Genom. 2009, 10, 564. [Google Scholar] [CrossRef] [PubMed]
  115. Turner, A.J.; Isaac, R.E.; Coates, D. The neprilysin (NEP) family of zinc metalloendopeptidases: Genomics and function. Bioessays 2001, 23, 261–269. [Google Scholar] [CrossRef]
  116. Pung, Y.F.; Wong, P.T.; Kumar, P.P.; Hodgson, W.C.; Kini, R.M. Ohanin, a novel protein from king cobra venom, induces hypolocomotion and hyperalgesia in mice. J. Boil. Chem. 2005, 280, 13137–13147. [Google Scholar] [CrossRef] [PubMed]
  117. Pung, Y.F.; Kumar, S.V.; Rajagopalan, N.; Fry, B.G.; Kumar, P.P.; Kini, R.M. Ohanin, a novel protein from king cobra venom: Its cDNA and genomic organization. Gene 2006, 371, 246–256. [Google Scholar] [CrossRef] [PubMed]
  118. Gopalakrishnakone, P.; Pu, X.C.; Wong, P.T.-H.; Gwee, M.C.E.; Kini, R.M. Therapeutic Molecules Derived from Snake Venom. U.S. Patent 6613745B1, 2 September 2003. [Google Scholar]
  119. Wordinger, R.J.; Clark, A.F. Growth factors and neurotrophic factors as targets. Ocul. Ther. 2008, 87–116. [Google Scholar] [CrossRef]
  120. Mori, T.; Shimizu, A.; Masuda, Y.; Fukuda, Y.; Yamanaka, N. Hepatocyte growth factor-stimulating endothelial cell growth and accelerating glomerular capillary repair in experimental progressive glomerulonephritis. Nephron Exp. Nephrol. 2003, 94, e44–e54. [Google Scholar] [CrossRef] [PubMed]
  121. Laemmli, J. SDS-polyacrylamide gel electrophoresis of protein. Nature 1970, 2, 680–685. [Google Scholar] [CrossRef]
Toxins 10 00434 g001 550
Figure 1. One-dimension electrophoretic profile of (a) C. rhodostoma and (b) O. hannah venoms. 12% SDS-PAGE gel picture showing the individual venoms prepared and run in various concentrations. M represents the molecular weight standard in kDa. The lanes (1–9) contain 100, 80, 60, 40, 20, 10, 5, 2.5, and 1 μg of venom proteins, respectively.
Figure 1. One-dimension electrophoretic profile of (a) C. rhodostoma and (b) O. hannah venoms. 12% SDS-PAGE gel picture showing the individual venoms prepared and run in various concentrations. M represents the molecular weight standard in kDa. The lanes (1–9) contain 100, 80, 60, 40, 20, 10, 5, 2.5, and 1 μg of venom proteins, respectively.
Toxins 10 00434 g001
Toxins 10 00434 g002 550
Figure 2. Superdex G75 chromatographic profile of Malaysian C. rhodostoma venom. Venom (10 mg) was dissolved in 10 mM ammonium acetate (250 μL), pH 7.0 and loaded onto a pre-equilibrated column. Pre-equilibration and isocratic elution was carried out using the sample buffer. The chromatographic separation yielded six prominent peaks (1–6) and two minor fractions (marked *). Individual runs were repeated (at least 10 times). Fractions were pooled, desalted, concentrated, and subjected to in-solution tryptic digestion before analysis by ESI-LC-MS/MS.
Figure 2. Superdex G75 chromatographic profile of Malaysian C. rhodostoma venom. Venom (10 mg) was dissolved in 10 mM ammonium acetate (250 μL), pH 7.0 and loaded onto a pre-equilibrated column. Pre-equilibration and isocratic elution was carried out using the sample buffer. The chromatographic separation yielded six prominent peaks (1–6) and two minor fractions (marked *). Individual runs were repeated (at least 10 times). Fractions were pooled, desalted, concentrated, and subjected to in-solution tryptic digestion before analysis by ESI-LC-MS/MS.
Toxins 10 00434 g002
Toxins 10 00434 g003 550
Figure 3. Superdex G75 chromatographic profile of Malaysian O. hannah venom. Venom (10 mg) was dissolved in 10 mM ammonium acetate (250 μL), pH 7.0 and loaded onto a pre-equilibrated column. Pre-equilibration and isocratic elution was carried out using the sample buffer. The chromatographic separation yielded six prominent peaks (1–6) and one minor fraction (marked *). Individual runs were repeated (at least 10 times). Fractions were pooled, desalted, concentrated, and subjected to in-solution tryptic digestion before analysis by ESI-LC-MS/MS.
Figure 3. Superdex G75 chromatographic profile of Malaysian O. hannah venom. Venom (10 mg) was dissolved in 10 mM ammonium acetate (250 μL), pH 7.0 and loaded onto a pre-equilibrated column. Pre-equilibration and isocratic elution was carried out using the sample buffer. The chromatographic separation yielded six prominent peaks (1–6) and one minor fraction (marked *). Individual runs were repeated (at least 10 times). Fractions were pooled, desalted, concentrated, and subjected to in-solution tryptic digestion before analysis by ESI-LC-MS/MS.
Toxins 10 00434 g003
Toxins 10 00434 g004 550
Figure 4. A two-way Venn diagram depicting the overlap of the number of proteins identified by two proteomic strategies for (a) C. rhodostoma and (b) O. hannah venoms. The strategies employed crude venom and gel filtration fractions to comprehensively characterize the entire venom proteome. The proteomic analyses were performed via a combination of nanoflow liquid chromatography coupled with electrospray tandem mass spectrometry (ESI-LC-MS/MS) and automated de novo sequencing. 32 and 50 proteins were common between the two data sets for C. rhodostoma and O. hannah, respectively. Both venoms demonstrated a higher number of proteins identified in the gel filtration strategy compared to the crude venom strategy. The identified proteins were matched using the UniProt database (Serpentes) and subsequently categorized according to individual toxin family.
Figure 4. A two-way Venn diagram depicting the overlap of the number of proteins identified by two proteomic strategies for (a) C. rhodostoma and (b) O. hannah venoms. The strategies employed crude venom and gel filtration fractions to comprehensively characterize the entire venom proteome. The proteomic analyses were performed via a combination of nanoflow liquid chromatography coupled with electrospray tandem mass spectrometry (ESI-LC-MS/MS) and automated de novo sequencing. 32 and 50 proteins were common between the two data sets for C. rhodostoma and O. hannah, respectively. Both venoms demonstrated a higher number of proteins identified in the gel filtration strategy compared to the crude venom strategy. The identified proteins were matched using the UniProt database (Serpentes) and subsequently categorized according to individual toxin family.
Toxins 10 00434 g004
Toxins 10 00434 g005 550
Figure 5. The proteome of Malaysian Calloselasma rhodostoma venom. Pie chart represents the relative abundance (expressed in percentage) of protein families identified in the whole crude venom. A total of 15 toxin families were determined by a combination of nanoflow liquid chromatography coupled with electrospray tandem mass spectrometry (ESI-LC-MS/MS) and automated de novo sequencing. Abbreviations: CRiSP: Cysteine-rich secretory proteins; CTL: C-type lectin; FMO: Flavin monoamine oxidase; NGF: Nerve-growth factor; PDE: Phosphodiesterase; PLA2: Phospholipase A2; PLB: Phospholipase B; QPCT: Glutaminyl-peptide cyclotransferase; SVMP: Snake venom metalloproteinase; SVSP: Snake venom serine proteinase; ANK: Ankyrin repeats.
Figure 5. The proteome of Malaysian Calloselasma rhodostoma venom. Pie chart represents the relative abundance (expressed in percentage) of protein families identified in the whole crude venom. A total of 15 toxin families were determined by a combination of nanoflow liquid chromatography coupled with electrospray tandem mass spectrometry (ESI-LC-MS/MS) and automated de novo sequencing. Abbreviations: CRiSP: Cysteine-rich secretory proteins; CTL: C-type lectin; FMO: Flavin monoamine oxidase; NGF: Nerve-growth factor; PDE: Phosphodiesterase; PLA2: Phospholipase A2; PLB: Phospholipase B; QPCT: Glutaminyl-peptide cyclotransferase; SVMP: Snake venom metalloproteinase; SVSP: Snake venom serine proteinase; ANK: Ankyrin repeats.
Toxins 10 00434 g005
Toxins 10 00434 g006 550
Figure 6. The proteome of Malaysian Ophiophagus hannah venom. Pie chart represents the relative abundance (expressed in percentage) of protein families identified in the whole crude venom. A total of 20 toxin families were determined by a combination of nanoflow liquid chromatography coupled with electrospray tandem mass spectrometry (ESI-LC-MS/MS) and automated de novo sequencing. Abbreviations: 3FTx: Three-finger toxins, AChE: Acetylcholinesterase; CTL: C-type lectin; CRiSP: Cysteine-rich secretory proteins; FMO: Flavin monoamine oxidase; NGF: Nerve-growth factor; PDGF/VEGF: Platelet-derived growth factor/vascular endothelial growth factor; PDE: Phosphodiesterase; PLA2: Phospholipase A2; PLB: Phospholipase B; QPCT: Glutaminyl-peptide cyclotransferase; SVGF: Snake venom-growth factor; SVMP: Snake venom metalloproteinase; SVSP: Snake venom serine protease.
Figure 6. The proteome of Malaysian Ophiophagus hannah venom. Pie chart represents the relative abundance (expressed in percentage) of protein families identified in the whole crude venom. A total of 20 toxin families were determined by a combination of nanoflow liquid chromatography coupled with electrospray tandem mass spectrometry (ESI-LC-MS/MS) and automated de novo sequencing. Abbreviations: 3FTx: Three-finger toxins, AChE: Acetylcholinesterase; CTL: C-type lectin; CRiSP: Cysteine-rich secretory proteins; FMO: Flavin monoamine oxidase; NGF: Nerve-growth factor; PDGF/VEGF: Platelet-derived growth factor/vascular endothelial growth factor; PDE: Phosphodiesterase; PLA2: Phospholipase A2; PLB: Phospholipase B; QPCT: Glutaminyl-peptide cyclotransferase; SVGF: Snake venom-growth factor; SVMP: Snake venom metalloproteinase; SVSP: Snake venom serine protease.
Toxins 10 00434 g006
Table
Table 1. List of proteins identified from Malaysian C. rhodostoma venom using the shotgun-ESI-LC-MS/MS technique. Refer to Supplementary File (Table S1) for complete proteomic data and abbreviations lists.
Table 1. List of proteins identified from Malaysian C. rhodostoma venom using the shotgun-ESI-LC-MS/MS technique. Refer to Supplementary File (Table S1) for complete proteomic data and abbreviations lists.
Protein FamilyProtein DescriptionNo. of ProteinsProtein Accession No.No. of PeptidesNo. of Unique Peptides
5′-NTD5′ nucleotidase (P. flavoviridis)5T2HRS942
5-nucleotidase (P. elegans)A0A077L6M542
5′ nucleotidase (P. flavoviridis)A0A077L7M942
Snake venom 5′-nucleotidase (A. piscivorus)A0A194APL932
BATXNUC1 (B. atrox)A0A1L8D66732
CRiSPCysteine-rich secretory protein LCCL 2 (C. adamanteus)1A0A0F7ZEB622
CTL/snaclecC-type lectin (A. piscivorus leucostoma)12G8FML644
C-type lectin 2 (A. piscivorus)A0A194APP266
C-type lectin APL (A. piscivorus piscivorus)P0DM3644
C-type lectin BjL (B. jararaca)Q9PRY744
C-type lectin PAL (B. arietans)Q9PSN044
Snaclec clone 2100755 (D. acutus)Q8JIV833
Snaclec rhodocetin subunit alpha (C. rhodostoma)P813971715
Snaclec rhodocetin subunit beta (C. rhodostoma)P813981514
Snaclec rhodocetin subunit gamma (C. rhodostoma)D2YW3964
Snaclec rhodocetin subunit delta (C. rhodostoma)D2YW4066
Snaclec rhodocytin subunit alpha (C. rhodostoma)Q9I84188
Snaclec rhodocytin subunit beta (C. rhodostoma)Q9I8401413
FMOAmine oxidase (O. okinavensis)2T2HQ57182
L-amino-acid oxidase (C. rhodostoma)P813825334
PDECadam10_PDE-1 (C. adamanteus)7A0A1W7RBB062
Ectonucleotide pyrophosphatase/phosphodiesterase 3 (C. horridus)T1DJT562
Phosphodiesterase (C. horridus)T1D6P763
Phosphodiesterase (C. adamanteus)A0A0F7Z2Q363
Phosphodiesterase (S. miliarius barbouri)A0A194AS0262
Venom phosphodiesterase 1 (C. adamanteus)J3SEZ363
Venom phosphodiesterase 2 (C. adamanteus)J3SBP363
PLA2Acidic phospholipase A2 H1E6 (C. rhodostoma)6Q9PVF21110
Acidic phospholipase A2 S1E6-b (C. rhodostoma)Q9PVF065
Acidic phospholipase A2 Ts-A4 (T. stejnegeri)Q6H3C71110
Basic phospholipase A2 W6D49 (C. rhodostoma)Q9PVF4216
Phospholipase A2 (C. rhodostoma)A0A0H3U266236
Phospholipase A2 (T. sabahi)A0A0H3U23232
PLBPhospholipase B (O. okinavensis)2T2HQ7532
Phospholipase B (P. regius)A0A098LY7453
QPCTGlutaminyl-peptide cyclotransferase (A. contortrix contortrix)1A0A1W7RH8822
SVMPMetalloprotease, mRNA (G. intermedius)5A0A0C4ZNF152
Snake venom metalloproteinase kistomin (C. rhodostoma)P0CB148686
Zinc metalloproteinase/disintegrin (C. rhodostoma)P304033836
Zinc metalloproteinase-disintegrin-like halysase (G. halys) Q8AWI542
Zinc metalloproteinase-disintegrin-like HV1 (P. flavoviridis)Q90ZI342
SVSPSerine protease 1 (P. elegans)4A0A077LA4632
Thrombin-like enzyme ancrod (C. rhodostoma)P263242120
Thrombin-like enzyme ancrod-2 (C. rhodostoma)P477971713
Thrombin-like enzyme elegaxobin-2 (P. elegans)P8478732
TrypsinogenTrypsinogen homolog (B. jararaca)1Q9PUF333
OthersAnkyrin repeat-containing protein (A. contortrix contortrix)1A0A1W7RJF122
TOTAL PROTEINS47
Abbreviations: 5′-NTD: 5′-nucleotidase; CRiSP: Cysteine-rich secretory protein; CTL: C-type lectin; FMO: Flavin monoamine oxidase; PDE: Phosphodiesterase; PLA2: Phospholipase A2; PLB: Phospholipase B; QPCT: Glutaminyl-peptide cyclotransferase; SVMP: Snake venom metalloproteinase; SVSP: Snake venom serine protease.
Table
Table 2. List of proteins identified in Malaysian O. hannah venom via the shotgun-ESI-LC-MS/MS technique. Refer to Supplementary File (Table S2) for complete proteomic data and abbreviations lists.
Table 2. List of proteins identified in Malaysian O. hannah venom via the shotgun-ESI-LC-MS/MS technique. Refer to Supplementary File (Table S2) for complete proteomic data and abbreviations lists.
Protein FamilyProtein DescriptionNo. of ProteinsProtein Accession No.No. of PeptidesNo. of Unique Peptides
3FTxAlpha-elapitoxin-Dpp2d (D. polylepis polylepis)23C0HJD722
Beta-cardiotoxin CTX27 (O. hannah)Q69CK0113
Haditoxin (O. hannah)A8N28633
Long neurotoxin 1 (O. hannah)P0138782
Long neurotoxin 1 (N. nivea)P0139022
Long neurotoxin 2 (O. hannah)P0138672
Long neurotoxin 4 (O. hannah)P8015622
Long neurotoxin LNTX1 (O. hannah)Q2VBP842
Long neurotoxin LNTX-2 homolog (O. hannah)A8N28582
Long neurotoxin LNTX37 (O. hannah)Q2VBP332
Long neurotoxin OH-5 (O. hannah)P8096522
Long neurotoxin OH-37 (O. hannah)Q53B5977
Long neurotoxin OH-55 (O. hannah)Q53B58136
Long neurotoxin OH-56 (O. hannah)Q53B5744
Long neurotoxin OH-57 (O. hannah)Q53B5633
Neurotoxin-like protein 1 (C. rhombeatus)P8471622
Short neurotoxin OH-35 (O. hannah)Q53B4955
Weak toxin DE-1 (O. hannah)P014123629
Weak toxin DE-1 homolog 1 (O. hannah)Q69CJ822
Weak neurotoxin OH-72 (O. hannah)Q53B6143
Weak neurotoxin WNTX33 (O. hannah)Q2VBN32316
Weak neurotoxin WNTX34 (O. hannah)Q2VBN254
Uncharacterized protein (O. hannah)—Toxin activityV8N3G955
5′-NTDEcto-5′-nucleotidase (M. tener)3A0A194AS9832
Ecto-5′-nucleotidase 1 (M. fulvius)U3FYP932
Ecto-5′-nucleotidase 1a (M. fulvius)A0A0F7YZM632
AChEAcetylcholinesterase (B. fasciatus)1Q9203555
CC3HO. venom factor (O. hannah)1I2C09086
CRiSPCysteine-rich venom protein (T. stejnegeri)4P6062322
Cysteine-rich venom protein ablomin (G. blomhoffii)Q8JI4022
Cysteine-rich venom protein mossambin (N. mossambica)P0DL1622
Cysteine-rich venom protein ophanin (O. hannah)Q7ZT981312
CTL/snaclecSnaclec rhodocetin subunit beta (C. rhodostoma)1P8139822
EndonucleaseEndonuclease domain-containing 1 protein (O. hannah)1V8NCE322
FMOAmine oxidase (O. hannah)3V8N3Q9292
L-amino-acid oxidase (O. hannah)P813831616
L-amino-acid oxidase (C. rhodostoma)P8138287
InsulinInsulin-like growth factor 1 (P. guttatus)1I6TIH676
Kunitz-typeKunitz-type serine protease inhibitor TCI (O. hannah)5B6RLX21311
Kunitz-type serine protease inhibitor Vur-Kin (V. renardi)P0DKL822
Kunitz trypsin inhibitor protein 2 (V. renardi)S4S37522
Venom Kunitz trypsin inhibitor protein 2 (V. renardi)S4S37422
Uncharacterized protein (O. hannah)—Peptidase inhibitor activityV8NAQ022
NeprilysinNeprilysin (B. irregularis)1A0A0B8RU8352
NGFVenom nerve growth factor (W. aegyptia)2V9I1F922
Uncharacterized protein (O. hannah)V8N4D822
Ohanin/VesprynOhanin (O. hannah)1P8323455
PDEEctonucleotide pyrophosphatase/phosphodiesterase 3 (M. fulvius)3U3FAB397
Phosphodiesterase (M. fulvius)A0A0F7YYZ897
Phosphodiesterase (M. tener)A0A194ARD797
PDGF/VEGFUncharacterized protein (O. hannah)1V8P6V355
PLA2Acidic phospholipase A2 2 (O. hannah)2Q9DF3344
Basic phospholipase A2 (C. rhodostoma)Q9PVF422
SVGFHepatocyte growth factor activator (O. hannah)1V8NIY222
SVMPMetalloproteinase (type III) 2 (Hypsiglena sp.)21A0A098M15632
Metalloproteinase type III 14a (Hypsiglena sp.)A0A098M22132
Metalloproteinase type III 14b (Hypsiglena sp.)A0A098M21932
Metalloproteinase (type III) 14c (Hypsiglena sp.)A0A098M22332
Metalloproteinase type III 2 (M. tener)A0A194AS4782
Scutellatease-1 (O. scutellatus)B5KFV543
Snake venom metalloproteinase BaP1 (B. asper)P8351232
Snake venom metalloproteinase BjussuMP-2 (B. jararacussu)Q7T1T432
Snake venom metalloproteinase BmooMPalpha-I (B. moojeni)P8531432
Snake venom metalloproteinase bothrojaractivase (B. jararaca)P0C7A932
Snake venom metalloproteinase kistomin (C. rhodostoma)P0CB141414
Snake venom metalloproteinase leucurolysin-A (B. leucurus)P8490732
Zinc metalloproteinase/disintegrin (B. asper)Q072L532
Zinc metalloproteinase/disintegrin (B. insularis)Q5XUW832
Zinc metalloproteinase/disintegrin (C. rhodostoma)P3040355
Zinc metalloproteinase-disintegrin-like atrase-A (N. atra)D5LMJ342
Zinc metalloproteinase-disintegrin-like atrase-B (N. atra)D6PXE852
Zinc metalloproteinase-disintegrin-like BjussuMP-1 (B. jararacussu)Q1PHZ432
Zinc metalloproteinase-disintegrin-like kaouthiagin (N. atra)D3TTC184
Zinc metalloproteinase-disintegrin-like MTP9 (D. coronoides)F8RKV962
Zinc metalloproteinase-disintegrin-like ohanin (O. hannah)A3R0T955
SVSPAlpha and beta fibrinogenase OhS1 (O. hannah)1A8QL5622
TOTAL PROTEINS76
Abbreviations: 3FTx: Three-finger toxin; 5′-NTD: 5′-nucleotidase; AChE: Acetylcholinesterase; CC3H: Complement C3 homolog; CRiSP: Cysteine-rich secretory protein; CTL: C-type lectin; FMO: Flavin monoamine oxidase; NGF: Nerve growth factor; PDE: Phosphodiesterase; PDGF/VEGF: Platelet-derived growth factor/vascular endothelial growth factor; PLA2: Phospholipase A2; SVGF: Snake venom growth factor; SVMP: Snake venom metalloproteinase; SVSP: Snake venom serine protease.
Table
Table 3. List of proteins identified from Malaysian C. rhodostoma venom using the GF-LC-MS/MS technique. Refer to Supplementary File (Table S3) for complete proteomic data and abbreviations lists.
Table 3. List of proteins identified from Malaysian C. rhodostoma venom using the GF-LC-MS/MS technique. Refer to Supplementary File (Table S3) for complete proteomic data and abbreviations lists.
Protein FamilyProtein DescriptionNo. of ProteinsProtein Accession No.No. of PeptidesNo. of Unique Peptides
5′-NTDSnake venom 5′-nucleotidase (B. irregularis)4A0A0B8RXZ9123
5′-nucleotidase (P. flavoviridis)T2HRS944
5′-nucleotidase (P. elegans)A0A077L6M554
5′-nucleotidase (P. flavoviridis)A0A077L7M944
AminopeptidaseAminopeptidase (C. horridus) 2T1DNX844
Aminopeptidase (C. horridus) A0A0K8S2L444
CRiSPCysteine-rich secretory protein LCCL … 2 (C. adamanteus)1A0A0F7ZEB699
CTL/snaclecSnaclec alboaggregin-D subunit alpha (T. albolabris)8P0DM3843
Snaclec rhodocetin subunit alpha (C. rhodostoma)P813975857
Snaclec rhodocetin subunit beta (C. rhodostoma)P813985453
Snaclec rhodocetin subunit gamma (C. rhodostoma)D2YW392523
Snaclec rhodocetin subunit delta (C. rhodostoma)D2YW402222
Snaclec rhodocytin subunit alpha (C. rhodostoma)Q9I8411211
Snaclec rhodocytin subunit beta (C. rhodostoma)Q9I840108
C-type lectin 12a (A. contortrix contortrix)A0A1W7RJZ777
EndonucleaseEndonuclease domain-containing protein (O. hannah)3U3FCT922
Endonuclease domain-containing protein (O. hannah)V8NCE322
Deoxyribonuclease-2-alpha-like (C. adamanteus)—DNaseJ3RZ1444
FMOAmine oxidase (O. okinavensis)8T2HQ57404
Amine oxidase (P. regius)A0A098LWS4123
Amine oxidase (P. guttatus)A0A098LX00112
L-amino-acid oxidase (C. rhodostoma)P8138214799
L-amino acid oxidase (B. schlegelii)A0A024BTN94610
L-amino-acid oxidase (T. stejnegeri)Q6WP39192
L-amino-acid oxidase (L. muta)J7H670152
L-amino-acid oxidase (B. fasciatus)A8QL52102
NGFVenom nerve growth factor (B. jararacussu)4Q90W3832
Nerve growth factor (P. flavoviridis)B1Q3K255
Nerve growth factor (P. flavoviridis)A0A077L85455
BATXNGF1 (B. atrox)A0A1L8D60832
PDECadam 10_PDE-1 (C. adamanteus)18A0A1W7RBB0122
Ectonucleotide pyrophosphatase/phosphodiesterase 3 (M. fulvius)U3FAB352
Ectonucleotide pyrophosphatase/phosphodiesterase 3 (C. horridus)T1DJT5123
Phosphodiesterase (P. flavoviridis)T2HQA0172
Phosphodiesterase (P. flavoviridis)T2HRT4172
Phosphodiesterase (P. flavoviridis)T2HPD6172
Phosphodiesterase (P. flavoviridis)T2HP62172
Phosphodiesterase (M. fulvius)A0A0F7YYZ852
Phosphodiesterase (M. tener)A0A194ARD752
Phosphodiesterase (M. fulvius)A0A0F7YYZ832
Phosphodiesterase (C. adamanteus)A0A0F7Z2Q3123
Phosphodiesterase (C. horridus)T1D6P7123
Phosphodiesterase (S. miliarius barbouri)A0A194AS02123
Phosphodiesterase (M. tener)A0A194ARD732
Phosphodiesterase (M. rudis)A0A141DWM143
Phosphodiesterase (M. lebetina)W8E7D142
Venom phosphodiesterase 1 (C. adamanteus) J3SEZ3122
Venom phosphodiesterase 2 (C. adamanteus)J3SBP3122
PLA2Acidic phospholipase A2 (C. rhodostoma)10Q9PVF043
Acidic phospholipase A2 (C. rhodostoma)Q9PVF1164
Acidic phospholipase A2 (C. rhodostoma)Q9PVF25754
Acidic phospholipase A2 (C. rhodostoma)Q9PVE962
Acidic phospholipase A2 (O. hannah)Q9DF3355
Acidic phospholipase A2 (T. stejnegeri)Q6H3C73432
Basic phospholipase A2 W6D49 (C. rhodostoma)Q9PVF42511
Phospholipase A2 (C. rhodostoma)A0A0H3U2664119
Phospholipase A2 (C. atrox)A0A1J0R08152
Phospholipase A2 1c (C. horridus)A0A0K8RZ1752
PLBPhospholipase B (C. adamanteus)6F8S101165
Phospholipase B (C. adamanteus)A0A0F7Z632165
Phospholipase B (P. regius)A0A098LY7482
Phospholipase B (P. guttatus)A0A098LWY942
Phospholipase B (P. flavoviridis)T2HP6875
Phospholipase B (P. elegans)A0A077L7E775
SVMPAAV1 protein (D. acutus) 23A0A0M4MEY542
BATXSVMPIII (B. atrox)A0A1L8D5X932
Disintegrin triflavin (P. flavoviridis)P2185986
Metalloprotease, mRNA (G. intermedius)A0A0C4ZNF182
Metalloproteinase type III 2 (M. tener)A0A194AS4732
Metalloproteinase type III 2 (M. tener)A0A194AR9173
Metalloproteinase type III 2 (Hypsiglena sp.)A0A098M15622
Metalloproteinase type III 14c (Hypsiglena sp.)A0A098M22322
Metalloproteinase type III 5 (Hypsiglena sp.)A0A098M21522
p-III snake venom metalloprotease (M. ikaheca)A0A024AXX752
Snake venom metalloproteinase kistomin (C. rhodostoma)P0CB14198197
SVMP-disintegrin-like mocarhagin (N. mossambica)Q1074984
SVMP-Hop-15 (H. bungaroides)R4FJZ452
Zinc metalloproteinase/disintegrin (C. rhodostoma)P304035858
Zinc metalloproteinase-disintegrin-like agkihagin (D. acutus)Q1PS4582
Zinc metalloproteinase-disintegrin-like atragin (N. atra)D3TTC242
Zinc metalloproteinase-disintegrin-like atrase-A (N. atra)D5LMJ342
Zinc metalloproteinase-disintegrin-like BfMP (B. fasciatus)A8QL4877
Zinc metalloproteinase-disintegrin-like BmMP (B. multicinctus)A8QL4977
Zinc metalloproteinase-disintegrin-like cobrin (N. kaouthia)Q9PVK752
Zinc metalloproteinase-disintegrin-like halysase (G. halys)Q8AWI582
Zinc metalloproteinase-disintegrin-like HV1 (P. flavoviridis)Q90ZI392
Zinc metalloproteinase-disintegrin-like ohanin (O. hannah)A3R0T91212
SVSPSnake venom serine protease catroxase (C. atrox)12Q8QHK232
Snake venom serine proteinase (C. adamanteus)F8S11632
Serine protease (P. flavoviridis)A0A077L6P3122
Serine proteinase 3 (A. piscivorus)A0A194APA832
Serine proteinase 3a (A. contortrix contortrix)A0A1W7RJX132
Serine proteinase 4 (A. contortrix contortrix)A0A1W7RJU232
Serine proteinase 9 (A. contortrix contortrix)A0A1W7RJU032
Serine proteinase 10 (C. horridus)T1DH1052
Serine proteinase 12a (A. piscivorus)A0A194APD332
Thrombin-like enzyme ancrod (C. rhodostoma)P263246459
Thrombin-like enzyme ancrod-2 (C. rhodostoma)P477973937
Thrombin-like protein (A. piscivorus leucostoma)E9NX1432
TOTAL PROTEINS99
Abbreviations: 5′-NTD: 5′-nucleotidase; CRiSP: Cysteine-rich secretory protein; CTL: C-type lectin; FMO: Flavin monoamine oxidase; NGF: Nerve growth factor; PDE: Phosphodiesterase; PLA2: Phospholipase A2; PLB: Phospholipase B; SVMP: Snake venom metalloproteinase; SVSP: Snake venom serine protease.
Table
Table 4. List of proteins identified from Malaysian O. hannah venom using the GF-LC-MS/MS technique. Refer to Supplementary File (Table S4) for complete proteomic data and abbreviations lists.
Table 4. List of proteins identified from Malaysian O. hannah venom using the GF-LC-MS/MS technique. Refer to Supplementary File (Table S4) for complete proteomic data and abbreviations lists.
Protein FamilyProtein DescriptionNo. of ProteinsProtein Accession No.No. of PeptidesNo. of Unique Peptides
3FTxAlpha-cobratoxin (N. kaouthia)51P0139142
Alpha-elapitoxin-Dpp2d (D. polylepis polylepis)C0HJD732
Alpha-elapitoxin-Dv2a (D. viridis)P0139532
Alpha-elapitoxin-Nno2a (N. oxiana)P0138254
Alpha-elapitoxin-Oh2b (O. hannah)P82662154
Haditoxin (O. hannah)A8N2861616
Kappa-6-bungarotoxin (B. multicinctus)Q9W72932
Long neurotoxin 1 (A. superbus)A8S6A832
Long neurotoxin 1 (N. annulata annulata)P3407492
Long neurotoxin 1 (N. haje haje)P2567494
Long neurotoxin 1 (N. nivea)P0139074
Long neurotoxin 1 (O. hannah)P01387266
Long neurotoxin 2 (O. hannah)P0138655
Long neurotoxin 2 (N. melanoleuca)P0138883
Long neurotoxin 3 (O. hannah)P0752652
Long neurotoxin 4 (O. hannah)P80156143
Long neurotoxin LlLong (L. laticaudata)Q7T2I332
Long neurotoxin OH-17 (O. hannah)Q53B5482
Long neurotoxin OH-37 (O. hannah)Q53B59209
Long neurotoxin OH-55 (O. hannah)Q53B584421
Long neurotoxin OH-56 (O. hannah)Q53B571815
Long neurotoxin OH-57 (O. hannah)Q53B5675
Long neurotoxin LNTX1 (O. hannah)Q2VBP8135
Long neurotoxin LNTX-2 (O. hannah)A8N285143
Long neurotoxin LNTX8 (O. hannah)Q2VBP6213
Long neurotoxin LNTX22 (O. hannah)Q2VBP553
Non-conventional three finger toxin isoform (B. flaviceps)D5J9P942
Oxiana weak toxin (N. oxiana)P8552033
Short neurotoxin OH-35 (O. hannah)Q53B4955
Three-finger toxin (B. multicinctus)E2IU1222
Three-finger toxin (B. multicinctus)E2IU1422
Weak neurotoxin (N. N.)P2918122
Weak neurotoxin OH-72 (O. hannah)Q53B611310
Weak neurotoxin WNTX33 (O. hannah)Q2VBN36555
Weak neurotoxin WNTX34 (O. hannah)Q2VBN21310
Weak neurotoxin 7 (N. N.)P2918143
Weak neurotoxin 10 (N. sputatrix)Q802B222
Weak toxin CM-11 (N. haje haje)P0140133
Weak toxin DE-1 (O. hannah)P01412117110
Weak toxin DE-1 homolog 1 (O. hannah)Q69CJ81212
Weak toxin S4C11 (N. melanoleuca)P0140033
Putative long neurotoxin (A. labialis)B2BRQ832
Putative long neurotoxin (A. labialis)B2BRQ932
Putative long neurotoxin (A. labialis)B2BRR032
Putative long neurotoxin (A. labialis)B2BRR132
Putative long neurotoxin (A. labialis)B2BRR232
Putative long neurotoxin (A. labialis)B2BRR332
Putative long neurotoxin (A. labialis)B2BRR632
Putative long neurotoxin (A. labialis)B2BRR832
Uncharacterized protein (O. hannah)—Toxin activityV8N38975
Uncharacterized protein (O. hannah)—Toxin activityV8N3G955
AChECarboxylic ester hydrolase (S. fasciata)1R4FKE622
CC3HComplement-depleting factor (O. hannah) 7A8QL5588
O. venom factor (O. hannah)I2C0903423
Venom factor (C. horridus)T1E3W822
Venom factor (Sistrurus catenatus tergeminus)A0A194APJ422
Venom factor (Agkistrodon piscivorus)A0A194ARG122
Venom factor (C. adamanteus)A0A0F7Z1I722
Venom factor (C. adamanteus)J3S83622
CRiSPCysteine-rich venom protein (P. olfersii)13Q09GJ932
Cysteine-rich venom protein mossambin (N. mossambica)P0DL1692
Cysteine-rich venom protein ophanin (O. hannah)Q7ZT984414
Cysteine-rich venom protein TRI1 (T. biscutatus)Q2XXP442
Cysteine-rich secretory protein (O. okinavensis)T2HPR872
Cysteine-rich secretory protein (P. chamissonis)A0A0B4SXI822
Cysteine-rich secretory protein Az-CRP (A. feae)F2Q6E362
Cysteine-rich secretory protein Bc-CRP (B. schlegelii)F2Q6E472
Cysteine-rich secretory protein Bc-CRPa (B. candidus)F2Q6G333
Cysteine-rich secretory protein Bc-CRPb (B. candidus)F2Q6G222
Cysteine-rich secretory protein 1b (B. irregularis)A0A0B8RZW732
Cysteine-rich secretory protein 1c (B. irregularis)A0A0B8RYV832
Uncharacterized protein (O. hannah) V8N8B4355
CTL/snaclecSnaclec rhodocetin subunit alpha (C. rhodostoma)1P8139722
EndonucleaseEndonuclease domain-containing protein (O. hannah)2V8NCE332
Endonuclease domain-containing protein (O. hannah)V8N4Y232
FMOAmine oxidase (O. hannah)2V8N3Q9528
L-amino-acid oxidase (O. hannah)P813835915
InsulinInsulin-like growth factor 1 (P. guttatus)5I6TIH61915
Insulin-like growth factor I (O. hannah)V8NHJ875
Insulin-like growth factor 1 (C. helena)I6S9K472
Insulin-like growth factor 1 (T. elegans)I6T66699
Uncharacterized protein (O. hannah)—Hormone activityV8N2J244
Kunitz-typeKunitz-type serine protease inhibitor (O. hannah)6P8296673
Kunitz-type serine protease inhibitor TCI (O. hannah)B6RLX23127
Kunitz-type serine protease inhibitor Vur-KIn (V. renardi)P0DKL843
Kunitz trypsin inhibitor protein 2 (V. renardi)S4S37543
Venom Kunitz trypsin inhibitor protein 2 (V. renardi)S4S37443
Uncharacterized protein (O. hannah)—Peptidase inhibitor activityV8NAQ044
NGFUncharacterized protein (O. hannah) 1V8N4D822
Ohanin/VesprynOhanin (O. hannah)2P832342217
Thaicobrin (N. kaouthia)P8288522
PDECadam 10_PDE-1 (C. adamanteus)11A0A1W7RBB032
Ectonucleotide pyrophosphatase/phosphodiesterase 3 (M. fulvius)U3FAB31111
Ectonucleotide pyrophosphatase/phosphodiesterase 3 (C. horridus)T1DJT532
Phosphodiesterase (M. tener)A0A194ARD71111
Phosphodiesterase (M. fulvius)A0A0F7YYZ81111
Phosphodiesterase (C. horridus)T1D6P732
Phosphodiesterase (C. adamanteus)A0A0F7Z2Q332
Phosphodiesterase (O. okinavensis)U3TBJ532
Phosphodiesterase (O. okinavensis)U3TAI732
Venom phosphodiesterase 1 (C. adamanteus)J3SEZ332
Venom phosphodiesterase 2 (C. adamanteus)J3SBP332
PDGF/VEGFUncharacterized protein (O. hannah) 1V8P6V377
PLA2Acidic phospholipase A2 2 (O. hannah)8Q9DF335449
Acidic phospholipase A2 (O. hannah)Q9DF562017
Acidic phospholipase A2 KBf-grIB (B. fasciatus)P0C55122
Phospholipase A2 12 (M. fulvius)U3FVG152
Phospholipase A2 23 (M. fulvius)U3FVF552
PLA2 (IB)-Tri1 (T. biscutatus)A7X41853
Secretory phospholipase A2 (Group 10) (O. hannah)V8PDP022
Uncharacterized protein (O. hannah) V8N1D822
PLBPhospholipase B-like 1 (O. hannah)1V8ND6899
SVGFHepatocyte growth factor activator (O. hannah)3V8NIY222
Macrophage colony-stimulating factor 1 (O. hannah)V8NW3533
Tumor necrosis factor receptor (O. hannah)V8P0T555
SVMPAsrin (A. superbus)33A6XJS7116
BATXSVMPII2 (B. atrox)A0A1L8D6A942
BATXSVMPII1 (B. atrox)A0A1L8D60042
BATXSVMPII3 (B. atrox)A0A1L8D5Z642
BATXSVMPIII9 (B. atrox)A0A1L8D64142
Metalloproteinase (A. andersonii)Q9PT4762
Metalloproteinase (D. vestigiata)B5G6F762
Metalloprotease BOJUMET III (B. jararacussu)Q7T1T342
Metalloproteinase PIII (T. gracilis)V5IWF422
Metalloproteinase type III 2 (Hypsiglena sp.) A0A098M15662
Metalloproteinase type III 14b (Hypsiglena sp.) A0A098M21974
Metalloproteinase type III 14c (Hypsiglena sp.) A0A098M22352
Metalloproteinase type III 1 (M. fulvius) U3EPC755
Metalloproteinase type III 2 (M. tener) A0A194AS47164
Metalloproteinase type III 2b (M. tener) A0A0F7YYV152
Scutellatease-1 (O. scutellatus)B5KFV5137
SVMP-disintegrin-like mocarhagin (N. mossambica)Q1074995
Snake venom metalloproteinase kistomin (C. rhodostoma)P0CB1444
Snake venom metalloprotease (P.olfersii)C6JUN342
Snake venom metalloprotease (P. olfersii)C6JUN442
Stephensease-1 (H. stephensii)B5KFV482
SVMP-Aca-4 (A. wellsi)R4G2D3114
SVMP-Den-9 (D. devisi)R4FIC4152
SVMP-Ech-32 (E. curta)R4FJM662
SVMP-Fur-1 (F. ornata)R4G2G342
SVMP-Hem-2 (H. signata)R4G2W983
SVMP-Hop-14 (H. bungaroides)R4FIM153
Zinc metalloproteinase-disintegrin-like atrase-A (N. atra)D5LMJ3143
Zinc metalloproteinase-disintegrin-like atrase-B (N. atra)D6PXE893
Zinc metalloproteinase-disintegrin-like cobrin (N. kaouthia)Q9PVK772
Zinc metalloproteinase-disintegrin-like kaouthiagin-like (N. atra)D3TTC193
Zinc metalloproteinase-disintegrin-like MTP4 (D. coronoides)F8RKW1113
Zinc metalloproteinase-disintegrin-like ohanin (O. hannah)A3R0T91910
SVSPAlpha and beta fibrinogenase OhS1 (O. hannah)2A8QL5664
Neuroendocrine convertase 1 (O. hannah)V8P1Y233
TOTAL PROTEINS150
Abbreviations: 3FTx: Three-finger toxin; AChE: Acetylcholinesterase; CC3H: Complement C3 homolog; CRiSP: Cysteine-rich secretory protein; CTL: C-type lectin; FMO: Flavin monoamine oxidase; NGF: Nerve growth factor; PDE: Phosphodiesterase; PDGF/VEGF: Platelet-derived growth factor/vascular endothelial growth factor; PLA2: Phospholipase A2; PLB: Phospholipase B; SVGF: Snake venom growth factor; SVMP: Snake venom metalloproteinase; SVSP: Snake venom serine protease.
Table
Table 5. Relative abundance of snake venom toxin families identified in Malaysian C. rhodostoma and O. hannah venoms. The relative abundance is presented in percentages (%) alongside the number of protein subtypes included in parentheses. Major species-specific families are highlighted in bold (CR, blue; OH, red).
Table 5. Relative abundance of snake venom toxin families identified in Malaysian C. rhodostoma and O. hannah venoms. The relative abundance is presented in percentages (%) alongside the number of protein subtypes included in parentheses. Major species-specific families are highlighted in bold (CR, blue; OH, red).
Snake Venom Protein FamiliesCalloselasma rhodostomaOphiophagus hannah
3FTxND31 (55)
Complement C3 HomologND4 (7)
Kunitz-typeND3 (6)
InsulinND3 (5)
SVGFND2 (3)
AChEND1 (2)
Ohanin/VesprynND1 (2)
NeprilysinND1 (1)
PDGF/VEGFND1 (1)
SVMP20 (23)25 (44)
PDE16 (18)6 (11)
SVSP12 (14)1 (2)
CTL/snaclec12 (14)1 (2)
PLA210 (11)5 (9)
FMO7 (8)2 (3)
PLB6 (7)1 (1)
5′-NTD5 (6)2 (3)
NGF4 (4)1 (2)
Endonuclease3 (3)1 (2)
CRiSP1 (1)9 (15)
Aminopeptidase2 (2)ND
QPCT1 (1)ND
Trypsinogen1 (1)ND
Others (ANK)1 (1)ND
TOTAL100 (114)100 (176)
Abbreviations: 3FTx: Three-finger toxin; 5′-NTD: 5′-nucleotidase; AChE: Acetylcholinesterase; ANK: Ankyrin repeat proteins; CTL: C-type lectin; CRiSP: Cysteine-rich secretory protein; FMO: Flavin monoamine oxidase; NGF: Nerve growth factor; PDE: Phosphodiesterase; PDGF/VEGF: Platelet-derived growth factor/vascular endothelial growth factor; PLA2: Phospholipase A2; PLB: Phospholipase B; QPCT: Glutaminyl-peptide cyclotransferase; SVGF: Snake venom growth factor; SVMP: Snake venom metalloproteinase; SVSP: Snake venom serine protease; ND: not determined.

Share and Cite

MDPI and ACS Style

Kunalan, S.; Othman, I.; Syed Hassan, S.; Hodgson, W.C. Proteomic Characterization of Two Medically Important Malaysian Snake Venoms, Calloselasma rhodostoma (Malayan Pit Viper) and Ophiophagus hannah (King Cobra). Toxins 2018, 10, 434. https://doi.org/10.3390/toxins10110434

AMA Style

Kunalan S, Othman I, Syed Hassan S, Hodgson WC. Proteomic Characterization of Two Medically Important Malaysian Snake Venoms, Calloselasma rhodostoma (Malayan Pit Viper) and Ophiophagus hannah (King Cobra). Toxins. 2018; 10(11):434. https://doi.org/10.3390/toxins10110434

Chicago/Turabian Style

Kunalan, Sugita, Iekhsan Othman, Sharifah Syed Hassan, and Wayne C. Hodgson. 2018. "Proteomic Characterization of Two Medically Important Malaysian Snake Venoms, Calloselasma rhodostoma (Malayan Pit Viper) and Ophiophagus hannah (King Cobra)" Toxins 10, no. 11: 434. https://doi.org/10.3390/toxins10110434

APA Style

Kunalan, S., Othman, I., Syed Hassan, S., & Hodgson, W. C. (2018). Proteomic Characterization of Two Medically Important Malaysian Snake Venoms, Calloselasma rhodostoma (Malayan Pit Viper) and Ophiophagus hannah (King Cobra). Toxins, 10(11), 434. https://doi.org/10.3390/toxins10110434

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