# Serine Protease Variants Encoded byEchis ocellatus Venom Gland cDNA: Cloning and Sequencing Analysis

Authors: S. S. Hasson; R. A. Mothana; T. A. Sallam; M. S. Al-balushi; M. T. Rahman; A. A. Al-Jabri
Journal: Journal of Biomedicine and Biotechnology
 (2010)

---

## Abstract

Envenoming byEchis saw-scaled viper is the leading cause of death and morbidity in Africa due to snake bite. Despite its medical importance, there have been few investigations into the toxin composition of the venom of this viper. Here, we report the cloning of cDNA sequences encoding four groups or isoforms of the haemostasis-disruptive Serine protease proteins (SPs) from the venom glands of Echis ocellatus. All these SP sequences encoded the cysteine residues scaffold that form the 6-disulphide bonds responsible for the characteristic tertiary structure of venom serine proteases. All the Echis ocellatus EoSP groups showed varying degrees of sequence similarity to published viper venom SPs. However, these groups also showed marked intercluster sequence conservation across them which were significantly different from that of previously published viper SPs. Because viper venom SPs exhibit a high degree of sequence similarity and yet exert profoundly different effects on the mammalian haemostatic system, no attempt was made to assign functionality to the new Echis ocellatus EoSPs on the basis of sequence alone. The extraordinary level of interspecific and intergeneric sequence conservation exhibited by the Echis ocellatus EoSPs and analogous serine proteases from other viper species leads us to speculate that antibodies to representative molecules should neutralise (that we will exploit, by epidermal DNA immunization) the biological function of this important group of venom toxins in vipers that are distributed throughout Africa, the Middle East, and the Indian subcontinent.


---

## Body

## 1. Introduction

Envenoming resulting from snake bites is an important public health hazard in many regions, particularly in tropical and subtropical countries [1, 2]. The saw-scaled viper Echis ocellatus is the most abundant [3] and medically important viper species in West Africa [4].Envenoming by saw-scaled viper (Viperidae: Echis) species is thought to be responsible for more snakebite deaths worldwide than any other snake genus [5]. In northern Nigeria, E. ocellatus is responsible for 95% of all envenoming by snakes [6], causing several hundred deaths annually [7]. The precise incidence of snakebite is difficult to determine and is often grossly underestimated, but in some areas of the Nigerian savannas, victims of E. ocellatus envenoming may occupy more than 10% of hospital beds [8]. In the Benue valley of Nigeria, for example, the estimated incidence is 497 per 100 000 population per year with 10%–20% untreated mortality [9]. Local effects of Echis viper envenoming include pain, swelling, blistering, and haemorrhage which, in severe cases, can lead to necrosis, permanent disfigurement, and even amputation of the affected limb [10]. Systemic effects include potentially lethal consumption coagulopathy, haemorrhage and hypovolaemic shock [10].Snake venoms contain a great variety of toxic proteases [11, 12]. Many of these components are proteases, for example, metalloproteases [13], serine proteases [14], phospholipases A2 [15] and C-type lectins [16] and mediate their toxicity by either stimulating or inhibiting the haemostatic system of human victims or experimental animals, resulting in clinical complications of blood clotting or uncontrolled haemorrhage [12, 17–19]. Several of these proteinases cleave plasma proteins of the victims in a specific manner with varying degrees of substrate specificity. Thus, while some serine proteases have both fibrinogenolytic and fibrinolytic activities, others have only fibrinogenolytic activity and are called “thrombin-like” proteases [19–25]. Approximately 100 snake venom toxins have been identified as “thrombin-like” enzymes activating the blood coagulation factor [26]. These “thrombin-like” proteases hydrolyze fibrinogen specifically and release either fibrinopeptide A or B or both [27] resulting in the disruption of the blood coagulation system by producing abnormal fibrin clots composed of short polymers that are rapidly dispersed and no longer cross-linked by activated factor XIII [28].Another group of serine proteases ofBatroxobin, Crotalase, and Ancrod venoms affect other substrates, for example, plasminogen [27] by cleaving fibrinogen in manner distinct from that of thrombin. Other venom serine proteases function like mammalian kallikrein (or kininogenase) releasing bradykinin from kininogen [29–31] and are called “kallikrein-like” proteases [29], an example of this is halystase [32], a kallikrein-like serine protease isolated from A. halys blomhoffii venom, which cleaves the β chain at Arg42 and slowly degrades the α chain of fibrinogen to generate a product that is no longer converted to normal fibrin clots by thrombin; this results in both reduction of blood pressure as well as inhibiting fibrinogen clotting in the victims. Another kallikrein-like serine protease with potent biological activity but with different physicochemical properties from those of halystase has been isolated from the venoms of  A. caliginosus, C. atrox, C. viridis, and Trimeresurus mucrosquamatus[29, 30, 32–34]. The latter showed both a strong β-fibrinogenolytic and kallikrein-like activities, cleaving β-chain of fibrinogen molecules specifically and releasing bradykinin from kininogen, respectively. Moreover, the purified enzymes indicated that they have specificities different from thrombin and thrombin-like proteases of snake venom reported previously by decreasing fibrinogen levels in plasma and prolonging bleeding without formation of fibrin clots. They also exhibit amidase activity against N-benzoyl-Pro-Phe-Arg-p-nitroanilide, which is a specific synthetic substrate for kallikrein-like proteases.In addition, there have been a few reports on venom serine proteases with a unique activity, such as ACC-C, a protein C activator isolated from theA. contortrix venom [35] (which inhibits blood coagulation by inactivating the activated forms of factor V and VIII), a plasminogen activator such as TSV-PA isolated from the T. stejnegeri venom [36, 37], PA-BJ, a platelet aggregating enzyme isolated from the B. jararaca and Trimeresurus mucrosquamatusvenoms [38],  and  RVV-V, a factor V-activating enzyme isolated from the V. russelli venom [39].These data indicate that snake venom serine proteases comprise an enzyme superfamily with multifunctional activities that may have diverged or have undergone gene duplication resulting in alteration of their biological properties during the process of evolution thus acquiring special functions [40, 41]. Although a considerable amount of data is now available, no standardised grouping of these venom serine proteases has yet been documented. However, in 2001 Wang et al. [27] compared sequences of 40 serine proteinases isolated from different snake venoms, using a constructed phylogram in which such sequences were clustered into three groups designated as coagulating enzymes, kininogenases, and plasminogen activators.No Serine proteinases have yet been purified from venom of the West African saw-scaled viperEchis ocellatus, in particular or for members of the Echis genus in general. However, the fact that the serine protease superfamily was important in the venom of the Viperidae suggested that such enzymes should be present in the venom of E. ocellatus and that serine protease-specific antibodies are likely to be an important factor in E. ocellatus envenoming. We therefore screened the E. ocellatus cDNA library in order to isolate and characterise different isoforms or variants of this enzyme superfamily.


