Expression, Docking, and Molecular Dynamics of Endo-β-1,4-xylanase I Gene of Trichoderma virens in Pichia stipitis

Wickramasinghe, Gammadde Hewa Ishan Maduka; Rathnayake, Pilimathalawe Panditharathna Attanayake Mudiyanselage Samith Indika; Chandrasekharan, Naduviladath Vishvanath; Weerasinghe, Mahindagoda Siril Samantha; Wijesundera, Ravindra Lakshman Chundananda; Wijesundera, Wijepurage Sandhya Sulochana

doi:https://doi.org/10.1155/2017/4658584

BioMed Research International

On this page

Abstract Introduction Materials and Methods Results and Discussion Conclusion Conflicts of Interest Acknowledgments References Copyright Related Articles

Research Article | Open Access

Volume 2017 | Article ID 4658584 | https://doi.org/10.1155/2017/4658584

Expression, Docking, and Molecular Dynamics of Endo-β-1,4-xylanase I Gene of Trichoderma virens in Pichia stipitis

Gammadde Hewa Ishan Maduka Wickramasinghe,¹Pilimathalawe Panditharathna Attanayake Mudiyanselage Samith Indika Rathnayake,¹Naduviladath Vishvanath Chandrasekharan,¹Mahindagoda Siril Samantha Weerasinghe,¹Ravindra Lakshman Chundananda Wijesundera,²and Wijepurage Sandhya Sulochana Wijesundera³

Academic Editor: Sherry L. Mowbray

Received19 Jan 2017

Accepted19 Jun 2017

Published10 Aug 2017

Abstract

It is essential that major carbohydrate polymers in the lignocellulosic biomass are converted into fermentable sugars for the economical production of energy. Xylan, the major component of hemicelluloses, is the second most naturally abundant carbohydrate polymer comprising 20–40% of the total biomass. Endoxylanase (EXN) hydrolyzes xylan into mixtures of xylooligosaccharides. The objective of this study was to genetically modify Pichia stipitis, a pentose sugar fermenting yeast species, to hydrolyze xylan into xylooligosaccharides via cloning and heterologous extracellular expression of EXNI gene from locally isolated Trichoderma virens species. Pichia stipitis was engineered to carry the EXNI gene of T. virens using pGAPZα expression vector. The open reading frame encodes 191 amino acids and SDS-PAGE analysis revealed a 24 kDA recombinant protein. The EXNI activity expressed by recombinant P. stipitis clone under standard conditions using 1% beechwood xylan was 31.7 U/ml. Molecular docking and molecular dynamics simulations were performed to investigate EXNI-xylan interactions. Free EXNI and xylan bound EXNI exhibited similar stabilities and structural behavior in aqueous medium. Furthermore, this in silico work opens avenues for the development of newer generation EXN proteins that can perform better and have enhanced catalytic activity.

1. Introduction

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 β-D-xylose units [7]. Xylanases, namely, β-1,4-xylosidases (XYL) and endo-xylanases (EXN), are responsible for depolymerization of xylan to xylose [8, 9]. Endo-β-1,4-xylanase (EXN) depolymerizes xylan to shorter xylobiose and xylooligomers [10] while β-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 industrial xylanase [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α integrative vector consisting of yeast α-mating factor (MFα) 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 [19]. The isolated species was cultivated for seven days on potato dextrose agar (PDA) at 30°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 (NH₄)₂SO₄ 1.4 g/L, KH₂PO₄ 2.0 g/L, CaCl₂ 0.3 g/L, MgSO₄·7H₂O 0.3 g/L, peptone 1.0 g/L, FeSO₄·7H₂O 5.0 mg/L, ZnSO₄·7H₂O 1.4 mg/L, MnSO₄·H₂O 1.6 mg/L, and CoCl₂ 2.0 mg/L, pH 5.5) supplemented with 2% beechwood xylan [20, 21]. Cultures were incubated at 30°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 6200 at 4°C for 10 minutes to obtain a cell free enzyme extract and stored at −20°C. Thereafter, the enzyme extracts were freeze dried and 100 μL of the crude enzyme extracts was tested for endoxylanase (EXN) activity [22] using 900 μL of 1% soluble beechwood xylan as the substrate. The reaction mixture was incubated at 50°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 (1 U) of endoxylanase activity is defined as the amount of 1 μmole of liberated hydrolysis product (equivalent to xylose) in 1 mL of enzyme volume per minute (μmol m1⁻¹ min⁻¹). 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 α-mating factor signal sequence of the pGAPZα vector.

