Next Article in Journal
Variability of Botulinum Toxins: Challenges and Opportunities for the Future
Previous Article in Journal
A New Concept to Secure Food Safety Standards against Fusarium Species and Aspergillus Flavus and Their Toxins in Maize
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Proteomic Deep Mining the Venom of the Red-Headed Krait, Bungarus flaviceps

1
Fundação Oswaldo Cruz-Ceará, Rua São José, 2º Pavimento, Precabura, Eusébio 61760-000, Brazil
2
Laboratório de Toxinologia, Instituto Oswaldo Cruz, Fiocruz, Rio de Janeiro 21045-900, Brazil
3
Computational Mass Spectrometry& Proteomics Group, Carlos Chagas Institute, Fiocruz, Paraná 81350-010, Brazil
4
Department of Biology, Stetson University, 421 N. Woodland Blvd, DeLand, FL 32723, USA
5
Department of Biological Sciences, National University of Singapore, 14 Science Drive 4, Singapore 117543, Singapore
6
School of Biological Sciences, University of Northern Colorado, 501 20th St., CB 92, Greeley, CO 80639-0017, USA
*
Authors to whom correspondence should be addressed.
Submission received: 13 June 2018 / Revised: 28 August 2018 / Accepted: 1 September 2018 / Published: 13 September 2018
(This article belongs to the Section Animal Venoms)

Abstract

:
The use of -omics technologies allows for the characterization of snake venom composition at a fast rate and at high levels of detail. In the present study, we investigated the protein content of Red-headed Krait (Bungarus flaviceps) venom. This analysis revealed a high diversity of snake venom protein families, as evidenced by high-throughput mass spectrometric analysis. We found all six venom protein families previously reported in a transcriptome study of the venom gland of B. flaviceps, including phospholipases A2 (PLA2s), Kunitz-type serine proteinase inhibitors (KSPIs), three-finger toxins (3FTxs), cysteine-rich secretory proteins (CRISPs), snaclecs, and natriuretic peptides. A combined approach of automated database searches and de novo sequencing of tandem mass spectra, followed by sequence similarity searches, revealed the presence of 12 additional toxin families. De novo sequencing alone was able to identify 58 additional peptides, and this approach contributed significantly to the comprehensive description of the venom. Abundant protein families comprise 3FTxs (22.3%), KSPIs (19%), acetylcholinesterases (12.6%), PLA2s (11.9%), venom endothelial growth factors (VEGFs, 8.4%), nucleotidases (4.3%), and C-type lectin-like proteins (snaclecs, 3.3%); an additional 11 toxin families are present at significantly lower concentrations, including complement depleting factors, a family not previously detected in Bungarus venoms. The utility of a multifaceted approach toward unraveling the proteome of snake venoms, employed here, allowed detection of even minor venom components. This more in-depth knowledge of the composition of B. flaviceps venom facilitates a better understanding of snake venom molecular evolution, in turn contributing to more effective treatment of krait bites.
Key Contribution: We describe the deep proteomic mining of Bungarus flaviceps venom that allowed us to identify the presence of 18 protein families, including several low abundance toxins. Utilization of this approach facilitates the understanding of the evolution of venom complexity and the differential expression of dominant venom toxins by allowing detection/identification of even trace venom components.

Graphical Abstract

1. Introduction

In the last decade, there has been a tremendous increase in the knowledge of snake venom composition and evolution, mainly because of the application of “omics” techniques, in particular, high-throughput transcriptomic investigations of venom gland tissue in combination with proteomic studies of venom [1,2,3,4,5,6]. The use of these highly sensitive technologies now makes it feasible to study the venom gland transcriptomes/proteomes of virtually any snake species in detail, even those producing low venom yields, such as kraits and many rear-fanged colubrid snakes. Recently, Viala et al. investigated the venom of the Australian Eastern Brown Snake (Pseudonaja textilis) by combining high throughput proteomics and transcriptomics techniques [7]. Their deep data-mining approach led to the identification of a new transcript variant of venom coagulation factor 5a. The authors also concluded that splicing variants represent an important source of toxin diversification. Likewise, a comparative deep venomics study including venom gland transcriptomes and proteomes of the Australian western brown snakes Pseudonaja aspidorhyncha, Pseudonaja nuchalis, and Pseudonaja textilis revealed not only accelerated evolution of different toxin families, but also a wide array of novel peptide sequences from six bioactive peptide families [8]. In addition, an integrated analysis of the transcriptome and proteome of the Amazon Puffing Snake Spilotes sulphureus identified two novel three-finger toxins: sulditoxin and sulmotoxin-1 [9]. These two toxins not only represent the most abundant venom proteins, but also showed distinct prey-specific toxicities. While sulditoxin exhibits high toxicity towards lizard prey, but is non-toxic towards mammalian prey, sulmotoxin 1 shows the reverse trend [9]. The authors reasoned that this is the first example in which a contrasting pattern of toxicity in relation to diet has been shown for proteins in the venom of a single snake. As these technologies are applied to poorly known venoms, new toxins may be discovered, as demonstrated by an investigation of the proteome and transcriptome of the venom of the Australian elapid snake Drysdalia coronoides, which revealed, among other protein families, novel three-finger toxins (3FTx) proteins [10]. These unique sequences could indicate novel biological functions related to specific types of prey consumed. For example, Brown Treesnake (Boiga irregularis), and Green Vinesnake (Oxybelis fulgidus), two rear-fanged colubrid snakes, both feed primarily on lizards and birds, and their venoms contain different novel 3FTxs that are specific for these taxa and are harmless toward mammals [11,12], strongly indicating selection for taxon-specific venom toxins. These examples point to the significance of deep mining venomics studies by, for example, combined transcriptome and proteome investigations to detect even less abundant venom sequences that might be of biological or pharmaceutical importance. In the present study, we investigated the venom proteome of B. flaviceps by a complementary approach of electrospray and MALDI mass spectrometry (MS) to probe deeply for even minor venom compounds.
The Bungarus genus (kraits) belongs to the family Elapidae. To date, 15 different species of this genus have been recognized. The genus Bungarus was notorious during the Vietnam War, and American troops referred to these snakes as the “two-step snake” in the mistaken belief that after being bitten, a person was only two steps away from death [13]. Investigations of venom gland transcriptomes [14] and proteomes [15] of B. multicinctus revealed that tandem duplications contributed significantly to the expansion of toxin multigene families in the transcriptome. Interestingly, the venom proteome [15] appeared to be less diverse than the corresponding transcriptome of B. multicinctus [14]. In addition, recent proteomic studies of the Malaysian kraits B. fasciatus and B. candidus led to the identification of natriuretic peptides, vespryns, and serine proteases, protein families previously not detected in these species. The venoms of the Vietnamese kraits B. multicinctus and B. fasciatus showed a pronounced difference in the amounts of three-finger toxins, which were highly abundant in the former venom (28%), but, surprisingly, were detected in very low amounts in the latter (1%). A recent study of venom variation in Bungarus caeruleus (Common Krait) concluded that Pakistani B. caeruleus should be incorporated into the Indian Vins Polyvalent Antivenom (VPAV), because krait bites at this locality showed only moderate neutralization by VPAV [16].
Bungarus flaviceps, the Red-headed Krait, is widely distributed across southeast Asia, from southern Myanmar and Thailand south to Borneo and Java [17] (Figure 1).
Initially, the transcriptome of B. flaviceps was investigated by sequencing 845 expressed sequence tags (ESTs) from a single venom gland of a B. flaviceps specimen. Roughly 75% of the transcripts were putative toxins and six different protein venom families were identified. However, the EST method has rather low resolution, especially when compared with next-generation high-throughput DNA sequencing techniques [8], which provide a deeper understanding of the venom composition.
Bungarus flaviceps feeds primarily on other snakes and on lizards [17], so its venom may also show novel toxins or specific compositional/structural motifs targeting these preferred prey. As B. flaviceps does not feed on mammals, but only on ectothermic reptiles, we might anticipate that specialized toxins may have evolved that more effectively incapacitate ectothermic prey. Venoms of several species of the genus Bungarus have been studied, and a general characteristic is a high content of potent neurotoxins, among which α-, β-, and κ-bungarotoxins are the most well-known [18,19,20]. The mass spectrometric data of the present study were analyzed by automated searches against a custom snake venom database, as well as by de novo sequencing of the tandem mass spectra followed by sequence-similarity analysis. These investigations revealed the presence of 18 snake venom protein families as evidenced by the detection of corresponding peptides, including protein families not previously detected in Bungarus venom. The results of this study have direct implications for envenomation treatment and for molecular venom evolution, particularly among elapid species.

2. Results and Discussion

The venom proteome of the elapid B. flaviceps was investigated by mass spectrometric analysis using complementary ionization techniques, including electrospray and MALDI. Based on the results of database searches and BLAST analyses of de novo sequenced tandem mass spectra, we were able to provide a comprehensive overview of the different protein families present in the venom. In addition, spectral counts of the tandem MS (MS/MS) data matched to known toxin sequences yielded a quantitative estimate of the different protein families in the venom. In total, 88 different proteins were identified that belong to 18 distinctive protein families; a similar number of proteins was found for the venoms of B. candidus and B. fasciatus [21]. Below, we discuss the different toxin families of B. flaviceps venom in detail.