## 2. Materials and Methods

### 2.1. Animals

AdultE. ocellatus(Nigeria) carpet viper used in this study was maintained in the herpetarium, Liverpool School of Tropical Medicine, Liverpool, UK.


### 2.2. Extraction of Total Venom Gland RNA and Construction of cDNA Libraries

Venom glands were dissected from threeEchis ocellatussnakes. The vipers were sacrificed 3 days after venom extraction when toxin gene transcription rates are at a peak. Glands were homogenized under liquid Nitrogen and total RNA extracted using guanidinium thiocyanate-phenol-chloroform as described previously [15]. Lambda phage cDNA libraries for E. ocellatuswere constructed by RT-PCR using the SMART cDNA library construction kit (Clontech, California, USA). The lambda phage of the E. ocellatuswas packaged using Gigapack III Gold Packaging Extract (Stratagene) and boiled for 5 min prior to being used as targets of polymerase chain reaction (PCR) amplification.


### 2.3. Isolation and Analysis of cDNA Sequences

A PCR strategy [42] was used to isolate sequences encoding serine proteinases from the cDNA libraries. A sense primer (5V-GGA-TCC-ATG-GTG-CTG-ATC-AGA-GTG-CTA-ATC-GCA-3V) and an antisense primer (5V-CTC-GAG-TGG-GGG-GCA-AGT-CGC-AGT-TGT-ATT-TCC-3V) complimentary to highly conserved amino-terminal signal peptide (M-V-L-I-R-V) and to the less conserved carboxy-terminal (T-T-A-T-C-P-P) domains of published serine proteinases DNA sequences of related viper species were synthesized commercially (Sigma-Genosys, UK). A TAG stop codon was inserted in the 3′ primer and BamH1 and Xho1 restriction endonuclease sites (bold) were included in the 5′  and 3′  primers, respectively, to facilitate future subcloning into mammalian expression plasmids. PCR was performed using an initial denaturation (95°C—6 minutes) and annealing (55°C—1 minute) step, followed by 35 cycles (1 minute each) of extension (74°C), denaturation (94°C), and annealing (55°C), and a terminal extension step (7 minutes) at 72°C in a thermal cycler (Gene Cycler, BioRad, Hercules, CA, USA). The inclusion of water-only controls with each PCR reaction allowed us to monitor and prevent cross-over contamination. The amplicons were subcloned into the TA cloning vector, pCR 2.1-TOPO, (Invitrogen, Groningen, The Netherlands) and used to transform chemically competent E. colicells (TOP10F’, Invitrogen) under ampicillin selection. Plasmid DNA was extracted (Mini-spin prep kit, Qiagen, Hilden, Germany) and digested with BamH1 and Xho1 at 37°C to select plasmids containing inserts of the predicted size for DNA sequencing. DNA sequencing was carried out by the dideoxy-nucleotide chain-termination method in a Beckman Coulter CEQk 2000 XL DNA Analysis System. To confirm that the cDNA sequences encoded CTLs, the predicted amino acid sequences were subjected to BLAST searches of the GenBank, PDB, SwissProt, PIR, and PRF databases. All the cDNAs exhibited significant sequence homology to Serine protienases of related vipers. The CLUSTALW program [43] with PAM 250 residue weight matrix was used to align deduced amino acid sequences representing each E. ocellatus Serine protienases isoforms with analogues in venoms from related Viperidae species as illustrated in Table 1. Serine proteinase (CAB62591) from V. lebetina [44], Serine protease 1 (AAR24534) from B. gabonica [45], Thrombin-like enzyme pre. (AAK12273) from D. acutus[46], Venom serine protease 5 (AAN52350) from T. stejnegeri[47], Serine proteinase 3 pre. (O13063) from gramineus[48], Serine proteinase A Precursor (Q9PTU8) from B. jararaca[46, 49], Serine proteinase 2A pre. (O13060) from T. gramineus[45, 48], Serine protease (AAP42416) from B. jararacussu[50], KN-BJ2 (BAA20283) from B. jararaca[51], Serine proteinase 1 pre. (AAG10788) from T. jerdonii[52], Thrombin-like serine protease (AAL68708) from G. ussuriensis[53], and, finally, Serine protease catroxase I pre. (AAL77226) from C. atrox[54]. The phylogenetic trees constructed from the above alignments were generated by a neighbour-joining [55] algorithm in Lasergene software (DNASTAR, USA). The predicted antigenic profile [56] of the published and new Echis ocellatusserine protease (EoSer) isoforms analysed here was determined using Protean Software (DNASTAR).Table 1
Percent sequence similarity betweenE. ocellatus serine proteases and analogous molecules from related viper species.
SpeciesAccession no.ReferencesSerine proteaseEoSP-1EoSP-3EoSP-7EoSP-17V. lebetinaCAB62591Siigur et al. [44]Serine proteinase65668064B. gabonicaAAR24534Francischetti et al. [45]Serine protease 162638265D. acutusAAK12273Liang et al. [46]Thrombin-like enzyme pre.67677170T. stejnegeriAAN52350Lee and Zhang [47]Venom serine protease 566666569T. gramineusO13063Deshimaru et al. [48]Serine proteinase 3 pre.71716176B. jararacaQ9PTU8Murayama, [49]Serine proteinase A pre.66686674T. gramineusO13060Deshimaru et al. [48]Serine proteinase 2A pre.65657370B. jararacussuAAP42416Kashima et al. [50]Serine protease63637268B. jararacaBAA20283Serrano et al. [51]KN-BJ269686269T. jerdoniiAAG10788Lu et al. [52]Serine proteinase 1 pre.65657170G. ussuriensisAAL68708Zhao et al. [53]Thrombin-like serine protease67677271C. atroxAAL77226Tsai et al. [54]Serine protease catroxase I pre.66656365


