Table of Contents Author Guidelines Submit a Manuscript
The Scientific World Journal
Volume 2014 (2014), Article ID 186029, 14 pages
http://dx.doi.org/10.1155/2014/186029
Research Article

Molecular Phylogeny and Predicted 3D Structure of Plant beta-D--Acetylhexosaminidase

1Department of Biochemistry and Molecular Biology, University of Rajshahi, Rajshahi 6205, Bangladesh
2Genetic Engineering Laboratory, Department of Molecular Biology, Faculty of Resource Science and Technology, Universiti Malaysia Sarawak, 94300 Kota Samarahan, Sarawak, Malaysia

Received 2 February 2014; Revised 14 June 2014; Accepted 21 June 2014; Published 20 July 2014

Academic Editor: Xiu-feng Wan

Copyright © 2014 Md. Anowar Hossain and Hairul Azman Roslan. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

beta-D--Acetylhexosaminidase, a family 20 glycosyl hydrolase, catalyzes the removal of -1,4-linked -acetylhexosamine residues from oligosaccharides and their conjugates. We constructed phylogenetic tree of -hexosaminidases to analyze the evolutionary history and predicted functions of plant hexosaminidases. Phylogenetic analysis reveals the complex history of evolution of plant -hexosaminidase that can be described by gene duplication events. The 3D structure of tomato -hexosaminidase (-Hex-Sl) was predicted by homology modeling using 1now as a template. Structural conformity studies of the best fit model showed that more than 98% of the residues lie inside the favoured and allowed regions where only 0.9% lie in the unfavourable region. Predicted 3D structure contains 531 amino acids residues with glycosyl hydrolase20b domain-I and glycosyl hydrolase20 superfamily domain-II including the (/)8 barrel in the central part. The and contents of the modeled structure were found to be 33.3% and 12.2%, respectively. Eleven amino acids were found to be involved in ligand-binding site; Asp(330) and Glu(331) could play important roles in enzyme-catalyzed reactions. The predicted model provides a structural framework that can act as a guide to develop a hypothesis for -Hex-Sl mutagenesis experiments for exploring the functions of this class of enzymes in plant kingdom.

1. Introduction

As a part of the study to elucidate the role of free -glycans and de--glycosylation mechanism working in plants, we have already characterized the PNGase, ENGase, -mannosidase and -hexosaminidase at molecular level [13]. The -D--acetylhexosaminidase (EC 3.2.1.52), a member of the glycosyl hydrolase family 20 (GH20), is an enzyme that hydrolyses nonreducing terminal monosaccharide residues of --acetylgalactosaminides and --acetylglucosaminides. It is widely distributed among the animals, insects, plants, fungus, and bacteria. Mammal lysosomal --acetyl-D hexosaminidases are mainly responsible for glycoconjugate degradation in lysosome. HexA is a heterodimer of subunits (encoded by the gene HexA) and (encoded by the gene HexB), whereas HexB is a homodimer of subunits. The subunits arose through a gene duplication event and the primary sequences are approximately 60% identical. Mutational defects that cause -hexosaminidase-A and B deficiency are responsible for Sandhoff and the Tay-Sachs diseases, respectively [4]. Recently, it has been reported that -hexosaminidase is a surrogate marker for renal function in autosomal dominant polycystic kidney disease [5]. In insects, it has been postulated to have specialized physiological functions, including posttranslational modification of -glycans, degradation of glycoconjugates, and egg-sperm recognition, suggesting that these enzymes have rather versatile physiological functions in the growth and development of insects [6]. Mammal --acetyl-D-hexosaminidases have been shown to be important for egg-sperm recognition [7], and the enzymes from Drosophila melanogaster sperm membrane also participate in the same process [8]. A fungal --acetyl-D-hexosaminidases has been expressed, characterized, and crystallized from Aspergillus oryzae, which has sequence similarity to bacterial and human enzymes ranges from 42% to 49% [9].

Recently, plant --acetyl-D-hexosaminidases has gained a lot of attention due to its presence in the ripening stages [2]. It has also been shown that the tomato fruit shelf life can be enhanced by the suppression of -glycan processing/degrading enzymes [10]. Plant glycoproteins contain substantial amounts of paucimannosidic -glycans lacking terminal GlcNAc residues at their nonreducing ends. It has been proposed that this is due to the action of -hexosaminidases during late stages of -glycan processing or in the course of -glycan turnover [11]. Although several -hexosaminidases have been reported from various parts of plants such as leaves, fruits, and seeds [1013], their physiological functions in plant biology are not yet fully understood. To elucidate the exact roles of this enzyme in plant kingdom, it is desirable to know about properties and behavior of the phylogenetically related enzymes from different species and their molecular evolutions. However, little is known about the phylogenetics and evolution of plant -hexosaminidases.