2.1. Major Venom Components

The most abundant components of the venom of B. flaviceps (22.3%; Figure 2 and Figure 3) are represented by 3FTxs, similar to that observed for the related B. fasciatus, in which 3FTxs make up 30% of the venom content [21]. As a non-enzymatic snake venom protein family, 3FTxs are structurally characterized by three-stranded loops protruding from a central core of the molecule. In spite of their similarity in structure and relatively small size (<10 kDa), 3FTxs have remarkable differences in function, including platelet aggregation, neurotoxicity, and cytotoxicity [22,23,24]. Furthermore, 3FTxs were identified as the major component of the B. flaviceps venom at the transcript level [25] and included both non-conventional and short-chain isoforms. We confirmed the presence of most of these 3FTxs in the proteome (Table 1), in addition to a protein similar to the κ-bungarotoxin from B. multicinctus (Table 1). However, no cytotoxins or cardiotoxins were encountered in the venom of B. flaviceps, consistent with the results of a recent analysis of B. caeruleus venom, in which these proteins were also not detected [16]. Interestingly, we were not able to demonstrate the presence of α-bungarotoxin, a 3FTx detected in the venoms of the congeneric species B. candidus, B. fasciatus [21], and B. multicinctus [26]. As it was also not detected in the transcriptome study of B. flaviceps [25], we speculate that its absence could indicate a differential need of suite of toxins for prey taken by B. flaviceps. Alternatively, it might indicate a random loss of the gene during the differentiation of B. flaviceps from a last common ancestor.
Serine proteinase inhibitors were the second most abundant toxin family in the venom of B. flaviceps and accounted for 19% of the total spectral count (Figure 3). Kunitz-type serine proteinase inhibitors (KSPI) carry a conserved Kunitz motif similar to that found in bovine pancreatic trypsin inhibitors [27,28]. They not only inhibit trypsin and chymotrypsin [29,30], but also cause different pharmacological effects, including anticoagulation and fibrinolysis. At the transcript level, six isoforms of KPSI were detected [25], while the proteome analysis revealed the presence of three different isoforms, suggesting that three isoforms are not expressed in the venom. However, the difference in the number of isoforms detected in the complementary studies might be related to poor ionization of unique peptides of the corresponding isoform rather than the absence of this protein in the venom. The B-chain of β-bungarotoxin exhibits high homology with KSPIs and, therefore, has been grouped together with this toxin family. At the protein level, we found four different isoforms of the B-chain (Table 1 and Table 2), but only one of those matched one of the four isoforms discovered in the transcriptome [25], perhaps because of differences in snakes sampled in the two studies or because of the lack of fine scale transcriptomic resolution used by Siang et al.
Acetylcholinesterases (AChEs) are important regulators of neurotransmission at the neuromuscular junction, and they rapidly hydrolyze acetylcholine. Acetylcholinesterases were the third most abundant toxin family in the venom (12.6%; Figure 3). The function of this protein family in snake venoms is still poorly understood, but it may potentiate the action of α-neurotoxins at the nAChR of the motor endplate. BfAChE was identified with relatively high sequence coverage, consistent with the view that only a single isoform of this protein family is present in the venom of B. flaviceps [25]. Of note is that AChEs in related Bungarus venoms are less abundant. For example, Malaysian B. candidus and B. fasciatus venoms revealed relative abundances of 4.9% and 4.7%, respectively [21]. In addition, AChEs were detected in very low abundances (1.1%) in the Vietnamese B. multicinctus venom [26]. This might indicate that AChEs play a more important role in envenomation in B. flaviceps than in congeneric venoms.
Phospholipases A2 (PLA2s) are broadly occurring esterolytic enzymes, and in venoms, they cause diverse pharmacological consequences including myotoxic, cardiotoxic, and anticoagulant pathologies, as well as pre- and postsynaptic neurotoxic effects [31]. While PLA2s represent the fourth most abundant component of the venom of B. flaviceps (11.9%; Figure 3), this protein family is considerably more abundant in B. fasciatus and B. candidus venoms, where PLA2s make up 44% and 25%, respectively, of the venom content [21]. The A chain of β-bungarotoxin, a heterodimeric protein with presynaptic neurotoxic effects [20], shows pronounced sequence similarity to PLA2s and we were able to detect two isoforms of the A chain in the venom (Table 1), while in the transcriptome, three isoforms were detected. Given that β-bungarotoxins are the main lethal components in B. multicinctus and B. candidus, the venom of B. flaviceps is most likely potently neurotoxic, with envenomation effects similar to those of other kraits. Additionally, we not only confirmed the presence of PLA2s also observed at the transcript level [25], but also found peptides homologous to those of PLA2s isolated from divergent elapid species, including B. fasciatus, B. candidus, Austrelaps superbus, Micropechis ikaheka, Laticauda colubrina, and L. semifasciata. This indicates that PLA2s are more diversified in B. flaviceps venom than was initially appreciated, and they are likely functionally important venom components.
Vascular endothelial growth factor (VEGF) is an important angiogenic factor in mammals and is produced by various cell types, including macrophages [32], platelets [33], and tumor cells [34]. VEGFs have also been identified in snake venoms (known as VEGF-F) and they have been found to strongly stimulate proliferation of vascular endothelial cells in vitro, as well as to induce hypotension in rats [35]. Four peptides that matched the sequence of a hypothetical venom protein of Ophiophagus hannah (King Cobra) were identified in the present study. As this hypothetical protein also shares about 80% sequence homology with a VEGF-F from the viper Vipera ammodytes [36], it is highly probable that VEGF-Fs are also part of the venom proteome of B. flaviceps.
Nucleotidases in snake venoms are believed to catalyze the production of adenosine, which could lead to paralysis and prey immobilization [37], but their role in envenomation is still unclear. Tryptic peptides derived from B. flaviceps venom matched two isoforms of nucleotidases originally isolated from venoms of the elapid Micrurus fulvius [38] and the colubrid Boiga irregularis [39]. Nucleotidases are typically found in low abundance in snake venoms, and it would be of interest to investigate their biological functions.
Snaclecs are mostly heterodimeric proteins that contain the same three-dimensional fold as classical C-type lectins [40]. Through binding to the von Willebrand factor or to different collagen receptors, they affect platelet activation. In this study, we detected peptides both identical to and similar to snaclec transcripts from B. fasciatus, Micrurus fulvius, and Bothrops jararaca. This appears to be the first identification of this protein family at the protein level in the genus Bungarus. Indeed, several reports found evidence of snaclecs in the venom glands of B. flaviceps [25], as well as B. fasciatus and B. multicinctus [41] at the transcript level, but to our knowledge, this is the first demonstration of snaclecs in the Bungarus proteome.
The venom of O. hannah was the original source of a novel snake venom family known as the vespryns [42]. Functionally, these proteins have been described to induce hypolocomotion and hyperalgesia in mice [42]. Vespryns matching those from several species, including O. hannah, Pseudechis australis, and Drysdalia coronoides, were encountered in the present study. Recently, vespryns were also found in the venom of B. candidus [21], but tryptic peptides of B. flaviceps origin did not match these proteins.
Phosphodiesterases (PDE) are basic enzymes with molecular masses in the range of 98 to 140 kDa [43], and they catalyze the hydrolysis of phosphodiester bonds from the 3’ terminus of polynucleotides. Though PDEs have long been known to be present in many venoms, the first complete primary structure of a PDE from snake venom was only recently published [44]. We encountered sequences similar to PDEs of the elapid Micrurus fulvius and of several different pitvipers (Table 1 and Table 2). These results further corroborate the presence of PDE in Bungarus venoms, as they were also reported in the venom of B. fasciatus [21].
Hyaluronidases from snake venoms are endo-β-N-acetyl-hexosaminidases with molecular masses in the range of 30 to 110 kDa [43]. These enzymes cleave hyaluronan, a major glycosaminoglycan constituent of the extracellular matrix, into N-acetylglucosamine and oligocarbohydrates. Therefore, hyaluronidases have been implicated as important factors for the distribution and dissemination of the venom in tissues and are often called “spreading factors” [43]. Peptides with pronounced homology to a hyaluronidase from the venom of the African Puff Adder, Bitis arietans, were identified in the venom of B. flaviceps. The low abundance (2.2%) of B. flaviceps hyaluronidase is consistent with observations in other venoms that this component is common but not abundant.
Nerve growth factors (NGFs) are relatively small proteins (up to 26 kDa) that induce growth and proliferation of certain neurons [45]. The NGFs encountered in snake venoms (vNGF) share significant sequence homology to their mammalian counterparts, but relatively little is known about their function and contribution to the envenomation process [46]. We detected three peptides related to vNGFs with homology to those from B. fasciatus and B. multicinctus, corroborating studies on the venoms of B. candidus [21], B. multicinctus, and B. fasciatus [26], in which vNGFs were also found in very low amounts.
The enzyme LAAO causes oxidative deamination of (primarily hydrophobic) L-amino acids and ultimately leads to the release of ammonia, α-keto acids and, importantly, hydrogen peroxide [47]. This compound (H2O2) is associated with various detrimental biological effects, including apoptosis induction and platelet aggregation [48]. It has also been speculated that the anti-bacterial effects of LAAOs, also likely dependent on peroxide production, may be related to maintenance and stabilization of the venom in the venom gland [43], but this seems unlikely, as several other stabilization mechanisms are known to be present [49].