## 2.1. Animals

AdultE. ocellatus(Nigeria) carpet viper used in this study was maintained in the herpetarium, Liverpool School of Tropical Medicine, Liverpool, UK.


## 2.2. Extraction of Total Venom Gland RNA and Construction of cDNA Libraries

Venom glands were dissected from threeEchis ocellatussnakes. The vipers were sacrificed 3 days after venom extraction when toxin gene transcription rates are at a peak. Glands were homogenized under liquid Nitrogen and total RNA extracted using guanidinium thiocyanate-phenol-chloroform as described previously [15]. Lambda phage cDNA libraries for E. ocellatuswere constructed by RT-PCR using the SMART cDNA library construction kit (Clontech, California, USA). The lambda phage of the E. ocellatuswas packaged using Gigapack III Gold Packaging Extract (Stratagene) and boiled for 5 min prior to being used as targets of polymerase chain reaction (PCR) amplification.


## 2.3. Isolation and Analysis of cDNA Sequences

A PCR strategy [42] was used to isolate sequences encoding serine proteinases from the cDNA libraries. A sense primer (5V-GGA-TCC-ATG-GTG-CTG-ATC-AGA-GTG-CTA-ATC-GCA-3V) and an antisense primer (5V-CTC-GAG-TGG-GGG-GCA-AGT-CGC-AGT-TGT-ATT-TCC-3V) complimentary to highly conserved amino-terminal signal peptide (M-V-L-I-R-V) and to the less conserved carboxy-terminal (T-T-A-T-C-P-P) domains of published serine proteinases DNA sequences of related viper species were synthesized commercially (Sigma-Genosys, UK). A TAG stop codon was inserted in the 3′ primer and BamH1 and Xho1 restriction endonuclease sites (bold) were included in the 5′  and 3′  primers, respectively, to facilitate future subcloning into mammalian expression plasmids. PCR was performed using an initial denaturation (95°C—6 minutes) and annealing (55°C—1 minute) step, followed by 35 cycles (1 minute each) of extension (74°C), denaturation (94°C), and annealing (55°C), and a terminal extension step (7 minutes) at 72°C in a thermal cycler (Gene Cycler, BioRad, Hercules, CA, USA). The inclusion of water-only controls with each PCR reaction allowed us to monitor and prevent cross-over contamination. The amplicons were subcloned into the TA cloning vector, pCR 2.1-TOPO, (Invitrogen, Groningen, The Netherlands) and used to transform chemically competent E. colicells (TOP10F’, Invitrogen) under ampicillin selection. Plasmid DNA was extracted (Mini-spin prep kit, Qiagen, Hilden, Germany) and digested with BamH1 and Xho1 at 37°C to select plasmids containing inserts of the predicted size for DNA sequencing. DNA sequencing was carried out by the dideoxy-nucleotide chain-termination method in a Beckman Coulter CEQk 2000 XL DNA Analysis System. To confirm that the cDNA sequences encoded CTLs, the predicted amino acid sequences were subjected to BLAST searches of the GenBank, PDB, SwissProt, PIR, and PRF databases. All the cDNAs exhibited significant sequence homology to Serine protienases of related vipers. The CLUSTALW program [43] with PAM 250 residue weight matrix was used to align deduced amino acid sequences representing each E. ocellatus Serine protienases isoforms with analogues in venoms from related Viperidae species as illustrated in Table 1. Serine proteinase (CAB62591) from V. lebetina [44], Serine protease 1 (AAR24534) from B. gabonica [45], Thrombin-like enzyme pre. (AAK12273) from D. acutus[46], Venom serine protease 5 (AAN52350) from T. stejnegeri[47], Serine proteinase 3 pre. (O13063) from gramineus[48], Serine proteinase A Precursor (Q9PTU8) from B. jararaca[46, 49], Serine proteinase 2A pre. (O13060) from T. gramineus[45, 48], Serine protease (AAP42416) from B. jararacussu[50], KN-BJ2 (BAA20283) from B. jararaca[51], Serine proteinase 1 pre. (AAG10788) from T. jerdonii[52], Thrombin-like serine protease (AAL68708) from G. ussuriensis[53], and, finally, Serine protease catroxase I pre. (AAL77226) from C. atrox[54]. The phylogenetic trees constructed from the above alignments were generated by a neighbour-joining [55] algorithm in Lasergene software (DNASTAR, USA). The predicted antigenic profile [56] of the published and new Echis ocellatusserine protease (EoSer) isoforms analysed here was determined using Protean Software (DNASTAR).Table 1
Percent sequence similarity betweenE. ocellatus serine proteases and analogous molecules from related viper species.
SpeciesAccession no.ReferencesSerine proteaseEoSP-1EoSP-3EoSP-7EoSP-17V. lebetinaCAB62591Siigur et al. [44]Serine proteinase65668064B. gabonicaAAR24534Francischetti et al. [45]Serine protease 162638265D. acutusAAK12273Liang et al. [46]Thrombin-like enzyme pre.67677170T. stejnegeriAAN52350Lee and Zhang [47]Venom serine protease 566666569T. gramineusO13063Deshimaru et al. [48]Serine proteinase 3 pre.71716176B. jararacaQ9PTU8Murayama, [49]Serine proteinase A pre.66686674T. gramineusO13060Deshimaru et al. [48]Serine proteinase 2A pre.65657370B. jararacussuAAP42416Kashima et al. [50]Serine protease63637268B. jararacaBAA20283Serrano et al. [51]KN-BJ269686269T. jerdoniiAAG10788Lu et al. [52]Serine proteinase 1 pre.65657170G. ussuriensisAAL68708Zhao et al. [53]Thrombin-like serine protease67677271C. atroxAAL77226Tsai et al. [54]Serine protease catroxase I pre.66656365


