Molecular cloning of a hyaluronidase fromBothrops pauloensisvenom gland

Abstract

Background

Hyaluronate is one of the major components of extracellular matrix from vertebrates whose breakdown is catalyzed by the enzyme hyaluronidase. These enzymes are widely described in snake venoms, in which they facilitate the spreading of the main toxins in the victim’s body during the envenoming. Snake venoms also present some variants (hyaluronidases-like substances) that are probably originated by alternative splicing, even though their relevance in envenomation is still under investigation. Hyaluronidases-like proteins have not yet been purified from any snake venom, but the cDNA that encodes these toxins was already identified in snake venom glands by transcriptomic analysis. Herein, we report the cloning and in silicoanalysis of the first hyaluronidase-like proteins from a Brazilian snake venom.

Methods

The cDNA sequence of hyaluronidase was cloned from the transcriptome of Bothrops pauloensisvenom glands. This sequence was submitted to multiple alignment with other related sequences by ClustalW. A phylogenetic analysis was performed using MEGA 4 software by the neighbor joining (NJ) method.

Results

The cDNA from Bothrops pauloensis venom gland that corresponds to hyaluronidase comprises 1175 bp and codifies a protein containing 194 amino acid residues. The sequence, denominated BpHyase, was identified as hyaluronidase-like since it shows high sequence identities (above 83%) with other described snake venom hyaluronidase-like sequences. Hyaluronidases-like proteins are thought to be products of alternative splicing implicated in deletions of central amino acids, including the catalytic residues. Structure-based sequence alignment of BpHyase to human hyaluronidase hHyal-1 demonstrates a loss of some key secondary structures. The phylogenetic analysis indicates an independent evolution of BpHyal when compared to other hyaluronidases. However, these toxins might share a common ancestor, thus suggesting a broad hyaluronidase-like distribution among venomous snakes.

Conclusions

This work is the first report of a cDNA sequence of hyaluronidase from Brazilian snake venoms. Moreover, the in silico analysis of its deduced amino acid sequence opens new perspectives about the biological function of hyaluronidases-like proteins and may direct further studies comprising their isolation and/or recombinant production, as well as their structural and functional characterization.

Alternative splicing; Hyaluronidase-like; Snake venom

Background

“Hyaluronidase” is a term introduced by Meyer [¹] to denote any enzyme that degrades hyaluronate. These enzymes are present in human testis, spleen, skin, eye, liver, kidney, uterus, placenta, tear, blood and sperm [²]. Hyaluronidases from many animal classes usually share common structural features, conserving all the critically important sites for enzyme activity [³]. Hyaluronidases were identified at first in bovine testis and in bacteria and were subsequently described as “spreading factors” [⁴,⁵].

Hyaluronidases are often found in a diversity of venoms, such as those from snakes, lizards and arthropods (scorpions, spiders, wasps and bees), in which they act as an immunogen. In snake venoms, this enzyme potentiates the toxicity and contributes to local damage at the bite site by affecting the extracellular matrix integrity due to hyaluronate degradation [⁶-⁸].

There are also some hyaluronidase variants, known as hyaluronidase-like proteins, which are products of alternative splicing [⁹]. These truncated forms of hyaluronidases have been recently identified in snake and bee venoms and in human serum as well [³,¹⁰,¹¹]. Ever since, some questions have been raised about their physiological role. A hyaluronidase-like isoform fromVespula vulgarisvenom was hypothesized to act as a lectin by binding to hyaluronate and/or other related substances without degrading them [¹¹].

The isolation and biological characterization of hyaluronidases from snake venoms, includingBothropsgenus, are usually difficult to achieve due to their instability and fast degradation, and because of their relatively low concentration in these biological samples. However, the purification of hyaluronidases fromNaja naja, Agkistrodon contortrix contortrix, Cerastes cerastesandCrotalus durissus terrificusvenoms has already been reported [¹²-¹⁵]. In the present work, we describe the molecular cloning andin silicoanalysis of a cDNA sequence that encodes a hyaluronidase-like protein from theBothrops pauloensisvenom gland. The sequence was compared to other known hyaluronidase-like sequences in order to screen conserved structural features, which may generate perspectives regarding its potential physiological functions and contributions to the envenoming.

Methods