2.2. Minor Venom Components

In the present study, there were also many toxin families present in very low abundance (up to 1.5%), some of which were detected by only a single peptide (Supplemental Table S2). The functions of some of the cysteine-rich secretory protein (CRISP or helveprin) family are related to the blockage of different ion channels, specifically, the inhibition of cyclic nucleotide gated ion channels and calcium channels, which could lead to the disruption of homeostasis in the prey [50]. While the predicted CRISP protein from the corresponding transcriptome study was not detected [25], we found peptides in the venom related to CRISP proteins from different snake families, including elapids, viperids, and colubrids, which might indicate an enhanced diversification of this toxin in B. flaviceps (Table 1 and Table 2). In addition, peptides related to a phospholipase B (PLB) from M. fulvius [38] were also found (Table 1 and Table 2). PLBs share PLA1 and PLA2 enzymatic activities and have been reported in venoms at the protein level in only a few species, including Sistrurus catenatus edwardsii [51], B. candidus [21], and B. fasciatus [26]. SVMPs are common components in viperid and colubrid venoms, in strong contrast to elapid venoms, where they often in much lower abundance [52]. SVMPs can be subdivided, based on in their domain structure, into different classes that are related to diverse biological activities [53]; P-III SVMPs were detected in the venom of B. flaviceps and accounted for 1.2% of the total spectral count (Figure 3). While most of the identified sequences were related to SVMPs of congeneric B. multicinctus and the elapid Notechis scutatus, we also identified peptides matching to SVMPs of viperid (e.g., Bothrops asper) and colubrid snakes, indicating an extensive diversification of this toxin that might point to accelerated evolution of this particular venom protein family in B. flaviceps [54]. Also, in this context, the sequence similarity search of the tandem mass spectra against a venom database proved to be especially useful, and most of the identified peptides belonged to SVMPs (Supplemental Table S2). Recent investigations of the closely related venoms B. candidus and B. fasciatus found SVMPs to account for approximately 5% of the total venom, while in the venom of B. multicinctus [26], the amount was rather low (0.8%). However, in these other Bungarus venoms, the SVMPs were much less diversified than in the B. flaviceps venom. In viperids, abundant SVMPs induce local and systemic hemorrhagic effects following envenomation, and we speculate that minor effects might occur following B. flaviceps envenomation.
Serine proteases are abundant components of many viperid venoms, and as a result of a detailed proteome analysis, they have recently been described in the venoms of B. fasciatus and B. candidus [21]. However, the first evidence of this protein family in the genus Bungarus stems from a study by Zhang and co-workers that identified a factor X activator in the venom of B. fasciatus as a serine protease [55]. In the present study, three different proteins similar to serine proteases from the pitvipers Crotalus adamanteus, Sistrurus catenatus edwardsii, and Deinagkistrodon acutus were detected (Table 1 and Table 2).
Natriuretic peptides play an important role in pressure-volume homeostasis, affecting blood pressure and renal filtration, and they have been isolated from several snake venoms. For example, the venom of the Green Mamba, Dendroaspis angusticeps, contains a natriuretic peptide (NP) known as Dendroaspis natriuretic peptide (DNP), which causes vasodilation and lowers blood pressure in victims [56]. A recent proteome study of the venoms of B. fasciatus and B. candidus revealed for the first time the presence of a natriuretic peptide (NP) at the protein level in the genus Bungarus [21], and the authors reasoned that it might be responsible for hyponatremia observed in victims of bites from B. candidus [21,57]. A full-length mRNA sequence of an NP was found in the transcriptome of B. flaviceps [25], and a protein called KNP was characterized [58]. We were able to confirm the translation of this protein by the detection of a single semi-tryptic peptide corresponding to the predicted sequence (Table 1). NPs isolated from snake venoms, especially DNP, have greater stability and increased potency compared with mammalian NPs and are considered potential drug lead candidates [59].
Cobra venom factor (CVF) is a glycoprotein found in different cobra venoms that activates the complement system in prey [60,61]. Complement activation by CVF eventually leads to the exhaustive consumption of many compounds of the innate immune system and the weakening of the bite victim. In this study, we confirmed the presence of a peptide identical to a predicted CVF sequence from B. fasciatus.

2.3. Quantification of Venom Components

In the present study, we have quantified the snake venom protein families based on spectral count of the corresponding tandem mass spectra followed by normalization in order to take into consideration the different length of the proteins. While this protocol is sound for those proteins also identified at the transcriptome level, as their exact amino acid sequence is known, it might cause a bias for proteins detected only by sequence identity or homology from other snake venoms, because the length of these proteins in B. flaviceps venom is not known and can only be estimated. For example, several different SVMPs have been identified in the venom proteome, but there is no sequence information on SVMPs from the corresponding transcriptome analysis of the B. flaviceps venom. However, it is reasonable to assume that they belong to the PIII class (approx. 600 amino acids) of SVMPs because in most elapid venoms characterized to date, this type of SVMP is the most abundant one [62]; additionally, we detected low metalloproteinase activity in assays of the venom (data not shown).
It is also interesting to compare the method of venom protein quantification used in the present study to other quantitative snake venom investigations. Quantitative analysis of the venoms of B. candidus and B. fasciatus, for example, was based on the number of the corresponding proteins of the different protein venom families identified [21]. To estimate the abundance of each protein family, the investigators expressed these numbers (number of proteins per family) in percentage compared to the total number of identified proteins. Using yet another method, Pla and coworkers estimated the relative protein abundances by separating the crude venom with reversed-phase HPLC [63]. Individual protein abundances were calculated by dividing the single peak area by the total area of the chromatogram [63]; in cases were more than one protein was present in a chromatographic fraction, their corresponding percentage was estimated by densitometry of Coomassie-stained SDS gels. These few examples illustrate how differently snake venom proteins are quantified by different workers, and it should be stressed that there is currently no consensus in the scientific community of how best to estimate protein venom abundances. This variety of the quantification methods makes it challenging to compare different snake venom protein analyses, and it indicates that official guidelines for the experimental protocols need to be delineated in order for results from different laboratories to generate comparable results. This would lead to more robust and more readily comparable datasets, which in turn will have a profound impact on investigations into the molecular mechanisms of snake venom evolution. We would like to stress that the present study is only the view of the authors to quantify protein families and we do not claim it represents the “best” solution for this task in general.

3. Conclusions

We have described one of the most exhaustive characterizations of a venom proteome of an elapid snake to date, revealing a venom complexity formerly unknown to the genus Bungarus [21,25]. Since the currently available venom databases are largely incomplete, the use of de novo sequencing as a tool to characterize protein families and isoforms has proven to be particularly useful in the present study, as evidenced by the detection of 58 additional peptides. Several protein families (3FTxs, serine protease inhibitors, acetylcholinesterases, and PLA2s) represent the most abundant toxin classes in the venom of B. flaviceps, but we also detected numerous less abundant snake venom proteins, including representatives of PLB, PDE, CRISP, vNGF, hyaluronidase, natriuretic peptide, nucleotidase, VEGF-F, acetylcholinesterase, SVSP, and vespryn venom protein families. In addition, we found protein families not yet reported in krait venoms (snaclecs, and complement depleting factors), indicating that Bungarus venoms are more complex and diverse than previously reported. However, our results represent the venom composition of only three individual specimens from a limited part of the species’ range, and further studies are needed to determine possible geographic, gender or age-related sources of variation in the venom composition of B. flaviceps. The results of this study have important implications for the most effective approaches to use for generating near-complete venom proteomes, and the approaches detailed here may shed light on divergence of snake venom gland genomes and how this has affected the molecular evolution of snake venoms.

4. Materials and Methods

4.1. Venom Extraction and Ethics Statement

Three specimens of Bungarus flaviceps were imported from southwest Sumatra (Bushmaster Reptile; Boulder, CO, USA) and maintained at the University of Northern Colorado Animal Resource Facility in accordance with UNC-IACUC protocol #1302D-SM-S-16 (approval date: 30 March 2016). Venom was extracted manually and then centrifuged at 10,000 rpm for 5 min, frozen at −80 °C, lyophilized, and stored at −20 °C until use; venom used in this study represents the combined extractions of these three individuals: (1) male, 1405 mm snout-vent length, 240 mm tail length, total length = 1645 mm; (2) female, 940 mm snout-vent length, 168 mm tail length, total length = 1108 mm; (3) female, 1060 mm snout-vent length, 150 mm tail length, total length = 1210 mm. Snakes were extracted four times each and the average yield was 24.7 L. Range = 5–75 L. No animals were sacrificed during this study.

4.2. Tryptic Digestion of Crude Venom

Lyophilized venom (100 µg) was dissolved in 40 L of 0.4 M ammonium bicarbonate and 8 M urea. After adding 10 L of 50 mM dithiothreitol (DTT), the solution was incubated for 3 h at 37 °C. Additionally, 10 µL of iodoacetamide (150 mM) was added at room temperature and the reaction was allowed to proceed for 15 min. The reaction was quenched with 50 mM of DTT (6 µL) for 15 min. The last two steps were performed in the dark. For digestion, the sample solution was diluted to 1 M urea (by adding 254 µL of water) and incubated with 10 µL (2 µg in 50 mM acetic acid) of trypsin (Promega, Madison, WI, USA) at 37 °C overnight. Quenching of the reaction occurred by adding trifluoroacetic acid (40 µL), and desalting of the sample was performed on Poros R2 microcolumns. As a final step, the peptides were dried in a vacuum centrifuge and brought up in 1% formic acid solution (approx. 50 µL).