## 3. Results

### 3.1. Isolation of cDNAs EncodingE. Ocellatus Serine Protease

PCR screening of theEchis ocellatusvenom gland cDNA libraries resulted in a total of 14 E. ocellatus(Eo) cDNAs whose sequences matched (BLAST searches) those of published Serine proteases. The cDNAs consisted of 822 nucleotides (Figures 1 and 2(a)) and were predicted to encode an open reading frame proteins of 264 amino acids (28.5 kDa) (Figure 2(b)). Alignment of the predicted amino acid sequences of the 14 specific cDNAs encoding the EoSPproteins (Figure 2(b)) revealed sequence variations. The sequence similarity between the EoSPvariants proteins was less than 60% for the mature protein-coding region but over 90% for regions coding both the signal peptide and the carboxyl-terminal end. Where two or more identical sequences were obtained from any one of these libraries, a single representative cDNA was used for subsequent analysis. Structural properties analysis (Emin algorithm-DNASTAR, USA) (Figure 3) was used to categorise the 14 Serine protease sequences into four distinct groups, based solely on sequence alignment.Figure 1
PCR product of theE. ocellatus serine proteases. Analysis of PCR amplification products by 0.7% agarose gel electrophoresis. Bands were visualised using the ultraviolet transillumination. Lane 2: represents the amplified PCR product (circled) of about 800 bp from E. ocellatus venom glands cDNA compared with Lane 1:1 kb ladder DNA-marker bands, of known molecular weight. Lanes 3 and 4 represent a H2O negative control and a SOD positive control, respectively.(a) The nucleotide sequence of the fourteenE. ocellatus venom gland cDNAs resulting from PCR amplification. (b) Deduced amino acid sequences of E. ocellatus venom gland serine protease cDNAs.
(a)(b)Figure 3
Differentiation of the fourteen cDNA-encodingE. ocellatus venom gland serine proteases.The predicted surface probabilities (Emin algorithm, DNASTAR, USA) of the 14 E. ocellatusserine protease cDNAs were aligned. The boxed areas indicate group specific structural motifs.