Isolation of cDNA hyaluronidase fromBothrops pauloensisvenom gland

A venom gland from aB. pauloensisadult snake was dissected three days after venom extraction, when transcription is most stimulated [¹⁶]. A cDNA library fromBothrops pauloensisvenom gland, formerly referred to asBothropoides pauloensis, was previously constructed by Rodrigueset al.[¹⁷]. Briefly, the pair of venom glands was homogenized by liquid nitrogen and the total RNA was extracted by the Trizol method (Invitrogen, UK). The mRNA was purified from total RNA by using PolyATract® mRNA Isolation kit (Promega, USA) and the cDNA library was obtained by CloneMiner cDNA Library Construction kit (Invitrogen, UK) with 3 μg of purified mRNA. First and second cDNA strands were synthesized as described by the manufacturer protocols whereas size fractioning of cDNA was carried out in a 1 mL column previously packed with Sephacryl S-500 resin. The cDNA was precipitated with ethanol and then resuspended in 50 mL of milli-Q water and submitted to Polymerase Chain Reaction (PCR). The PCR products were purified and sequenced using DYEnamic ET Terminator Cycle Sequencing kit (GE Healthcare, UK) on a MEGA-BACE 1000 automated DNA sequencer (GE Healthcare, UK).

After the cDNA sequencing, specific primers were designed in order to certify that the full open reading frame (ORF) was obtained. The specific primers were Hyase internal forward (5′-TTGGTGAAACAGCGGCCATG-3′) and Hyase internal reverse (5′-CTTTTCATCCAGCACAATAC-3′). After amplification, the PCR products were analyzed by electrophoresis on 1% agarose gel. The bands containing the PCR products were purified from gel using the Wizard SV Gel and PCR clean up system kit (Promega, Brazil), according to the manufacturer’s specifications. The Ins T/A clone PCR Product kit (Fermentas, Lithuania) was used for rapid cloning of PCR products in pTZ57R/T plasmids. Bacteria colonies were selected on a medium containing ampicillin, IPTG and X-Gal. The recombinant colonies were analyzed by PCR and gel electrophoresis. PCR products were purified and submitted to sequencing using DYEnamic ET Terminator Cycle Sequencing Kit (GE Healthcare, UK) on a MEGA-BACE 1000 automated DNA sequencer (GE Healthcare, UK). The software Base Caller Cimarron 3.12 was used to analyze the electropherograms and generate sequences, which were then aligned in the software Bioedit version 7.0.5.3.

In silicoanalysis of cDNA sequences