4.3. Chromatography and Mass Spectrometry

4.3.1. One-Dimensional Sodium Dodecyl Sulfate Polyacrylamide Gel Electrophoresis (SDS-PAGE)

Crude venom was reduced with 2.5% β-mercaptoethanol in Laemmli sample buffer (Bio-Rad, Hercules, CA, USA.) by incubation at 95 °C for 10 min. Samples (20 µg and 50 µg amounts) were run on a 12% Bis-Tris acrylamide gel, stained with 0.1% Coomassie Brilliant Blue R-250 overnight, and destained in 10% acetic acid: 40% methanol: 50% ddH20. A Novex Mark 12 unstained mass standard (Life Technologies, Grand Island, NY, USA) was also run for band mass estimation. Protein bands were excised and submitted to the Protein and Proteomics Centre in the Department of Biological Sciences, National University of Singapore, for LC-MS/MS analysis.

4.3.2. Liquid Chromatography-Tandem Mass Spectrometry (LC-MS/MS) Analysis of Excised Gel Bands

Bands were reduced with 10 mM DTT (dithiothreitol) and alkylated with 55 mM IAA (iodoacetamide), then digested with trypsin (13 ng/µL) overnight in 25 mM ammonium bicarbonate, 10% ACN (acetonitrile). Samples were desalted using a Sep-Pak tC18 μ Elution Plate (Waters, Milford, MA, USA), and reconstituted with 20 L of diluent (97.5% H2O, 2% ACN, 0.05% formic acid). Peptide separation was carried out on an Eksigent nanoLC Ultra and ChiPLC-nanoflexLC-MS (Eksigent, Dublin, CA, USA) in Trap Elute configuration. A total of 5 L of the sample was loaded onto a 200 μm × 0.5 mm trap column and eluted on an analytical 75 μm × 150 mm column. Trap and analytical columns were made of ChromXP C18-CL, 3 μm (Eksigent, Dublin, CA, USA). Peptides were separated by a gradient formed by 2% ACN, 0.1% FA (mobile phase A) and 98% ACN, 0.1% FA (mobile phase B); 5 to 7% of mobile phase B in 0.1 min, 7 to 30% of mobile phase B in 10 min, 30 to 60% of mobile phase B in 4 min, 60 to 90% of mobile phase B in 1 min, 90 to 90% of mobile phase B in 5 min, 90 to 5% of mobile phase B in 1 min and 5% of mobile phase B for 10 min, at a flow rate of 300 nL/min. The MS analysis was performed on a TripleTOF 5600 system (AB SCIEX, Redwood City, CA, USA) in Information Dependent Mode. MS spectra were acquired across the mass range of 400–1250 m/z in high resolution mode (>30,000) using 250 ms accumulation time per spectrum. A maximum of 10 precursors per cycle were chosen for fragmentation from each MS spectrum with 100 ms minimum accumulation time for each precursor and dynamic exclusion for 8 s. Tandem mass spectra were recorded in high sensitivity mode (resolution > 15,000) with rolling collision energy on adjustment. Survey-IDA experiment with charge states 2 to 4, which exceeds 125 cps, was selected. Peptide identification was achieved with ProteinPilot 5.0 software Revision 4769 (AB SCIEX, Redwood City, CA, USA) using the Paragon database search algorithm (5.0.0.0.4767) for peptide identification and the integrated false discovery rate (FDR) analysis function. The data were searched against a database consisting of SerpentesDB database (total 345,092 entries). The search parameters are as follows: Sample Type: Identification; Cys Alkylation: Iodoacetamide; Digestion: Trypsin; Special Factors: None; Species: None. The processing was specified as follows: ID Focus: Biological Modifications; Search Effort: Thorough; Detected Protein Threshold: 0.05 (10.0%); and Competitor Error Margin (ProtScore): 2.00.

4.3.3. Electrospray Mass Spectrometry of Digested Venom

Approximately 300 ng of trypsin-digested venom was applied to a nano-chromatography system (Ultimate 3000 system, Thermo Scientific Dionex, Sunnyvale, CA, USA) connected to an electrospray mass spectrometer (Q Exactive™ Plus Hybrid Quadrupole-Orbitrap™, Thermo Scientific, Waltham, MA, USA). The chromatography consisted of an initial desalting and concentration step of the sample using a capillary “trap” column (2 cm × 100 m i.d.) packed with 5 m, 200 Å Magic C18 AQ matrix (Michrom Bioresources, Auburn, CA, USA). The fractionation of the peptides was performed on a 30 cm analytical column (75 m i.d., packed with 1.9 m ReproSil-Pur 120 C18-AQ; Dr. Maisch) including a laser-pulled tip (~5 m). Using standard mobile phases (solvent A: 0.1% formic acid in water; solvent B: 0.1% formic acid in acetonitrile), we applied a shallow gradient of 2–40% of B over 160 min for maximum peptide resolution. To elute possible hydrophobic peptides, the gradient was increased to 80% B over 4 min at the end of the run. The capillary temperature was set to 250 °C and a spray voltage of 1.9 kV with no sheath or auxiliary gas flow was applied. During the acquisition of the mass spectrometric data, up to 12 of the most intense precursor ions were selected for fragmentation using higher energy collisional dissociation (HCD). Normalized collision energy of 30% was applied, and singly-charged peptides or peptides with unresolved charge states were excluded. The acquisition parameters for full MS were as follows: resolution 70,000 (FWHM, full width at half-maximum, at m/z 200), automatic gain control pre-scan (AGC) target 1 × 106, maximum ion injection time (IT) of 100 ms, scan range 300–1500 m/z, and profile spectrum data type. Data-dependent (dd) MS/MS acquisition parameters were as follows: resolution 17,500, AGC target 5 × 104, maximum injection time 50 ms, isolation window 2.0 m/z (with 0.5 m/z offset). Dynamic exclusion time was set to 60 s. For external calibration, a mixture of caffeine, peptide MRFA (Met–Arg–Phe–Ala) and Ultramark 1621 was used following the manufacturer’s instructions.

4.3.4. MALDI Mass Spectrometry of Digested Venom