### 3.2. BLAST Search of the Predicted Amino Acid Sequence

Accession numbers assigned to the newEchis ocellatusSerine protease sequences are as follows: “group 1” EoSer-1 (GU562413), “group 2” EoSer-3 (GU592440), “group 3” EoSer-17 (GU592441), and “group 4” EoSer-7 (GU592439). The predicted amino acid sequences of the EoSP-01, EoSP-03, EoSP-07,and EoSP-17were submitted to BLAST searches of the genetic data bases and their similarity to published viper serine protease (Table 1) confirmed that the EoSPcDNAs encoded serine proteases.


### 3.3. Comparison ofE. Ocellatus cDNAs with Analogous Serine Proteases from Other Viper Species

All theEoSer-variants contained the serine protease-consensus 24 amino acid signal peptide sequence (Figure 4, arrows), including the six-amino acids-activated motif. The signal peptide residues were followed by a protease domain of 236 residues. The deduced primary structures of all EoSPcDNA clones include the requisite, highly conserved, 12 cysteine residues that form the 6-disulphide bonds responsible for the characteristic tertiary structure of venom serine proteases. The complete amino acid sequences of the EoSPvariants were aligned with those of other venom serine proteases (Figure 4). Viper venom SP sequences in the genetic databases were compared with the E. o groups (Table 2 and Figure 4) by BLAST. Groups 1–4 represent novel, highly similar, SP isomers with less than 65% sequence similarity to analogues in related viper species. Group 4 showed the greatest sequence similarity (80% and 82%) to the Serine protease of the African V. lebetina and B. gabonicavipers, respectively. Of all the EoSPclusters seemed to represent a SP sequence which showed the highest sequence similarity range between 62% to 70% to the SP of the vipers. None of the clusters showed more than 72% sequence similarity to the partial peptide sequences for the Thrombin-like serine protease isolated from the venom of the G. ussuriensisviper [52]. Similarly, the Serine protease catroxase I pre. of C. atroxvenom showed no greater than 65% sequence similarity to any of the EoSPsequences.Table 2
Comparison of amino acid motifs which are responsible for the potent effects and characterisation of some published venom serine proteases with the fourEoSP cDNAs.
Amino acidTSV-PABatroxobinAncrodEoSP-1EoSP-17EoSP-3EoSP-7ReferencesH/RH57H57H57H67H67H67R67Braud et al. [57]DD102D102D102D112D112D112D112SS195S195S195S208S208S208T208HH192G192N192K205L205K205K205FF193G193S193G206G206G206A206Guinto et al. [58]DD189D189D189G202D202G202D202PP225P225P225P235P235P235P235PP219P219P219P228V228P228P228Braud et al. [57]DD96N96R96Y106Y106Y106Y106Lee and Zhang [47]DD97V97T97T107T107T107T107EE98I98S98L108L108K108R108HDS: Catalytic Traid; H/F: substrate specificity; D & P: Architecture of water channel; P: Evolutionary region to kallikrein; DDE: substrate specificity to plasminogen.Figure 4
Amino acid sequence similarity betweenEoSPVariants and serine proteases from related vipers. The residues shaded in black correspond to residues that are identical to EoSP-01.The asteriks [*] represented the tweleve conserved cysteine residues. The catalytic traid His/Arg(67), Asp(110) and Ser(208) are represented in red circules. Activated peptide where the mature proteine cleaved is represented by green rectangle.


### 3.4. Predicted Antigenic Profile Analysis ofE. ocellatus Serine Proteases with Analogous Molecules