Genomic DNA was extracted from T. virens using a simple extraction method [25] and EXN1 was amplified using the above designed primers. EXN1FP: 5′-ATCGTGAATTCCAGACGATTGGCCCCGGCACT-3′ and EXN1 RP: 5′-TTGATTCTAGATTAGCTYACGTTAATGTTGGCGTTACCAGAGCT-3′. PCR cycling conditions were optimized (initial denaturation at 94°C for 2 min followed by 35 cycles of denaturation at 94°C for 30 sec, annealing at 65°C for 30 sec, and extension at 72°C for 30 sec. This was followed by a final extension at 72°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 [26]. 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 μg/mL), 5-bromo-4-chloro-indolyl-D-galactoside (0.2 mM), and isopropyl-thio-β-D-galactopyranoside (40 μ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α vector. Ligated products were transformed into E. coli JM109 and transformants were selected in zeocin (25 μg/mL) containing low salt LB agar plates. The resulting clone was designated as pGAPZα/gEXNI.

2.3. Transformation into Pichia stipitis

The pGAPZα/gEXNI vector construct was linearized using BglII and purified. The plasmid concentration was determined and adjusted to ~1 μg/μ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°C overnight in a rotary shaker at 150 rpm. A volume of 500 μL from the above grown culture was inoculated into 50 mL broth in a 250 mL conical flask and incubated at 30°C at 150 rpm until the OD₆₀₀ reached 1.4. Yeast electrocompetent cells were prepared according to the procedure given in pGAPZα vector manual (Invitrogen, USA). A volume of 80 μL P. stipitis competent cells was mixed with 5–10 μg of linearized pGAPZα/gEXNI plasmid DNA. The mixture was subjected to electroporation under the optimized conditions of 1.5 Kw, 200 mA, and 25 μ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 1 M sorbitol) plates with 100 μg/mL zeocin as the selection marker. The plates were then incubated for 3 days at 30°C to obtain positive transformants. Twenty yeast colonies were selected and streaked on fresh YPDS plates containing zeocin (100 μ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α/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α/gEXNI was sequenced. Thereafter it was cultured in YPD broth with 100 μ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°C.

2.5. SDS-PAGE and Expression Analysis of EXNI in Recombinant P. stipitis

Y-pGAPZα/gEXNI recombinant P. stipitis was inoculated into beechwood xylan broth culture and incubated overnight at 30°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°C. The enzyme extract was concentrated by freeze drying. SDS-PAGE analysis [27] was performed and the results were analyzed by comparing with native P. stipitis.

Enzyme activity assay was carried out on the Y-pGAPZα/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°C. The enzyme harvest was freeze dried and a volume of 900 μL of solubilized 1% beechwood xylan was treated with 100 μ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 [32] 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 [33], PROCHECK [34], and ERRAT [35] 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 [38]. Protein and above optimized ligand were prepared using UCSF chimera [39], 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 [42].

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 [43]. The GROMOS54a7 united atom force field was assigned for the model protein. The force field parameters of the ligand were obtained from PRODRG server [44]. The protein-ligand complex was inserted in the center of a box with 9 × 9 × 9 nm³ volume. Electroneutrality of the system was maintained by adding Na⁺ ions to the box and solvated with SPC/E water molecules [45]. Electrostatic interactions were modeled by particle mesh Ewald (PME) with a short-range cutoff of 1.2 nm [46]. Berendsen’s weak coupling algorithm was employed to maintain the temperature and pressure of the system at 300 K and 1 bar [47]. Using LINCS algorithm [48] 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® Core™ 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 [49]. 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 387 bp). 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 EXNI 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 ~24 kDA (Figure 1).