Hyaluronidase sequences were searched against the NCBI database (http://www.ncbi.nlm.nih.gov/webcite). BLASTp 2.2.19 was used for scoring the sequence alignments and the maximum e-value obtained was 2e^-100[¹⁸]. The predicted sequence of hyaluronidase fromBothrops pauloensisvenom and other full-length hyaluronidases and hyaluronidase-like sequences from the database were aligned by ClustalW (available in http://www.ebi.ac.uk/Tools/clustalw2/index.htmlwebcite). The putative glycosylation sites were identified by the glycosylation prediction software NetNGlyc (http://www.cbs.dtu.dk/services/NetNGlyc/webcite). Predicted signal peptide cleavage site was determined by the SignalP algorithm [¹⁹]. The theoretical isoelectric point and molecular mass were computed using the tool ProtParam [²⁰].

Phylogenetic tree

Sequences were aligned by ClustalW algorithm and the phylogenetic analysis was performed using the software MEGA 4 by the neighbor joining (NJ) method [²¹]. The evolutionary distances were computed by the JTT matrix-based method. The reliability of NJ trees was evaluated by analyzing 1000 bootstrap replicates. Human hyaluronidase (NP009296.1) was employed as an out-group.

Results and discussion

Since snakes need to kill their prey quickly and efficiently, a systemic delivery of the main venom toxins is required in order to potentiate the lethal effects. Thus, these toxins enter into the circulatory system of the victim with the aid of toxins that degrade the extracellular matrix (ECM) (metalloproteinases, myotoxic phospholipases A₂and hyaluronidases) [¹⁰]. Hyaluronidases have been identified in some snake venoms, such as those fromAgkistrodon acutus, Naja naja, Vipera russelli siamensis, Trimeresurus flavoviridis, Trimeresurus popeorum, Trimeresurus macrops, Trimeresurus albolabris, Agkistrodon contortrix contortrixandCrotalus durissus terrificus[⁶,¹⁵,²²-²⁴].

In this study, we present the amino acid sequence of a hyaluronidase-like (BpHyase) protein deduced from a cDNA obtained fromB. pauloensisvenom gland transcriptome [¹⁷]. Interestingly, the identification of a single truncated hyaluronidase-encoding EST was achieved in an attempt to clone true hyaluronidase, which may reflect its low representation in the venom when compared to other toxin classes. On the other hand, most snake venom-gland transcriptomes reveal the presence of transcripts corresponding to hyaluronidase [¹⁰,²⁵].

The cDNA sequence of hyaluronidase fromB. pauloensisgland, denominated BpHyase, is composed of 1175 bp and codifies 194 amino acid residues for the mature protein, including eight cysteine residues (Figure 1). The full-length sequence of BpHyase comprises an ORF of 582 bp, flanked by a 5′ UTR of 100 bp and a 3′ UTR of 493 bp. The initiating methionine of BpHyase is followed by a predicted signal peptidase I (SPase I) cleavage site at FNG20-VH, which is consistent with the secreted nature of toxins. This prepeptide is believed to initiate the transport of preBpHyase into the endoplasmatic reticulum for glycosylation and is characterized by an N-terminal basic region (Met1-Lys7), a hydrophobic region (Cys8-Phe14) and a polar C-terminal (Leu15-Gly20) [²⁶]. The presence of a prepeptide followed by a start codon indicates that BpHyase is probably translated into protein. Moreover, this prepeptide probably releases a mature protein that, after cleavage, possesses 174 amino acids with a theoretical pI of 9.60 and molecular mass of 19,892.3 Da. Therefore, this prepeptide may act as a signal sequence that directs the protein to the secretory pathway of venom gland cells [²⁵].

Figure 1
Full-length nucleotide and deduced amino acid sequence of BpHyase. The underlined amino acids indicate the inferred sequence of signal peptide which is characterized by an N-terminal basic region, marked in pink, a hydrophobic region (green) and a polar C-terminal (light blue) The nucleotide and amino acid sequences reported herein are available in GenBank with accession numbers GR955246 and FJ654998.1, respectively.

The presence of N-linked glycans is supposed to be necessary for the stabilization of intramolecular folding and the consequent retention of enzymatic activity. Furthermore, changes in glycosylation are likely responsible for the diversity of biological functions exhibited by protein isoforms [²⁷]. In relation to BpHyase, several asparagine residues identified in its sequence could potentially constitute glycosylation sites, thus influencing some physical and chemical parameters of the molecule. The glycosylation prediction algorithm (NetNGlyc) indicated the following glycosylation sites for BpHyase: N101, V102, T103 and N146, A147 and T148 (Figure 2). The glycosylation consensus triad is NXS or T, where X represents any amino acid, except proline[²⁸,²⁹]. However, further structural analyses are of great importance to reveal the residues truly involved in glycosylation.

Figure 2
BpHyase was aligned with truncated hyaluronidases fromEchis carinatus sochureki(HyL1000_Ec; Genbank: ABI33950.1 and ABI33949.1),Echis pyramidum leakeyi(HyL500_Epl; Genbank: ABI33949.1 and ABI33942.1.), andBitis arietans(HyL1000_Ba, HyL750_Ba and Hy-L-500; Genbank: ABI33947.1, ABI33946.1 and ABI33948.1) snake venoms. The putative points of glycosylation are marked in the boxes.

Three cDNA variants of truncated hyaluronidase fromEchis pyramidum leakeyi, Echis carinatus sochurekiandBitis arietansvenom glands were already identified: Hy-L-1000 that encodes the consensus amino- and carboxy-termini with a central deletion of 256 residues, Hy-L-750 that lacks the consensus amino-terminus and Hy-L-500 that lacks the amino-terminus and encodes a shorter carboxy-terminal segment [¹⁰]. Hy-L-1000 is probably translated into a protein without enzymatic activity, while Hy-L-750 and Hy-L-500 represent non-translated transcripts due the absence of an essential translation initiating motif. The inferred protein-coding sequence of BpHyase was classified into the Glycol-Hydro-56 superfamily by protein BLAST analysis, and the highest identity (88%) was presented by truncated hyaluronidase fromEchis carinatus sochurekivenom (GenBank: ABI33950) (Figure 2). In order to confirm its identity, BpHyase was aligned by ClustalW algorithm against other reported hyaluronidase-like sequences from snake venoms, in which the highest sequence identities (above 86%) were observed for Hy-L-1000 truncated hyaluronidases, revealing that BpHyase presents the same residue deletion pattern as these molecules.

It would be tempting to speculate that partial hyaluronidases or hyaluronidases-like proteins represent vestigial enzymes with no activity, since some authors affirmed they lack catalytic residues because of deletions of central residues during their evolutionary history [¹⁰,³⁰,³¹]. The predicted BpHyase amino acid sequence was aligned with other full-length and truncated hyaluronidases from snake venoms, as well as human hyaluronidase (hHyal-1), in order to investigate its deletion pattern (Figure 3). The multiple alignment revealed a substantial deletion of 255 amino acids, starting at residue 52, resulting in the loss of two cysteines, the catalytic (Glu135) and positional residues (Asp133, Tyr206, Tyr253) from full-length viper hyaluronidases. Structural data, site-directed mutagenesis and steady state enzyme kinetics allowed the determination of some important residues for human Hyal-1 catalysis [³²].

Figure 3
Structure-based sequence alignment of mature truncated and complete hyaluronidases. Full-length hyaluronidases fromEchis pyramidum leakeyi(E_pyramidum; Genbank: ABI33941.1), Hy-L-1000 truncated hyaluronidase fromEchis carinatus sochureki(E_carinatus; Genbank: ABI33950.1), full-length hyaluronidase fromBitis arietans(B_arietans; Genbank: ABI33945.1), Hy-L-1000 truncated hyaluronidase fromBitis arietans(B_arietans; Genbank: ABI33947.1) and human hyaluronidase (Human_hyal; Genbank: 2PE4). Key catalytic residue is shown in red and positional residues appear in blue. Cysteine residues are marked in black. Secondary structures were based on the human hyaluronidase crystal structure [³⁸]. Secondary structure elements for human hyaluronidase are shown below the sequences: pink arrows represent β-strands, blue cylinders α-helices, and red cylinders 310 helices.

An essential direct role in chemical catalysis was suggested for Glu131 and a supporting role for Asp129, which was also observed by Arminget al.[³³]. In these cases, the acidic character of the residues is critical for enzymatic activity while Glu131 acts as a proton donor to the hydroxyl group in glycosidic cleavage. These acid residues are also conserved into hyaluronidases fromTrichoderma reesei(Glu212 and Glu217),Bacillus agaradherans(Glu139 and Glu228),Echis ocellatusvenom (Asp133 and Glu 135), as well as in those from bovine testis PH20 (Asp147 and Glu 149) andApis melliferavenom (Asp111 and Glu113) [³⁴-³⁷,¹⁰]. Moreover, Tyr202 and Tyr247 are also essential residues for catalytic activity, since Tyr202 probably binds the substrate and Tyr247 is suggested as coordinating and stabilizing the oxidation during transition state formation[³²].

Figure 3 also indicates that BpHyase encodes a protein containing the consensus amino and carboxi-termini, as well as Hy-L-1000 truncated hyaluronidase variants fromE. c. sochurekiandBitis arietansvenoms. hHyal-1 human hyaluronidase showed a cross-generic sequence conservation of 38% when compared to BpHyase, representing the most similar hyaluronidase that had its structure solved. The alignment of hHyal-1 to BpHyase allowed the mapping of the secondary structures lost by deletions and demonstrated the lack of regions that are probably involved in the formation of β-strands 2 to 7 and α-helices 2 to 6. The crystal structures of bee venom (bvHyal) and human (hHyal-1) hyaluronidases reveal a classical (β/α)₈TIM barrel fold, which is common to many hydrolases [³⁶,³⁸,³⁹]. In bvHyal, the barrel is formed by only seven strands (β1- β7), in contrast to hHyal-1, which presents all the eight strands [³⁶,³⁸]. The alignment of BpHyase and other hyaluronidases and hyaluronidases-like prteins from snake venoms with hHyal-1 demonstrates that they lack the residues responsible for forming the β-strands 2 to 7 (Figure 3). This fact indicates that BpHyase and other snake venom hyaluronidases-like proteins have a different folding pattern than that described for hHyal-1 and bvHyal, although this does not necessarily imply a loss of biological activity during envenoming.

Alternative splicing is a molecular mechanism by which different combinations of exons can be alternatively linked in order to produce different mRNA isoforms. In the globular enzymatic proteins, such as hyaluronidase, a negative selection pressure operates against gene duplication and diversification, since their correct folding is generally more sensitive to mutations [⁴⁰]. Therefore, alternative splicing is the most important source of functional diversity for globular proteins in eukaryotes [⁴¹,⁴²]. Alternative splicing variants of hyaluronidases were previously described in hHyal-1 and hHyal-3 from human prostate cancer cell lines andVespula vulgarisvenom [¹¹,⁴³]. Moreover, researchers hypothesized that spliced variants would have their function silenced and demonstrated that some hyaluronidase splicing variants from bladder tumor tissues would form a complex with true hyaluronidases, displaying the spreading activity and, consequently, regulating the functional aspects of these true hyaluronidases by alternative mRNA splicing [⁴⁴,⁴⁵]. Therefore, spliced variants of hyaluronidases would provide a valuable tool for modeling a metastasis inhibitor. Nevertheless, the impact of alternative splicing needs to be further investigated in hyaluronidases from human and venoms, in order to screen a biotechnological application for these recently discovered enzymes. In the present work, we suggest that BpHyase would also operate in association with a true hyaluronidase fromBothrops pauloensisvenom in order to potentiate its activity as a “spreading factor” during the envenoming.

In addition, it was supposed that hHyal-1 splice variants may play regulatory roles by binding to partner proteins via interaction with the carboxi-terminal HyalEGF-like domain, resulting in modulation of its enzymatic activity [³⁸]. The EGF domain is present in many extracellular proteins and is involved in cell adhesion and cell-cell communication [⁴⁶,⁴⁷]. In mammalian hyaluronidases, the HyalEGF-like domain is characterized by an EGF disulfide bond signature sequence. The three disulfide bonds of the HyalEGF-like domain are possibly responsible for the maintenance of its fold even when the catalytic domain unfolds [³⁸]. Mapping deletions on the BpHyase structure demonstrates that in BpHyal-1 the segment correspondent to the HyalEGF-like domain remains intact (Figure 3), which may confer regulatory roles on these hyaluronidase-like molecules. This result is in agreement with our hypothesis that BpHyase contributes to the spreading of the main toxins into the envenomed body of the victim. On the other hand, these findings have not previously been reported for snake venoms. In this conception, further studies concerning the isolation and/or heterologous expression of hyaluronidase-like toxins would be of great interest for testifying their biological role during envenoming.

The phylogenetic analysis of hyaluronidase-like sequence and true hyaluronidases from the Viperidae family (Figure 4) showed that these molecules form a monophyletic group, indicating a recent divergence among them. The branches amongst BpHyase and other hyaluronidases may be due to point mutations and/or gene duplication, which would result in new amino acid sequences. This approach is based on the description by Futuyma [⁴⁸] in which the homologous genes from different species evolve at much higher rates than others, and also, when a single gene pool changes, it can evolve into new species. BpHyase was visualized in a separate branch, which indicates an independent evolution of this toxin when compared to hyaluronidases and hyaluronidases-like proteins fromBitis, EchisandCerastisgenders [⁴⁹]. The most interesting observation is that hyaluronidases-like proteins are, until now, described only forBothrops pauloensis, Echis carinatusandBitis arietansvenoms, which inhabit different continents. All these observations suggest that hyaluronidases-like proteins may share a common ancestor, thus presenting a broad distribution among venomous snakes.

Figure 4
Evolutionary relationships of snake venom hyaluronidases inferred using the neighbor-joining method.Human hyaluronidase was employed as the out-group.

Conclusions

Herein we have reported the first hyaluronidase-like cDNA sequence from a Brazilian snake venom. Thein silicoanalysis of its deduced amino acid sequence opens new perspectives concerning its biological functions, suggesting a regulatory role that needs to be investigated. Moreover, this work may guide further studies comprising their isolation and/or recombinant production, as well as their detailed structural and functional characterization.

Acknowledgments

The authors gratefully acknowledge the financial support of Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq), Universidade Federal de Uberlândia (UFU), Fundação de Amparo à Pesquisa dos Estados de Minas Gerais (FAPEMIG) and Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES), Brazil. We also thank the technical support of Genetic and Evolution Department from Federal University of São Carlos, Brazil.