The tryptic digest was separated on a C18 reversed phase column (Agilent Zorbax 300SB-C18 1.0 × 150 mm × 3.5 m) by running a linear gradient from 0% to 50% acetonitrile over 120 min applying a flow rate of 40 L/min (solvent A 0.1% TFA, solvent B contains 80% ACN with 0.1% TFA) on an Eksigent system (Dionex, Vernon Hills, IL, USA). During the chromatographic run, approximately 200 fractions were manually collected directly onto the MALDI plate and 0.3 µL of a saturated matrix solution (10 mg/mL α-cyano-4-hydroxycinnamic acid (Aldrich, Milwaukee, WI, USA) in 50% acetonitrile with 0.1% TFA was added. Finally, the solution was dried at room temperature (dried-droplet method). Raw data for protein identification were obtained on an AB Sciex 5800 (AB Sciex, Foster City, CA, USA). Up to 15 of the most intense ion signals with a signal-to-noise ratio above 2 were selected as precursors for MS/MS acquisition (excluding common trypsin autolysis and keratin peaks). External calibration in MS mode was performed using a mixture of six singly charged peptides: des-Arg1-Bradykinin (m/z = 904.468), angiotensin I (m/z = 1296.685), Glu1-fibrinopeptideB (m/z = 1570.677), ACTH (1–17 clip; m/z = 2093.087), ACTH (18–39 clip; m/z = 2465.199), and ACTH (7–38 clip; m/z = 3657.929). MS/MS spectra were externally calibrated using known fragment ion masses observed in the tandem mass spectrum of Glu1-fibrinopeptideB.

4.4. Data Analysis of Digested Venom

Database searches of the mass spectra acquired by MALDI or electrospray MS were searched against a database compiled of venom sequences deposited on NCBI (71798 entries, 8 June 2015) using Mascot version 2.1 (Matrix Science Inc., Boston, MA, USA). The following search parameters were applied: no restrictions on species of origin or protein molecular weight, two tryptic missed cleavages allowed, variable modifications of cysteine (carbamidomethylation) and methionine (oxidation), and pyro-glutamate formation at N-terminal glutamine of peptides. The enzyme conditions were set to “semi-tryptic” including the option to search for non-tryptic cleavage sites either at the N-terminus or C-terminus of the corresponding peptide. The database search of electrospray data (Orbitrap Q Exactive™, ThermoFisher Scientific, Waltham, MA, USA) applying a p-value of 0.05% showed a false discovery rate (FDR) of 1.82%, while the database search of spectra acquired with MALDI mass spectrometry revealed a FDR of 0.49% (p-value 0.05%).

4.5. De Novo Sequencing and Similarity-Driven Analysis of Digested Venom

In addition to automated database searches, all tandem mass spectra were sequenced by de novo analysis using the program PEAKS version 6 (Bioinformatics Solutions Inc., Waterloo, KW, Canada). The following search parameters were applied: cysteine modification (carbamidomethylation of cysteine), deamidation of asparagine (N), oxidation of methionine, and pyroglutamate modification at N-terminal glutamate; 20 ppm precursor mass tolerance and 0.1 Da fragment mass tolerance. High quality de novo sequences (Peaks ALC score higher than 50 and sequences with six or more amino acid residues) were submitted to the algorithm PepExplorer [64] to identify possible snake venom peptides not detected by automated database searches (against a snake venom database including 71,798 entries). The search included the PAM (point accepted mutation) 30MS substitution matrix to score the alignments. Initially, peptides with at least 65% of sequence identity compared with a sequence from the database were accepted. Finally, tandem mass spectra of sequences that satisfied the above mentioned criteria were manually confirmed

4.6. Protein Quantitation

Relative protein abundances were estimated by normalized spectral-counts of the corresponding MS/MS spectra, referred to as the normalized spectral abundance factor (NSAF) [65].

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2072-6651/10/9/373/s1, Table S1: Proteins identified in different Coomassie-stained bands of a 1D gel of the crude venom of B. flaviceps, Table S2: Proteins identified in the venom of B. flaviceps including detailed information on peptide sequences, m/z values, and ppm values.

Author Contributions

Conceptualization, A.C., S.P.M. and R.M.K.; Methodology, A.C., S.P.M. and R.J.R.M..; Software, P.C.; Validation, A.C., S.P.M., R.M.K. and R.J.R.M. Formal Analysis, A.C., S.P.M., C.M.M., P.C., A.S. and J.P.; Resources, A.C. and S.P.M.; Writing-Original Draft Preparation, A.C. and S.P.M.; Writing-Review & Editing, A.C., S.P.M., R.M.K., C.M.M., R.J.R.M., P.C., A.S. and J.P.; Project Administration, A.C., S.P.M. and R.M.K.; Funding Acquisition, S.P.M.

Funding

Funding for this study was provided in part by a UNC Provost Grant (QR088) to SPM.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Brahma, R.K.; McCleary, R.J.; Kini, R.M.; Doley, R. Venom gland transcriptomics for identifying, cataloging, and characterizing venom proteins in snakes. Toxicon 2015, 93, 1–10. [Google Scholar] [CrossRef] [PubMed]
  2. Ducancel, F.; Durban, J.; Verdenaud, M. Transcriptomics and venomics: Implications for medicinal chemistry. Future Med. Chem. 2014, 6, 1629–1643. [Google Scholar] [CrossRef] [PubMed]
  3. Calvete, J.J. Snake venomics: From the inventory of toxins to biology. Toxicon 2013, 75, 44–62. [Google Scholar] [CrossRef] [PubMed]
  4. Zelanis, A.; Tashima, A.K. Unraveling snake venom complexity with ‘omics’ approaches: Challenges and perspectives. Toxicon 2014, 87, 131–134. [Google Scholar] [CrossRef] [PubMed]
  5. Calvete, J.J.; Sanz, L.; Pla, D.; Lomonte, B.; Gutierrez, J.M. Omics meets biology: Application to the design and preclinical assessment of antivenoms. Toxins 2014, 6, 3388–3405. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Pla, D.; Petras, D.; Saviola, A.J.; Modahl, C.M.; Sanz, L.; Pérez, A.; Juárez, E.; Frietze, S.; Dorrestein, P.C.; Mackessy, S.P.; et al. Transcriptomics-guided bottom-up and top-down venomics of neonate and adult specimens of the arboreal rear-fanged Brown Treesnake, Boiga irregularis, from Guam. J. Proteom. 2018, 174, 71–84. [Google Scholar] [CrossRef] [PubMed]
  7. Viala, V.L.; Hildebrand, D.; Trusch, M.; Fucase, T.M.; Sciani, J.M.; Pimenta, D.C.; Arni, R.K.; Schlüter, H.; Betzel, C.; Mirtschin, P.; et al. Venomics of the Australian eastern brown snake (Pseudonaja textilis): Detection of new venom proteins and splicing variants. Toxicon 2015, 107, 252–265. [Google Scholar] [CrossRef] [PubMed]
  8. Reeks, T.; Lavergne, V.; Sunagar, K.; Jones, A.; Undheim, E.; Dunstan, N.; Fry, B.; Alewood, P.F. Deep venomics of the Pseudonaja genus reveals inter- and intra-specific variation. J. Proteom. 2016, 133, 20–32. [Google Scholar] [CrossRef] [PubMed]
  9. Modahl, C.M.; Frietze, S.; Mackessy, S.P. Adaptive evolution of distinct prey-specific toxin genes in rear-fanged snake venom. Proc. Biol. Sci. 2018, 285. [Google Scholar] [CrossRef]
  10. Chatrath, S.T.; Chapeaurouge, A.; Lin, Q.; Lim, T.K.; Dunstan, N.; Mirtschin, P.; Kumar, P.P.; Kini, R.M. Identification of novel proteins from the venom of a cryptic snake Drysdalia coronoides by a combined transcriptomics and proteomics approach. J. Proteome Res. 2011, 10, 739–750. [Google Scholar] [CrossRef] [PubMed]
  11. Pawlak, J.; Mackessy, S.P.; Sixberry, N.M.; Stura, E.A.; Le Du, M.H.; Menez, R.; Foo, C.S.; Menez, A.; Nirthanan, S.; Kini, R.M. Irditoxin, a novel covalently linked heterodimeric three-finger toxin with high taxon-specific neurotoxicity. FASEB J. 2009, 23, 534–545. [Google Scholar] [CrossRef] [PubMed]
  12. Heyborne, W.H.; Mackessy, S.P. Identification and characterization of a taxon-specific three-finger toxin from the venom of the Green Vinesnake (Oxybelis fulgidus; family Colubridae). Biochimie 2013, 95, 1923–1932. [Google Scholar] [CrossRef] [PubMed]
  13. James, J. The Snake Charmer: A Life and Death in Pursuit of Knowledge; Hyperion Press: New York, NY, USA, 2008; p. 288. [Google Scholar]
  14. Jiang, Y.; Li, Y.; Lee, W.; Xu, X.; Zhang, Y.; Zhao, R.; Zhang, Y.; Wang, W. Venom gland transcriptomes of two elapid snakes (Bungarus multicinctus and Naja atra) and evolution of toxin genes. BMC Genom. 2011, 12, 1. [Google Scholar] [CrossRef] [PubMed]
  15. Shan, L.L.; Gao, J.F.; Zhang, Y.X.; Shen, S.S.; He, Y.; Wang, J.; Ma, X.M.; Ji, X. Proteomic characterization and comparison of venoms from two elapid snakes (Bungarus multicinctus and Naja atra) from China. J. Proteom. 2016, 138, 83–94. [Google Scholar] [CrossRef] [PubMed]
  16. Oh, A.M.F.; Tan, C.H.; Ariaranee, G.C.; Quraishi, N.; Tan, N.H. Venomics of Bungarus caeruleus (Indian krait): Comparable venom profiles, variable immunoreactivities among specimens from Sri Lanka, India and Pakistan. J. Proteom. 2017, 164, 1–18. [Google Scholar] [CrossRef] [PubMed]
  17. Das, I. A Field Guide to the Reptiles of South-East Asia; New Holland Publishers: London, UK, 2010; p. 376. [Google Scholar]
  18. Chang, C.C.; Lee, C.Y. Isolation of neurotoxins from the venom of Bungarus multicinctus and their modes of neuromuscular blocking action. Arch. Int. Pharmacodyn. Ther. 1963, 144, 241–257. [Google Scholar] [PubMed]
  19. Chu, C.C.; Li, S.H.; Chen, Y.H. Resolution of isotoxins in the beta-bungarotoxin family. J. Chromatogr. A 1995, 694, 492–497. [Google Scholar] [CrossRef]
  20. Su, M.J.; Chang, C.C. Presynaptic effects of snake venom toxins which have phospholipase A2 activity (beta-bungarotoxin, taipoxin, crotoxin). Toxicon 1984, 22, 631–640. [Google Scholar] [CrossRef]
  21. Rusmili, M.R.A.; Yee, T.T.; Mustafa, M.R.; Hodgson, W.C.; Othman, I. Proteomic characterization and comparison of Malaysian Bungarus candidus and Bungarus fasciatus venoms. J. Proteom. 2014, 110, 129–144. [Google Scholar] [CrossRef] [PubMed]
  22. McDowell, R.S.; Dennis, M.S.; Louie, A.; Shuster, M.; Mulkerrin, M.G.; Lazarus, R.A.; Ma, Z.Q. Mambin, a potent glycoprotein IIb-IIIa antagonist and platelet aggregation inhibitor structurally related to the short neurotoxins. Biochemistry 1992, 31, 4766–4772. [Google Scholar] [CrossRef] [PubMed]
  23. Chaki, S.; Muramatsu, M.; Ushiyama, Y.; Otomo, S. Purification and partial characterization of K+ channel blockers from the venom of Dendroaspis angusticeps. Neurochem. Int. 1992, 20, 553–558. [Google Scholar] [CrossRef]
  24. Kini, R.M.; Doley, R.S. Structure, function and evolution of three-finger toxins: Mini proteins with multiple targets. Toxicon 2010, 56, 855–867. [Google Scholar] [CrossRef] [PubMed]
  25. Siang, A.S.; Doley, R.; Vonk, F.J.; Kini, R.M. Transcriptomic analysis of the venom gland of the red-headed krait (Bungarus flaviceps) using expressed sequence tags. BMC Mol. Biol. 2010, 11, 11–15. [Google Scholar] [CrossRef] [PubMed]
  26. Ziganshin, R.H.; Kovalchuk, S.I.; Arapidi, G.P.; Starkov, V.; Hoang, A.N.; Thi Nguyen, 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]
  27. Laskowski, M.J.; Kato, I. Protein inhibitors of proteinases. Annu. Rev. Biochem. 1980, 49, 593–626. [Google Scholar] [CrossRef] [PubMed]
  28. Laskowski, M.J. Protein inhibitors of serine proteinases—Mechanism and classification. Adv. Exp. Med. Biol. 1986, 199, 1–17. [Google Scholar] [PubMed]
  29. Chang, L.S.; Chung, C.; Huang, H.B.; Lin, S.R. Purification and characterization of a chymotrypsin inhibitor from the venom of Ophiophagus hannah (King Cobra). Biochem. Biophys. Res. Commun. 2001, 283, 862–867. [Google Scholar] [CrossRef] [PubMed]
  30. Cheng, Y.C.; Yan, F.R.; Chang, L.S. Taiwan cobra chymotrypsin inhibitor: Cloning, functional expression and gene organization. Biochim. Biophys. Acta 2005, 1747, 213–220. [Google Scholar] [CrossRef] [PubMed]
  31. Kini, R.M. Structure-function relationships and mechanism of anticoagulant phospholipase A2 enzymes from snake venoms. Toxicon 2005, 45, 1147–1161. [Google Scholar] [CrossRef] [PubMed]
  32. Sunderkotter, C.; Steinbrink, K.; Goebeler, M.; Bhardwaj, R.; Sorg, C. Macrophages and angiogenesis. J. Leukocyte Biol. 1994, 55, 410–422. [Google Scholar] [CrossRef] [PubMed]
  33. Verheul, H.M.W.; Hoekman, K.; Bakker, S.L.D.; Eekman, C.A.; Folman, C.C.; Broxterman, H.J.; Pinedo, H.M. Platelet: Transporter of vascular endothelial growth factor. Clin. Cancer Res. 1997, 3, 2187–2190. [Google Scholar] [PubMed]
  34. Itakura, J.; Ishiwata, T.; Shen, B.; Kornmann, M.; Korc, M. Concomitant over-expression of vascular endothelial growth factor and its receptors in pancreatic cancer. Int. J. Cancer 2000, 85, 27–34. [Google Scholar] [CrossRef] [Green Version]
  35. Yamazaki, Y.; Takani, K.; Atoda, H.; Morita, T. Snake venom vascular endothelial growth factors (VEGFs) exhibit potent activity through their specific recognition of KDR (VEGF receptor 2). J. Biol. Chem. 2003, 278, 51985–51988. [Google Scholar] [CrossRef] [PubMed]
  36. Yamazaki, Y.; Matsunaga, Y.; Tokunaga, Y.; Obayashi, S.; Saito, M.; Morita, T. Snake venom vascular endothelial growth factors (VEGF-Fs) exclusively vary their structures and functions among species. J. Biol. Chem. 2009, 284, 9885–9891. [Google Scholar] [CrossRef] [PubMed]
  37. Aird, S.D. Ophidian envenomation strategies and the role of purines. Toxicon 2002, 40, 335–393. [Google Scholar] [CrossRef]
  38. Margres, M.J.; Aronow, K.; Loyacano, J.; Rokyta, D.R. The venom-gland transcriptome of the eastern coral snake (Micrurus fulvius) reveals high venom complexity in the intragenomic evolution of venoms. BMC Genom. 2013, 14, 531. [Google Scholar] [CrossRef] [PubMed]
  39. McGivern, J.J.; Wray, K.P.; Margres, M.J.; Couch, M.E.; Mackessy, S.P.; Rokyta, D.R. RNA-seq and high-definition mass spectrometry reveal the complex and divergent venoms of two rear-fanged colubrid snakes. BMC Genom. 2014, 15, 1061. [Google Scholar] [CrossRef] [PubMed]
  40. Lu, Q.M.; Clemetson, J.M.; Clemetson, K.J. Snake venom C-type lectins interacting with platelet receptors. Toxin Rev. 2007, 26, 77–93. [Google Scholar] [CrossRef]
  41. Zha, H.G.; Lee, W.H.; Zhang, Y. Cloning of cDNAs encoding C-type lectins from Elapidae snakes Bungarus fasciatus and Bungarus multicinctus. Toxicon 2001, 39, 1887–1892. [Google Scholar] [CrossRef]
  42. 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. Biol. Chem. 2005, 280, 13137–13147. [Google Scholar] [CrossRef] [PubMed]
  43. Fox, J.W. A brief review of the scientific history of several lesser-known snake venom proteins: L-amino acid oxidases, hyaluronidases and phosphodiesterases. Toxicon 2013, 62, 75–82. [Google Scholar] [CrossRef] [PubMed]
  44. Trummal, K.; Aaspollu, A.; Tonismagi, K.; Samel, M.; Subbi, J.; Siigur, J.; Siigur, E. Phosphodiesterase from Vipera lebetina venom—Structure and characterization. Biochimie 2014, 106, 48–55. [Google Scholar] [CrossRef] [PubMed]
  45. Lavin, M.F.; Earl, S.; Birrell, G.; St. Pierre, L.; Guddat, L.; de Jersey, J.; Masci, P. Snake venom nerve growth factors. In Handbook of Venoms and Toxins of Reptiles; Mackessy, S.P., Ed.; CRC Press: Boca Raton, FL, USA, 2009; pp. 377–391. [Google Scholar]
  46. 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]
  47. Chen, H.S.; Wang, Y.M.; Huang, W.T.; Huang, K.F.; Tsai, I.H. Cloning, characterization and mutagenesis of Russell’s viper venom L-amino acid oxidase: Insights into its catalytic mechanism. Biochimie 2012, 94, 335–344. [Google Scholar] [CrossRef] [PubMed]
  48. Zuliani, J.P.; Kayano, A.M.; Zaqueo, K.D.; Neto, A.C.; Sampaio, S.V.; Soares, A.M.; Stabeli, R.G. Snake venom L-amino acid oxidases: Some consideration about their functional characterization. Protein Pept. Lett. 2009, 16, 908–912. [Google Scholar] [CrossRef] [PubMed]
  49. Mackessy, S.P.; Baxter, L.M. Bioweapons synthesis and storage: The venom gland of front-fanged snakes. Zool. Anz. 2006, 245, 147–159. [Google Scholar] [CrossRef]
  50. Yamazaki, Y.; Brown, R.L.; Morita, T. Purification and cloning of toxins from elapid venoms that target cyclic nucleotide-gated ion channels. Biochemistry 2002, 41, 11331–11337. [Google Scholar] [CrossRef] [PubMed]
  51. Chapeaurouge, A.; Abu Reza, M.D.; Mackessy, S.P.; Carvalho, P.C.; Valente, R.H.; Teixeira-Ferreira, A.; Jonas 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. [Google Scholar] [CrossRef] [PubMed]
  52. Markland, F.S.; Swenson, S. Snake venom metalloproteinases. Toxicon 2013, 62, 3–18. [Google Scholar] [CrossRef] [PubMed]
  53. Fox, J.W.; Serrano, S.M. Insights into and speculations about snake venom metalloproteinase (SVMP) synthesis, folding and disulfide bond formation and their contribution to venom complexity. FEBS J. 2008, 275, 3016–3030. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Doley, R.; Mackessy, S.P.; Kini, R.M. Role of accelerated segment switch in exons to alter targeting (ASSET) in the molecular evolution of snake venom proteins. BMC Evol. Biol. 2009, 9, 146. [Google Scholar] [CrossRef] [PubMed]
  55. Zhang, Y.; Xiong, Y.L.; Bon, C. An activator of blood-coagulation factor-X from the venom of Bungarus fasciatus. Toxicon 1995, 33, 1277–1288. [Google Scholar] [CrossRef]
  56. Johns, D.G.; Ao, Z.H.; Heidrich, B.J.; Hunsberger, G.E.; Graham, T.; Payne, L.; Elshourbagy, N.; Lu, Q.; Aiyar, N.; Douglas, S.A. Dendroaspis natriuretic peptide binds to the natriuretic peptide clearance receptor. Biochem. Biophys. Res. Commun. 2007, 358, 145–149. [Google Scholar] [CrossRef] [PubMed]
  57. Trinh, K.X.; Khac, Q.L.; Trinh, L.X.; Warren, D.A. Hyponatraemia, rhabdomyolysis, alterations in blood pressure and persistent mydriasis in patients envenomed by Malayan kraits (Bungarus candidus) in southern Viet Nam. Toxicon 2010, 56, 1070–1075. [Google Scholar] [CrossRef] [PubMed]
  58. Sridharan, S.; Kini, R.M. Snake venom natriuretic peptides: Potential molecular probes. BMC Pharmacol. Toxicol. 2015, 16. [Google Scholar] [CrossRef] [Green Version]
  59. Vink, S.; Jin, A.H.; Poth, K.J.; Head, G.A.; Alewood, P.F. Natriuretic peptide drug leads from snake venom. Toxicon 2012, 59, 434–445. [Google Scholar] [CrossRef] [PubMed]
  60. Rehana, S.; Kini, R.M. Molecular isoforms of cobra venom factor-like proteins in the venom of Austrelaps superbus. Toxicon 2007, 50, 32–52. [Google Scholar] [CrossRef] [PubMed]
  61. Vogel, C.W.; Bredehorst, R.; Fritzinger, D.C.; Grunwald, T.; Ziegelmuller, P.; Kock, M.A. Structure and function of cobra venom factor, the complement-activating protein in cobra venom. Adv. Exp. Med. Biol. 1996, 391, 97–114. [Google Scholar] [PubMed]
  62. Guan, H.H.; Goh, K.S.; Davamani, F.; Wu, P.L.; Huang, Y.W.; Jeyakanthan, J.; Wu, W.; Chen, C.J. Structures of two elapid snake venom metalloproteases with distinct activities highlight the disulfide patterns in the D domain of ADAMalysin family proteins. J. Struct. Biol. 2010, 169, 294–303. [Google Scholar] [CrossRef] [PubMed]
  63. Pla, D.; Sanz, L.; Sasa, M.; Acevedo, M.E.; Dwyer, Q.; Durban, J.; Pérez, A.; Rodriguez, Y.; Lomonte, B.; Calvete, J.J. Proteomic analysis of venom variability and ontogeny across the arboreal palm-pitvipers (genus Bothriechis). J. Proteom. 2017, 152, 1–12. [Google Scholar] [CrossRef] [PubMed]
  64. Leprevost, F.V.; Valente, R.H.; Lima, D.B.; Perales, J.; Melani, R.; Yates, J.R., III; Barbosa, V.C.; Junqueira, M.; Carvalho, P.C. PepExplorer: A similarity-driven tool for analyzing de novo sequencing results. Mol. Cell. Proteom. 2014, 13, 2480–2489. [Google Scholar] [CrossRef] [PubMed]
  65. Florens, L.; Carozza, M.J.; Swanson, S.K.; Fournier, M.; Coleman, M.K.; Workman, J.L.; Washburn, M.P. Analyzing chromatin remodeling complexes using shotgun proteomics and normalized spectral abundance factors. Methods 2006, 40, 303–311. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Approximate distribution of the Red-headed Krait (Bungarus flaviceps). Inset A shows the extraction of venom from B. flaviceps, a low-yielding front-fanged snake. Distributions are adapted from the literature ([17], 2010, New Holland Publishers) and the map is from Google Earth.
