Unique N-glycosylation of a recombinant exo-inulinase from Kluyveromyces cicerisporus and its effect on enzymatic activity and thermostability
Journal of Biological Engineering volume 13, Article number: 81 (2019)
Inulinase can hydrolyze polyfructan into high-fructose syrups and fructoligosaccharides, which are widely used in food, the medical industry and the biorefinery of Jerusalem artichoke. In the present study, a recombinant exo-inulinase (rKcINU1), derived from Kluyveromyces cicerisporus CBS4857, was proven as an N-linked glycoprotein, and the removal of N-linked glycan chains led to reduced activity.
Five N-glycosylation sites with variable high mannose-type oligosaccharides (Man3–9GlcNAc2) were confirmed in the rKcINU1. The structural modeling showed that all five glycosylation sites (Asn-362, Asn-370, Asn-399, Asn-467 and Asn-526) were located at the C-terminus β-sandwich domain, which has been proven to be more conducive to the occurrence of glycosylation modification than the N-terminus domain. Single-site N-glycosylation mutants with Asn substituted by Gln were obtained, and the Mut with all five N-glycosylation sites removed was constructed, which resulted in the loss of all enzyme activity. Interestingly, the N362Q led to an 18% increase in the specific activity against inulin, while a significant decrease in thermostability (2.91 °C decrease in Tm) occurred, and other single mutations resulted in the decrease in the specific activity to various extents, among which N467Q demonstrated the lowest enzyme activity.
The increased enzyme activity in N362Q, combined with thermostability testing, 3D modeling, kinetics data and secondary structure analysis, implied that the N-linked glycan chains at the Asn-362 position functioned negatively, mainly as a type of steric hindrance toward its adjacent N-glycans to bring rigidity. Meanwhile, the N-glycosylation at the other four sites positively regulated enzyme activity caused by altered substrate affinity by means of fine-tuning the β-sandwich domain configuration. This may have facilitated the capture and transfer of substrates to the enzyme active cavity, in a manner quite similar to that of carbohydrate binding modules (CBMs), i.e. the chains endowed the β-sandwich domain with the functions of CBM. This study discovered a unique C-terminal sequence which is more favorable to glycosylation, thereby casting a novel view for glycoengineering of enzymes from fungi via redesigning the amino acid sequence at the C-terminal domain, so as to optimize the enzymatic properties.
Inulinases can hydrolyze β-2,1-D-fructosidic linkages in polyfructan, which are found in a wide range of plants and microorganisms, among which yeasts, especially Kluyveromyces, generally produces higher levels of inulinase than other filamentous fungi or bacterial . The hydrolysis products, high-fructose syrups and fructoligosaccharides, are widely used in food, medical and bioenergy industries. Inulinases are divided into endo-inulinase (EC 184.108.40.206) and exo-inulinase (EC 220.127.116.11) according to the mode of action on inulin . Exo-inulinase as the key enzyme in Jerusalem artichoke biorefinery can hydrolyze inulin from the non-reducing end, obtaining the high-fructose syrups with purity up to 95% within one step, which can be further transformed into fuel ethanol, single cell oil, 2,3-butanediol, lactic acid or other chemical production . Therefore, exo-inulinase has attracted much attention in scientific research for its potential application in biomass processing.
Exo-inulinases belong to the family of glycoside hydrolase 32 (GH32), a family including other inulinases, invertases, β-fructofuranosidases and fructosyltransferases, which share a common structural feature: a N-terminal five-blade β-propeller catalytic domain with four antiparallel β-strands exhibiting ‘W’ topology for each blade, surrounding the catalytic active center, followed by a C-terminal β-sandwich domain constituting with two antiparallel six-stranded β-sheets . At present, though with much more penetration in the research about the N-terminal catalytic module, little is known about the function of the C-terminal domain. However, researchers found that it might be involved in substrate binding [5, 6] and regulating substrate specificity . For example, the β-sandwich domain of a dimeric β-fructofuranosidase from Schwanniomyces occidentalis was found to be involved in substrate binding through the interaction between β-sandwich domain of the adjacent subunit within the dimer with the substrate ; and the C-terminal domain (BsCBM66) of the exo-acting levanase SacC from Bacillus subtilis proved to belong to CBM66 with the function of identifying and binding substrate, as well as facilitating the orienting of the catalytic domain to the substrate [8, 9], thereby enhancing enzymatic activity through increasing the concentration of the appended enzymes in the vicinity of the substrate. It was just the extensive interactions of BsCBM66 with the terminal fructose moiety (Fru-3) of levantriose that conferred SacC the substrate specificity .
In recent years, many exo-inulinases were cloned in various hosts and characterized [2, 10,11,12], and it was found that exo-inulinases expressed by yeast strains, such as P. pastoris, were usually modified by N-linked glycosylation [1, 13], one of the most common and important post-translational modifications in eukaryotic cells. N-glycosylation often affects the folding and conformation of proteins [14, 15], secretory expression [16,17,18], enzyme activity  and stability . The effects of N-glycosylation modification on the enzyme activity are quite complicated, which are including glycosylation sites, glycoforms and the length of glycan chains, etc. Glycosylation at different sites sometimes could lead to opposite effects on enzyme activity. As for Cellobiohydrolases (CBHs) PvCel6A and PvCel7A, both from Penicillium verruculosum belonging to GH6  and GH7 family  respectively, the removal of N-linked glycans close to the entrance to the enzyme active center turned out to be an increase in the rate of cellulose hydrolysis likely due to the elimination of steric hindrance and the ease of threading the cellulose molecule into the tunnel. Meanwhile the deletion of N-linked glycans at the positions locating at the “bottom” of their catalytic domains resulted in a partial decrease of enzyme activity, because the glycan chains attaching to these residues could help the enzyme orient appropriately to the surface of cellulose microfibril. In most cases, the loss of N-linked glycosylation would reduce the enzyme activity, which was probably caused by the change of secondary structure via circular dichroism spectroscopy analysis . Meanwhile it has been reported that the enzyme activity is increased by introducing glycosylation sites recently [18, 24]. In addition, enzymes which are from different expression systems often bring different activities due to variable glycoforms or the length of glycan chains. Wei et al. have reported that the hydrolytic activity of Aspergillus terreus β-glucosidase expressed by P.pastoris was lower than that by Trichoderma reesei, with a conclusion that longer N-linked glycan chains could weaken the substrate affinity .
In general, N-glycosylation plays diverse roles in regulating enzyme activity, though its mechanism is still not completely understood. The impact of N-glycosylation on the conformation of enzyme [14, 26,27,28] and the substrate affinity [25, 29] are probably the main factors for the changed activity. Gao et al. raised the hypothesis of “synergism” based on the research on CBHI from Penicillium decumbens, i.e. N-glycosylation of CBHI-A, which has no detectable hydrolysis activity, could favor for other cellulases to hydrolyze substrates due to the disruption of the hydrogen bonds from cellulose and the convertion of the high crystalline cellulose into amorphous types . In addition, Dotsenko and Gusakov et al. pointed out the novel concept of “lifting factor” for the function of N-linked glycosylation . They believed N-glycosylation located at the “bottom” of the catalytic domains of CBHs belonging to GH6 and GH7 families could help the catalytic domain of the enzymes positioning along the microcrystalline chains through non-specific dynamic interaction, followed by the detachment of a single cellulose chain from the crystal, i.e. the N-linked glycans served as a chain “lifting actor”, and then threaded substrates into the active tunnel.
However, little is known about the effects of N-glycosylation on exo-inulinase activity compared to the cellulose, and still with no report about the effect of N-glycosylation on the activity of inulinase at present. For GH32 family, limited research has reported about the influence of N-glycosylation on substrate specificity of fructan 1-exohydrolase , invertase activity or stability [32, 33] and β-fructofuranosidase activity as well as oligomerization .
In the present study, the recombinant N-glycoprotein rKcINU1 was expressed in P. pastoris X-33. Meanwhile single-site glycosylation mutants and the Mut with five N-glycosylation sites removed were constructed to investigate the effect of N-glycans on enzymatic activity, kinetics properties, the second structure and thermostability. Specifically, its five N-glycosylation sites were found to be located at the C-terminal β-sandwich domain and likelihood mechanism was given based on the experimental results and 3D modeling, which was the N-linked glycan chains endowed the β-sandwich region with the functions of CBM.
Materials and methods
Strains, plasmids and culture conditions
Escherichia coli TOP10 (TaKaRa, Dalian, China) was used as the host strain for plasmid amplification. P. pastoris X-33 (Invitrogen, Carlsbad, USA) served as the eukaryotic expression host for different enzymes. The plasmid pPICZαA (Invitrogen, Carlsbad, USA) was used to construct expression vectors. E. coli strain was cultivated in low salt Luria-Bertani (LB) medium containing 25 μg/mL zeocin at 37 °C. P. pastoris was grown in YPDSZ medium (w/v): 1% yeast extract, 2% peptone, 2% glucose, 1 M sorbitol and 100 μg/mL zeocin at 30 °C for selecting transformants. Protein expression was performed in BMGY and BMMY media (w/v) at 28 °C: 1% yeast extract, 2% peptone, 1.34% YNB, 4 × 10− 5% biotin, 100 mM potassium phosphate, pH 6.0 and 1% glycerol or 0.5% methanol as carbon sources, respectively.
Site-directed mutagenesis, protein expression and purification
The exo-inulinase gene kcINU1 from K. cicerisporus CBS4857 (GenBank™ accession number AF178979) was amplified and cloned into the pPICZαA expression vector to yield the plasmid pPICZαA-rkcINU1 as described previously . Restriction-free cloning was carried out to construct the N-glycosylation mutations N362Q, N370Q, N399Q, N467Q and N526Q. The Asn codons were mutated to Gln using specific primers listed in Additional file 1: Table S1 with pPICZαA-rkcINU1 as a template following the method published previously . A variant designated as Mut with all five glycosylated Asn sites mutated into Gln was also obtained. The protein expression and purification followed the method described by Ma , with details described in Additional file 1.
Glycosidase treatment, Zymogram analysis and enzymes activity assays
The deglycosylation reaction of purified wild-type rKcINU1 using PNGase F (New England Biolabs, Beverly, MA, USA) and EndoF1 (Sigma, Saint Louis, MO, USA) were carried out according to the specification. The zymogram analysis and activity assay by DNS method were described as referred previously . The details above were given in Additional file 1.
For kinetics study, the initial velocities were measured by incubating 0.4 μg of purified rKcINU1 and its variants with 1–20 mg/mL sucrose at 55 °C for 20 min in 100 mM acetate buffer, pH 4.5. The initial rates versus substrate concentrations were hyperbola fitted based on the Michaelis-Menten equation using Origin program (version 8). Kinetics parameters Km and Vmax were calculated from the curve.
Mass spectrometry analysis
Glycosylation sites of the purified wild-type rKcINU1 and Mut were confirmed by using the mass spectrometry analysis described previously by Jiang et al. . The details involved were given in Additional file 1.
The three-dimensional structure of the wild-type rKcINU1 was built by the Swiss-Model server (http://swissmodel.expasy.org/) based on the crystal structure of S.cerevisiae invertase (PDB code 4EQV) which sharing 54.82% sequence identity with the rKcINU1. Docking of the N-linked glycan chains to the model mentioned above was then followed using the program GLYCAM, whereas the (Man)7(GlcNAc)2 glycan chains were selected from the high-mannose oligosaccharide library.
Circular Dichroism (CD) spectroscopy analysis
Far-UV CD spectra of the purified wild-type rKcINU1 and its mutants were measured in 20 mM phosphoric acid (pH 4.5) at a protein concentration of 0.2 mg/mL with a MOS-450 (Bio-Logic, Seyssinet-Pariset, France) spectropolarimeter using 2 mm path length quartz cuvettes. The protein spectra were recorded over a wavelength range from 190 nm to 260 nm at 25 °C with a 1 nm bandwidth. The speed of scanning was 50 nm/min. For each sample, the average of three spectra was collected by subtraction of the equivalent buffer spectrum. The secondary structure of each sample was analyzed with CDpro software.
Differential scanning Calorimetry (DSC) analysis
The thermodynamic parameters of wild-type rKcINU1 and its variants were determined on a Nano-Differential Scanning Calorimeter (Nano-DSC, TA Instruments, New Castle, DE, USA). The purified recombinants mentioned above were prepared in 100 mM sodium acetate buffer at pH 4.5 with a final protein concentration of 1 mg/mL. The samples were heated in the temperature ranging from 25 °C to 110 °C with the acceleration rate of 1 °C/min. All samples were degassed in vacuum prior to subjecting to the system, with a baseline of 100 mM sodium acetate buffer, pH 4.5. The protein spectra were recorded by subtraction of the baseline. Data analysis was performed on the DSC run software.
Results and discussion
The rKcINU1 modified by N-glycosylation
The recombinant rKcINU1, which corresponds to the mature enzyme sequence with an additional N-terminal His-tag, consists of 538 amino acid residues, with its theoretical molecular mass deduced from its primary sequence as 60.54 kDa. However, the band of purified rKcINU1 smeared between 70 kDa and 100 kDa (Fig. 1 a, Lane 1), which is much higher than the calculated mass of 60.54 kDa. This is likely due to the glycosylation modification of the enzyme by P. pastoris, which has the ability to post-translationally modify recombinant proteins [1, 37]. An analysis of the amino acid sequence of rKcINU1 produced 1 potential N-glycosylation sites using the website named NetNGlyc 1.0 Server (http://www.cbs.dtu.dk/services/NetNGlyc/), which were then used to predict the N-glycosylation sites (Fig. 1 b). Next, rKcINU1 was digested with PNGase F, which was known to cleave N-linked glycans between the innermost acetylglucosamine (GlcNAc) and the Asn residue, a single band instead of smear appeared in the SDS-PAGE (Fig. 1 a, Lane 2), which was almost exactly consistent with the theoretical molecular weight of 60.54 kDa. Later, MALDI-TOF mass spectrometry was used to obtain the precise molecular mass of rKcINU1 79.5 kDa (Fig. 1 c). All of the data indicated that rKcINU1 was an N-glycoprotein, and its average degree of glycosylation was 23.9%. Furthermore, rKcINU1 was digested with PNGase F under natural conditions (Fig. 1 a, Lane 3), and the results showed a band of about 65 kDa on the gel, hinting at the fact that rKcINU1 possessed more than one glycosylation site, some of which were likely buried deep inside the enzyme, where they inhibited the accessibility of PNGase F to those sites.
Effect of Deglycosylation by Endo F1 on enzyme activity
In order to investigate the effects of N-linked oligosaccharide chains on the hydrolytic activity of rKcINU1, purified rKcINU1 was subjected to Endo F1 treatment under non-denaturing conditions, followed by zymogram analysis and DNS assay to detect enzymatic activity. Endo F1 is not sensitive to protein conformation, and is thus suitable for the cleaving of oligosaccharide chains, even those deeply buried under natural conditions, though it leaves a single GlcNAc unit at the Asn residue. The hydrolytic activity of purified rKcINU1 with glycosylation modification was shown to be higher with both inulin and sucrose as the reaction substrate, compared to that of the deglycosylated rKcINU1 as presented in situ (Fig. 2). In other words, glycosylation modification allowed rKcINU1 to maintain high hydrolytic activity, whereas the removal of oligosaccharide moieties dramatically decreased the enzymatic activity, while still carrying the inner GlcNAc residue. These results confirmed not only the glycosylation modification, but also the fact that the length of oligosaccharide chains contributed to the positive effect on the enzymatic activity.
The DNS method was applied to quantitatively measure the decrement extent of enzyme activity, due to the removal of the N-linked glycan chains. After treatment by Endo F1 for 4 h (Additional file 1: Table S2), the hydrolytic activity of rKcINU1 could only keep around 70% against inulin and 63% toward sucrose; and when more glycan chains were removed, the data for both substrates fell to between 30 and 40% when treated for 16 h. These results further strengthened the conclusion that N-glycosylation modification was necessary to maintain the enzyme activity of rKcINU1.
Determination of the glycosylation sites of the wild-type rKcINU1 by mass spectrometry analysis
In order to identify N-linked glycans sites, the purified wild-type rKcINU1 was subjected to trypsin hydrolysis, followed by PNGase F digestion, prior to mass analysis on the LTQ-Orbitrap Elite hybrid mass spectrometer. Mass spectrometry assigned a single deamidation modification to each Asn residue within the Asn-X-Ser/Thr (X ≠ Pro) motif, producing the five sites of Asn-362, Asn-370, Asn-399, Asn-467 and Asn-526 (Additional file 1: Figure S1). While no mass alteration was found on the PNGase F treated trypsin peptides containing Asn-9, Asn-99, Asn-100, Asn-106, Asn-153, Asn-203 or Asn-267, the absence of N-glycosylation modification of those putative sites in the wild-type rKcINU1 was demonstrated. Quite interestingly, the authentic glycosylation modification site Asn-399, found in the VFNGSP sequence, was not included in the list of putative sites based on prediction in silicon. The ESI-MS profile from the wild-type rKcINU1 proved the variable high mannose-type oligosaccharides (Man3–9GlcNAc2) carrying up to nine mannose residues of the N-linked glycan chains (Fig. 3). Furthermore, the MALDI-TOF mass spectra of the wild-type rKcINU1 suggested different glycoform ensembles caused by glycoform heterogeneity , each with variable numbers of mannose residues, as indicated by arrows for the typical clusters in Fi. 1C. These results demonstrated that only 4 of the 11 predicted N-glycosylation sites, plus the Asn-399 site in the wild-type rKcINU1 sequence, were finally modified by the heterogeneous high mannose-type glycan chains.
Three-dimensional modeling of the wild-type rKcINU1
Next, the glycoprotein 3D model structure of the wild-type rKcINU1 with (Man)7(GlcNAc)2 glycans was built (Fig. 4). As illustrated, the wild-type rKcINU1 adopted the bimodular structure typical of the GH32 family. It folded into a five-blade β-propeller catalytic domain at the N-terminus. This accommodated for the active site of the enzyme (with each blade composed of four antiparallel β-strands with a ‘W’ topology) and a C-terminus β-sandwich domain assembled from two six-stranded antiparallel β-sheets, with a short loop linking the two domains . Three key acidic amino acid residues were located in the active site, and were responsible for substrate binding and hydrolysis. These included two Asp and one Glu residues in the conserved sequences of GH32 family, namely NDPNG (D acted as anucleophile) and RDP (D acted as a stabilizer of the transient state), as well as EC (E acted as an acid base catalyst).
The priority of N-glycosylation for the C-terminus domain with a unique sequence
A unique feature of the wild-type rKcINU1 model was that the five glycosylation sites mentioned above were located at the C-terminus β-sandwich domain, all of which were placed on the same six-stranded β-sheet located opposite the catalytic active center. As shown in the built three-dimensional model (Fig. 4), Asn-362 at the C-terminus domain near the linker was quite close to the entrance for the substrate into the catalytic active cavity, in the upper open pocket of the ‘W’ topology. Meanwhile, the other N-glycosylation sites were located on the β-strands, or turns in the C-terminus β-sandwich domain at the reverse side of the catalytic cavity, isolated from the catalytic active pocket. The multiple sequence alignment of KcINU1 (native sequence), along with the other three members of the family GH32 whose glycosylation sites have been identified in previous research from fungi, was used (Fig. 5). However, no conservative N-glycosylation sites were found in any of the 4 enzymes, even between KcINU1 and ScINV, which shared the highest sequence homogeneity of up to 52%, with only two N-glycosylation sites conserved: Asn-370 and Asn-399 in KcINU1 equivalent to Asn-349 and Asn-378 in ScINV, respectively. Meanwhile, Asn-362, Asn-467 and Asn-526 in KcINU1 represent “unique” N-glycosylation sites, which are absent in all other enzymes analyzed. The most remarkable feature was that no authentic N-glycosylation modification in KcINU1 was located in the N-terminal catalytic domain, despite its seven putative glycosylation sites in this sequence observed from the prediction results (Fig. 1 b). This is unlike three other enzymes possessing N-linked glycan sites in the catalytic domain and the β-sandwich domain simultaneously. All of these results showed that the C-terminal sequence of the KcINU1 may be unique among the GH32 family in fungi, thus ensuring the priority of glycosylation modification in the C-terminus β-sandwich domain, instead of the N-terminal catalytic domain. Based on previous research, an Asn residue in an Asn-X-Thr sequon was more readily glycosylated than that in an Asn-X-Ser sequon . In addition, the middle residue with small hydrophobic and positively charged side chain of an Asn-X-Ser consensus motif was more efficiently N-glycosylated . Recently, Huang et al. discovered the enhanced aromatic sequon (Phe-X-Asn-X-Thr and Phe-X-X-Asn-X-Thr) could be efficiently N-glycosylated, namely aromatic residues, especially Trp and sulfur-containing residues at the X1 position in five-residue sequon X1-X2-N-X3-S/T, improved N-glycosylation efficiency, while positively charged residues such as Arg suppressed N-glycosylation . However, none of the theories described above could rationally explain the high efficiency of glycosylation at Asn-399 within the VFNXSP sequence. As shown in Fig. 1 b, Asn-399 actually possessed little likelihood of being glycosylated according to the prediction result compared to Asn-9, Asn-99, Asn-203 and Asn-267 (all in NXT motifs), which in the end were not glycosylated for unknown reasons. Interestingly, Asn-399 in the VFNXSP sequence was located right between two consecutive NXT motifs (NXT-NXT-NXS-NXT-NXT), whose binding to oligosaccharyltransferase (OST) may hinder the potential skipping of OST through Asn-399, also for unknown reasons .
In any case, the unique sequence above was favorable to N-glycosylation modification, and thus provided a novel insight of artificially designing the glycosylation modification in the C-terminal β-sandwich domain to fine-tune the enzyme function. Furthermore, in previous research regarding the bimodular CBH enzymes, it was the N-terminus domain incorporating the catalytic cavity that was more ready to be N-glycosylated, while the C-terminal domain and the linker between the two were sensitive to O-glycosylation [21, 22, 41, 42]. However, in this case it was observed that all five N-glycosylation modification sites were located in the C-terminus non-catalytic domain, far from the active cavity, while still efficiently regulating the enzymatic activity. This implied a novel regulation mechanism for the enzyme properties, as well as a model of bimodular enzymes, by which to study the effects of N-glycosylation in non-catalytic domains on enzyme properties.
Mutagenesis and expression of rKcINU1
To elucidate the precise role which N-glycosylation modification played in maintaining the activity and stability of rKcINU1, each glycosylation site was mutated to Gln by means of site-directed mutagenesis. In addition, a mutant without any N-glycosylation sites (mut) was obtained by mutating all five Asn residues mentioned above to Gln. The DNA sequencing of the entire gene fragment confirmed that the amino acid substitution and no other undesired mutations. Then, the six mutated genes were successfully expressed in P. pastorisX-33 and purified by nickel affinity chromatography, yielding six recombinant proteins, i.e. N362Q, N370Q, N399Q, N467Q, N526Q and Mut (Fig. 6 a). Among these the variant Mut was expected to be expressed without N-glycosylation. However, to some extent N-glycosylation was still unexpectedly detected in the Mut from the smear band between 70 and 100 kDa in SDS-PAGE as the wild-type. In order to determine the exact glycosylation sites, the Mut was hydrolyzed by trypsin, followed by PNGase F treatment, then applied to the LTQ-Orbitrap Elite hybrid mass spectrometer. Astonishingly, mass spectrometry revealed six new N-linked glycan sites which were all located at the N-terminal catalytic domain (Fig. 6 b), i.e. Asn-9, Asn-147, Asn-153, Asn-197, Asn-203, Asn-232 (Additional file 1: Figure S2). This once again clearly illustrated the priority of N-glycosylation modification at the C-terminal sequence.
In order to elucidate the glycan pattern, purified proteins with single-site mutation digested with PNGase F were analyzed via mass spectrometry. The ESI-MS profiles revealed that the N-linked glycans of mutants N362Q, N370Q, N399Q, N467Q and N526Q respectively represented Man4–9GlcNAc2, Man4, 6–9GlcNAc2, Man4, 6–9GlcNAc2, Man4, 7–9GlcNAc2, and Man4, 7–9GlcNAc2 (Additional file 1: Figure S3).
N-linked Glycans are crucial for enzyme activity
The specific activities of the wild-type rKcINU1 and its mutants against inulin and sucrose are shown in Fig. 6 c. In general, mutants with a single N-glycosylation site mutated (N362Q, N370Q, N399Q, N467Q and N526Q) kept activity to some extent toward inulin and sucrose, whereas the Mut lost all hydrolytic activity, likely due to defective protein folding caused by incorrect glycosylation. This fact suggested the conclusion that the N-glycosylation modification at the C-terminal β-sandwich domain was necessary for the correct folding of rKcINU1, and essential for keeping an active form. Among the single mutants, N362Q led to an 18% increase in the specific activity against inulin in comparison with the wild-type, without any effect of hydrolysis activity toward sucrose on the other hand. This demonstrated that the N-glycosylationat Asn-362 negatively affected the rKcINU1 activity. On the contrary, other single mutations, such as N370Q, N399Q, N467Q and N526Q, resulted in a decrease of the specific activity to various degrees toward inulin and sucrose compared to the wild-type. Among these, the most dramatic decline of specific activity was observed in the case of N467Q, respectively to 40.71 and 54.66% of enzyme activity against inulin and sucrose compared to the wild-type, thus indicating that the glycosylation at Asn-467 played the most essential role in sustaining the enzymatic activity of the wild-type. In addition, N526Q presented the least significant reduction in specific activity, both for inulin and sucrose as a substrate. These results implied the potential regulation of N-glycosylation in maintaining enzymatic activity with contribution to various extents, according to the glycosylated site. Meanwhile, the glycosylation at Asn-467 exhibited the most outstanding effect on maintaining enzyme activity, while the glycosylation at Asn-362 led to a negative effect in the hydrolysis activity of the enzyme toward inulin.
Next, the kinetic properties with sucrose as the substrate for both the wild-type rKcINU1 and its mutants were investigated (Table 1). The data showed that the Km values of single-site mutants were higher than that of the wild-type. They also had lower apparent catalytic efficiency (kcat/Km), such as N467Q with the lowest activity among all single-site mutants, presenting a 46-fold lower kcat/Km value and 2-fold higher Km value compared to the wild-type. These results indicated that N-glycosylation modification could positively affect the substrate-enzyme affinity, which may be a key reason for the decline in the activities of the mutants.
Effect of N-glycosylation on the second structure of rKcINU1
To elucidate whether N-glycosylation modification had an effect on the protein conformation, the secondary structures of the wild-type rKcINU1 and its mutants were determined using CD spectra in the far-UV spectrum region (Fig. 7, Table 2). The far-UV CD spectrum of the wild-type exhibited a strong negative peak around 218 nm, a typical feature for the protein rich in β-sheet , further followed by a quantitative analysis giving the ration of β-sheet of 42.11%, which is consistent with the feature of the GH32 family with abundant β-sheet. The CD spectra of N362Q, N370Q, N467Q and N526Q showed similar profiles to that of the wild-type, with an increase in the negative peak approximately at 218 nm. This corresponded to an increase in content of both the β-sheet and random coil conformation, but a decline in the ratio of α-helix. For N399Q, an additional positive peak around 192 nm was observed, except for the regular negative peak around 218 nm, which was reflected with the highest increase (by 7%) in the ratio of β-helix from 42.11 to 44.98%, as well as a decrease in random coil (by 5%) and β-turn (by 6%). These results suggested that glycosylation mutations altered the secondary structure of the rKcINU1. Therefore, it was assumed that the change in the secondary structure may have led to altered enzyme activity through fine-tuning substrate affinity.
Effect of N-glycosylation on thermal stability of rKcINU1
Glycosylation is closely related to thermal stability in glycoproteins. Therefore, the impact of N-glycosylation on the thermostability of rKcINU1 was evaluated by DSC, and the results are presented in Table 3. In general, the melting temperature values (Tm) of single-site glycosylation mutants were lower than that of the wild-type (88.89 °C), whereas N362Q outweighed the other mutants with 2.91 °C compared with the wild-type. In other words, the protein with removal of the N-linked glycan chains at Asn-362 would lead to significant decrease in thermostability, thereby implying that the glycosylation at Asn-362 contributed to stabilizing the structure of rKcINU1.
In addition, only the Tm of N526Q was slightly increased, as shown in Table 3, partially due to its location too close to glycosylation site Asn-362. This might owe to the intensive distribution of glycan chains, which contributed to the interactions between glycans instead of the protein and glycan, which is not conducive to protein thermostability. In contrast, more dispersed glycosylation sites would favor protein-glycan interactions, thus enhancing thermostability .
Furthermore, it was concluded that the glycan position, rather than the glycan pattern, had important effects on enzyme thermostability , just as the ESI-MS analysis (Additional file 1: Figure S3) indicated that the glycan patterns of N370Q were similar to those of N399Q, and the glycan patterns of mutant N467Q were consistent with those of N526Q. However, each mutant exhibited varying thermostability.
Mechanism of the effect of N-glycosylation on enzyme properties
From all the data above, it was concluded that the N-glycosylation at C-terminus β-sandwich domain had the following functions. First of all, the N-linked glycans could stabilize the C-terminal domain by keeping the β-strands in appropriate conformation, which was responsive to both enzymatic activity and thermostability. For example, Asn-467 was located at the loop region between two antiparallel β-sheets in the C-terminal domain. It was concluded that the glycan chains at this site were involved in stabilizing the conformation of the β-sandwich domain, similarly to those at Asn-297 of the immunoglobulin G1 Fc, which played a role in stabilizing the C’E loop through intramolecular interactions between carbohydrate and amino acid residues. Without these the Cγ2 domain would collapse, further rendering Fc incapable of binding FcγRІІІA [44, 45]. Another case was the N-glycosylation of the ectodomain (ECD) in the epidermal growth factor receptor (EGFR), which was crucial for the conformational arrangement of the ECD thus effecting its orientation to the membrane and the ligand-binding function of EGFR .
In this research, the N-glycan chains at Asn-467 likely kept the two antiparallel β-sheets positioned in the relative stable orientation. The N-linked glycan chains could stringently limit the motion of the amide nitrogen of the glycosylation site, as well as rigidifying the residues proximal to the glycosylation site, thereby stabilizing the region surrounding the glycosylation site with increased overall protein stability simultaneously . At the same time, the N-linked glycans could also induce the formation of the secondary structure of the nascent polypeptide proximal to the glycosylation site, thereby enabling the protein to be folded into the desired conformation [48,49,50]. When this glycosylation site was removed via mutagenesis, the two antiparallel β-sheets topology may collapse due to the loss of the support of glycan chains, thereby leading to a dramatic decrement of enzyme activity, as proven by the 2-fold higher Km and 46-fold lower kcat/Km.
Another function of N-linked glycans at the site close to the linker between the catalytic domain and β-sandwich domain could enhance the overall protein stability, albeit with deterioration in enzymatic activity caused by steric hindrance to adjacent N-glycans. For Asn-362, the glycosylation site was located adjacent to the linker between the catalytic domain and β-sandwich domain, which may maintain the relative conformation of the two domains mentioned above via increased rigidity of peptides nearby. Meanwhile, the glycosylation performed as a steric hindrance to the adjacent N-glycans whose interaction with the substrate was interfered with, thereby leading to delayed access of the substrate into the catalytic center. This in turn negatively regulated the enzyme activity, especially for high polymerization substrates like inulin, but without efficiently affecting the activity with low-molecular-weight substrates like sucrose.
In general, all five glycosylation sites were located in the β-sandwich domain, away from the catalytic active center with Asp-30, Asp-159 and Glu-215 at the bottom of the active pocket, thus making it difficult for the glycan chains to gain access to the core amino acid residues in the active pocket, where they could produce direct interaction. Consequently, the regulation of enzyme activity by the N-linked glycan chains originated from the stabilization of the configuration of the β-sandwich domain and the ease of capturing or transferring the substrate to the active center, which altered the micro-environment closely around the catalytic pocket with higher substrate concentration. This is quite similar to the functions of CBM, i.e. the N-linked glycan chains endowed the β-sandwich domain with the functions of CBM. This discovery provided a novel aspect by which to fine-tune the enzyme with desired characters for industrial enzymes.
In this study, the effects of glycosylation on the rKcINU1 activity, secondary structure and thermal stability was studied, uncovering a unique feature of the enzyme rKcINU1 that five glycosylation sites were all located at the C-terminal β-sandwich domain, away from the catalytic domain, which could still efficiently regulate the enzyme activity by stabilizing the configuration of the β-sandwich region and favoring the substrate into the catalytic center, which increased the substrate concentration in the micro-environment around the active pocket. Therefore it was believed the N-linked glycan chains in the non-catalytic domain endowed the C-terminal β-sandwich domain with the functions of CBM. This discovery provided a new manipulation strategy for glycoengineering the inulinase from fungi as well as other industrial enzymes from GH32 family.
Availability of data and materials
We provide all the necessary data for the publication of this article. All additional data is present in the additional material.
Buffered Glycerol-complex Medium
Buffered Methanol-complex Medium
Yeast Nitrogen Base
Yeast Extract Peptone Dextrose Medium with Sorbitol
Zhang S, Yang F, Wang Q, Hua Y, Zhao ZK. High-level secretory expression and characterization of the recombinant Kluyveromyces marxianus inulinase. Process Biochem. 2012;47:151–5.
Arjomand MR, Ahmadian G, Habibi-Rezaei M, Hassanzadeh M, Karkhane AA, Moosavi-Movahedi AA, Amanlou M. The importance of the non-active site and non-periodical structure located histidine residue respect to the structure and function of exo-inulinase. Int J Biol Macromol. 2017;98:542–9.
Chi ZM, Zhang T, Cao TS, Liu XY, Cui W, Zhao CH. Biotechnological potential of inulin for bioprocesses. Bioresour Technol. 2011;102:4295–303.
Sainz-Polo MA, Ramirez-Escudero M, Lafraya A, Gonzalez B, Marin-Navarro J, Polaina J, Sanz-Aparicio J. Three-dimensional structure of Saccharomyces invertase: role of a non-catalytic domain in oligomerization and substrate specificity. J Biol Chem. 2013;288:9755–66.
Alvaro-Benito M, Polo A, Gonzalez B, Fernandez-Lobato M, Sanz-Aparicio J. Structural and kinetic analysis of Schwanniomyces occidentalis invertase reveals a new oligomerization pattern and the role of its supplementary domain in substrate binding. J Biol Chem. 2010;285:13930–41.
Alvaro-Benito M, Sainz-Polo MA, Gonzalez-Perez D, Gonzalez B, Plou FJ, Fernandez-Lobato M, Sanz-Aparicio J. Structural and kinetic insights reveal that the amino acid pair Gln-228/Asn-254 modulates the transfructosylating specificity of Schwanniomyces occidentalis beta-fructofuranosidase, an enzyme that produces prebiotics. J Biol Chem. 2012;287:19674–86.
Cuskin F, Flint JE, Gloster TM, Morland C, Basle A, Henrissat B, Coutinho PM, Strazzulli A, Solovyova AS, Davies GJ, Gilbert HJ. How nature can exploit nonspecific catalytic and carbohydrate binding modules to create enzymatic specificity. PNAS. 2012;109:20889–94.
Venditto I, Najmudin S, Luis AS, Ferreira LM, Sakka K, Knox JP, Gilbert HJ, Fontes CM. Family 46 carbohydrate-binding modules contribute to the enzymatic hydrolysis of xyloglucan and beta-1,3-1,4-glucans through distinct mechanisms. J Biol Chem. 2015;290:10572–86.
Hettle A, Fillo A, Abe K, Massel P, Pluvinage B, Langelaan DN, Smith SP, Boraston AB. Properties of a family 56 carbohydrate-binding module and its role in the recognition and hydrolysis of beta-1,3-glucan. J Biol Chem. 2017;292:16955–68.
Cao TS, Wang GY, Chi Z, Wang ZP, Chi ZM. Cloning, characterization and heterelogous expression of the INU1 gene from Cryptococcus aureus HYA. Gene. 2013;516:255–62.
Zhang LL, Chi Z, Tan MJ, Wang GY, Liu GL, Chi ZM. Cloning and characterization of an inulinase gene from the marine yeast Candida membranifaciens subsp. flavinogenie W14-3 and its expression in Saccharomyces sp. W0 for ethanol production. Mol Biotechnol. 2015;57:337–47.
Zhou J, Lu Q, Peng M, Zhang R, Mo M, Tang X, Li J, Xu B, Ding J, Huang Z. Cold-active and NaCl-tolerant exo-inulinase from a cold-adapted Arthrobacter sp. MN8 and its potential for use in the production of fructose at low temperatures. J Biosci Bioeng. 2015;119:267–74.
Hong SJ, Kim HJ, Kim JW, Lee DH, Seo JH. Optimizing promoters and secretory signal sequences for producing ethanol from inulin by recombinant Saccharomyces cerevisiae carrying Kluyveromyces marxianus inulinase. Bioprocess Biosyst Eng. 2014;38:263–72.
Kaszubaa K, Grzybek M, Orłowskia A, Dannea R, Róga T, Simonsd K, Coskun Ü, Vattulainen I. N-glycosylation as determinant of epidermal growth factor receptor conformation in membranes. PNAS. 2015;112:4334–9.
Subedi GP, Barb AW. The structural role of antibody N-glycosylation in receptor interactions. Structure. 2015;23:1573–83.
Qi F, Zhang W, Zhang F, Chen G, Liu W. Deciphering the effect of the different N-glycosylation sites on the decretion, activity, and stability of cellobiohydrolase I from Trichoderma reesei. Applied Environ Microb. 2014;80:3962–71.
Yang M, Yu XW, Zheng H, Sha C, Zhao C, Qian M, Xu Y. Role of N-linked glycosylation in the secretion and enzymatic properties of Rhizopus chinensis lipase expressed in Pichia pastoris. Microb Cell Factories. 2015;14:40.
Yu XW, Yang M, Jiang C, Zhang X, Xu Y. N-glycosylation engineering to improve the constitutive expression of Rhizopus oryzaelipase in Komagataella phaffii. J Agric Food Chem. 2017;65:6009–15.
Mendez-Yanez A, Beltran D, Campano-Romero C, Molinett S, Herrera R, Moya-Leon MA, Morales-Quintana L. Glycosylation is important for FcXTH1 activity as judged by its structural and biochemical characterization. Plant Physiol Biochem. 2017;119:200–10.
Fonseca-Maldonado R, Vieira DS, Alponti JS, Bonneil E, Thibault P, Ward RJ. Engineering the pattern of protein glycosylation modulates the thermostability of a GH11 xylanase. J Biol Chem. 2013;288:25522–34.
Gusakov AV, Dotsenko AS, Rozhkova AM, Sinitsyn AP. N-linked glycans are an important component of the processive machinery of cellobiohydrolases. Biochimie. 2017;132:102–8.
Dotsenko AS, Gusakov AV, Volkov PV, Rozhkova AM, Sinitsyn AP. N-linked glycosylation of recombinant cellobiohydrolase I (Cel7A) from Penicillium verruculosum and its effect on the enzyme activity. Biotechnol Bioeng. 2016;113:283–91.
Qin Y, Qu Y. Asn124 of Cel5A from Hypocrea jecorina not only provides the N-glycosylation site but is also essential in maintaining enzymatic activity. BMB Rep. 2014;47:256–61.
Han M, Wang X, Yan G, Wang W, Tao Y, Liu X, Cao H, Yu X. Modification of recombinant elastase expressed in Pichia pastoris by introduction of N-glycosylation sites. J Biotechnol. 2014;171:3–7.
Wei W, Chen L, Zou G, Wang QF, Yan X, Zhang J, Wang CS, Zhou ZH. N-glycosylation affects the proper folding, enzymatic characteristics and production of a fungal β-glucosidase. Biotechnol Bioeng. 2013;110:3075–84.
Maestre-Reyna M, Liu WC, Jeng WY, Lee CC, Hsu CA, Wen TN, Wang AH, Shyur LF. Structural and functional roles of glycosylation in fungal laccase from Lentinus sp. PLoS One. 2015;10:e0120601.
Cai X, Thinn AMM, Wang Z, Shan H, Zhu J. The importance of N-glycosylation on β3 integrin ligand binding and conformational regulation. Sci Rep. 2017;7:4656.
Kar B, Verma P, den Haan R, Sharma AK. Effect of N-linked glycosylation on the activity and stability of a β-glucosidase from Putranjiva roxburghii. Int J Biol Macromol. 2018;112:490–8.
Lee SM, Booe JM, Gingell JJ, Sjoelund V, Hay DL, Pioszak AA. N-glycosylation of asparagine 130 in the extracellular domain of the human calcitonin receptor significantly increases peptide hormone affinity. Biochemistry. 2017;56:3380–93.
Gao L, Gao F, Wang L, Geng C, Chi L, Zhao J, Qu Y. N-glycoform diversity of cellobiohydrolase I from Penicillium decumbens and synergism of nonhydrolytic glycoform in cellulose degradation. J Biol Chem. 2012;287:15906–15.
Le Roy K, Verhaest M, Rabijns A, Clerens S, Van Laere A, Van den Ende W. N-glycosylation affects substrate specificity of chicory fructan 1-exohydrolase: evidence for the presence of an inulin binding cleft. New Phytol. 2007;176:317–24.
Tauzin AS, Sulzenbacher G, Lafond M, Desseaux V, Reca IB, Perrier J, Bellincampi D, Fourquet P, Leveque C, Giardina T. Functional characterization of a vacuolar invertase from Solanum lycopersicum: post-translational regulation by N-glycosylation and a proteinaceous inhibitor. Biochimie. 2014;101:39–49.
Plascencia-Espinosa M, Santiago-Hernández A, Pavón-Orozco P, Vallejo-Becerra V, Trejo-Estrada S, Sosa-Peinado A, Benitez-Cardoza CG, Hidalgo-Lara ME. Effect of deglycosylation on the properties of thermophilic invertase purified from the yeast Candida guilliermondii MpIIIa. Process Biochem. 2014;49:1480–7.
Ramírez-Escudero M, Gimeno-Pérez M, González B, Linde D, Merdzo Z, Fernández-Lobato M, Sanz-Aparicio J. Structural analysis of β-Fructofuranosidase from Xanthophyllomyces dendrorhous reveals unique features and the crucial role of N-glycosylation in oligomerization and activity. J Biol Chem. 2016;291:6843–57.
Ma JY, Cao HL, Tan HD, Hu XJ, Liu WJ, Du YG, Yin H. Cloning, expression, characterization, and mutagenesis of a thermostable exoinulinase from Kluyveromyces cicerisporus. Appl Biochem Biotechnol. 2016;178:144–58.
Jiang H, Yuan H, Qu Y, Liang Y, Jiang B, Wu Q, Deng N, Liang Z, Zhang L, Zhang Y. Preparation of hydrophilic monolithic capillary column by in situ photo-polymerization of N-vinyl-2-pyrrolidinone and acrylamide for highly selective and sensitive enrichment of N-linked glycopeptides. Talanta. 2016;146:225–30.
Ge F, Zhu L, Aang A, Song P, Li W, Tao Y, Du G. Recent advances in enhanced enzyme activity, thermostability and secretion by N-glycosylation regulation in yeast. Biotechnol Lett. 2018;40:847–54.
Lizak C, Fan YY, Weber TC, Aebi M. N-linked glycosylation of antibody fragments in Escherichia coli. Bioconjug Chem. 2011;22:488–96.
Malaby HL, Kobertz WR. The middle X residue influences cotranslational N-glycosylation consensus site skipping. Biochemistry. 2014;53:4884–493.
Huang YW, Yang HI, Wu YT, Hsu TL, Lin TW, Kelly JW, Wong CH. Residues comprising the enhanced aromatic sequon influence protein N-glycosylation efficiency. J Am Chem Soc. 2017;139:12947–55.
Chen CD. Specificity of O-glycosylation in enhancing the stability and cellulose binding affinity of family 1 carbohydratebinding modules. PNAS. 2014;111:7612–7.
Badino SF, Bathke JK, Sorensen TH, Windahl MS, Jensen K, Peters GHJ, Borch K, Westh P. The influence of different linker modifications on the catalytic activity and cellulose affinity of cellobiohydrolase Cel7A from Hypocrea jecorina. Protein Eng Des Sel. 2017;30:495–501.
Ma L, Ma Q, Cai R, Zong Z, Du L, Guo G, Zhang Y, Xiao D. Effect of β-mannanase domain from Trichoderma reesei on its biochemical characters and synergistic hydrolysis of sugarcane bagasse. J Sci Food Agric. 2017;98:2540.
Kajiwara H, Toda M, Mine T, Nakada H, Wariishi H, Yamamoto T. Visualization of sialic acid produced on bacterial cell surfaces by lectin staining. Microbes Environ. 2010;25:152–5.
Jenkins GA, Figueira M, Kumar GA, Sweetman WA, Makepeace K, Pelton SI, Moxon R, Hood DW. Sialic acid mediated transcriptional modulation of a highly conserved sialometabolism gene cluster in Haemophilus influenzae and its effect on virulence. BMC Microbiol. 2010;10:48.
Kaushik S, Mohanty D, Surolia A. Molecular dynamics simulations on pars intercerebralis major peptide-C (PMP-C) reveal the role of glycosylation and disulfide bonds in its enhanced structural stability and function. J Biomol Struct Dyn. 2012;29:905–20.
Imperiali B, O’Connor SE. Effect of N-linked glycosylation on glycopeptide and glycoprotein structure. Curr OpinChem Biol. 1999;3:643–9.
Imperiali B, Rickert KW. Conformational implications of asparagine-linked glycosylation. PNAS. 1995;92:97–101.
O’Connor SE, Imperiali B. A molecular basis for glycosylation induced conformational switching. Chem Biol. 1998;5:427–37.
Hanson SR, Culyba EK, Hsu TL, Wong CH, Kelly JW, Powers ET. The core trisaccharide of an N-linked glycoprotein intrinsically accelerates folding and enhances stability. PNAS. 2009;106:3131–6.
The authors gratefully acknowledge Prof. Zhongfu Wang and Prof. Ying Zhang (Northwest University, China) for the ESI-MS analysis of N-linked glycan chains.
This work was supported by the Natural Science Foundation of China (31670803, 31500040), the National Key Research and Development Project of China (2017YFD0200900) Subject 2 (2017YFD0200902), the Natural Science Foundation of Liaoning Province (20170540027).
Ethics approval and consent to participate
Consent for publication
All authors consent for publication.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Materials and Methods about Site-Directed Mutagenesis, Protein Expression and Purification, Glycosidase Treatment, Zymogram Analysis, Enzymes Activity Assays and Mass Spectrometry Analysis. Figure S1. LC-MS/MS spectrum analysis of the glycosylation sites within the wild-type rKcINU1. The sample was treated by trypsin followed with PNGase F digestion and then subjected to LC-MS/MS. The peak with m/z of 976.95, 694.66, 942.49, 1188.59 and 817.90 confirmed the presence of deamidation modification on (A) Asn-362, (B) Asn-370, (C) Asn-399, (D) Asn-467 and (E) Asn-526, respectively (substitution of Asn with Glu residue). Each peptide sequence was shown in the panel and the arrows indicated the N-glycosylation sites. Figure S2. LC-MS/MS spectrum analysis of the glycosylation sites within the variant Mut. The sample was treated by trypsin followed with PNGase F digestion and then subjected to LC-MS/MS. The peak with m/z of 621.32, 1298.60, 1298.60, 990.45, 990.45 and 874.37 confirmed the presence of deamidation modification on (A) Asn-9, (B) Asn-147, Asn-153 (C) Asn-197, Asn-203 and (D) Asn-233, respectively (substitution of Asn with Glu residue). Each peptide sequence was shown in the panel and the arrows indicated the N-glycosylation sites. Figure S3. Positive ion ESI-MS analysis of N-linked glycan chains for the single-sited mutants of (A) N362Q, (B) N370Q, (C) N399Q, (D) N467Q and (E) N526Q, respectively. The peak of m/z 1905.17 confirmed the presence of the high mannose oligosaccharide (Man)7(GlcNAc)2. Series of peaks by 162 Da (the mass of anhydrohexose) demonstrated the glycoform heterogeneity. Solid squares and circles represented N-acetylglucosamine (GlcNAc) and Mannose (Man), respectively. Table S1. Primers of the site-directed mutagenesis. Table S2. The enzyme activity of the wildtype rKcINU1 treated with Endo F1
About this article
Cite this article
Ma, J., Li, Q., Tan, H. et al. Unique N-glycosylation of a recombinant exo-inulinase from Kluyveromyces cicerisporus and its effect on enzymatic activity and thermostability. J Biol Eng 13, 81 (2019). https://doi.org/10.1186/s13036-019-0215-y