Since the main focus of our research is to develop toxin neutralising antibodies by immunisation with DNA encoding specific toxins in venoms of the most medically important African vipers [15, 59, 60], we next compared the algorithm-predicted immunogenicity of the E. ocellatusserine protease cluster cDNA sequences with those of all the published SPs from vipers of African origin (Figure 5). The predicted antigenic profiles of the published and new E. ocellatusserine proteases were analysed as shown in Figure 5 using Protean Software (DNASTAR, USA) [53]. The deduced signal peptide domains of the EoSPvariants are separated by a vertical dotted line, as these would normally be cleaved from the native proteins during posttranslational. The thin vertical boxes depict the residues comprising the catalytic traid, H/R/N, D/G/N, and S/P/N/T (67, 110, and 208), that show the greatest immunogenic domains conservation common to all the new and published African viper venom SPs sequences as demonstrated in Figure 5.Figure 5
Comparison of antigenic profile of theEoSer variants with analogous serine proteases used in Figure 4. The top horizontal scale represents the number of amino acid residues. The conserved signal peptide is separated from the mature protein by a vertical dotted line. The three vertical boxes were drawn to indicate the conserved catalytic traid regions described in the text.


## 3.1. Isolation of cDNAs EncodingE. Ocellatus Serine Protease

PCR screening of theEchis ocellatusvenom gland cDNA libraries resulted in a total of 14 E. ocellatus(Eo) cDNAs whose sequences matched (BLAST searches) those of published Serine proteases. The cDNAs consisted of 822 nucleotides (Figures 1 and 2(a)) and were predicted to encode an open reading frame proteins of 264 amino acids (28.5 kDa) (Figure 2(b)). Alignment of the predicted amino acid sequences of the 14 specific cDNAs encoding the EoSPproteins (Figure 2(b)) revealed sequence variations. The sequence similarity between the EoSPvariants proteins was less than 60% for the mature protein-coding region but over 90% for regions coding both the signal peptide and the carboxyl-terminal end. Where two or more identical sequences were obtained from any one of these libraries, a single representative cDNA was used for subsequent analysis. Structural properties analysis (Emin algorithm-DNASTAR, USA) (Figure 3) was used to categorise the 14 Serine protease sequences into four distinct groups, based solely on sequence alignment.Figure 1
PCR product of theE. ocellatus serine proteases. Analysis of PCR amplification products by 0.7% agarose gel electrophoresis. Bands were visualised using the ultraviolet transillumination. Lane 2: represents the amplified PCR product (circled) of about 800 bp from E. ocellatus venom glands cDNA compared with Lane 1:1 kb ladder DNA-marker bands, of known molecular weight. Lanes 3 and 4 represent a H2O negative control and a SOD positive control, respectively.(a) The nucleotide sequence of the fourteenE. ocellatus venom gland cDNAs resulting from PCR amplification. (b) Deduced amino acid sequences of E. ocellatus venom gland serine protease cDNAs.
(a)(b)Figure 3
Differentiation of the fourteen cDNA-encodingE. ocellatus venom gland serine proteases.The predicted surface probabilities (Emin algorithm, DNASTAR, USA) of the 14 E. ocellatusserine protease cDNAs were aligned. The boxed areas indicate group specific structural motifs.


## 3.2. BLAST Search of the Predicted Amino Acid Sequence

Accession numbers assigned to the newEchis ocellatusSerine protease sequences are as follows: “group 1” EoSer-1 (GU562413), “group 2” EoSer-3 (GU592440), “group 3” EoSer-17 (GU592441), and “group 4” EoSer-7 (GU592439). The predicted amino acid sequences of the EoSP-01, EoSP-03, EoSP-07,and EoSP-17were submitted to BLAST searches of the genetic data bases and their similarity to published viper serine protease (Table 1) confirmed that the EoSPcDNAs encoded serine proteases.


## 3.3. Comparison ofE. Ocellatus cDNAs with Analogous Serine Proteases from Other Viper Species

All theEoSer-variants contained the serine protease-consensus 24 amino acid signal peptide sequence (Figure 4, arrows), including the six-amino acids-activated motif. The signal peptide residues were followed by a protease domain of 236 residues. The deduced primary structures of all EoSPcDNA clones include the requisite, highly conserved, 12 cysteine residues that form the 6-disulphide bonds responsible for the characteristic tertiary structure of venom serine proteases. The complete amino acid sequences of the EoSPvariants were aligned with those of other venom serine proteases (Figure 4). Viper venom SP sequences in the genetic databases were compared with the E. o groups (Table 2 and Figure 4) by BLAST. Groups 1–4 represent novel, highly similar, SP isomers with less than 65% sequence similarity to analogues in related viper species. Group 4 showed the greatest sequence similarity (80% and 82%) to the Serine protease of the African V. lebetina and B. gabonicavipers, respectively. Of all the EoSPclusters seemed to represent a SP sequence which showed the highest sequence similarity range between 62% to 70% to the SP of the vipers. None of the clusters showed more than 72% sequence similarity to the partial peptide sequences for the Thrombin-like serine protease isolated from the venom of the G. ussuriensisviper [52]. Similarly, the Serine protease catroxase I pre. of C. atroxvenom showed no greater than 65% sequence similarity to any of the EoSPsequences.Table 2
Comparison of amino acid motifs which are responsible for the potent effects and characterisation of some published venom serine proteases with the fourEoSP cDNAs.
Amino acidTSV-PABatroxobinAncrodEoSP-1EoSP-17EoSP-3EoSP-7ReferencesH/RH57H57H57H67H67H67R67Braud et al. [57]DD102D102D102D112D112D112D112SS195S195S195S208S208S208T208HH192G192N192K205L205K205K205FF193G193S193G206G206G206A206Guinto et al. [58]DD189D189D189G202D202G202D202PP225P225P225P235P235P235P235PP219P219P219P228V228P228P228Braud et al. [57]DD96N96R96Y106Y106Y106Y106Lee and Zhang [47]DD97V97T97T107T107T107T107EE98I98S98L108L108K108R108HDS: Catalytic Traid; H/F: substrate specificity; D & P: Architecture of water channel; P: Evolutionary region to kallikrein; DDE: substrate specificity to plasminogen.Figure 4
Amino acid sequence similarity betweenEoSPVariants and serine proteases from related vipers. The residues shaded in black correspond to residues that are identical to EoSP-01.The asteriks [*] represented the tweleve conserved cysteine residues. The catalytic traid His/Arg(67), Asp(110) and Ser(208) are represented in red circules. Activated peptide where the mature proteine cleaved is represented by green rectangle.