So far eight crystal structures of GH20 --acetyl-D-hexosaminidases have been reported including two humans, one insect, and six bacterial enzymes. Both the human HexA and HexB are the --acetyl-D-hexosaminidases that degrade glycoconjugate in the lysosome [4, 14]. OfHex1, the enzyme from the Asian corn borer Ostrinia furnacalis (one of the most destructive pests), has been reported to function merely in chitin degradation [6]. The bacterial enzymes include SpHex and SmCHB, which are found in the chitinolytic bacteria Streptomyces plicatus and Serratia marcescens, respectively [1517]. AaDspB, which is isolated from Aggregatibacter actinomycetemcomitans, is involved in the degradation of biofilm (polymeric -1,6-linked GlcNAc) [18]. The enzyme, PsHex from Paenibacillus sp. TS12, can efficiently degrade various glycosphingolipids [19]. PgGcnA, the enzyme found in the endocarditis pathogen, Streptococcus gordonii, is involved in the release of dietary carbohydrates [20]. Recently, it has been found that a novel --acetylhexosaminidase, StrH protein from Streptococcus pneumoniae R6, is involved in the catalytic specificity towards the (1,2)-linked --acetylglucosides and key residues in the active site are Trp-443 and Tyr-482 [21]. Thus, it is interesting to know how these enzymes could carry out their specialized functions in terms of their structural features. To our knowledge, no crystal structure of plant --acetyl-D-hexosaminidase has yet been reported. Therefore, comparative homology modeling of tomato --acetyl-D-hexosaminidase is desirable to elucidate the functional prediction, active site information, and mechanism of action.

In the present work, first we identified the 83 homologous sequences of --acetyl-D-hexosaminidase in GenBank by the NCBI BLAST-PSI search. We did multiple sequence alignments and reconstructed the phylogenetic tree. Secondly, in order to initiate structural studies of this enzyme, we performed sequence alignment and 3D-structure homology modeling and constructed a molecular model of this enzyme and of its complex with the natural substrate. We also performed molecular docking of the enzymes and predicted the active site residues responsible for catalytic activity. The predicted 3D structural information will be useful to study the site-directed mutagenesis wet lab experiments as well as the physiological functions of tomato --acetyl-D-hexosaminidase in the plant kingdom.

2. Material and Methods

2.1. Data Retrieval