References

• ¹
  Meyer K: Hyaluronidases. In The Enzymes. Volume 5 Edited by Boyer PD. New York: Academic Press; 1971:307-320.
• ²
  Menzel EJ, Farr C: Hyaluronidase and its substrate: biochemistry, biological activities and therapeutic uses. Cancer Lett 1998, 131(1):3-11.
• ³
  Csóka AB, Scherer SW, Stern R: Expression analysis of six paralogous human hyaluronidase genes clustered on chromosomes 3p21 and 7q31. Genomics1999, 60(3):356-361.
• ⁴
  Duran-Reynalds F: Exaltation de l’activité de virus vaccinal par les extraits de certains organs. Compt Rend Soc Biol 1928, 9:6-7.
• ⁵
  Duran-Reynalds F: Studies on a certain spreading factor existing in bacteria and its significance for bacterial invasiveness. J Exp Med 1933, 58(2):161-181.
• ⁶
  Girish KS, Jagadeesha DK, Rajeev KB, Kemparaju K: Snake venom hyaluronidase: an evidence for isoforms and extracellular matrix degradation. Mol Cell Biochem 2002, 240:105-110.
• ⁷
  Shashidharamurthy R, Jagadeesha DK, Girish KS, Kemparaju K: Variation in biochemical and pharmacological properties of Indian cobra (Naja naja naja) venom due to geographical distribution. Mol Cell Biochem2002, 229:93-101.