Figure 1. Approximate distribution of the Red-headed Krait (Bungarus flaviceps). Inset A shows the extraction of venom from B. flaviceps, a low-yielding front-fanged snake. Distributions are adapted from the literature ([17], 2010, New Holland Publishers) and the map is from Google Earth.
Toxins 10 00373 g001
Figure 2. One-dimensional gel of the crude venom of B. flaviceps. Prominent venom protein families (boxed) are based on their molecular masses and confirmed by liquid chromatography (LC)-tandem mass spectrometry (MS/MS) (see also Table S1). Bands were excised for further analysis. PDE—phosphodiesterases; AChE—acetylcholinesterases; LAAO—L-amino-acid oxidase; 5’-NUC—5’inucleotidase; SVMP—snake venom metalloproteinase; PLB—phospholipase B; CRISP—cysteine-rich seceretory proteins; PLA2—phospholipase A; KUN—Kunitz-type serine proteinase inhibitor; 3FTx—three-finger toxin.
Figure 2. One-dimensional gel of the crude venom of B. flaviceps. Prominent venom protein families (boxed) are based on their molecular masses and confirmed by liquid chromatography (LC)-tandem mass spectrometry (MS/MS) (see also Table S1). Bands were excised for further analysis. PDE—phosphodiesterases; AChE—acetylcholinesterases; LAAO—L-amino-acid oxidase; 5’-NUC—5’inucleotidase; SVMP—snake venom metalloproteinase; PLB—phospholipase B; CRISP—cysteine-rich seceretory proteins; PLA2—phospholipase A; KUN—Kunitz-type serine proteinase inhibitor; 3FTx—three-finger toxin.
Toxins 10 00373 g002
Figure 3. (A) Abundances of the venom protein families of B. flaviceps as evidenced by normalized mass spectrometric spectral count. (B) Comparison of abundances of venom protein families (B. flaviceps) by transcriptomic (red, adapted from [25], 2010, Springer Nature) and proteomic analysis. KSPI—Kunitz-type serine proteinase inhibitors; VEGF—vascular endothelial growth factor; VNGF—venom nerve growth factor; SVSP—snake venom serine proteinase; CVF—cobra venom factor.
Figure 3. (A) Abundances of the venom protein families of B. flaviceps as evidenced by normalized mass spectrometric spectral count. (B) Comparison of abundances of venom protein families (B. flaviceps) by transcriptomic (red, adapted from [25], 2010, Springer Nature) and proteomic analysis. KSPI—Kunitz-type serine proteinase inhibitors; VEGF—vascular endothelial growth factor; VNGF—venom nerve growth factor; SVSP—snake venom serine proteinase; CVF—cobra venom factor.
Toxins 10 00373 g003
Table 1. Snake venom protein families of B. flaviceps identified by automated database search. SVMP—snake venom metalloproteinase; SVSP—snake venom serine proteinase.
Table 1. Snake venom protein families of B. flaviceps identified by automated database search. SVMP—snake venom metalloproteinase; SVSP—snake venom serine proteinase.
Protein FamilyProteinAccession No.SpeciesNumber of Peptides Matched
3FTxNon-conventional three finger toxin isoform 1294961050Bungarus flaviceps6
3FTxNon-conventional three finger toxin isoform 6294961060Bungarus flaviceps5
3FTxShort-chain three finger toxin isoform 4294961042Bungarus flaviceps4
3FTxShort-chain three finger toxin isoform 7294961048Bungarus flaviceps9
3FTxShort-chain three finger toxin isoform 6294961046Bungarus flaviceps7
3FTxShort-chain three finger toxin isoform 1294961036Bungarus flaviceps2
3FTxκ-bungarotoxin809178Bungarus multicinctus1
3FTxShort-chain three finger toxin isoform 3294961040Bungarus flaviceps1
3FTxκ-flavitoxin128938Bungarus flaviceps9
3FTxMuscarinic toxin-like protein294961066Bungarus flaviceps9
Serine protease inhibitorβ-bungarotoxin B chain precursor31745053Bungarus flaviceps7
Serine protease inhibitorKunitz-type serine proteinase inhibitor isoform 5 294961076Bungarus flaviceps8
Serine protease inhibitorKunitz-type serine proteinase inhibitor isoform 1294961068 Bungarus flaviceps4
AcetylcholinesteraseAcetylcholinesterase1389604Bungarus flaviceps20
AcetylcholinesteraseAcetylcholinesterase DEN476538388Denisoni adevisi4
PLA2β-bungarotoxin A2 chain precursor31745049Bungarus flaviceps16
PLA2β-bungarotoxin A1 chain precursor 31745051Bungarus flaviceps3
PLA2Phospholipase A2II precursor 31745057Bungarus flaviceps21
PLA2Phospholipase A2263083 Bungarus fasciatus8
PLA2Phospholipase A2 precursor 31745055Bungarus flaviceps23
PLA2Phospholipase A2 isoform 3 294961092Bungarus flaviceps5
PLA2Phospholipase A2Kbf-III 110559306Bungarus fasciatus4
PLA2Phospholipase A2 isozyme 124638470Laticauda semifasciata1
PLA2Phospholipase A229422777Bungarus candidus2
PLA2Phospholipase A25924345Austrelaps superbus3
PLA2Phospholipase A2152032644Bungarus fasciatus3
PLA2Phospholipase A2 precursor 156257593Bungarus fasciatus1
PLA2Phospholipase A2 48425218Bungarus caeruleus1
PLA2Phospholipase A2 129428Laticauda colubrina1
VEGFHypothetical protein L345_04144565318860Ophiophagus hannah4
NucleotidaseEcto-5′-nucleotidase 1 537444870Micrurus fulvius11
SnaclecSnaclec factor IX/factor X-binding protein B chain398488Bothrops jararaca1
SnaclecC-type lectin-like protein 1 13876735Bungarus fasciatus3
VesprynOhanin precursor70907886Ophiophagus hannah2
VesprynVespryn22336042222Drysdalia coronoides1
PhosphodiesterasePhosphodiesterase 1537444868Micrurus fulvius20
PhosphodiesterasePhosphodiesterase 1 338855302Crotalus adamanteus1
HyaluronidaseHyaluronidase 113203681Bitis arietans4
VNGFVenom nerve growth factor precursor266299Bungarus multicinctus4
LAAOL-amino-acid oxidase126035653Bungarus fasciatus5
LAAOL-amino acid oxidase126035649Bungarus multicinctus2
LAAOL-amino-acid oxidase426205815Crotalus durissus cumanensis2
PLBPhospholipase B 537444729Micrurus fulvius3
CRISPCysteine-rich seceretory protein190195343Bungarus candidus3
CRISPOpharin precursor225547744Ophiophagus hannah2
SVMPScutatease-1 (PIII)145982766Notechis scutatus8
SVMPMetalloproteinase (PIII)126035640Bungarus multicinctus6
SVMPMetalloproteinase MTP9 (PIII)336042214Drysdalia coronoides4
SVMPMetalloproteinase (PIII)126035635Bungarus fasciatus4
SVMPP-III633276509Micropechis ikaheka4
SVMPMTP4 (PIII)537463069Micrurus fulvius3
SVMPAtragin precursor(PIII)224482347Naja atra3
SVMPMetalloproteinase isoform 3 (PIII)109254964Sistrurus catenatus edwardsi2
SVMPSVMP-Hop-14, partial (PIII)476539284Hoplocephalus bungaroides2
SVMPSVMP-Hop-46, partial(PIII)476539268Hoplocephalus bungaroides2
SVMPSVMP 1537444726Micrurus fulvius2
SVMPMetalloproteinase (PII)82466485Bothrops asper1
SVMPFur-1, partial (PI)476538467Furinaor nata1
SVMPjararhagin (PIII)62468Bothrops jararaca1
SVMPMetalloproteinase (PIII)241995585Philodrya solfersii1
SVMPLeucurolysin-B (PIII)223635807Bothrops leucurus1
SVMPEch-32 (PIII)476538400Echiopsis curta1
SVMPCobrin precursor(PIII)6006966Naja naja1
SVMPMetalloproteinase (PII)297594122Echis pyramidum leakeyi1
SVMPCohPH-3 (PII)522802426Crotalus oreganus helleri1
SVSPSerine proteinase isoform 2 109254940Sistrurus catenatus edwardsi2
SVSPSVSP 11387014258Crotalus adamanteus1
Natriuretic peptideNatriuretic peptide294961100Bungarus flaviceps1
Complement-depleting factorComplement-depleting factor126035660Bungarus fasciatus1
Table 2. Snake venom protein families of B. flaviceps identified by de novo sequencing of tandem mass spectra followed by sequence-similarity analysis. The corresponding sequences are indicated.
Table 2. Snake venom protein families of B. flaviceps identified by de novo sequencing of tandem mass spectra followed by sequence-similarity analysis. The corresponding sequences are indicated.
Protein FamilyProteinAccession No.SpeciesDatabase SequenceAlignment ScoreDe Novo SequenceALC (%)m/zppm to De Novo Derived Sequence
3FTxMuscarinic toxin-like protein 117606606Ophiophagus hannahTFTCPELTPN64CGTYTCPELTPDRR52575.935413.5
PLA2Phospholipase A29453880Laticauda semifasciataVHDDCYGEAEK87VHDDCYGEAEK58563.5938−1.8
LVQFTYLIQCANKGSRASYHYADYGCY84RLAQFACLLQCADEEVGVDLHYADYGCYNCK46906.9115.8
NucleotidaseSnake venom 5′-nucleotidase 752779244Boiga irregularisVVQFMNSLR 69EAAAVVQFMNSLR51479.25377.1
HGQGTGELLQVSGIKVVYDLS99HPQGTGELLQVSLGKVVYDLSKGAPVNK52587.72924.2
SnaclecC-type lectin 2 537444936Micrurus fulviusDHHCPSDWYSFDKFCYKFI105DNHCNADWSSFDKFCYQFLR52653.27992
VesprynVespryn 23 336042224Drysdalia coronoidesHFFEVK49HFFEVK60403.7129−1.1
IVVFLDYSEGK68LVVFLDYKEGK60437.5828−1.2
IVVFLDY51LVVFLDYK71498.791−1.7
PhosphodiesterasePhosphodiesterase675402318Protobothrops elegansFYTLYIEEPDTTGHK106QFLPVFFTLYLEEPDTTGHK57592.05337
RPDFYTLYIEEPDTTGHK130RPDFFTLYLEEPDTTGHK59542.2665−2.6
VNGFVenom nerve growth factor335892638Bungarus fasciatusYFFETK51AGYFFETK51481.7344−0.2
LAAOL-amino acid oxidase 126035653Bungarus fasciatusSYVTADYVIVCAT76LTSYVTADYLLVSATMGR51981.0001−6.1
TADYVIVCATSR80TLSAAHTADYLLVCATSR52650.672314.1
LAAOL-amino acid oxidase 401021343Lachesis mutaLSAAYVLAEAGHQVTVLEASER105LSAAYVLSEKHQVTVLQASER52583.3127−2.5
LAAOL-amino acid oxidase 537444909Micrurus fulviusSAAYVLAEAGHKVTLLE103QLLVVQGMQDSAAYVLSEAGHKVTVLESNNK50669.762920.9
AGMAGLSAAYVLAEAGHK116VKHHSVPPRVAGMAGLSAAYVLSEAGHK52574.716910.5
LAAOL-amino acid oxidase 3426324Crotalus adamanteusRVVIVGAGMAGLSAAYVL104RLALVGAGMAGLSAAYVLNPSSPGVLSQLLWSR53671.7708−3.1
LAAOL-amino acid oxidase 1a 537444909Micrurus fulviusGMAGLSAAYVLAEAGHK110RTPVVKGMAGLSAAYVLSEAGHK48590.3197−1.2
Phospholipase BPhospholipase B 537444729Micrurus fulviusKGYWPSYNIPFHK94KPYWPSYNLPFHK56559.6232−1.9
CRISPCRISP precursor158262802Austrelaps superbusVDKHNALR58VAAVDKHNALR45398.5624−2
YLYVCQYCPAGNIRGSIATPYK84YLYVCQYCFRAFGAPYNNKHGSLATPYK43847.8989−8.6
CRISPCysteine-rich secretory protein A 698375481Opheodrys aestivusWNSNAAQNAK76EMWNSNAAQNAK44682.3055−1.5
PAGNIVGSIATPYK95MENTQAAVCSDPAGNLVGSLATPYK44846.41147.4
PAGNIVGSIATPYK95LYYCMNEGQMAPAGNLVGSLATPYK44897.7457−14.4
ASCFCH50HCASCFCHGGK48632.2381−10.1
SVSPVenom thrombin-like enzyme 118500915Deinagkistrodon acutusDECDINEHR 72VVNDECDLNEHR51500.5556−1.3
SVMPMetalloproteinase 13b 699236656Hypsiglena sp. JMG-2014DPDNGMVEPGTK82LNMMSSSPNACADPDEGMVEPGTK50846.6896−2.9
SVMPMetalloproteinase 4a752782853Boiga irregularisYLVVDNR53LYLVVDNR72496.282−0.4
AYEMVNILNVIFR81AYEMVNLLNTMFR61801.3932−1.5
YLVVDNR53TNTARCLTLYLKLKYQPPAYLVVDNR51622.94474
SPDYGMVEPGTK64DYPGHHACALEDDQLLMCSEDYGMVEPASK51850.118211.6
SVMPMetalloproteinase172046653Naja mossambicaTSMVAITMAHQMGHNLGMNDDR122NTDMVAGTMAHQMGHNLGLNDHR52841.04857.2
AHQMGHNLGMNDDR81EFVRQAFANALYPPWWMPKYMVPQKCKAHEMGHNLGLNDHR50986.093912.5
SVMPMetalloproteinase (type III) 2a 752782867Boiga irregularisTEGMVEPGTK 75QTEGMVEPGTK60580.26880.3
TESFVASTMAHELGHNLGINHDR136TNPFVGSTMAHEMGHNLGLNHDR60634.55047.2
SVMPAtrase B, partial289655973Naja atraFGEWRETVLLPR 92AFGEWRETVLLPR63525.28780.2
SVMPMetalloproteinase 11a 699236668Hypsiglena sp. JMG-2014QPCQNNQGYCYNGK94GPGLQPCGNNQGYCYDGK44964.4064−0.3
KYIEFYIVVDH78SMTKPCLMRVCDLVKYLQFYLVVDHVKSYVLCYR48823.632211.2

