Expression, Docking, and Molecular Dynamics of Endo-[beta]-1,4-xylanase I Gene of Trichoderma virens in Pichia stipitis.
Production of alternate fuels is crucial for economic development and therefore second and third generation bioethanol production have become key areas of research in many countries [1, 2]. Lignocellulosic biomass consists of cellulose (40-60%), hemicellulose (20-40%), and lignin (10-30%). Much research has already been focused on biofuel generation using cellulose. However, for economically feasible production of ethanol, the hemicellulose of plant biomass also should be utilized in an effective manner [3, 4].
Hemicelluloses are linear or branched heteropolysaccharides comprising different types of sugar residues, namely, xylose, mannose, galactose, arabinose, or glucose [5, 6]. The major component in hemicellulose is xylan, a polymer of 1,4-linked [beta]-D-xylose units . Xylanases, namely, [beta]-1,4-xylosidases (XYL) and endo-xylanases (EXN), are responsible for depolymerization of xylan to xylose [8, 9]. Endo-[beta]-1,4-xylanase (EXN) depolymerizes xylan to shorter xylobiose and xylooligomers  while [beta]-1,4-xylosidases (XYL) hydrolyzes xylobiose and xylooligomers to D-xylose generated by the action of endo-xylanases (EXN). Isolation and purification of EXN enzymes from cellulolytic microorganisms have been strongly directed towards research and development due to its incomparable commercial applications such as in paper pulping and biobleaching and in food and animal feed industries [6, 11]. Among many different microbes, filamentous fungi including Trichoderma species have been identified as better xylanase producers and they have been used for more than 50 years in the production of industrialxylanase [12,13].
Pentose fermentation pathway is restricted to a few yeast species such as Pichia stipitis, Candida shehatae, and Pachysolen tannophilus [14-16]. Therefore, P. stipitis has been genetically modified to improve xylose like pentose fermentation [17,18]. The present study describes the characterization, cloning, and expression of EXNI from locally isolated Trichoderma virens species with the long-term objective of producing second generation biofuel from hemicellulose using recombinant Pichia stipitis. In this study, the pGAPZ[alpha] integrative vector consisting of yeast a-mating factor (MFa) signal peptide with a glyceraldehyde 3-phosphate dehydrogenase (GAP) promoter driven expression system was used.
A three-dimensional (3D) structure of EXNI was built using homology modeling with quality assessments for the understanding of the structural features and properties. A molecular docking study was carried out to investigate the docking of substrates to the catalytic site of the enzyme. Molecular dynamics (MD) simulations provide evidence for conformational changes in the molecular system and the time-dependent behavior of biomolecules and study how the enzyme functions are affected by xylan.
2. Materials and Methods
2.1. Culturing of Trichoderma virens and Xylanase Activity Assay. Locally isolated Trichoderma species was confirmed as Trichoderma virens using PCR based internal transcribed spacer (ITS) analysis . The isolated species was cultivated for seven days on potato dextrose agar (PDA) at 30[degrees]C. An equal amount of inocula was selected from the above grown fungal culture using a sterile cork borer (5 mm diameter) and inoculated into conical flasks (100 mL) containing 25 mL of Mandel's medium (MM) (constituents of MM include [(N[H.sub.4]).sub.2]S[O.sub.4] 1.4g/L, K[H.sub.2]P[O.sub.4] 2.0g/L, Ca[Cl.sub.2] 0.3 g/L, MgS[O.sub.4]-7[H.sub.2]O 0.3 g/L, peptone 1.0 g/L, FeS[O.sub.4]-7[H.sub.2]O 5.0 mg/L, ZnS[O.sub.4]-7[H.sub.2]O 1.4 mg/L, MnS[O.sub.4]-[H.sub.2]O 1.6 mg/L, and Co[Cl.sub.2] 2.0 mg/L, pH 5.5) supplemented with 2% beechwood xylan [20, 21]. Cultures were incubated at 30[degrees]C, in a rotary shaker at 150 rpm. The enzyme extracts were harvested by filtration at 24 hour intervals for nine days. They were centrifuged at 6200g at 4[degrees]C for 10 minutes to obtain a cell free enzyme extract and stored at -20[degrees]C. Thereafter, the enzyme extracts were freeze dried and 100 [micro]L of the crude enzyme extracts was tested for endoxylanase (EXN) activity  using 900 [micro]L of 1% soluble beechwood xylan as the substrate. The reaction mixture was incubated at 50[degrees]C in a water bath for 5 minutes. Three milliliters of 3,5-dinitrosalicylic acid (DNS) solution was added to each reaction mixture and placed in a boiling water bath for 15 minutes for color development. The color intensity is proportional to the amount of reducing sugar (xylose) produced in the reaction mixture. Absorbance was measured at 540 nm to determine enzyme activity with reference to the substrate and enzyme blanks [23, 24]. One unit (1U) of endoxylanase activity is defined as the amount of 1 [micro]mole of liberated hydrolysis product (equivalent to xylose) in 1 mL of enzyme volume per minute ([micro]mol [ml.sup.-1] [min.sup.-1]). Enzyme activity was calculated with reference to the xylose standard graph.
2.2. Cloning and Sequence Analysis of EXN1 Gene of Trichoderma virens. Gene specific primers were designed for PCR amplification of EXNI gene prior to cloning. Sequences homologous to EXNI were identified using the nucleotide blast search in NCBI. Open reading frame (ORF) of the consensus sequence of EXNI was determined using ORF finder (https://www.ncbi.nlm.nih.gov/orffinder/). The native signal sequence, predicted by the SignalP3.0 server (http://www.cbs .dtu.dk/services/SignalP/), was omitted in primer designing of EXNI to accommodate the a-mating factor signal sequence of the pGAPZ[alpha] vector.
Genomic DNA was extracted from T. virens using a simple extraction method  and EXN1 was amplified using the above designed primers. EXN1FP: 5'-ATCGT[GAATTCC.bar]AGACGATTGGCCCCGGCACT-3' and EXN1RP: 5'-TTGA[TTCTAGA.bar]TTAGCTYACGTTAATGTTGGCGTTACC- AGAGCT-3'. PCR cycling conditions were optimized (initial denaturation at 94[degrees] C for 2 min followed by 35 cycles of denaturation at 94[degrees]C for 30 sec, annealing at 65[degrees]C for 30 sec, and extension at 72[degrees]C for 30 sec. This was followed by a final extension at 72[degrees]C for 7 min). The PCR product (~700) bp was electrophoresed on a 0.8% agarose gel, purified, and ligated into pGEM-T vector (Promega, USA). Thereafter, it was transformed into E. coli JM109 competent cells (Promega, USA) using the heat-shock method . The transformants were spread on to low salt Luria-Bertani (LB) agar medium (1% tryptone, 0.5% yeast extract, 0.5% NaCl, and 1.5% agar, pH 7.5) containing ampicillin (100 [micro]g/mL), 5-bromo4-chloro-indolyl-D-galactoside (0.2 mM), and isopropylthio-[beta]-D-galactopyranoside (40 [micro]g/mL). Plasmid extractions were carried out on selected white colonies and custom sequenced (Macrogen, Korea). A recombinant EXNI clone was designated as pGEM-T/gEXNI. After sequence confirmation, pGEM-T/gEXNI was digested with EcoRI and XbaI restriction enzymes and the insert purified and ligated into pGAPZ[alpha] vector. Ligated products were transformed into E. coli JM109 and transformants were selected in zeocin (25 [micro]g/mL) containing low salt LB agar plates. The resulting clone was designated as pGAPZ[alpha]/gEXNI.
2.3. Transformation into Pichia stipitis. The pGAPZ[alpha]/gEXNI vector construct was linearized using BglII and purified. The plasmid concentration was determined and adjusted to ~1 [micro]g/[micro]L. Pichia stipitis (ATCC 58376) was inoculated into 0.5 mL YPD broth (1% yeast extract, 2% peptone, and 2% glucose) in a 1.5 mL microcentrifuge tube and incubated at 30[degrees]C overnight in a rotary shaker at 150 rpm. A volume of 500 [micro]L from the above grown culture was inoculated into 50 mL broth in a 250 mL conical flask and incubated at 30[degrees]C at 150 rpm until the O[D.sub.600] reached 1.4. Yeast electrocompetent cells were prepared according to the procedure given in pGAPZ[alpha] vector manual (Invitrogen, USA). A volume of 80 [micro]L P stipitis competent cells was mixed with 5-10 pg of linearized pGAPZ[alpha]/gEXNI plasmid DNA. The mixture was subjected to electroporation under the optimized conditions of 1.5 Kw, 200 mA, and 25 [micro]F (pulse time of 5 ms) in a 0.2 cm electroporation cuvette. The resulting transformation mixture was spread on to YPDS (1% yeast extract, 2% peptone, 2% glucose, 2% agar, and 1M sorbitol) plates with 100 [micro]g/mL zeocin as the selection marker. The plates were then incubated for 3 days at 30[degrees]C to obtain positive transformants. Twenty yeast colonies were selected and streaked on fresh YPDS plates containing zeocin (100 [micro]g/mL).
2.4. Screening of Recombinant Yeast. Colony PCR was performed for the above selected colonies to confirm the presence of the integrated EXNI gene in the P. stipitis genome. Native P. stipitis and pGAPZ[alpha]/gEXNI plasmid DNA were used as the negative and positive controls, respectively. PCR amplified product was subjected to agarose (0.8%) gel electrophoresis. A putative clone designated Y-pGAPZ[alpha]/gEXNI was sequenced. Thereafter it was cultured in YPD broth with 100 [micro]g/mL zeocin as the selection marker and incubated overnight. These cultures were used for expression analysis and long-term storage purposes in 50% glycerol at -80[degrees]C.
2.5. SDS-PAGE and Expression Analysis of EXNI in Recombinant P. stipitis. Y-pGAPZ[alpha]/gEXNI recombinant P. stipitis was inoculated into beechwood xylan broth culture and incubated overnight at 30[degrees]C in a rotary shaker at 200 rpm. Native P. stipitis was used as the control. The enzyme was harvested by centrifugation of culture medium at 12000 rpm for 2 min at 4[degrees]C. The enzyme extract was concentrated by freeze drying. SDS-PAGE analysis  was performed and the results were analyzed by comparing with native P. stipitis.
Enzyme activity assay was carried out on the YpGAPZ[alpha]/gEXN recombinant P. stipitis using native P. stipitis as the control. They were separately inoculated from YPD streak plates into 0.5 mL of YPD broth cultures in 1.5 mL microcentrifuge tubes and incubated for 3 days in a rotary shaker at 200 rpm at 30[degrees]C. The enzyme harvest was freeze dried and a volume of 900 [micro]L of solubilized 1% beechwood xylan was treated with 100 [micro]L of above enzyme extract. The standard assay procedure was performed and recombinant EXNI enzyme activity was quantitatively determined according to the method described in Section 2.1. Both enzyme and substrate controls were maintained throughout the assay procedure.
2.6. Structure Prediction and Evaluation of EXNI Protein. The tertiary structure of EXNI protein was generated with the aid of MODELLER (version 9.13) program [28, 29]. The most similar X-ray crystallographic structures were identified in RSCB PDB protein databank [30, 31] using the Blast Protein tool  against the deduced amino acid sequence of EXNI as the target. Multiple sequence alignment method was used for homology modeling and the generated model was based on the templates of PDB IDs: 2JIC, 4S2H, and 4HKW for EXNI in the RSCB PDB protein databank. The constructed model was validated using structure validation tools of VERIFY3D , PROCHECK , and ERRAT  to analyze the compatibility of the model with its amino acid sequence, to verify the geometrical and stereochemical constraints of the model and to determine the overall quality factor respectively. The COACH server (https://zhanglab.ccmb.med.umich.edu/COACH) [36, 37] was used to identify the binding domain of the above generated EXNI model.
2.7. Molecular Docking. The 3D structure of the oligoxylose ligand was constructed and geometrically optimized before docking, with 6-31g ** basis set using Gaussian 09 (linux version) software . Protein and above optimized ligand were prepared using UCSF chimera , and the ligand docked into the active site of the model structure of EXNI using DOCK6 software [40,41] using flexible docking method. The binding strength of the protein and the ligand was ranked using the grid score energies .
2.8. Molecular Dynamics Simulation. Molecular dynamics simulation was performed in two phases. In the first phase, the docked protein-ligand complex obtained from above docking process with lowest binding energy was selected as the initial configuration for the molecular dynamics (MD) simulations using the GROMACS v4.6.5 . The GROMOS54a7 united atom force field was assigned for the model protein. The force field parameters of the ligand were obtained from PRODRG server . The protein-ligand complex was inserted in the center of a box with 9 x 9 x 9n[m.sup.3] volume. Electroneutrality of the system was maintained by adding [Na.sup.+] ions to the box and solvated with SPC/E water molecules . Electrostatic interactions were modeled by particle mesh Ewald (PME) with a short-range cutoff of 1.2 nm . Berendsen's weak coupling algorithm was employed to maintain the temperature and pressure of the system at 300 K and 1 bar . Using LINCS algorithm  all bonds were constrained at their equilibrium distances while allowing other internal motions of bending and torsion during molecular dynamics simulation. The system was subjected to 2000 steps of energy minimization with steepest descent algorithm followed by 200 ps long MD simulation to equilibrate the simulation system. After the equilibration step, 15 ns MD simulation was carried out using a desktop computer with Intel[R] Core[TM] i7-950 Processor. Configurations of the system at every 2 ps intervals were stored for further analysis. At the end of the simulation, the noncovalent interactions between ligand and the protein were analyzed by the LigPlot+v.14.5 software . The same protocol was used in the second phase, to simulate the bare protein, obtained from protein homology modeling process.
The first phase molecular dynamics simulation was carried out to investigate the stability of the docked protein-ligand complex and the second phase molecular dynamics simulation was employed to investigate the structural comparison between docked complex and bare protein during the simulations in the aqueous medium.
3. Results and Discussion
The maximum EXN activity of Trichoderma virens was 114.2 U/mL on day seven. PCR amplification of EXNI gene using the genomic DNA of Trichoderma virens yielded a 689 bp fragment (GeneBank: KJ882380). Sequence analysis revealed a single intron of 115 bp (from 274 bp to 387bp). The coding region consists of a 573 bp open reading frame that encodes 191 amino acids. A sequence similarity search (NCBI) indicated that the above amplified EXNI was 99% identical to the DNA sequence of endo-1,4-beta-xylanase I gene of Trichoderma asperellum strain T-1 (KM277356.1) and 99% identical to the mRNA sequence of endo-1,4-beta-xylanase of Trichoderma viride strain YNUCC0183 (AY320048.1). Coding sequence of the EXN I was translated and was 100% similar to the Trichoderma viride amino acid sequence (AAP83925.1). Further, the amino acid residues ranging from 9 to 189 were identified to belong to the glycosyl hydrolases family 11 by the InterProScan server (EMBL). Many characterized xylanases are classified under this family. Three potential N-glycosylation sites at positions 61, 97, and 188 of the amino acid sequence were identified using NetOGlyc 4.0 Server and GlycoEP server. Search for O-glycosylation sites using the NetNGlyc 1.0 Server did not reveal any O-glycosylation sites. The theoretical molecular weight of the EXNI protein was calculated as 20.7 kDA and isoelectric pH was 8.75. SDS-PAGE analysis confirmed that the EXNI recombinant enzyme expressed by P. stipitis was ~24kDA (Figure 1).
The EXNI enzyme activity expressed by Y-pGAPZ[alpha]/gEXNI P stipitis clone was 31.7 U/mL indicating the successful approach of genetic engineering in the heterologous extracellular expression of EXNI using the (GAP) promoter driven expression system. However, the EXNI activity of recombinant Y-pGAPZ[alpha]/gEXNI was 3.6 times less than the native T. virens. It is reported that yeast species hyperglycosylated secretory recombinant proteins [50, 51]. The experimental molecular weight of the EXNI was comparatively higher (~24kDA) than the molecular weight of the EXNI deduced from the amino acid sequence. There is thus a~ 3 kDA difference in the experimental and the theoretical molecular mass of EXNI. A possible explanation for the low enzyme activity could be due to hyperglycosylation of recombinant EXNI at N-glycosylation sites by P. stipitis compared to the native fungus. It has been reported that N-linked hyperglycosylation can have a significant effect on protein expression and function [52-54]. According to these studies N-linked hyperglycosylation usually occurs at the sequence Asn-Xaa-Ser/Thr (Xaa, any amino acid). As mentioned above, there are three potential glycosylation sites at 61 (Asn-Phe-Ser), 97 (Asn-Pro-Ser), and 188 (Asn-ValSer) which can be hyperglycosylated to support the above possibility. Recent mutational studies on glycosylation sites have confirmed that the presence of a proline residue between Asn and Ser/Thr will inhibit N-linked hyperglycosylation [55, 56]. At position 97 there is a proline residue that may inhibit hyperglycosylation. However, the other two positions are possible sites for hyperglycosylation.
The recombinant EXNI activity (31.7 U/mL) observed in the present study was considerably lower than the heterologous expression of different xylanases reported for Pichia pastoris in several recent studies where the highest activity ranged from 120U/ml to 746U/ml [8, 57-59]. These high enzyme activities are attributed to the low hyperglycosylation of recombinant proteins expressed by P.pastoris compared to expression of heterologous secretory proteins in other yeast species as the host [8, 60]. However, P. pastoris does not harbour pentose utilization enzymes needed for fermentation of xylose to ethanol. It is reported that the pentose utilization ability is much more efficient in P. stipitis compared to all other Pichia species [61-63]. To our knowledge, this is the only recorded study, for the heterologous expression of Trichoderma EXNI in P. stipitis. The long-term goal of the present study is to express both xylanases, EXN andXYL , in a yeast species to ferment xylan, the major component of hemicellulose in biomass. To achieve this goal the expression of recombinant EXNI in P. stipitis should be much higher than the present expression. This could be achieved by (i) introducing proline residues between Asn and Ser/Th or by introduction of specific amino acids, such as Trp, Asp, Glu, or Leu, residues between Asn and Ser/Thr at 61 and 188 positions of EXNI to convert the sequon to a poor oligosaccharide acceptor, to omit the hyperglycosylation. However, at the same time the integrity of the tertiary structure and the catalytic activity of the protein should be maintained with the aid of in silico strategies including homology modeling and molecular dynamics simulation techniques [65, 66]. (ii) According to the in silico strategies of this study the active site residues of EXNI were identified at positions 126-131. Random mutagenesis could be used to increase the acidity or the basicity of the above amino acid residues which can facilitate the catalytic activity under "retention" and "inversion" mechanisms described for the glycosyl hydrolases [67, 68]. (iii) An alternative approach would be to introduce genes encoding xylose utilizing pathway enzymes into P.pastoris by advanced genomic shuffling [69,70]. This would involve recombination of entire genome of Y-pGAPZ[alpha]/gEXNI recombinant P. stipitis with that of P. pastoris.
Five probable models were obtained in homology modeling from the MODELLER 9.13 software and they were ranked according to their normalized Discrete Optimized Protein Energy (zDOPE) and GA341 score. The model comprising the best scores was selected as the theoretical model for EXNI protein.
According to the characterization of the EXNI model by DSSP program  the secondary structure of EXNI is composed of 1 [alpha]-helix and 14 [beta]-sheets as given below.
Alpha Helix. VAL152-LEU162
[beta]-Sheets. (GLY6-ASN10), (TYR13-ASN19), (VAL25-ASN29), (SER34-TRP39), (ASN44-TRP51), (VAL59-ASN69), (ASN71-SER80), (ILE85-PHE93), (THR103-SER110), (SER113-ARG122), (PH134-ARG141), (GLY148-THR151), (GLN172-TYR179), (SER182-SER190).
Results from the PROCHECK analysis of EXNI are given in Table 1 and the Ramachandran plot generated by the same program is represented in Figure 2. According to the statistical score of the Ramachandran plot none of amino acids are in the disallowed region.
VERIFY3D profile of EXNI shows that all the residues have an averaged 3D-1D score greater than zero (Figure 3). Furthermore, to pass the VERIFY3D test, it is essential to show that at least 80% of the amino acids have scored more than 0.2  in the 3D/1D profile. In this study, the results indicate that 96.32% of the residues had an averaged 3D-1D score of greater than or equal to 0.2. Moreover, the ERRAT program evaluated the overall quality factor as 86.81 for the modeled 3D structure of EXNI. Above evaluations concluded a highly reliable 3D structure of EXNI.
The EXNI-xylan complex with the lowest binding energy was selected for the molecular docking method. Initially the binding residues of EXNI (TRP18, VAL46, TYR77, TRP79, GLU86, T YR88, ARG122, PRO126, SER127, PHE134, GLN136, TYR171, ASN44, TYR77, ILE173, and SER16) were identified using I-TASSER-COACH server. Above binding residues were interestingly matched with the results obtained from molecular docking. The recorded best grid score for EXNI-xylan was -80.47 kcal mol-1. Above grid score is a summation of van der Waals dispersive and electrostatic interaction energy which indicates the approximate binding energy of the ligand. Figure 4 represents the best protein-ligand complex from molecular docking.
The general catalytic mechanisms of glycosyl hydrolases follow either retaining or inverting mechanisms with the assistance of acidic amino acid residues in the active site [67, 68]. According to the docking results the presence of glutamic and aspartic acids in the active site may perform any of the above mechanisms.
Two MD simulations of 15 ns each were carried out for the protein-ligand complex and the bare protein in aqueous medium with 21923 SPC/E water molecules. The noncovalent interaction (H bond) of the final configuration (after 15 ns) of protein-ligand complex identified from LigPlot+v.145 software is presented in Figures 5 and 6. The LigPlot analysis of the protein structure indicates that the ligand forms three strong hydrogen bonds with SER16 (3.39 A), ILE173 (3.47 [Angstrom]), and TYR77 (2.94 [Angstrom]). Detailed information is presented in Table 2. The stability of all these H bonds was studied using gdist tool in the GROMACS program.
Throughout the simulation time the distance between the Centers of mass of the two groups of atoms which was involved in H bond formation was maintained nearly at a constant value confirming the continuance, stability, and effectiveness of the H bonding.
The top pane of Figure 7 compares the root mean square deviation (RMSD) of the EXNI backbone of bare protein and the protein-ligand (EXNI-xylan) complex system separately. Both systems indicate stable structures by maintaining a steady RMSD of about 0.2 nm with the simulation time. The bottom pane of Figure 7 gives the variation of radius of gyration (Rg) as a function of simulation time which indicates the compactness of the protein. As seen in the figure, Rg of both systems were maintained approximately the same value with minimum deviations. These results suggest that EXNI maintains its tertiary structure even after forming a complex with the ligand (xylan).
The root mean square fluctuations (RMSF) of amino acid residues of EXNI (Figure 8) represent the stability of 3D structures for both systems of bare protein and protein-ligand complex. Most of the fluctuations are concentrated in 96, 100, and 128 amino acid residues in both systems. Significant fluctuations of amino acid residues have been observed in 126-131 regions in the free protein and in 58 and 98 positions in the protein-ligand complex. Further, none of the high fluctuating residues of the protein of the complex were in the predicted active site. Thus it can be postulated that EXNI can initiate its catalytic activity with xylan (ligand) via either retention or inversion mechanisms.
Five probable models were obtained in homology modeling from the MODELLER 9.13 software and they were ranked according to their normalized Discrete Optimized Protein Energy (zDOPE) and GA341 score. The model comprising the best scores was selected as the theoretical model for EXNI protein.
EXNI gene was successfully characterized, cloned, and expressed from locally isolated T. virens in P stipitis. Three potential N-linked glycosylation sites were identified in EXNI. The increase in molecular weight observed from SDS-PAGE analysis of EXNI compared to the theoretically calculated value can be attributed to hyperglycosylation. Recombinant P. stipitis containing EXNI has the potential to degrade xylan the major component of hemicellulose fraction in plant biomass. Moreover it can be utilized to produce bioethanol via the combinatory simultaneous action with XYL for the hydrolysis and fermentation of xylan into ethanol.
Molecular docking studies provided major information on ligand binding domain and the active site residues of the protein. MD simulation results indicated almost similar and higher stability of EXNI in the form of protein-ligand complex and also in the free form of the protein. Enzyme-substrate association was steadily maintained by three hydrogen bonds in the ligand binding domain of EXNI. The predicted model of EXNI is realistic and therefore it provides more information of EXNI-xylan interaction for the designing of mutagenic experiments aimed at improving the catalytic action.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
The authors are grateful to National Science Foundation (NSF), Sri Lanka, for funding this research (Grant no. RG/2012/BT/02).
 IEA, "Technology Roadmap Biofuels for Transport," International Energy Agency: France, pp. 19-145, 2011.
 J. Goldemberg, "Biomass and energy," Quimica Nova, vol. 32, no. 3, pp. 582-587, 2009.
 Y. Sun and J. Cheng, "Hydrolysis of lignocellulosic materials for ethanol production: a review," Bioresource Technology, vol. 83, no. 1, pp. 1-11, 2002.
 O. J. Sanchez and C. A. Cardona, "Trends in biotechnological production of fuel ethanol from different feedstocks," Bioresource Technology, vol. 99, no. 13, pp. 5270-5295, 2008.
 M. P. Coughlan and G. P. Hazlewood, "[beta]-1,4-D-Xylan-degrading enzyme systems: biochemistry, molecular biology and applications," Biotechnology and Applied Biochemistry, vol. 17, pp. 259-289, 1993.
 Q. K. Beg, M. Kapoor, L. Mahajan, and G. S. Hoondal, "Microbial xylanases and their industrial applications: a review," Applied Microbiology and Biotechnology, vol. 56, no. 3-4, pp. 326-338, 2001.
 B. Lindberg, K. G. Rossell, and S. Svensson, "Position of O-acetyl groups in birch xylan," Acta Chemica Scandinavica, vol. 76, pp. 30-32, 1973.
 J. He, B. Yu, K. Zhang, X. Ding, and D. Chen, "Expression of endo-1,4-beta-xylanase from Trichoderma reesei in Pichia pastoris and functional characterization of the produced enzyme," BMC Biotechnology, vol. 9, article no. 56, 2009.
 M. Sharma and A. Kumar, "Xylanases: an overview," British Biotechnology Journal, vol. 3, no. 1, pp. 1-28, 2013.
 T. Collins, C. Gerday, and G. Feller, "Xylanases, xylanase families and extremophilic xylanases," FEMS Microbiology Reviews, vol. 29, no. 1, pp. 3-23, 2005.
 S. Subramaniyan and P. Prema, "Biotechnology of microbial xylanases: enzymology, molecular biology, and application," Critical Reviews in Biotechnology, vol. 22, no. 1, pp. 33-64, 2002.
 H. Dalboge, "Expression cloning of fungal enzyme genes: a novel approach for efficient isolation of enzyme genes of industrial relevance," FEMS Microbiology Reviews, vol. 21, no. 1, pp. 29-42, 1997
 M. Naveen and K. G. Siddalingeshwara, "Molecular confirmation, identification and influence of carbon source for the production of xylanase from Penicillium citrinum," Journal of Drug Delivery and Therapeutics, vol. 5, no. 6, pp. 63-67, 2015.
 J. C. du Preez, M. Bosch, and B. A. Prior, "The fermentation of hexose and pentose sugars by Candida shehatae and Pichia stipitis," Applied Microbiology and Biotechnology, vol. 23, no. 34, pp. 228-233, 1986.
 J. C. du Preez, M. Bosch, and B. A. Prior, "Temperature profiles of growth and ethanol tolerance of the xylose-fermenting yeasts Candida shehatae and Pichia stipitis," Applied Microbiology and Biotechnology, vol. 25, no. 6, pp. 521-525, 1987
 P Basaran, N. Basaran, and Y. D. Hang, "Isolation and characterization of Pichia stipitis mutants with enhanced xylanase activity," World Journal of Microbiology and Biotechnology, vol. 16, no. 6, pp. 545-550, 2000.
 J. Shi, M. Zhang, L. Zhang, P Wang, L. Jiang, and H. Deng, "Xylose-fermenting Pichia stipitis by genome shuffling for improved ethanol production," Microbial Biotechnology, vol. 7, no. 2, pp. 90-99, 2014.
 C.-F. Huang, T-H. Lin, G.-L. Guo, and W.-S. Hwang, "Enhanced ethanol production by fermentation of rice straw hydrolysate without detoxification using a newly adapted strain of Pichia stipitis," Bioresource Technology, vol. 100, no. 17, pp. 3914-3920, 2009.
 N. Sadfi-Zouaoui, S. Chaabani, M. Rouaissi, A. Hedi, M. R. Hajlaoui, and A. Boudabous, "Analysis of the diversity of Trichoderma spp. in soil horizons using digested ITS regions," Annals of Microbiology, vol. 59, no. 3, pp. 459-463, 2009.
 M. Mandels and D. Sternburg, "Recent advances in cellular technology," Journal of Fermentation Technology, vol. 54, pp. 267-286, 1976.
 N. H. Bhavsar, B. V. Raol, G. G. Raol, and P R. Bhatt, "Isolation, Screening and Optimisation of Xylanase Producing Fungi from Compost Pit," Imperial Journal of Interdisciplinary Research, vol. 2, no. 11, pp. 782-788, 2016.
 M. J. Bailey, P. Biely, and K. Poutanen, "Interlaboratory testing of methods for assay of xylanase activity," Journal of Biotechnology, vol. 23, no. 3, pp. 257-270, 1992.
 M. Nummi, J. M. Perrin, M. L. Niku-Paavola, and T. M. Enari, "Measurement of xylanase activity with insoluble xylan substrate," Biochemical Journal, vol. 226, no. 2, pp. 617-620, 1985.
 T. K. Ghose, "Measurement of cellulase activities," Pure and Applied Chemistry, vol. 59, no. 2, pp. 257-268, 1987
 T. H. Al-Samarrai and J. Schmid, "A simple method for extraction of fungal genomic DNA," Letters in Applied Microbiology, vol. 30, no. 1, pp. 53-56, 2000.
 J. Sambrook, E. F. Fritsch, and T. Maniatis, Molecular Cloning: A Laboratory Manual, Cold Spring Harbor Laboratory Press, New York, USA, 1989.
 J. Sambrook and D. Russell, Molecular Cloning: A Laboratory Manual, Cold Spring Harbor Laboratory Press, New York, USA, 3rd edition, 2001.
 B. Webb and A. Sali, "Comparative protein structure modeling using MODELLER," in Current Protocols in Bioinformatics, pp. 5.6.1-5.6.32, John Wiley & Sons, New York, NY, USA, 2014.
 M. A. Marti-Renom, A. C. Stuart, A. Fiser, R. Sanchez, F. Melo, and A. Sali, "Comparative protein structure modeling of genes and genomes," Annual Review of Biophysics and Biomolecular Structure, vol. 29, pp. 291-325, 2000.
 H. M. Berman, J. Westbrook, Z. Feng et al., "The protein data bank," Nucleic Acids Research, vol. 28, no. 1, pp. 235-242, 2000.
 H. Berman, K. Henrick, and H. Nakamura, "Announcing the worldwide Protein Data Bank," Nature Structural Biology, vol. 10, no. 12, p. 980, 2003.
 T. L. Madden, R. L. Tatusov, and J. Zhang, "Applications of network BLAST server," Methods in Enzymology, vol. 266, pp. 131-141, 1996.
 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.
 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.
 C. Colovos and T. O. Yeates, "ERRAT: an empirical atom-based method for validating protein structures," Protein Science, vol. 2, pp. 1511-1519, 1993.
 J. Yang, A. Roy, and Y. Zhang, "Protein-ligand binding site recognition using complementary binding-specific substructure comparison and sequence profile alignment," Bioinformatics, vol. 29, no. 20, pp. 2588-2595, 2013.
 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.
 M. J. E. A. Frisch, Gaussian 09 Revision A 02, Gaussian Inc, Wallingford, Conn, USA, 2009.
 E. F. Pettersen, T. D. Goddard, C. C. Huang et al., "UCSF Chimera--a visualization system for exploratory research and analysis," Journal of Computational Chemistry, vol. 25, no. 13, pp. 1605-1612, 2004.
 S. R. Brozell, S. Mukherjee, T. E. Balius, D. R. Roe, D. A. Case, and R. C. Rizzo, "Evaluation of DOCK 6 as a pose generation and database enrichment tool," Journal of Computer-Aided Molecular Design, vol. 26, no. 6, pp. 749-773, 2012.
 W. J. Allen, T. E. Balius, S. Mukherjee et al., "DOCK 6: impact of new features and current docking performance," Journal of Computational Chemistry, vol. 36, no. 15, pp. 1132-1156, 2015.
 A. C. Wallace, R. A. Laskowski, and J. M. Thornton, "LIGPLOT: a program to generate schematic diagrams of protein-ligand interactions," Protein Engineering, vol. 8, no. 2, pp. 127-134, 1995.
 H. J. C. Berendsen, D. van der Spoel, and R. van Drunen, "GROMACS: a message-passing parallel molecular dynamics implementation," Computer Physics Communications, vol. 91, no. 1-3, pp. 43-56, 1995.
 A. W. Schuttelkopf and D. M. F. van Aalten, "PRODRG: a tool for high-throughput crystallography of protein-ligand complexes," Acta Crystallographica, Section D: Biological Crystallography, vol. 60, no. 8, pp. 1355-1363, 2004.
 H. J. C. Berendsen, J. R. Grigera, and T. P Straatsma, "The missing term in effective pair potentials," Journal of Physical Chemistry, vol. 91, no. 24, pp. 6269-6271, 1987
 U. Essmann, L. Perera, M. L. Berkowitz, T. Darden, H. Lee, and L. G. Pedersen, "A smooth particle mesh Ewald method," The Journal of Chemical Physics, vol. 103, no. 19, pp. 8577-8593, 1995.
 H. J. C. Berendsen, J. P. M. Postma, W. F. Van Gunsteren, A. Dinola, and J. R. Haak, "Molecular dynamics with coupling to an external bath," The Journal of Chemical Physics, vol. 81, no. 8, pp. 3684-3690, 1984.
 B. Hess, "P-LINCS: a parallel linear constraint solver for molecular simulation," Journal of Chemical Theory and Computation, vol. 4, no. 1, pp. 116-122, 2008.
 R. A. Laskowski and M. B. Swindells, "LigPlot+: multiple ligand-protein interaction diagrams for drug discovery," Journal of Chemical Information and Modeling, vol. 51, no. 10, pp. 2778-2786, 2011.
 M. A. Romanos, C. A. Scorer, and J. J. Clare, "Foreign gene expression in yeast: a review," Yeast, vol. 8, no. 6, pp. 423-488, 1992.
 R. Conde, R. Cueva, G. Pablo, J. Polaina, and G. Larriba, "A search for hyperglycosylation signals in yeast glycoproteins," Journal of Biological Chemistry, vol. 279, no. 42, pp. 43789-43798, 2004.
 W. H. Shakin-Eshleman and S. L. Spitalnik, "Role of Individual iV-Linked Glycans in Cell Surface Expression and Secretion of Glycoproteins," Trends in Glycoscience and Glycotechnology, vol. 5, no. 25, pp. 355-368, 1993.
 P J. Gallagher, J. M. Henneberry, J. F. Sambrook, and M.J. H. Gething, "Glycosylation requirements for intracellular transport and function of the hemagglutinin of influenza virus," Journal of Virology, vol. 66, no. 12, pp. 7136-7145, 1992.
 C. E. Machamer and J. K. Rose, "Influence of new glycosylation sites on expression of the vesicular stomatitis virus G protein at the plasma membrane," Journal of Biological Chemistry, vol. 263, no. 12, pp. 5948-5954, 1988.
 E. Bieberich, "Synthesis, Processing, and Function of N-glycans in N-glycoproteins," in Glycobiology of the Nervous System, vol. 9 of Advances in Neurobiology, pp. 47-70, Springer New York, New York, NY, USA, 2014.
 Y. Gavel and G. V. Heijne, "Sequence differences between glycosylated and non-glycosylated Asn-X-Thr/ser acceptor sites: Implications for protein engineering," Protein Engineering, Design and Selection, vol. 3, no. 5, pp. 433-442, 1990.
 V. Ruanglek, R. Sriprang, N. Ratanaphan et al., "Cloning, expression, characterization, and high cell-density production of recombinant endo-1,4-[beta]-xylanase from Aspergillus niger in Pichia pastoris," Enzyme and Microbial Technology, vol. 41, no. 1-2, pp. 19-25, 2007.
 Y. Cao, J. Qiao, Y. Li, and W. Lu, "De novo synthesis, constitutive expression of Aspergillus sulphureus [beta]-xylanase gene in Pichia pastoris and partial enzymic characterization," Applied Microbiology and Biotechnology, vol. 76, no. 3, pp. 579-585, 2007
 M. C. T. Damaso, M. S. Almeida, E. Kurtenbach et al., "Optimized Expression of a Thermostable Xylanase from Thermomyces lanuginosus in Pichia pastoris," Applied and Environmental Microbiology, vol. 69, no. 10, pp. 6064-6072, 2003.
 P. Fickers, "Pichia pastoris: a workhorse for recombinant protein production," Current Research in Microbiology and Biotechnology, vol. 2, no. 3, pp. 354-363, 2014.
 T. W. Jeffries, I. V. Grigoriev, J. Grimwood et al., "Genome sequence of the lignocellulose-bioconverting and xylose-fermenting yeast Pichia stipitis," Nature Biotechnology, vol. 25, no. 3, pp. 319-326, 2007.
 F. K. Agbogbo and G. Coward-Kelly, "Cellulosic ethanol production using the naturally occurring xylose-fermenting yeast, Pichia stipitis," Biotechnology Letters, vol. 30, no. 9, pp. 1515-1524, 2008.
 S. Fernandes and P. G. Murray, "Metabolic engineering for improved microbial pentose fermentation," Bioengineered Bugs, vol. 01, no. 06, pp. 424-428, 2010.
 G. H. I. M. Wickramasinghe, N. V. Chandrasekharan, S. Rathnayake, S. Weerasinghe, R. L. C. Wijesundera, and W. S. S. Wijesundera, "[beta]-1,4-xylosidase (XYL) Containing Recombinant Pichia stipitis: Cloning, Expression and Homology Modeling," in Proceedings of the International Conference on Environmental Engineering, Transport and Renewable Energy (EETRE-2016), 2016.
 X. Ji, W. Wang, Y. Zheng, J. Hao, and M. Sun, "Homology modeling and molecular dynamics simulation studies of a marine alkaline protease," Bioinformatics and Biology Insights, vol. 6, pp. 255-263, 2012.
 H. A. Carlson, "Protein flexibility and drug design: How to hit a moving target," Current Opinion in Chemical Biology, vol. 6, no. 4, pp. 447-452, 2002.
 G. Davies and B. Henrissat, "Structures and mechanisms of glycosyl hydrolases," Structure, vol. 3, no. 9, pp. 853-859, 1995.
 D. E. Koshland, "Stereochemistry and the mechanism of enzymatic reactions," Biological Reviews of the Cambridge Philosophical Society, vol. 28, no. 4, pp. 416-436, 1953.
 W. Zhang and A. Geng, "Improved ethanol production by a xylose-fermenting recombinant yeast strain constructed through a modified genome shuffling method," Biotechnology for Biofuels, vol. 5, article no. 46, 2012.
 L. Hou, "Improved production of ethanol by novel genome shuffling in saccharomyces cerevisiae," Applied Biochemistry and Biotechnology, vol. 160, no. 4, pp. 1084-1093, 2010.
 W. Kabsch and C. Sander, "DSSP: definition of secondary structure of proteins given a set of 3D coordinates," Biopolymers, vol. 22, pp. 2577-2637, 1983.
Gammadde Hewa Ishan Maduka Wickramasinghe, (1) Pilimathalawe Panditharathna Attanayake Mudiyanselage Samith Indika Rathnayake, (1) Naduviladath Vishvanath Chandrasekharan, (1) Mahindagoda Siril Samantha Weerasinghe, (1) Ravindra Lakshman Chundananda Wijesundera, (2) and Wijepurage Sandhya Sulochana Wijesundera (3)
(1) Department of Chemistry, Faculty of Science, University of Colombo, Colombo, Sri Lanka
(2) Department of Plant Sciences, Faculty of Science, University of Colombo, Colombo, Sri Lanka
(3) Department of Biochemistry and Molecular Biology, Faculty of Medicine, University of Colombo, KynseyRoad, Colombo 08, SriLanka
Correspondence should be addressed to Wijepurage Sandhya Sulochana Wijesundera; firstname.lastname@example.org
Received 19 January 2017; Accepted 19 June 2017; Published 10 August 2017
Academic Editor: Sherry L. Mowbray
Caption: Figure 1: SDS-PAGE analysis of P. stipitis extracts. Lane (1): enzyme extract of nonrecombinant P. stipitis. Lane (2): recombinant EXNI enzyme (~24kDA) secreted by recombinant Y-pGAPZ[alpha]/gEXNI P. stipitis clone. Lane (3): broad range protein molecular weight marker.
Caption: Figure 2: Ramachandran map of modeled EXNI protein.
Caption: Figure 3: VERIFY3D score profile that shows more than 96.32% of residues having average 3D-1D score of greater than or equal to 0.2.
Caption: Figure 4: Protein-ligand complex from molecular docking.
Caption: Figure 5: Three-dimensional view of H bonds between ligand and the protein residues.
Caption: Figure 6: H bonds between ligand and the protein residues from LigPlot program.
Caption: Figure 7: Root mean square deviations (RMSD) of the backbone and radius of gyration (Rg) of the protein from 15 ns long MD trajectory.
Caption: Figure 8: Root mean square fluctuation (RMSF) of the residues in the protein over 15 ns long MD trajectory.
Table 1: Statistics of the 3D model of EXNI from the Ramachandran plot. Ramachandran plot statistics EXNI Amino acid in most favoured regions 143 92.3% Amino acid in additional allowed regions 11 7.1% Amino acid in generously allowed regions 1 0.6% Amino acid in disallowed regions 0 0.0% Number of nonglycine and nonproline residues 155 Number of end residues 2 Number of glycine residues 27 Number of proline residues 6 Total number of residues 190 Table 2: Detailed information of H bonds formed between ligand and the protein (A-Acceptor and D-Donor). Residue Amino acid Distance H-A Distance D-A 16 SER 2.79 3.39 173 ILE 2.84 3.47 77 TYR 3.14 2.94 Residue Donor angle Protein donor Side chain 16 119.03 No Yes 173 157.74 Yes No 77 69.35 Yes Yes
|Printer friendly Cite/link Email Feedback|
|Title Annotation:||Research Article|
|Author:||Wickramasinghe, Gammadde Hewa Ishan Maduka; Rathnayake, Pilimathalawe Panditharathna Attanayake Mudi|
|Publication:||BioMed Research International|
|Date:||Jan 1, 2017|
|Previous Article:||Antibiotic Elution from Hip and Knee Acrylic Bone Cement Spacers: A Systematic Review.|
|Next Article:||Protective and Predisposing Morphological Factors in Suprascapular Nerve Entrapment Syndrome: A Fundamental Review Based on Recent Observations.|