• ⁸
  Girish KS, Kemparaju K: Inhibition of Naja naja venom hyaluronidase: Role in the management of poisonous bite. Life Sci 2006, 78(13):1433-1440.
• ⁹
  de Sá VK, Canavez FC, Silva IA, Srougi M, Leite KRM: Isoforms of hyaluronidases can be a predictor of a prostate cancer of good prognosis. Urol Oncol 2009, 27(4):377-381.
• ¹⁰
  Harrison RA, Ibison F, Wilbrahan D, Wagstaff SC: Identification of cDNA enconding viper venom hyaluronidases: cross-generic sequence conservation of full-length and unusually short variant transcripts. Gene 2007, 392(1–2):22-33.
• ¹¹
  Kolarich D, Loos A, Léonard R, Mach L, Marzban G, Hemmer W, Altmann F: A proteomic study of the major allergens from yellow jacket venoms.Proteomics 2007, 7(10):1615-1623.
• ¹²
  Girish KS, Shashidharamurthy R, Nagaraju S, Gowda TV, Kemparaju K: Isolation and characterization of hyaluronidase a “spreading factor” from Indian cobra (Naja naja) venom. Biochimie 2004, 86(3):193-202.
• ¹³
  Wagstaff SC, Sanz L, Juárez P, Harrison RA, Calvete JJ: Combined snake venomics and venom gland transcriptomic analysis of the ocellated carpet viper. Echis ocellatus. J Proteomics 2009, 71(6):609-623.