The EXNI enzyme activity expressed by Y-pGAPZα/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α/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 (~24 kDA) 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-Val-Ser) 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 120 U/ml to 746 U/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 and XYL [64], 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α/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 [71] the secondary structure of EXNI is composed of 1 α-helix and 14 β-sheets as given below.

Alpha Helix. VAL152-LEU162

β-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 [33] 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, TYR88, 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⁻¹. 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 Å), ILE173 (3.47 Å), and TYR77 (2.94 Å). Detailed information is presented in Table 2. The stability of all these H bonds was studied using g_dist 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.

4. Conclusion

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.

Acknowledgments

The authors are grateful to National Science Foundation (NSF), Sri Lanka, for funding this research (Grant no. RG/2012/BT/02).

References

IEA, “Technology Roadmap Biofuels for Transport,” International Energy Agency: France, pp. 19–145, 2011.
View at: Google Scholar
J. Goldemberg, “Biomass and energy,” Química Nova, vol. 32, no. 3, pp. 582–587, 2009.
View at: Publisher Site | Google Scholar
Y. Sun and J. Cheng, “Hydrolysis of lignocellulosic materials for ethanol production: a review,” Bioresource Technology, vol. 83, no. 1, pp. 1–11, 2002.
View at: Publisher Site | Google Scholar
Ó. J. Sánchez and C. A. Cardona, “Trends in biotechnological production of fuel ethanol from different feedstocks,” Bioresource Technology, vol. 99, no. 13, pp. 5270–5295, 2008.
View at: Publisher Site | Google Scholar
M. P. Coughlan and G. P. Hazlewood, “β-1,4-D-Xylan-degrading enzyme systems: biochemistry, molecular biology and applications,” Biotechnology and Applied Biochemistry, vol. 17, pp. 259–289, 1993.
View at: Google Scholar
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.
View at: Publisher Site | Google Scholar
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.
View at: Google Scholar
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.
View at: Publisher Site | Google Scholar
M. Sharma and A. Kumar, “Xylanases: an overview,” British Biotechnology Journal, vol. 3, no. 1, pp. 1–28, 2013.
View at: Publisher Site | Google Scholar
T. Collins, C. Gerday, and G. Feller, “Xylanases, xylanase families and extremophilic xylanases,” FEMS Microbiology Reviews, vol. 29, no. 1, pp. 3–23, 2005.
View at: Publisher Site | Google Scholar
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.
View at: Publisher Site | Google Scholar
H. Dalbøge, “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.
View at: Publisher Site | Google Scholar
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.
View at: Google Scholar
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. 3-4, pp. 228–233, 1986.
View at: Publisher Site | Google Scholar
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.
View at: Publisher Site | Google Scholar
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.
View at: Publisher Site | Google Scholar
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.
View at: Publisher Site | Google Scholar
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.
View at: Publisher Site | Google Scholar
N. Sadfi-Zouaoui, S. Châabani, 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.
View at: Publisher Site | Google Scholar
M. Mandels and D. Sternburg, “Recent advances in cellular technology,” Journal of Fermentation Technology, vol. 54, pp. 267–286, 1976.
View at: Google Scholar
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.
View at: Google Scholar
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.
View at: Publisher Site | Google Scholar
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.
View at: Publisher Site | Google Scholar
T. K. Ghose, “Measurement of cellulase activities,” Pure and Applied Chemistry, vol. 59, no. 2, pp. 257–268, 1987.
View at: Publisher Site | Google Scholar
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.
View at: Publisher Site | Google Scholar
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.
View at: Publisher Site | Google Scholar
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.
View at: Publisher Site | Google Scholar
H. M. Berman, J. Westbrook, Z. Feng et al., “The protein data bank,” Nucleic Acids Research, vol. 28, no. 1, pp. 235–242, 2000.
View at: Publisher Site | Google Scholar
H. Berman, K. Henrick, and H. Nakamura, “Announcing the worldwide Protein Data Bank,” Nature Structural Biology, vol. 10, no. 12, p. 980, 2003.
View at: Publisher Site | Google Scholar
T. L. Madden, R. L. Tatusov, and J. Zhang, “Applications of network BLAST server,” Methods in Enzymology, vol. 266, pp. 131–141, 1996.
View at: Publisher Site | Google Scholar
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 Site | Google Scholar
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: Publisher Site | Google Scholar
C. Colovos and T. O. Yeates, “ERRAT: an empirical atom-based method for validating protein structures,” Protein Science, vol. 2, pp. 1511–1519, 1993.
View at: Google Scholar
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.
View at: Publisher Site | Google Scholar
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 Site | Google Scholar
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.
View at: Publisher Site | Google Scholar
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.
View at: Publisher Site | Google Scholar
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.
View at: Publisher Site | Google Scholar
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.
View at: Publisher Site | Google Scholar
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.
View at: Publisher Site | Google Scholar
A. W. Schüttelkopf 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.
View at: Publisher Site | Google Scholar
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.
View at: Publisher Site | Google Scholar
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.
View at: Publisher Site | Google Scholar
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.
View at: Publisher Site | Google Scholar
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.
View at: Publisher Site | Google Scholar
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.
View at: Publisher Site | Google Scholar
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.
View at: Publisher Site | Google Scholar
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.
View at: Publisher Site | Google Scholar
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.
View at: Publisher Site | Google Scholar
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.
View at: Google Scholar
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.
View at: Google Scholar
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.
View at: Publisher Site | Google Scholar
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.
View at: Publisher Site | Google Scholar
V. Ruanglek, R. Sriprang, N. Ratanaphan et al., “Cloning, expression, characterization, and high cell-density production of recombinant endo-1,4-β-xylanase from Aspergillus niger in Pichia pastoris,” Enzyme and Microbial Technology, vol. 41, no. 1-2, pp. 19–25, 2007.
View at: Publisher Site | Google Scholar
Y. Cao, J. Qiao, Y. Li, and W. Lu, “De novo synthesis, constitutive expression of Aspergillus sulphureus β-xylanase gene in Pichia pastoris and partial enzymic characterization,” Applied Microbiology and Biotechnology, vol. 76, no. 3, pp. 579–585, 2007.
View at: Publisher Site | Google Scholar
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.
View at: Publisher Site | Google Scholar
P. Fickers, “Pichia pastoris: a workhorse for recombinant protein production,” Current Research in Microbiology and Biotechnology, vol. 2, no. 3, pp. 354–363, 2014.
View at: Google Scholar
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.
View at: Publisher Site | Google Scholar
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.
View at: Publisher Site | Google Scholar
S. Fernandes and P. G. Murray, “Metabolic engineering for improved microbial pentose fermentation,” Bioengineered Bugs, vol. 01, no. 06, pp. 424–428, 2010.
View at: Publisher Site | Google Scholar
G. H. I. M. Wickramasinghe, N. V. Chandrasekharan, S. Rathnayake, S. Weerasinghe, R. L. C. Wijesundera, and W. S. S. Wijesundera, “β-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.
View at: Google Scholar
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.
View at: Publisher Site | Google Scholar
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.
View at: Publisher Site | Google Scholar
G. Davies and B. Henrissat, “Structures and mechanisms of glycosyl hydrolases,” Structure, vol. 3, no. 9, pp. 853–859, 1995.
View at: Publisher Site | Google Scholar
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.
View at: Publisher Site | Google Scholar
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.
View at: Publisher Site | Google Scholar
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.
View at: Publisher Site | Google Scholar
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.
View at: Google Scholar

Copyright

Copyright © 2017 Gammadde Hewa Ishan Maduka Wickramasinghe et al. 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.

PDF Download Citation

Download other formats

Order printed copies

Views

2391

Downloads

1208

Citations