Research article | Open | Published:
Purification, biochemical characterization and structural modelling of alkali-stable β-1,4-xylan xylanohydrolase from Aspergillus fumigatus R1 isolated from soil
BMC Biotechnologyvolume 16, Article number: 11 (2016)
Aspergillus fumigatus R1 produced xylanase under submerged fermentation which degrades the complex hemicelluloses contained in agricultural substrates. Xylanases have gained considerable attention because of their tremendous applications in industries. The purpose of our study was to purify xylanase and study its biochemical properties. We have predicted the secondary structure of purified xylanase and evaluated its active site residues and substrate binding sites based on the global and local structural similarity.
Various microorganisms were isolated from Puducherry soil and screened by Congo-red test. The best isolate was identified to be Aspergillus fumigatus R1. The production kinetics showed the highest xylanase production (208 IU/ml) by this organism in 96 h using 1 % rice bran as the only carbon source. The purification of extracellular xylanase was carried out by fractional ammonium sulphate precipitation (30–55 %), followed by extensive dialysis and Bio-Gel P-60 Gel-filtration chromatography. The enzyme was purified 58.10 folds with a specific activity of 38196.22 IU/mg. The biochemical characterization of the pure enzyme was carried out for its optimum pH and temperature (5.0 and 500C), pH and temperature stability, molecular mass (Mr) (24.5 kDa) and pI (6.29). The complete sequence of protein was obtained by mass spectrometry analysis. Apparent Km and Vmax values of the xylanase for birchwood xylan were 11.66 mg/ml and 87.6 μmol min−1 mg−1 respectively.
Purified xylanase was analyzed by mass-spectrometry which revealed 2 unique peptides. Xylanase under current study showed significant production using agricultural residues and a broad range of pH stability in the alkaline region. Xylanase produced by Aspergillus fumigatus R1 could serve as the enzyme of choice in industries.
Xylan, is a major structural polysaccharide in plant cells and the second most abundant polysaccharide in nature accounting for approximately one-third of all renewable organic carbon on earth  Hemicellulose is a complex of polymeric carbohydrates including xylan, xyloglucan, glucomannan, galactoglucomannan and arabinogalactan . This together with cellulose and lignin constitute the major polymeric constituents of plant cell walls [3, 4]. Within the cell wall structure, all three constituents interact via covalent and non-covalent linkages with the xylan being found at the surface between lignin and cellulose where it is believed to be significant for fibre cohesion and plant cell wall integrity. Due to the heterogeneous and complex structure, the complete hydrolysis of xylan requires a number of enzymes which carry out its hydrolysis by a phenomenon known as cooperativity exhibited by endo-1,4-β-xylanase (EC 220.127.116.11), β-D-xylosidase (18.104.22.168) and a series of enzymes that hydrolyze the side chain groups . Based on the mode of action and mechanism of catalysis xylanases have been assigned to glycoside hydrolase (GH) families 10 and 11 however some xylanases have also been categorized to families 5, 7, 8 and 43 .
Agricultural substrates contain hemicelluloses which have been globally generated. The use of agricultural residues in bioprocesses not only helps to reduce the pollution problems but also serves as the alternate substrate for the production of secondary metabolites . Fungi such as Aspergillus spp. and Trichoderma spp. are of great significance as they produce high levels of xylanase than yeast and bacteria . In this study, xylanase was produced using agricultural waste by Aspergillus fumigatus R1 which was isolated and identified in our laboratory. The xylanase was purified and its biochemical properties were studied to establish its relationship with the other characterized xylanases. Mass-spectrometric analysis reported that the purified xylanase is putative. Moreover, purified xylanase demonstrated significant biochemical characteristics which make it a suitable industrial enzyme.
Identification and Phylogenetic analysis of the culture R1
The xylanase-producing fungus was isolated on agar plate supplemented with the birchwood xylan as the sole source of the carbon at 37 0C. The isolate was identified based on the highest xylanase activities as shown by the clear zone of hydrolysis when flooded with 0.1 % Cong-red stain. The best isolate was selected for further studies. Morphological characteristics indicated that the isolate is Aspergillus. Phylogenetic analysis of 18S rRNA gene sequence indicated that this strain shared 99 % sequence identity with 18S rRNA gene sequence of Aspergillus fumigatus, hence it was assigned as Aspergillus fumigatus strain R1 (Fig. 1).
Xylanase production kinetics
A submerged fermentation process was carried out with the optimized production medium as mentioned above. It was observed that Aspergillus fumigatus R1 produced xylanase and the highest yield reached 208 U/ml after 96 h of cultivation with rice bran as the only carbon source and decreased after 120 h (Fig. 2).
Purification of Xylanase
The xylanase activity increased with the increase of incubation time and reached the peak of activity about 208 U/ml after 96 h of incubation. So, the enzyme was harvested for the purification. Aspergillus fumigatus R1 xylanase was purified by fractional ammonium sulphate precipitation, dialysis and Bio-Gel P-60 gel-filtration chromatography. The result of the purification is summarized in Table 1. The crude supernatant was first subjected to 30–55 % fractional ammonium salt precipitation. The pellet obtained was then dissolved in minimum volume of phosphate buffer (pH 7) and subjected to extensive dialysis. After ammonium sulphate precipitation, the protein with a specific activity of 6300.28 U/mg was obtained whilst the specific activity of the enzyme preparation increased 9.04-fold. The dialysed enzyme sample was then subjected to Bio-Gel P-60 gel-filtration chromatography and the elution profile displayed more peaks of protein whilst a single elution peak of xylanase was obtained which indicated that xylanase of higher purity was achieved (Fig. 3). After the final step of purification, the enzyme was purified 58.10-fold with a specific activity of 38196.22 U/mg with a 3.43 % yield.
SDS-Polyacrylamide gel electrophoresis and activity analysis
SDS-PAGE of the purified xylanase was performed with linear (5–15 %) acrylamide gradient revealed single band when stained by silver staining method (supplementary data). The relative molecular mass of xylanase was found to be approximately 25 kDa. Wong et al.  illustrated a conserved relationship between molecular weight and pI of multiple xylanases from various microorganisms. For instance they are classified as low MW (16–22 kDa)/ basic (pI, 8.3–10.0) xylanases and high MW (43–50 kDa)/ acidic (pI, 3.6-4.5) xylanases. However, this system does not classify all xylanases from different microorganisms satisfactorily. The pure xylanase exhibited relatively clear activity under nondenaturating conditions using 2 % birchwood xylan in 12 % (w/v) acrylamide in gels. The activity was detected by Congo-red staining which indicated that the xylanase was active (supplementary data).
Mass-spectrometry and Sequence analysis
Mass spectrometry data showed 228 amino acids with 2 unique peptides of protein (supplementary data). The molecular weight and isoelectric point of the protein were estimated to be 24.5 kDa and 6.29 respectively. Putative proteins are identified using the bioinformatics tools while they exhibit sequence similarity to an extent with characterized proteins. Usually, the similarity is found only in their conserved amino acid residues and no other significant region is similar to annotated proteins . The amino acid sequence similarity and molecular structure of the protein define the enzyme classification system. 133 different glycoside hydrolase families have been reported [http://www.cazy.org] to date. Enzymes within a specific family have similar three-dimensional structure  and similar molecular mechanism . Xylanases have been classified to families 10 and 11 exclusively where family 11 contains monospecific xylanases.
Homology modeling of Xylanase
To predict the secondary structure of the protein sequence of xylanase, profile-profile alignment was performed with 5798 non-redundant protein structures in I-TASSER (Iterative Threading ASSEmbly Refinement) employing PSI-BLAST (Position-Specific Iterative Basic Local Alignment Search Tool). The secondary structure was hence predicted using PSI-PRED. All α-helices, β-sheets and coiling were identified with their confidence scores in the predicted structure (Fig. 4). The query sequence was threaded through PDB library using LOMETS after the prediction of secondary structure by PSI-PRED. Initially to obtain the best performance from the ten state-of-art treading programs in LOMETS, top ten models were elected from each program solely based on sequence identity coverage and combined energy Z-score values (Table 2). Fold 2vgdA, the top ranked threading fold was identified by three individual threading programs in LOMETS with 66 % sequence identity of the templates in the threading alignment region with the query sequence. The threading folds were considered based on the Z-scores. More the Z-score value, more the significance level of the predicted structure. As it’s evident from the Table 2 that many folds have Z-score values, for instance 3zseA (5.63) or 1xndA (3.76) have more than the Z-score value of 2vgdA which creates ambiguity to predict the true positive structural fold of the protein. Hence, to identify the false positive matches and highlight true positive folds, another programming template modeling align (TM-align) was performed. The result of this structural alignment was ranked according to the TM-score (Table 3). Fold 2vgdA was thus identified and ranked in the first order (TM-score = 0.833; RMSD = 0.61). Hence, 2vgdA was identified as the possible native structure of the query sequence.
Prediction of 3D structure
Five top ranking 3D models were generated by the I-TASSER server by employing the top ten threading folds as the native structures (Fig. 5). Each model was validated based on C-score (Confidence-Score), TM-score, RMSD (The root-mean-square deviation) and cluster density. The model 1 had the highest C-score (-0.93) value reflecting a model of better validation with the structural similarity between the predicted model and the native structure (TM-score = 0.60 ± 0.14 and RMSD = 7.6 ± 4.3Ao).
Prediction of EC number and Ligand Binding site
Model 1 which has the highest C-score (−0.93) was then analyzed by some benchmarking tests on the I-TASSER server for potential enzyme analogs prediction. The analysis showed that the model 1 predicted by I-TASSER matched with five enzyme analogs, all of which exhibited identical EC nos.: 22.214.171.124 (Table 4) whereas the EC-scores (Enzyme Commission-Score) of those analogs differed. Highest EC-score was recorded in the PDB hit 2vgdA was 0.591 whereas the lowest was 0.582 for PDB hits 1hixB and 1te1B. As per ExPaSy enzyme database, the name of the protein analog 2vgdA is Endo-1,4-beta-xylanase whose biological function is endohydrolysis of 1,4-beta-D-xylosidic linkages in xylan. The ligand binding sites of the query protein were predicted using I-TASSER for further verification of protein function. This is based on determination of binding site scores (BS-score) which depend on the structural similarity and local sequence between the query structure and the binding sites of the templates. Based on the BS-score (more than 1.1 indicates prediction with high confidence) residues Asn-109, Tyr-111, Tyr-134 and Tyr-217 were identified as the active binding residues of the query protein (data not shown).
Biochemical characterization of purified xylanase
Effect of pH on xylanase activity
The enzyme activity is greatly affected by pH. For the substrate binding and catalysis depend on charge distribution of both substrate and the enzyme molecules. The reaction pHs were adjusted to 4.0–10.0 with various buffers as described above. Xylanase was active over wide range of pH values from 4.0–10.0 with an optimum pH of 5.0; the optimal acetate buffer pH was 5.0 and retained more than 90 % of its activity (Fig. 6a).
Effect of pH on xylanase stability
The stability of the purified xylanase from Aspergillus fumigatus R1 was determined at 37 0C for 1 h when tested over a wide pH range of 4.0–10.0 (Fig. 6a). The xylanase was found to be highly stable at pH values of 8.0 and 9.0. The enzyme retained 70 % of its activity at pH 9.0 for 1 h at 370C.
Optimum temperature for xylanase activity
The xylanase was active over a broad temperature range of 40–60 0C and showed its optimal temperature at 500C (Fig. 6b).
Thermal stability of xylanase
The enzyme was found to be more stable at temperature below 500C. The enzyme retained about 40 % of its activity at 600C after 30 min of incubation (Fig. 6b). Thermal stability of the enzyme is considered to be one of the significant characteristics in the industry; the thermostability of the enzyme can be increased by different methods like addition of metal ions or by protein engineering.
Determination of kinetic parameters
Km is the Michaelis constant, can be found by measuring the substrate concentration when half the max velocity is achieved. Km is a constant that remains same for a given set of enzyme and substrate. Therefore, low Km increases the affinity of enzyme with the substrate. The rate dependence of the enzyme reaction on different birchwood xylan concentrations at pH 5.0 and 50 0C followed Michaelis-Menten reaction kinetics. Reciprocal plots (Fig. 7) displayed apparent Km and Vmax values of 11.66 mg ml−1 and 87.6 μmol min−1 mg−1. In this study, the purified enzyme showed significantly low Km towards birchwood xylan.
In spite of advanced knowledge of microbial xylanases in the past decades, several terms are still significantly considered for choosing a microorganism for the production of xylanases. Pulping industries require cost effective ways to produce xylanases which must have significant biochemical characteristics to withstand the hostile catalytic conditions. Agricultural residues serve as the prominent source of hemicelluloses which are both environment friendly and cost effective. In the present study, the xylanase has been produced using agricultural residue (rice bran) as the only source of carbon from Aspergillus fumigatus R1 which yielded 208U/ml of xylanase in 96 h of incubation. Aspergillus species have been the target organisms for the production of xylanases mainly because of the production capacity using the submerged fermentation as compared to the bacteria where the xylanase production is found to be low . Xylanase production has been reported from many Aspergillus species. The xylanase activities of Aspergillus fumigatus Z5  and Aspergillus fumigatus MA28  had been found to be 15.02U/ml and 8.45U/ml respectively whereas the xylanase activity of Aspergillus fumigatus R1 was reported to be 208U/ml in 96 h of incubation. Thus, xylanase production by Aspergillus fumigatus R1 is found to be energy efficient and less time consuming which is an important criterion for selection of strain in industry for cost-effective production.
For in silico studies of the purified xylanase, the protein sequence containing 228 amino acids obtained after mass-spectrometry was analyzed by using the I-TASSER server . I-TASSER protein structural modelling is based on the programme known as secondary structure enhanced Profile-Profile threading alignment. Secondary structure of the target protein was predicted by employing the I-TASSER and then a 3D structural model of the target protein indicating substrate binding residues was constructed. Nuclear membrane resonance (NMR) and/or protein crystallography have been used to construct 3D structural model of the target protein  but these techniques are more laborious, expensive and time consuming whereas, bioinformatics analyses offer comparatively less laborious, inexpensive and relatively faster experimentation. However, there are some shortcomings in these analyses; for instance, it is inevitably difficult to select most appropriate candidate when several variable parameters are to be considered in determination of structural models. In the current study, the threading fold 2vgdA (Z-Score 3.37) was selected as a native structure for further analysis of structural modelling (Table 2) where the other folds 3zseA and 1xndA shown higher Z-Score values (5.63 and 3.55 respectively Table 2). In this particular case the decision making factors were the sequence identities (Iden 1 and Iden 2). Thus, there is no problem in selection of a model when all required parameters are significantly higher. In this experiment, model 1 was selected as the top ranked model by I-TASSER which was based on the highest confidence score (C-score = -0.93; Fig. 5a) which was further validated by highest TM-score (0.60 ± 0.14) and RMSD (7.6 ± 4.3Ao). Therefore, this model was selected for further studies. Hence, in the in silico studies the final model of purified xylanase predicted consists of residues Asn-109, Tyr-111, Tyr-134 and Tyr-217 which were identified as the active binding residues of the query protein.
Xylanase from Aspergillus ficuum AF-98 had optimum pH of 5.0 which retained 50 % and 30 % of its relative activities at pH 5.0 and pH 3.0 respectively . The optimum pH of purified xylanase from A. fumigatus R1 is found to be close to other species of Aspergillus such as A. kawachii pH 5.5 , A. nidulans pH 6.0 , A. foetidus pH 5.3 . The change in behaviour of the enzyme at pH 5 and 6 may be due to change in the composition of the buffer. The pH stability of xylanase between pH 7.0–9.0 (370C) has been reported from other Aspergillus species such as A. carneus M34  and xylanase stability at pH 2.0–11.0 (300C), 5.2–5.7 (550C) has also been investigated from different genera such as Streptomyces olivaceoviridis E-86  and Streptomyces olivaceoviridis A1  respectively. Thus, stability of the purified xylanase at broad pH range, especially to alkaline conditions is significant characteristic for applications in industry.
The optimum temperature of xylanase from Aspergillus fumigatus R1 was found to be similar to the xylanase from Aspergillus fumigatus MA28  which showed optimal temperature of 500C. The optimal temperature of purified xylanase from various strains of Aspergillus species is different such as A. versicolor (550C) , A. caespitosus (50–550C) , A. ficuum AF-98 (450C)  and Penicillium species such as Penicillium occitanis Pol6 (450C) . The temperature stability of xylanase in the current study is similar to another xylanase from Penicillium occitanis Pol6  which retained 50 % of its activity at 600C after 1 h incubation. However, incubation at higher temperature decreased the xylanolytic activity. At 700C, it lost 28 % of its activity in 30 min.
Kinetic parameter analysis of purified xylanase from A. fumigatus R1 showed that it had a Km of 11.66 mg ml−1 and Vmax of 87.6 μmol min−1 mg−1. Bakira et al.  studied the kinetic parameter of xylanase from Rhizopus oryzae and showed that Km value for xylanase towards birchwood xylan was 18.5 mg ml−1. The Km value of xylanase from A. fumigatus R1 was found to be similar to the Km value of xylanase from Streptomyces cyaneus SN 32  which was found to be 11.1 mg ml−1 towards birchwood xylan. The Km values of xylanase II from A. sydowii SBS 45  towards birchwood xylan and oat spelt xylan were reported to be 6.51 mg ml−1 and 7.69 mg ml−1 respectively whereas the Vmax of xylanase II for birchwood xylan and oat spelt xylan was found to be 1587 μmol min−1 mg−1 and 2381 μmol min−1 mg−1 respectively. Penicillium occitanis Pol6  showed Km of 14.13 mg ml−1 and Vmax of 806.3 μmol min−1 mg−1 .
In summary, we show that xylanase produced by Aspergillus fumigatus R1 under submerged fermentation is novel enzyme. The xylanase was analyzed by mass-spectrometry, which indicated that the xylanase under study has never been reported earlier. The purified xylanase showed acidic optimum pH and broad range of pH stability mainly in the alkaline region. These characteristics coupled with the xylanase production using agricultural residue as the only source of carbon make the enzyme suitable for many applications, especially in food industry and a probable alternative to the use of chemicals in bioleaching in paper and pulp industries. Taken together, we propose that the 3D model of xylanase established could be used for rational design of xylanase-specific inhibitors and to further understand the role of each fold of this protein to determine its function. Moreover, the predicted 3D structural model can be used as a basal structure for obtaining point mutation to improve the catalytic efficiency of this protein and may be to increase thermostability.
Birchwood xylan and dialysis tubing (M.W. cut off, 10 kDa; average flat width, 25 mm; capacity, 60 ml per foot) were purchased from Sigma (Sigma Aldrich Co. Ltd., Germany). Bio-Gel P-60 for gel filtration was purchased from Bio-Rad. Rice bran was obtained from a local mill. All other chemicals were of reagent grade unless otherwise mentioned and obtained from suppliers, HiMedia, Sigma, GeNei, Merck and SRL.
Fungal strain isolation and culture conditions
Aspergillus fumigatus R1 was isolated in our laboratory from the soil samples (Puducherry, India) collected from the decaying paper and wood material. The strain was identified by lactophenol cotton blue staining and 18S rRNA sequence analysis. For isolation of xylanase producer, the preculture medium (g/L) contained: Birchwood xylan, 1; KH2PO4, 0.5; MgSO4, 0.25; NH4Cl, 1; Yeast-extract, 0.1. 1 g Soil sample was inoculated into the broth medium and incubated in shaker for 5–7 days at 70 rpm at 37 0C. For xylanase production, the basal medium for the flask culture contained (g/L): Rice bran, 10; KH2PO4, 1; MgSO4.7H2O, 0.2; K2HPO4, 1.2; (NH4)2SO4, 5; CaCl2.2H2O, 0.03; FeSO4.7H2O, 0.01. Shake flasks, prepared each containing liquid basal medium, were inoculated with the fungal spores and incubated at 370C.70 rpm for 96 h.
Culture from pre-culture medium was spread inoculated on sterile xylan-agar medium and incubated at 370C. After 48 h, the plate was flooded with 0.1 % (w/v) Congo-red for 15 min followed by 10 min treatment with 1 M NaCl and subsequently applied with 5 % (v/v) acetic acid .
The 18S rRNA partial sequence of the isolate was compared with other closely related strains using BLAST and NCBI GeneBank data base. Alignment and the phylogenetic tree were constructed using MEGA 5.05 software . The neighbor-joining (NJ) tree of the Aspergillus fumigatus R1 xylanase was evaluated using 1000 bootstrap replications . The nucleotide sequence of the culture Aspergillus fumigatus R1 was submitted in NCBI GeneBank database with an accession number KJ001801.
Xylanase production kinetics
The fungal culture was grown in basal medium containing rice-bran as the only source of carbon at 370C for 24 to 96 h. At regular time lapse of 24 h, the samples were withdrawn under aseptic conditions to check the enzyme activity.
Xylanase activity was assayed according to Bailey et al. . The reaction mixture contained equal volumes of 1 % Birchwood xylan and the suitably diluted enzyme solution with 50 mM phosphate buffer (pH 7) at 370C for 30 min. The amount of reducing sugar liberated was determined by Dinitrosalicylic acid method using xylose as the standard. One unit of the xylanase activity was defined as the amount of enzyme that catalyzes the release of 1millimol of xylose equivalent in 1 min. under standard assay conditions.
The protein concentration was determined according to the method of Lowry et al.  using bovine serum albumin as a standard.
Purification of xylanase
All purification steps were performed at 16 0C unless otherwise mentioned. The crude extracellular xylanase was obtained by centrifuging the culture broth at 8000 rpm for 10 min. The crude supernatant was subjected to 30–55 % fractional ammonium salt precipitation, stirred for 2–3 h and centrifuged at 10000 rpm for 10 min at 4 0C. The pellet obtained was subjected to extensive dialysis for removal of salts. The dialyzed sample was further concentrated by lyophilisation. The lyophilized dialyzed sample, 1 ml, was applied to Bio-Gel P-60 gel-filtration column (5 cm × 52 cm). The column was equilibrated with 50 mM phosphate buffer (pH 7) and then washed with the same buffer. The packing of column was confirmed by constant flow rate of 0.2 ml/min. The void volume of the column was determined by blue dextran. Elution of proteins was carried out with the same buffer and A 280 was used to monitor the protein in the purification steps. Fractions showing high specific activity were checked for their purity by SDS-PAGE. Those showing homogeneity of fractions were pooled and concentrated by lyophilisation and their purity and homogeneity was checked by sodium dodecyl sulphate-polyacrylamide gel electrophoresis (SDS-PAGE).
SDS-polyacrylamide gel electrophoresis (SDS-PAGE) and activity analysis
SDS-PAGE was carried out with vertical slab unit (Genaxy scientific Pvt. Ltd., US) using 15 % (w/v) acrylamide in gels under denaturing conditions according to the method described by Laemmli . The gels were stained by silver staining method . The relative molecular mass (Mr) of the pure xylanase was determined by using protein molecular weight markers (GeNei PMW-L) as standard: (Ovalbumin 43 kDa, Carbonic Anhydrase 29 kDa, Trypsin Inhibitor 20 kDa, Lysozyme 14.3 kDa and Insulin 3.5 kDa). Activity analysis was performed under nondenaturating conditions using birchwood xylan with 1 % final concentration in 15 % (w/v) native PAGE according the method of Breccia et al.  after non-denaturing electrophoresis, the gel was incubated in acetate buffer pH 5 for 30 min and thereafter at 60 0C for 15 min in the same buffer. Subsequently, the gel was stained with 0.1 % Congo-red to localize the band of protein with xylanase activity.
The purified xylanase band from native-PAGE was cut and sent to National Centre for Biological Sciences (NCBS), Bangalore, India for mass-spectrometric analysis. The purified protein was subjected to in-gel digestion with additional reduction and alkylation to bring in better sequence coverage according to Shevchenko et al. . Digested peptides were reconstituted in 15 μL of the 0.1 % formic acid and 1 μL was injected onto agilent 1200 nano flow HPLC system in-line coupled though advion nanomate to orbitrap discovery. Peptides were separated on Agilent ZORBAX SB-C18 nano column by a gradient developed from 1 % [v/v] acetonitrile, 0.1 % [v/v] formic acid to 80 % [v/v] acetonitrile, 0.1 % [v/v] formic acid in water over 70 min at a flow rate of 300 nl/min. Full MS in a mass range between m/z 300–2,000 was performed in the Orbitrap mass analyzer with a resolution of 30,000 at m/z 400 and an AGC target of 5 × 105. The strongest five signals were selected for CID (collision induced dissociation)-MS/MS in the LTQ ion trap at a normalized collision energy of 35 % using an AGC target of 3 × 104 and two microscans. Dynamic exclusion was enabled with one count during 30 s and an exclusion period of 180 s. The exclusion mass width was set to 0.01. For protein/peptide identification MS/MS data was searched against both uniprot swissprot amino acid sequence database (non-redundant database with reviewed proteins) and Uniprot TrEMBL database (database with unreviewed proteins) downloaded in August 2013 using an in-house mascot server (version 2.4) through Proteome Discoverer 1.3 software. The search was set up for full tryptic peptides with a maximum of three missed cleavage sites. Carbamidomethyl on cysteine, and oxidized methionine were included as variable modifications. The precursor mass tolerance threshold was 10 ppm and the maximum fragment mass error was 0.8 Da. The significance threshold of the ion score was calculated based on a false discovery rate of < 1 % estimated by the peptide validator node of the Proteome Discoverer software. Minimum of two high confident peptides were used as a prerequisite to identify the proteins. 12.5fmoles of Standard BSA digest was analyzed at the beginning and end of sequence to check the performance of the instrument.
Threading of protein structure
Two programs were employed on the I-TASSER server  to predict the secondary structure of the target protein xylanase: position specific iterated prediction (PSI-PRED) and position specific iterated-BLAST (PSI-BLAST). Initially, PSI-BLAST was performed for alignment of the query sequence against a non-redundant sequence database and then the secondary structure of the protein was predicted using PSI-PRED on the basis of the sequences generated by multiple alignments of the sequence homologs. LOMETS  is a locally installed meta-threading program on I-TASSER server which combines ten state-of-art treading programs (MUSTER, FFAS-3D, SPARKS-X, HHSEARCH2, HHSEARCH1, Neff-PPAS, HHSEARCH, pGenTHREADER, wdPPAS and cdPPAS). The secondary structure of the query sequence predicted was then threaded through a library of PDB structures using LOMETS. The templates were ranked on the basis of sequence-based profile-profile alignment, structure-based profile- profile alignment, secondary structure match, backbone torsion angle match, solvent accessibility match and general hydrobhobicity scoring matrix. Top-ranked template hits were selected from each threading program for further analysis.
Assembly and Refinement of protein structure
Two programs were used Support Vector Machine SVMSEQ  and SPICKER  for assembly and refinement of the target protein structure. SPICKER is clustering program which performs clustering using representative set of decoy confirmations. Unified knowledge-based force field has three components which guide simulations of structural assembly: (1) PDB Cα/ side chain correlations , H-bonds and hydrobhobicity ; (2) spatial restraints of threading templates ; and (3) SVMSEQ sequence based contact predictions. There are two sets of simulations for determination of structural assembly. (1) Generation of initial structures for threading templates and (2) SPICKER generates cluster centroids from trajectories obtained in the first set of simulations. Finally, TM-align identifies the PDB structures which are structurally close to the cluster centroids .
REMO: Modeling and 3D structure ranking for protein function prediction
The program REMO  depends on construction of full atomic models by optimizing H-bond network from C-alpha traces. Secondary structure specific backbone isomer was constructed using the program PSI-PRED. Full atomic model of our protein (xylanase) was constructed with reduced modeling simulation. Cα and side-chain centers play as the main simulation points for construction of models in reducing representation. These simulation points must be around to obtain better global topology of full atomic model with regular bond angle and bond length. In current study, Cα-RMSD and TM-scores parameters were considered for the verification of estimation of movement of Cα atoms relative to the native structure. Cα, C, N and O basic protein backbone fragments were matched from a secondary structure-specific backbone isomer library of non-redundant isomers from high-resolution PDB structures. H-bonding, clash/break-amendment, CHARMM22  and I-TASSER were the driving parameters in REMO refinement protocol. CHARMM22 identifies the potential energy functions by mathematical equations and respective force field dimensions of biomolecules. The models constructed by I-TASSER were ranked on the basis of structural density in SPICKER clusters. Eventually, the 3D models which had the highest scores were compared with the proteins of known structure and function in the PDB. Structural analogs of query protein in enzyme commission (EC) numbers, substrate binding site libraries and active site residues were evaluated on the basis of global and local structural similarity.
Biochemical characterisation of Xylanase
Effect of pH on xylanase activity
The optimum pH for xylanase was determined under standard assay conditions with 1 % birchwood xylan dissolved in different 50 mM pH buffers, including citrate buffer (pH 4), acetate buffer (pH 5), sodium-phosphate buffer (pH 6–7), Tris buffer (pH 8), glycine- NaOH buffer (pH 9) and carbonate-bicarbonate buffer (pH 10). 0.786 μg of xylanase was incubated with 1 % xylan at different pH and the enzyme activity was carried out under standard assay conditions.
Effect of pH on xylanase stability
For the pH stability test, 0.786 μg pure enzyme was placed in the above buffers with different pH and incubated at room temperature for 1 h and the relative activity was then assayed as described above at its optimum pH.
Effect of temperature on xylanase activity
The optimum temperature for xylanase was determined under standard assay conditions, by incubating the enzyme at different temperatures (30–70 0C) in 50 mM acetate buffer pH 5 (optimum pH) for 30 min and the relative activity was then assayed as described above.
Effect of temperature on xylanase stability
To determine the temperature stability of xylanase, the enzyme was incubated at different temperatures (30–70 0C) in acetate buffer (pH 5) for 1 h. The enzyme samples were cooled on ice for 10 min and relative xylanase activities were measured by the standard assay protocol as described above.
Determination of kinetic parameters
Purified enzyme was incubated with various concentrations of birchwood xylan (2.5–40 mg/ml) in 50 mM acetate buffer (pH 5) at 50 0C for 30 min. Kinetic parameters, Km and Vmax were calculated by linear regression from Lineweaver-Burk plots.
Ethical approval and consent from patients/ participants
This work did not involve any studies on human volunteers.
collision induced dissociation
high performance liquid chromatography
Iterative Threading ASSEmbly Refinement
Position-Specific Iterative Basic Local Alignment Search Tool
sodium dodecyl sulphate-polyacrylamide gel electrophoresis
support vector machine
Template Modeling align
Prade RA. Xylanases: from biology to biotechnology. Biotechnol Genet Eng Rev. 1995;13:100–31.
Shallom D, Shoham Y. Microbial hemicellulases. Curr Opin Microbiol. 2003;6:219–28.
Kulkarni N, Shendye A, Rao M. Molecular and biotechnological aspects of xylanases. FEMS Microbiol Rev. 1999;23:411–56.
Beg QK, Kapoor M, Mahajan L, Hoondal GS. Microbial xylanases and their industrial applications: a review. Appl Microbiol Biotechnol. 2001;56:326–38.
Collins T, Gerday C, Feller G. Xylanases, xylanase families and extremophilic xylanases. FEMS Microbiol Rev. 2005;29:3–23.
Paës G, Berrin JG, Beaugrand J. GH11 xylanases: structure, function, properties, relationships and applications. Biotechnol Adv. 2012;30:564–92.
Botella C, Diaz A, Ory I, Webb C, Blandino A. Xylanase and pectinase production by Aspergillus awamori on grape pomace in solid state fermentation. Process Biochem. 2007;42:98–101.
Haltrich D, Nidetzky B, Kulbe KD, Steiner W, Zupancic S. Production of fungal xylanases. Bioresour Technol. 1996;58:137–61.
Wong KKY, Tan LUL, Saddler JN. Multiplicity of beta-1,4-xylanase in microorganisms: functions and applications. Microbiol Rev. 1988;52:305–17.
Junker VL, Apweiler R, Bairoch A. Representation of functional information in the SWISS-PROT data bank. Bioinformatics. 1999;12:1066–7.
Henrissat B, Coutinho PM. Classification of glycoside hydrolases and glycosyltransferases from hyperthermophiles. Methods Enzymol. 2001;330:183–201.
Gebler J, Gilkes NR, Claeyssens M, Wilson DB, Beguin P, Wakarchuk WW, et al. Stereoselective hydrolysis catalyzed by related beta-1,4-glucanases and beta-1,4-xylanases. J Biol Chem. 1992;267:12559–61.
Liu D, Li J, Zhao S. Secretome diversity and quantitative analysis of cellulolytic Aspergillus fumigatus Z5 in the presence of different carbon sources. Biotechnol biofuels. 2013;6:149.
Bajaj BK, Abbass M. Studies on an alkali-thermostable xylanase from Aspergillus fumigatus MA28. 3 Biotech. 2011;3:161–71.
Yang J, Yan R, Roy A, Xu D, Poisson J, Zhang Y. The I-TASSER Suite: Protein structure and function prediction. Nat Methods. 2015;12:7–8.
Pittelkow M, Tschapek B, Smits SH, Schmitt L, Bremer E. The crystal structure of the substrate-binding protein OpuBC from Bacillus subtilis in complex with choline. J Mol Biol. 2011;411:53–67.
Lu FX, Lu M, Lu ZX, Bie XM, Zhao HZ, Wang Y. Purification and characterization of xylanase from Aspergillus ficuum AF-98. Bioresour Technol. 2008;99:5938–41.
Ito K, Ogasawara H, Sugimoto T, Ishikawa T. Purification and properties of acid stable xylanases from Aspergillus kawachii. Biosci Biotechnol Biochem. 1992;56:547–50.
Fernandez-Espinar MT, Pinaga F, de Graaff L, Visser J, Ramon D, Valles S. Purification, characterization and regulation of the synthesis of an Aspergillus nidulans acidic xylanase. Appl Microbiol Biotechnol. 1994;42:555–62.
Amita RS, Shah RK, Madamwar D. Improvement of the quality of whole wheat bread by supplementation of xylanase from Aspergillus foetidus. Bioresour Technol. 2006;97:2047–53.
Fang HY, Chang SM, Lan CH, Fang TJ. Purification and characterization of xylanase from Aspergillus carneus M34 and its potential use in photoprotectant preparation. Process Biochem. 2008;43:49–55.
Fujimoto Z, Kuno A, Kaneko S, Yoshida S, Kobayashi H, Kusakabe I, et al. Crystal structure of Streptomyces olivaceoviridis E-86 1,4-β-xylanase containing xylan-binding domain. J Mol Biol. 2000;300:575–85.
Wang YR, Zhang HL, He YZ, Luo HY, Yao B. Characterization, gene cloning and expression of a novel xylanase XYNB from Streptomyces olivaceoviridis A1. Aquaculture. 2007;267:328–34.
Carmona EC, Fialho MB, Buchgnani ÉB, Coelho GD, Brocheto-Braga MR, Jorge JA. Production, purification and characterization of a minor form of xylanase from Aspergillus versicolor. Process Biochem. 2005;40:359–64.
Sandrim VC, Rizzattib ACS, Terenzib HF, Jorgeb JA, Milagresc AMF, Polizelib MLTM. Purification and biochemical characterization of two xylanases produced by Aspergillus caespitosus and their potential for kraft pulp bleaching. Process Biochem. 2005;40:1823–8.
Driss D, Bhiri F, Elleuch L, Bouly N, Stals I, Miled N, et al. Purification and properties of an extracellular acidophilic endo-1,4-β xylanase, naturally deleted in the “thumb”, from Penicillium occitanis Pol6. Process Biochem. 2011;46:1299–306.
Bakira U, Yavascaoglu S, Guvenc F, Ersayin A. An endo-β-1,4-xylanase from Rhizopus oryzae: production, partial purification and biochemical characterization. Enzym Microb Tech. 2001;29:328–34.
Ninawe S, Kapoor M, Kuhad RC. Purification and characterization of extracellular xylanase from Streptomyces cyaneus SN32. Bioresour Technol. 2008;99:1252–8.
Nair SG, Sindhu R, Shashidhar S. Purification and biochemical characterization of two xylanases from Aspergillus sydowii SBS 45. Appl Biochem Biotechnol. 2008;149:229–43.
Teather RM, Wood PJ. Use of Congo red polysaccharide interactions in enumeration and characterization of cellulolytic bacteria from the bovine rumen. Appl Environ Microbiol. 1982;43:777–80.
Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S. MEGA5: Molecular Evolutionary Genetics Analysis using Maximum Likelihood, Evolutionary Distance and Maximum Parsimony Methods. Mol Biol Evol. 2011;28:2731–9.
Saitou N, Nei M. The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987;4:406–25.
Bailey MJ, Biely P, Poutanen K. Interlaboratory testing of methods for assay of xylanase activity. J Biotechnol. 1992;23:257–70.
Lowry OH, Rousebrough NJ, Farr AL, Randall RJ. Protein measurement with the Folin phenol reagent. J Biol Chem. 1951;193:265–75.
Laemmli UK. Cleavage of structural proteins during assembly of head of bacteriophage T4. Nature. 1970;227:680–5.
Chevallet M, Luche S, Rabilloud T. Silver staining of proteins in Polyacrylamide gels. Nat Protoc. 2006;4:1852–8.
Breccia JD, Baigori MD, Castro GR, Sineriz F. Detection of endo-xylanase activities in electrophoretic gels with Congo Red Staining. Biotechnol Tech. 1995;9:145–8.
Shevchenko A, Tomas H, Havlis J, Olsen JV, Mann M. In-gel digestion for mass spectrometric characterization of proteins and proteomes. Nat Protoc. 2006;6:2856–60.
Wu S, Zhang Y. LOMETS: A local meta-threading-server for protein structure prediction. Nucleic Acids Res. 2007;35:3375–82.
Wu S, Zhang Y. A comprehensive assessment of sequence-based and template-based methods for protein contact prediction. Bioinformatics. 2008;24:924–31.
Zhang Y, Skolnick J. SPICKER: A clustering approach to identify near-native folds. J Comput Chem. 2004;25:865–71.
Holm L, Sander C. Database algorithm for protein backbone and side-chains coordinates from a C-alpha trace application to model building and detection of coordinates error. J Biol. 1991;218:183–94.
McDonald IK, Thornton JM. Satisfying hydrogen bonding potential in proteins. J Mol Biol. 1994;238:777–93.
Bennett MJ, Schlunegger MP, Eisenberg D. 3D domain swapping: a mechanism for oligomer assembly. Protein Sci. 1995;4:2455–68.
Zhang Y, Skolnick J. TM-align: A protein structure alignment algorithm based on TM-Score. Nucleic Acids Res. 2005;33:2302–9.
Li Y, Zhang Y. REMO: A new protocol to refine full atomic protein models from C-alpha traces by optimizing hydrogen-bonding networks. Proteins. 2009;76:665–76.
Buck M, Bouguet-Bonnet S, Pastor RW, MacKerell AD. Importance of the CMAP correction to the CHARMM22 protein force field: Dynamics of hen lysozyme. Biophys J. 2006;90:36–8.
Financial Assistance from UGC, India to Dr. Sharmili Jagtap is duly acknowledged. The authors thank Mass-spectrometry-Proteomics facility, C-CAMP, NCBS, Bangalore, India.
I confirm that I have read BioMed Central's guidance on competing interests and we have no financial or non-financial competing interests.
RAD: Performed the experiments of purification and characterization and carried out the structural modelling and contributed to article preparation. MKM: Identification of the isolate, phylogenetic tree analysis, sequence alignment. SKM: Production of enzyme, and proofreading of the manuscript. SJ: conception and Experimental design, Isolation of microorganism, preparation and editing of manuscript, analysis and interpretation of data. All authors read and approved the final manuscript