• ¹⁴
  Wahby AF, el SM M, El-Mezayen HA, Salama WH, Abdel-Aty AM, Fahmy AS: Egyptian horned viper Cerastes cerastes venom hyaluronidase: Purification, partial characterization and evidence for its action as a spreading factor. Toxicon 2012, 60(8):1380-1399.
• ¹⁵
  Bordon KC, Perino MG, Giglio JR, Arantes EC: Isolation, enzymatic characterization and antiedematogenic activity of the first reported rattlesnake hyaluronidase from Crotalus durissus terrificus venom.Biochimie 2012, 94(12):2740-2748.
• ¹⁶
  Rotenberg D, Bamberger ES, Kochva E: Studies on ribonucleic acid synthesis in the venom glands of Vipera palaestinae Biochem J1971, 121(4):609-612.
• ¹⁷
  Rodrigues RS, Boldrini-França J, Fonseca FP, de la Torre P, Henrique-Silva F, Sanz L, Calvete JJ, Rodrigues VM: Combined snake venomics and venom gland transcriptomic analysis of Bothropoides pauloensis J Proteomics 2012, 75(9):2707-2720.
• ¹⁸
  Altschul SF, Madden TL, Schäffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs.Nucleic Acids Res 1997, 25(17):3389-3402.
• ¹⁹
  Bendtsen D, Nielsen H, von Heijne G, Brunak S: Improved prediction of signal peptides: SignalP 3.0. J Mol Biol 2004, 340(4):783-795.