## 3.4. Predicted Antigenic Profile Analysis ofE. ocellatus Serine Proteases with Analogous Molecules

Since the main focus of our research is to develop toxin neutralising antibodies by immunisation with DNA encoding specific toxins in venoms of the most medically important African vipers [15, 59, 60], we next compared the algorithm-predicted immunogenicity of the E. ocellatusserine protease cluster cDNA sequences with those of all the published SPs from vipers of African origin (Figure 5). The predicted antigenic profiles of the published and new E. ocellatusserine proteases were analysed as shown in Figure 5 using Protean Software (DNASTAR, USA) [53]. The deduced signal peptide domains of the EoSPvariants are separated by a vertical dotted line, as these would normally be cleaved from the native proteins during posttranslational. The thin vertical boxes depict the residues comprising the catalytic traid, H/R/N, D/G/N, and S/P/N/T (67, 110, and 208), that show the greatest immunogenic domains conservation common to all the new and published African viper venom SPs sequences as demonstrated in Figure 5.Figure 5
Comparison of antigenic profile of theEoSer variants with analogous serine proteases used in Figure 4. The top horizontal scale represents the number of amino acid residues. The conserved signal peptide is separated from the mature protein by a vertical dotted line. The three vertical boxes were drawn to indicate the conserved catalytic traid regions described in the text.


## 4. Discussion