Share and Cite

MDPI and ACS Style

Chapeaurouge, A.; Silva, A.; Carvalho, P.; McCleary, R.J.R.; Modahl, C.M.; Perales, J.; Kini, R.M.; Mackessy, S.P. Proteomic Deep Mining the Venom of the Red-Headed Krait, Bungarus flaviceps. Toxins 2018, 10, 373. https://0-doi-org.brum.beds.ac.uk/10.3390/toxins10090373

AMA Style

Chapeaurouge A, Silva A, Carvalho P, McCleary RJR, Modahl CM, Perales J, Kini RM, Mackessy SP. Proteomic Deep Mining the Venom of the Red-Headed Krait, Bungarus flaviceps. Toxins. 2018; 10(9):373. https://0-doi-org.brum.beds.ac.uk/10.3390/toxins10090373

Chicago/Turabian Style

Chapeaurouge, Alex, Andreza Silva, Paulo Carvalho, Ryan J. R. McCleary, Cassandra Marie Modahl, Jonas Perales, R. Manjunatha Kini, and Stephen P. Mackessy. 2018. "Proteomic Deep Mining the Venom of the Red-Headed Krait, Bungarus flaviceps" Toxins 10, no. 9: 373. https://0-doi-org.brum.beds.ac.uk/10.3390/toxins10090373

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