In this study, we retrieved all of the sequences from the National Center for Biotechnology Information (NCBI) GenBank database as described by Gonzalez and Jordan [22]. Shortly, an initial dataset of the previously published and functionally characterized --acetyl-D-hexosaminidase amino acid sequences was retrieved manually from Entrez (http://www.ncbi.nlm.nih.gov/entrez). The representative sequences including the -Hex-Sl were isolated from a wide phylogenetic range of eukaryotes and prokaryotes, which possessed a variety of biochemical activities. A CD-Hit clustering program was used to group these sequences by amino acid identities into clusters [23]. Divergent --acetyl-D-hexosaminidase amino acid sequences with representatives from each cluster were used as queries in a series of PSI-BLAST (Position-Specific Iterated BLAST) searches of the protein database throughout all organisms at NCBI [24]. The representative sequences were Solanum lycopersicum beta-hexosaminidase sequence [gi:350540008], Arabidopsis thaliana AtHex1[gi:30694211], Homo sapiens protein sequences HexA[gi:4261632] and HexB[gi:867691], Drosophila melanogaster hexosaminidase sequences Hexo1[gi:17647501] and Hexo2[gi:17933586], Aspergillus oryzae HexA[gi:169766420], and Streptomyces plicatus HexA[gi:13786695]. We chose the sequences from BLAST results based on the high similarities of amino acids (>60% identities) with the query representative sequences. The picked sequences were checked manually to exclude incomplete and redundant sequences. For the feature analysis and construction of phylogenetic tree we took a total of 83 sequences, which are already characterized as predicted or true --acetyl-D-hexosaminidase from the GenBank, to reduce computational burden. An archea sequence was also retrieved from GenBank that was used as an outgroup in the construction of phylogenetic tree.

2.2. Multiple Sequence Alignments and Construction of Phylogenetic Tree

MUSCLE program [25] was used to align all 83 amino acid sequences of --acetyl-D-hexosaminidases and the alignments were checked manually. Unambiguously aligned regions were identified using GBlocks program [26]. The phylogenetic relationships between the genes were analyzed using the maximum-likelihood (ML) method. For the ML analyses, we used the PROTML program of PHYLIP version 3.6 [27].

We employed the WAG model of amino acid substitution with gamma distribution site rate and invariable site category for phylogenetic analysis [28]. All indels were counted as missing. We performed ten random sequence addition searches using the J option and global branch swapping using the G option to isolate the ML tree with the best log likelihood. In addition, we performed bootstrap analysis with 100 replications.

2.3. Comparative Homology Modeling of Tomato --Acetyl-D-hexosaminidase

Amino acid sequence of Solanum lycopersicum --acetyl-D-hexosaminidase (-hex-Sl) composed of 575 residues was retrieved from NCBI GenBank (GI: 350540008 and Accession no. NP_001234608.1). The SWISS-MODEL web server [29] was used to identify the template structure, 1now, and also used for homology modeling. The online ModWeb Comparative Modeling Server version SVN.r1340:1348M and I-TASSER [30] were also used for further modeling to compare which is the most correct model. The DFire [31], QMEAN [32], PROCHECK [33], WHAT_CHECK [34], and VERIFY_3D [35] methods and ModEval model evaluation server [36] were used to check the validity of the modeled structures. UCSF Chimera and Swiss-PdbViewer were used to view the models and images preparation. The COFACTOR, a structure-based method for biological function annotation of protein molecules, was used to identify the functional insights including ligand-binding site, gene-ontology terms, and enzyme classification [3739].

3. Results and Discussion

3.1. Sequence Analysis of -Hex-Sl

The -Hex-Sl protein sequence was analyzed by NCBI CD-search tool (CDD V3.0-44354 PSSMs) to identify the conserved domains (CD). The sequence contains a Glyco_hydro_20b (46~149 aa), GH20_HexA_HexB-like domain (167~549 aa), and a glycosyl hydrolase family 20, catalytic domain (167~522 aa) belongs to the GH20_hexosaminidase superfamily proteins (Figure 1). Based on CD database available and three-dimensional structure-activity relationship, the amino acid residues Arg(178), Asp(207), His(261), Asp(330), Glu(331), Trp(378), Trp(404), Tyr(430), Asp(432), Trp(494), and Glu(496) were predicted to be present in the active site of -Hex-Sl with other sequences (Figure 2). The online tool NetNGlyc 1.0 server was used to identify the -glycosylation site present in the protein sequence. The predicted -glycosylation sites were position at 50 (NFTI), 86 (NLTS), 112 (NESY), 151 (NPTR), 299 (NPSI), 350 (NGTL), 362 (NNTL), 372 (NRTV), 390 (NPSL), 409 (NNTK), and 441 (NDSR) (data not shown). The software SignalP 4.1 server was used to predict the signal peptide cleavage site that was found to be in between positions 23 and 24 in the amino acid sequence.

186029.fig.001
Figure 1: Conserved domains for tomato -hexosaminidase, analyzed using Conserved Domain Database search in NCBI-BLAST.
186029.fig.002
Figure 2: Sequence alignment of -Hex-Sl with nine other sequences by CD search. The amino acid residues Arg(178), Asp(207), His(261), Asp(330), Glu(331), Trp(378), Trp(404), Tyr(430), Asp(432), Trp(494), and Glu(496) were predicted to be responsible for the activity of -Hex-S. The conserved amino acids are shown as yellow color.

3.2. Phylogenetic Analysis of -Hexosaminidase Sequences

In order to know the evolutional history and properties of plant beta-hexosaminidases, we reconstructed the phylogenetic tree. We aimed to collect the sequence data of the beta-hexosaminidases from a wide range of organisms so that we could get a lot of information including their physicochemical, structural, and biological functions. A total of 83 amino acid sequences were retrieved from the GenBank database by previously characterized representative sequences. These sequences used in the analysis include 23 experimentally characterized --acetyl-D-hexosaminidase enzymes as well as 60 novel predicted or putative --acetyl-D-hexosaminidase sequences (Table 1). MUSCLE program was also used to align the sequences, whereas maximum likelihood method was used in phylogenetic reconstruction. Our phylogenetic analysis shows that --acetyl-D-hexosaminidases are widely distributed among plant, animal, insects, fungi, and bacteria, belonging to the glycosyl hydrolase 20 superfamily (Figure 3). It reveals the complex history of evolution of --acetyl-D-hexosaminidases that can be described by multiple gene duplication events.

tab1
Table 1: Proteins sequences used for construction of phylogenetic studies.
186029.fig.003
Figure 3: The phylogenetic tree based on beta-hexosaminidase amino acid sequences obtained by the maximum likelihood method. Thermofilum (Archea) was used as an outgroup to reconstruct the phylogenetic tree. The percentage of replicate trees in which the associated taxa clustered together in the bootstrap test (100 replicates) is shown next to the branches. All analyses were performed with the WAG amino acid substitution model and 1 invariable and 4 gamma distributed site rate categories. Detailed information about the sequences is shown in Table 1.

Eukaryotic -hexosaminidases might be originated from common bacterial ancestor through multiple gene duplications. Bacteria and fungi clades mostly contain one gene for hexosaminidase in each species albeit few have two genes. Bacteria clade consists of -hexosaminidases that have the peptidoglycan degradation and chitinolytic activities. Those bacterial species, which contain two genes of hexaminidases, might acquire their last copies either by horizontal gene transfer or gene duplication. Fungi sequences clearly showed its own clade and only few species have more than one gene and might be originated either lineage specific mutation and/or gene duplication. Insects clade-I and clade II and plants clade-I and II also contain at least one hexosaminidase gene in each species. Insects (I and II) clades hexosaminidases are chitinolytic enzymes, which separately form paraphylactic groups that could be evolved by gene duplication. Plants clade-I and clade-II also constitute paraphylactic group and also split into monocotyledons and dicotyledons that have functional divergences. Plant -hexosaminidases are involved in -glycan processing of cell walls. Animal clade clearly splits into two clades, A and B, that contain the isoenzymes, HexA and HexB, respectively.

Gene duplication is considered a major driving force for evolution of genetic novelty, thereby facilitating functional divergence and organismal diversity, including the process of speciation. It can be generated by several mechanisms, including tandem duplication, transposition, and large-scale duplication (e.g., segmental/whole genome duplication (WGD)). Also, segmental duplications (SDs) are increasingly recognized as frequent phenomena, especially in primate genomes; for example, approximately 5% of the human genome consists of duplicated segments [40]. More than 300 gene duplication events have been detected by phylogenetic analysis of plant, animal, and fungi before the separation of three major eukaryotic lineages [41]. Specifically, copy numbers for genes with highly conserved functions seem to be more stable than the number of genes with more divergent functions. beta-Hexosaminidases from each kingdom (plant, animal, and insect) are separated into two clades (clusters) and each clade contains at least one member. Human genome data analysis showed that both genes, HexA and HexB, are located in different locus in the chromosomes-15q23-24 and -5q13, respectively. They are originated by gene duplication [42]. Most of the higher eukaryotes contain two or more genes for the hexosaminidases. For example, Arabidopsis thaliana contains Hex1, Hex2, and Hex3 [11]. Likewise Drosophila melanogaster has three genes, Hexo1, Hexo2, and fdl for hexosaminidase isoenzymes [8]. Even these proteins are also located in different organelles. It has been reported that some legume species have at least two Adh gene loci and resulted from relatively ancient duplication events [43]. From the accumulated evidences and phylogenetic topology, it can be speculated that eukaryotic hexosaminidases might be originated by multiple gene duplication, although more experimental evidences are required to establish our hypothesis.

Most of the prokaryotic and eukaryotic -hexosaminidases reported so far play an important physiological role in chitin recycling, a structural components of cell walls [6, 17, 44]. Plant -hexosaminidases have been investigated in a variety of tissues including seeds and leaves suggesting a role in the storage of glycoproteins [4547]. They have also been proposed to be involved in plants defense mechanisms and reported as chitin-degrading enzymes [46, 48]. A molecular study of Arabidopsis -hexosaminidases has shown that HEXO1 participates in -glycan trimming in the vacuole, whereas HEXO2 and/or HEXO3 could be responsible for the processing of -glycans present on secretory glycoproteins [11]. The -Hex is also present at high levels during the ripening of many fruits, including the climacteric fruit tomato [49] and mango [50]. Recently, it has been reported that suppression of -glycan processing enzymes increases the shelf life of tomato fruits and capsicum [10, 51]. The -Hex, a cell wall enzyme, cleaves the terminal -acetyl-D-hexosamine residues and generates the paucimannosidic -glycans present in most plant glycoproteins which in turn downregulate the genes that encode for certain cell wall degrading proteins, such as pectin methylesterase, glucan endo-1,3--D-glucosidase, -1,3-glucanase, endoxyloglucan transferase, pectinesterase, expansin, pectinacetylesterase, -galactosidase, pectate lyase, (1-4)--mannan endohydrolase, and -galactosidase [10]. Therefore, suppression of -Hex activity in transgenic fruits not only inhibited -glycoprotein degradation but also affects cellulose, hemicellulose, and pectin degradation. Altogether, our phylogenetic analysis of various GH20 -Hexosaminidases with their comparative functional properties suggests that plant -Hexosaminidases are cell wall bound enzymes derived from common bacterial ancestor through multiple gene duplications and are involved in -glycan degradation or processing.

3.3. Resolved Predicted 3D Structure and Function

The SWISS-MODEL web server [29] was used to identify the 1now as template structure for homology modeling with 38.41% the target-template sequence identity. Another online server ModWeb Comparative Modeling Server version SVN.r1340:1348 M and I-TASSER [30] were also used for further modeling for appropriate model selection. To obtain an accurate homology model, it is very important that appropriate steps are built into the process to assess the quality of the model. Therefore, the accuracies of the predicted models were checked through a series of tests such as DFire [31], QMEAN [32], PROCHECK [33], WHAT_CHECK [34], VERIFY_3D [35], and also ModEval Model evaluation server [36]. A high quality predicted model was obtained from ModWeb comparative modeling web server through the analysis of predicted structures when compared with each other. However, the data for the rest of modeled structures are not shown. The Dfire energy and QMEAN score of best model were −716.03 and 0.511, respectively. The Ramachandran plot showed 88.1% of the residues in the most favoured region, 10.4% in the additional allowed region, 0.7% in the generally allowed region, and only 0.9% in the unfavourable region (Figure 4). Ramachandran -score is −0.669 indicating how well the backbone conformations of all residues are corresponding to the known allowed areas in the Ramachandran plot and within expected ranges for a well-refined structure. None of the individual amino acid residues was in a bad packaging region. The structural average for the second-generation quality control value is within the normal range. All contacts average is −0.484 and -score is −2.49, which were within the normal ranges. The Anolea, QMean graph and DSSP (define secondary structure of protein) of modeled -Hex-Sl obtained from the structural assessment by Swiss-model workplace are shown in Figure 5.

186029.fig.004
Figure 4: Ramachandran plot of the modeled structure of tomato --acetyl hexosaminidase provided by PROCHECK.
186029.fig.005
Figure 5: Anolea, Qmean, and DSSP (define secondary structure of protein) obtained from the structural assessment by SWISS-MODEL workplace online software.

The X-ray crystal structure of human -hexosaminidase started at position 55 of its gene-translated protein sequence [4]. However, the 3D modeled structure of -Hex-Sl started at 38 position of its amino acid sequence as N-terminal. An overall structural model of -Hex-Sl is shown in Figure 6(a), which contains 531 residues in structural parts, glycosyl hydrolase 20b domain-I, and glycosyl hydrolase 20 superfamily domain-II including the (/)8 barrel in the middle part. The (/)8 barrel structure houses the active site within loops extending from the C termini of the strands that constitute the -barrel. The homologous domains are found in the crystal structure of S. plicatus (SpHEX) and S. marcescens (SmCHB) [16, 17]. An important secondary-structural motif comprised 19 helices and 13 strands. The - and -contents of the modeled protein were found to be 33.3% and 12.2%, respectively, as predicted by the program PROMOTIF (Figure 6(a)). Structural similarity was further compared by superimposition of modeled structure with template. The modeled structure -Hex-Sl closely resembled the template structure (1nowB) and it had good similarity with the template upon superimposition (Figure 6(b)). The online 3D ligand site prediction software [52] was used to identify the ligand-binding site of the modeled structure -Hex-Sl. The amino acid residues Arg(178), Asp(207), His(261), Asp(330), Glu(331), Trp(378), Trp(404), Tyr(430), Asp(432), Trp(494), and Glu(496) were predicted to be present in the ligand-biding site of -Hex-Sl modeled structure (Figure 6(c)). The space filled view of ligand-biding site of -Hex-Sl with docking substrate -acetyl--D-glucosamine (NAG) is shown in Figure 6(d). The COFACTOR online software was used to identify the functional motifs including ligand-binding site, gene-ontology terms, and enzyme classification. The top 10 structural analogs of -Hex-Sl modeled structure were identified in the protein data bank (Table 2). The 1nowB, which had the TM-score 0.785 and RMSD 2.21, was found to be the top ranked among the various the homologous proteins analyzed (Table 2). The results indicated that our predicted model structure of -Hex-Sl was good, accurate, and reliable.

tab2
Table 2: Top 10 identified structural analogs in PDB by COFACTOR.
fig6
Figure 6: The molecular 3D modeling of tomato beta--acetyl hexosaminidase (-Hex-Sl). SPDB viewer and Chimera were used to prepare the images. (a) The predicted 3D modeled structure is shown as ribbon diagram. The structure contains two fold domains (I and II) including -helix (red), -pleated sheets (purple), and coils (gray) The catalytic domain II is a (/)8 barrel with the active site located at the C terminus of the barrel. Template used for building this structure was 1now_B(PDB). (b) Superimposition magic fit image of the modeled structure -Hex-Sl (blue) with template structure human 1now, human --acetyl-hexosaminidase (red), and human -hexosaminidase B-subunit. (c) The predicted ligand-binding site (active site) residues identified are depicted by as blue color. (d) Space filled view of ligand biding site of -Hex-Sl with docking substrate -acetyl--D-glucosamine (NAG).

The COFACTOR identified -Hex-Sl with the classification EC3.2.1.52 and predicted that amino acid residues Asp(330) and Glu(331) could play important role in enzymatic reaction (Table 3). It was also used to search other known homologous binding to compare the consensus binding with predicted ligand binding site. The three proteins (3lmyA, 2gk1G, 2gjx1) were found to have similar consensus binding sites that were identical to the previously predicted ligand-binding sites (Table 4). To predict the functions of modeled structure of -Hex-Sl, we used COFACTOR and identified 19 gene ontology (GO) terms. The consensus prediction of GO terms and their GO-scores are shown in Table 5. Table 5 shows a consistence of function (GO terms) amongst top scoring templates. The GO score associated with each prediction is defined as the average weight of the GO term, where the weights are assigned based on of the template from which the GO term is derived. The most striking features for -Hex-Sl described by GO terms are homodimerization activities and localization in cell membrane. In humans, two major -hexosaminidase isoenzymes exist: Hex A and Hex B. Hex A is a heterodimer of subunits and (60% identity), whereas Hex B is a homodimer of subunits [4]. The molecular weight of purified -Hex-Sl as determined by gel-filtration (native condition) also showed about four times greater value than that determined by SDS-PAGE (denaturation condition) [10]. This happened due to the dissociation of four subunits from each other by denaturing agent like SDS. The -Hex-Sl modeled 3D structure is a single chain protein containing 531 amino acids but it does not have any other-Hex-subunit like animals. Taken altogether our studies suggested that -Hex-Sl may need to exist as a homotetrameric structure during its functional state and be located at the plant cell wall. Although an involvement of -Hex-Sl in plant cell wall or fruit ripening has been reported recently [10], depending on the properties and behaviour of hexosaminidase homologues we could not exclude the possibilities of their involvements in the other physiological processes such as pathogenic resistance and abiotic stress tolerance in plants.

tab3
Table 3: Top 5 enzyme homologs in PDB by COFACTOR.
tab4
Table 4: Template proteins with similar binding sites searched by COFACTOR.
tab5
Table 5: Consensus prediction of gene ontology terms searched by COFACTOR.

4. Conclusion

We used the 23 previously characterized -hexosaminidases and the 60 novel putative -hexosaminidase amino acid sequences to reconstruct the phylogenetic tree. Phylogenetic analysis placed -Hex-Sl into the plant group, which might originate from the common bacterial ancestral origin by multiple gene duplications. Predicted 3D structure of -Hex-Sl contains 531 amino acids with glycosyl hydrolase 20b domain-I and glycosyl hydrolase 20 superfamily domain-II including the barrel (/)8 in the central part. An important secondary-structural motif comprised 19 helices and 13 strands. The - and -contents of the modeled protein were found to be 33.3% and 12.2%, respectively. Eleven amino acids were found to be involved in ligand-binding site of -Hex-Sl. The amino acid residues Asp(330) and Glu(331) could play important role in enzyme-catalyzed reaction. The fully functional state of -Hex-Sl needs to exist as a tetrameric structure and be located at the plant cell wall. The predicted model provides a structural framework that can act as a guide to develop a functional hypothesis to interpret experimental data of --acetyl-D-hexosaminidases. They may also facilitate efforts to design further site-directed mutagenesis to explore the ligand recognition and the downstream signaling mechanisms for the fruit ripening. The presented modeling approach can be extended to other proteins as well.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgment

The authors are grateful to Abdul Wahab Abideen, a graduate student, Department of Molecular Biology, FRST, UNIMAS, for his help in correction of spelling, grammar, and punctuation of this paper.

References

  1. M. A. Hossain, R. Nakano, K. Nakamura, and Y. Kimura, “Molecular identification and characterization of an acidic peptide:N-glycanase from tomato (Lycopersicum esculentum) fruits,” Journal of Biochemistry, vol. 147, no. 2, pp. 157–165, 2010. View at Publisher · View at Google Scholar · View at Scopus
  2. K. Nakamura, M. Inoue, M. Maeda et al., “Molecular cloning and gene expression analysis of tomato endo-β-n-acetylglucosaminidase, an endoglycosidase involved in the production of high-mannose type free N-glycans during tomato fruit ripening,” Bioscience, Biotechnology and Biochemistry, vol. 73, no. 2, pp. 461–464, 2009. View at Publisher · View at Google Scholar · View at Scopus
  3. M. A. Hossain, R. Nakano, K. Nakamura, M. T. Hossain, and Y. Kimura, “Molecular characterization of plant acidic α-mannosidase, a member of glycosylhydrolase family 38, involved in the turnover of N-glycans during tomato fruit ripening,” Journal of Biochemistry, vol. 148, no. 5, pp. 603–616, 2010. View at Publisher · View at Google Scholar · View at Scopus
  4. B. L. Mark, D. J. Mahuran, M. M. Cherney, D. Zhao, S. Knapp, and M. N. G. James, “Crystal structure of human β-hexosaminidase β: understanding the molecular basis of sandhoff and tay-sachs disease,” Journal of Molecular Biology, vol. 327, no. 5, pp. 1093–1109, 2003. View at Publisher · View at Google Scholar · View at Scopus
  5. H. C. Park, J. H. Hwang, A.-Y. Kang et al., “Urinary N-acetyl-β-D glucosaminidase as a surrogate marker for renal function in autosomal dominant polycystic kidney disease: 1 year prospective cohort study,” BMC Nephrology, vol. 13, no. 1, article 93, 2012. View at Publisher · View at Google Scholar · View at Scopus
  6. T. Liu, H. Zhang, F. Liu, Q. Wu, X. Shen, and Q. Yang, “Structural determinants of an insect β-N-acetyl-D-hexosaminidase specialized as a chitinolytic enzyme,” Journal of Biological Chemistry, vol. 286, no. 6, pp. 4049–4058, 2011. View at Publisher · View at Google Scholar · View at Scopus
  7. P. V. Miranda, F. González-Echeverría, J. A. Blaquier, D. J. Mahuran, and J. G. Tezón, “Evidence for the participation of β-hexosaminidase in human sperm-zona pellucida interaction in vitro,” Molecular Human Reproduction, vol. 6, no. 8, pp. 699–706, 2000. View at Publisher · View at Google Scholar · View at Scopus
  8. F. Cattaneo, J. Intra, M. Matsumoto, F. Briani, M. Hoshi, and M. E. Perotti, “Identification and expression analysis of Drosophila melanogaster genes encoding β-hexosaminidases of the sperm plasma membrane,” Glycobiology, vol. 16, no. 9, pp. 786–800, 2006. View at Publisher · View at Google Scholar · View at Scopus
  9. O. Vaněk, J. Brynda, K. Hofbauerová et al., “Crystallization and diffraction analysis of β-N-acetylhexosaminidase from Aspergillus oryzae,” Acta Crystallographica F: Structural Biology and Crystallization Communications, vol. 67, no. 4, pp. 498–503, 2011. View at Publisher · View at Google Scholar · View at Scopus
  10. V. S. Meli, S. Ghosh, T. N. Prabha, N. Chakraborty, S. Chakraborty, and A. Datta, “Enhancement of fruit shelf life by suppressing N-glycan processing enzymes,” Proceedings of the National Academy of Sciences of the United States of America, vol. 107, no. 6, pp. 2413–2418, 2010. View at Publisher · View at Google Scholar · View at Scopus
  11. R. Strasser, J. S. Bondili, J. Schoberer et al., “Enzymatic properties and subcellular localization of arabidopsis β-N-acetylhexosaminidases,” Plant Physiology, vol. 145, no. 1, pp. 5–16, 2007. View at Publisher · View at Google Scholar · View at Scopus
  12. B. H. Jagadeesh, T. N. Prabha, and K. Srinivasan, “Activities of β-hexosaminidase and α-mannosidase during development and ripening of bell capsicum (Capsicum annuum var. variata),” Plant Science, vol. 167, no. 6, pp. 1263–1271, 2004. View at Publisher · View at Google Scholar · View at Scopus
  13. Y. L. Jin, Y. Y. Jo, K. Y. Kim, J. H. Shim, Y. W. Kim, and R. D. Park, “Purification and characterization of β-N-acetylhexosaminidase from rice seeds.,” Journal of biochemistry and molecular biology, vol. 35, no. 3, pp. 313–319, 2002. View at Publisher · View at Google Scholar · View at Scopus
  14. T. Maier, N. Strater, C. G. Schuette, R. Klingenstein, K. Sandhoff, and W. Saenger, “The X-ray crystal structure of human β-hexosaminidase B provides new insights into Sandhoff disease,” Journal of Molecular Biology, vol. 328, no. 3, pp. 669–681, 2003. View at Publisher · View at Google Scholar · View at Scopus
  15. B. L. Mark, D. J. Vocadlo, S. Knapp, B. L. Triggs-Raine, S. G. Withers, and M. N. G. James, “Crystallographic evidence for substrate-assisted catalysis in a bacterial β-hexosaminidase,” The Journal of Biological Chemistry, vol. 276, no. 13, pp. 10330–10337, 2001. View at Publisher · View at Google Scholar · View at Scopus
  16. S. J. Williams, B. L. Mark, D. J. Vocadlo, M. N. G. James, and S. G. Withers, “Aspartate 313 in the Streptomyces plicatus hexosaminidase plays a critical role in substrate-assisted catalysis by orienting the 2-acetamido group and stabilizing the transition state,” Journal of Biological Chemistry, vol. 277, no. 42, pp. 40055–40065, 2002. View at Publisher · View at Google Scholar · View at Scopus
  17. G. Prag, Y. Papanikolau, G. Tavlas, C. E. Vorgias, K. Petratos, and A. B. Oppenheim, “Structures of chitobiase mutants complexed with the substrate di-N-acetyl-D-glucosamine: the catalytic role of the conserved acidic pair, aspartate 539 and glutamate 540,” Journal of Molecular Biology, vol. 300, no. 3, pp. 611–617, 2000. View at Publisher · View at Google Scholar · View at Scopus
  18. N. Ramasubbu, L. M. Thomas, C. Ragunath, and J. B. Kaplan, “Structural analysis of dispersin B, a biofilm-releasing glycoside hydrolase from the periodontopathogen Actinobacillus actinomycetemcomitans,” Journal of Molecular Biology, vol. 349, no. 3, pp. 475–486, 2005. View at Publisher · View at Google Scholar · View at Scopus
  19. T. Sumida, R. Ishii, T. Yanagisawa, S. Yokoyama, and M. Ito, “Molecular cloning and crystal structural analysis of a novel β-N- acetylhexosaminidase from Paenibacillus sp. TS12 capable of degrading glycosphingolipids,” Journal of Molecular Biology, vol. 392, no. 1, pp. 87–99, 2009. View at Publisher · View at Google Scholar · View at Scopus
  20. D. B. Langley, D. W. S. Harty, N. A. Jacques, N. Hunter, J. M. Guss, and C. A. Collyer, “Structure of N-acetyl-β-D-glucosaminidase (GcnA) from the endocarditis pathogen Streptococcus gordonii and its complex with the mechanism-based inhibitor NAG-thiazoline,” Journal of Molecular Biology, vol. 377, no. 1, pp. 104–116, 2008. View at Publisher · View at Google Scholar · View at Scopus
  21. Y. L. Jiang, W. L. Yu, J. W. Zhang et al., “Structural basis for the substrate specificity of a novel β-N-acetylhexosaminidase StrH protein from Streptococcus pneumoniae R6,” The Journal of Biological Chemistry, vol. 286, no. 50, pp. 43004–43012, 2011. View at Publisher · View at Google Scholar · View at Scopus
  22. D. S. Gonzalez and I. K. Jordan, “The α-mannosidases: phylogeny and adaptive diversification,” Molecular Biology and Evolution, vol. 17, no. 2, pp. 292–300, 2000. View at Publisher · View at Google Scholar · View at Scopus
  23. Y. Huang, B. Niu, Y. Gao, L. Fu, and W. Li, “CD-HIT Suite: a web server for clustering and comparing biological sequences,” Bioinformatics, vol. 26, no. 5, pp. 680–682, 2010. View at Publisher · View at Google Scholar · View at Scopus
  24. D. L. Wheeler, T. Barrett, D. A. Benson et al., “Database resources of the National Center for Biotechnology Information,” Nucleic Acids Research, vol. 35, no. 1, pp. D5–D12, 2007. View at Publisher · View at Google Scholar · View at Scopus
  25. R. C. Edgar, “MUSCLE: a multiple sequence alignment method with reduced time and space complexity,” BMC Bioinformatics, vol. 5, article 113, 2004. View at Publisher · View at Google Scholar · View at Scopus
  26. J. Castresana, “Selection of conserved blocks from multiple alignments for their use in phylogenetic analysis,” Molecular Biology and Evolution, vol. 17, no. 4, pp. 540–552, 2000. View at Publisher · View at Google Scholar · View at Scopus
  27. J. Felsenstein, PHYLIP: Phylogeny Inference Package, Version 3.6 (Alpha), University of Washington, Seattle, Wash, USA, 2000.
  28. S. Whelan and N. Goldman, “A general empirical model of protein evolution derived from multiple protein families using a maximum-likelihood approach,” Molecular Biology and Evolution, vol. 18, no. 5, pp. 691–699, 2001. View at Publisher · View at Google Scholar · View at Scopus
  29. K. Arnold, L. Bordoli, J. Kopp, and T. Schwede, “The SWISS-MODEL workspace: a web-based environment for protein structure homology modelling,” Bioinformatics, vol. 22, no. 2, pp. 195–201, 2006. View at Publisher · View at Google Scholar · View at Scopus
  30. Y. Zhang, “I-TASSER server for protein 3D structure prediction,” BMC Bioinformatics, vol. 9, article 40, 2008. View at Publisher · View at Google Scholar · View at Scopus
  31. H. Zhou and Y. Zhou, “Distance-scaled, finite ideal-gas reference state improves structure-derived potentials of mean force for structure selection and stability prediction,” Protein Science, vol. 11, no. 11, pp. 2714–2726, 2002. View at Google Scholar
  32. P. Benkert, M. Biasini, and T. Schwede, “Toward the estimation of the absolute quality of individual protein structure models,” Bioinformatics, vol. 27, no. 3, Article ID btq662, pp. 343–350, 2011. View at Publisher · View at Google Scholar · View at Scopus
  33. R. A. Laskowski, M. W. MacArthur, D. S. Moss, and J. M. Thornton, “PROCHECK: a program to check the stereochemical quality of protein structures,” Journal of Applied Crystallography, vol. 26, pp. 283–291, 1993. View at Google Scholar
  34. R. W. W. Hooft, G. Vriend, C. Sander, and E. E. Abola, “Errors in protein structures,” Nature, vol. 381, no. 6580, p. 272, 1996. View at Google Scholar · View at Scopus
  35. R. Luthy, J. U. Bowie, and D. Eisenberg, “Assesment of protein models with three-dimensional profiles,” Nature, vol. 356, no. 6364, pp. 83–85, 1992. View at Publisher · View at Google Scholar · View at Scopus
  36. F. Melo, R. Sánchez, and A. Sali, “Statistical potentials for fold assessment,” Protein Science, vol. 11, no. 2, pp. 430–448, 2002. View at Publisher · View at Google Scholar · View at Scopus
  37. A. Roy, J. Yang, and Y. Zhang, “COFACTOR: an accurate comparative algorithm for structure-based protein function annotation,” Nucleic Acids Research, vol. 40, no. 1, pp. W471–W477, 2012. View at Publisher · View at Google Scholar · View at Scopus
  38. A. Roy and Y. Zhang, “Recognizing protein-ligand binding sites by global structural alignment and local geometry refinement,” Structure, vol. 20, no. 6, pp. 987–997, 2012. View at Publisher · View at Google Scholar · View at Scopus
  39. J. Yang, A. Roy, and Y. Zhang, “BioLiP: a semi-manually curated database for biologically relevant ligand-protein interactions,” Nucleic Acids Research, vol. 41, no. 1, pp. D1096–D1103, 2013. View at Publisher · View at Google Scholar · View at Scopus
  40. T. Marques-Bonet, S. Girirajan, and E. E. Eichler, “The origins and impact of primate segmental duplications,” Trends in Genetics, vol. 25, no. 10, pp. 443–454, 2009. View at Publisher · View at Google Scholar · View at Scopus
  41. X. Zhou, Z. Lin, and H. Ma, “Phylogenetic detection of numerous gene duplications shared by animals, fungi and plants,” Genome Biology, vol. 11, no. 4, article r38, 2010. View at Publisher · View at Google Scholar · View at Scopus
  42. K. Zwierz, A. Zalewska, and W. Zoch-Zwierz, “Isoenzymes of N-acetyl-β-hexosaminidase,” Acta Biochimica Polonica, vol. 46, no. 3, pp. 739–751, 1999. View at Google Scholar · View at Scopus
  43. T. Fukuda, J. Yokoyama, T. Nakamura et al., “Molecular phylogeny and evolution of alcohol dehydrogenase (Adh) genes in legumes,” BMC Plant Biology, vol. 5, article 6, 2005. View at Publisher · View at Google Scholar · View at Scopus
  44. C. E. Bulawa, “Genetics and molecular biology of chitin synthesis in fungi,” Annual Review of Microbiology, vol. 47, pp. 505–534, 1993. View at Publisher · View at Google Scholar · View at Scopus
  45. M. S. Barber and J. P. Ride, “Purification and properties of a wheat leaf N-acetyl-β-d-hexosaminidase,” Plant Science, vol. 60, no. 2, pp. 163–172, 1989. View at Publisher · View at Google Scholar · View at Scopus
  46. A. Oikawa, E. Itoh, A. Ishihara, and H. Iwamura, “Purification and characterization of β-N-acetylhexosaminidase from maize seedlings,” Journal of Plant Physiology, vol. 160, no. 9, pp. 991–999, 2003. View at Publisher · View at Google Scholar · View at Scopus
  47. N. Harris and M. Chrispeels, “Histochemical and biochemical observations on storage protein metabolism and protein body autolysis in cotyledons of germinating mung beans,” Plant Physiology, vol. 56, no. 2, pp. 292–299, 1975. View at Google Scholar
  48. P. F. Dowd, E. T. Johnson, and T. S. Pinkerton, “Oral toxicity of β-N-acetyl hexosaminidase to insects,” Journal of Agricultural and Food Chemistry, vol. 55, no. 9, pp. 3421–3428, 2007. View at Publisher · View at Google Scholar · View at Scopus
  49. B. H. Jagadeesh, T. N. Prabha, and K. Srinivasan, “Activities of glycosidases during fruit development and ripening of tomato (Lycopersicum esculantum L.): implication in fruit ripening,” Plant Science, vol. 166, no. 6, pp. 1451–1459, 2004. View at Publisher · View at Google Scholar · View at Scopus
  50. M. A. Hossain, M. M. Ranam, Y. Kimura, and H. A. Roslan, “Changes in biochemical Character istics and activities of ripening associated enzymes in mango fruit during storage at different temperatures,” BioMed Research Internal. In press.
  51. S. Ghosh, V. S. Meli, A. Kumar et al., “The N-glycan processing enzymes α-mannosidase and β-D-N-acetylhexosaminidase are involved in ripening-associated softening in the non-climacteric fruits of capsicum,” Journal of Experimental Botany, vol. 62, no. 2, pp. 571–582, 2011. View at Publisher · View at Google Scholar · View at Scopus
  52. M. N. Wass, L. A. Kelley, and M. J. E. Sternberg, “3DLigandSite: predicting ligand-binding sites using similar structures,” Nucleic Acids Research, vol. 38, no. 2, pp. W469–W473, 2010. View at Publisher · View at Google Scholar · View at Scopus