• ²⁰
  Gasteiger E, Hoogland C, Gattiker A, Duvaud S, Wilkins MR, Appel RD, Bairoch A: Protein Identification and Analysis Tools on the ExPASy Server. In The Proteomics Protocols Handbook Edited by Walker JM. New Jersey: Humana Press; 2005:571-607.
• ²¹
  Tamura K, Dudley J, Nei M, Kumar S: MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Mol Biol Evol 2007, 24(8):1596-1599.
• ²²
  Xu X, Wang XS, Xi XT, Liu J, Huang JT, Lu ZX: Purification and partial characterization of hyaluronidase from five pace snake (Agkistrodon acutus) venom. Toxicon 1982, 20(6):973-981.
• ²³
  Pukrittayakamee S, Warrell DA, Desakorn V, McMichael AJ, White NJ, Bunnag D: The hyaluronidase activities of some Southeast Asian snake venoms.Toxicon 1988, 26(7):629-637.
• ²⁴
  Kudo K, Tu AT: Characterization of hyaluronidase isolated from Agkistrodon contortrix contortrix (Southern copperhead) venom.Arch Biochem Biophys 2001, 386(2):154-162.
• ²⁵
  Fox JW: 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.
• ²⁶
  Ng HC, Ranganathan S, Chua KL, Khoon EK: Cloning and molecular characterization of the first aquatic hyaluronidase, SFHYA1, from the venom of stonefish (Synanceja horrida). Gene 2005, 346:71-81.
• ²⁷
  Li MW, Yudin AI, Robertson KR, Cherr GN, Overstreet JW: Importance of glycosylation and disulfide bonds in hyaluronidase activity of macaque sperm surface PH-20. J Androl 2002, 23(2):211-219.
• ²⁸
  Ronin C, Granier C, Caseti C, Bouchilloux S, Van Rietschoten J: Synthetic substrates for thyroid oligosaccharide transferase. Effects of peptide chain length and modifications in the Asn-Xaa-Thr-region. Eur J Biochem 1981, 118(1):159-164.
• ²⁹
  Cherr GN, Yudin AI, Overstreet JW: The dual functions of GPI-anchored PH-20: hyaluronidase and intracellular signaling. Matrix Biol 2001, 20(8):515-525.
• ³⁰
  Fry BG: From genome to “venome”: molecular origin and evolution of the snake venom proteome inferred from phylogenetic analysis of toxin sequences and related body proteins. Genome Res 2005, 15(3):403-420.
• ³¹
  Fry BG, Scheib H, van der Weerd L, Young B, McNaughtan J, Ramjan SF, Vidal N, Poelmann RE, Norman JA: Evolution of an arsenal: structural and functional diversification of the venom system in the advanced snakes (Caenophidia). Mol Cell Proteomics 2008, 7(2):215-246.
• ³²
  Zhang L, Bharadwaj AG, Casper A, Barkley J, Barycki JJ, Simpson MA: Hyaluronidase activity of human hyal1 requires active site acidic and tyrosine residues. J Biol Chem 2009, 284(14):9433-9442.
• ³³
  Arming S, Strobl B, Wechselberger C, Kreil G: In vitromutagenesis of PH-20 hyaluronidase from human sperm. Eur J Biochem 1997, 247(3):810-814.
• ³⁴
  Davies GJ, Dauter M, Brzozowski AM, Bjornvad ME, Andersen KV, Schülein M: Structure of the Bacillus agaradherans family 5 endoglucanase at 1.6 a and its cellobiose complex at 2.0 Å resolution. Biochemistry 1998, 37(7):1926-1932.