Serine proteases are a major component of viper venoms and are thought to disrupt several distinct elements of the blood coagulation system of envenomed victims. A detailed understanding of the functions of these enzymes is important for both acquiring a full understanding of the pathology of envenoming and because these venom proteins have shown a vital role in treating blood coagulation disorders.In general, serine proteinases including fibrinogenolytic enzymes are very abundant inViperidae venoms in which they may account for 20% of their total protein content [61].The unique specificity of snake venom proteinases makes them potentially useful in research of fibrinogen-depletion and limited proteolysis [62, 63]. This may be due to the existence of multiple forms of serine proteases in the venom of a single viper species which is likely to contribute to the diverse biological effects exerted by the whole venom.Therefore, screening the E. ocellatuscDNA library to isolate different isoforms or variants of serine proteases was the aim of this research work. The results obtained in this work provide the first molecular sequence data for E. ocellatusserine proteases they also reveal that the serine protease composition of E. ocellatusis as complex as that of the better characterised Viperidaespecies. The utilization of PCR amplification of E. ocellatusvenom gland cDNA with the new viper serine protease-specific primers was successful and produced fourteen cDNAs sequences that were identified (BLAST) as belonging to the serine protease enzyme family. All EoSPcDNAs were of similar total length (approximately 0.80 kb, Figure 1) and encoded 260 amino acids (Figure 2(b)) with a predicted molecular weight of 28.5 kDa. To differentiate between the isolated EoSPclones a surface probability algorithm was used to assign the 14 E. ocellatusserine protease cDNAs into four main groups (Figure 3). A single representative clone from each group was chosen for further analyses as described earlier. The sequence similarity between the EoSPvariants proteins was less than 60% for the mature protein-coding region but over 90% for regions coding both the signal peptide and the carboxyl-terminal end. Thus the latter two regions are highly conserved, which explains why the PCR experiment to amplify the cDNAs-encoding EoSPclones was successful.TheEoSPcDNA sequences were confirmed by BLAST searches as encoding serine proteases (Figure 4). The greatest sequence similarity was between EoSer-7 and B. gabonicaand V. labetina(80% and 85%) with the remaining EoSPcDNAs showing 60%–76% sequence similarity with other snake venom serine proteinases as illustrated in Table 1. From the proteins with known biological activity, sequence similarities of the EoSPvariants (i.e., EoSer-01, EoSer-03, EoSer-07and EoSer-17) were 62%–69% with the kinin-releasing and fibrinogen-clotting serine protease (KN-BJ) from venom of B. jararaca[51] (Table 1). The putative 18 amino acid signal-peptide of the EoSPvariants was as conserved (over 90% sequence similarity) as that in the serine proteases of other viper species (Figure 4, arrows). Following the signal peptide all the EoSPvariants contained the predicted six-amino acid cleavage (activation) site Q-K/T/M/E-S-S-E-L/P (Figure 4 in green) as proposed for batroxobin [64]; thus cleavage generates a hydrophilic zymogen peptide, based on the processing site of pre-peptides of mammalian serine proteinases [65–67]. Comparison of the EoSPvariants with analogous members of the serine protease family revealed that all EoSPvariants encoded the presumed catalytic triad, which is common to venom serine proteases H67, D110 and S208 as shown in Figure 4. Such residues were highly conserved in groups 1–3, except proteins of group 4 (Figures 2(b) and 4) which contain R instead of H at the same position (Figure 4). Furthermore, comparison of the EoSPamino acid sequence alignment with analogous venom serine proteases (Figure 4) revealed a conserved consensus active site of L-T/S-A-A-H/R/N-C corresponding to position 63–68, as previously determined [68]. Most SVSPs are likely to be glycoproteins showing a variable number of N- or O-glycosylation sites in sequence positions that differ from one SVSP to the other [69]. Using the primary structure of EoSPvariants (Figure 4) the putative N-linked glycosylation sites, Asn-X-Thr/Ser [45], were found and are located at two different positions. EoSer-01, EoSer-03, and EoSer-17[N44-X45-S46  andN257-X258-T259] and EoSer-07 [N124-R125-T126  and N257-T258-T258]. Although such motifs are thought to be needed for protein stabilization rather than for the catalytic function of the venom enzymes [30], confirmation of the roles of such motifs in venom proteases remain to be investigated. All serine proteases have a common pattern of 6-disulfide bridges [69, 70]. They contain twelve cysteine residues, ten of which form five disulfide bonds, based on the homology with trypsin [64]; the remaining two cysteines form a unique and conserved bridge among SVSPs, involving Cys245e (chymotrypsinogen numbering), found in the C-terminal extension [35].From the results obtained this was found in allEoSPclones (Figure 2(b)) that encoded the common 12 cysteine residues in which are strongly conserved forming putative disulphide bridges which are located atCys31 Cys52, C68, C100, C145, C165, C176, C204, C214, C229,andC260 (Figure 4). This suggests that the EoSPproteins possess a similar tertiary structure to that of other serine proteases which are well characterized.Despite such sequence and structural conservation, viper venom serine proteases show very divergent effects on haemostasis as previously stated. In some cases certain amino acid sequences have been shown to be responsible for such effects as demonstrated in Table2. Although such table gives a preliminary prediction of the functional characterization of the EoSPcDNAs in comparison with well-known characterized venom serine proteases, it cannot be considered as a functional confirmation or even a categorization strategy to differentiate between the four EoSPcDNAs. However, from Figure 4 and Table 2 it can be generally concluded that such comparison demonstrates that the enzymes encoded by the four EoSPcDNAs confer multiple haemostasis-disruptive activities to E. ocellatusvenom. Furthermore, the sequence and predicted structural similarities of these four EoSPgroups suggest that an antibody generated to one group may be capable of neutralizing the other group of EoSPs. To examine this permeability the sequences of EoSPgroups were subjected to a more specific algorithm that predicted amino acid motifs of high immunogenicity. A protein structure-predicting algorithm [56] has been used (i) to identify domains of strong antigenic potential in the toxin gene product and (ii) to determine whether these domains are conserved in analogous venom toxin gene products of related vipers. The signal peptide was separated from the mature protein by dotted line as would be cleaved posttranslationally. The peaks shown by the EoSPs profile indicate the numerous domains predicted to have a surface location and potential for antibody induction. Although the antigenic peaks of the catalytic traid of the EoSPs showed less similarity with that of the analogous venom SPs particularly those at residues 67 and 110, many antigenic residue similarities of EoSPs are shared with other SVSPs of related vipers. Therefore, it is likely that antibodies raised by EoSPDNA immunisation are likely to possess considerable cross-reactivity and might competitively inhibit the function of these domains in the similar venom toxins of related vipers. However, binding of antibodies specific to conserved antigenic domains without a known function are equally as likely to disrupt protein function by virtue of steric hindrance. The veracity of these speculations need to be confirmed experimentally and thus is a focus of our current research.In conclusion, the predicted Jameson-Wolf antigenic profiles (DNASTAR, USA) of theEoSPvariants aligned with very low identity to their (BLAST) analogous serine proteases. This observation strongly suggests that an antibody raised by immunisation with group one EoSPDNA is likely to be less effective against the gene products of groups 2, 3, or 4. Therefore additional antibodies generated against antigenic index that showed less conservation will be required.



---

*Source: 134232-2010-09-29.xml*