• ³⁵
  Divne C, Stahlberg J, Teeri TT, Jones TA: High-resolution crystal structures reveal how a cellulose chain is bound in the 50 A long tunnel of cellobiohydrolase I from Trichoderma reesei J Mol Biol1998, 275(2):309-325.
• ³⁶
  Marković-Housley K, Miglierini G, Soldatova L, Rizkallah PJ, Müller U, Schirmer T: Crystal structure of hyaluronidase, a major allergen of bee venom.Structure 2000, 8(10):1025-1035.
• ³⁷
  Botzki A, Rigden DJ, Braun S, Nukui M, Salmen S, Hoechstetter J, Bernhardt G, Dove S, Jedrzejas MJ, Buschauer A: L-Ascorbic acid 6-hexadecanoate, a potent hyaluronidase inhibitor. X-ray structure and molecular modeling of enzyme-inhibitor complexes. J Biol Biochem 2004, 279(44):45990-45997.
• ³⁸
  Chao KL, Muthukumar L, Herzberg O: Structure of human hyaluronidase-1, a hyaluronan hydrolyzing enzyme involved in tumor growth and angiogenesis.Biochemistry 2007, 46(23):6911-6920.
• ³⁹
  Henrissat B, Callebaut I, Fabrega S, Lehn P, Mornon JP, Davies G: Conserved catalytic machinery and the prediction of a common fold for several families of glycosyl hydrolases. Proc Natl Acad Sci U S A 1995, 92(15):7090-7094.
• ⁴⁰
  Fry BG, Roelants K, Champagne DE, Scheib H, Tyndall JD, King GF, Nevalainen TJ, Norman JA, Lewis RJ, Norton RS, Renjifo C, de la Vega RC: The toxicogenomic multiverse: convergent recruitment of proteins into animal venoms. Annu Rev Genomics Hum Genet 2009, 10:483-511.
• ⁴¹
  Valenzuela A, Talavera D, Orozco M, de la Cruz X: Alternative splicing mechanisms for the modulation of protein function: conservation between human and other species. J Mol Biol 2004, 335(2):495-502.
• ⁴²
  Stamm S, Ben-Ari S, Rafalska I, Tang Y, Zhang Z, Toiber D, Thanaraj TA, Soreq H:Function of alternative splicing. Gene 2005, 344:1-20.
• ⁴³
  Junker N, Latini S, Petersen LN, Kristjansen PE: Expression and regulation patterns of hyaluronidases in small cell lung cancer and glioma lines. Oncol Rep 2003, 10(3):609-616.
• ⁴⁴
  Lokeshwar VB, Schroeder GL, Carey RI, Soloway MS, Iida N: Regulation of hyaluronidase activity by alternative mRNA splicing. J Biol Chem 2002, 277(37):33654-33663.
• ⁴⁵
  Lokeshwar VB, Estrella V, Lopez L, Kramer M, Gomez P, Soloway MS, Lokeshwar BL:HYAL1-v1, an alternatively spliced variant of HYAL1 hyaluronidase: a negative regulator of bladder cancer. Cancer Res 2006, 66(23):11219-11227.
• ⁴⁶
  Harris RC, Chung E, Coffey RJ: EGF receptor ligands. Exp Cell Res 2003, 284(1):2-13.
• ⁴⁷
  Wouters MA, Rigoutsos I, Chu CK, Feng LL, Sparrow DB, Dunwoodie SL: Evolution of distinct EGF domains with specific functions. Protein Sci2005, 14(4):1091-1103.
• ⁴⁸
  Futuyma DJ: Biologia Evolutiva 2nd edition. Ribeirão Preto: Sociedade Brasileira de Genética/CNPq; 1992.
• ⁴⁹
  Feng L, Gao R, Meng J, Gopalakrishnakone P: Cloning and molecular characterization of BmHYA1, a novel hyaluronidase from the venom of Chinese red scorpionButhus martensi karsch Toxicon 2010, 56(3):474-479.

• Ethics committee approval
  Bothrops pauloensis was donated by the Reptile Sector of the Federal University of Uberlândia, Minas Gerais state, Brazil. The serpentarium is registered in the Brazilian Institute of Environment and Renewable Natural Resources – IBAMA (n. 301286).

Publication Dates

• Publication in this collection
  2014

History

• Received
  17 Jan 2014
• Accepted
  30 May 2014
• Reviewed
  10 June 2014